![FIR数字滤波器设计与使用_第1页](http://file3.renrendoc.com/fileroot_temp3/2022-5/5/00a04f34-cebc-46db-b4f5-18f655b20878/00a04f34-cebc-46db-b4f5-18f655b208781.gif)
![FIR数字滤波器设计与使用_第2页](http://file3.renrendoc.com/fileroot_temp3/2022-5/5/00a04f34-cebc-46db-b4f5-18f655b20878/00a04f34-cebc-46db-b4f5-18f655b208782.gif)
![FIR数字滤波器设计与使用_第3页](http://file3.renrendoc.com/fileroot_temp3/2022-5/5/00a04f34-cebc-46db-b4f5-18f655b20878/00a04f34-cebc-46db-b4f5-18f655b208783.gif)
![FIR数字滤波器设计与使用_第4页](http://file3.renrendoc.com/fileroot_temp3/2022-5/5/00a04f34-cebc-46db-b4f5-18f655b20878/00a04f34-cebc-46db-b4f5-18f655b208784.gif)
![FIR数字滤波器设计与使用_第5页](http://file3.renrendoc.com/fileroot_temp3/2022-5/5/00a04f34-cebc-46db-b4f5-18f655b20878/00a04f34-cebc-46db-b4f5-18f655b208785.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、刖户a货实验报告课程名称:数字信号处理指导老师:刘英成绩:实验名称:FIR数字滤波器设计与使用同组学生姓名:一、实验目的和要求设计和应用FIR低通滤波器。掌握FIR数字滤波器的窗函数设计法,了解设计参数(窗型、窗长)的影响。二、实验内容和步骤编写MATLAB程序,完成以下工作。2-1设计两个FIR低通滤波器,截止频率6c=0.5几(1)用矩形窗,窗长N=41o得出第一个滤波器的单位抽样响应序列hi(n)0记下hi(n)的各个抽样值,显示h1(n)的图形(用stem(.)。求出该滤波器的频率响应(的N个抽样)H1(k),显示|H1(k)|的图形(用plot(.)。(2)用汉明窗,窗长N=41o得
2、出第二个滤波器的单位抽样响应序列h2(n)0记下h2(n)的各个抽样值,显示h2(n)的图形。求出滤波器的频率响应H2(k),显示|H2(k)|的图形。(3)由图形,比较h1(n)与h2(n)的差异,|H1(k)|与|H2(k)|的差异。2-2产生长度为200点、均值为零的随机信号序列x(n)(用rand(1,200)-0.5)。显示x(n)。求出并显示其幅度谱|X(k)|,观察特征。2-3滤波(1)将x(n)作为输入,经过第一个滤波器后的输出序列记为y1(n),其幅度谱记为|Y1(k)|。显示|X(k)|与|Y1(k)|,讨论滤波前后信号的频谱特征。(2)将x(n)作为输入,经过第二个滤波器
3、后的输出序列记为y2(n),其幅度谱记为|Y2(k)|。比较|Y1(k)|与|Y2(k)|的图形,讨论不同的窗函数设计出的滤波器的滤波效果。2-4设计第三个FIR低通滤波器,截止频率切C=0.5几用矩形窗,窗长N=127O用它对x(n)进行滤波。显示输出信号y3(n)的幅度谱|Y3(k)|,并与|Y1(k)|比较,讨论不同的窗长设计出的滤波器的滤波效果。三、主要仪器设备自行编程。四、操作方法和实验步骤(参见“二、实验内容和步骤”)五、实验数据记录和处理5.1列出MATLAB程序清单,加注释。%2-1clear;closeall;clc%thelengthofwindow=41wc=0.5*pi
4、;N=41;n=0:N-1;h1=fir1(N-1,wc/pi,boxcar(N);H1,w1=freqz(h1);h2=fir1(N-1,wc/pi);H2,w2=freqz(h2);figure;subplot(1,2,1);stem(n,h1,'f川ed');gridon;xlabel('n','FontSize',14);ylabel('h1(n)','FontSize',14);title('矩形窗得到的低通滤波器h1(n)的时域','FontSize',14);subpl
5、ot(1,2,2);stem(n,h2,'filled');gridon;axis(030-0.20.6);xlabel('n','FontSize',14);ylabel(title('矩形窗得到的低通滤波器h2(n)figure;subplot(1,2,1);plot(w1/pi,abs(H1);gridon;xlabel('omega/pi','FontSize'title('矩形窗得到的低通滤波器h1(n)subplot(1,2,2);plot(w2/pi,abs(H2);gridon;xl
6、abel('omega/pi','FontSize'title('矩形窗得到的低通滤波器h2(n)figure;subplot(1,2,1);'h2(n)','FontSize',14);的时域,FontSize',14);,14);ylabel('|H1(k)|','FontSize'的幅度谱','FontSize',14);,14);ylabel('|H2(k)|','FontSize'的幅度谱','Font
7、Size',14);,14);,14);plot(w1/pi,20*log10(abs(H1);gridon;axis(01-10020);/pi','FontSize',14);ylabel('20lg|H1(k)|/dB','FontSize',14);title('h1(n)的幅度谱,FontSize',14);subplot(1,2,2);plot(w2/pi,20*log10(abs(H2);gridon;xlabel('omega/pi','FontSize',14);y
8、label('20lg|H2(k)|/dB','FontSize',14);title('h2(n)的幅度谱,FontSize',14);xlabel('omegapause;%2-2N=200;n=0:N-1;x=rand(1,N)-0.5;X,w=freqz(x);figure;subplot(2,1,1);stem(n,x,filled');gridon;xlabel('n','FontSize',14);ylabel('x(n)','FontSize'titl
9、e('随机序列x(n)的时域,FontSize',14);subplot(2,1,2);plot(w/pi,abs(X);gridon;xlabel('omega/pi','FontSize',14);ylabel('|X(k)|'title('随机序列x(n)的幅度谱','FontSize',14);pause;,14);,'FontSize',14);%2-3y1=filter(h1,1,x);Y1,w1=freqz(y1);y2=filter(h2,1,x);Y2,w2=fre
10、qz(y2);figure;subplot(2,1,1)plot(w1/pi,abs(Y1);gridon;axis(01012);set(gca,xlabel('omega/pi','FontSize',14);ylabel(title('输出信号y1(n)的幅度谱,FontSize',14);subplot(2,1,2);plot(w2/pi,abs(Y2);gridon;axis(01012);set(gca,xlabel('omega/pi','FontSize',14);ylabel(title('
11、;输出信号y2(n)的幅度谱,FontSize',14);'ytick',0:4:12);'|Y1(k)|','FontSize',14);'ytick',0:4:12);'|Y2(k)|','FontSize',14);pause;%2-4N=127;n=0:N-1;h3=fir1(N-1,wc/pi,boxcar(N);H3,w3=freqz(h3);y3=filter(h3,1,x);丫3,wy3=freqz(y3);figure;plot(wy3/pi,abs(Y3);gridon;
12、axis(01012);set(gca,'ytick',0:4:12);xlabel('omega/pi','FontSize',14);ylabel('|Y3(k)|','FontSize',14);title('输出信号y3(n)的幅度谱,FontSize',14);figure;subplot(1,2,1);plot(w3/pi,abs(H3);gridon;xlabel('omega/pi','FontSize',14);ylabel('|H3(k)|
13、','FontSize',14);title('矩形窗得到的低通滤波器h3(n)的幅度谱,FontSize',14);subplot(1,2,2);plot(w3/pi,20*log10(abs(H3);gridon;axis(01-10020);xlabel('omega/pi','FontSize',14);ylabel('20lg|H3(k)|/dB','FontSize',14);title('h3(n)的幅度谱,FontSize',14);5.2列出计算结果,包括h
14、1(n)和h2(n)的各个抽样值,|H1(k)卜|H2(k)|、|X(k)|、|Y1(k)|、|Y2(k)|和|Y3(k)|的图形。(1) h1(n)和h2(n)世明宙百乳逋注衣器同组M城百刎的旗器M的时现-n?(2) |H1(k)卜|H2(k)|语用审号上的位如更稗履H1时的甯注,14109口。国明智将到的跟通淖疵糖的幅度港1J口30字a1DOHDZ0.3M05D.B£>Znh北加的幅度港0n.1D203O.d.0.5D607O.B39-ahi30h2,巾的徜叟借(3)|X(k)|、|Yi(k)卜|Y2(k)|和|Y3(k)|051Dmwsrlx晋N10D-DodajB.BC
15、M口号a.E:/0605也反B6IL序制川川的时域40140|£0tti100nSD机序列即的幅度璃BD10203DI063B'田氐范出.:强E鼠号春君六、实验结果与分析观察结果,进行讨论,最后总结:滤波器的频率响应中的过渡带宽度取决于哪些设计参数?有什么规律?阻带最小衰减取决于哪些设计参数?有什么规律?(1)观察实验结果,逐一进行讨论拒悔窗得到的底通滤被嚣M附的时战演用多四氢时代通声旅图h2旧药M医0.50<QL1-Q.I-I必13店-02E1口恰为X期招的结论:从上面的序列图和记录的h1(n)与h2(n)的序列值中,我们可以看出h1(n)与h2(n)的最大值都相同,
16、在主瓣上取样点数相同且取样值几乎相等,但旁瓣上h1(n)的取样值在0的周围波动较大,h2(n)的取样值在0的周围波动较小。W典皆若到的跟通鹿液幅的旧境谱从|H1(k)|与|H2(k)|的图中(即h1(n)与h2(n)的幅度谱)的对比中,我们可以看出两者在w=0时,值都为1,ftw=wc=0.5W,值都为1/2H(0)=0.5,在g处,值都为0,说明两者都是满足设计要求的低通滤波器。但在通带内,|H1(k)|的波动较大,而|H2(k)|几乎没有波动;但|H1(k)|的过渡带宽度明显小于|H2(k)|的过渡带宽度。h割地和睡里迳30用BP/IWE工Fsri从dB发减图的对比上,我们可以看出,|H2
17、(k)|的最小我减为-53dB左右,而|H1(k)|的最小衰减为-20dB左右,|H2(k)|的最小衰减明显更小。(2)观察x(n)和|X(k)|第机的时域从上图我们可以看出随机序列x(n)既存在高频分量也存在低频分量(3)不同窗型对于滤波效果的比较从上图我们可以看出经过滤波之后,大于截止频率wc=0.5冗的高频部分几乎为0,保留了小于截止频率wc=0.5冗的低频部分。从上图|X(k)|、|Y1(k)|、|Y2(k)|的对比,我们可以看出用矩形窗和汉明窗设计的低通滤波器都过滤掉了大于截止频率wc=0.5冗的高频部分,保留了小于截止频率wc=0.5冗的低频部分;但矩形窗设计的低通滤波器在高频部分
18、还存在微小的量,而汉明窗设计的低通滤波器在高频部分则几乎为0,高频部分滤波效果更好;但在刚大于截止频率wc=0.5冗的部分通过矩形窗设计的低通滤波器输出信号很快衰减为0,而通过汉明窗设计的低通滤波器输出信号则还存在部分较大的分量;可以证明通过矩形窗设计的低通滤波器的。(4)窗长的影响MS)的福比¥ham的斑iii心的的礴星通从上图可以看出,窗长为127设计出的低通滤波器通带内波动起伏变密,但相对振荡幅度却几乎不改变,证明了吉布斯效应的存在;同时窗长为127设计出的低通滤波器的过渡带更窄更陡,但窗长为127设计出的低通滤波器和窗长为31设计出的低通滤波器的阻带最小衰减都为-20dB左右,说明滤波器的阻带最小衰减与窗长无关。总结:(1)从实验中可以看出:滤波器的频率响
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2023年医用制氧机(系统)项目综合评估报告
- 【正版授权】 IEC 62443-2-4:2015+AMD1:2017 CSV EN Security for industrial automation and control systems - Part 2-4: Security program requirements for IACS service providers
- 2024年校园低碳生活演讲稿
- 2024年创业交流会演讲稿
- 联醇技术问答
- 2024年医院安全生产月活动总结(5篇)
- 2025高考物理复习探究两个互成角度的力的合成规律课件教案练习题
- 2024体育教研组个人年终工作总结(34篇)
- 临床急诊值班阅片技巧
- 电气控制与PLC应用技术(FX5U) 课件全套 贺建华 模块1-6 继电控制电路装调 - PLC组网通信和模拟量输入输出的应用
- 云南省民办高校教师薪酬激励机制研究
- 学习宣传贯彻中国工会十八大精神试题附答案
- 《邮轮乘务员职业道德与素养》教学课件-01职业道德
- 辽宁省大连市甘井子区2021-2022学年五年级下学期期末科学试卷
- DINEN1706铝和铝合金铸件化学成分和机械性能(中文版)
- 游乐园安全责任书
- 4升5奥数思维特训:三角形综合试题-数学五年级上册人教版
- 智能电能表功能课件
- 2022-2023学年广东省深圳市龙岗区四年级数学第二学期期末综合测试试题含答案
- 外国民商案例选读-东北大学中国大学mooc课后章节答案期末考试题库2023年
- 跨京开高速公路中承式钢箱混凝土拱桥监控监测方案
评论
0/150
提交评论