切比雪夫逼近法设计FIR滤波器

概念

切比雪夫逼近法，是在所需要的区间[a,b]内，使误差函数E(x)=|p(x)-f(x)|较均匀一致，并且通过合理选择p(x)，使E(x)的最大值En达到最小。切比雪夫最佳一致逼近的基本思想是，对于给定区间[a，b]上的连续函数f(x)，在所有n次多项式的集合中，寻找一个多项式，使其在[a，b]上对f(x)的偏差和其他一切属于集合的多项式对f(x)的偏差相比是最小的。 En=max|p(x)-f(x)| E(x)=p(x)-f(x) p(x)是f(x)最佳一致逼近多项式的充要条件是E(x)在[a，b]上至少存在n+2个交错点，使得E(xi)=En或者-En，同时要满足E(xi)= -E(xi+1)。 该n+2个点为交错点组，也为E(x)的极值点。 利用切比雪夫逼近理论设计FIR数字滤波器时，满足最佳一致逼近的滤波器又具有等纹波性质。存在通带纹波以及阻带纹波。 在实际求解过程中，利用数值分析中的Remez算法，靠一次次迭代来求得一组交错点组。

设计举例

%-------------------------------------------------------------------------- % FIR滤波器设计要求 % 切比雪夫逼近法：通带边缘频率0.6pi，阻带边缘频率0.7pi % to test remez.m and to design Low-pass FIR filter; %------------------------------------------------------------------------- clear; close; clc; f=[0 .6 .7 1]; % 给定频率轴分点； A=[1 1 0 0]; % 给定在这些频率分点上理想的幅频响应； weigh=[1 10]; % 给定在这些频率分点上的加权； b=remez(32,f,A,weigh); % 设计出切比雪夫最佳一致逼近滤波器； % [h,w]=freqz(b,1,256,1); h=abs(h); h=20*log10(h); figure(1) stem(b,'.');grid; figure(2) plot(w,h);grid; 左图是单位冲击响应，右图为其频响。

---

%-------------------------------------------------------------------------- % FIR滤波器设计要求 % 切比雪夫逼近法：Fs=500hz, 陷波 50 100 150hz % to test remez.m and to design multi-band FIR filter; %------------------------------------------------------------------------- clear; close; clc; % 用切比雪夫最佳一致逼近设计线性相位多带FIR滤波器； f=[0 .14 .18 .22 .26 .34 .38 .42 .46 .54 .58 .62 .66 1]; A=[1 1 0 0 1 1 0 0 1 1 0 0 1 1]; weigh=[8 1 8 1 8 1 8]; b=remez(64,f,A,weigh); % [h,w]=freqz(b,1,256,1); hr=abs(h); h=abs(h); h=20*log10(h); figure(1) stem(b,'.');grid; figure(2) plot(w,h);grid; 左图是单位冲击响应，右图为其频响。

滤波器设计总结

下面列出胡广书版的信号处理对滤波器设计的总结：IIR可取到非常好的通带和阻带的衰减，得到精准的通带阻带的边缘频率，且计算量小。缺点是不具有线性相位，且存在稳定性的问题。如果我们强调的是最大限度地去除噪声而没有别的限制，最佳选择是IIR滤波器。FIR优点是可以取得线性相位，且没有稳定性的问题。如果不要求实时实现，可实现零相位滤波。缺点是计算量大，不易实时实现。另外，切比雪夫逼近法是公认的最佳FIR滤波器的设计方法。梳状滤波器是针对信号中含有周期的噪声，或噪声中有周期的信号。除了传统意义上的滤波器外，还有广泛用于多抽样率信号处理的滤波器组、基于统计意义的最优现代滤波器（维纳滤波及自适应滤波等）等各种滤波器理论和设计方法。