《数字信号处理_第九题带阻_》由会员分享，可在线阅读，更多相关《数字信号处理_第九题带阻_（12页珍藏版）》请在金锄头文库上搜索。
1、课程设计报告(2013- 2014年度第一学期)名 称: 信号的分析与处理 题 目：带阻FIR数字滤波器设计 院 系 动力工程系 班 级： 测控11K2 学 号： 111912020229 学生姓名： 杨睿 指导教师： 白康 设计周数： 一周 成 绩：日期：2014年01月12日数字信号处理课程设计任 务 书一、 目的与要求能够运用本课程中学到的知识，设计基于窗口函数法的FIR 数字滤波器。要求掌握数字信号处理的基本方法；FIR滤波器的设计步骤和方法；能够熟练采用C 语言或MATLAB语言进行计算机辅助设计和仿真验证设计内容的合理性。二、 主要内容1） 掌握设计数字滤波器的基本步骤；2） 重点
2、掌握利用窗口函数法设计FIR数字滤波器的设计方法，不同窗函数对滤波器滤波性能的影响以及滤波器单位冲激响应长度对滤波器延时特性的影响；3） 能够利用傅立叶变换在理论上分析滤波器的频率响应；4） 掌握计算机C语言或MATLAB的编程技巧；三、 进度计划序号设计(实验)内容完成时间备注1熟悉利用窗口函数法设计FIR带阻滤波器的方法和步骤，针对具体题目进行设计2天2编制程序2天3调试， 答辩1天四、 设计（实验）成果要求设计报告五、 考核方式答辩+平时表现学生姓名：杨睿指导教师：白康 2014 年1月12日数字信号处理课程设计一、数字信号处理课程设计的目的与要求能够运用本课程中学到的知识，设计基于窗口
3、函数法的FIR 数字滤波器。要求掌握数字信号处理的基本方法；FIR滤波器的设计步骤和方法；能够熟练采用C 语言或MATLAB语言进行计算机辅助设计和仿真验证设计内容的合理性。二、正文 1设计题目：采用窗口函数法设计一个带阻FIR数字滤波器，要求设计的频率响应为2设计步骤（1）根据设计要求，确定FIR滤波器的频率响应，包括过渡带宽、通带和阻带阻尼限制： 由设计题目可得知：滤波器的通带:0.55;通带宽度:0.45和0.通带宽度0.2;过渡带:,宽度0.1和，过渡带宽0.05；阻带:,宽度0.2.阻带最小衰减:。（2）根据前一步的过渡带宽、通带和阻带阻尼限制，确定选择的窗口函数类型 和的列长度：
4、由于阻带最小衰减要求为-33.98dB，根据窗函数基本参数可查的汉宁窗的阻带最小衰减为-44dB,符合要求且在所有符合要求中最简,因此选择汉宁窗.同样可以查的汉宁窗的过渡带宽为8/N.所以0.05=8/N,N=8/(0.05)=160,取N=160. 同样可以查得汉宁窗（3）根据第一步确定理想的频率响应，通过求的逆傅里叶变换得到，对进行右移、截短和窗口函数处理，求得因果FIR数字滤波器的单位冲激响应 : 由已知的条件可以得出理想得频率响应 。右移、截短 .因此因果FIR数字滤波器的单位冲激响应* （4）求得傅里叶变换，检验是否符合要求； .各部分MATLAB的仿真曲线如下:（1）.理想 b).
5、汉宁窗:3). 实际， .d)幅度响应（db）: 由此可以得知设计的滤波器阻带最小衰减为-44db,验证了汉宁窗的阻带最小衰减为-44db,同时说明该设计的滤波器在阻带最小衰减的性能指标满足设计要求的-33.98db.滤波器的性能检验:1).f1=3*sin(0.15*pi*n)2).f2=3*sin(0.4*pi*n)3).f3=3*sin(0.75*pi*n)4)叠加后的图像f4=+f1+f2+f35).通过滤波器后的图像f4*h(n)6).f5=f1+f3分析：通过设计的带阻滤波器后，只有f1和f3的图形能通过滤波器，所以f5=f1+f3的图像应该和f4*h(n)的图像差不多。由上面的波
6、形图比较可知，该滤波器基本满足要求。三、收获和体会： 在本次课程设计之前，我已经在老师的指导下初步熟悉了MATLAB这个非常好用的软件。通过本次的课程设计，使我对上课学到的东西有了更深刻的认识，加深了对各种数字滤波器的了解，同时对其的设计方法及各种性能指标都有一定的了解。可以用MATLAB编写一些简单的程序，得到理想的波形。在课程设计的初期，由于对MATLAB软件的不熟悉，在编完程序后进行编辑运行的时候，会出现很多的问题，例如程序的错误，有些定义软件识别不了，少了加一个括号之类的，导致程序不能运行，不能输出波形；有时候坐标轴长度选取不合适，导致输出波形之后无法完全显示或者出现重叠，最后通过查找
7、相关资料，请教老师和同学，我一一解决了这些困难，终于比较完满的地完成了此次课程设计。这次课程设计使我学习到了很多知识，更加深刻地掌握了数字信号处理的相关知识，加深的对MATLAB软件的理解和应用。以后会更加的加强联系，掌握更多的技能。四、参考文献：1、数字信号处理原理及其MATLAB实现 从玉良 电子工业出版社 2、Matlab信号处理与应用 董长虹 国防工业出版社 3、数字信号处理及MATLAB实现 余成波 清华大学出版社五、源程序：MATLAB程序代码:clear all; Wpl=0.2*pi;Wph=0.55*pi;Wsl=0.3*pi;Wsh=0.5*pi;tr_width=min(
8、Wsl-Wpl),(Wph-Wsh); N=ceil(8*pi/tr_width)+1 n=0:N-1;Wcl=(Wsl+Wpl)/2; Wch=(Wsh+Wph)/2;for i=0:N-1 hd1(i+1)=sin(Wcl*(i-(N-1)/2+eps)./(pi*(i-(N-1)/2+eps); hd2(i+1)=(sin(pi*(i-(N-1)/2+eps)-sin(Wch*(i-(N-1)/2+eps)./(pi*(i-(N-1)/2+eps); hd(i+1)=hd1(i+1)+hd2(i+1);endw_hann=(hanning(N); w0=w_hann;h=hd.*w0;h
9、1,w=freqz(h,1); f1=3*sin(0.15*pi*n)subplot(321)stem(n,f1)title(f1=3*sin(0.15*pi*n)axis(0 N-1 -4 4);ylabel(f1);text(N,-0.1,n);hold onf2=3*sin(0.4*pi*n)subplot(322)stem(n,f2)title(f2=3*sin(0.4*pi*n)axis(0 N-1 -4 4);ylabel(f2);text(N,-0.1,n);hold onf3=3*sin(0.75*pi*n)subplot(323)stem(n,f3)title(f3=3*si
10、n(0.75*pi*n)axis(0 N-1 -4 4);ylabel(f3);text(N,-0.1,n);hold onf4=f1+f2+f3subplot(324)stem(n,f4)title(f4=f1+f2+f3)axis(0 N-1 -10 10);ylabel(f4);text(N,-0.1,n);hold onf=conv(f4,h)subplot(325)k=(0:320)stem(k,f)title()axis(0 320 -8.0 8.0);ylabel(f4*hn);text(N,-0.1,k);hold onf5=f1+f3subplot(326)stem(n,f5
11、)title(f5=f1+f3)%axis(-20 200 -8 8);ylable(f5);text(N,-0.1,n);hold on figure(2)subplot(221)stem(n,hd)title(hd(n)axis(0 N-1 -0.5 1);ylabel(hd(n);text(N,-0.1,n);hold onsubplot(222)stem(n,w_hann)title(w(n)axis(0 N-1 0 1);ylabel(w(n);text(N,-0.1,n);hold onsubplot(223)stem(n,h)title(h(n)axis(0 N-1 -0.5 1
12、);ylabel(h(n);text(N,-0.1,n);hold onsubplot(224)plot(w/pi,20*log(abs(h1)/abs(h1(1);title(dB)ylabel(db);text(N,-0.1,t);hold on 三、分析总结： 通过本次高通FIR数字滤波器的设计，加深了对数字滤波器的了解，同时对其的设计方法及各种性能指标都有一定的了解。由阻带最小衰减确窗函数的类型，由过渡带宽确定窗函数的点数N，由理想频率响应通过傅立叶逆变换得到和hd(n)截断右移，由窗函数特性确定的w(n),加窗后得到实际h(n)=hd(n)*w(n) 。最后验证各项指标（包括阻带最小衰减，过渡带宽，通带阻带以及截至频率）。 12
