基于MATLAB的IIR数字滤波器的设计与DSP实现开题报告

《基于MATLAB的IIR数字滤波器的设计与DSP实现》开题

报告

一、课题的目的和意义

1、研究目的

我选择了《基于MATLAB的IIR滤波器的设计》这个课题作为毕业设计其主要目的是通过此次课程设计进一步学习和巩固数字信号处理和DSP及其相关知识，并学会利用所学的知识技能，在设计过程中能综合运用所学知识内容，进一步熟悉和掌握MATLAB和DSP的使用方法；对IIR滤波器的原理及其实现有较深的了解；为即将进入社会参加工作打下坚实的基础；掌握收集资料、消化资料和综合资料的能力等等。

2、研究意义

随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。

随着电子工业的发展，对滤波器的性能要求越来越高，功能也越来越多，并且要求它们向集成方向发展。我国滤波器研制和生产与上述要求相差甚远，为缩短这个差距，电子工程和科技人员负有重大的历史责任。

从事电子通信业而不能熟练操作使用MATLAB电子线路设计软件，在工作和学习中将是寸步难行的。在数学、电子、金融等行业，使用MATLAB等计算机软件对产品进行设计、仿真在很早以前就已经成为了一种趋势，这类软件的问世也极大地提高了设计人员在通信、电子等行业的产品设计质量与效率。

二、文献综述

数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高

通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运行在0与1两个电平状态)、灵活性强等优点。数字滤波器按单位脉冲响应的性质可分为无限长单位脉冲响应滤波器IIR和有限长单位脉冲响应滤波器(FIR)两种。IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法。FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。

目前数字滤波器的设计有许多现成的高级语言设计程序，但他们都存在设计效率较低，不具有可视图形，不便于修改参数等缺点，而MATLAB为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。他以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。尤其是MATLAB工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。其中的信号处理工具箱、图像处理工具箱、小波工具箱等更是为数字滤波研究的蓬勃发展提供了可能。

三、研究（设计）内容和拟定解决的关键问题

1．研究内容

（1）研究数字滤波器及其传统设计方法。

（2）研究基于MATLAB的IIR滤波器设计算法。

（3）研究利用FDATool工具来设计IIR滤波器及MATLAB语言设计IIR滤波器。（4）研究DSP实现IIR滤波器的算法及DSP芯片TMS320VC5416的结构和使用方法。2．研究途径

收集资料和文献，系统学习数字滤波器的原理、MATLAB软件和DSP原理及应用。按照IIR滤波器的设计理论建立数学模型，建立完整仿真模型并仿真。用MS320VC5416DSP芯片作为实现工具来实现MATLAB仿真出来的IIR滤波器效果，然后不断调试程序直至得到最佳实现效果。

3．技术路线

（1）根据IIR滤波器的原理，利用MA TLAB进行IIR数字低通滤波器设计,首先

是在给定滤波器性能情况下,确定通带和阻带的截止频率,还有通带和阻带区的波纹系数。然后通过But terworth滤波器buttord函数计算出满足指定性能的模拟But terworth滤波器的阶和截止频率。有了上述的前提和基础,我们可利用MA TLAB语言中Butter函数设计出模拟域的Butterworth滤波器。然后再利用impinvar函数,用脉冲响应不变法实现模拟到数字滤波器的转换。

（2）借助于计算机强大的仿真功能，在MATLAB环境下设计一个IIR数字滤波器变得十分容易。但是要在DSP上实现此滤波器功能，则需要对滤波器的结构、系数等参数进行改进。从理论上说，可以用高阶IIR数字滤波器实现良好的滤波效果。但由于DSP本身有限字长和精度的因素，加上IIR滤波器在结构上存在反馈回路，是递归型的，再者高阶滤波器参数的动态范围很大。这样一来造成两个后果：结果溢出和误差增大，从而导致算法无法在DSP上实现。解决此问题的有效方法是把高阶IIR数字滤波器简化成几个2阶滤波器来设计，即采用级联结构。

4．拟解决的关键问题

所要解决本课题的主要问题是：IIR滤波器的设计和仿真，以及在DSP芯片上进行调试和实现。

四、研究（设计）方案与进度计划安排

1．研究（设计）方案：

在MATLAB上对IIR滤波器进行设计设计时，有两种设计方案，基于此，对两种设计方案进行如下说明。

方案一：程序设计法。MATLAB信号处理工具箱提供了几个直接设计IIR数字滤波器的函数，直接调用就可以设计滤波器，这为设计通用滤波器提供了方便。

方案二：FDATool界面设计。FDATool（Filter Design&Analysis Tool）是MATLAB信号处理工具箱里专用的滤波器设计分析工具。FDATool可以设计几乎所有的常规滤波器，包括FIR和IIR的各种设计方法。它操作简单，方便灵活。

