数字信号处理-有限脉冲响应数字滤波器设计.ppt

* 6.4 切比雪夫最佳逼近法 设希望设计的滤波器幅度特性为Hd(ω)，实际设计的滤波器幅度特性为Hg(ω)，其加权误差E(ω)用下式表示： E(ω)=W(ω)[Hd(ω)-Hg(ω)] 什么是最佳逼近？ 设计具有线性相位的FIR滤波器，使E(ω)最优 其单位脉冲响应序列h(n) 必须满足一定条件，例如：h(n)=h(n-N-1)且N为奇数 * M=(N-1)/2 最佳逼近问题： 选择M+1个系数a(n),使加权误差E(ω)的最大值为最小 W(ω)可以在不同频段选择不同权重 * 该定理指出最佳一致逼近的充要条件是E(ω)在A上至少呈现M+2个“交错”(alternations)，使得 Alternation Theorem 属于多项式逼近问题，证明已超出要求 * 对于低通滤波器 性能指标 Using Matlab 设计带阻滤波器，通带[0, 0.2pi]、[0.8pi, pi]，ap=1dB；阻带[0.35pi, 0.65pi]，as=60dB。 %BS_Equirip f=[0.2,0.35,0.65,0.8]; m=[1,0,1]; ap=1;as=60; delta1=(10^(ap/20)-1)/(10^(ap/20)+1); delta2=10^(-as/20); rip=[delta1,delta2,delta1]; [M,fo,mo,w]=remezord(f,m,rip); hn=remez(M,fo,mo,w); figure(1);stem(0:(length(hn)-1),hn); fk=0:1/200:1; hk=freqz(hn,1,fk*pi); figure(2);plot(fk,20*log10(abs(hk))); xlabel(Frequency(pi));ylabel(Magnitude(dB)); * M=28, N=29, fo=[0, .2, .35, .65, .8, 1], mo=[1, 1, 0, 0, 1, 1], w=[1, 57.5011,1] 设计数字低通滤波器, Ωp=1.5kHz, ap=1dB; Ωs=2.5kHz, as=40dB; 采样频率10kHz %LP_Equirip Fs=10000; f=[1500,2500]; m=[1,0]; ap=1;as=40; delta1=(10^(ap/20)-1)/(10^(ap/20)+1); delta2=10^(-as/20); rip=[delta1,delta2]; [M,fo,mo,w]=remezord(f,m,rip,Fs); M=M+1; hn=remez(M,fo,mo,w); figure(1);stem(0:(length(hn)-1),hn); fk=0:1/200:1; hk=freqz(hn,1,fk*pi); figure(2);plot(fk,20*log10(abs(hk))); xlabel(Frequency(pi));ylabel(Magnitude(dB)); * M=15, N=16, w=[1, 5.7501] * 6.5 IIR和FIR数字滤波器的比较 性能： IIR滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的 对于同样的性能指标，FIR滤波器所要求的阶数一般比IIR滤波器高5-10倍，成本高，信号时延大，但是为线性相位 结构： IIR滤波器必须采用递归结构，极点位置必须在单位圆内，否则系统将不稳定；另外，由于运算过程中的舍入处理，有限字长效应可能会引起寄生振荡 FIR滤波器采用非递归结构，在理论和实际中都不存在稳定性问题，运算误差引起的噪声较小；此外，FIR可以采用FFT算法实现，运算速度可以大大提高 * 设计工具： IIR滤波器可以借助于模拟滤波器的成果，因此一般都有有效的封闭形式的设计公式可供准确计算，计算工作量比较小，对计算工具的要求不高。 FIR滤波器计算通带、阻带衰减等无显式表达式，其边界频率不易精确控制；FIR的设计依赖计算工具 灵活性： IIR滤波器虽然设计简单，但主要用于低通、高通、带通和带阻等滤波器，脱离不了典型滤波器的频响特性 FIR滤波器则要灵活很多，易于适应某些特殊的应用 * 6.6 FDATool Matlab：FDATool * 图6.2.5 常用窗函数的幅度特性 (a)矩形窗；(b)巴特利特窗(三角形窗)；(c)汉宁窗； (d)哈明窗；(e)布莱克曼窗 * 图6