数字信号处理课程设计技术总结.docVIP

河北科技大学 课程设计报告 学生姓名： 张亚洲 学 号： 12L0751130 专业班级： 电信L123 课程名称： 数字信号处理课程设计 学年学期： 2 014 —2 015 学年第2学期 指导教师： 张秀清 2 015 年 7 月 课程设计成绩评定表 学生姓名张亚洲学 号12L0751130成绩专业班级电信L123起止时间2015.7.6-2015.7.10设计题目1. 窗函数法设计FIR数字低通滤波器 2. DFT在信号频谱分析中的应用指 导 教 师 评 语 指导教师： 年 月 日 目 录 （一）窗函数法设计FIR数字低通滤波器 一、设计题目……………………………………………………………1 二、设计目的……………………………………………………………1 三、设计原理……………………………………………………………1 四、实现方法……………………………………………………………1 五、设计内容及要求……………………………………………………1 六、运行结果……………………………………………………………6 七、心得体会……………………………………………………………9 （二）DFT在信号频谱分析中的应用 一、设计题目………………………………………………………………10 二、设计目的………………………………………………………………10 三、设计原理………………………………………………………………10 四、设计源程序……………………………………………………………10 五、运行结果………………………………………………………………14 六、参考文献………………………………………………………………20 PAGE \* MERGEFORMAT21 一、设计题目：窗函数法设计FIR数字低通滤波器 二、设计目的 1.掌握用窗函数法设计FIR数字滤波器的原理和方法。 2.熟悉线性相位FIR数字滤波器的特性。 3.了解各个窗函数对滤波器特性的影响。 三、设计原理 窗函数设计FIR DF的基本思想是根据给定的滤波器技术指标，选择滤波器长度N和窗函数w(n),使其具有最窄宽度的主瓣，和最小的旁瓣，设计步骤为： 1.由指标要求，确定滤波器长度N（窗长N）和窗函数类型。 2.在进而确定理想低通hd(n),所以h(n)=hd(n)×w(n) 3.由h(n)求H(ejw)验证所求滤波器是否满足给定指标要求. 四、实现方法\ 用malab编程实现给定指标要求的滤波器设计 五、设计内容及要求 1.熟悉各种窗函数，在matlab命令窗口键入一下命令行，浏览各种窗函数， 将图写入。 2.编写求理想低通滤波器中冲激响应hd(n)的MATLAB程序： hd(n)= function[hd]=ideal(Wc,N) %Wc= 低通截频 %N为滤波器长度，且X=(N-1)/2 %hd为0~N-1 范围 内的理想响应。 q=(N-1)/2; n=0:N-1; m=n-q+eps; hd=sin(wc*m)./(pi*m);存盘，ideal.m function[H]=fr(b,a,w); m=0:length(b)-1; l=0:length(a)-1; num=b*exp(-j*m*w); den=a*exp(-j*l*w); H=num./den; 存盘，fr.m 3.调用编写的函数文件ideal.m及fr.m设计给定指标要求的线性相位FIR滤波器指标：Wp=0.2π，Rp=0.25dB(通带最大波动) Ws=0.3π， As=50dB(阻带最小衰减) 要求：①根据指标要求选择窗函数类型 ②编程序计算窗长Ni 0≤n≤N-1范围的理想低通冲激响应hd(n)实际滤波器的冲激响应h(n)，实际的滤波器的频幅特性|H(ejw)|并用分贝表示。 ③验证所用设计的DF是否满足指标要求，计算通带内的最大波动和阻带内的最小衰减。 窗函数选择： 由As=5dB(阻带最小衰减)查《数字信号处理》知hanming与blackman都可以选择设计，此文选hanming窗，旁瓣峰值-41dB,主瓣宽度8π/N，过滞带宽,阻带最小衰减为-53dB. 冲激响应; hd(n)= 六、运行结果 N=110 Rp1=0.0036dB (通带