物理与电子信息工程学院
实验报告
实验课程名称：数字信号处理
实验名称： FIR数字滤波器设计与软件实现 班 级： 1012341 姓 名： 严娅 学 号： 101234153 成 绩：_______
实验时间： 2012年12月20 日
一、实验目的
（1）掌握用窗函数法设计FIR数字滤波器的原理和方法。
（2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。（3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。
二、实验原理
1、用窗函数法设计FIR数字滤波器的原理和方法。
如果所希望的滤波器的理想频率响应函数为 Hd(ej?)，则其对应的单位脉冲响应为
hd(n)=
12?????Hd(ej?)ej?d? （2-1）
窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼近hd(n)。由于hd(n)往往是无限长序列，且是非因果的，所以用窗函数?(n)将hd(n)截断，并进行加权处理，得到：
h(n)＝hd(n)?(n) (2-2)
h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，
其频率响应函数Hd(ej?)为：
Hd(e)＝?h(n)e?j? (2-3)
j?n?0N?1式中，N为所选窗函数?(n)的长度。
由第七章可知，用窗函数法设计的滤波器性能取决于窗函数?(n)的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。
这样选定窗函数类型和长度N后，求出单位脉冲响应h(n)＝hd(n)·?(n)，并按式（2-3）求出H(ej?)。H(ej?)是否满足要求，要进行验算。一般在h(n)尾部加零使长度满足于2的整数次幂，以便用FFT计算H(ej?)。如果要观察细节，补零点数增多即可。如果H(ej?)不满足要求，则要重新选择窗函数
类型和长度N，再次验算，直至满足要求。
如果要求线性相位特性，则h(n)还必须满足 h(n)??h(N?1?n) (2-4)
根据上式中的正负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据设计的滤波特性正确选择其中一类。例如，要设计线性低通特征，可选择h(n)?h(N?1?n)一类，而不能选h(n)??h(N?1?n)一类。 2、等波纹最佳逼近法
等波纹最佳逼近法是一种优化设计法，它克服了窗函数设计法和频率采样法的缺点，使最大误差（即波纹的峰值）最小化，并在整个逼近频段上均匀分布。用等波纹最佳逼近法设计的FIR数字滤波器的幅频响应在通带和阻带都是等波纹的，而且可以分别控制通带和阻带波纹幅度。这就是等波纹的含义。最佳逼近是指在滤波器长度给定的条件下，使加权误差波纹幅度最小化。与窗函数设计法和频率采样法比较，由于这种设计法使滤波器的最大逼近误差均匀分布，所以设计的滤波器性能价格比最高。阶数相同时，这种设计法使滤波器的最大逼近误差最小，即通带最大衰减最小，阻带最小衰减最大；指标相同时，这种设计法使滤波器阶数最低。实现FIR数字滤波器的等波纹最佳逼近法的MATLAB信号处理工具函数为remez和remezord。Remez函数采用数值分析中的remez多重交换迭代算法求解等波纹最佳逼近问题，求的满足等波纹最佳逼近准则的FIR数字滤波器的单位脉冲响应h（n）。
三、 实验内容及步骤
（1）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示；