电磁波大学课件第7章-2.pptVIP

7.2窗函数设计方法 *理论学习及滤波器的设计 FIR 滤波器的EDA设计 7.2.1 理想低通滤波器 理想低通的单位脉冲响应为： 是无限长、非因果的。 用窗函数法设计FIR滤波器的思路。 7.2.2 时域加窗 实际的滤波器长度为N，故： 有限项近似必然产生截断效应。  设 则由复卷积定理得： 其中矩形窗的频谱为： 理想低通滤波器的频谱为： 其幅度特性为： 冲激响应截短后的低通滤波器频谱： 幅度特性等于理想低通滤波器的幅度 特性与矩形窗幅度特性的卷积。 卷积结果 矩形窗对理想低通幅度特性的影响 对 加矩形窗处理后 其频率响应将产生以下几点影响： （1）在理想特性不连续点 附近形成过 渡带。过渡带宽近似等于 主瓣宽 度 ； （2）在通带和阻带内出现波动，并在截止频率 的两边 出现最大尖峰值； （3）主瓣附近窗的频率响应为： 随着N的加大，振荡变密，主瓣变窄； 主瓣与旁瓣的幅度亦有所加大，但主 瓣与旁瓣的相对比例不变（吉布斯现象?）。 7.2.3 矩形窗 （Rectangular Window） 矩形窗函数： 其频率响应为： 矩形窗设计法的基本步骤： ?1 确定窗口的宽度N； ?2 对理想低通滤波器的冲激响应进 行截断； ?3 进行延时，变为因果序列。 （ 矩形窗频谱的主瓣宽度为4?/N, 第一旁瓣比主瓣低13dB。） Example 1 试用矩形窗设计一个FIR数字低通滤波器，其通带截止频率ωc =0.3π，矩形窗的长度分别取N=41 和 N=121。 Solution: （1）N=41时，有 a=(N-1)/2=20 （2）低通滤波器的冲激响应: h(n)=hd(n-20) n=0,1,…,40 h(20)=hd (0)= ωc/π=0.3 % n-a=0 （3）N = 121 时 a=(N-1)/2=60, h(60)=had(0)=0.3 其低通滤波器的冲激响应为： h(n)=had(n-60) n=0,1,…,120 fir_lp.m % 低通滤波器设计 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); %画出脉冲响应 fir_lp.m title(Impulse Response,N=41); xlabel(time sample n); ylabel(h(n)); X=fft(h,128); subplot(1,2,2); plot((2*(-64:63))/128, 10*log10(abs(fftshift(X)))); %画幅频特性 title(Magnitude Response,N=41); xlabel(ω in units of\pi ); ylabel(H(ω)); Grid on 7.2.4 其它常用窗函数（wintool） （1）三角窗（Bartlett Window）： *用窗函数设计FIR滤波器的步骤： （1）按性能指标要求给定希望的逼近的频率响应函数Hd(ejω)。 （2）求出 （3）由过渡带宽和阻带最小衰减的要求，确定窗的形状和阶数N（试探多次）一般 ， A决定于窗口形状。 表7.2.2 六种窗函数的基本参数 7.2.5 用窗函数法设计其它通带滤波器 （1）FIR数字高通滤波器的设计 理想高通频响特性： 其中 单位脉冲响应为 一般在 处， 则不能做为高通滤波器，故N只能为奇数。 FIR数字高通单位脉冲响应 （2）FIR数字带通滤波器设计 理想带通的频率响应 （3）FIR数字带阻滤波器设计 理想带阻频率特性： FIR数字带阻： Example 3 用哈明窗设计一个数字高通滤波器， 选取 N=81，截止频率 ωc=0.3π。 解： 哈明窗函数 理想脉冲响应滤波器的冲激响应 Example 根据下列技术指标，设计一个FIR数字低