FIR数字滤波器设计

FIR数字滤波器设计

摘要 数字滤波器的输入输出均为数字信号，信号通过数字滤波器后，可以改变

频率成分的相对比例或滤除某些频率成分。数字滤波器可以分为IIR数字滤波器

和FIR数字滤波器。文中只介绍FIR数字滤波器的设计，可以根据所给定的频率

特性直接设计FIR数字滤波器。FIR数字滤波器在保证幅度特性满足要求的同时，

能够做到严格的线性特性。本文采用了窗函数法、频率采样法以及基于firls函

数和remez函数的最优化方法设计FIR滤波器。在文中对FIR滤波器进行了详细

的理论分析，并且对应于每种方法都给出了设计实例。通过编写MATLAB语言程序，

运行程序，得到幅频和相频特性图。对于窗函数和firls函数设计的滤波器，还

通过建立Simulink系统模块进行仿真，观察滤波器滤波情况。

关键词 窗函数法；频率采样法；最优化方法；仿真

I

Design of FIR Digital Filter

Abstract Firstly, the paper gives a brief description of digital filter and MATLAB

language.The digital filter can be divided into IIR digital filter and FIR digital filter.The paper only give the digital filter of FIR, 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，frequency sampling method and optimization method.As for all filter design method above, the paper all describes in detail. And for each method, the paper also gives some examples. At last, through the preparation of the relevant procedures in MATLAB simulation software, operating procedures and records, observation of frequency and phase of plans to FIR digital filter in the MATLAB software in the implementation process.

Keywords

II

目 次

引言 ............................................... 1 1数字滤波器 ........................................ 2 1.1数字滤波器的概述 ............................. 2 1.2数字滤波器的分类 ............................. 2 1.3数字滤波器设计指标 ........................... 3 2 FIR数字滤波器设计基础 ............................ 4 2.1 FIR数字滤波器的特点 ......................... 4 2.2 FIR数字滤波器的线性相位条件： ............... 4 2.3 FIR数字滤波器的基本结构 ..................... 5 3 FIR数字滤波器设计 ................................ 9 3.1 窗函数法设计FIR数字滤波器 .................. 9 3.2频率采样法设计FIR数字滤波器 ................ 15 3.3最优化方法设计FIR数字滤波器 ................ 17 结 论 ............................................. 21 参考文献 .......................................... 23 致 谢 ............................ 错误！未定义书签。

III

引言

随着信息时代与数字技术的发展，数字信号处理己逐渐发展成为当今极其重要的学科与技术领域之一。数字信号处理在通信、语音、图像、自动控制雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理的基本方法中，通常会涉及到变换、滤波、频谱分析、调制解调和编码解码等处理。其中，滤波是应用非常广泛的一个环节，数字滤波器的相关理论也一直都是人们研究的重点之一。数字滤波器是数字信号处理的重要基础，在对信号的滤波、检测及参数的估计等信号应用中，数字滤波器是使用最为广泛的一种线性系统。数字滤波器根据其单位冲击响应函数的时域特性可分为两类：无限冲击响应(IIR)数字滤波器和有限冲击响应(FIR)数字滤波器。与IIR数字滤波器相比,FIR数字滤波器的实现是非递归的，稳定性好，精度高；更重要的是FIR数字滤波器在满足幅度响应要求的同时，可以获得严格的线性相位。因此，它在高保真的信号处理中，如数字音频、图像处理、数据传输和生物医学等领域得到广泛应用。

1

1数字滤波器

1.1数字滤波器的概述

所谓数字滤波器，是指输入输出均为数字信号，通过一定的运算关系，改变输入信号中所含频率成分的相对比例，或则滤除某些频率成分的器件[3]。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。对于数字滤波器而言，若系统函数为H(z),其脉冲响应为h(n),输入时间序列为x(n),则它们在时域内的关系式如下：

y(n)?h(n)?x(n) (1-1) 在Z域内，输入和输出存在如下关系：

Y(z)?H(z)X(z) (1-2)

式中， X(z)、Y(z)分别为x(n)和y(n)的Z变换。 在频域内，输入和输出则存在如下关系：

Y(j?)?H(j?)X(j?) (1-3) 式中，H(j?)是数字滤波器的频率特性；X(j?)、Y(j?)分别为x(n)和y(n)的频谱，而?为数字角频率。

1.2数字滤波器的分类

