7第七章-FIR数字滤波器的设计.pptVIP

无穷级数 第七章FIR数字滤波器的设计方法 一、 为何设计FIR数字滤波器 2.应用需求 3.FIR DF 优点 二、 线性相位FIR数字滤波器 2.频率响应 3. H(ejω)线性相位 4.FIR滤波器具有线性相位的条件 5. 线性相位FIR滤波器的类型 三、FIR DF 设计 一窗函数设计方法 2.加窗影响——吉布斯效应 3. 几种典型窗函数 4.改善函数形状标准 矩形窗设计法的基本步骤： Example 1 试用矩形窗设计一个FIR数字低通滤波器，其通带截止频率ωc =0.3π，矩形窗的长度分别取N=41 和 N=121。 （2）低通滤波器的冲激响应: h(n)=hd(n-20) n=0,1,…,40 （3）N = 121 时 Program: % 低通滤波器设计 N=41; a=(N-1)/2 ; n=0: (N-1); h=sin(0.3*pi*(n-a))./(pi*(n-a+eps)); h(a+1)=0.3; %序号必须从1开始 subplot(1,2,1); plot(n,h); %画出脉冲响应 title(Impulse Response,N=41); axis([0,130,-0.1,0.4]); xlabel(time sample n); ylabel(h(n)); X=fft(h,128); subplot(1,2,2); plot(abs(X)); %画幅频特性 title(Magnitude Response,N=41); axis([0,60,0,1.1]); xlabel(ω in units of ?/N ); ylabel(H(ω)); grid 用窗函数设计FIR滤波器的步骤： （1）按性能指标要求给定理想频率特性 （2）求出 （3）由过渡带宽和阻带最小衰减的要求， 确定窗的形状和阶数N（试探多次） 一般 ， A决定于窗口 形状。 （4）求得FIR滤波器的单位脉冲响应 （5）求 验证结果是否满足设计要求。 Example 2 用哈明窗设计一个数字高通滤波器， 选取 N=81，截止频率 ωc=0.3π。 解： 哈明窗函数 理想脉冲响应滤波器的冲激响应 数字高通滤波器的脉冲响应 Example 5 根据下列技术指标，设计一个FIR数字低通滤 波器。 ωp=0.2? ，Ap=0.25dB， ωs=0.3? ，As=50dB。 查表7-3，根据过渡带宽的要求可计算N，得 N=6.6 ?/ Δω，其中Δω为过渡带的宽度。 程序： %根据技术指标，设计FIR数字低通滤波器 %选哈明窗 wp=0.2*pi;ws=0.3*pi; deltaw=ws-wp; %Δω为过渡带的宽度 N1=ceil(6.6*pi/deltaw); %计算N N=N1+mod(N1+1,2); %确保N的长度为奇数,实现第一类线性相位 滤波 wdham=(hamming(N)); %求哈明窗 wc=(wp+ws)/2; %截止频率取两边频率的平均值 hd=ideallp(wc,N); h=hd.*wdham; %检验 H=fft(h,128); HdB=20*log10(abs(H+eps)/max(abs(H))); %计算相对振幅 Ap=-min(HdB(1:wp/(2*pi/128)+1)); %检验通带波动 plot(HdB); legend(哈明窗,FIR数字低通滤波幅频特性) 二频率抽样设计法 2. 线性相位的约束 3. 逼近误差及其改进措施 三最优化设计法 最优化设计的matlab函数 2. cremez() 1、从性能上进行比较 2、从结构上看 3、从设计工具看 另外 总结 误差的大小与理想频率响应曲线的形状有关，越平滑的特性部分误差越小，内插值越接近理想值；在特性曲线间断点附近，插值与理想值的误差最大，并产生尖峰和波纹，形成振荡。振荡的形成会减小阻带衰减，破坏滤波器性能。 改进措施： （1）增加N，减少逼近误差，但间断点 附近的误差仍很大； （2）在频响间断点附近区间增加一个或 多个过渡采样点，使不连续点变成 缓慢过渡。结果是加大了过渡带， 但增大了阻带衰减。 将所有抽样值作为变量，选取某一函数作为逼近模型，在某一优化准则下，通过迭代运算，以得到最优结