Matlab数字滤波器设计

1 基于Matlab的FIR数字滤波器设计1郑伟鑫 ( 1. 华南理工大学机械与汽车工程学院机械电子工程1班 ) 摘要：数字滤波器不是用物理元器件构成，而是按照某种算法编写一段程序或模块，对数字信号进行处理加工，从而达到滤波的要求。数字滤波器是实现滤波过程的一种数字信号处理系统，具有离散时间系统的基本特征。鉴于滤波器在数字信号处理中的广泛应用以及在Matlab环境下设计滤波器的快捷方便等特性，本文研究了Matlab环境下，数字滤波器的设计与实现。数字滤波器包括IIR（无限脉冲响应）数字滤波器和FIR（有限脉冲响应）数字滤波器。其实现方式主要有软件方式实现（即Matlab编程实现）和硬件方式实现。本文主要介绍FIR数字滤波器的软硬件方式实现。Matlab强大的数字信号处理工具箱为我们提供了易用、有力的设计工具。通过正确的使用函数，很容易设计出性能良好的数字滤波器。关键词：数字滤波器 Matlab FIR Design of FIR digital filter based on Matlab 1 Zheng Weixin(1.South China University of Technology, Mechanical and Automotive institute, Class 1 of Mechanical and Electronic Engineering)：Abstract Digital filter is not composed of physical components, but writing a program or module according to some algorithms. The digital signal processing is aimed to achieve the filtering requirements. The digital filter is a digital signal processing system implemented in filtering process, which is also the basic feature of the discrete time system. In view of the extensive applications of filters in digital signal processing and the fast and convenient features of filter

design in the Matlab environment, the paper mainly introduces the design and implementation of digital filters under it. The digital filter includes IIR (infinite impulse responsibility) digital filter and FIR (finite impulse responsibility) digital filter. There’re two methods to achieve it, including the software (i.e. programming by Matlab) and the hardware mode. This paper mainly introduces the software and the hardware achievement of FIR digital filter. The powerful digital signal processing toolbox in Matlab provides us with a powerful and easy design tool. By using the function correctly, we are easy to design the digital filter with good performance. ：Matlab FIR Key

words Digital filter 高、稳定、体积小、重量轻、灵活、不要求阻抗匹0 前言配以及可以避免模拟滤波器所无法克服的电压漂移和噪声问题。

数字滤波器的设计，其实质是数学逼数字信号处理在通信、语音、图像,

自动控制、近理论的应用，通过计算使物理可实现的实际滤波雷达、军事、航空航天、医疗和家用电器等众多领器频率特性逼近理想的或给定的频率特性，以达到域得到了广泛的应用。在数字信号处理中，滤波器去除干扰提取有用信号的目的。占有非常重要的地位。数字滤波器在语音和图像处Matlab是由美国MathWorks公司于1982 年推出理、HDTV、模式识别、频谱分析等方面得到广泛的一套高性能的数值计算和可视化软件,它集数值应用。分析、矩阵运算、信号处理和图形显示于一体。在所谓滤波器，是指对输入信号进行滤波的软件或硬件。所谓数字滤波器，是指输入、输出均为数Matlab内部配备了涉及到自动控制、信号处理和计字信号，通过一定运算关系改变输入信号所含频率算机仿真等种类繁

多的工具箱,所以Matlab的应用成分的相对比例或者滤除某些频率成分的

硬件或软非常广泛,它可涉足于数值分析、控制、信号分析和件。由于数

字滤波信号形式与实现滤波方法与模拟通信等多种领域。Matlab不仅可

完成基本代数运算滤波方法不同，数字滤波器具有比模拟滤波器精度

