|
0
|
1 |
#!/usr/bin/env python
|
|
|
2 |
# 9.16
|
|
|
3 |
import os, sys
|
|
|
4 |
sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]
|
|
|
5 |
|
|
|
6 |
import scipy as sp
|
|
|
7 |
from zpowk import zpowk
|
|
|
8 |
from pp_im import pp_im
|
|
|
9 |
|
|
|
10 |
B = 1
|
|
|
11 |
A = [1 -1]
|
|
|
12 |
k = 1
|
|
|
13 |
zk = zpowk(k)
|
|
|
14 |
Ts = 1
|
|
|
15 |
phi = sp.array([1, -0.5])
|
|
|
16 |
|
|
|
17 |
Delta = 1 # Choice of internal model of step
|
|
|
18 |
Rc, Sc, Tc, gamm = pp_im(B, A, k, phi, Delta)
|
|
|
19 |
|
|
|
20 |
# simulation parameters for stb_disc
|
|
|
21 |
st = 1 # desired step change
|
|
|
22 |
t_init = 0 # simulation start time
|
|
|
23 |
t_final = 20 # simulation end time
|
|
|
24 |
xInitial = [0, 0]
|
|
|
25 |
C = 0
|
|
|
26 |
D = 1
|
|
|
27 |
N_var = 0
|