方案三：模拟原型法。采用经典低通滤波器作为连续域上的设计模型，通过频域变换得到IIR数字滤波器，最后还要进行离散化处理。MATLAB提供的低通模拟滤波器原型函数包括：besselap，buttap，cheb1lp，cheb2ap，ellipap；频域变换函数包括：lp2bp，lp2bs，lp2hp，lp2lp;离散化处理函数有bilinear 和impinvar。

方案四：直接设计法。直接设计方法的思想是基于给定的滤波器参数直接在离散域上寻找合适的数字滤波器，他不限于常规的滤波器类型，如低通、高通、带通和带阻等。这种方法甚至可以设计多带的频率响应，MATLAB提供yulewalk 函数用于辅助设计。

方案五：通用Butterworth设计方法。使用这种方法设计的Butterworth 数字滤波器可以有不同数目的零点和极点，MATLAB提供的maxflat函数实现了这一功能。这个函数与butter函数很相似，但他可以指定两个阶参数，其中归一化和非归一化各一个。如果这两个参数的值相同，那么他与butter函数的结果就是相同的。

方案六：参数建模法。寻找接近于所需要设计的滤波器的通用模型，时域上的建模函数为lpc，prony，Stmcb;频域上的建模函数有invfreqs和invfreqz。比较上述方案，在该课题里，我将采用方案一设计一个低通滤波器，采用方案二设计一个带阻滤波器。

2．进度计划安排

阶段起始日期终止日期进度技术指标

第一阶段2008.03.302008.4.18完成选题和资料收集，

撰写开题报告

全面丰富课题

内容

第二阶段2008.04.192008.05.4IIR设计仿真及训练符合精度要求

第三阶段2008.05.52008.05.25完成设计和软件编程设计合理机构化、模块化

第四阶段2008.05.262008.06.21进行在线调试与仿真，

撰写毕业论文

实现制动控制

第五阶段2008.06.222008.06.29细化论文，准备答辩详细全面，提交论文

五、预期结果和创新成果

设计生成IIR低通和带通数字滤波器的幅频和相频图，然后调用程序和SPTool工具箱对所设计的滤波器进行仿真，并做出待处理噪音信号的滤波前后的相关图形。

TMS320VC5416是16位DSP。在CCS下,它的浮点数格式占用两个字，即32bit。通过CCS的数据输入输出功能和RTDX功能，可以随时观察DSP结果和理论仿真结果之间的差异，及时发现数据的溢出和计算误差。数据采集时要先通过DSP的DMA功能将输入数据保存在外部SRAM中，采样结束后再将要滤波的数据通过DMA载入内部RAM中，从而提高处理速度。考虑到IIR数字滤波器的运算特点，某一时刻的输入数据参加一次运算后不再参与下次运算，只有各级滤波器输出结果需要参与多次运算，因此滤波结果和输入数据可以共用同一存储地址，从而节省存储空间。

除此之处，还要注意以下问题：

(1)IIR数字滤波器是针对采样率固定的系统设计的，更改系统或系统采样率改变时应重新设计。

(2)用到的数据存储单元，程序开始时应对这些单元进行初始化。如果要进行连续滤波，应保存上一次滤波的结果。

(3)IIR数字滤波器的滤波结果会引起相位的延迟，故对相位严格要求的场合，需加全通网络进行较正，或选用FIR滤波器实现滤波功能。

六、成文时间和提交成果形式

2008年6月22日完成论文“基于MATLAB的IIR滤波器的设计”。

提交成果：开题报告、设计论文、有关设计的英文翻译、毕业设计手册。参考文献目录

1．戴明桢周建江编著，TMS320C54x DSP结构、原理及应用，北京航空航天大学出版社，2004

2．郑红吴冠编著，TMS320C54x DSP应用系统设计，北京航空航天大学出版社，2003 3．程佩青编著，数字信号处理教程（第二版），清华大学出版社，2002

4．蒋志凯．数字滤波器与卡尔曼滤波[M]．北京：中国科学技术出版社，1993

5．楼顺天，李博菡．基于MATLAB的系统分析与设计一信号处理[M]．西安：西安电子科技大学出版社，1998．

6．胡广书．数字信号处理、理论、算法与实现[M]．北京：清华大学出版社，1997

7．倪养华,等．数字信号处理与实现．上海：上海交通大学出版社,1998

8．施阳．等．MATLAB语言工具箱．西安：西北工业大学出版社,1999

9．桑国明,刘智.数字滤波器的DSP实现[J].大连大学学报,2001,22(4):69-72. 10．苏涛,庄德靖,吴顺君.一种新的滤波器快速优化算法[J].雷达科学与技术,2004, 2(6):364-369.

11．成跃乐，李东兴，王文畅.MATLAB6.1环境下IIR数字滤波器的设计[J].山东理工大学学报(自然科学版)，2003；17（5）：82~85

matlab滤波器设计

