1 #!/usr/bin/python

     2 # 7.4

     3 import os, sys

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

     5 

     6 import pylab as pl

     7 from scipy import signal

     8 from bode import bode

     9 

    10 w = pl.linspace(0.001,pl.pi,10000)

    11 a = pl.linspace(0.001,0.999,100)

    12 

    13 zero = a

    14 pole = 0.9*a

    15 A = pl.vstack((pl.ones_like(zero), -zero))

    16 A = A.T

    17 B = pl.vstack((pl.ones_like(pole), -pole))

    18 B = B.T

    19 sys = map(signal.lti, A, B)

    20 

    21 

    22 

    23 #pl.show()