数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx

上传人:王** 文档编号:1040242 上传时间:2024-03-16 格式:DOCX 页数:7 大小:225.04KB
下载 相关 举报
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第1页
第1页 / 共7页
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第2页
第2页 / 共7页
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第3页
第3页 / 共7页
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第4页
第4页 / 共7页
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第5页
第5页 / 共7页
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第6页
第6页 / 共7页
数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx_第7页
第7页 / 共7页
亲,该文档总共7页,全部预览完了,如果喜欢就下载吧!
资源描述

《数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx(7页珍藏版)》请在优知文库上搜索。

1、实验四用窗函数法设计FIR数字滤波器一、实验目的(1)掌握用窗函数法设计FlR数字滤波器的原理和方法。(2)熟悉线性相位FIR数字滤波器特性。(3)了解各种窗函数对滤波特性的影响。二、实验原理、滤波器的理想频率响应函数为Hd(),那么其对应的单位脉冲响应为hd(n)二白瓦32咽窗函数设计法的根本原理是用有限长单位脉冲响应序列h(n)逼hd(n)t由于hd(n)往往是无限长序列,且是非因果的,所以用窗函数。w(n)将hd(n)截断,并进行加权处理:h(n)=hcj(n)w(n)h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数HG)为N-IH(ej)=Yjh(n)ejn=

2、G用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择适宜的窗函数类型和窗口长度N。三、实验内容1 .MATALAB程序任务一Nl=15;N2=33;bl=firl(14,l4,hanning(15);b2=fir1(32,1/4,hanning(33);Hl,W=freqz(bl,l);HLdb=20*1Og10(abs(Hl);magHl=abs(Hl);phaHl=angle(Hl);H2,W=freqz(b2,l);H2.db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(

3、H2);figure(l);subplot(2,1,1);stem(b1);title(N=15时,汉宁窗得到的FIR滤波器脉冲响应subplot(2,1,2);stem(b2);title(N=33时,汉宁窗得到的FIR滤波器脉冲响应);figure(2);subplot(3,1,1);plot(W/pi,magH1);title(N=15时汉宁窗得到的FIR滤波器的幅频特性)subplot(3,1,2);plot(W/pi,Hl_db);title(N=15时汉宁窗得到的FIR滤波器的对数幅频特性)SUbPIOt(3,3);PIot(Wpi,phaH1);title(N=15时汉宁窗得到的

4、FIR滤波器的相频特性)figure(3);subplot(3,1,1);plot(W/pi,magH2);title(N=33时汉宁窗得到的FIR滤波器的幅频特性)subplot(3,1,2);plot(W/pi,H2_db);title(N=33时汉宁窗得到的FIR滤波器的对数幅频特性)subplot(3,1,3);plot(Wpi,phaH2);title(N=33时汉宁窗得到的FIR滤波器的相频特性)任务二N=33;b1=firl(32,l4,boxcar(33);b2=firl(32,l4,hanning(33);b3=fir1(32,1/4,bartlett(33);b4=firl

5、(32,l4,blackman(33);Hl,W=freqz(bl,l);HLdb=20*1Og10(abs(H1);magHl=abs(Hl);phaHl=angle(Hl);H2,W=freqz(b2,l);H2_db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(H2);H3,W=freqz(b3,l);H3_db=20*log10(abs(H3);magH3=abs(H3);phaH3=angle(H3);H4,W=freqz(b4,l);H4_db=20*log10(abs(H4);magH4=abs(H4);phaH4=angle(H4);f

6、igure(l);subplot(2,2,1);stem(b1);titleC矩形窗得到的FIR滤波器脉冲响应)subplot(2,2,2);stem(b2);title,汉宁窗得到的FIR滤波器脉冲响应)SUbPIOt(2,2,3);Stem(b3);title。三角窗得到的FIR滤波器脉冲响应)SUbPlOt(2,2,4);Stem(b4);title,布拉克曼窗得到的FIR滤波器脉冲响应)figure(2);subplot(2,2,1);plot(W/pi,magH1);titleC矩形窗得到的FIR滤波器幅频特性)SUbPIOt(2,2,2);PIOt(Wpi,magH2);title

7、,汉宁窗得到的FIR滤波器幅频特性)SUbPlOt(2,2,3);PlOt(Wpi,magH3);titleC三角窗得到的FIR滤波器幅频特性)SUbPIOt(2,2,4);PIOt(Wpi,magH4);title,布拉克曼窗得到的FIR滤波器幅频特性)figure(3);subplot(2,2,1);plot(W/pi,Hl_db);出世矩形窗得到的FIR滤菠器相频特性)subplot(2,2,2);plot(W/pi,H2_db);title,汉宁得到的FIR滤波器相频特性)subplot(2,2,3);plot(W/pi,H3_db);出世三角窗得到的FIR滤菠器相频特性)SUbPIO

8、t(2,2,4);Plot(W/pi,H4_db);title,布拉克曼得到的FIR禧波器相频特性)figure(4);subplot(2,2,1);plot(Wpi,phaH1);title,矩形窗得到的FIR滤波器相频特性)SUbPIOt(2,2,2);PIOt(Wpi,phaH2);titleC汉宁窗得到的FIR滤波器相频特性)SUbPIOt(2,2,3);PIOt(Wpi,phaH3);titleC三角窗得到的FIR滤波器相频特性)SUbPIOt(2,2,4);PlOt(Wpi,phaH4);title,布拉克曼窗得到的FIR滤波器相频特性)2 .实验波形图任务一0.30.250.20

9、.150.10050N=15时,汉宇苜得到的Fl磁波器肽冲响应0250.2015N=33时双宁窗得到的Fl磁波器脉冲响应0.33035N=15时汉宁亩得到的FIR滤波器的幅频特性任务二布拽克曼由得到的Fl礴源器脉冲雨应0.30250.20.150.10.05四、实验结论1.N的大小决定了窗谱的主瓣宽度,N越大,窗谱的主瓣宽度越大2 .最小阻带衰减只有窗行决定,不受N的影响,过渡带宽度与N和窗形都有关,N越大,过渡带宽越小3 .由实验可知滤波特性:布拉克曼窗汉宁窗三角窗矩形窗五、思考题(1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?写出设计步骤。答:根据A二s一%计算出所要设计的滤波器的过渡带查表找到满足的阻带最小衰减AS所需的窗函数利用已选的窗函数的过渡带宽除以AG得到低通滤波器单位脉冲响应的长度N根据/)二吊()求出h(n)Af-I再根据H(ei)=2%求出响应函数H(ei)n=0如果要求用窗函数法设计带通滤波器,且给定上、下边带截止频率为必和0)2,试求理想带通的单位脉冲响应hd(n)o答:由=竺言算出纭,求出滤波器过渡带宽Ag=q-牡和滤波器长度N,最后由hd(n)=7-_cSinl.(一理声c或()=幺,=7求出儿(),其中丁=-n-)2

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 高等教育 > 工学

copyright@ 2008-2023 yzwku网站版权所有

经营许可证编号:宁ICP备2022001189号-2

本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!