数字滤波器可以有很多种分类方法，但总体上可分为两大类。一类称为经典滤波器，即一般的滤波器，其特点是输入信号中的有用成分和希望滤除的成分占用不同的频带，通过合适的选频滤波器可以实现滤波[4]。例如，若输入信号中有干扰，信号和干扰的频带互不重叠，则可滤出信号中的干扰得到纯信号。但是，如果输入信号中信号和干扰的频带相互重叠，则干扰就不能被有效的滤除。另一类称为现代滤波器，如维纳滤波器、卡尔曼滤波器等，其输入信号中有用信号和希望滤除的频带成分重叠。对于经典滤波器，从频域上也可以分为低通、高通、带通和带阻滤波器。从时域特性上看，数字滤波器还可以分为有限脉冲响应（FIR，finite impulse response）数字滤波器和无限脉冲响应（IIR, infinite impulse response）数字滤波器[5]。

2

对于有限脉冲响应（FIR）数字滤波器，其输出y(n)只取决于有限个过去和现在的输入，x(n),x(n-1),…,x(n-m),滤波器的输入输出关系可表示为

My(n)??brx(n?r) (1-4)

r?0对于无限脉冲响应（IIR）数字滤波器，它的输出不仅取决于过去和现在的输入，而且还取决于过去的输出，其差分方程为

NMy(n)??aky(n?k)??brx(n?r) (1-5)

k?1r?0该差分方程的单位冲激响应是无限延续的。

1.3数字滤波器设计指标

设数字滤波器的传输函数用下式表示：

H(ej?)?H(ej?)ej?(?)

(1-6)

式中，｜H(ej?)｜为幅频特性，?(?)为相频特性[6]。幅频特性表示信号通过滤波器后各频率成分的衰减情况，相频特性则反映各频率成分通过滤波器后在时间上的延时情况。通常，选频滤波器的指标要求都以幅频特性给出，对相频特性不作要求，如果需要对输出波形有严格要求，如语音合成、波形传输等，则要求设计线性相位数字滤波器[7]。

数字滤波器的参数指标是?p、?s、?p和?s。?p和?s分别称为通带截止频率和阻带截止频率。通带和阻带内允许的衰减一般用分贝数表示，通带内允许的最大衰减用?p表示，阻带内允许的最小衰减用?s表示，?p和?s分别定义为

