《实验四 窗函数法设计FIR数字滤波器》由会员分享，可在线阅读，更多相关《实验四 窗函数法设计FIR数字滤波器（16页珍藏版）》请在金锄头文库上搜索。
1、实验四 窗函数法设计FIR数字滤波器一、实验目的1、掌握窗函数法设计FIR数字滤波器的原理及具体方法。 2、掌握频率取样法设计FIR数字滤波器的原理和基本方法. 3、学习利用窗函数法和频率取样法设计低通、带通、高通、带阻数字滤波器。二、实验环境计算机、MATLAB软件三、实验基础理论窗函数设计FIR滤波器1。基本原理窗函数设计法的基本思想为，首先选择一个适当的理想的滤波器，然后用窗函数截取它的单位脉冲响应,得到线性相位和因果的FIR滤波器。这种方法的重点是选择一个合适的窗函数和理想滤波器，使设计的滤波器的单位脉冲响应逼近理想滤波器的单位脉冲响应。2.设计步骤（1）给定理想滤波器的频率响应，在通
2、带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽为的低通滤波器由下式给定： 其中为采样延迟，其作用是为了得到一个因果系统。（2）确定这个滤波器的单位脉冲响应为了得到一个长度为N的因果的线性相位FIR滤波器，我们令 （3）用窗函数截取得到所设计FIR数字滤波器： 3.窗函数的选择常用的窗函数有矩形(Rectangular）窗，汉宁（Hanning）窗，海明(Hamming）窗、布莱克曼（Blackman）窗、凯瑟(Kaiser）窗等表4-1 MATLAB中产生窗函数的命令MATLAB函数窗函数MATLAB函数窗函数Boxcar矩形窗函数Blackman布莱克曼窗Hanning汉宁窗函数K
3、aiser凯瑟窗函数Hamming海明窗表42 常用窗函数的特性窗函数窗函数频率特性加窗后滤波器指标旁瓣峰值dB主瓣宽度过渡带宽最小阻带衰减dB矩形窗134/N1。8/N-21汉宁窗318/N6.2/N-44海明窗-418/N6。6/N-53布莱克曼窗5712/N11/N-74凯瑟窗是一种广泛在实际中广泛应用的窗函数，它由下式给定：其中是修正的零阶贝塞尔函数,参数控制最小阻带衰减，这种窗函数对于相同的N可以提供不同的过渡带宽.由于贝塞尔函数比较复杂,这种窗函数的设计方程很难推导，然而幸运的是，有一些经验设计方程可以直接使用.已知给定的指标 ,滤波器长度N和凯瑟窗参数可以按如下凯瑟窗方程给出过渡
4、带带宽：频率取样设计FIR滤波器1.基本原理频率取样法从频域出发，把理想的滤波器等间隔采样得到,将作为实际设计滤波器的：得到以后可以由来确定唯一确定滤波器的单位脉冲响应,可以由求得:其中为内插函数：有求得的频率响应将逼近。如果我们设计的是线性相位FIR滤波器，则的幅度和相位满足线性相位滤波器的约束条件。我们将表示为如下形式当为实数,则由此得到即为中心偶对称.在利用线性相位条件可知,对于1型和2型线性相位滤波器：对于3型和4型线性相位滤波器2.设计步骤(1）由给定的理想滤波器给出和.（2)由求得(3）根据求得或四、实验内容1、设计一个数字低通FIR滤波器，其技术指标如下：分别采用矩形窗、汉宁窗、
5、海明窗、布莱克曼窗、凯瑟窗设计该滤波器。结合实验结果，分别讨论采用上述方法设计的数字滤波器是否都能满足给定指标要求。（1） 矩形窗程序代码：wp=0。2pi;wst=0。3*pi；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)(nalpha);w_boxcar=boxcar（N)；h=hd。*w_boxcar；subplot(221）；stem（n，hd，filled);axis tight；xlabel(n）;ylabel（hd（n）；Hr,w1
6、=zerophase（h）；subplot(222）;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);plot(w/pi,20log10(abs(H）/max(H）)）;axis tight;xlabel（omega/pi)；ylabel(dB)；grid on;MATLAB图形：(2）汉宁窗程序代码:wp=0.2*pi；wst=0.3*pi;tr_wid
7、th=wstwp;N=ceil（6.2pi/tr_width）n=0:N1；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（H(omega)；subplot（223)；stem（
8、n，h，filled);axis tight;xlabel(n）；ylabel（h(n）;H,w=freqz(h,1）;subplot(224);plot（w/pi,20log10(abs(H)/max(H）)）；axis tight；xlabel(omega/pi）;ylabel（dB)；grid on；MATLAB图形：(3)海明窗程序代码：wp=0。2pi；wst=0。3*pi;tr_width=wstwp；N=ceil（6。6*pi/tr_width)n=0：N1;wc=(wst+wp)/2;alpha=（N1）/2;hd=(wc/pi)sinc(wc/pi)(n-alpha）；w_b
9、oxcar=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,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);plot（w/pi，20log10(abs(H）/ma
10、x(H）;axis tight；xlabel(omega/pi);ylabel(dB)；grid on;MATLAB图形:（4）布莱克曼窗程序代码：wp=0.2*pi；wst=0。3pi；tr_width=wst-wp;N=ceil(11pi/tr_width）n=0:N1;wc=(wst+wp)/2;alpha=（N1)/2;hd=(wc/pi）*sinc(wc/pi）（nalpha);w_boxcar=blackman（N）;h=hd。w_boxcar;subplot（221);stem（n，hd，filled）;axis tight；xlabel(n）；ylabel（hd（n）;Hr,w
11、1=zerophase（h）；subplot（222)；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）；plot（w/pi，20log10(abs(H）/max(H）);axis tight；xlabel（omega/pi)；ylabel（dB)；grid on；MATLAB图形为:（5）凯瑟窗程序代码：wp=0。2*pi；wst=0.3pi；tr_wid
12、th=wstwp;As=50；N=ceil（As-7.95）/（2。285tr_width）+1;beta=0。1102(As8。7）;n=0：N1;wc=(wst+wp)/2;alpha=(N-1）/2;hd=（wc/pi)sinc（(wc/pi)*(nalpha)）;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(
13、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);plot（w/pi，20log10(abs(H）/max（H）)；axis tight;xlabel(omega/pi)；ylabel（dB）;grid on；MATLAB图形：2、设计一个数字带通FIR滤波器，其技术指标如下：下阻带边缘：下通带边缘:上通带边缘:上阻带边缘：程序代码:wp1=0。2*pi；Rp1=1；wst1=0。35pi；A1=60；widt
14、h1=wst1wp1；N1=ceil(11pi/width1)+1;n1=0:（N1-1）；wc1=(wp1+wst1)/2；alpha=（N11)/2；wp2=0。65*pi;Rp2=1；wst2=0.8*pi;A2=60；width2=wst2wp2;N2=ceil(11pi/width2)+1；n2=0:（N2-1);wc2=(wp2+wst2）/2;alpha=（N21）/2；hd=(wc2/pi)*sinc(wc2/pi）*（n2-alpha)）（wc1/pi）*sinc(（wc1/pi）（n1-alpha）；w_w=blackman(N1)； h=hd.w_w; subplot（22