2 同之处，IIR滤波器设计方法有两类，经常用到的操作,而且还可完成矩阵函数运算,提供丰富的实用一类设计方法是借助于模拟滤波器的设计方法进行函数命令。另外,用户还可以根据自己的需要编写函的。其设计思路是：先设计模拟滤波器得到传输函数。数H(s),然后将H(s)按某种方法转换成数字滤波器的系统函数H(Z)。这类方法是基于模拟滤波器的设计 1 数字示波器方法相对比较成熟。它不仅有完整的设计公式，也有完整的图表供查阅。更可以直接调用MATLAB 1.1 数字示波器的分类中的对应的函数进行设计。另一类是直接在频域或数字滤波器从功能上分类：可分为低通滤波者时域中进行设计的，设计时必须用计算机作辅助器、高通滤波器、带通滤波器、带阻滤波器。设计，直接调用MatLab中的一些程序或者函数可从滤波器的网络结构或者从单位脉冲响应分以很方便地设计出所需要的滤波器。类：可分为IIR滤波器(即无限长单位冲激响应滤波FIR滤波器不能采用由模拟滤波器的设计进行器)和FIR滤波器(即有限长单位冲激响应滤波器)。转换的方法。经常用的是窗函数法和频率采样法。IIR数字滤波器和FIR数字滤波

器各有优缺点，也也可以借助计算机辅助设计软件采用切

比雪夫等波各有不同的应用场合。逼近法进行设计。下面

主要介绍FIR数字滤波器的1.2 数字滤波器的表示方法软

件及硬件方式实现。数字滤波器包括以下四种表示方法：

(1) 差分方式表示MN 2 FIR数字滤波器的软件设计

其中输入信号为x(t),

y(n) ax(n i) by(n i)ii i 0i 1

输出信号为y(t) ,a,b为常数。FIR滤波器在保证幅度特性满

足技术要求的同ii(2) 利用传递函数表示时，很容易做到有

严格的线性相位特性。设FIR滤M波器单位脉冲响应h(n)

长度为N，其系统函数为

rbzr1N r 0 nH(z) h(n)z H(z)

是的（N-1）次多项式，在Z平面上

Nn 0 k1 azH(z)

有（N-1）kk 1个零点，原点z=0是（N-1）阶重极点。因

此，H(z)在差分方程的基础上，方程取Z变换，化简即可得

永远稳定。稳定和线性相位特性是FIR滤波器突出上式。

的优点。(3) 按单位响应冲击响应 2.1 窗函数法设计FIR

数字滤波器假设输入单位冲击时间序列的滤波器输出序列

jH(e)理想低通滤波器的传输函数为 为：，则任

意输入信号时的数字滤h(n) T*(n)+波器输出为： -ja

e, c

j H(e)= y(n) x(k)h(n k) x(n)*h(n)

0, k

c x(k n)h(n) h(n)*x(n)相应的单位取样响应h(n)为

-jajnc(4) 按频率特性表示k sin((n ))1

c h(n)=ee

d - 2(n a)c频率特性是数字滤波器对正弦输入序列的响应，利用频率特性与传递函数之间的关系可得：可以看出，理想低通滤波器的单位取样响应

H(j) H(z) H(jT)e限长，且为非因果序是无 j()

列。为了构造一个长度为N的sjwTsz e 其中，H(jwTs)和频率f构成幅频特性，ψ(w)线性相位滤波器，只有将h(n)截取一段，并保证截和频率f构成相频特性。取的一段对（N-1）/2对称。这就是截断函数的一 1.3 数字滤波器设计方法概述 j H(e)个例子。我们知道是一个以2为周期的函IIR滤波器和FIR滤波器的设计方法有很多不

2 通滤波器的频率响应，带通滤波器的通带频率(归一数，可以展为傅氏级数。即 化后的频率)为0.25到0.60，阶数为38。 j n j H(e)=h(n)e，设定截止频率和阶数：n=- 设计FIR滤波器就是根据要求找到有限个傅氏Wn=[0.25,0.60]; n=38; 级数系数，来代替并近似无限项傅氏级数，这样一使有默认汉明窗的带通滤波器: 来，在频率不连续点附近会产生误差，即截断效应。b1=fir1(n, Wn); 截断效应，指的是采取截断函数时，截取的有使用切比雪夫窗函数的带通

