数字信号处理实验-IIR滤波器设计

《数字信号处理实验-IIR滤波器设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验四IIR数字滤波器的设计(1)，，，,;设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。解：程序：9clear;9fc=300;fr=200;fs=1000;rp=0.8;rs=20;wc=2*fs*tan(2*pi*fc/(2*fs));wt=2*fs*tan(2*pi*fr/(2*fs));[N,wn]=cheb1ord(wc,wt,rp,rs,'s');[B,A]=cheby1(N,rp,wn,'high','s');[bz,az]=bilinear(B,A,fs);[h,w]=fr

2、eqz(bz,az);f=w*fs/(2*pi);plot(f,20*log10(abs(h)));axis([0,fs/2,-80,10]);grid;xlabel('频率/Hz');ylabel('幅度/dB');99分析：f=200Hz时阻带衰减大于30dB，通过修改axis([0,fs/2,-80,10])为axis([200,fs/2,-1,1])发现通带波动rs满足<0.8。9bz=[0.0262-0.10470.1570-0.10470.0262]az=[1.00001.52891.65370.

3、94520.2796]系统函数为：（2），，，,;分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。解：程序：9clear;fs=1000;fc=200;fr=300;rp=1;rs=25;%脉冲响应不变法wp=2*pi*fc;ws=2*pi*fr;[N,wn]=buttord(wp,ws,rp,rs,'s');[b1a1]=butter(N,wn,'s');[bz1,az1]=impinvar(b1

4、,a1,fs);[h1,w]=freqz(bz1,az1);%双线性变换法wp=2*fs*tan(2*pi*fc/fs/2);ws=2*fs*tan(2*pi*fr/fs/2);[N,wn]=buttord(wp,ws,rp,rs,'s');[b2a2]=butter(N,wn,'s');[bz2,az2]=bilinear(b2,a2,fs);[h2,w]=freqz(bz2,az2);f=w/(2*pi)*fs;figure;plot(f,abs(h1),'-.r',f,abs(h2),'-b');gri

5、d;xlabel('频率/Hz');ylabel('幅度');legend('脉冲响应不变法','双线性变换法');title('巴特沃思低通滤波器，线性幅度谱');99bz1=[0.00000.00020.01530.09950.14440.06110.00750.00020.00000]az1=[1.0000-1.91992.5324-2.20531.3868-0.63090.2045-0.04500.0060-0.0004]因此脉冲响应不变法的系统函数为:bz2=[0.01790.10720.26810

6、.35750.26810.10720.0179]az2=[1.0000-0.60190.9130-0.29890.1501-0.02080.0025]因此双线性变换法的系统函数为:分析：脉冲响应不变法的N=9，双线性变换法的N=6，由图知它们都满足要求，但脉冲响应的衰减较快，双线性变换的过渡带窄一些，且阶数比脉冲小，容易实现。（3）利用双线性变换法分别设计满足下列指标的巴特沃思型、切比雪夫型和椭圆型数字低通滤波器，并作图验证设计结果：，，，,。解：程序：9clear;fs=8000;fc=1200;fr=20

7、00;rp=0.5;rs=40;%巴特沃思低通滤波器wp=2*fs*tan(2*pi*fc/fs/2);ws=2*fs*tan(2*pi*fr/fs/2);[N,wn]=buttord(wp,ws,rp,rs,'s');[b1a1]=butter(N,wn,'s');[bz1,az1]=bilinear(b1,a1,fs);[h1,w]=freqz(bz1,az1);H1=20*log10(abs(h1));f=w/(2*pi)*fs;figure;plot(f,H1);%对数幅度谱axis([0,fs/2,

8、-100,10]);grid;xlabel('频率/Hz');ylabel('幅度');title('巴特沃思低通滤波器，对数幅度谱');%切比雪夫低通滤波器wc=2*fs*tan(2*pi*fc/(2*fs));wt=2*fs*tan(2*pi*fr/(2*fs));[N,wn]=cheb1ord(wc,wt,rp,rs,'s');[b2,a2]=cheby1(N,rp,wn,'low','s')