实验四数字滤波器的设计实验报告.docx
- 文档编号:7708099
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:8
- 大小:21.57KB
实验四数字滤波器的设计实验报告.docx
《实验四数字滤波器的设计实验报告.docx》由会员分享,可在线阅读,更多相关《实验四数字滤波器的设计实验报告.docx(8页珍藏版)》请在冰豆网上搜索。
实验四数字滤波器的设计实验报告
数字信号处理
实验报告
实验四HR数字滤波器的设计
学
生姓
名
张志翔
班
级
电于信息工程1203班
学
号
指
导教
师
实验四
IIR数字滤波器的设计
一、实验目的:
1.掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设
计方法及其原理,熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的MATLAB编程。
2.观察双线性变换及脉冲响应不变法设计的滤波器的频域特性,了解双
线性变换法及脉冲响应不变法的特点。
3.熟悉Butterworth滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。
二、实验原理:
1.脉冲响应不变法
用数字滤波器的单位脉冲响应序列模仿模拟滤波器的冲激响应,让正好等于的采样值,即,其中为采样间隔,如果以及分别表示的拉式变换及的Z变换,则
2.双线性变换法
S平面与z平面之间满足以下映射关系:
S平面的虚轴单值地映射于Z平面的单位圆上,S平而的左半平而完全映射到Z平而的单位圆内。
双线性变换不存在混叠问题。
双线性变换是一种非线性变换,这种非线性引起的幅频特性畸变可通过预畸而得到校正。
三、实验内容及步2实验中有关变量的定义:
fc通带边界频率;fr阻带边界频率;6通带波动;At最小阻带
衰减;fs采样频率;T采样周期
(1)=0.3KHz,6=0.8Db,=0.2KHz,At=20Db,T=lms;
设计一个切比雪夫高通滤波器,观察其通带损耗和阻带衰减是否满足要求。
MATLAB源程序:
wp=2*1000*tan(2*pi*300/(2*1000));
ws=2*1000*tan(2*pi*200/(2*1000));
[N,wn]=cheblord(wp,ws,0.8,20,'s');%给定通带(wp)和阻带(ws)边界角频率,通带波动波动0.8,阻带最小衰减20dB,求出最低阶数和通带滤波器的通带边界频率Wn
[B,A]=chebyl(N,0.5,wn,'high',‘s');%给定通带(wp)和阻带(ws)边界角频率,通带波动
[num,den]=bilinear(B,A,1000);
[h,w]=freqz(num,den);
f二w/(2*pi)*1000;
plot(f,20*logl0(abs(h)));
axis([O,500,-80,10]);
grid;xlabel('频率');ylabel('幅度/dB')
程序结果
num=0.0304
-0.1218
0.1827-0.12180.0304
den=11
.38341.4721
0.80120.2286
系统函数:
0.0304-0.1218z_1+0.1827z-2-0.1218z-3+0.0304Z"4
Z)~1.0000+1.3834z_,+1.4721z-2+0.8012z-3+0.2286Z-4
幅频响应图:
分析:
由图可知,切比雪夫滤波器幅频响应是通带波纹,阻带单调衰减的。
6=0.8,fr=O.2kHz,At=30Db,满足设计要求
(2)fc=O.2kHz,6二ldB,fr=O.3kHz,At=25dB,T二Ims;分别用脉冲响应不变法及双线性变换法设计一Butterworth数字低通滤波器,观察所设计数字滤波器的幅频特性曲线,记录带宽和衰减量,检查是否满足要求。
比较这两种方法的优缺点。
MATLAB源程序:
T=0.001;fs=1000;fc=200;fr=300;
wpl=2*pi*fc;wrl=2*pi*ft;
[Nl,wnl]=buttord(wpl,wrl,1,25,'s')
[Bl,Al]=butter(Nl,wnl,'s');
[numl,deni]=impinvar(Bl,Al,fs);%脉冲响应不变法
[hl,w]=freqz(numl,deni);
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*logl0(abs(hl)),f,20*logl0(abs(h2)),'-;
axis([0,500,TOO,10]);grid;xlabel('频率/Hz');ylabel('幅度
/dBj
titleC巴特沃思数字低通滤波器');
legend脉冲相应不变法','双线性变换法',1);
结果分析:
脉冲响应不变法的低通滤波器系统函数:
numl-2.36470.00020.01530.09950.14440.0611
0.00750.00023.65690
deni1-1.91992.5324-2.20531.3868-0.63090.2045
-0.04500.0060-0.0004
双线性变换法设计的低通滤波器系统函数:
num20.01790.10720.26810.35750.2681
0.10720.0179
den21-0.60190.9130-0.29890.1501-0.0208
0.0025
分析:
脉冲响应不变法的频率变化是线性的,数字滤波器频谱响应岀现了混叠,影响了过渡带的衰减特性,并且无传输零点;双线性变化法的频率响应是非线性的,因而消除了频谱混叠,在f=500Hz出有一个传输零点。
脉冲响应不变法的一个重要特点是频率坐标的变换是线性的,3=QT,3与Q是线性关系:
在某些场合,要求数字滤波器在时域上能模仿模拟滤波器的功能时,如要实现时域冲激响应的模仿,一般使用脉冲响应不变法。
脉冲响应不变法的最大缺点:
有频谱周期延拓效应,因此只能用于带限的频响特性,如衰减特性很好的低通或带通,而高频衰减越大,频响的混淆效应越小,至于高通和带阻滤波器,由于它们在高频部分不衰减,因此将完全混淆在低频响应中,此时可增加一保护滤波器,滤掉高于的频带,再用脉冲响应不变法转换为数字滤波器,这会增加设计的复杂性和滤波器阶数,只有在一定要满足频率线性关系或保持网络瞬态响应时才米用。
双线性变换法的主要优点是S平面与Z平面一一单值对应,S平面的虚轴(整个jQ)对应于Z平面单位圆的一周,S平面的Q二0处对应于Z平而的3二0处,Q二8处对应于z平面的3二n处,即数字滤波器的频率响应终止于折叠频率处,所以双线性变换不存在混迭效应。
双线性变换缺点:
Q与3成非线性关系,导致:
a.数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变,(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。
b.线性相位模拟滤波器经双线性变换后,得到的数字滤波器为非线性相位。
c.要求模拟滤波器的幅频响应必须是分段恒定的,故双线性变换只能用于设计低通、高通、带通、带阻等选频滤波器。
(3)利用双线性变换法分别设计满足下列指标的Butterworth型、Chebyshev型和椭圆型数字低通滤波器,并作图验证设计结果:
fc=l.2kHz,§WO.5dB,fr二2kHz,At^40dB,fs=8kHz,比较这种滤波器的阶数。
MATLAB源程序:
clearall;
wc=2*pi*1200;wr二2*pi*2000;rp=O.5;rs=40;fs=8000;
wl=2*fs*tan(wc/(2*fs));
w2=2*fs*tan(wi7(2*fs));
[Nb,wn]二buttord(wl,w2,rp,rs,'s')%巴特沃思
[B,A]=butter(Nb,wn,'s');
[numl,denl]=bilinear(B,A,fs);
[hl,w]=freqz(numl,deni);
[Nc,wn]=cheblord(wl,w2,rp,rs,'s,)%切比雪夫
[B,A]=chebyl(Nc,rp,wn,'s');
[num2,den2]=bilinear(B,A,fs);
[h2,w]=freqz(num2,den2);
[Ne,wn]=ellipord(wl,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*logl0(abs(hl)),'-',f,20*logl0(abs(h2)f,20*log
10(abs(h3)),':
;
axis([0,4000,TOO,10]);grid;
xlabel(*FrequencyinHz');ylabelCGainindB');titleC三种数字低通滤波器');
legend(,巴特沃思数字低通滤波器',’切比雪夫数字低通滤波器',’椭圆数字低通滤波器',3);
巴特沃思数字低通滤波器的系统函数系数:
numl=
0.0032
0.01290.0302
0.0453
0.0453
0.0302
0.0129
0.0032
0.0003
deni二
-2.7996
4.4581—4.5412
3.2404
-1.6330
0.5780
-0.13700.0197-0.0013
切比雪夫数字低通滤波器的系统函数系数:
num2=0.00260.01320.02640.02640.01320.0026
den2=1-2.97754.2932-3.51241.6145-0.3334
椭圆数字低通滤波器的系统函数系数:
num3=0.038870.03630.06650.03630.0389
den3=1-2.14442.3658-1.32500.3332
程序结果图:
分析:
设计结果表明,巴特沃思数字低通滤波器、切比雪夫数字低通滤波器、椭圆数字低通滤波器的阶数分别是9、5、4阶。
可见,对于给定的阶数,椭圆数字低通滤波器的阶数最少(换言之,对于给定的阶数,过渡带最窄),就这一点来说,他是最优滤波器。
由图表明,巴特沃思数字低通滤波器过渡带最宽,幅频响应单调下降;椭圆数字低通滤波器过渡带最窄,并具有等波纹的通带和阻带响应;切比雪夫数字低通滤波器的过渡带介于两者之间。
(4)分别用脉冲响应不变法及双线性变换法设计一Butterworth型数字带通滤波器,己知,其等效的模拟滤波器指标6<3dB,2kHz MATLAB源程序: wpl=2*pi*2000;wp2=2*pi*3000; wsl=2*pi*1500;ws2=2*pi*6000; [Nl,wnlj=buttord([wplwp2],[wslws2],3,20,‘s');%求巴特沃思滤波器的阶数 [Bl,Al]=butter(Nl,wnl/s^;%给定阶数和边界频率设计滤波器 [numl,deni]=impinvar(Bl,Al,30000);%脉冲相应不变法 [hl,w]=freqz(numl,deni); wl二2*30000*tan(2*pi*2000/(2*30000)); w2=2*30000*tan(2*pi*3000/(2*30000)); wrl=2*30000*tan(2*pi*1500/(2*30000)); wr2=2*30000*tan(2*pi*6000/(2*30000)); [N,wn]=buttord([wlw2],[wrlwr2],3,20,'s');%求巴特沃思滤波器的阶数 [B,A]=butter(N,wn,'s'); [num,den]=bilinear(B,A,30000);%双线性变化法 [h2,w]=freqz(num,den); f二w/(2*pi)*30000; plot(f,20*logl0(abs(hl)),'-,f,20*logl0(abs(h2)),'-');axis([0,15000,-60,10]); xlabel(*FrequencyinHz');ylabel(,GainindB,); grid; titleC巴特沃思数字带通滤波器'); legendC脉冲相应不变法',’双线性变换法',1); 脉冲相应不变法设计的巴特沃思数字带通滤波器系统函数的分子、分母多项式系数: numl=-1.51580.0057-0.01220.00250.0089-0.0049 deni二1-4.805610.2376-12.26258.7012-3.47190.6145双线性变换法设计的巴特沃思数字带通滤波器系统函数的分子、分母多项式的系数: num=0.00140-0.00427.10540.00425.7732-0.0014 den=1-4.807110.2473-12.28388.7245-3.48490.6176 (5)利用双线性变换法设计满足下列指标的Chebyshev型数字带阻滤波器,并作图验证设计结果: 当lMzv/«2Mz时,At>\SdB;当/5500Hz以及f>3kHz时,553dB;采样频率fs=\OkHz0 MATLAB源程序: wl=2*10000*tan(2*pi*1000/(2*10000)); w2=2*10000*tan(2*pi*2000/(2*10000));wrl=2*10000*tan(2*pi*500/(2*10000)); wr2=2*10000*tan(2*pi*3000/(2*10000)); [N,wn]=cheblord([wrlwr2],[wlw2],3,l&'s');%计算阶数[B/A]=chebyl(N/3/wn;stop7s');%给定阶数和参数设计滤波器[numzden]=bilinear(B,A,10000);%双线性变化法 [h,w]=freqz(num,den);%频率响应 f=w/(2*pi)*1OOOO; plot(f,20*logl0(abs(h))); axis([0,5000,-120,10]); grid;xlabel(濒率/Hz');ylabel('幅度/dB') title/切比雪夫数字带阻滤波器匕 程序结果图: 四、实验思考题 1.双线性变换法中Q和3之间的关系是非线性的,在实验中你注意到这种非线性关系了吗从哪几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系 答: 在双线性变化法中,模拟频率与数字频率不再是线性关系,所以一个线性相位模拟滤波器经双线性变换后,得到的数字滤波器不再保持原有的线性相位了,在每一幅使用了双线性变换的图中,可以看到在采样频率一半处,幅度为零,这显然不是线性变换能够产生的,这是由于双线性变换将模拟域中的无穷远点映射到了改点处。 2.能否利用公式刃⑵「孑]完成脉冲响应不变法的数字滤波器设计为什么 答: HR数字滤波器的设计实际上是求解滤波器的系数线它是数学上的一种逼近问题,即在规定意义上(通常采用最小均方误差准则)去逼近系统的特性。 如果在S平而上去逼近,就得到模拟滤波器;如果在z平而上去逼近,就得到数字滤波器。 但是它的缺点是,存在频率混叠效应,故只适用于阻带的模拟滤波器。 五、实验总结 数字滤波器的设计是信号处理方而的重要内容,通过运用MATLAB软件来设计IIR数字滤波器,使我熟悉了MATLAB的强大功能,同时也对数字滤波器的特点、作用有更深入的理解。 了解了熟练的利用MATLAB这一功能强大的软件来设计数字滤波器对数字信号处理这一领域有着重要的意义与价。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 数字滤波器 设计 报告