def _plotArrays(self):
import gist
for array in self._getArrays():
gist.plg(*array)
gist.logxy(self.xlog, self.ylog)
def bode(w,H,win=0,frame=0,lcolor='blue',color='black',tcolor='black',freq='rad'):
"""Plot a bode plot of the transfer function H as a function of w.
"""
if freq == 'Hz':
w = w /2.0 / pi
subplot(2,1,win,hsep=120,frame=frame,color=color)
gist.plsys(1)
gist.plg(20*numpy.log10(abs(H)),w,type='solid',color=lcolor,marks=0)
gist.logxy(1,0)
gist.gridxy(1,1)
if freq == 'Hz':
xlabel('Frequency (Hz)',color=tcolor,deltay=-0.005)
else:
xlabel('Frequency (rad/s)',color=tcolor,deltay=-0.005)
ylabel('Magnitude (dB)',color=tcolor,deltax=-0.005)
title("Bode Plot",color=tcolor)
gist.plsys(2)
gist.plg(180/pi*numpy.unwrap(MLab.angle(H)),w,type='solid',color=lcolor,marks=0)
gist.logxy(1,0)
gist.gridxy(1,1)
if freq == 'Hz':
xlabel('Frequency (Hz)',color=tcolor,deltay=-0.005)
else:
xlabel('Frequency (rad/s)',color=tcolor,deltay=-0.005)
ylabel('Phase (deg.)',color=tcolor,deltax=-0.005)
def _plotArrays(self):
import gist
for array in self._getArrays():
gist.plg(*array)
gist.logxy(self.xlog, self.ylog)
def bode(w,H,win=0,frame=0,lcolor='blue',color='black',tcolor='black',freq='rad'):
"""Plot a bode plot of the transfer function H as a function of w.
"""
if freq == 'Hz':
w = w /2.0 / pi
subplot(2,1,win,hsep=120,frame=frame,color=color)
gist.plsys(1)
gist.plg(20*numpy.log10(abs(H)),w,type='solid',color=lcolor,marks=0)
gist.logxy(1,0)
gist.gridxy(1,1)
if freq == 'Hz':
xlabel('Frequency (Hz)',color=tcolor,deltay=-0.005)
else:
xlabel('Frequency (rad/s)',color=tcolor,deltay=-0.005)
ylabel('Magnitude (dB)',color=tcolor,deltax=-0.005)
title("Bode Plot",color=tcolor)
gist.plsys(2)
gist.plg(180/pi*numpy.unwrap(numpy.angle(H)),w,type='solid',color=lcolor,marks=0)
gist.logxy(1,0)
gist.gridxy(1,1)
if freq == 'Hz':
xlabel('Frequency (Hz)',color=tcolor,deltay=-0.005)
else:
xlabel('Frequency (rad/s)',color=tcolor,deltay=-0.005)
ylabel('Phase (deg.)',color=tcolor,deltax=-0.005)
