北京邮电大学数字信号处理习题答案第5章

《北京邮电大学数字信号处理习题答案第5章》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、FIR数字滤波器设计本章知识点：对于一个离散时间系统,若分母多项式中系数，则此系统就变成一个FIR系统，其中系数即为该系统的单位取样响应h(0),h(1),…h(N-1),且当n>N-1时，h(n)=0。FIR系统函数H(z)在Z平面上有N-1个零点，在原点z=0处有N-1个重极点。这类系统不容易取得较好的通带和阻带特性，要想得到与IIR系统类似的衰减特性，则要求较高的H(z)阶次。相比于IIR系统来说，FIR系统主要有三大突出优点：1）系统永远稳定；2）易于实现线性相位系统；3）易于实现多通带（或多组带）系统。线性相位FIR滤波器实现的充要条件是：对于任意给定的数值N

2、（奇数或偶数），冲激响应h[n]相对其中心轴必须成偶对称或奇对称，此时滤波器的相位特性是线性的，且群延时均为常数。由于h(n)有奇对称和偶对称两种情况，h(n)的点数N有奇数、偶数之分。因此，h（n）可以有4种不同的类型，分别对应于4种线性相位FIR数字滤波器：h[n]偶对称N为奇数、h[n]偶对称N为偶数、h[n]奇对称N为奇数、h[n]奇对称N为偶数。四种线性相位FIR滤波器的特性归纳对比于表5.１中。表5.1线性相位FIR滤波器特性偶对称h(n)=h(N-1-n)奇对称h(n)=-h(N-1-n)相位函数ωθ(ω)ωθ(ω)相N为奇数n=1,2,…,(N-1)/2

3、相N为偶数b(n)=，n=1,2,…N/2n=1,2,…,N/2一．FIRDF设计方法FIRDF的设计实现不能像IIRDF设计那样借助于模拟滤波器的设计方法来实现，其设计方法主要是建立在对理想滤波器频率特性进行不同程度逼近的基础上，主要的逼近方法有三种：窗函数法；频率抽样法；最佳一致逼近法。1．窗函数法窗函数法是设计FIR滤波器的最直接方法，它通过采用不同时宽的窗函数，对理想滤波器的无限长冲激响应hd(n)进行截短，从而得到系统的有限长冲激响应h(n)，这一过程可用式5－1来描述：（5.1）其中WR(n)是时宽为N的窗函数。由分析可知，加窗处理使所得滤波器的频响与理想滤

4、波器频响之间产生差异，具体表现在过渡带和波动的出现。我们希望所设计的滤波器尽量逼近理想滤波器，就要设法减小波动的幅值，同时使过渡带变窄。因此为了改善滤波器的性能，要求窗函数尽量具有以下特性：1）主瓣宽度尽可能地窄，以获得尽量陡的过渡带。2）最大旁瓣相对于主瓣尽可能地小，即能量尽可能集中于主瓣内，以使肩峰和波动减小。对于窗函数，以上两个要求是相互矛盾的，不可能同时达到最佳，要根据需要进行折衷的选择，通常是以增加主瓣宽度来换取对旁瓣的抑制。下面是几种常用的窗函数及其频谱特性：¨矩形窗（5.2）(5.3)¨三角窗（或巴特利特Bartlett窗）(5.4)(5.5)¨汉宁（Ha

5、nning）窗——升余弦窗(5.6)（5.7）¨汉明（Hamming）窗——改进的升余弦窗(5.8)(5.9)¨布莱克曼（Blackman）窗——又称二阶升余弦窗(5.10)(5.11)¨凯瑟（Kaiser）窗(5.12)表5.2列出了以上各窗函数的综合性能指标。表5.3给出了不同β值下的凯瑟窗性能总结。表5.2常用窗函数的性能指标窗函数旁瓣峰值衰减（dB）窗函数主瓣宽度加窗后滤波器过渡带宽（△ω）加窗后滤波器阻带最小衰减（dB）矩形窗-134π/N1.8π/N-21汉宁窗(升余弦窗)-318π/N6.2π/N-44汉明窗(改进升余弦窗)-418π/N6.6π/N-53

6、布莱克曼窗(二阶升余弦窗)-5712π/N11π/N-74凯塞窗（β=7.865）-5710π/N10π/N-80表5.3不同β值下的凯瑟窗性能β加窗后滤波器过渡带（△ω）加窗后滤波器阻带最小衰减（dB）2.1203.00π/N-303.3844.46π/N-404.5385.86π/N-505.6587.24π/N-606.7648.64π/N-707.86510.0π/N-808.96011.4π/N-9010.05612.8π/N-100用窗函数法设计FIR滤波器的实现步骤如下：(1)根据所要设计滤波器的性能指标（阻带最小衰减、过渡带宽），通过查表来选定窗函数w(

7、n)，通过计算求得窗的宽度N：根据所要设计线性相位FIR滤波器的类型来决定最终N取奇数或偶数，一般情况下取奇数。(2)根据相应的频率响应函数，求傅立叶反变换，得到hd(n)。(3)按所得窗函数求出FIR滤波器的冲激响应：h(n)=hd(n)w(n)n=0,1,2,…,N–1。(4)利用h(n)计算FIR滤波器的频响，并检验各项指标，如不符合指标，则重新修改N及w(n)。2．频率取样法窗函数法是在时域内，以有限长冲激响应h(n)去近似所要求的理想冲激响应hd(n)，从而实现FIR滤波器的设计。而频率取样法则是在频域内，以有限个频率响应取样，