长安大学数字信号处理综合设计专业_______电子信息工程_______ 班级__24030602___________ 姓名_______张舒_______ 学号2403060203 指导教师陈玲日期_______2008-12-27________

一、课程设计目的： 1. 进一步理解数字信号处理的基本概念、基本理论和基本方法； 2.熟悉在Windows环境下语音信号采集的方法； 3.学会用MATLAB软件对信号进行分析和处理； 4.综合运用数字信号处理理论知识，掌握用MATLAB软件设计FIR和IIR数字滤波器的方法； 5. 提高依据所学知识及查阅的课外资料来分析问题解决问题的能力。二、课程设计内容： 1.语音信号的采集利用windows下的录音机录制一段自己的话音，时间控制在1秒左右；并对语音信号进行采样，理解采样频率、采样位数等概念。 2.语音信号的频谱分析利用函数fft对采样后语音信号进行快速傅里叶变换，得到信号的频谱特性。 3.设计数字滤波器采用窗函数法和双线性变换法设计下列要求的三种滤波器，根据语音信号的特点给出有关滤波器的性能指标： 1）低通滤波器性能指标，fp=1000Hz，fc=1200Hz， As=100dB，Ap=1dB； 2）高通滤波器性能指标，fc=4800Hz，fp=5000Hz ，As=100dB，Ap=1dB； 3）带通滤波器性能指标，fp1=1200Hz，fp2=3000Hz，fc1=1000Hz，fc2=3200Hz，As=100dB，Ap=1dB。 4.对语音信号进行滤波比较用两种方法设计的各滤波器的性能，然后用性能好的滤波器分别对采集的语音信号进行滤波；并比较滤波前后语音信号的波形及频谱，分析信号的变化。 5.回放语音信号，感觉滤波前后的声音变化。三、实验原理（一）基于双线性Z变换法的IIR数字滤波器设计由于的频率映射关系是根据推导的，所以使jΩ轴每隔2π/Ts便映射到单位圆上一周，利用冲激响应不变法设计数字滤波器时可能会导致上述的频域混叠现象。为了克服这一问题，需要找到由s平面到z平面的另外的映射关系，这种关系应保证： 1) s平面的整个jΩ轴仅映射为z平面单位圆上的一周； 2) 若G(s)是稳定的，由G(s)映射得到的H(z)也应该是稳定的； 3) 这种映射是可逆的，既能由G(s)得到H(z)，也能由H(z)得到G(s)； 4) 如果G(j0)=1，那么。双线性Z变换满足以上4个条件的映射关系，其变换公式为

matlab中关于数字滤波器的函数介绍

MATLAB下的数字信号处理实现示例一信号、系统和系统响应 1、理想采样信号序列（1）首先产生信号x(n),0<=n<=50 n=0:50; %定义序列的长度是50 A=444.128; %设置信号有关的参数 a=50*sqrt(2.0)*pi; T=0.001; %采样率 w0=50*sqrt(2.0)*pi; x=A*exp(-a*n*T).*sin(w0*n*T); %pi是MATLAB定义的π，信号乘可采用“.*”close all %清除已经绘制的x(n)图形 subplot(3,1,1);stem(x); %绘制x(n)的图形 title(‘理想采样信号序列’); （2）绘制信号x(n)的幅度谱和相位谱 k=-25:25; W=(pi/12.5)*k; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X); %绘制x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’); angX=angle(X); %绘制x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)（3）改变参数为：1,0734.2,4.0,10==Ω==TAα n=0:50; %定义序列的长度是50 A=1; %设置信号有关的参数 a=0.4; T=1; %采样率 w0=2.0734; x=A*exp(-a*n*T).*sin(w0*n*T); %pi是MATLAB定义的π，信号乘可采用“.*”close all %清除已经绘制的x(n)图形 subplot(3,1,1);stem(x); %绘制x(n)的图形 title(‘理想采样信号序列’); k=-25:25; W=(pi/12.5)*k; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X); %绘制x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’); angX=angle(X); %绘制x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’) 2、单位脉冲序列在MatLab中，这一函数可以用zeros函数实现： n=1:50; %定义序列的长度是50 x=zeros(1,50); %注意：MATLAB中数组下标从1开始

基于matlab的数字滤波器设计

淮北煤炭师范学院 2009届学士学位论文基于MA TLAB的数字滤波器设计学院、专业物理与电子信息学院电子信息科学与技术研究方向基于MATLAB的数字滤波器设计学生姓名耿博学号200513432024 指导教师姓名邹锋指导教师职称讲师 2009 年4 月18

