ms20_4.m

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

%脉冲响应不变法用于高通电路产生混叠
%数字滤波器指标
wp=0.4*pi;            %滤波器的通带截止频率
ws=0.25*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和截止频率
[z0,p0,k0]=cheb2ap(n,As);   %设计归一化的模拟原型滤波器
ba1=k0*real(poly(z0));          %求原型滤波器系数b
aa1=real(poly(p0));             %求原型滤波器系数a
[ba,aa]=lp2hp(ba1,aa1,Omgc);    %变换为模拟滤波器
[bd,ad]=impinvar(ba,aa,Fs); %用脉冲响应不变法计算数字滤波器系数

t=0:T:(30*T);nt=length(t);
subplot(2,1,1),
ha=impulse(ba,aa,t);plot(t,ha*T,'r'),hold on
h=impz(bd,ad,nt);stem(t,h);

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