ms22_1n3.m

《数字信号处理实验(MATLAB版)》这

%双线性变换法设计数字低通滤波器
%数字滤波器指标
fp=150;fs=250;Fs=800;T=1/Fs;
wp=fp/Fs*2;             %数字滤波器的通带截止频率
ws=fs/Fs*2;             %数字滤波器的阻带截止频率
Rp=3;As=20;                %输入滤波器的通阻带衰减指标
[n,wc]=buttord(wp,ws,Rp,As) %计算阶数n和截止频率
%方法3:直接求数字滤波器系数
[bd3,ad3]=butter(n,wc)  
[H3,w3]=freqz(bd3,ad3);   
dbH3=20*log10(abs(H3)/max(abs(H3))); %化为分贝值
%
subplot(1,2,1),plot(w3/2/pi*Fs,dbH3,'k');
title('直接法  幅度响应(dB)');axis([0,Fs/2,-40,5]); 
ylabel('dB');xlabel('频率(Hz)');
set(gca,'XTickMode','manual','XTick',[0,fp,fs,Fs/2]);
set(gca,'YTickMode','manual','YTick',[-50,-20,-3,0]);grid
subplot(1,2,2),plot(w3/2/pi*Fs,angle(H3)/pi*180,'k');
title('相位响应');axis([0,Fs/2,-180,180]);
ylabel('\phi');xlabel('频率(Hz)');
set(gca,'XTickMode','manual','XTick',[0,fp,fs,Fs/2]);
set(gca,'YTickMode','manual','YTick',[-180,0,180]);grid