资源描述：
《实验四 IIR数字滤波器的设计 数字信号处理 DSP》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第二次实验一.实验目的:(1)掌握双线形变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线形变换法及脉冲响应不变法设计低通，高通和带通IIR数字滤波器的计算机编程。(2)观察双线形变换法及脉冲响应不变法设计的滤波器的频域特性，了解双线形变换法及脉冲响应不变法的特点。(3)熟悉巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器的频率特性。二.实验原理:(1)实验中有关变量的定义：fc通带边界频率，fr阻带边界频率，tao通带波动，at最小阻带衰减，fs采样频率，t采样周期。(2)设计一个数字滤波器一般包括以
2、下两步：a.按照任务要求，确定滤波器性能指标b.用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求(3)数字滤波器的实现：对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。用双线形变换法设计IIR数字滤波器的过程：a.将设计性能指标中的关键频率点进行“预畸”b.利用“预畸”得到的频率点设计一个模拟滤波器。c.双线形变换，确定系统函数三.实验内容:(1)fc=0.3kHz,δ=0.8dB,fr=0.2kHz,At=2
3、0dB,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
4、(f,20*log10(abs(h)));axis([0,500,-80,10]);gridon;xlabel('频率/Hz');7ylabel('幅度/dB');图形为:由图易得，该滤波器完全满足要求。δ趋近于0，f=200Hz时，幅度约为-30dB，满足At=20dB的要求。(1)fc=0.2kHz,δ=1dB,fr=0.3kHz,At=25dB,T=1ms;分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。程
5、序为:clc;clear;fs=1000;wp=2*pi*200;ws=2*pi*300[n,wn]=buttord(wp,ws,1,25,'s');[B1,A1]=butter(n,wn,'s');[b1,a1]=impinvar(B1,A1,fs);[B2,A2]=butter(n,2000*tan(pi/5),'s');[b2,a2]=bilinear(B2,A2,fs);[h1,w1]=freqz(b1,a1,100,'whole');c1=20*log10(abs(h1)+0.0001);[h2,w2]=freq
6、z(b2,a2,100,'whole');c2=20*log10(abs(h2)+0.0001);plot(w1(1:50),c1(1:50),'-.');gridon;holdon;plot(w2(1:50),c2(1:50),'-');7gridon;图形为:图中实线为双曲不变法，虚线为脉冲响应不变法由图可以看出，利用脉冲响应不变法得到的滤波器在高频部分（阻带）衰减特性不好；而双线性变换法将模拟频率的无穷大映射到数字频率的pi,这样使得其对应的滤波器在高频部分衰减很快。(1)利用双线性变换法分别设计满足下列指标的巴特沃
7、思滤波器，切比雪夫滤波器和椭圆滤波器，并作图验证设计结果：fc=1.2kHz,δ<=0.5dB,fr=2kHz,At>=40dB,fs=8kHz。程序:clc;clear;fs=8000;wp=2*pi*1200;ws=2*pi*2000[n,wn]=buttord(wp,ws,0.5,40,'s');[B1,A1]=butter(n,wn,'s');[b1,a1]=impinvar(B1,A1,fs);[B2,A2]=butter(n,16000*tan(pi*0.3),'s');[b2,a2]=bilinear(B2,
8、A2,fs);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');c1=20*log10(abs(h1)+0.000
