nivs

Untitled

Jan 23rd, 2013
148
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. def analitStepDis((A, B, C, D), T):
  2.     T = arange(T[0], T[1], T[2])
  3.     J, V = eig(A)
  4.     J = mt(J * eye(3, 3))
  5.     V = mt(V)
  6.     y = ([], [])
  7.     for i in range(len(T)):
  8.         y[0].append(T[i])
  9.         yy = mt(((0, 0, 0))).transpose()
  10.         for j in range(i - 1):
  11.             yy = yy + matrix_power(A, i - 1 - j) * B
  12.         yy = C * yy
  13.         y[1].append(array(yy[0])[0][0])
  14.     plot(y[0], y[1])
  15.     title('Step construction')
  16.     xlabel('t, s.')
  17.     ylabel('y(t)')
  18.     grid(color='g', linestyle='--')
  19.     show()
Advertisement
Add Comment
Please, Sign In to add comment