ms20_3_2.m

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

%脉冲响应不变法用于低通电路
%数字滤波器指标
wp=0.25*pi;           %滤波器的通带截止频率
ws=0.4*pi;            %滤波器的阻带截止频率
Rp=1;As=40;           %输入滤波器的通阻带衰减指标
%转换为模拟原型滤波器指标
Fs=2000;T=1/Fs;
Omgp=wp*Fs;Omgs=ws*Fs;
%模拟滤波器计算
[n,Omgn]=cheb2ord(Omgp,Omgs,Rp,As,'s')%计算阶数n和截止频率
[ba,aa]=cheby2(n,As,Omgn,'s');          %求模拟滤波器系数
[bd,ad]=impinvar(ba,aa,Fs);    %用脉冲响应不变法计算数字滤波器系数
%模拟滤波器与数字滤波器的冲激响应
t=0:T:(30*T);nt=length(t);
ha=impulse(ba,aa,t);
h=impz(bd,ad,nt);
subplot(2,1,1),plot(t,ha*T,'r'),hold on
stem(t,h,'k')
title('模拟与数字滤波器的冲激响应');

%模拟滤波器与数字滤波器的幅频响应
wb=[0:Fs]*2*pi;        %为作图建立频率向量 
Ha=freqs(ba,aa,wb);    %计算模拟频率响应
H=freqz(bd,ad,wb/Fs);  %计算数字频率响应
subplot(2,1,2),
plot(wb/(2*pi),abs(Ha)/max(abs(Ha)),'r'),hold on
plot(wb/(2*pi),abs(H)/max(abs(H)),'k');   
title('脉冲响应不变法用于低通滤波');
axis([0,Fs,-0.1,1.1]);
ylabel('幅度');xlabel('f(Hz)');