课程设计：窗函数法设计fir滤波器

《课程设计：窗函数法设计fir滤波器》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、目录目录目录0第一章课程设计的目的1第二章主要设计内容2第三章设计原理及思想33.1设计FIR滤波器的原理33.2各种窗函数的具体参数4第四章主程序框图5第五章参考程序6第六章FIR滤波器图像对比96.1低通滤波器设计结果图96.2带通滤波器设计结果图96.3高通滤波器设计结果图10第七章滤波器的应用117.1应用程序117.2滤波结果图12第八章心得体会14第九章参考文献1518第一章课程设计的目的语音信号处理的好坏将影响并导致语音信号的好坏，只有通过将这些语音信号进行一系列的的数字处理，才能将那些非必要

2、的噪声杂波妥善滤除，得到纯净的单纯的语音信号，因此语音信号去噪处理是具有现实意义的研究课程。通过课设可以达到以下目的：1、通过课程设计把自己在大学中所学的知识应用到实践当中。2、深入了解利用Matlab设计FIR数字滤波器的基本方法。3、在课程设计的过程中掌握程序编译及软件设计的基本方法。4、提高自己对于新知识的学习能力及进行实际操作的能力。5、锻炼自己通过网络及各种资料解决实际问题的能力。18第二章主要设计内容利用窗函数法设计FIR滤波器，绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进

3、行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。1、学习用窗函数法设计FIR数字滤波器的原理及其设计步骤；2、用Matlab方法对FIR数字滤波器数学模型的系数进行求解；并用窗函数法完成低通、高通、带通滤波器的设计；采用频率为16KHZ、截止频率为200HZ，阶数为30阶；3、用三角形窗、哈明窗、布莱克曼窗完成FIR数字滤波器的设计，并对结果进行比较和分析。18第三章设计原理及思想3.1设计FIR滤波器的原理如果所希望的滤波器的理想的频率响应函数为，则其对应的单位脉冲响应为（3.1）窗函数设计法的

4、基本原理是用有限长单位脉冲响应序列逼近。由于往往是无限长序列，而且是非因果的，所以用窗函数将18截断，并进行加权处理，得到：（3.2）就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为（3.3）式中，N为所选窗函数的长度。我们知道，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见表1。这样选定窗函数类型和长度N之后，求出单位脉冲响应，并按照式

5、（3.3）求出。是否满足要求，要进行演算。一般在尾部加零使长度满足2的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果不满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。如果要求线性相位特性，则还必须满足根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类，例如，要设计线性相位低通特性，可以选择这一类，而不能选择这一类。3.2各种窗函数的具体参数18表1各种窗函数的基本参数窗函数旁瓣峰值幅度/dB过渡带宽阻带最小衰减/

6、dB矩形窗-134π/N-12三角形窗-258π/N-25汉宁窗-318π/N-44哈明窗-418π/N-53不莱克曼窗-5712π/N-74凯塞窗(α=7.865)-5710π/N-8018第四章主程序框图开始读入窗口长度N计算hd(n)调用窗函数子程序求w(n)调用子程序（函数）计算H(k)=DFT[h(n)]调用绘图子程序（函数）绘制H(k)幅度相位曲线结束图1主程序框图计算h(n)=hd(n)w(n)18第五章参考程序主要程序：近似计算截止频率wc=0.5*(wp+ws)/pi；过渡带宽Bt=ws-

7、wp或wp-ws；三角形窗滤波器阶数N1=ceil(6.1*pi/Bt);哈明窗滤波器阶数N1=ceil(6.6*pi/Bt);布莱克曼窗滤波器阶数N1=ceil(11*pi/Bt);实现此题的程序如下：%低通滤波closeall;clc;clearall;%清零，清屏，关闭其他窗口wpl=pi/4;%通带截止频率wsl=pi/2;%阻带截止频率Btl=wsl-wpl;%过渡带宽wcl=0.5*(wpl+wsl)/pi;%近似计算截止频率N1l=ceil(6.1*pi/Btl);%根据过渡带宽度求滤波器阶数

8、window=bartlett(N1l+1);%三角窗b1l=fir1(N1l,wcl,window);%三角窗低通滤波器设计N2l=ceil(6.6*pi/Btl);%根据过渡带宽度求滤波器阶数window=hamming(N2l+1);%哈明窗b2l=fir1(N2l,wcl,window);%哈明窗低通滤波器设计18N3l=ceil(11*pi/Btl);%根据过渡带宽度求滤波器阶数window=black