《数字信号处理》实验6 FIR数字滤波器设计.pptxVIP

《数字信号处理实验与课程设计教程--面向工程教育》实验六 FIR数字滤波器设计实验目的：本实验包含3个实验子项目，通过matlab编程仿真，掌握FIR数字滤波器设计方法，包括：1.掌握用窗函数法设计FIR数字滤波器，并了解各种窗函数对滤波特性的影响。2.熟悉线性相位滤波器的特点。3.掌握用频率采样法设计FIR数字滤波器。4.熟悉FIR数字滤波器的应用实例。 实验六 FIR数字滤波器设计6.1 窗函数法设计FIR数字滤波器6.2 频率采样法设计FIR数字滤波器6.3 FIR数字滤波器的应用实例6.1 窗函数法设计FIR数字滤波器一、实验原理1. 窗函数法设计FIR滤波器的步骤(1)根据技术指标设计线性相位的理想滤波器hd(n)。理想低通滤波器： hd(n)= (6.1)其中:ωc=(ωp+ωs)/2---- hd(n)的截止频率;τ=(N-1)/2--- hd(n)的对称中心；N--- hd(n)的阶数。(2) 根据技术指标：阻带最小衰减as选择适当的窗函数，对hd(n)进行加窗处理得到实际FIR滤波器单位脉冲响应：h(n)=hd(n)w(n)。其中:w(n)---窗函数。常用窗函数有：矩形窗、三角窗、汉宁窗、哈明窗、布莱克曼窗、凯塞窗等。Matlab函数：1)矩形窗： w=boxcar(N)其中:N---窗函数的长度；w---窗函数w(n)。2)三角窗： w=bartlett(N)3)汉宁窗： w=hanning(N)4)哈明窗： w=hamming(N)5)布莱克曼窗：w=blackman(N)6)凯塞窗：w=kaiser(N,beta)其中：beta---凯塞窗的调整参数。2.高通/带通/带阻FIR滤波器设计(1)理想高通FIR滤波器 hd(n)= (6.2)(2) 理想带通FIR滤波器 hd(n)= (6.3)其中:ωc1=(ωp1+ωs1)/2, ωc2=(ωp2+ωs2)/2---- 带通滤波器hd(n)的截止频率。(3) 理想带阻FIR滤波器 (6.4) 3.线性相位滤波器(1)第一类线性相位滤波器 (6.5)(2)第二类线性相位滤波器 (6.6)4. FIR滤波的实现方法FIR滤波的输出是滤波器输入信号与滤波器单位脉冲响应的卷积： y(n)=x(n)*h(n) (6.7)其中：x(n)---输入信号；h(n)---FIR滤波器的单位脉冲响应。Matlab函数：y=conv(x,h)或y=filter(h,1,x)其中：x---输入信号；h---FIR滤波器的单位脉冲响应h(n);y---输出信号。二、实验环境1.计算机1台。2. Windows 7 以上操作系统。3. Matlab 7.0 以上版本软件。三、实验参考与实验内容1.窗函数法设计FIR低通滤波器 采用窗函数法设计一个线性相位低通滤波器，通带截止频率wp=0.3πrad,阻带截止频率ws=0.5πrad,要求阻带衰减as=40dB,选择适当的窗函数及长度，求：滤波器的单位脉冲响应h(n)。理论分析：由于阻带衰减as=40dB,可选择汉宁窗、哈明窗、布莱克曼窗。%实验参考程序：参见: 实验项目参考源码\ch6\ch6prog1.m实验要求：(1)运行程序，画出理想低通滤波器hd(n)、汉宁窗、实际滤波器h(n)及其幅频特性H(w)波形。(2)观察hd(n)与h(n),这两个波形有何不同？在程序中添加语句，画出hd(n)的幅频特性Hd(w)，比较Hd(w)与H(w),这两个波形又有何不同？（提示：hd(n)其实是理想低通滤波器加矩形窗）(3)在命令窗口中读出滤波器实际的通带衰减ap及阻带衰减as,分别等于？是否符合设计指标要求？(4)在程序中修改N的表达式及窗函数，改为哈明窗及布莱克曼窗，重新运行程序，问此时，设计出的滤波器实际阻带衰减as=?由此得出什么结论？图6.1 窗函数法设计FIR数字低通滤波器2.窗函数法设计FIR高通/带通/带阻滤波器(1)FIR高通滤波器设计 采用窗函数法设计一个线性相位高通滤波器，通带截止频率wp=0.5πrad,阻带截止频率ws=0.3πrad,要求阻带衰减as=50dB,选择哈明窗，求：滤波器的单位脉冲响应h(n)。理论分析：高通滤波器=全通滤波器-低通滤波器，可根据公式(6.2)设计其h(n)。实验要求：(1)编写程序：ch6prog2.m,设计此FIR高通滤波器，画出其理想低通滤波器hd(n)、哈明窗、实际滤波器h(n)及其幅频特性H(w)波形，并列出程序清单。(实验参考程序：参见: 实验项目参考源码\ch6\ch6prog2.m)(2)在命令窗口中读出滤波器实际的通带衰减ap及阻带衰减as,分别等于？是否符