equal
deleted
inserted
replaced
|
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) |