数字信号处理第二次实验报告.doc

数字信号处理实验报告第二次实验IIR数字滤波器的设计姓名:陈桐 学号:04004316 实验日期:2006年11月14日一. 实验目的:(1) 掌握双线形变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线形变换法及脉冲响应不变法设计低通，高通和带通IIR数字滤波器的计算机编程。(2) 观察双线形变换法及脉冲响应不变法设计的滤波器的频域特性，了解双线形变换法及脉冲响应不变法的特点。(3) 熟悉巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器的频率特性。二. 实验原理:(1) 实验中有关变量的定义：fc通带边界频率，fr阻带边界频率，tao通带波动，at最小阻带衰减，fs采样频率，t采样周期。(2) 设计一个数字滤波器一般包括以下两步：a.按照任务要求，确定滤波器性能指标 b.用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求(3) 数字滤波器的实现：对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。用双线形变换法设计IIR数字滤波器的过程：a.将设计性能指标中的关键频率点进行“预畸”b.利用“预畸”得到的频率点设计一个模拟滤波器。c.双线形变换，确定系统函数三. 实验内容:(1) fc=0.3kHz,=0.8dB,fr=0.2kHz,At=20dB,T=1ms;设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。程序: clc;clear;figure(2);%预畸变模拟边界频率wc=2*1000*tan(2*pi*300/(2*1000);wr=2*1000*tan(2*pi*200/(2*1000);%阻带单调，通带波纹，方法为切比雪夫N,wn=cheb1ord(wc,wr,0.8,20,s);%高通B,A=cheby1(N,0.8,wn,high,s);num,den=bilinear(B,A,1000);%频谱h,w=freqz(num,den);f=w/pi*500;%画图plot(f,20*log10(abs(h);axis(0,500,-80,10);grid on;xlabel(频率/Hz);ylabel(幅度/dB);%图形为:由图易得，该滤波器完全满足要求。趋近于0，f=200Hz时，幅度约为-30dB，满足At=20dB的要求。(2) fc=0.2kHz, =1dB,fr=0.3kHz,At=25dB,T=1ms;分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。程序为: clc;clear;fs=1000; wp=2*pi*200;ws=2*pi*300n,wn=buttord(wp,ws,1,25,s);% 设计Butterworth低通模拟滤波器% 通带衰减小于1分贝、阻带（0.6pi）衰减大于25分贝 % n: 阶数% wn:低通模拟滤波器的3分贝通带截止频率B1,A1 = BUTTER(n,wn,s);% 设计Butterworth低通模拟滤波器% B,A：模拟滤波器传递函数分子、分母多项式系数b1,a1=impinvar(B1,A1,fs);% 用脉冲响应不变法将模拟低通滤波器转换为数字滤波器b,a：分子、分母多项式系数% 双曲变换的预畸变B2,A2=butter(n,2000*tan(pi/5),s);% 用双曲变换法将模拟低通滤波器转化为数字滤波器b,a:分子,分母多项式系数b2,a2=bilinear(B2,A2,fs); h1,w1=freqz(b1,a1,100,whole);% h: 在单位圆上采样100点的频率响应% w：频率矢量c1=20*log10(abs(h1)+0.0001); h2,w2=freqz(b2,a2,100,whole);% h: 在单位圆上采样100点的频率响应% w：频率矢量c2=20*log10(abs(h2)+0.0001); plot(w1(1:50), c1(1:50),-.);grid on;hold on;plot(w2(1:50), c2(1:50),-);grid on;%图形为:由图可以看出，利用脉冲响应不变法得到的滤波器在高频部分（阻带）衰减特性不好；而双线性变换法将模拟频率的无穷大映射到数字频率的pi,这样使得其对应的滤波器在高频部分衰减很快。(3) 利用双线性变换法分别设计满足下列指标的巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器，并作图验证设计结果：fc=1.2kHz, =40dB,fs=8kHz。程序: clc;clear;fs=8000; wp=2*pi*1200;ws=2*pi*2000n,wn=buttord(wp,ws,0.5,40,s);% 设计Butterworth低通模拟滤波器% 通带衰减小于1分贝、阻带（0.6pi）衰减大于25分贝 % n: 阶数% wn:低通模拟滤波器的3分贝通带截止频率B1,A1 = BUTTER(n,wn,s);% 设计Butterworth低通模拟滤波器% B,A：模拟滤波器传递函数分子、分母多项式系数b1,a1=impinvar(B1,A1,fs);% 用脉冲响应不变法将模拟低通滤波器转换为数字滤波器b,a：分子、分母多项式系数% 双曲变换的预畸变B2,A2=butter(n,16000*tan(pi*0.3),s);% 用双曲变换法将模拟低通滤波器转化为数字滤波器b,a:分子,分母多项式系数b2,a2=bilinear(B2,A2,fs);% 椭圆变换法,先对fs/2进行归一化wp=1200/(fs/2);ws=2000/(fs/2);n,wn=ellipord(wp,ws,0.5,40);% 椭圆参数数字化b3,a3 = ellip(n,0.5,40,wn);% 频率h1,w1=freqz(b1,a1,100,whole);% h: 在单位圆上采样100点的频率响应% w：频率矢量c1=20*log10(abs(h1)+0.0001); h2,w2=freqz(b2,a2,100,whole);% h: 在单位圆上采样100点的频率响应% w：频率矢量c2=20*log10(abs(h2)+0.0001);h3,w3=freqz(b3,a3,100,whole);% h: 在单位圆上采样100点的频率响应% w：频率矢量c3=20*log10(abs(h3)+0.0001);% 画图plot(w1(1:50), c1(1:50),-.);grid on;hold on;plot(w2(1:50), c2(1:50),-);grid on;hold on;plot(w3(1:50), c3(1:50),-);%图形为:由图可以看出，满足设计要求。(4) 分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字带通滤波器，已知fs=3kHz，其等效的模拟滤波器指标为3kHz, 2kHzf=5dB, f=6kHz, At=20dB,f=1.5kHz。程序: %参数定义fs=30000;tao=3;f1=2000;f2=3000;at1=5;f3=6000;f4=1500;at2=20;%预畸变w11=2*fs*tan(pi*f1/fs);w12=2*fs*tan(pi*f2/fs);w13=2*fs*tan(pi*f3/fs);w14=2*fs*tan(pi*f4/fs);%数字化n,wn=buttord(w11 w12,w14,w13,5,20,s);B,A=butter(n,wn,s);num,den=bilinear(B,A,fs);h,w=freqz(num,den);f=w/(pi*2)*fs;subplot(2,1,1);plot(f,20*log10(abs(h);grid;axis(0,15000,-100,0);title(双线性变换法);%脉冲响应的频率变化w1=2*pi*f1;w2=2*pi*f2;w3=2*pi*f3;w4=2*pi*f4;%数字化n,wn=buttord(w11 w12,w14,w13,5,20,s);B,A=butter(n,wn,s);num1,den1=impinvar(B,A,fs);h1,w1=freqz(num1,den1);f=w1/(pi*2)*fs;subplot(2,1,2),plot(f,20*log10(abs(h1);title(脉冲响应不变法);grid;axis(0,15000,-100,0);%图形为:由上图可以看出，用脉冲响应不变法由于滤波器的混叠作用在过度带和阻带都衰减的较双线性变换法慢。(5) 利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波器，并作图验证设计结果：当1kHz=f=18dB；当f=3kHz 时，=3dB ；采样频率fs=10kHz。程序: f3=1000;f4=2000;ar=18;tao=3;f1=500;f2=3000;fs=10000;%参数变化fp1=pi*f1/fs;fp2=pi*f2/fs;fc1=pi*f3/fs;fc2=pi*f4/fs;w1=2*fs*tan(fp1) 2*fs*tan(fp2);w2=2*fs*tan(fc1) 2*fs*tan(fc2);%数字化n,wn=cheb1ord(w1,w2,tao,ar,s);B,A=cheby1(n,tao,wn,stop,s);num,den=bilinear(B,A,fs);h,w=freqz(num,den);plot(w*fs/(2*pi),20*log10(abs(h);grid;%图形为:由上图可以看出，设计满足要求。四. 思考题：（1）双线性变换法中模拟频率和数字频率的关系是非线性的，在实验中你注意到这种非线性关系了吗？从哪几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系？答：观察到了。从（2）、（