SHOW:
|
|
- or go back to the newest paste.
| 1 | - | def analitMihailovDis(): |
| 1 | + | def analitStepDis((A, B, C, D), T): |
| 2 | T = arange(T[0], T[1], T[2]) | |
| 3 | - | for i in map(lambda x: x / 100.0, range(0, 314)): |
| 3 | + | J, V = eig(A) |
| 4 | - | yy = exp(3 * i * 1j) - 2.887 * exp(2 * i * 1j) + 2.775 * exp(i * 1j) - 0.89 |
| 4 | + | J = mt(J * eye(3, 3)) |
| 5 | - | y[0].append(yy.real) |
| 5 | + | V = mt(V) |
| 6 | - | y[1].append(yy.imag) |
| 6 | + | |
| 7 | - | plot(y[0], y[1]) |
| 7 | + | for i in range(len(T)): |
| 8 | - | title('Godograf Mihailova')
|
| 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() |
