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 ejw = exp(-1j*omega)

     9 G1 = 1.5*(1-0.4*ejw)

    10 mag1 = abs(G1)

    11 angle1 = angle(G1)*180/pi

    12 G2 = -0.6*(1-2.5*ejw)

    13 mag2 = abs(G2)

    14 angle2 = angle(G2)*180/pi

    15 

    16 subplot(2,1,1)

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

    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)