第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） 确定相应的数字滤波器指标: 通带截止频率为 阻带截止频率为 阻带最小衰减为 ? s =40 dB （2） 用窗函数法设计FIR数字低通滤波器，为了降低阶数选择凯塞窗。根据式（7.2.16）计算凯塞窗的控制参数为  指标要求过渡带宽度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）的凯塞窗函数。 　　实现本例设计的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) M=1024; hk=fft(hn,M); n=0:N; subplot(2,2,1);stem(n,hn,.);line([0,30