|
0
|
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)
|