?p?20lgH(ej0)H(ep)j???20lgH(ep) dB (1-7）

j??s?20lgH(ej0)H(ej?s)??20lgH(ej?s) dB (1-8)

式中均假定H(ej0)已被归一化为1

3

2 FIR数字滤波器设计基础

2.1 FIR数字滤波器的特点

FIR滤波器在保证幅度特性的同时，很容易做到严格的线性相位特性。在数字滤波器中，FIR滤波器的最主要特点是没有反馈回路，故不存在不稳定的问题；同时，在幅度特性可以任意设置的同时，保证了精确的线性相位。稳定和线性相位是FIR滤波器的突出优点。另外还有以下特点：设计方式是线性的；硬件容易实现；滤波器过渡过程具有有限区间；相对IIR滤波器而言，阶次较高，其延迟也要比同样性能的IIR滤波器大得多[8]。

2.2 FIR数字滤波器的线性相位条件：

设滤波器单位脉冲响应的长度为N，系统函数为

H(z)??h(n)z?n (2-1)

n?0N?1由此式可见，H(z)是z?1的(N-1)次多项式,它在Z平面上有(N-1)个零点，原点z=0是(N-1)阶重极点，位于r =1的单位圆内，系统永远稳定。稳定性和线性相位特性是FIR滤波器的突出优点。

FIR滤波器的设计任务是选择有线长度的h(n)，使传输函数H(ej?)满足要求。

线性相位条件：

对于长度为N的h(n)，传输函数为

N?1n?0H(e)??h(n)e?j?n (2-2)

jwH(ejw)?Hg(?)e?j?(w) (2-3)

式中,Hg(?)称为幅度特性，?(?)称为相位特性。线性相位是指相位函数?(?)满

4

足如下特性：?(w)???w或?(?)??0??w, ?0是起始相位，?为常数，一般称第一种情况为第一类线性相位，称第二种情况为第二类线性相位。

满足第一类线性相位的充要条件是：h(n)为实序列，并且对(N-1)/2偶对称，即h(n)?h(N?n?1)；满足第二类线性相位的充要条件是：h(n)为实序列，并且对(N-1)/2奇对称，即h(n)??h(N?n?1)。

2.3 FIR数字滤波器的基本结构

FIR滤波器的基本结构有以下几种：直接型、级联型、线性相位型、频率采样型。

1. 直接型

设FIR滤波器的单位冲击响应h(n)为一个长度为N的序列，则滤波器系统函数为：

N?1n?0H(z)??h(n)z?n （2-4）

表示这一系统输入输出关系的差分方程为

N?1y(n)??h(m)x(n?m) (2-5)

m?0直接由差分方程可得出对应的网络结构如图2-1所示：

x(n)z?1h(0)z?1h(1)h(2)z?1h(N-2)h(N-1)y(n)

图2-1 FIR滤波器的直接型结构 直接型结构的优点：简单直观，乘法运算量较少。 缺点：调整零点较难。 2．级联型

当需要控制滤波器的传输零点时，可将H(z)分解为实系数二阶因子的乘积形式：

N/2k?1H(z)??(?0k??1kz?1??2kz?2) (2-6)

5

式中，H(z)为h(n)的z变换，?0k，?1k，?2k为实数。级联型结构如图2-2所示：

x(n)?01z?1z?1?02z?1z?1?N??0???2?y(n)?11?12z?1?1?N???2??Nz?1?2????2??

?21?22图2-2 FIR滤波器的级联型结构

该结构的优点：调整零点比直接型方便。

缺点：H(z)中的系数比直接型多，因而需要的乘法器多。当H(z)的阶次高时，也不易分解。

3.线性相位型结构

FIR滤波器的线性相位结构有偶对称和奇对称，不论h(n)为偶对称还是奇对称都有：

当N为偶数时，系统函数为

N?12n?0H(z)??h(n)?z?n?z?(N?1?n)? (2-7)

当N为奇数时，系统函数为

H(z)??h(n)?zn?0N?12?n?z?(N?1?n)??h???N2?N?1????1??2?? (2-8) ?z?对这两种情况，都可以用FIR直接型实现，其信号流图如图2-3所示。

x(n)z?1?z?1z?1y(n)b?0?z?1h?1?z?1?H?h??1??2?

(a)N为偶数

6

x(n)z??z?1z?1z?1y(n)b?0?z?1h?1??H?h??1??2?

(b)N为奇数 图2-3 线性相位型结构

这种结构在本质上是直接型，但乘法次数比直接型省了一半。

4．频率采样型

频率采样型结构是一种用系数将滤波器参数化的一种实现结构。一个有限长序列可以由相同长度频域采样值惟一确定。

系统函数在单位圆上作N等分取样就是单位取样相应h(n)的离散傅里叶变换

H(k)。H(k)与系统函数之间的关系可用内插公式表示：

H(k) (2-9) H(z)?(1)(1?z?N)?k?1N1?WZN式中 Hc(z)?1?z?N

Hk(z)?H(k) ?k?11?WNZ这样，H(z)是由梳状滤波器Hc(z)和N个一阶网络Hk(z)的并联结构进行级联而成的，其网络结构（信号流图）如图2-3所示。Hc(z)是一个梳妆网络，其零点为

?kWN?exp(jk2?), k= 0, 1，2…，N-1 N刚好和极点一样，等间隔地分布在单位圆上。理论上，极点和零点相互抵消，保证了网络的稳定性。

7

x(n)H(0)y(n)1N?z?N0WNz?1H(1)?1WNz?1H(N-1)?N?1WNz?1

图2-5 FIR滤波器的频率采样结构

频率采样结构的优点：

1)在频率采样点?k，H(ej?k)?H(k)，只要调整H(k)就可以有效地调整频响特性。

2）只要h(n)长度N相同，对于任何频响，其梳状滤波器部分和N个一阶网络部分完全相同，只是各支路增益H(k)不同。相同部分便于标准化、模块化。

缺点：

1）寄存器长度都是有限的，零、级点可能不能正好抵消，造成系统不稳定。2）当N很大时，其结构很复杂，需要的乘法器和延时单元很多。

8

3 FIR数字滤波器设计

FIR数字滤波器的设计方法有窗函数法、频率采样法和基于firls函数和

