数字信号处理实验 fFIR IIR 滤波器的设计.doc

数字信号处理实验报告 FIR和IIR滤波器的设计 实验目的 分别利用脉冲响应不变法、频率采样法设计IIR滤波器，观察滤波器形状，并进行比较。 双线性变化法设计FIR滤波器。 二、实验条件 计算机、MATLAB2013a 实验内容 、IIR滤波器的设计。 IIR滤波器的特征是，具有无限持续时间冲激响应。这种滤波器也称之为递归型滤波器。这些滤波器产生新的输出，不但需要过去的输入和现在的输入，还需要过去的输出。其相应的差分方程为 ，假设，整理之后便得到滤波器的方程，由所学知识很明显看出该方程为递归的，系数的存在意味着递归滤波器有无限项的脉冲响应。相较与非递归滤波器，递归滤波器设计时所需的系数要少的多，更容易实现，一般会选择具有待求特性的原型模拟滤波器，然后将其转化为数字滤波器。IIR 数字滤波器的设计常采用两种方法 ：冲激响应不变法、双线性变换法。 利用脉冲响应不变法设计IIR滤波器。 、原理：每个模拟滤波器都有冲击响应h(t)，就像每个数字滤波器都有脉冲响应h[n]一样，脉冲响应不变法选择的数字脉冲响应h[n]是满足设计要求的模拟滤波器冲击响应h(t)的采样值。即h[N]=h(nT),其中T为所用的采样间隔。如果令 H a( s)是h ( t)的拉普拉斯变换，H (z)为h( n )的z变换，利用采样序列z变换和模拟信号的拉氏变换之间的关系，即，从而将模拟滤波器的s面转换为数字滤波器的z平面。 、设计步骤： 、确定滤波器各个参数，截止频率fs、阶数、采样频率Fs等； 、将数字滤波器频率指标转换成相应的模拟滤波器频率指标 将数字频率和模拟频率联系起来； 、设计一个阶数为1的巴特沃斯滤波器，其截止频率为750HZ，采样频率为4000HZ； 、利用impinvar和butter函数，带入滤波器系数，具体代码如下： 其中butter的具体使用方法如注释，当设计高通滤波器时，将‘s’改为‘high’即可。运行结果如图： 分析结果发现与实际课本上的滤波器略有不同，判断可能是由于采样频率过大导致，修改采样频率后结果如下： 利用双线性变换设计IIR滤波器。 原理分析：双线性变换的定义为 ，fs为采样频率，用代替z，利用欧拉公式，最终变换得 ，称为预扭曲方程，从而将模拟角频率与数字角频率联系起来，此时w与是非线性关系，只有当趋于0时，为线性关系,因此对于相同参数的滤波器，其对应的模拟滤波器和数字滤波器并不能达到完全吻合。 、滤波器的具体参数同上，1阶巴特沃斯滤波器，采样频率为4HZ（由上个实验发现采样频率不能超过10HZ，否则会出现严重的失真）。 、利用bilinear函数，对模拟滤波器实现双线性变换，代码如下： 结果如下： 将两种设计方法进行对比，并之后与对应得模拟滤波器相比较： 和与之对应的模拟滤波器对比如下： 结果分析： 其中在采样频率为4HZ的时候，分别用双线性和脉冲响应不变法设计IIR滤波器，红色曲线为脉冲不变法，蓝色为双线性法，可知 ①、脉冲响应不变方法设计滤波器在通频带的增益要小但是其阻带频率较高，衰减幅度大，滤波性相对较好； ②、双线性变换法在通频带其增益较高但阻带频率高，在实际的应用中可能不能很好地实现滤除噪声的功能。 ③、由于采样间隔较大，导致采样点数很少，对比模拟滤波器很容易发现数字滤波器不能很好地趋近于模拟滤波器。 结合所学知识，对于脉冲响应不变法和双线性变换法设计滤波器的比较如下： 脉冲不变： ①、模拟频率到数字频率的转换时线性的，不需要借助预扭曲方程； ②、数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应，因此时域特性逼近好； ③、采样间隔越小，越能逼近模拟滤波器的形状； ④、会产生频谱混叠现象，只适合带限滤波器；双线性变换法： ①、克服多值映射得关系，可以消除频率的混叠； ②、表示模拟传输函数的滤波器必须进行预扭曲； ③、是非线性的，在高频处有较大的失真。 、FIR的数字滤波器的设计 有限脉冲响应滤波器又称为非递归滤波器，其输出仅取决于过去的输入，与过去的输出无关，与之对应的差分方程为，对比之前的IIR不难发现，其脉冲响应为有限项，，包括M+1个脉冲函数。一般来说，滤波器的滚降特性越好，其需要的系数就越多，在FIR滤波器的设计中，如果要想达到满意的性能，一般需要100~200个系数，对比前面提到的IIR滤波器系数较多，但是非递归滤波器具有递归滤波器所没有的一些优点，它可以保证稳定，消除通带内的相位失真，这使得它更适合音频应用。在实际应用中，FIR滤波器可以利用窗函数或者频率采样法设计，窗函数设计FIR滤波器，上个实验已经做了具体的观察，此次实验主要采用的是频率采样法设计滤波器。 原理分析：频率采样法的基本思想是使所设计的FIR数