1 #!/usr/bin/python

     2 # 7.3

     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 pol1 = [1, -0.9]

    11 pol2 = [1, -0.8]

    12 G1 = signal.lti(pol1,[1, 0])

    13 G2 = signal.lti(pol2,[1, 0])

    14 w = pl.linspace(0.001,0.5,1000)

    15 bode(G1, w)

    16 bode(G2, w)

    17 pl.figure()

    18 G = signal.lti(pol1,pol2)

    19 bode(G,w)

    20 

    21 pl.show()