实验8：FIR数字滤.ppt

实验8：FIR数字滤波器设计及应用 实验目的 1、掌握用窗函数法设计FIR数字滤波器的原 理和方法。 2、熟悉线性相位FIR数字滤波器特性。 3、了解各种窗函数对滤波器特性的影响。 实验原理 1、有限冲激响应数字滤波器的基础理论(请 参考《数字信号处理教程》)。 2、窗函数法设计FIR滤波器的原理。 3、数字滤波器系数的确定方法。 4、根据要求设计低通FIR滤波器。 要求：通带边缘频率10kHz，阻带边缘频率22kHz，阻带衰减75dB，采样频率50kHz。 设计过程： (1)过渡带宽度=阻带边缘频率-通带边缘频率=22-10 =12kHz (2)理想低通滤波器的截止频率： fc=通带边缘频率+(过渡带宽度)/2 =10000+12000/2=16kHz Ωc=2πf1/fs=0.64π (3)理想低通滤波器脉冲响应： h1[n]=sin(nΩ1)/n/π=sin(0.64πn)/n/π (4)根据阻带衰减要求，选择布莱克曼窗，窗函数长度为： N=5.98fs/过渡带宽度=5.98*50/12=24.9 选择N=25，窗函数为： w[n]=0.42-0.5cos(2πn/24)+0.08cos(4πn/24) (5)滤波器脉冲响应为： h[n]=h1[n]w[n] |n|≤12 h[n]=0 |n|＞12 (6)根据上面计算，各式计算出h[n]，然后将脉冲响应值移位为因果序列。 b = 0.0000 -0.0000 0.0008 -0.0016 -0.0019 0.0100 -0.0087 -0.0180 0.0492 -0.0204 -0.1095 0.2800 0.6400 0.2800 -0.1095 -0.0204 0.0492 -0.0180 -0.0087 0.0100 -0.0019 -0.0016 0.0008 -0.0000 0.0000 (7)完成的滤波器的差分方程为： y[n]=0.001x[n-2]-0.002x[n-3]-0.002x[n-4] +0.01x[n-5] -0.009x[n-6]-0.018x[n-7] +0.049x[n-8]-0.02x[n-9]-0.11x[n-10] +0.28x[n-11]+0.64x[n-12]+0.28x[n-13] -0.11x[n-14]-0.02x[n-15]+0.049x[n-16] -0.018x[n-17]-0.009x[n-18]+0.01x[n-19] -0.002x[n-20]-0.002x[n-21]+0.001x[n-22] 程序流程图 实验步骤 1、实验准备： 设置软件仿真模式； 启动CCS。 2、打开工程，浏览程序：工程目录为 C:\ICETEK\VC5509Ae\VC5509Ae\ Lab0501-FIR\Fir.pjt。 3、编译并下载程序。 4、打开观察窗口。 选择菜单View-GraphTime/Frequency…， 进行如下设置： 5、设置断点：在有注释“break point”的语句设置软件断点。 6、运行并观察结果： ⑴ 选择“Debug”菜单的“Animate”项，或按Alt+F5键运行程序。 ⑵ 观察“Input”、“Output”窗口中时域图形；观察 滤波效果。 ⑶ 鼠标右键单击“Input”和“Output”窗口，选择“Properties…”项，设置“Display Type”为“FFT Magitude”，再单击“OK”按钮结束设置。 ⑷ 观察“Input”、“Output”窗口中频域图形；理解滤波效果。 7、退出CCS。 实验结果 * *