remez函数的最优化方法。MATLAB语言中的数字信号处理工具箱，提供了一些滤波器的函数，使FIR滤波器的运算更加方便和快捷。在MATLAB中提供的滤波函数有fir1(),此函数以经典的方法实现加窗线性相位FIR数字滤波器设计，可以设计出低通、高通、带通和带阻滤波器；fir2函数设计的FIR滤波器，其滤波的频率特性由矢量f和m决定，f和m分别为滤波器的期望幅频响应的频率相量和幅值相量。Firls()和remez()的基本格式用于设计I型和II型线性相位FIR滤波器，I型和II型的区别是偶函数还是奇函数。freqz()用于求数字滤波器的频率响应。并且提供了各种窗函数的函数，比如，hamming()是海明窗函数，hanning()是汉宁窗函数，kaiser()是凯泽窗函数，使在设计的过程中，不用自己重新设计窗函数。

Simulink是MATLAB众多工具包中的一员，对于建模，Simulink提供了一个图形化的用户界面（GUI）。Simulink包括一个复杂的由接收器、信号源、线性和非线性组件以及连接件组成的模块库。定义完一个模型后，就可以通过Simulink的菜单或者在MATLAB的命令窗口输入命令对它进行仿真。使用Scopes或者其他的显示模块，可以在运行仿真时观察到仿真的结果。另外，还可以在仿真时改变参数，并且立即就可看到变化。

3.1 窗函数法设计FIR数字滤波器

设我们所要设计的FIR滤波器的传输函数是Hd(ej?),hd(n)是与其对应的单位脉冲响应，因此

?Hd(ej?)?n????j?nh(n)e （3-1） ?d1hd(n)?2?????Hd(ej?)ej?nd? （3-2)

如果我们能够在Hd(ej?)已知的情况下，求出hd(n)，经过Z变换可得到滤波器的系统函数。通常情况下理想数字滤波器的单位脉冲相应hd(n)是无限长的，且是非

9

因果序列。获得有限脉冲响应滤波器的一种可能方法是对hd(n)截取一段h(n)来近似代替hd(n)，可是这样会改变原来的滤波器指标，出现吉布斯效应误差。

窗函数法就是用被称为窗函数的有限加权序列w(n)来修正式(1)的傅里叶基数以求得要求的有限脉冲响应序列h(n)，即

h(n)?hd(n)w(n) （3-3）

w(n)是有限长序列，当n<0或n>N-1时，w(n)=0。

这种方法的重点在于选择某种合适的窗函数。要求窗函数主瓣宽度尽可能窄，以获得最小的过渡带；旁瓣相对值尽可能小，以使得通带波纹小，并且阻带衰减大。

下面介绍几种常用的窗函数：

1.矩形窗(Rectangle Window)

wR(n)?RN(n)（3-4）

其频率函数为：

N?1?j?sin(N?/2)2WR(ej?)?e （3-5）

sin(?/2)

2.三角形窗(Bartlett Window)

1?2n,0?n?(N?1)??N?12wBr(n)?? （3-6）

2n1?2?,(N?1)?n?N?1?N?12?其频率函数为：

N?1?j(???)??2sin(N?/4)j?2WBr(e)??e （3-7) ?N?sin(?/2)?23.汉宁(Hanning)窗，又称升余弦窗

1?2n??wHn(n)??1?cos()?RN(n) （3-8）

2?N?1?利用傅里叶变换得到频率函数为：

?2?2????j?N2?1?W(e)??0.5WR(?)?0.25?WR(??)?WR(??)??eN?1N?1????j??WHn(?)eN?1?j?2 （3-9）

当N?1时，N?1?N，所以窗函数的幅度函数为

2?2???WHn(?)?0.5WR(?)?0.25?WR(??)?WR(??)? (3-10)

NN??4.汉明(Hamming)窗，又称改进的升余弦窗

10

2n???wHm(n)??0.54?0.46cos()?RN(n) （3-11）

N?1??其幅度函数为：

2?2???WHm(?)?0.54WR(?)?0.23?WR(??)?WR(??)? （3-12）

N?1N?1??5.布莱克曼(Blankman)窗，又称二阶升余弦窗

2n?4n???wBl(n)??0.42?0.5cos()?0.08cos()?RN(n) （3-13)

N?1N?1??其幅度函数为：