基于MATLAB的数字滤波器设计耿博（淮北煤炭师范学院物理与电子信息学院235000）摘要随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。数字滤波是数字信号处理的重要内容，数字滤波器可分为IIR和FIR两大类。对于IIR数字滤波器的设计，需要借助模拟原型滤波器，再将模拟滤波器转化为数字滤波器，文中采用的设计方法是脉冲响应不变法、双向性变换法和完全函数设计法；对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，文中采用的设计方法是窗函数法。本文根据IIR滤波器和FIR滤波器的特点，在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并对采集的语音信号进行分析，最后给出了IIR和FIR对语音滤波的效果。关键词数字滤波器；IIR ；FIR ；MATLAB

The Design of Digital Filter based on MATLAB Geng Bo School of Physics and Electronics Information, Huaibei Coal Industry Teachers? College, 235000 ABSTRACT Along with the information age and the digital world arrival, the digital signal processing has become a now extremely important discipline and the area of technology.The digital signal processing in the correspondence, the multitudinous domains the pronunciation such as the image, the automatic control, the radar, the military, the aerospace, the medical service and the domestic electric appliances and so on have obtained the widespread application.In the digital signal processing application, the digital filter are extremely important and have obtained the widespread application. The digital filter are the digital signal processing important content, the digital filter may divide into IIR and the FIR two main kinds. As for the IIR digital filter design, we need the help of analog prototype filter, and then transform analog filter into digital filter. In the paper we use the design of the pulse response invariable method, the bilinear method and full function design; as for the FIR filter, we can design it directly based on the giving frequency, in the paper it uses the design of the window function.This article according to the IIR filter and the FIR filter characteristic, uses the bilinearity method of transformation under the MATLAB bad boundary to design IIR and to design the FIR numeral filter separately with the window box number, and carries on the analysis to the gathering pronunciation signal, and finally gives IIR and FIR to the pronunciation filter effect. Keywords Digtial Filter;IIR;FIR;MATLAB

数字滤波器的MATLAB设计与DSP上的实现

数字滤波器的MAT LAB设计与 DSP上的实现数字滤波器的MATLAB 设计与DSP上的实现公文易文秘资源网佚名2007-11-15 11:56:42我要投稿添加到百度搜藏摘要：以窗函数法设计线性相位FIR数字滤波器为例，介绍用MATLAB工具软件设计数字滤波器的方法和在定点DSP上的实现。实现时，先在CCS5000仿真开发，然后将程序加载到TMS320VC5409评估板上实时运行，结果实现了目标要求。文中还讨论了定标、误差、循环寻址等在DSP上实现的关键问题。关键词摘要：以窗函数法设计线性相位 FIR数字滤波器为例，介绍用 MATLAB工具软件设计数字滤波器的方法和在定点DSP上的实现。实现时，先在 CCS5000仿真开发，然后将程序加载到 TMS320VC5 409评估板上实时运行，结果实现了目标要求。文中还讨论了定标、误差、循环寻址等在DSP上实现的关键问题。关键词：数字滤波器MATLAB DSP 引言随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应

