samp5_9.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 15 行

%Samp5_9
[z,p,k]=buttap(6);   %设计阶数为6的Butterworth模拟原型滤波器
[b,a]=zp2tf(z,p,k);  %将零点极点增益形式转换为传递函数形式
[H,w]=freqs(b,a,0:0.01:2); %计算指定频率点的复数频率响应
subplot(2,1,1),plot(w,abs(H).^2); %绘制平方幅频响应
xlabel('w/wc');ylabel('Butterworth |H(jw)|^2');
title('Btterworth低通原型滤波器(wc=1)')
w1=0.7;w2=1.5;   %带阻滤波器的下限及上限频率
w0=sqrt(w1*w2);bw=w2-w1;   %计算中心点频率和阻带宽度
[bt,at]=lp2bs(b,a,w0,bw);   %频率转换
[Ht,wt]=freqs(bt,at,0:0.01:2);  %计算带阻滤波器指定点的复数频率响应
subplot(2,1,2),plot(wt,abs(Ht).^2);%绘制平方幅频响应
xlabel('w/wc');ylabel('Butterworth |H(jw)|^2');
title('Butterworth带阻滤波器(wc=0.7~1.5) ')