SHOW:
|
|
- or go back to the newest paste.
| 1 | - | def analitStep((A, B, C, D), T): |
| 1 | + | A = [-10 0 -100; 2 -2 0; 0 0.1 0.2]; |
| 2 | - | ss = StateSpace(A, B, C, D) |
| 2 | + | B = [10; 0; 0]; |
| 3 | - | |
| 3 | + | C = [0 0 1]; |
| 4 | - | T = arange(T[0], T[1], T[2]) |
| 4 | + | D = [0]; |
| 5 | - | J, V = eig(A) |
| 5 | + | syms t |
| 6 | - | J = mt(J * eye(3, 3)) |
| 6 | + | wn = C * expm(A * t) * B; |
| 7 | - | V = mt(V) |
| 7 | + | tn = [0:0.01:7]; |
| 8 | - | y = ([], []) |
| 8 | + | plot(tn, yn), grid on; |
