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