实验四IIR数字滤波器的设计实验报告.docx
《实验四IIR数字滤波器的设计实验报告.docx》由会员分享，可在线阅读，更多相关《实验四IIR数字滤波器的设计实验报告.docx（11页珍藏版）》请在冰豆网上搜索。
实验四IIR数字滤波器的设计实验报告
数字信号处理
实验报告
实验四IIR数字滤波器的设计
学
生
姓
名
张志翔
班
级
电子信息工程1203班
学
号
12401720522
指
导
教
师
2015429
实验四IIR数字滤波器的设计
一、实验目的：
1.掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体
设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低
通、高通和带通IIR数字滤波器的MATLA编程。
2.观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了
解双线性变换法及脉冲响应不变法的特点。
3.熟悉Butterworth滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。
二、实验原理：
1.脉冲响应不变法用数字滤波器的单位脉冲响应序列模仿模拟滤波器的冲激响应，让正好等于的采样值，即，其中为采样间隔，如果以及分别表示的拉式变换及的Z变换，则
12
H（z）zesT〒Ha（sj—m）
ImI
2.双线性变换法
S平面与z平面之间满足以下映射关系:
1Is
1Ts
s平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完
全映射到z平面的单位圆内双线性变换不存在混叠问题。
双线性变换是一种非线性变换，这种非线性引起的幅频特性畸变
可通过预畸而得到校正。
三、实验内容及步骤：
实验中有关变量的定义：
fc通带边界频率；fr阻带边界频率；3通带波动；At最小阻
带衰减；fs采样频率；T采样周期
（1）=0.3KHz,S=0.8Db,=0.2KHz,At=20Db,T=1ms;
设计一个切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。
MATLA源程序：
wp=2*1000*tan（2*pi*300/（2*1000））;
ws=2*1000*tan（2*pi*200/（2*1000））;
[N,wn]=cheb1ord（wp,ws,0.8,20,'s'）;%给定通带（wp）和阻带（ws）
边界角频率，通带波动波动0.8，阻带最小衰减20dB,求出最低阶数和通带滤波器的通带边界频率Wn
[B,A]=cheby1（N,0.5,wn,'high','s'）;%给定通带（wp）和阻带（ws）边
界角频率，通带波动
[num,den]=bilinear（BA1000）;
[h,w]=freqz（num,den）;
f=w/（2*pi）*1000;
plot（f,20*log10（abs（h）））;
axis（[0,500,-80,10]）;
grid;xlabel（'频率'）;ylabel（'幅度/dB'）
程序结果
num=0.0304-0.1218
0.1827-0.1218
0.0304
den=11.38341.47210.80120.2286
=0.0304-0.1218Z10.1827z2-0.1218z30.0304z4
系统函数：
=1.0000+1.3834z1+1.4721z2+0.8012z3+0.2286z4
幅频响应图:
Figure1
FileEditViewInsertToolsDesktopWindowHelp
□曰£空|脅®熄貳TD|口图|・口
m蚩■£
分析：
由图可知，切比雪夫滤波器幅频响应是通带波纹，阻带单调衰
减的。
3=0.8，fr=0.2kHz,At=30Db,满足设计要求
（2）fc=0.2kHz,S=1dB,fr=0.3kHz,At=25dB,T=1ms;分别用脉冲响应不变法及双线性变换法设计一Butterworth数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。
比较这两种方法的优缺点。
MATLA源程序：
T=0.001;fs=1000;fc=200;fr=300;
wp1=2*pi*fc;wr1=2*pi*fr;
[N1,wn1]=buttord（wp1,wr1,1,25,'s'）
[B1,A1]=butter（N1,wn1,'s'）;
[num1,den1]=impinvar（B1,A1,fs）;%脉冲响应不变法
[h1,w]=freqz（num1,den1）;
wp2=2*fs*tan（2*pi*fc/（2*fs））
wr2=2*fs*tan（2*pi*fr/（2*fs））
[N2,wn2]=buttord（wp2,wr2,1,25,'s'）
[B2,A2]=butter（N2,wn2,'s'）;
[num2,den2]=bilinear（B2,A2,fs）;%双线性变换法
[h2,w]=freqz（num2,den2）;
f=w/（2*pi）*fs;
plot（f,20*log10（abs（h1））,'-.',f,20*log10（abs（h2））,'-'）;
axis（［0,500,-100,10］）;grid;xlabel（'频率/Hz'）;ylabel（'幅度
/dB'）
title（'巴特沃思数字低通滤波器'）;
legend（'脉冲相应不变法','双线性变换法',1）;
结果分析：
脉冲响应不变法的低通滤波器系统函数：
num1-2.36470.00020.01530.09950.14440.0611
0.00750.00023.65690
1.3868-0.6309
den11-1.91992.5324-2.2053
0.2045-0.04500.0060-0.0004
2.36470.0002Z0.0153z0.0995z0.1444z0.061传0.0075z0.0002z3.6569z
LJ（
11.9199z2.5324z2.2053z1.3869z0.6309z0.2045z0.0450z0.0060z0.0004z
双线性变换法设计的低通滤波器系统函数
0.0025
0.01790.107N
分析：
脉冲响应不变法的频率变化是线性的，数字滤波器频谱响应出现了混叠，影响了过渡带的衰减特性，并且无传输零点；双线性变化法的频率响应是非线性的，因而消除了频谱混叠，在f=500Hz出有一个传输零点。
脉冲响应不变法的一个重要特点是频率坐标的变换是线性的，w=QT,®与Q是线性关系：
在某些场合，要求数字滤波器在时域上能模仿模拟滤波器的功能时，如要实现时域冲激响应的模仿，一般使用脉冲响应不变法。
脉冲响应不变法的最大缺点：
有频谱周期延拓效应，因此只能用于带限的频响特性，如衰减特性很好的低通或带通，而高频衰减越大，频响的混淆效应越小，至于高通和带阻滤波器，由于它们在高频部分不衰减，因此将完全混淆在低频响应中，此时可增加一保护滤波器，滤
掉高于的频带，再用脉冲响应不变法转换为数字滤波器，这会增
加设计的复杂性和滤波器阶数，只有在一定要满足频率线性关系或保持网络瞬态响应时才采用。
双线性变换法的主要优点是S平面与Z平面一一单值对应，s平面的虚轴（整个jQ）对应于Z平面单位圆的一周，S平面的Q=0处对应于Z平面的3=0处，Q二乂处对应于Z平面的3=n处，即数字滤波器的频率响应终止于折叠频率处，所以双线性变换不存在混迭效应。
双线性变换缺点：
Q与3成非线性关系，导致：
a.数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，
（使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变）。
b.线性相位模拟滤波器经双线性变换后，得到的数字滤波器为非线性相位。
c.要求模拟滤波器的幅频响应必须是分段恒定的，故双线性变换只
能用于设计低通、高通、带通、带阻等选频滤波器。
（3）利用双线性变换法分别设计满足下列指标的Butterworth型、
Chebyshev型和椭圆型数字低通滤波器，并作图验证设计结果：
fc=1.2kHz,0.5dB,fr=2kHz,At>40dB,fs=8kHz，比较这
种滤波器的阶数。
MATLA源程序：
clearall;
wc=2*pi*1200;wr=2*pi*2000;rp=0.5;rs=40;fs=8000;w1二2*fs*tan（wc/（2*fs））;
w2=2*fs*tan（wr/（2*fs））;
[Nb,wn]二buttord（w1,w2,rp,rs,'s'）%巴特沃思
[B,A]=butter（Nb,wn,'s'）;
[num1,den1]=bilinear（B,A,fs）;
[h1,w]=freqz（num1,den1）;
[Nc,wn]=cheb1ord（w1,w2,rp,rs,'s'）%切比雪夫
[B,A]=cheby1（Nc,rp,wn,'s'）;
[num2,den2]=bilinear（B,A,fs）;
[h2,w]=freqz（num2,den2）;
[Ne,wn]=ellipord（w1,w2,rp,rs,'s'）%椭圆型
[B,A]=ellip（Ne,rp,rs,wn,'low','s'）;
[num3,den3]=bilinear（B,A,fs）;
[h3,w]=freqz（num3,den3）;
f=w/（2*pi）*fs;
plot（f,20*log10（abs（h1））,'-',f,20*log10（abs（h2））,'--
：
f,20*log10（abs（h3））,':
'）;
axis（[0,4000,-100,10]）;grid;
xlabel（'FrequencyinHz'）;ylabel（'GainindB'）;
title（'三种数字低通滤波器'）;
legend（'巴特沃思数字低通滤波器','切比雪夫数字低通滤波器椭圆数字低通滤波器',3）;
巴特沃思数字低通滤波器的系统函数系数:
-0.13700.0197-0.0013
切比雪夫数字低通滤波器的系统函数系数:
椭圆数字低通滤波器的系统函数系数:
num3=0.038870.03630.06650.03630.0389
den3二1-2.14442.3658-1.32500.3332
程序结果图:
三种数字低通滤液器
10
050010001500200Q2500300035004000
FrequencyinHz
分析：
设计结果表明，巴特沃思数字低通滤波器、切比雪夫数字低通滤波器、椭圆数字低通滤波器的阶数分别是9、5、4阶。
可见，对于给定的阶数，椭圆数字低通滤波器的阶数最少（换言之，对于给定的阶数，过渡带最窄），就这一点来说，他是最优滤波器。
由图表明，巴特沃思数字低通滤波器过渡带最宽，幅频响应单调下降；椭圆数字低通滤波器过渡带最窄，并具有等波纹的通带和阻带响应；切比雪夫数字低通滤波器的过渡带介于两者之间。
（4）分别用脉冲响应不变法及双线性变换法设计一Butterworth型数字带通滤波器，已知，其等效的模拟滤波器指标S<3dB,2kHzvfw3kHz;At>5dB,f>6kHz;At>20dB,f<1.5kHz。
MATLA源程序：
wp1=2*pi*2000;wp2=
升级会员 