SHOW:
|
|
- or go back to the newest paste.
| 1 | - | def eqS((A, B, C, D), (An, Bn, Cn, Dn)): |
| 1 | + | def Initial((A, B, C, D), x0, T): |
| 2 | - | A = mt(A); B = mt(B); C = mt(C); D = mt(D) |
| 2 | + | ss = StateSpace(A, B, C, D) |
| 3 | - | An = mt(An); Bn = mt(Bn); Cn = mt(Cn); Dn = mt(Dn) |
| 3 | + | T = arange(0, 7, 0.001) # шаг построения |
| 4 | - | b = B[:, 0] |
| 4 | + | x0 = array([1.0, 1.0, 1.0]) # начальные условия |
| 5 | - | bn = Bn[:, 0] |
| 5 | + | |
| 6 | - | n = len(A) |
| 6 | + | iss = initial(ss, T, x0) |
| 7 | - | Sy = b |
| 7 | + | |
| 8 | - | Sny = bn |
| 8 | + | plot(iss[1], iss[0]) |
| 9 | - | for i in range(1, n): |
| 9 | + | title('initial() construction')
|
| 10 | - | Apow = copy(A) |
| 10 | + | xlabel('t, s.')
|
| 11 | - | Anpow = copy(An) |
| 11 | + | ylabel('y(t)')
|
| 12 | - | if not i == 1: |
| 12 | + | grid(color='g', linestyle='--') |
| 13 | - | for j in range(i-1): |
| 13 | + | show() |
