第五章 IIR滤波器 设计 - 重庆大学电气工程学院.ppt

 3）增大N 如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数N。 例如，同样边界频率ωc=0.5π, 以N=65采样，并在k=17和k=48插入由阻带衰减最优化计算得到的采样值H17=H48=0.5886,在k=18、47处插入经阻带衰减最优化计算获得的采样值H17=H48=0.1065 , 这时得到的 H(ejω)，过渡带为6π/65，而阻带衰减增加了20多分贝，达-60dB以上，当然，代价是滤波器阶数增加，运算量增加。 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; axis([0 1 -100 20]); xlabel(\pi); title(过渡带设两个采样点，总采样点数提高一倍); 六、IIR和FIR滤波器的比较 ? FIR IIR 设计方法 一般无解析的设计公式，要借助计算机程序完成 利用AF的成果，可简单、有效地完成设计 ? 设计结果 可得到幅频特性（可以多带）和线性相位（最大优点） 只能得到幅频特性，相频特性未知（一大缺点），如需要线性相位，须用全通网络校准，但增加滤波器阶数和复杂性 稳定性 极点全部在原点（永远稳定）无稳定性问题 有稳定性问题 ? 阶 数 高 ? 结 构 非递归 递归系统 运算误差 一般无反馈，运算误差小 有反馈，误差大 快速算法 可用FFT实现，减少运算量 无快速运算方法 低 从上面的简单比较我们可以看到IIR与FIR滤波器各有所长，所以在实际应用时应该从多方面考虑来加以选择。例如，从使用要求上来看，在对相位要求不敏感的场合，如语言通讯等，选用IIR较为合适，这样可以充分发挥其经济高效的特点，而对于图像信号处理，数据传输等以波形携带信息的系统，则对线性相位要求较高，如果有条件，采用FIR滤波器较好，当然，在实际应用中应考虑经济上的要求以及计算工具的条件等多方面的因素。 wn=kaiser(30,4.55); nn=[0:1:29]; alfa=(30-1)/2; hd=sin(0.4*pi*(nn-alfa))./(pi*(nn-alfa)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-80,10]); grid; xlabel(归一化频率/p) ylabel(幅度/dB) 利用矩形窗设计结果 利用凯塞窗设计结果 各窗函数的时域表示 窗函数 表达式W[n] ， -0.5(N-1)?n?0.5(N-1) Boxcar 1 Bartlett Hanning Hamming Blackman Kaiser 以下给出上面各窗函数的幅度函数 六种窗函数的基本参数 四、FIR窗函数法的设计步骤 1、根据技术要求确定逼近理想滤波器的单位冲激响应hd(n). 如果给出通带阻带衰减和边界频率的要求，可选用理想滤波器作为逼近函数，从而用理想滤波器的特性作付里叶反变换，求出hd(n)。 2、根据对过渡带及阻带衰减的要求，选择窗函数w(n)的形式和估计窗口长度N。 设待求滤波器的过渡带用?ω表示，它近似等于窗函数主瓣宽度。 N?A/? ω A决定窗口形式。 例如：矩形窗A=4?，汉宁窗A=8?等 按照过渡带及阻带衰减情况，根据P342页表7-3，选择窗函数形式。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。 3、对理想滤波器进行加窗，计算滤波器的单位冲激响应h(n)。 式中w(n)是上面选择好的窗函数。如果要求线性相位，则要求hd(n)和w(n)均对(N-1)/2对称。 4、验算技术指标是否满足要求。设计出的滤波器频率响应用下式计算： 计算上式时可用FFT算法。如果H(ejw)不满足要求，可根据具体情况重复2、3、4步，直到满足要求。 窗函数设计法的特点： 设计过程简单、方便实用。但边界频率不易精确控制。所以设计完以后，必须检验结果。 FIR DF设计的窗函数法不但可以用来设计普通的L