数字信号处理第4章有限长单位脉冲响应(FIR)滤波器的设计方法分析报告.ppt

N=65;k=0:(N-1)/2; Wm=2*pi*k./N; Ad(1:(N+1)/2)=1; Ad(18)=0.5886;Ad(19)=0.1065;Ad(20:33)=0; Hd=Ad.*exp(-j*0.5*(N-1)*Wm); Hd=[Hd conj(fliplr( Hd(2:(N+1)/2) ) )]; h=real(ifft(Hd)); w=linspace(0,pi-0.1,1000); H=freqz(h,[1],w); plot(w/pi,20*log10(abs(H)));grid; 小结： 频率采样设计法优点： ①???直接从频域进行设计，物理概念清楚，直观方便； ②?适合于窄带滤波器设计，这时频率响应只有少数几个非零值。 缺点：截止频率难以控制。 因频率取样点都局限在2π/N的整数倍点上，所以在指定通带和阻带截止频率时，这种方法受到限制，比较死板。充分加大N，可以接近任何 给定的频率，但计算量和复杂性增加。 4.4 FIR数字滤波器的最优化设计 最优化设计的前提是最优准则的确定，在FIR滤波器最优化设计中，常用的准则有 ①均方误差最小准则 ②最大误差最小化准则。 1)??均方误差最小化准则， 若以E(ejω)表示逼近误差，则 均方误差为  根据Parseval定理 所以，矩形窗窗口设计法是一个最小均方误差FIR设计，根据前面的讨论，我们知道其优点是过渡带较窄，缺点是局部点误差大,或者说误差分布不均匀。 2)??最大误差最小化准则（也叫最佳一致逼近准则） 其中F是根据要求预先给定的一个频率取值范围，可以是通带，也可以是阻带。最佳一致逼近即选择N个频率采样值 （ 或时域 h(n) 值 ），在给定频带范围内使频响的最大逼近误差达到最小。也叫等波纹逼近。 优点：可保证局部频率点的性能也是最优的，误差分布均匀，相同指标下，可用最少的阶数达到最佳化。  例如，我们提到的频率采样最优化设计，它是从已知的采样点数N、预定的一组频率取样和已知的一组可变的频率取样（即过渡带取样）出发，利用迭代法（或解析法）得到具有最小的阻带最大逼近误差（即最大的阻带最小衰减）的FIR滤波器。但它只是通过改变过渡带的一个或几个采样值来调整滤波器特性。如果所有频率采样值（或FIR时域序列h(m)）都可调整，显然，滤波器的性能可得到进一步提高。 4.4.1 非线性最优化 为了简化分析，假定滤波器的单位脉冲响应为对称于n=0的单位脉冲响应，且N为奇数。 令N=2M+1 低通滤波器的误差分配  切比雪夫最佳一致逼近 如图，用等波纹逼近法设计滤波器需要确定五个参数： M、ωc、ωr、δ1、δ2 按上图所示的误差容限设计低通滤波器，就是说要在通带 0? ω?ωp 内以最大误差 δ1 逼近1，在阻带ωr ? ω?? 内 以最大误差δ2逼近零。 要同时确定上述五个参数较困难。常用的两种逼近方法： 1）给定M、δ1、δ2，以ωc和ωr为变量。 缺点：边界频率不能精确确定。 2）给定M、ωc和ωr，以δ1和δ2为变量，通过迭代运算 ，使逼近误差δ1和δ2 最小，并确定h(n)——切比雪 夫最佳一致逼近。 特点：能准确地指定通带和阻带边界频率。?  等波动逼近的低通滤波器 ?c ?r  一.误差函数 定义逼近误差函数： 为所设计的滤波器与理想滤波器的幅频特性在通带和阻带内的误差值， 是已知的权函数，在不同频带可取不同的值， 所要设计的滤波器的幅频特性 理想滤波器的幅频特性 例如，希望在固定 M, ?c, ?r 的情况下逼近一个低通滤波器，这时有 对于表4.1中的第一种滤波器， 于是 切比雪夫逼近问题变为，寻求一组系数 使逼近误差的最大值达到最小，即 给定后等效于求 最小。  二.交替定理（最佳逼近定理） 令F表示闭区间 的任意闭子集，为了使 在 F 上唯一最佳地逼近于 ,其充分必要条件是误差函数