第六章 无限脉.pptVIP

线性相位FIR 的特点 窗函数法设计数字滤波器 频率取样法设计数字滤波器 FIR滤波器和IIR滤波器的比较 数字滤波器的应用 6.1线性相位FIR的特点 例如：FIR线性相位滤波器的h(n)是实数且n0和n6时，h(n)＝0，如果h(0)＝1且系统函数在z=3和各 有一个零点，求H(Z) 二、窗函数对频响的影响 二、窗函数对频响的影响 4） 时，主瓣全部在通带外， 出现负的肩峰。 6.3.1采样值H(k)具有线性相位的约束 性能上：IIR用较少的阶数获得较高的选择特性，非线性相位；FIR运算量大，线性相位； 结构上：IIR递归型结构，有舍入误差引起寄生振荡；FIR非递归型结构，不存在稳定性问题，误差小，可以采用FFT提高运算速度； 设计工作：IIR计算量小，FIR计算工具要求高，只有程序可循。 适用范围：IIR拘泥于模拟滤波器的格局；FIR灵活 N=61时，最小的阻带衰减为18dB，仍不满足指标。 结论：频率采样法设计FIR滤波器，不能用增大采样点数的方法提高阻带的衰减。 在Hd(ejω)的间断点附近区间加入若干个过渡采样点，能否使H(ejω)与 Hd(ejω)的误差减小？ 例如：增加两个采样点T1=0.7和T2=0.2。 N=61 两个过渡点 从幅度响应曲线可见，N=61时，最小的阻带衰减为45dB，达到设计要求。 结论：为改善滤波器的特性，可以在频响间断点附近插入一个或几个过渡采样点，适当增加采样点数提高阻带的衰减。 6.5 FIR滤波器和IIR滤波器的比较 例如：设计一个截止频率为 的数字低通滤波器 ★矩形窗过渡带最窄，布拉克曼窗过渡带最宽，但是对于阻带最小衰减性矩形窗最差，只有－21dB，布莱克曼窗最好，达到－75dB。汉明窗与汉明窗的过渡带稍宽，阻带衰减性也较好，是较为常用的窗函数。 改变N只能改变过渡带宽窄，不能改变阻带的衰减性； 阻带的衰减性由窗函数类型决定，矩形窗最差，布莱克曼窗最好; 设计数字滤波器时，一般希望滤波器既要有很好过渡性又要有很好的阻带衰减性。过渡带窄时，但阻带衰减不好；阻带衰减好时，但过渡带宽；因此，实际中窗函数的选择往往是它们的折衷。 三. 选择窗函数的原则 四 窗函数设计FIR滤波器的步骤 1.给定希望逼进的 2.确定待求滤波器的单位取样响应 3.根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。 4.计算滤波器的单位取样响应 5.确定H(Z) 例1．根据下列技术指标，设计一个数字FIR低通滤波器：ωp=0.15π,Ap=0.25dB，ωs=0.4π,As=50dB，选择一个恰当的窗函数，确定单位脉冲响应，绘出所设计的滤波器的幅度函数。 解：根据窗函数最小阻带衰减的特性，只有海明窗和布莱克曼窗可提供大于50dB的衰减。故选择海明窗，它提供较小的过渡带，因此具有较小的阶数。 ωc=(ωp+ωs)/2=(0.4π+0.15π)/2=0.275π 当N=34 时： α=(N-1)/2=(34-1)/2=16.5 Matlab 中FIR数字滤波器设计相关函数 　窗函数 ? ? ? ? ? ? ? ? ? ? 矩形窗 三角窗 宁汉窗 MATLAB 函数 MATLAB 函数 　窗函数 boxcar(N) triang(N) hanning(N) hamming(N) blackman(N) Kaiser(N,BETA) 汉明窗 布莱克曼 凯泽窗 1.窗函数 b = fir1(n,Wn)　　　　　　 b = fir1(n,Wn,ftype) b = fir1(n,Wn,window)　　　　 b = fir1(n,Wn,ftype,window) n滤波器的阶数：对于高通、带阻滤波器n为偶数； ftype滤波器类型：低通low、高通high、带通bandpass、带阻stop Wn截止频率：低通、高通时Wn，取值范围为（0.0~1.0），其中1对应0.5fs，fs为采样频率，带通、带阻时Wn ＝[W1　W2]且W1W2，默认为低通或高通。 window窗函数类型：window可取boxcar、hanning、blankman及kaiser等，缺省为hamming； b为FIR滤波器系数向量：为n+1。 2.FIR数字滤波器设计函数 例：用窗函数法设计FIR低通滤波器，通带截止频率Wp＝0.2π，阻带截止频率Ws＝0.3 π，阻带衰减不小于40dB，通带衰减不大于3dB wp=0.2*pi;ws=0.3*pi; wdelta=ws-wp; N=ceil(8*pi/wdelta); Wn=(0.2+0.3)*pi/2; b=fir1(N,Wn/pi,h