滤波器，先输入切限长信号不能完全反映原信号的频率特性。具体地比雪夫窗函数: 说,会增加新的频率成分,并且使谱值大小发生变window=chebwin(n+1,40); 化,这种现象称为频率泄漏。从能量角度来讲,频b2=fir1(n, Wn, window); 率泄漏现象相当于原信号各种频率成分处的能量渗从图1可以看出，汉明窗函数具有较低的旁瓣透到其他频率成分上,所以又称为功率泄漏。泄漏幅度，尤其是第一旁瓣幅度；而切比雪夫窗函数具是由于无限长信号的突然截断造成的。因此,自然有比较窄的主瓣，相比之下，过渡带也较窄。从分想到,如果能改变这种突然截断方式,泄漏会得到析结果来看，汉明窗函数适用于对阻带衰减有要求改善。选择适当的窗函数,对所取样本函数进行不的场合，而切比雪夫窗适用于要求较窄的过渡带的等权处理,便是一种有效的措施。需要注意的是，场合。在使用窗函数法时，选取傅氏级数的基数愈多，引起的误差就愈小，但同时项数增多也使成本、体积增加。 2.2 Matlab中的窗函数法软件实现Matlab中，函数fir1和fir2就是基于窗函

数方法的。firl函数实现了加窗线性相位FIR 数字滤波器设计的经典方法。主要用于常用的标准通带滤波器设计，包括:低通、带通、高通和带阻数字滤波器。 b=fir1(n, Wn)可得到n阶低通FIR滤波器，滤图1 用汉明窗和切比雪夫窗设计的FIR滤波器波器系数包括在b中，它是一个截止频率为Wn的函数fir2()也设计加窗的FIR滤波器，但它针对加汉明窗的线性相位FIR滤波器。如果Wn是一个任意形状的分段线性频率响应，这一点在fir1中是包含两个元素的向量，Wn=[W1 W2],firl返回一个n受到限制的。阶的带通滤波器，其通带为W1

率的一半。它们必须按递增的顺序从0至‘DC-1’)，使第一频带为通带；b=fir1(n,Wn,‘DC-0’)1结束。缺省情况下，函数fir2()默认使用的是使第一频带为阻带。Hamming窗，可在后面参数中指定其它的窗函数。对于在Fs/2附近为通带的滤波器如高通或带阻关于利用函数fir2()设计的FIR滤波器这里不再敷滤波器，n必须是偶数。述。缺省情况下，firl()使用Hamming窗。可以在参数window中指定其它窗，包括矩形窗、Harming窗、Bartlett窗、Blackman窗、Kaise窗等。 3 FIR数字滤波器的硬件设计缺省情况下，滤波器被归一化，以使经加窗后的第一个通带的中心幅值刚好是1。使用参数’在Matlab的toolbox中有现成的滤波器模块，noscale’可以阻止这样做。调用时只需设置相应参数，便可使用，十分方便，如图1是分别用切比雪夫窗和汉明窗设计的带是一种很实用的数字滤波硬件设计。

2 在Matlab中Simulink命令下有Digital Filter 0.01，采用凯塞窗设计。Design模块可以设计用户所需要的数字滤波器。通通过simulink绘制如图3所示的正弦波形框过fdatool 工具将所设计好的封装导入到Digital 图，输入信号为 Filter Design模块即可

n(30t) sin(60t) 3.1 FDATOOL工具箱观察滤波结果。 f sin(10t) si

FDATOOL是Matlab软件中提供滤波器设计的专门工具箱，可以设计几乎所有的常规滤波器，其界面如图2所示。图3 FIR数字滤波器仿真框图设计窗函数幅频

响应如图4所示，图2 fdatool界面界面上半部分是特性区，用来显示滤波器的各种特性。界面下半部分是Design Filter，用来设置滤波器的设计参数，主要分为：ResoponseType（滤波器类型）选项、Design Method（设计方法）选项、Filter Order （滤波器阶数）选项、Frenquency Specifications（频率特性）选项和Magnitude Specifications（幅度特性）选项。图4 kaiser窗函数幅频响应3.2 Simulink简介将封装好的滤波器文件kaiser.fda导入Digital Filter Design模块中，可得如下波形比较：Simulink是一个进行动态系统建模、仿真和综合分析的集成软件包，它可以处理线性、非线性系统；离散系统、连续系统以及混合系统；单任务、多任务离散事件系统等，目前已经覆盖通信、控制、信号处理、电力系统等诸多领域。

