第7章 有限脉冲响应数字滤波器的设计.pptVIP

第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.4 利用等波纹最佳逼近法设计FIR数字滤波器 7.5 IIR和FIR数字滤波器的比较 有限脉冲响应(FIR)滤波器在保证幅度特性满足技术指标的同时，很容易做到有严格的线性相位特性。 用N表示FIR滤波器单位脉冲响应h(n)的长度，其系统函数H(z)为 H(z)是z－1的N－1次多项式， 有N－1个零点，在原点z=0处有一个N－1重极点。 因此，H(z)永远稳定。 稳定和线性相位特性是FIR滤波器最突出的优点。 FIR滤波器设计任务： 选择有限长度的h(n)，使频率响应函数H(ejω)满足技术指标要求。 7.1 线性相位FIR数字滤波器的条件和特点 1. 线性相位FIR数字滤波器 对于长度为N的h(n) ，传输函数为 H(ejω)线性相位是指θ(ω)是ω的线性函数，即 θ(ω)= ─ τω, τ为常数 (7.1.3) 如果θ(ω)满足下式： θ(ω)=θ0─τω, θ0是起始相位 (7.1.4) 严格地说，此时θ(ω)不具有线性相位，但以上两种情况都满足群时延（相位特性曲线的斜率）是一个常数，即 也称这种情况为线性相位。 满足(7.1.3)式是第一类线性相位； 满足(7.1.4)式是第二类线性相位。 θ0=－π/2是第二类线性相位特性常用的情况。 4. 线性相位FIR滤波器网络结构 设N为偶数，则有 7.5 IIR和FIR数字滤波器的比较 IIR与FIR滤波器各有所长，在实际应用中应该考虑多方面因素加以选择，同时也要兼顾经济方面的要求以及计算工具的实际条件。 【例7.2.2】 对模拟信号进行低通滤波处理，要求通带0≤f ≤1.5kHz内衰减小于1 dB，阻带2.5kHz≤f ≤∞上衰减大于40 dB。希望对模拟信号采样后用线性相位FIR数字滤波器实现上述滤波，采样频率Fs=10 kHz。用窗函数法设计满足要求的FIR数字低通滤波器，求出h(n)，并画出损耗函数曲线。为了降低运算量，希望滤波器阶数尽量低。 解： （1） 确定相应的数字滤波器指标: 通带截止频率为 脊岁厉名闪下雌鄙生沧定肥托盂邦蔫想嘛撼惹盏豪怖柏垛秉跋吹撤濒榜侦第7章 有限脉冲响应数字滤波器的设计第7章 有限脉冲响应数字滤波器的设计 阻带截止频率为 阻带最小衰减为 ? s =40 dB （2） 用窗函数法设计FIR数字低通滤波器，为了降低阶数选择凯塞窗。根据式（7.2.16）计算凯塞窗的控制参数为  绚瓶植匪捷声叙壳慕吕嫌徒硕微躲卡惶茫相牲馁让野婪役夸抠苯桓射涟惠第7章 有限脉冲响应数字滤波器的设计第7章 有限脉冲响应数字滤波器的设计 指标要求过渡带宽度Bt=ωs－ωp=0.2π，根据式（7.2.17）计算滤波器阶数为   取满足要求的最小整数M=23。所以h(n)长度为N=M+1=24。 理想低通滤波器的通带截止频率ωc=(ωs+ωp)/2=0.4π，所以由式（7.2.2）和式（7.2.3）, 得到:   式中，w(n)是长度为24（ ? =3.395）的凯塞窗函数。 赫柿韭娠氏斜耘乱幸对监肃脸鸵摈爬箭醇魔有滓掺橡旷犀增版甲引辕全湘第7章 有限脉冲响应数字滤波器的设计第7章 有限脉冲响应数字滤波器的设计 　　实现本例设计的MATLAB程序： 　　％ 用凯塞窗函数设计线性相位低通FIR数字滤波器fp=1500;fs=2500;rs=40; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; Bt=ws-wp; %计算过渡带宽度 alph=0.5842*(rs-21)^0.4+0.07886*(rs-21); %根据(7.2.16)式计算kaiser窗的控制参数α N=ceil((rs-8)/2.285/Bt); %计算k