《实验四 窗函数法设计fir数字滤波器》由会员分享，可在线阅读，更多相关《实验四 窗函数法设计fir数字滤波器（16页珍藏版）》请在金锄头文库上搜索。
1、数字信号处理实验报告1实验四 窗函数法设计 FIR 数字滤波器一、实验目的1、掌握窗函数法设计 FIR 数字滤波器的原理及具体方法。2、掌握频率取样法设计 FIR 数字滤波器的原理和基本方法。3、学习利用窗函数法和频率取样法设计低通、带通、高通、带阻数字滤波器。二、实验环境计算机、MATLAB 软件三、实验基础理论窗函数设计 FIR 滤波器1.基本原理窗函数设计法的基本思想为，首先选择一个适当的理想的滤波器 ，然后()jdHe用窗函数截取它的单位脉冲响应 ，得到线性相位和因果的 FIR 滤波器。这种方法的(n)dh重点是选择一个合适的窗函数和理想滤波器，使设计的滤波器的单位脉冲响应逼近理想滤波
2、器的单位脉冲响应。2.设计步骤（1）给定理想滤波器的频率响应 ，在通带上具有单位增益和线性相位，()jdHe在阻带上具有零响应。一个带宽为 的低通滤波器由下式给定：c |,0)(,|,)( cjdjajd eeH其中 为采样延迟，其作用是为了得到一个因果系统。（2）确定这个滤波器的单位脉冲响应 )(sin)(ahcd为了得到一个 长度为 N 的因果的线性相位 FIR 滤波器，我们令21Na（3）用窗函数截取 得到所设计 FIR 数字滤波器： (n)dh )()(nRhnNd3.窗函数的选择常用的窗函数有矩形（Rectangular）窗，汉宁（Hanning）窗，海明（Hamming）窗、布莱克
3、曼（Blackman）窗、凯瑟（Kaiser）窗等表 4-1 MATLAB 中产生窗函数的命令MATLAB 函数 窗函数 MATLAB 函数 窗函数Boxcar 矩形窗函数 Blackman 布莱克曼窗Hanning 汉宁窗函数 Kaiser 凯瑟窗函数Hamming 海明窗数字信号处理实验报告2表 4-2 常用窗函数的特性窗函数频率特性 加窗后滤波器指标窗函数 旁瓣峰值dB 主瓣宽度 过渡带宽 最小阻带衰减 dB矩形窗 -13 4/N 1.8/N -21汉宁窗 -31 8/N 6.2/N -44海明窗 -41 8/N 6.6/N -53布莱克曼窗 -57 12/N 11/N -74凯瑟窗是一
4、种广泛在实际中广泛应用的窗函数，它由下式给定： 2001()()nINn其中 是修正的零阶贝塞尔函数，参数 控制最小阻带衰减，这种窗函数对于相0Ix同的 N 可以提供不同的过渡带宽。由于贝塞尔函数比较复杂，这种窗函数的设计方程很难推导，然而幸运的是，有一些经验设计方程可以直接使用。已知给定的指标 ，滤波器长度 N 和凯瑟窗参数 可以按如下凯瑟窗,RpstpsA和 方程给出过渡带带宽： stp7.95128sAN0.40.(),5786(21),50ssssAA频率取样设计 FIR 滤波器1.基本原理频率取样法从频域出发，把理想的滤波器 等间隔采样得到 ,将()jdHe()dHk作为实际设计滤波
5、器的 ：()dHk()k2()|0,1jdkNkHe L得到 以后可以由 来确定唯一确定滤波器的单位脉冲响应 , 可()() ()hn)jHe以由 求得：Hk数字信号处理实验报告310()()2)NjkhnIDFTHek其中 为内插函数：()x12sin/)(NjeN（有 求得的频率响应 将逼近 。()Hk()jH()jde如果我们设计的是线性相位 FIR 滤波器，则 的幅度和相位满足线性相位滤波器k的约束条件。我们将 表示为如下形式()k()()=jkjkrHeHe当 为实数，则()hn*()kNk由此得到 ()()rrH即 为中心偶对称。在利用线性相位条件可知，对于 1 型和 2 型线性/
6、2rk以相位滤波器： 11()0,2()2,NkNkL对于 3 型和 4 型线性相位滤波器 11()0,22() ,NkNk Lm2.设计步骤（1）由给定的理想滤波器给出 和 。()rHk（2）由 求得()()()=jkjrHkee()k数字信号处理实验报告4（3）根据 求得 或()Hk()hn)je四、实验内容1、设计一个数字低通 FIR 滤波器，其技术指标如下：0.2,.5ppRdB3stsA分别采用矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗设计该滤波器。结合实验结果，分别讨论采用上述方法设计的数字滤波器是否都能满足给定指标要求。（1）矩形窗程序代码：wp=0.2*pi;wst=0.3*p
7、i;tr_width=wst-wp;N=ceil(1.8*pi/tr_width)n=0:N-1;wc=(wst+wp)/2;alpha=(N-1)/2;hd=(wc/pi)*sinc(wc/pi)*(n-alpha);w_boxcar=boxcar(N);h=hd.*w_boxcar;subplot(221);stem(n,hd,filled);axis tight;xlabel(n);ylabel(hd(n);Hr,w1=zerophase(h);subplot(222);plot(w1/pi,Hr);axis;xlabel(omega/pi);ylabel(H(omega);subplo
8、t(223);stem(n,h,filled);axis tight;xlabel(n);ylabel(h(n);H,w=freqz(h,1);subplot(224);plot(w/pi,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);grid on;MATLAB 图形：数字信号处理实验报告50 5 10 1500.10.2nhd(n)0 0.5 1-0.500.511.5/H()0 5 10 1500.10.2nh(n)0 0.2 0.4 0.6 0.8-80-60-40-200/dB（2）汉宁窗程序代码：wp=0
9、.2*pi;wst=0.3*pi;tr_width=wst-wp;N=ceil(6.2*pi/tr_width)n=0:N-1;wc=(wst+wp)/2;alpha=(N-1)/2;hd=(wc/pi)*sinc(wc/pi)*(n-alpha);w_boxcar=hanning(N);h=hd.*w_boxcar;subplot(221);stem(n,hd,filled);axis tight;xlabel(n);ylabel(hd(n);Hr,w1=zerophase(h);subplot(222);plot(w1/pi,Hr);axis;xlabel(omega/pi);ylabel
10、(H(omega);subplot(223);stem(n,h,filled);axis tight;xlabel(n);ylabel(h(n);H,w=freqz(h,1);subplot(224);plot(w/pi,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);数字信号处理实验报告6grid on;MATLAB图形：0 20 40 6000.10.2nhd(n)0 0.5 1-0.500.511.5/H()0 20 40 6000.050.10.150.2nh(n)0 0.2 0.4 0.6 0.8-120-1
11、00-80-60-40-200/dB（3）海明窗程序代码：wp=0.2*pi;wst=0.3*pi;tr_width=wst-wp;N=ceil(6.6*pi/tr_width)n=0:N-1;wc=(wst+wp)/2;alpha=(N-1)/2;hd=(wc/pi)*sinc(wc/pi)*(n-alpha);w_boxcar=hamming(N);h=hd.*w_boxcar;subplot(221);stem(n,hd,filled);axis tight;xlabel(n);ylabel(hd(n);Hr,w1=zerophase(h);subplot(222);plot(w1/pi
12、,Hr);axis;xlabel(omega/pi);ylabel(H(omega);subplot(223);stem(n,h,filled);axis tight;xlabel(n);ylabel(h(n);H,w=freqz(h,1);数字信号处理实验报告7subplot(224);plot(w/pi,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);grid on;MATLAB图形：0 20 40 6000.10.2nhd(n)0 0.5 1-0.500.511.5/H()0 20 40 6000.10.2nh(
13、n)0 0.2 0.4 0.6 0.8-100-500/dB（4）布莱克曼窗程序代码：wp=0.2*pi;wst=0.3*pi;tr_width=wst-wp;N=ceil(11*pi/tr_width)n=0:N-1;wc=(wst+wp)/2;alpha=(N-1)/2;hd=(wc/pi)*sinc(wc/pi)*(n-alpha);w_boxcar=blackman(N);h=hd.*w_boxcar;subplot(221);stem(n,hd,filled);axis tight;xlabel(n);ylabel(hd(n);Hr,w1=zerophase(h);subplot(2
14、22);plot(w1/pi,Hr);axis;xlabel(omega/pi);ylabel(H(omega);subplot(223);stem(n,h,filled);axis tight;xlabel(n);ylabel(h(n);H,w=freqz(h,1);subplot(224);数字信号处理实验报告8plot(w/pi,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);grid on;MATLAB图形为：0 50 10000.10.2nhd(n)0 0.5 1-0.500.511.5/H()0 50 10
15、000.10.2nh(n)0 0.2 0.4 0.6 0.8-150-100-500/dB（5）凯瑟窗程序代码：wp=0.2*pi;wst=0.3*pi;tr_width=wst-wp;As=50;N=ceil(As-7.95)/(2.285*tr_width)+1;beta=0.1102*(As-8.7);n=0:N-1;wc=(wst+wp)/2;alpha=(N-1)/2;hd=(wc/pi)*sinc(wc/pi)*(n-alpha);w_boxcar=kaiser(N,beta);h=hd.*w_boxcar;subplot(221);stem(n,hd,filled);axis tight;xlabel(n);ylabel(hd(n);Hr,w1=zerophase(h);subplot(222);plot(w1/pi,Hr);axis;xlabel(omega/pi);ylabel(H(omega);subplot(223);stem(n,h,filled);ax
