1 #!/usr/bin/python

     2 # 5.2

     3 

     4 import scipy as sp

     5 from pylab import *

     6 

     7 omega = linspace(0, pi, 100)

     8 g1 = 0.5 / (cos(omega)-0.5+1j*sin(omega))

     9 mag1 = abs(g1)

    10 angle1 = angle(g1) * 180/pi

    11 g2 = (0.5+0.5*cos(omega)-1.5*1j*sin(omega)) \

    12      * 0.25 / (1.25-cos(omega))

    13 mag2 = abs(g2)

    14 angle2 = angle(g2) * 180/pi

    15 subplot(2,1,1)

    16 plot(omega,mag1,omega,mag2,'--')

    17 

    18 ylabel('Magnitude')

    19 

    20 subplot(2,1,2) 

    21 plot(omega,angle1,omega,angle2,'--')

    22 xlabel('w (rad/s)')

    23 ylabel('Phase')

    24 

    25 show()

    26 #axis tight, label('',18,'w (rad/s)','Phase',18)