资源描述：
《东南大学信息学院DSP第二次实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验三IIR数字滤波器的设计04013222张嘉俊一、实验目的（1）掌握双线形变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线形变换法及脉冲响应不变法设计低通，高通和带通IIR数字滤波器的计算机编程。（2）观察双线形变换法及脉冲响应不变法设计的滤波器的频域特性，了解双线形变换法及脉冲响应不变法的特点。（3）熟悉巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器的频率特性。（3）实验中有关变量的定义（4）设计一个数字滤波器一般包括以下两步a.按照任务要求，确定滤波器性能指标b.用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求（5）数字滤波器的
2、实现对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。用双线形变换法设计IIR数字滤波器的过程：a.将设计性能指标中的关键频率点进行“预畸”b.利用“预畸”得到的频率点设计一个模拟滤波器。c.双线形变换，确定系统函数二、实验内容（1）fc=0.3kHz，δ=0.8dB，fr=0.2kHz，At=20dB，T=1ms；设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。源程序：clc,clear;wc=2000*tan(2*pi*300*0.001/2);wr
3、=2000*tan(2*pi*200*0.001/2);[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/(2*pi)*1000;plot(f,20*log10(abs(h)));axis([0,500,-80,10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('ChebyshevHigh-passFilt
4、er');实验结果：分析与结论：观察图形可知，δ趋近于0；f=200Hz时，幅度约为-30dB，满足At=20dB的要求，故其通带损耗和阻带衰减满足要求。（2）fc=0.2kHz，δ=1dB，fr=0.3kHz，At=25dB，T=1ms；分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。clc,clear;wc=2*pi*200;wr=2*pi*300;wc2=2*1000*tan((2*pi)*200*0.001/2);wr2=2*1000*tan((2*p
5、i)*300*0.001/2);%脉冲响应不变法[N1,wn]=buttord(wc,wr,1,25,'s');[B1,A1]=butter(N1,wn,'s');[num1,den1]=impinvar(B1,A1,1000);[h1,w1]=freqz(num1,den1);%双线性变换法[N2,wn2]=buttord(wc2,wr2,1,25,'s');[B2,A2]=butter(N2,wn2,'s');[num2,den2]=bilinear(B2,A2,1000);[h2,w2]=freqz(num2,den2);f1=w1/pi*500;f2=w2/pi
6、*500;plot(f1,20*log10(abs(h1)),'-',f2,20*log10(abs(h2)),'--');legend('脉冲响应不变法','双线性变换法');axis([0,500,-80,10]);grid;xlabel('Frequency/Hz');ylabel('Amplitude/dB');title('ButterworthLow-passFilter');分析与结论：通带边界和阻带边界分别为200Hz,300Hz,最小阻带衰减量满足25dB。脉冲响应不变法优点：模拟频率到数字频率的转换是线性的，ω＝ΩΤ，时域特性逼近。缺点：有频谱周期延
7、拓效应，会产生频谱混叠现象，因此只适合带限滤波器双线性变换法优点：通过严重非线性关系将S平面与Z平面一一对应，不存在频谱混迭效应。缺点：时域到频域是非线性变换，数字滤波器的响应有畸变，在高频处有较大失真。（3）利用双线性变换法分别设计满足下列指标的巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器，并作图验证设计结果：fc=1.2kHz，δ<=0.5dB，fr=2kHz，At>=40dB，fs=8kHz。比较这三种滤波器的阶数。源程序：clc,clear;wr=2*8000*tan(2*pi*2000/(2*8000));wc=2*8000*ta
