信号分析与处理课程设计67713045.docx

课程设计报告( 2014 - 2015年度第一学期)名 称：信号分析与处理a课程设计 题 目： fir数字滤波器的设计 院 系： 自动化 班 级： 测控1201 学 号： 201202030102 学生姓名： 蔡文斌 指导教师： 金秀章 设计周数： 第22 周 成 绩： 日期：2015年 01月23日信号分析与处理a课程设计任务书一、 目的与要求能够运用本课程中学到的知识，设计基于窗口函数法的fir 数字滤波器。要求掌握数字信号处理的基本方法；fir滤波器的设计步骤和方法；能够熟练采用c 语言或matlab语言进行计算机辅助设计和仿真验证设计内容的合理性。二、 主要内容2. 请采用窗口函数法设计一个低通fir数字滤波器，要求设计的频率响应为1-0.02hej1+0.02 ,0.3 hej0.02 ,0.4 三、进度计划序号设计(实验)内容完成时间备注1领取设计任务 查阅相关资料2015.1.192制定滤波器设计的方案2015.1.203编程实现2015.1.20-2015.1.224答辩2015.1.23四、设计（实验）成果要求课程设计报告五、考核方式设计报告+答辩学生姓名：蔡文斌指导教师：金秀章2015年1 月19日课程设计正文一、 课程设计的目的与要求能够运用本课程中学到的知识，设计基于窗口函数法的fir 数字滤波器。要求掌握数字信号处理的基本方法；fir滤波器的设计步骤和方法；能够熟练采用c 语言或matlab语言进行计算机辅助设计和仿真验证设计内容的合理性。二、 设计正文滤波器的设计1、所需设计的滤波器要求采用窗口函数法设计一个低通fir数字滤波器，要求设计的频率响应为1-0.02hej1+0.02 ,0.3 hej0.02 ,0.4 2、分析所需设计的滤波器并确定序列长度本题要求我们设计一个低通fir滤波器，其要求的波动范围为0.02，由20log0.02-34db，查表，根据最小阻带衰减，所以选择海宁窗。该滤波器的通带截止频率p=0.3，阻带截止频率s=0.4，所以截止频率c=12p+s=0.35rad。过渡带宽=s-p=0.1rad。根据海宁窗的过渡带宽=8n，可以确定序列的长度为80，由于n必须为奇数，所以取n=81。3、构造理想的频率响应将此滤波器化为一个理想低通数字滤波器，其频率特性为hdej=1， 0.35 0, 0.35 4、计算数字滤波器的单位采样响应并对其进行加窗处理通过对理想低通数字滤波器的频率特性做傅里叶变换可得理想单位采样响应即hdn=12-ccejnd=sincnn=sin0.35nn所选海宁窗窗函数为hnn=0.51+cos(2（n-(n-1)2)n-1),0n(n-1)0, 其他 运用matlab对所选窗函数进行仿真验证： 对海宁窗移位并施加海宁窗，得：hn=sinn-40*0.35(n-40)*0.51+cos2n-4080 0n80运用matlab仿真可得单位采样响应如下：5、对单位冲激响应做傅里叶变换，检验其是否合格单位采样响应的频率响应为hej=n=0n-1hne-jn=n=080sinn-40*0.35(n-40)*0.51+cos2n-4080e-jn利用matlab进行仿真可得：在图中可清楚看出，通带截止频率p=247800=0.309，阻带截止频率s=314800=0.393，其通带为（00.309），阻带为（0.393），所设计的滤波器通带阻带都包含题目所要求的。由程序可知，在通带（00.3 ）中，max=1.00640.98，通带满足设计要求。由程序可知，在阻带（0.4）中，max= 0.0064-0.02,阻带满足设计要求。综上，所设计的滤波器满足要求。6、构造通带和阻带内不同频率和不同幅值的正弦序列组合而成的输入信号序列，用所设计的滤波器对其进行滤波，验证设计滤波器的合理性。通带内取值得x1=12sin(0.12n)，阻带内取值得x2=23(sin0.6n)，x=x1+x2。仿真如图：分析:由上图可知，滤波后的信号y(n)与信号x(n)波形相同。证明所设计的滤波器性能良好，满足设计要求。三、 课程设计的总结或结论通过本次低通fir数字滤波器的设计，加深了对数字滤波器的理解，以及更加清楚了数字滤波器的在实际生活中的应用，也让我明白了信号分析与处理这门课的意义所在。在这次课设中，运用了在课堂上学到的有限冲激响应数字滤波器设计的知识，以及运用在课下学到的通过matlab绘图分析的知识。另外，这次课设让我学到了一些课本上没有的知识，以及加深了对数字滤波器设计的步骤，收获颇丰。四、 参考文献1崔翔主编.信号分析与处理. 中国电力出版社, 第二版. 2011年2 张志涌，杨祖樱.matlab教程. 北京航空航天大学出版社, 2010年附录：matlab程序clear;n=81;wc=0.35*pi;for n=1:n w(n)=0.5*(1+cos(pi/40*(n-41); if n=41 hd(n)=wc/pi; h(n)=wc/pi; else hd(n)=sin(n-41)*wc)/(pi*(n-41); h(n)=sin(n-41)*wc)/(pi*(n-41)*w(n); endendm=0:pi/800:pi;for n=1:801; h(n)=0; for k=1:n; h(n)=h(n)+h(k)*exp(-j*m(n)*(k-1); end h(n)=abs(h(n);endfor n=1:200 x1(n)=12*cos(0.12*pi*n) x2(n)=23*cos(0.6*pi*n); x(n)=x1(n)+x2(n);endfor m=2:281 y(m)=0; for k=1:n if mk if km-201 y(m)=y(m)+x(m-k)*h(k); end end endendfigure(1)plot(w);xlabel(n);ylabel(w(n);title(海宁窗)figure(2)plot(h);title(实际单位采样响应)xlabel(n);ylabel(h(n)figure(3)plot(h);title(单位冲激响应幅度特性)xlabel(w（*pi/800）)ylabel(h(w)max=max(h(1:241)min=min(h(1:241)max=max(h(321:800)min=min(h(321:800)figure(4)subplot(311)plot(x1);title(曲线x1)xlabel(n);ylabel(x1);axis(60 160 -20 2