freqtransform.m

Signal and Sytem Using Matlab 书的源代码 注释详尽 m文件

%频率变换，低通、高通、带通、带阻函数分别为lp2lp,lp2hp,lp2bp,lp2bs
%page_353 upper
%怎样才能和书上作出一样的图呢？？？？？？？？？？？？？？？？？？
[z,p,k]=buttap(3);%3 pole filter
[num,den]=zp2tf(z,p,k);%convert to polynoials
wb=5*2*pi;  %new bandwidth in rad/sec
[num,den]=lp2lp(num,den,wb);%频率转换
f=0:15/200:15;
w=2*pi*f;
H=tf(num,den);
bode(H,w);

% 以下的画图方法是错误的，因为H是一个tf对象，不是一个可求值的表达式
% subplot(2,1,1);plot(w,abs(H));      
% subplot(2,1,2);plot(w,angle(H)*180/pi);