--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/freq/nmp.py	Fri May 27 14:24:59 2011 +0530
@@ -0,0 +1,26 @@
+#!/usr/bin/python
+# 5.2
+
+import scipy as sp
+from pylab import *
+
+omega = linspace(0, pi, 100)
+ejw = exp(-1j*omega)
+G1 = 1.5*(1-0.4*ejw)
+mag1 = abs(G1)
+angle1 = angle(G1)*180/pi
+G2 = -0.6*(1-2.5*ejw)
+mag2 = abs(G2)
+angle2 = angle(G2)*180/pi
+
+subplot(2,1,1)
+plot(omega,mag1,omega,mag2,'--')
+ylabel('Magnitude')
+
+subplot(2,1,2) 
+plot(omega,angle1,omega,angle2,'--')
+xlabel('w (rad/s)')
+ylabel('Phase')
+
+show()
+#axis tight, label('',18,'w (rad/s)','Phase',18)