信号与信号处理实验九FIR数字滤波器的设计.pptVIP

实验九 FIR数字滤波器的设计 一、 窗函数设计法 例：用窗函数法设计线性相位FIR低通数字滤波器，通带截止频率Wp=0.5π，阻带截止频率Ws=0.6π，阻带衰减不小于40dB，通带衰减不大于3dB. 根据阻带衰减，采用汉宁窗设计滤波器。MATLAB程序： Wp=0.5*pi; Ws=0.6*pi; wdel=Ws-Wp;%过度带宽度 N=ceil(6.2*pi/wdel);%计算滤波器长度 Wn=(Wp+Ws)/2;%截止频率 window=hanning(N+1); b=fir1(N,Wn/pi,window); freqz(b,1,512) figure；n=0:N;stem(n,b) * 1.根据滤波器技术要求直接设计一个理想数字滤波器的频率响应 理想频率响应的滤波器具有矩形频率特性，因而其时域响应为无限长的 （一）设计原理 2.求具有上述频率响应的滤波器对应的冲激响应序列 3.将求得的冲激响应序列截短，并使之具有对称性 以低通滤波器为例 1.理想数字滤波器 表示其群时延 2.冲激响应序列 3. 截短 相当于 原序列 加窗(乘窗函数） FIR冲激响应 FIR的频率响应 1. 主瓣宽度窄，过渡带陡。 窗函数要求： 2. 旁瓣小，改善通带平稳度和增大阻带衰减。 (β＝7.865) -80 5 -57 凯泽窗 -74 5.5 6 -57 不拉克曼窗 -53 3.3 4 -41 海明窗 -44 3.1 4 -31 汉宁窗 -25 2.1 4 -25 三角形窗 -21 0.9 2 －13 矩形窗 阻带最小衰减/dB 过渡带宽 /(2π/N) 主瓣宽度 /(2π/N) 旁瓣峰值 /dB 加窗后滤波器性能指标 窗谱性能指标 窗函数 几种窗函数的主要性能： MATLAB工具箱提供了采用窗函数法设计FIR滤波器的函数fir1 其调用格式为： b=fir1(n,Wn) b=fir1(n,Wn,’ftype’) b=fir1(n,Wn,window) b=fir1(n,Wn,’ftype’,window) 说明：n表示滤波器的阶数，‘ftype’ 参数代表所设计的滤波器类型：‘high’——高通， ‘stop’——带阻滤波器，默认时为低通和带通滤波器；‘window’为窗函数，是长度为n+1的列向量，默认时窗函数自动取Hamming窗。 clear all; Wp=0.2*pi; Ws=0.4*pi; tr_width=Ws-Wp; %过渡带宽度 N=ceil(6.6*pi/tr_width)+1 %滤波器长度（选择海明窗） Wc=(Ws+Wp)/2; %理想低通滤波器的截止频率 h=fir1(N,Wc/pi,hamming(N+1)) subplot(211) stem(h) [H,w]=freqz(h,1,512); mag=abs(H); db=20*log10((mag+eps)/max(mag)); %相对幅值响应 subplot(212) plot(w/pi,db) *