基于MATLAB的FIR低通滤波器的设计讲述.doc

吉首大学信息科学与工程学院 课程设计报告书 课 程 通信系统设计与仿真 课 题：基于MATLAB的低通FIR滤波器设计 姓 名： 何其兵 杨冰清 岳裕旭 学 号： 20124055066 20124055067 20124055068 专 业： 通信工程（物联网方向） 年 级： 2012级 指导教师： 雷可君 2015年 11 月 项目介绍与设计目的 1.项目介绍 FIR数字滤波器的设计方法有很多，比较常用的有窗函数设计法、频率采样设计法、等波纹逼近法等。本课题通过运用窗函数设计FIR数字低通滤波器，并实现对给定的信号进行滤波。窗函数设计法是最基本的数字滤波方法，是利用傅里叶反变换(IDTFT)计算给定的频响的理想单位脉冲响应，再加以窗函数进行截断和平滑。Matlab软件的信号处理工具箱提供了FIR数字滤波器设计的子函数，运用Matlab软件设计可以避免繁杂的数学运算，而且具有丰富的绘图功能，可以方便地查看所设计的数字滤波器的幅度响应和相位响应是否满足设计要求。因此，本课题在理论分析各种FIR数字滤波器设计方法的基础上，运用Matlab软件进行仿真分析。 随着信息科学和计算机技术的不断发展，数字信号处理(DSP，Digital Signal Processing)的理论和技术也得到了飞速的发展，并逐渐成为一门重要的学科，它的重要性在日常通信、图像处理、遥感、声纳、生物医学、地震、消费电子、国防军事、医疗方面等显得尤为突出。在我们面临的信息革命中，数字信号处理几乎涉及了所有的工程技术领域。 二窗函数法的Matlab实现 Matlab信号处理工具箱提供了基于加窗的线性相位FIR滤波器设计函数fir1。 fir1函数介绍 fir1调用格式为：b=fir1(n, Wc,ftype,window) 函数参数说明如下： 1.n表示滤波器的阶数 2. Wc为滤波器的归一化截止频率，它是一个大于0小于1的一个数 3.ftype表示所设计滤波器的类型，如果ftype=high，则表示高通滤波器，如果ftype=stop，则表示带阻滤波器，如果此时没有参数，就表示低通滤波器 4.window表示的是指定的窗函数，如矩形窗为(n)，三角窗为bartlett(n)，如果缺省window参数，则fir1默认为是哈明窗hamming(n) 基于fir1函数的窗函数法FIR滤波器设计 下面给出利用fir1函数的窗函数数字低通滤波器利用fir1函数进行设计，这种设计方法只需要给出滤波器的阶数、截止频率、窗函数等参数，Matlab即可自行完成设计，并可通过freqz函数查看滤波器的幅频响应和相频响应，以验证滤波器是否满足设计要求。n=27; %设滤器的阶数为27 Wn=0.2; %归一化截止频率为0.2 window=(28); %矩形窗 b=fir1(n,Wn,window); freqz(b); 运行结果：从幅频响应上看，通带，阻带波纹较大，因而阻带特性较不理想，相频响应曲线在通带内为直线，效果较好，信号失真小。 图矩形窗滤波器幅频和相频曲线 利用三角窗进行设计时，原理与窗基本相同，只不过生成窗函数时采用triang（n）函数生成三角窗，程序运行结果： 图三角窗设计的滤波器幅频和相频曲线 从幅频响应上看，通带，相频响应曲线在为直线，。利用窗进行设计时，原理与窗相同，只不过生成窗函数时采用hamming（n）函数生成三角窗，程序运行结果如下： 图-3哈明窗设计的滤波器幅频和相频曲线 利用布莱克曼窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用blackman（n）函数生成窗，程序运行结果如下： 图-4布莱克曼窗设计的滤波器幅频和相频曲线 利用窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用kaiser（n,beta）函数生成窗，程序运行结果： 图-5凯塞窗设计的滤波器幅频和相频曲线 窗函数设计法对信号加窗之后会使不连续点处边沿加宽形成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主宽度。在处出现肩峰值，两侧形成振荡，振荡的幅度和多少取决于旁瓣的幅度和多少，改变N只能改变窗谱的主宽度，不能改变主瓣旁瓣的相对比例，其相对比例由窗函数形状决定。 窗函数设计法设计FIR数字滤波器是傅里叶变换的典型运用，而频率采样法设计的指导思想是通过频域采样点实现的，同时为了保证过渡带宽的不变，滤波器的采样点数也要相应