第7章FIR数字滤波器的设计剖析.ppt

第7章 有限脉冲响应数字滤波器(Finite Impulse Response Filters——FIR) 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.4 利用等波纹最佳逼近法设计FIR滤波器 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 四、线性相位FIR低通滤波器的设计 五、线性相位FIR高通滤波器的设计 六、线性相位FIR带通滤波器的设计 七、线性相位FIR带阻滤波器的设计 八、窗函数法的MATLAB设计函数 7.3 利用频率采样法设计FIR滤波器 7.4 利用等波纹最佳逼近法设计FIR滤波器 7.4 利用等波纹最佳逼近法设计FIR滤波器 （1）理想线性相位带通滤波器的频率响应为 （2）求单位脉冲响应 （3）根据对通带、阻带衰减的要求以及对过渡带宽的要求，选择窗函数w(n)及窗的点数N。由此可求得所需线性相位带通滤波器的单位脉冲响应h(n) （4）求FT[h(n) ]，验证各项指标是否满足要求，如不满足要求，则要改变N，或改变窗形状（或者两者都改变），然后重新计算。 低通 - 低通 （1）理想线性相位带通滤波器的频率响应为 （2）求单位脉冲响应 （3）根据对通带、阻带衰减的要求以及对过渡带宽的要求，选择窗函数w(n)及窗的点数N。由此可求得所需线性相位带通滤波器的单位脉冲响应h(n) （4）求FT[h(n) ]，验证各项指标是否满足要求，如不满足要求，则要改变N，或改变窗形状（或者两者都改变），然后重新计算。 低通+高通 b=fir1(M, wc, ’ftype’, window) M：FIR滤波器阶数，M=N-1 window：窗口类型 window=boxcar(N)，矩形窗 window=bartlett(N)，三角形窗 window=hanning(N)，汉宁窗 wc:对pi归一化的数字频率 wc=(wp+ws)/2，ftype不填，设计FIR低通滤波器 wc=(wp+ws)/2，ftype=high，设计FIR高通滤波器 wc=[wcl,wcu]，ftype不填，设计FIR带通滤波器 wc=[wcl,wcu]，ftype=stop，设计FIR带阻滤波器 window=hamming(N)，哈明窗 window=blackman(N)，布莱克曼窗 window=kaiser(N,alph)，凯塞窗 【例7.2.2】对模拟信号进行低通滤波处理，要求通带0≤f ≤1.5kHz内衰减小于1dB，阻带2.5kHz≤f ≤∞上衰减大于40dB。采用线性相位FIR数字滤波器实现滤波，采样频率Fs=10kHz。用窗函数法设计FIR滤波器，求出h(n)，并画出损耗函数曲线。为降低计算量，希望滤波器阶数尽量低。 Fs=10*1000; fp=1500; fs=2500; rs=40; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; Bt=ws-wp; alph=0.5842*(rs-21)^0.4+0.07886*(rs-21); %计算kaiser窗的参数α，P210 N=ceil((rs-8)/2.285/Bt); %根据(7.2.17)式计算kaiser窗所需阶数N wc=(wp+ws)/2/pi; %计算理想高通滤波器通带截止频率(关于π归一化） hn=fir1(N,wc,kaiser(N+1,alph)); %调用kaiser计算低通FIRDF的h(n) Fs=10*1000; fp=1500;fs=2500;rs=40; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; Bt=ws-wp; alph=0.5842*(rs-21)^0.4+0.07886*(rs-21); %根据(7.2.16)式计算kaiser窗的参数α N=ceil((rs-8)/2.285/Bt); %根据(7.2.17)式计算kaiser窗所需阶数N wc=(wp+ws)/2/pi; %计算理想高通滤波器通带截止频率(关于π归一化） hn=fir1(N,wc,kaiser(N+1,alph)); %调用kaiser计算低通FIRDF的h(n) %以下是绘图部分 M=1024; hk=fft(hn,M); n=0:M; subplot(2,2,1);stem(n,hn,.);line([0,30],[0,0]) xlabel(n); ylabel(h(n)); k = 1:M/2; w = (0:M/2-1) / (M/2); subplot(2,2,2); plot(w,20*log10(abs(hk(k)))); axis([0,1,-80,5]); xla