Simulink和Matlab的良好结合使得用户可以利用Matlab丰富的资源，建立仿真模型，监控仿真过程，并且可以在线修改参数，并观察改变后的结果。 3.3 FIR带通滤波器设计此处举例一个硬件实现带通滤波器的电路及其实现。针对含有5Hz、15Hz 和30Hz的混合正弦波信号，设计一个FIR带通滤波器。参数要求：采样频率Fs=100Hz，通带下限截止频率Fc1=10Hz，通带上限截止频率Fc2=20Hz，过渡带宽6Hz，通阻带波动

2

图5 导入窗函数后的波形与之前比较图5中横坐标表示时间（t），纵坐标表示振幅。滤波后，信号周期约为0.067s，即频率为15Hz，通过对滤波前后的离散波形进行对比可以很明显的看出，复合信号通过设计的带通滤波器后，所需要的有用信号（15Hz）分离了出来，基本达到预期目的。

5 总结参考文献 Matlab是一个强大的数学计算/仿真工具，它[1]丁玉美，高西全.数字信号处理（第二版）. 西自带的信号处理工具箱（signal processing 安:西安电子科技大学出版社,2000：151-317 toolbox）可以在数字信号处理领域很好的设计和处[2]黄建明，

隋燕. 基于数字滤波器设计的讨论. 理数字滤波，为数字滤波，信号处理提供了有力，南

开大学学报(自然科学版),2003,26(3):15-18 易用，强大的工具。[3]李茂，任学恒.基于Matlab/Simulink的数字本文研究了基于MATLAB下的FIR滤波器的设计滤波器的设计

与仿真.磁性材料与器实现、进行了数字信号

处理、数字图像处理中的滤件.2007.02:65-68

波器设计。在设计过程中，Matlab强大性，

易用性[4]聂祥飞.基于Matlab的IIR数字滤

波器设计得到了充分的体现。从本次设计中

我学会了独立完研究.信息技术.2002.26(12):38-40 成任务，这次设计中学

会了很多东西包括软件编程，[5]Robert J. Schilling, Sandra L. Harris. 硬件设计等有效的滤

波设计方案。但是对于Matlab Fundamentals

of Digital Signal Processing Using 的强大功能还

只是点到为止，以后将会深入了解MATLAB （影印版）西安:西安交通大学出版社 Matlab

这个强大的编程软件。2005：332-613 [6]Louise.Schultz.Digital Processing: A System Orientation. Prentice Hall,1963:101-153

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

实验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的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

