的fir滤波器设计.docVIP

《DSP原理与技术》 实 验 报 告 姓 名：XXX 学 号：XXXX 院 系：航天学院控制科学与工程系 指导老师：XXX 2011年5月25日 1 FIR滤波器的设计与DSP实现 1.1 FIR滤波器的理论基础 FIR滤波器即有限长冲击响应滤波器的单位冲击响应是有限长的（）,其Z变换为 在有限Z平面（）内有（）个零点，在Z平面的原点处有（）阶极点。 的频率响应为 当单位冲击响应序列为实序列时，表示成 当具有线性相位即 或 当满足条件可得出 上式即为FIR滤波器具有线性相位条件的充要条件，它要求单位冲击响应序列是以为对称中心的偶对称序列。当满足条件可得出 此时要求单位冲击响应序列是以为对称中心的奇对称序列。 鉴于工程实际应用，具有线性相位的FIR滤波器具有更加广阔的应用。对于非线性相位的FIR滤波器，一般可以用IIR滤波器来代替。同样幅度特性，IIR滤波器所需阶数比FIR滤波器所需阶数少得多。 1.2 FIR滤波器的设计 常见的FIR滤波器设计方法包括窗函数设计法和频率采样设计法。这里讨论用窗函数设计方法设计数字低通滤波器。 理想的数字低通滤波器的频率响应如图1.2.1所示。根据傅立叶逆变换，由图1.2.1可得其单位脉冲响应为 此时是一个无限长的偶对称序列。 窗函数的设计即是用有限长度的窗去截断，使之变成有限长度序列。假设长度为M的窗函数的单位冲击响应为，则截断后的序列为 为保证系统的因果性，须将序列右移M个采样周期即得到物理可实现的实际数字滤波器的单位脉冲响应。 下面以借助于Matlab的滤波器设计工具FDATool设计一个截至频率为30KHz、系统采样频率为300KHz的数字低通滤波器为例来说明数字滤波器的一般设计步骤。 在Matlab命令行输入位置输入 fdatool 即进入FDATool对话框界面，如图1.2.2所示。 图1.2.2 FDATool设计界面 在滤波器响应类型Response Type栏选择Lowpass，设计方法Design Method选择FIR和Window，在Options一栏选择汉宁窗Hann，然后在频率设置一栏中Fs代表系统采样频率,这里填入，Fc代表截止频率填入30000后点击Designed Filter按钮即完成设计。图1.2.3即为设计的滤波器的幅频响应。 图1.2.3 幅频响应 FDATool工具可以将所设计的滤波器参数导出成Matlab能识别的m语言文件和Simulink Model文件用于仿真验证，也可以导出成C语言头文件（.H）用于其他平台的开发。图1.2.4是将所设计的滤波器导出Simulink Model后建立的验证性仿真结构框图。 图1.2.4 滤波器验证框图 图1.2.5a为加入频率为80KHz噪声后的信号功率谱，b为滤波后的功率谱，比较二者可知噪声在通过滤波器时完全被滤掉，基本达到滤波的功能。 a)滤波前 b)滤波后 图1.2.5 信号功率谱 1.3 FIR滤波器的DSP实现 将在Matlab的FDATool中设计好的滤波器参数导出成C语言头文件，从而在DSP开发环境中引用该参数即可完成在DSP中实现数字滤波器的设计。 附录一为FIR实验的源代码。实验中给出的信号是250Hz的有用信号和4.5KHz的噪声信号的叠加，采样率为10KHz。 1.3.1 FIR低通滤波器的设计实验 根据以上参数在FDATool中设计的FIR低通滤波器的幅频响应如图1.3.1所示。 图1.3.1 FIR低通滤波器，fc=300Hz 其中，该低通滤波器阶数为11，系数为：0, 0.392, 0.67, 0.33,0.6,0.02,0.6,0.33,0.67, 0.392, 0。 将附录一中的源代码中的滤波器阶数变量FIRNUMBER和滤波器系数变量fHn改成以上参数后，运行之，在CCS的Graph工具中观察滤波前后的效果，如图1.3.2所示。 a)低通滤波前的输入信号 b）低通滤波后的输入信号 图1.3.2 低通滤波前后的信号对比 通过观察以上两个波形并对比可得，原先输入信号中含有的高频分量被完全滤除了，说明所设计的低通滤波器各项指标符合实验要求。 1.3.2 FIR高通滤波器的设计实验 根据信号的各项参数在Matlab的FDATool工具中设计的高通滤波器幅频响应如图1.3.3所示。 图1.3.3 FIR高通滤波器，fc=4KHz 其中，该低通滤波器阶数为11，系数为：0, 0.9511,-0.214, 0.74,-0.87,0.46,-0.87, 0.74, -0.214, 0.9511, 0。 将附录一中的源代码中的滤波器阶数变量FIRNUMBE