数字信号处理吴镇扬DSPCH4分析.ppt

例 用频率采样法设计一线性相位低通FIR数字滤波器 0≤ω≤ωc 其他 ωc=0.5π，采样点数为奇数N=33。 试求各采样点的幅值Hk及相位θk，也即求采样值H(k)。 解 N=33, 且低通幅度特性H(0)=1。第一类线性相位滤波器 0≤k≤8, 25≤k≤32 9≤k≤24 0≤k≤32 频率采样法小结 优点是可以在频域直接设计，并且适合最优化设计 缺点是采样频率只能等于2π/N的整数倍，因而不能确保截止频率ωc的自由取值，要想实现自由地选择截止频率，必须增加采样点数N，但这又使计算量加大。 4.5 FIR滤波器和IIR滤波器的比较 首先，从性能上来说，IIR滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的。 从结构上看，IIR滤波器必须采用递归结构，极点位置必须在单位圆内，否则系统将不稳定。 从设计工具看，IIR滤波器可以借助于模拟滤波器的成果，因此一般都有有效的封闭形式的设计公式可供准确计算，计算工作量比较小，对计算工具的要求不高。 To Learn Chapter 5! To Learn Chapter 5! To Learn Chapter 5! 设计方法 0≤n≤N-1 其他 h(n)=hd(n)w(n) 窗函数法设计FIR数字滤波器是在时域进行 0≤n≤N-1 其他 窗函数为矩形窗 ---------窗口设计法 例如，要求设计一个FIR低通数字滤波器， 假设理想低通滤波器的频率响应为  |ω|≤ωc ωc|ω|≤π 理想低通的单位脉冲响应及矩形窗 矩形窗对理想低通幅频特性的影响 h(n)=hd(n)wR(n) 矩形窗 理想滤波器的频率响应 实际设计的FIR滤波器频率响应 矩形窗影响 加窗函数处理后， 对理想频率响应产生以下几点影响 （1）H(ω)将Hd(ω)在截止频率处的间断点变成了连续曲线， 使理想频率特性不连续点处边沿加宽，形成一个过渡带，过渡带的宽度等于窗的频率响应WR(ω)的主瓣宽度Δω=4π/N，即正肩峰与负肩峰的间隔为 4π/N。窗函数的主瓣越宽，过渡带也越宽。 （2）在截止频率ωc的两边即ω=ωc±(2π/N)的地方，H(ω)出现最大的肩峰值，肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少，则取决于旁瓣的多少。 （3）改变N，只能改变窗谱函数的主瓣宽度，改变ω的坐标比例以及改变WR(ω)的绝对值大小。例如，在矩形窗情况下， 当截取长度N增加时，只会减小过渡带宽度（4π/N），但不能改变主瓣与旁瓣幅值的相对比例; 同样，也不会改变肩峰的相对值。这个相对比例是由窗函数形状决定的，与N无关。换句话说，增加截取窗函数的长度N只能相应的减少过渡带，而不能改变肩峰值。 五种常用的窗函数 各种窗函数的傅里叶变换（N=51），A=20 lg|W(ω)/W(θ)| (a) 矩形窗 (b) 巴特利特窗（三角形窗） (c) 汉宁窗 (d) 海明窗 (e) 布拉克曼窗 理想低通滤波器加窗后的幅度响应（N=51）, A=20lg|H(ω)/H(0)| (a) 矩形窗 (b) 巴特利特窗（三角形窗） (c) 汉宁窗 (d) 海明窗 (e) 布拉克曼窗 凯塞（Kaiser）窗 窗函数法的设计步骤归纳 （1） 给定希望逼近的频率响应函数Hd(ejω)。 （2） 根据下式求单位脉冲响应hd(n)。 如果Hd(ejω)很复杂或不能直接计算积分，则必须用求和代替积分，也就是要计算离散傅里叶反变换， 一般都采用FFT来计算。将积分限分成M段，也就是令采样频率为ωk=2πk/M，k=0, 1, 2, …, M-1，则有 由于hd(n)有可能是无限长的序列，因此严格说，必须当M→∞时，hＭ(n)才能等于hd(n)而不产生混叠现象，即 。实际上，由于hd(n)随n的增加衰减很快，一般只要M足够大，即MN，近似就足够了。 （3） 按照过渡带及阻带衰减情况，选择窗函数形式。原则是在保证阻带衰减满足要求的情况下， 尽量选择主瓣窄的窗函数。 （4） 最后，计算所设计的FIR滤波器的单位脉冲响应。 h(n)=hd(n)w(n) 0≤n≤N-1 （5）由h(n)求FIR滤波器的系统函数H(z) 例 用矩形窗设计一个线性相位带通滤波器 -ωc≤ω-ω0≤ωc 0≤ωω0-ωc, ω0+ωcω≤π （1） 设计N为奇数时的h(n)。 （2） 设计N为偶数时的h(n)。 （