基于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的数字滤波器设计(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很快成为应用学科等领域不可或缺的基础软件。它可以快速有效地实现数字滤波器的设计、分析和仿真，极大地减轻了工作量,有利于滤波器设计的最优化。

IIR数字滤波器设计matlab代码

大连理工大学实验报告 学院（系）：电子信息与电气工程学部专业：电子信息工程班级： 姓名：学号：组：实验时间：实验室：大黑楼C221实验台：指导教师签字：成绩：实验名称 实验三IIR 数字滤波器设计 一、实验目的和要求 (1)掌握双线性变换法设计IIR 数字滤波器。(2)掌握脉冲响应不变法设计IIR 数字滤波器。 二、实验原理和内容 1.巴特沃思低通滤波器 巴特沃思低通滤波器幅度平方函数定义为 |H a (jΩ)|2=N c 211 ???? ??ΩΩ+式中N 为正整数，代表滤波器的阶次，c Ω称为截止频率。当c Ω=Ω时，有21 )(2 = Ωc a j H 即dB j H j H c a a 3) ()0(lg 201=Ω=δ所以又称c Ω为巴特沃思低通滤波器的3分贝带宽。 巴特沃思低通滤波器的幅度特性如图所示。

图1巴特沃思低通滤波器的幅度特性 2.切贝雪夫低通滤波器 切贝雪夫I 型滤波器的幅度平方函数为 ??? ? ??ΩΩ+=Ωc N a C j H 22211 )(ε其中，ε是表示通带波纹大小的一个参数，ε越大，波纹也越大。 c ΩΩ为Ω对c Ω的归一化频率，c Ω为截止频率，)(x C N 是N 阶切贝雪夫多项式。 切贝雪夫I 型低通滤波器的幅度特性如图所示 图2切贝雪夫I 型低通滤波器的幅度特性 三、主要仪器设备 Matlab 四、实验内容及结果 1.用双线性变换法设计一个巴特沃斯IIR 低通数字滤波器。设计指标参数为：在通带内频率低于0.4*pi 时，最大衰减小于2dB ；在阻带内[0.6*pi,pi]频率区间上，最小衰减大于15dB 。要求： ?(1)用buttord ，buttap ，zp2tf ，lp2lp ，bilinear 等指令实现 ?(2)用buttord ，butter 指令实现

实验4 基于MATLAB的FIR数字滤波器设计

实验4 基于MATLAB 的FIR 数字滤波器设计 实验目的：加深对数字滤波器的常用指标和设计过程的理解。 实验原理：低通滤波器的常用指标： P P P for H Ω≤Ω+≤Ω≤-,1)(1δδ πδ≤Ω≤Ω≤ΩS S for H ,)( 通带边缘频率P Ω，阻带边缘频率S Ω ，通带起伏 P δ， 通带峰值起伏] )[1(log 2010dB p p δα--=， 阻带起伏s δ，最小阻带衰减])[(log 2010dB s S δα-=。 数字滤波器有IIR 和FIR 两种类型，它们的特点和设计方法不同。 在MATLAB 中，可以用b=fir1(N,Wn,’ftype’,taper) 等函数辅助设计FIR 数字滤波器。N 代表滤波器阶数；Wn 代表滤波器的截止频率(归一化频率)，当设计带通和带阻滤波器时，Wn 为双元素相量；ftype 代表滤波器类型，如’high ’高通，’stop ’带阻等；taper 为窗函数类型，默认为海明窗，窗系数需要实现用窗函数blackman, hamming,hanning chebwin, kaiser 产生。 S P P S Passband Stopband Transition band Fig 1 Typical magnitude specification for a digital LPF

例1 用凯塞窗设计一FIR低通滤波器，通带边界频率π3.0 ，阻带边界频 Ω = p ，阻带衰减不小于 率π5.0 Ω = s 50dB。

解首先由过渡带宽和阻带衰减来决定凯塞窗的N和 π2.0 = Ω - Ω = ?Ω p s ， ，

用MATLAB设计FIR数字滤波器

实验八 用MATLAB 设计FIR 数字滤波器（二） 一、实验目的： 1、加深对窗函数法设计FIR 数字滤波器的基本原理的理解。 2、学习用MATLAB 语言的窗函数法编写设计FIR 数字滤波器的程序。 3、了解MATLAB 语言有关窗函数法设计FIR 数字滤波器的常用函数用法。 二、实验原理： 1、用窗函数法设计FIR 数字滤波器 FIR 数字滤波器的系统函数为 N-1 -n n=0H(z)=h(n)z ∑ 这个公式也可以看成是离散LSI 系统的系统函数 M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 分母a 0为1，其余a k 全都为0时的一个特例。由于极点全部集中在零点，稳定和线性相位特性是FIR 滤波器的突出优点，因此在实际中广泛使用。 FIR 滤波器的设计任务是选择有限长度的h(n)，使传输函数H(e j ω)满足技术要求。主要设计方法有窗函数法、频率采样法和切比雪夫等波纹逼近法等。本实验主要介绍窗函数法。 用窗函数法设计FIR 数字滤波器的基本步骤如下： （1）根据过渡带和阻带衰减指标选择窗函数的类型，估算滤波器的阶数N 。 （2）由数字滤波器的理想频率响应H(e j ω)求出其单位脉冲响应h d (n)。

可用自定义函数ideal_lp实现理想数字低通滤波器单位脉冲响应的求解。 程序清单如下： function hd=ideal_lp(wc,N) %点0到N-1之间的理想脉冲响应 %wc=截止频率（弧度） %N=理想滤波器的长度 tao=(N-1)/2; n=[0:(N-1)]; m=n-tao+eps; %加一个小数以避免0作除数 hd=sin(wc*m)./(pi*m); 其它选频滤波器可以由低通频响特性合成。如一个通带在ωc1～ωc2之间的带通滤波器在给定N值的条件下，可以用下列程序实现： Hd=ideal_lp(wc2,N)-ideal_lp(wc1,N) （3）计算数字滤波器的单位冲激响应h(n)=w(n)h d(n)。 （4）检查设计的滤波器是否满足技术指标。 如果设计的滤波器不满足技术指标，则需要重新选择或调整窗函数的类型，估算滤波器的阶数N。再重复前面的四个步骤，直到满足指标。 常用的窗函数有矩形窗、三角形窗、汉宁窗、哈明窗、切比雪夫窗、布莱克曼窗、凯塞窗等，MATLAB均有相应的函数可以调用。另外，MATLAB信号处理工具箱还提供了fir1函数，可以用于窗函数法设计FIR滤波器。 由于第一类线性相位滤波器（类型Ⅰ）能进行低通、高通、带通、带阻滤波器的设计，因此，本实验所有滤波器均采用第一类线性相位滤波器。 三、实验容：

实验报告基于MATLAB的数字滤波器设计

实验7\8 基于MATLAB 的数字滤波器设计 实验目的：加深对数字滤波器的常用指标和设计过程的理解。 实验原理：低通滤波器的常用指标： P P j P for e G ωωδδω ≤+≤≤-,1)(1 π ωωδω ≤≤≤S S j for e G ,)( 通带边缘频率：P ω，阻带边缘频率：S ω ， 通带起伏：P δ，通带峰值起伏： ])[1(log 2010 dB p p δα --=，阻带起伏： s δ， 最小阻带衰减： ] )[(log 2010dB s S δα-=。 数字滤波器有IIR 和FIR 两种类型，它们的特点和设计方法不同。 在MATLAB 中，可以用[b ，a]=butter （N,Wn ）等函数辅助设计IIR 数字滤波器, 也可以用b=fir1(N,Wn,’ftype ’) 等函数辅助设计FIR 数字滤波器。 实验内容： 利用MATLAB 编程设计一个数字带通滤波器，指标要求如下： 通带边缘频率：πω45.01=P ，πω65.02=P ，通带峰值起伏：] [1dB p ≤α 。 阻带边缘频率： π ω3.01=S ， π ω75.02=S ，最小阻带衰减： ] [40dB S ≥α。 分别用IIR 和FIR 两种数字滤波器类型进行设计。 实验要求：给出IIR 数字滤波器参数和FIR 数字滤波器的冲激响应，绘出它们的幅度和相位 频响曲线，讨论它们各自的实现形式和特点。 实验内容： IRR 代码： wp=[0.45*pi,0.65*pi]; ws=[0.3*pi,0.75*pi]; Ap=1; As=40; [N,Wc]=buttord(wp/pi,ws/pi,Ap,As) ; [b,a]=butter(N,Wc)%[b,a] = butter(n,Wn,'ftype') S ω-P ω-P ωS ωPassband Stopband Transition band Fig 7.1 Typical magnitude specification for a digital LPF

Matlab数字滤波器设计

1 基于Matlab的FIR数字滤波器设计1郑伟鑫 ( 1. 华南理工大学机械与汽车工程学院机械电子工程1班 ) 摘要：数字滤波器不是用物理元器件构成，而是按照某种算法编写一段程序或模块，对数字信号进行处理加工，从而达到滤波的要求。数字滤波器是实现滤波过程的一种数字信号处理系统，具有离散时间系统的基本特征。鉴于滤波器在数字信号处理中的广泛应用以及在Matlab环境下设计滤波器的快捷方便等特性，本文研究了Matlab环境下，数字滤波器的设计与实现。数字滤波器包括IIR（无限脉冲响应）数字滤波器和FIR（有限脉冲响应）数字滤波器。其实现方式主要有软件方式实现（即Matlab编程实现）和硬件方式实现。本文主要介绍FIR数字滤波器的软硬件方式实现。Matlab强大的数字信号处理工具箱为我们提供了易用、有力的设计工具。通过正确的使用函数，很容易设计出性能良好的数字滤波器。关键词：数字滤波器 Matlab FIR Design of FIR digital filter based on Matlab 1 Zheng Weixin(1.South China University of Technology, Mechanical and Automotive institute, Class 1 of Mechanical and Electronic Engineering)：Abstract Digital filter is not composed of physical components, but writing a program or module according to some algorithms. The digital signal processing is aimed to achieve the filtering requirements. The digital filter is a digital signal processing system implemented in filtering process, which is also the basic feature of the discrete time system. In view of the extensive applications of filters in digital signal processing and the fast and convenient features of filter