SHOW:
|
|
- or go back to the newest paste.
| 1 | - | def analitMihailovDis(): |
| 1 | + | def Nyquist((A, B, C, D), T): |
| 2 | - | y = ([], []) |
| 2 | + | sys = ss(A, B, C, D) |
| 3 | - | for i in map(lambda x: x / 100.0, range(0, 314)): |
| 3 | + | real, imag, freq = nyquist(sys, omega=arange(T[0], T[1], T[2]), Plot=False) |
| 4 | - | yy = exp(3 * i * 1j) - 2.887 * exp(2 * i * 1j) + 2.775 * exp(i * 1j) - 0.89 |
| 4 | + | plot(real, imag) |
| 5 | - | y[0].append(yy.real) |
| 5 | + | title('Nyquist construction')
|
| 6 | - | y[1].append(yy.imag) |
| 6 | + | xlabel('Real')
|
| 7 | - | plot(y[0], y[1]) |
| 7 | + | ylabel('Imag')
|
| 8 | - | title('Godograf Mihailova')
|
| 8 | + | |
| 9 | - | xlabel('t, s.')
|
| 9 | + |
