数字滤波器理论与设计.分析报告.ppt

H(ω)与原理想低通 Hd(ω)差别有以下2点： H(ω)在 ω=ωC附近形成过渡带，过渡带宽度B=4?/N，近似于矩形序列幅度谱RN(ω)的主瓣宽度； 通带内增加了波动，最大的峰值在 ω=ωC?2?/N 处，阻带内产生了余振，最大的负峰值在 ω=ωC+2?/N 处。幅度谱RN(ω)波动越快(N加大)，通带、阻带内波动越快，其旁瓣的大小直接影响H(ω)波动的大小。 Hd(ω)在加窗后在频域中的现象称为吉布斯效应 影响: (1)通带内的波动影响滤波器通带的平稳性； (2)阻带内波动影响阻带的衰减，可使最小衰减不满足技术要求； 2、改善窗函数的形状 减少带内波动以及加大阻带的衰减只能从窗函数的形状找出解决方法，主要考虑以下2点因素： 尽量减小主瓣宽度，以获得较窄的过渡带； 尽量使窗函数的最大副瓣相对于主瓣要小，使设计出来的滤波器幅度特性中肩峰和余振较小，阻带衰减较大。 三、窗函数法的设计步骤 1．确定希望逼近的滤波器的频响函数Hd(ejω) 2、根据Hd(ejω)确定其对应的单位脉冲响应hd(n) IDFT RM(n) n (1) Hd(ejω)可封闭求解，则： (2) Hd(ejω) 不可封闭求解，对Hd(ejω) 从 ω=0?2?采样M点，采样值为 Hd(ej2?k/M)，k=0,1,…,M-1，用2?/M 代替上式中dω ，则： (3) 如果已知通带(或阻带)衰减和边界截止频率ωc，选理想滤波器作为逼近函数，对理想滤波器频响函数作 IFT ，求出 hd(n)。 2、选择窗函数 根据过渡带与阻带衰减的要求，选择满足条件的窗函数形式，并估计窗口长度N。原则是保证阻带衰减的前提下，尽量选主瓣窄的窗函数。 3、计算所要设计的滤波器的单位采样响应h(n) 计算滤波器的单位取样响应 h(n)= hd(n)?w(n)。其中 w(n)是上面选择好的窗函数，hd(n)与w(n)都应满足线性相位要求。 4、验证技术指标是否满足要求 已设计出的滤波器的频率响应 。验算H(ejω)是否满足设计要求，若不满足要求，重复上面2，3，4过程。 窗函数法优点： 从时域出发的一种设计方法，设计简单，方便，实用。 缺点是： 要求用计算机实现，边界频率不易控制。 Hd(ejω) h(n) hd(n) IFT 加窗截断 FT 比较，满足设计要求，则设计完毕，不合格则修改窗函数 H (ejω) 一、频率采样法基本原理 设 Hd(ejω) 为所要设计的数字滤波器的频率响应，则: 1、在?=???? 范围内对 Hd(ejω) 进行N点等间隔采样，得到Hd(k) 用频率采样法设计FIR滤波器 2、对N点Hd(K)进行IDFT，得所设计的滤波器的单位 脉冲响应h(n) 3．对求出的h(n)进行z变换，得到滤波器的系统函数H(z) 或利用频域内插公式 Hd(ejω) Hd(k) h(n) H(ejω) 等间隔采样 IDFT 内插函数 是否满足指标 （2）h(n)偶对称，N为偶数时： 设计方法： 用理想滤波器作为逼近滤波器，截止频率为ωc，采样点数为N，则Hk和θk的计算公式为： （1）h(n)偶对称，N为奇数时： Hk = HN-k = 1， k=0，1，…，kc； Hk = 0， k= kc+1，kc+2，…，N-kc-1； k=0，1，…，N-1； Hk = 1， k=0，1，…，kc； Hk = 0， k= kc+1，kc+2，…，N-kc-1； k=0，1，…，N-1； HN-k = -1， k=0，1，…，kc； kc取小于等于 ωcN /2?的最大整数 对 H(k) 进行 IDFT变换，求出 h(n)； 由 h(n)可求出所设计滤波器的频响 H(ejw) 由前面计算出的Hk和θk的值可构造出H(k) 等效于在[0，2π] 上的N个采样值 k=0，1，…，N-1； 总结：频率采样法设计线性相位FIR滤波器步骤 1、根据ωc 及N的奇偶性，确定滤波器 Hk和 ?k及kC值； 2、由H(k) = Hkej?k，求出 H (k)； 3、对 H (k) 进行 IDFT变换，求出 h(n)； 4、由 h(n) 求出所设计的滤波器的频率响应 H(ejw)，并分析误差，优化设计。 主要内容 数字