2?2???WBl(?)?0.42WR(?)?0.25?WR(??)?WR(??)?N?1N?1?? (3-14）

4?4????0.04?WR(??)?WR(??)?N?1N?1??6.凯泽(Kaiser)窗

I0(?1?[1?2n/(N?1)]2)wk(n)?,I0(?)上述窗函数的基本参数如下表

窗函数 矩形窗 三角形窗 汉宁窗 汉明窗 布莱克曼窗 凯泽窗 旁瓣峰值幅度/db -13 -26 -31 -40 -57 -57 0?n?N?1 （3-15）

其中：β是一个可自由选择的参数，I0( x)是第一类修正零阶贝塞尔函数[10].

过渡带宽 4?/N 8?/N 8?/N 8?/N 12?/N 10?/N 阻带最小衰减/db -21 -25 -44 -53 -74 -80 窗函数法设计滤波器的步骤：

1）根据技术要求确定待求滤波器的单位取样响应hd(n)。

2）根据对过渡带和阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。 3）计算滤波器的单位取样响应h(n)：

h(n)?hd(n)?(n) （3-16)

式中,?(n)是前面所选择好的窗函数。

4）检验技术指标是否满足要求。根据下式计算：

11

H(e)??h(n)e?j?nj?n?0N?1 (3-17)

如果H(ej?)不满足要求，根据具体情况重复步骤(2)(3)(4)步，直到满足要求为止。

本文以一个FIR滤波器的设计为例说明如何使用MATLAB设计数字滤波器 设计实例：用窗函数法设计线性相位FIR低通数字滤波器，要求通带截止频率Wp=0.4*?,阻带截止频率Ws=0.5*?, 通带衰减不大于3db，阻带衰减不小于40db。 1.程序如下： Wp=0.4*pi; Ws=0.5*pi; Wdel=Ws-Wp; N=ceil(8*pi/Wdel); Wn=(0.4+0.5)*pi/2; window=hanning(N+1); b=fir1(N,Wn/pi,window); freqz(b,1,512)

程序执行后得幅频和相频如下图所示：

50Magnitude (dB)0-50-100-15000.10.20.30.40.50.60.70.8Normalized Frequency (?? rad/sample)0.910Phase (degrees)-1000-2000-3000-400000.10.20.30.40.50.60.70.8Normalized Frequency (?? rad/sample)0.91

12

2．用Simulink仿真对窗函数设计的低通滤波器进行滤波：

FIR低通滤波器的Simulink仿真模块图如下所示：

输入信号：s(t)=sin(30πt)+sin(56πt)+sin(70πt)

数字滤波器参数设置：已知滤波器的阶数n=80。本例中，首先在Filter Type中选择lowpass；在Design Method选项中选择FIR Window，接着在Window选项中选取hann；指定Filter Order项中的Specify order为80；采样频率Fs=100Hz,截止频率Fc=22.5HZ。设置完以后点击窗口下方的Design Filter，在窗口上方就会看到所设计滤波器的幅频响应和相频响应曲线。设计界面如下图所示：

13

信号源参数设置：三个信号源的频率分别是15Hz、28Hz、35Hz ：

运行仿真，得波形如下图所示：

14

15HZ信号波形 28HZ信号波形 35HZ信号波形

滤波前的叠加波形 滤波后波形

由上图可知，15Hz、28Hz、35Hz叠加后通过截止频率为22.5HZ的FIR低通滤波器，输出与15Hz信号波形相似，其它的高频被阻隔。

3.2频率采样法设计FIR数字滤波器

1.对理想滤波器的系统函数Hd(z)进行频率采样以得到系统的理想频响Hd(ejw)的等间隔采样值H(k)。H(k)实际上是所要求的滤波器的单位采样响应(h(n))的离散傅里叶变换（DFT），如下试：

