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