实验八IIR滤波器的算法

《实验八IIR滤波器的算法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验八IIR滤波器的算法一、实验目的和任务1．理解并掌握IIR数字滤波器的机理；2.掌握设计IIR数字滤波器的原理和方法；3.观察IIR滤波器的频域特性。二、实验原理1、设计原理（1）脉冲响应不变法所谓脉冲相应不变法就是使数字滤波器的单位脉冲响应序列h(n)等于模拟滤波器的单位冲击响应和(t)的采样值。（2）双向性变换法：由于s平面和z平面的单值双线性映射关系为s＝其中T为采样周期。因此，若已知模拟滤波器的传递函数，将上式代入即可得到数字滤波器的系统函数H(z)。在双线性变换中，模拟角频率和数字角频率的变换关系为可见，Ω和w之间的变换关系为非线性的。3   三、IIR滤波器的

2、设计参考实例一：IIR低通滤波（Butterworth）%IIRLowpassUseButterworthclear;fs=20;fpass=4;fstop=5;Ap=0.5;As=10;wp=2*pi*fpass/fs;ws=2*pi*fstop/fs;omegap=tan(wp/2);omegas=tan(ws/2);ep=sqrt(10^(Ap/10)-1);es=sqrt(10^(As/10)-1);N=ceil(log(es/ep)/log(omegas/omegap));omega0=omegap/ep^(1/N);K=floor(N/2);fori=1:K   

3、theta(i)=pi*(N-1+2*i)/(2*N);endfori=1:K   G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2);endfori=1:K   a1(i)=2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2);endfori=1:K   a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);endifK<(N/2)   G0=omega0/(omega0+1);

4、a0=(omega0-1)/(omega0+1);endw=0:pi/300:pi;Hw2=1./(1+(tan(w/2)/omega0).^(2*N));plot(w/pi,Hw2);grid;参考实例二：IIR高通滤波（Butterworth）%IIRHightpassUseButterworthclear;fs=20;fpass=5;fstop=4;Ap=0.5;As=10;wp=2*pi*fpass/fs;ws=2*pi*fstop/fs;omegap=cot(wp/2);omegas=cot(ws/2);ep=sqrt(10^(Ap/10)-1);es=sqrt(1

5、0^(As/10)-1);N=ceil(log(es/ep)/log(omegas/omegap));omega0=omegap/ep^(1/N);K=floor(N/2);fori=1:K   theta(i)=pi*(N-1+2*i)/(2*N);endfori=1:K   G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2);endfori=1:K   a1(i)=-2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2);endfori=1:K   a2(i)=(1+2*omeg

6、a0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2);endifK<(N/2)   G0=omega0/(omega0+1);a0=-(omega0-1)/(omega0+1);endw=(0+eps):pi/300:pi;Hw2=1./(1+(cot(w/2)/omega0).^(2*N));plot(w/pi,Hw2);grid; 参考实例三：IIR带通滤波（Butterworth）%IIRBandpassUseButterworthclear;fs=20;fpa=2;fpb=4;fsa=1.5;fs

7、b=4.5;Ap=0.0877;As=16.9897;wpa=2*pi*fpa/fs;wpb=2*pi*fpb/fs;wsa=2*pi*fsa/fs;wsb=2*pi*fsb/fs;c=sin(wpa+wpb)/(sin(wpa)+sin(wpb));omegap=abs((c-cos(wpb))/sin(wpb));omegasa=(c-cos(wsa))/sin(wsa);omegasb=(c-cos(wsb))/sin(wsb);omegas=min(abs(omegasa),abs(omegas