2?1?jwN2?1N?1H(e)??H(k)?(w?k)?eH(k)e?NNk?0k?0jwN?12??jNNw)2w?k (3-18) sin(?)2Nsin(1h(n)?N?H(k)ek?0N?1j2?knN (3-19)

15

为了减小H（k）的通带边缘由于抽样点的变化而引起的起伏振荡，可以增加过渡点，加宽过渡带以减小通带的起伏。

每一个抽样值产生一个与sin(

Nww)/sin()成正比，并位移（2?k）/N的频率22响应，而H（k）与内插函数的线性组合就是FIR滤波器的频率响应，增加一点过渡可以使阻带衰减提高到-44~54dB,二点过渡衰减-65~75,三点过渡衰减-85~95dB.

如果不能使过渡带太宽，同时要求增大阻带衰减，可以增加取样点数N，但这样会增加计算量、延时和误差。

频率取样型FIR滤波器设计步骤： (1)给定理想滤波器频率响应Hd(ejw)。

(2)根据过渡带宽和阻带衰减确定过渡点数和h(n)的长度N。

H(k)?Hd(z)|z?e2?jk/N (3-20)

(3)由IFFT计算IDFT得到：

1H(n)?N?H(k)ek?0N?1j2?knN (3-21)

设计实例：率采样法设计一个带通滤波器，满足：低阻带边缘：w1s=0.2*?;低通带边缘：w1p=0.35*?;高通带边缘：w2p=0.65*?;高阻带边缘：w2s=0.8*?。设计过渡带中的频率样本值为t1和t2，取t1=0.109021,t2=0.59417456。 设计程序如下： M=40; al=(M-1)/2; l=0:M-1; t1=0.109021; t2=0.59417456;

Hrs=[zeros(1,5),t1,t2,ones(1,7),t2,t1,zeros(1,9),t1,t2,ones(1,7),t2,t1,zeros(1,4)]; k1=0:floor((M-1)/2); k2=floor((M-1)/2)+1:M-1;

angh=[-al*(2*pi)/M*k1,al*(2*pi)/M*(M-k2)]; H=Hrs.*exp(j*angh); h=real(ifft(H,M)); freqz(h,1,512,1000)

16

实验得幅频相频特性如下图所示：

200Magnitude (dB)0-200-400050100150200250300Frequency (Hz)3504004505001000Phase (degrees)0-1000-2000050100150200250300Frequency (Hz)350400450500

3.3最优化方法设计FIR数字滤波器

MATLAB7.0信号处理工具箱提供了通用的工具函数firls和remez，采用不同的优化方法可以设计最优滤波器和多频带滤波器。

firls是fir1和fir2函数的扩展，其基本设计准则是利用最小二乘法使期望的频率响应和实际的频率响应之间的整体误差最小。Remez函数则利用Parks-McClellan算法设计线性相位FIR滤波器。该算法能使期望频率响应和实际频率响应之间的最大误差最小。Firls函数设计的滤波器在频率响应上出现等波纹，因此也叫等波纹滤波器。函数firls、remez调用的语法规则相同。

(1)firls。其调用格式为 Ⅰ. b=firls(n,f,a) Ⅱ. b=firls(n,f,a,w) Ⅲ. b=firls(n,f,a,’ftype’) Ⅳ. b=firls(n,f,a,w,’ftype’)

该函数用于设计n阶FIR滤波器，其幅频特性由f和a向量确定，f是频率点向量，其范围为0~1，是递增向量；a是指定频率点的幅度响应，与f必须同长度；

17

b为函数返回的滤波器系数，长度为n+1,且具有偶对称关系：

b(k)=-b(n+2-k), k=1,2,…,n+1 b=firls(n,f,a,w)则使用权系数w给误差加权。

b=firls(n,f,a,w,’ftype’)调用形式中，参数’ftype’用于指定所设计的滤波器类型，ftype=Hilbert,为奇对称的线性相位滤波器，返回的滤波器系数满足b(k)=-b(n+2-k), k=1,2,…,n+1；ftype=differentiator,则采用特殊加权技术，生成奇对称的线性相位滤波器，使低频段误差大大小于高频段误差。

(2)函数remez的基本调用格式为 b=remez(n,f,a) 其中各项意义同函数firls所述。

设计实例：分别使用函数firls和remez函数设计20阶带通滤波器，阻带频率为ws1=0.3*?,ws2=0.7*?,通带频率为wp1=0.4*?,wp2=0.6*?,并画出幅频特性曲线。

1.程序如下所示： f=[0 0.3 0.4 0.6 0.7 1]; m=[0 0 1 1 0 0]; n=20;

b=firls(n,f,m); %firls函数设计FIR数字滤波器 [h,w]=freqz(b);

bb=remez(n,f,m); %remez函数设计FIR数字滤波器 [hh,w]=freqz(bb);

axes('position',[0.2 0.2 0.5 0.5]); plot(w/pi,abs(h),'b:',w/pi,abs(hh),'r-'); xlabel('frequency'); ylabel('magnitude'); grid on;

实验的幅频特性曲线如下：

18

1.41.21magnitude0.80.60.40.2000.10.20.30.40.50.6frequency0.70.80.91

2.用Simulink仿真对用firls函数设计的带通滤波器进行滤波： 仿真模块如下：

输入信号：s(t)=sin(20πt)+sin(50πt)+sin(90πt)

数字滤波器参数设置：已知滤波器的阶数n=20。本例中，首先在Filter Type中选择bandpass；在Design Method选项中选择FIR least-squares;指定Filter Order项中的Specify order为20；采样频率Fs=100Hz,Fstop1=15HZ，Fpass1=20HZ,

19

Fpass2=30HZ,Fstop2=35HZ。设置完以后点击窗口下方的Design Filter，在窗口上方就会看到所设计滤波器的幅频响应曲线。。设计界面如下图所示：

信号源参数设置：三个信号源的频率分别是10Hz、25Hz、45Hz ：

运行仿真，得波形如下图所示：

20

（1）10HZ信号波形 （2） 25HZ信号波形 （3）45HZ信号波形

(4） 滤波前波形 (5) 滤波后波形

由以上波形可知，10Hz、25Hz、45Hz叠加后通过通频带为20Hz~30Hz的FIR带通滤波器，输出与25Hz信号相似，低频和高频部分被阻隔。

21

结 论

本文在数字滤波器的设计过程中，采用的设计方法是基于MATLAB的数字滤波器的设计。设计过程中充分发挥了MATLAB语言编程效率高的特点以及Simulink仿真能力强的特点，完成各种计算与图形绘制工作，使得数字滤波器的设计变得简便快捷。本文通过FIR数字滤波器的设计实例，说明如何利用MATLAB来完成数字滤波器的设计，设计结果的各项性能指标均达到指定要求，Simulink仿真验证了设计方法的正确性。本文所述设计过程简便易行，具有很好的实际运用意义。

21

参考文献

【1】陈怀琛，吴大正, 高西全.MATLAB及在电子信息课程中的应用.北京：电子工业出版社.2009年2月第三版

【2】张亮，郭仕剑等.MATLAB7.X系统建模与仿真.人民邮电出版社.2006.11 【3】程佩青.数字信号处理教程.北京:清华大学出版社，2003

【4】丁吉，姜涛.基于Matlab 的FIR数字滤波器的设计.长春工业大学学报(自然科学版).2006，9

【5】罗敏，刘嵩，刘忠兵.基于MATLAB的IIR数字滤波器的设计.山西电子术.2007年第3期

【6】姜献忠.MATLAB在数字滤波器设计中的应用.南京工业职业技术学院学报.2003年12月第3卷第4期

【7】丁玉美，高西全.数字信号处理.西安：西安电子科技大学出版社，2005 【8】廖颀，刘良福，肖士斌.基于MATLAB的数字滤波器设计实现[J].赣南师范

学院学报.2006年第6期

【9】丁玉美，高西全.数字信号处理. 西安：西安电子科技大学出版社，2005第二版

【10】罗军辉，罗勇江，白义臣等.MATLAB7.0在数字信号处理中的应用. 北京：机械工业出版社，2005年5月第一版

23

参考文献

【1】陈怀琛，吴大正, 高西全.MATLAB及在电子信息课程中的应用.北京：电子工业出版社.2009年2月第三版

【2】张亮，郭仕剑等.MATLAB7.X系统建模与仿真.人民邮电出版社.2006.11 【3】程佩青.数字信号处理教程.北京:清华大学出版社，2003

【4】丁吉，姜涛.基于Matlab 的FIR数字滤波器的设计.长春工业大学学报(自然科学版).2006，9

【5】罗敏，刘嵩，刘忠兵.基于MATLAB的IIR数字滤波器的设计.山西电子术.2007年第3期

【6】姜献忠.MATLAB在数字滤波器设计中的应用.南京工业职业技术学院学报.2003年12月第3卷第4期

【7】丁玉美，高西全.数字信号处理.西安：西安电子科技大学出版社，2005 【8】廖颀，刘良福，肖士斌.基于MATLAB的数字滤波器设计实现[J].赣南师范

学院学报.2006年第6期

【9】丁玉美，高西全.数字信号处理. 西安：西安电子科技大学出版社，2005第二版

【10】罗军辉，罗勇江，白义臣等.MATLAB7.0在数字信号处理中的应用. 北京：机械工业出版社，2005年5月第一版

23

本文来源：https://www.bwwdw.com/article/u2q.html