第FIR数字滤波器设计.ppt

* 图 7-17 内插后矩形理想特性与梯形特性频率采样的响应 由图可见： 误差大小取决于理想频率响应曲线形状。 理想频率响应特性变化越平缓，则内插（采样）值越接 近理想值，逼近误差越小。 * 7.4.4 过渡带采样的最优设计 为了提高逼近质量，在理想特性不连续点处人为加入过 渡采样点（1～3个），虽然加宽了过渡带，但缓和了边缘上 两采样点之间的突变，将有效的减少起伏振荡，提高阻带衰 减。 过渡点的值，可由计算机通过线性最优化技术来确定。 例7-8 利用频率采样法设计一个线性相位FIR低通滤波 器，已知 （1）采样点数N=33，wc=p/2； （2）采样点数N=34，wc=p/2 ，设置两个过渡点 |H1(k)|=0.5925， |H2(k)|=0.1099。 解：(1)首先选择滤波器的种类。由于要设计的是低 通，且N为奇数，故选择1型滤波器，于是有 由 ，确定通带内的采样点数 取整数k=8，应在通带内设置9个采样点（k=0～8），第 10个采样点已在通带截止频率之外，处于阻带内。 * 根据Hk=HN-k，可得Hk为 因此 将H(k)代入 内插公式即得所设计滤波器的频率响 应。 (2)选择滤波器的种类。由于要设计的是低通，且N为 偶数，故选择第2类滤波器，于是有 由 ，确定通带内的采样点数 取整数k=8，应在通带内设置9个采样点（k=0～8），第 10个采样点已在通带截止频率之外，可将第10、11点设为过 渡点。 根据Hk=-HN-k，可得Hk为 Hk=1(k=0～8)，0.5925(k=9)，0.1099(k=10), -0.1099(k=24), -0.5925(k=25)，-1(k=26～33)，0(其它) * * 因此 将H(k)代入 内插公式即得所设计滤波器的频率响 应。 * 例7-9 要求编程实现例题7-8中所要求的内容。 解：Matlab主要程序 %过渡带不设置采样点 N=33;alpha=(N-1)/2;k=0:N-1;wk=(2*pi/N)*k; Hk=[ones(1,9),zeros(1,16),ones(1,8)]; Hd=[1,1,0,0];wd=[0,0.5,0.5,1]; angH=-alpha*(2*pi)/N*k;H=Hk.*exp(i*angH); h=real(ifft(H,N)); [db,mag,pha,grd,w]=freqz_m(h,1); [Hr,ww,a,L]=Hr_Type1(h); subplot(2,2,1);plot(ww/pi,Hr,wk(1:17)/pi,Hk(1:17),o,wd,Hd); axis([0,1,-0.5,1.2]); subplot(2,2,2);plot(w/pi,db);axis([0,1,-100,10]);grid 仿真曲线 * 7.5 FIR滤波器和IIR滤波器的比较 前面讨论了IIR和FIR两种滤波器传输函数的设计方法。 这两种滤波器究竟各自有什么特点？在实际运用时应该怎样 去选择它们呢？ 为此对这两种滤波器作一简单的比较。 (1).从性能上比较 (2).从结构上比较 (3).从设计工具上比较 * (1).从性能上进行比较 IIR滤波器 传输函数的极点可位于单位圆内的任何地方，因此可用 较低的阶数获得高的选择性，所用的存贮单元少，所以经济 而效率高。但是这个高效率是以相位的非线性为代价的。选 择性越好，则相位非线性越严重。 FIR滤波器 可以得到严格的线性相位，然而由于FIR滤波器传输函 数的极点固定在原点，所以只能用较高的阶数达到高的选择 性；对于同样的滤波器设计指标，FIR滤波器所要求的阶数 可以比IIR滤波器高5-10倍，所以，成本较高，信号延时也 较大。 * (2).从结构上进行比较 IIR滤波器 必须采用递归结构，极点位置必须在单位圆内，否则系 统将不稳定。另外，在这种结构中，由于运算过程中对序列 的舍入处理，这种有限字长效应有时会引入寄生振荡。 FIR滤波器 主要采用非递归结构，不论在理论上还是在实际的有限 精度运算中都不存在稳定性问题，运算误差也较小。此外， FIR滤波器可以采用快速傅里叶变换算法，在相同阶数的条 件下，运算速度可以快得多。 * (3).从设计工具进行比较 IIR滤波器 可以借助于模拟滤波器的成