用窗函数法设计FIR数字低通滤波器

《用窗函数法设计FIR数字低通滤波器》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、河北科技大学课程设计报告学生姓名：学号：专业班级：课程名称：学年学期指导教师：20年月-课程设计成绩评定表学生姓名学号成绩专业班级起止时间设计题目指导教师评语指导教师：年月日目录1.窗函数设计低通滤波器1.1设计目的……………………………………………………………………11.2设计原理推导与计算………………………………………………………11.3设计内容与要求……………………………………………………………21.4设计源程序与运行结果……………………………………………………31.5思考题……………………………………………………………………10

2、1.6心得体会…………………………………………………………………14参考文献………………………………………………………………………151.窗函数设计低通滤波器1.1设计目的1.熟悉设计线性相位数字滤波器的一般步骤。2.掌握用窗函数法设计FIR数字滤波器的原理和方法。3.熟悉各种窗函数的作用以及各种窗函数对滤波器特性的影响。4.学会根据指标要求选择合适的窗函数。1.2设计原理推导与计算如果所希望的滤波器的理想的频率响应函数为，则其对应的单位脉冲响应为（4.1）窗函数设计法的基本原理是设计设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函

3、数，即,其中用有限长单位脉冲响应序列逼近。由于往往是无限长序列，而且是非因果的，所以用窗函数将截断，并进行加权处理，得到：（4.2）就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为（4.3）式中，N为所选窗函数的长度。用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N-15-的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见表（一）。窗函数旁瓣峰值幅度/dB过渡带宽阻带最小衰减/dB过渡带带宽矩形窗-134π/N-120.9三

4、角形窗-258π/N-252.1汉宁窗-318π/N-443.1哈明窗-418π/N-533.3不莱克曼窗-5712π/N-745.5凯塞窗(=7.865)-5710π/N-805表（一）各种窗函数的基本参数这样选定窗函数类型和长度N之后，求出单位脉冲响应，并按照式（4.3）求出。是否满足要求，如果不满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。1.3设计内容与要求（一）设计要求：1.学会计算滤波器各项性能指标及如何来满足给定的指标要求。2.用MATLAB语言编程实现给定指标要求的滤波器的设计。3.熟悉MATLAB语言，独

5、立编写程序。4.设计低通FIR滤波器的指标：-15-通带最大波动，阻带最小衰减，（二）、设计内容：1.熟悉各种窗函数,在MATLAB命令窗下浏览各种窗函数，绘出（或打印）所看到的窗函数图。2.编写计算理想低通滤波器单位抽样响应hd(n)的m函数文件ideal.m。3.编写计算N阶差分方程所描述系统频响函数的m函数文件fr.m。4.根据指标要求选择窗函数的形状与长度N。（至少选择两种符合要求的窗函数及其对应的长度）。5.编写.m程序文件，通过调用ideal.m和fr．m文件，计算你设计的实际低通FIR滤波器的单位抽样响应h(n)和频率响应，打印

6、在频率区间[O，π]上的幅频响应特性曲线，幅度用分贝表示。6.验证所设计的滤波器是否满足指标要求。1.4设计的源程序及运行结果：1、利用MATLAB窗口观察各种窗函数：-15-%巴特利特窗w=bartlett(20);subplot(3,2,1);plot(w);stem(w,'y');%'y'表示黄色%stem表示以离散图输出title('巴特利特床窗');xlabel('n');%横坐标为nylabel('w(n)');%纵坐标为w(n)%布莱克曼窗w=blackman(20);subplot(3,2,2);plot(w);stem(w,

7、'b');%'b'表示蓝色title('布莱克曼窗');xlabel('n');ylabel('w(n)');%矩形窗w=boxcar(20);subplot(3,2,3);plot(w);stem(w,'r');title('矩形窗');xlabel('n');-15-ylabel('w(n)');%海明窗w=hamming(20);plot(w);stem(w,'m');%'m'表示紫色title('海明窗');xlabel('n');ylabel('w(n)');%汉宁窗w=hanning(20);subplot(3,2,5);plot

8、(w);stem(w,'g');%'g'表示绿色title('汉宁窗');xlabel('n');ylabel('w(n)');%凯泽窗beta=5.6533;w=