1 #!/usr/bin/python

     2 ## 4.4

     3 import os, sys

     4 sys.path += [os.getcwdu() + os.sep + ".." + os.sep + "python"]

     5 

     6 import scipy as sp

     7 from scipy import signal

     8 from dscr import dscr

     9 

    10 F = sp.array([[0,0],[1,-0.1]])

    11 G = sp.array([[0.1],[0]])

    12 C = sp.array([0,1])

    13 D = 0

    14 Ts = 0.2

    15 

    16 H = signal.lti(F, G, C, D)

    17 

    18 H = dscr(H,Ts)

    19 n, d = H.num, H.den

    20 

    21 print "Numerator: ", n

    22 print "Denominator: ", d

    23 print "Sampling Time: ", Ts