用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 1数字滤波器的设计 1.1数字滤波器设计的基本步骤数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR ）滤波器和有限长冲激响应（FIR ）滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。随着 MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。数字滤波器设计的基本步骤如下： (1确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给岀幅度和相位响应。幅度指标主要以两种方式给岀。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于FI R滤波器的设计。第二种指标是相对指标。它以分贝值的形式给岀要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中人有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点：①只包含实数算法，不涉及复数运算；②不存在延迟失真，只有固定数量的延迟；③长度为N的滤波器(阶数为N-1)，计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。 (2)逼近

数字滤波器matlab的程序

数字滤波器matlab的源代码 function lvbo(Ua,Ub,choise) %参考指令：lvbo(2*pi,10*pi,1/0/-1) U1=min(Ua,Ub); U2=max(Ua,Ub); Us=16*U2; T=2*pi/Us; T_sum=4*max(2*pi/Ua,2*pi/Ub); sum=T_sum/T; t=T:T:T_sum; x=sin(U1*t)+0.8*sin(U2*t); X=DFT(x); figure(1); subplot(221) U=Us/sum:Us/sum:Us; stem(U,abs(X));grid on axis([Us/sum,Us/2,0,1.2*max(abs(X))]) title('原模拟信号采样频谱图') Ucd=U1+(U2-U1)*1/5;Usd=U2-(U2-U1)*1/5; switch choise case 1 Hz_ejw=IIR_DF_BW(Ucd,1,Usd,30,T,sum); case -1 Hz_ejw=IIR_DF_CF(Ucd,1,Usd,30,T,sum); case 0 Hz_ejw=FIR_DF_HM(U1,U2,T,sum); otherwise Hz_ejw=IIR_DF_BW(Ucd,1,Usd,30,T,sum); end Y=X.*Hz_ejw; y=1/sum*conj(DFT(conj(Y))); figure(1); subplot(224) plot(t,real(y)); title('模拟信号滤波后');grid on axis([0,T_sum,-max(real(y))*1.5,max(real(y))*1.5]) subplot(222); plot(t,x); hold on

基于矩形窗设计FIR数字滤波器

郑州航空工业管理学院《电子信息系统仿真》课程设计 14级专业班级题目基于矩形窗设计FIR数字滤波器姓名学号二О一六年十一月二十五日

第一章 FIR 滤波器的设计原理及方法 FIR 滤波器通常采用窗函数方法来设计。窗设计的基本思想是，首先选择一个适当的理想选频滤波器（它总是具有一个非因果，无限持续时间脉冲响应），然后街区（加窗）它的脉冲响应得到线性相位和因果FIR 滤波器。我们用Hd （e^jw ）表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽wc

实验11 用MATLAB设计FIR数字滤波器

实验11 用MATLAB 设计FIR 数字滤波器一、实验目的： 1、加深对窗函数法设计FIR 数字滤波器的基本原理的理解。 2、学习用MA TLAB 语言的窗函数法编写设计FIR 数字滤波器的程序。 3、了解MATLAB 语言有关窗函数法设计FIR 数字滤波器的常用函数用法。二、实验内容及步骤 2、选择合适的窗函数设计FIR 数字低通滤波器，要求： w p =0.2π，R p =0.05dB ； w s =0.3π，A s =40dB 。描绘该滤波器的脉冲响应、窗函数及滤波器的幅频响应曲线和相频响应曲线。分析：根据设计指标要求，并查表11-1，选择汉宁窗。程序清单如下： function hd=ideal_lp(wc,N) wp=0.2*pi;ws=0.3*pi;deltaw=ws-wp; tao=(N-1)/2; n=[0:(N-1)]; m=n-tao+eps; hd=sin(wc*m)./(pi*m); function[db,mag,pha,grd,w]=freqz_m(b,a); [H,w]=freqz(b,a,1000,'whole'); H=(H(1:501))';w=(w(1:501))'; mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); grd=grpdelay(b,a,w); wp=0.2*pi;ws=0.3*pi;deltaw=ws-wp; wc=(ws+wp)/2; 课程名称：数字信号处理实验成绩：指导教师：实验报告院系：信息工程学院班级：电信二班学号：姓名：日期：

matlab数字滤波器设计程序

%要求设计一butterworth低通数字滤波器，wp=30hz，ws=40hz，rp=0.5，rs=40，fs=100hz。>>wp=30;ws=40;rp=0.5;rs=40;fs=100; >>wp=30*2*pi;ws=40*2*pi; >> [n,wn]=buttord(wp,ws,rp,rs,'s'); >> [z,p,k]=buttap(n); >> [num,den]=zp2tf(z,p,k); >> [num1,den1]=impinvar(num,den); Warning: The output is not correct/robust. Coeffs of B(s)/A(s) are real, but B(z)/A(z) has complex coeffs. Probable cause is rooting of high-order repeated poles in A(s). > In impinvar at 124 >> [num2,den2]=bilinear(num,den,100); >> [h,w]=freqz(num1,den1); >> [h1,w1]=freqz(num2,den2); >>subplot(1,2,1); >>plot(w*fs/(2*pi),abs(h)); >>subplot(1,2,2); >>plot(w1*fs/(2*pi),abs(h1)); >>figure(1); >>subplot(1,2,1); >>zplane(num1,den1); >>subplot(1,2,2); >>zplane(num2,den2);

高级数字滤波器设计及Matlab实现

高级数字滤波器设计及Matlab 实现利用Parks-McClellan 算法设计线性相位FIR 滤波器一、算法原理长度为2n+1的线性相位数字滤波器的传输函数为：20 ()n k k k G z h Z -== ∑。当 Z=exp(j T ω)=exp(2j F π)时，可得到频率响应： ()exp(2)cos 2n k k G F j nF d k F ππ==-∑ exp(2)()j F H F π=- 其中2n k k d h -=，k=0，……，n-1，0n d h = max ()()()W F D F H F - 对于一个理想的低通滤波器上式中的H(F)可以表示为：误差加权函数：允许设计者自己给定通带和阻带内的误差范围。令p s B B A = ，设计长度为2n+1的线性相位低通滤波器只要找到k d 并使得m a x ()()()W F D F H F -最小。设空间A 是[0,1/2]的封闭子空间，为了使0 ()cos 2n k k H F d kF π== ∑是D(F) 在A 上唯一的最佳逼近，加权误差方程()()[()()]E F W F D F H F =-在A 上至少要有n+2个交错点。因此1()()i i E F E F E -=-=±，011,n i F F F F A +<∈ ，

max ()E E F =。算法的流程如图1所示。对于给定的n+2个频率点，需要计算n+2个方程： ()(()())(1)k k k k W F H F D F ρ-=-- 写成矩阵的形式就是：图 1.

