实验四-窗函数法设计FIR数字滤波器(共16页).docx
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《实验四-窗函数法设计FIR数字滤波器(共16页).docx》由会员分享,可在线阅读,更多相关《实验四-窗函数法设计FIR数字滤波器(共16页).docx(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布莱克曼窗
3、Hanning汉宁窗函数Kaiser凯瑟窗函数Hamming海明窗表4-2 常用窗函数的特性窗函数窗函数频率特性加窗后滤波器指标旁瓣峰值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.2*pi;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)*(n-alpha);w_boxcar=boxcar(N);h=hd.*w_boxcar;subplot(221);stem(n,hd,filled);axis tight;xlabel(n
6、);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,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);grid on;MATLAB图形:(2)汉宁窗程序代码:wp=0.2*pi
7、;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(H(om
8、ega);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);grid on;MATLAB图形:(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/p
9、i)*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,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)
10、;plot(w/pi,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);grid on;MATLAB图形:(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);ax
11、is 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,20*log10(abs(H)/max(H);axis tight;xlabel(omega/pi);ylabel(dB);grid on;MATLAB图形为:(5
12、)凯瑟窗程序代码: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);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 函数 设计 FIR 数字滤波器 16
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内