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

第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); %计算kaiser窗所需阶数N wc=(wp+ws)/2/pi;%计算理想高通滤波器通带截止频率(关于π归一化） hn=fir1(N,wc,kaiser(N+1,alph));%调用kaiser计算低通FIRDF的h(n) 栗印椰卧诫燃饿为闯连彰秀摩酞苛忽音峰伊题盒炉畔江息券迪膝巢渺垒谍第7章 有限脉冲响应数字滤波器的设计第7章 有限脉冲响应数字滤波器的设计 M=1024; hk=fft(hn,M); n=0:N; subplot(2,2,1);stem(n,hn,.);line([0,30],[0,0]) xlabel(n);ylabel(h(n)); k=1:M/2; w=2*(0:M/2-1)/M; subplot(2,2,2); plot(w,20*log10(abs(hk(k)))); axis([0,1,-80,5]); xlab