0000001 1 1 1 1011 1 1 1 11 1cos 2cos 4cos 2()()11cos 2cos 4cos 2()()()(1)1cos 2cos 4cos 2()n n n n n n n F F nF W F d D F d F F nF D F W F d D F F F nF W F ππππππρπππ++++++?? ???? ?? ?? ??-????????????=??????????????????-?? ????? ? 通过该方程组可得：其中：利用拉格朗日插值公式可得：这里利用求得的H(F)求出误差函数E(F)。如果对所有的频率都有()E F ρ≤，说明ρ是纹波极值，交错频率点121,n F F F + 是交错频率点。若存在某些频率使得()E F ρ>，说明初始交错点组中的某些点需要交换。对于上次确定的121,n F F F + 中每一点，都检查其附近是否存在某一频率 ()E F ρ>，如果有再在该点附近找出局部极值点，并用该，点代替原来的

基于matlab的数字滤波器的设计与仿真开题报告.pdf

一、选题背景、目的和意义背景信号与信号处理时信息科学中近几十年来发展最为迅速的学科之一。长期以来?信号处理技术一直用于转换、产生模拟或数字信号?其中最为频繁应用的领域就是信号的滤波。数字滤波是语音、图像处理、模式识别和谱分析等应用中的一个基本处理部件?它可以满足滤波器对幅度和相位特性的严格要求?避免模拟滤波器无法克服的电压漂移、温度漂移和噪声等问题。目的我的毕业设计的课题是《基于matlab的数字滤波器的设计与仿真》?其主要目的是通过此次课程设计进一步学习和巩固数字信号处理及其相关知识?并学会利用所学的知识能?在设计过程中能综合运用所学知识内容?进一步熟悉和掌握matlab的使用方法;对数字滤波器的原理有较深的了解;为即将进入社会参加工作打下坚实的基础;掌握收集资料、消化资料和综合资料的能力等等。意义从事电子通信业而不能熟练操作使用matlab电子线路设计软件?在工作和学习中将是寸步难行的。在数学、电子、金融等行业?使用matlab等计算机软件对产品进行设计、仿真在很早以前就已经成为了一种趋势?这类软件的问世也极大地提高了设计人员在通信、电子等行业的产品设计质量与效率。众所周知?实际过程中信号传输都要经过调制与解调这一过程?由于消息传过来的原始信号即调制信号具有频谱较低的频谱分量?这种信号在许多信道中不宜传输。因而?在通信系统的发送端通常需要有调制过程?反之在接收端则需要有解调过程。二、国内外的研究综述数字滤波在DSP?数字信号处理?中占有重要地位。数字滤波器按实现的网络结构或者从单位脉冲响应?分为IIR(无限脉冲响应?和FIR?有限脉冲响应?滤波器。如果IIR滤波器和FIR滤波器具有相同的性能?那么通常IIR滤波器可以用较低的阶数获得高的选择性?执行速度更快?所有的储存单元更少?所有既经济又高效。

FIR数字滤波器设计论文

毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名：日期：

学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名：日期：年月日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期：年月日导师签名：日期：年月日

基于matlab的FIR数字滤波器设计(多通带,窗函数法)

数字信号处理课程设计报告设计名称：基于matlab的FIR数字滤波器设计彪

一、课程设计的目的 1、通过课程设计把自己在大学中所学的知识应用到实践当中。 2、深入了解利用Matlab设计FIR数字滤波器的基本方法。 3、在课程设计的过程中掌握程序编译及软件设计的基本方法。 4、提高自己对于新知识的学习能力及进行实际操作的能力。 5、锻炼自己通过网络及各种资料解决实际问题的能力。二、主要设计内容利用窗函数法设计FIR滤波器，绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。三、设计原理 FIR 滤波器具有严格的相位特性，对于信号处理和数据传输是很重要的。目前 FIR滤波器的设计方法主要有三种：窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求高的时候是比较灵活方便的。如果 FIR 滤波器的 h(n)为实数, 而且满足以下任意条件,滤波器就具有准确的线性相位：第一种：偶对称，h(n)=h(N-1-n)，φ (ω)=-(N-1)ω/2 第二种：奇对称，h(n)=-h(N-1-n), φ(ω)=-(N-1)ω/2+pi/2 对称中心在n=(N-1)/2处四、设计步骤 1.设计滤波器 2.所设计的滤波器对多个频带叠加的正弦信号进行处理 3.比较滤波前后信号的波形及频谱五、用窗函数设FIR 滤波器的基本方法基本思路：从时域出发设计 h(n)逼近理想 hd(n)。设理想滤波器的单位响应在时域表达为hd(n),则Hd(n) 一般是无限长的，且是非因果的，不能

FIR数字滤波器设计及MATLAB使用要点

数字信号处理课程设计《数字信号处理》课程设计报告 FIR数字滤波器设计及MATLAB实现专业：通信工程班级：通信1101班组次：第9组姓名及学号：姓名及学号：

目录一、设计目的 (3) 二、设计任务 (3) 三、设计原理 (3) 3.1窗函数法 (3) 3.2频率采样法 (4) 3.3最优化设计 (5) 3.3.1等波纹切比雪夫逼近准则 (5) 3.3.2仿真函数 (6) 四、设计过程 (7) 五、收获与体会 (13) 参考文献 (13)

FIR数字滤波器设计及MATLAB实现一、设计目的 FIR滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。滤波器设计是根据给定滤波器的频率特性，求得满足该特性的传输函数。二、设计任务 FIR滤波器设计的任务是选择有限长度的() H e满足一定 h n，使传输函数()jw 的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个基本问题：（1）根据实际要求确定数字滤波器性能指标；（2）用一个因果稳定的系统函数去逼近这个理想性能指标；（3）用一个有限精度的运算去实现这个传输函数。三、设计原理 FIR滤波器设计的任务是选择有限长度的() H e满足一定 h n，使传输函数()jw 的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个基本问题：（1）根据实际要求确定数字滤波器性能指标；（2）用一个因果稳定的系统函数去逼近这个理想性能指标；（3）用一个有限精度的运算去实现这个传输函数。 3.1窗函数法设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应()jw H e,设计 d

自适应滤波器的设计开题报告

自适应滤波器的设计与应用一、题目来源来源于其他二、研究目的和意义滤波技术在当今信息处理领域中有着极其重要的应用。滤波是从连续的或离散的输入数据中除去噪音和干扰以提取有用信息的过程，相应的装置就称为滤波器。滤波器实际上是一种选频系统，他对某些频率的信号予以很小的衰减，使该部分信号顺利通过。而对其他不需要的频率信号予以很大的衰减，尽可能阻止这些信号通过。滤波器研究的一个目的就是：如何设计和制造最佳的(或最优的)滤波器。在数字信号处理中，数字滤波是语音和图像处理、模式识别、频谱分析等应用中的一个基本处理算法。在许多应用场合，由于无法预先知道信号和噪声的特性或者它们是随时间变化的，仅仅用 FIR 和 IIR两种具有固定滤波系数的滤波器无法实现最优滤波。在这种情况下，必须设计自适应滤波器，以跟踪信号和噪声的变化。自适应滤波器是利用前一时刻已获得的滤波器参数，自动地调节、更新现时刻的滤波器参数，以适应信号和噪声未知的统计特性，从而实现最优滤波。当在未知统计特性的环境下处理观测信号时，利用自适应滤波器可以获得令人满意的效果，其性能远超过通用方法所设计的固定参数滤波器。

三、阅读的主要参考文献及资料名称 1、《数字信号处理》刘益成（第二版）西安电子科技出版社 2、《数字信号处理》张小虹（第二版）机械工业出版社 3、自适应信号处理[M]．西安：西安电子科技大学出版社，2001． 4．邹理和,数字信号处理, 国防工业出版社,1985 5．丁玉美等, 数字信号处理,西安电子科技大学出版社,1999 6．程佩青, 数字信号处理,清华大学出版社,2001 7. The MathWorks Inc, Signal Processing Toolbox For Use with MATLAB, Sept. 2000 8. vinay K.Ingle, John G.Proakis,数字信号处理及MATLAB实现,陈怀琛等译,电子工业出版社,1998.9 9、《MATLAB编程参考手册》 10、中国期刊网的相关文献 11、赫金，自适应滤波器原理第四版，西安工业出版社，2010-5-1 四、国内外现状和发展趋势与主攻方向自适应滤波器的理论与技术是50年代末和60年代初发展起来的。它是现代信号处理技术的重要组成部分，对复杂信号的处理具有独特的功能。自适应滤波器在数字滤波器中试属于随机数字信号处理的范畴。对于随机数字信号的滤波处理，通常有维纳滤波，卡尔曼滤波和自适应滤波，维纳滤波的权系数是固定的，适用于平稳随机信号；卡尔曼滤波器的权系数是可变的，适用于非平稳随机信号中。但是，只有在对信号和噪声的统计特性先验

基于MATLAB的数字滤波器设计

基于matlab的数字滤波器设计摘要：本文介绍的是数字滤波器在MATLAB环境下的设计方法。数字滤波是数字信号处理的重要内容，在实际应用中有非常大的作用。我们研究的数字滤波器可分为IIR和FIR两大类。对于IIR数字滤波器的设计，我们需要借助模拟原型滤波器，然后再将模拟滤波器转化为数字滤波器，文中采用的设计方法是脉冲响应不变法、双向性变换法和完全函数设计法；对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，文中采用的设计方法是窗函数法。根据IIR 滤波器和FIR滤波器的特点，本文在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并让这两种滤波器对采集的语音信号进行分析和比较，经过分析，最后给出了IIR和FIR对语音滤波的效果，并总结这两种滤波器在MATLAB环境下设计方法的优缺点。关键词：数字滤波器；IIR ；FIR ；MATLAB

The Design of Digital Filter based on MATLAB Abstract：This article describes a digital filter in the MATLAB environment design. Digital filtering is an important part of digital signal processing which is playing a very big role in practice .The digital filter we studied can be divided into two categories——IIR and FIR. For the IIR digital filter design, we will need the help of simulation prototype filter, analog filters and then converted it into digital filter For the IIR digital filter design .The design methods used in the text is the same impulse response method, bi-sexual transformation and full function design ;We can based on the frequency characteristics of the given direct design, design method used in the text is the window function for FIR digital filter design. Based on the characteristics of IIR filter and FIR filters ,the bad paper in the MATLAB environment under the bilinear transformation method were used to design IIR and FIR with window function digital filter design and filters to capture both the voice signal analysis and compare. Through analysis of IIR and FIR Finally, the effect of filtering on the speech, and concluded the advantages and disadvantages in the two filter design methods in the MATLAB environment. Key words: Digital Filter ；IIR；FIR；MATLAB

基于Matlab的FIR低通滤波器设计[开题报告]2011-11-07

电子信息学院本科毕业设计(论文)开题报告论文题目基于Matlab的FIR低通滤波器设计学生姓名专业班级性别一、课题研究意义及现状所谓数字信号处理技术是指运用数值计算的方法对信号进行分析、变换、综合、估值与识别等的处理。随着电子技术和电子计算机技术的发展，数字信号处理技术受到了越来越广泛的关注，数字信号处理的理论和技术也在不断丰富和完善，新的理论和新技术层出不穷。可以说，数字信号处理是发展最快、应用最广泛、成效最显著的新科学之一，目前已广泛地应用在语音、雷达、声纳、地震、图像、通信、控制、生物医学、遥感遥测、地质勘探、航空航天、故障检测、自动化仪表等领域。数字滤波器是数字信号处理中的重要组成部分之一，在数字信号处理中，数字滤波占有极其重要的地位。滤波器可以分为IIR滤波器和FIR滤波器两大类，IIR滤波器并不能得到严格的线性相位特性，因此在许多实际应用中为了得到线性相位特征，还必须另外增加相位校正网络，使滤波器设计变得复杂，成本也高，而FIR滤波器在保证幅度特性满足技术要求的同时，还可以很容易的做到有严格的线性相位特性。它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。目前对数字滤波器的设计有多种方法，其中著名的Matlab软件包，功能强大、使用方便。传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。利用Matlab信号处理工具箱可以快速有效的设计由软件组成的常规数字滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。本课题在基于Matlab的FIR低通滤波器的设计方法、编程以及实际应用上有着积极的意义。

基于MATLAB的数字滤波器设计(DOC)

目录 1 引言 (1) 2 设计任务 (2) 2.1设计内容 (2) 2.2设计要求 (2) 3 语音信号的采集及时频分析 (3) 3.1语音信号的采集 (3) 3.2语音信号的时频分析 (3) 4 基于MATLAB的数字滤波器的设计 (5) 4.1数字滤波器的设计 (5) 4.1.1数字滤波器的基本概念 (5) 4.1.2 IIR滤波器设计思想 (5) 4.2IIR数字滤波器设计 (5) 4.2.1 IIR低通滤波器设计 (5) 4.2.2 IIR带通滤波器设计 (7) 4.2.3 IIR带通滤波器设计 (9) 5 合成信号及其滤波 (12) 5.1合成信号 (12) 5.2合成信号滤波 (13) 6 设计系统界面 (15) 6.1系统界面设计工具—GUI概述 (15) 6.2界面设计及使用说明 (15) 7 心得体会 (18) 参考文献 (19) 附录 (20)

1 引言数字信号处理是利用计算机或专用处理设备，以数值计算的方法对信号进行采集、抽样、变换、综合、估值与识别等加工处理，借以达到提取信息和便于应用的目的。它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等优点。数字滤波器, 是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应( FIR，Finite Impulse Response)滤波器和无限冲激响应( IIR，Infinite Impulse Response)滤波器。 IIR滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。同时，IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。 MATALB 可以创建图形用户界面GUI (GraphicalUser Interface) ,它是用户和计算机之间交流的工具。MATLAB 将所有GUl 支持的用户控件都集成在这个环境中并提供界面外观、属性和行为响应方式的设置方法,随着版本的提高,这种能力还会不断加强。而且具有强大的绘图功能,可以轻松的获得更高质量的曲线图。滤波器的设计可以通过软件或设计专用的硬件两种方式来实现。随着MATLAB软件及信号处理工具箱的不断完善，MATLAB很快成为应用学科等领域不可或缺的基础软件。它可以快速有效地实现数字滤波器的设计、分析和仿真，极大地减轻了工作量,有利于滤波器设计的最优化。