窗函数法设计FIR数字低通滤波器课程设计文档格式.docx

窗函数法设计FIR数字低通滤波器课程设计文档格式.docx

《窗函数法设计FIR数字低通滤波器课程设计文档格式.docx》由会员分享，可在线阅读，更多相关《窗函数法设计FIR数字低通滤波器课程设计文档格式.docx（16页珍藏版）》请在冰豆网上搜索。

FIR滤波器在截止频率的边沿陡峭性能虽然不及IIR滤波器，但是，考虑到FIR滤波器严格的线性相位特性和不像IIR滤波器存在稳定性的问题，FIR滤波器能够在数字信号处理领域得到广泛的应用。

FIR是有限冲激响应（FiniteImpulseResponse）的简称。

由线性系统理论可知，在某种适度条件下，输入到线性系统的一个冲击完全可以表征系统。

当我们处理有限的离散数据时，线形系统的响应（包括对冲击的响应）也是有限的。

若线性系统仅是一个空间滤波器，则通过简单地观察它对冲击的响应，我们就可以完全确定该滤波器。

通过这种方式确定的滤波器称为有限冲击响应（FIR）滤波器。

FIR滤波器是在数字信号处理（DSP）中经常使用的两种基本的滤波器之一。

FIR滤波器具有严格的相位特性，对于信号处理和数据传输是很重要的。

目前FIR滤波器的设计方法主要有三种：

窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。

常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。

用窗函数设FIR滤波器的基本思路：

从时域出发设计h（n）逼近理想hd（n）。

设理想滤波器的单位响应在时域表达为hd（n）,则Hd（n）一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。

要想得到一个因果的有限长的滤波器单位抽样响应h（n）最直接的方法是先将hd（n）往右平移，再进行截断，即截取为有限长因果序列：

h（n）=hd（n）w（n），并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。

按照线性相位滤波器的要求，线性相位FIR数字低通滤波器的单位抽样响应h（n）必须是偶对称的。

对称中心必须等于滤波器的延时常数，即用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，这个现象称为吉布斯（Gibbs）效应。

为了消除吉布斯效应，一般采用其他类型的窗函数。

MATLAB设计FIR滤波器有多种方法和对应的函数。

窗函数设计法不仅在数字滤波器的设计中占有重要的地位，同时可以用于功率谱的估计，从根本上讲，使用窗函数的目的就是消除由无限序列的截短而引起的Gibbs现象所

带来的影响。

一、课程设计的目的和意义

1、通过课程设计把自己在大学中所学的知识应用到实践当中。

2、深入了解利用Matlab设计FIR数字滤波器的基本方法。

3、在课程设计的过程中掌握程序编译及软件设计的基本方法。

4、提高自己对于新知识的学习能力及进行实际操作的能力。

5、锻炼自己通过网络及各种资料解决实际问题的能力。

二、课程设计题目描述和要求

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

题目:

对模拟信号进行低通滤波处理，要求通带0≤f≤2kHZ内衰减小于1dB，阻带3kHZ≤f≤∞上衰减大于30dB。

希望对模拟信号采样后用线性相位FIR数字滤波器实现上述滤波，采样频率Fs=10kHZ。

要求:

用窗函数法设计满足要求的FIR数字低通滤波器，为了降低运算量，滤波器阶数尽量低。

三、设计原理

1.基本原理

设计低通FIR数字滤波器，寻求一系统函数

，使其频率响应

逼近滤波器要求的理想频率响应

，其对应的单位脉冲响应

（1）

其中

如果所希望的滤波器的理想的频率响应函数为

，则其对应的单位脉冲响应为

（2）

窗函数设计法的基本原理是用有限长单位脉冲响应序列

逼近

。

由于

往往是无限长序列，而且是非因果的，所以用窗函数

将

截断，并进行加权处理，得到：

（3）

就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数

为

（4）

式中，N为所选窗函数

的长度。

2.典型的窗函数

（1）矩形窗（RectangleWindow）

（5）

其频率响应和幅度响应分别为：

，

（6）

（2）三角形窗（BartlettWindow）

（7）

其频率响应为：

（3）汉宁（Hanning）窗，又称升余弦窗

（8）

（4）汉明（Hamming）窗，又称改进的升余弦窗

（9）

其幅度响应为：

（5）布莱克曼（Blankman）窗，又称二阶升余弦窗

（10）

其幅度响应：

（6）凯塞（Kaiser）窗

（11）

其中：

β是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的交换关系，一般说来，β越大，过渡带越宽，阻带越小衰减也越大。

I0（·

）是第一类修正零阶贝塞尔函数。

若阻带最小衰减表示为

，β的确定可采用下述经验公式：

若滤波器通带和阻带波纹相等即δp=δs时，滤波器节数可通过下式确定：

式中：

我们知道，用窗函数法设计的滤波器性能取决于窗函数

的类型及窗口长度N的取值。

设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。

各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见下表1

表1各种窗函数的基本参数

窗函数

旁瓣峰值幅度/dB

过渡带宽

阻带最小衰减/dB

矩形窗

-13

4π/N

-12

三角形窗

-25

8π/N

汉宁窗

-31

-44

哈明窗

-41

-53

布莱克曼窗

-57

12π/N

-74

凯塞窗（α=7.865）

10π/N

-80

这样选定窗函数类型和长度N之后，求出单位脉冲响应

，并按照式（4）求出

是否满足要求，要进行演算。

一般在

尾部加零使长度满足2的整数次幂，以便用FFT计算

如果要观察细节，补零点数增多即可。

如果

不满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。

如果要求线性相位特性，则

还必须满足

（12）

根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。

要根据所设计的滤波特性正确选择其中一类，例如，要设计线性相位低通特性，可以选择

这一类，而不能选择

这一类。

我们在设计滤波器时，希望窗谱主瓣尽可能窄，以获得较陡的过渡带，同时尽可能减小最大旁瓣的相对幅度，我们发现这是不可能的，因为这本身就是一个矛盾体，所以在设计滤波器时只是根据实际情况来选择合适的窗函数

主程序框图如图下图所示。

其中幅度特性要求用dB表示。

四、FIR滤波器的窗函数设计基本方法

1.FIR滤波器简介

FIR数字滤波器设计最简单的方法是窗函数法，通常也称为傅立叶级数法。

它是在时域进行的，因而必须由理想滤波器的频率响应

推导出其单位冲激响应

，在设计一个FIR数字滤波器的单位冲激响应

去逼近

根据冲激响应的时域特性，数字滤波器可分为无限长冲激响应（IIR）和有限长冲激响应滤波器（FIR），FIR的突出优点是：

系统总是稳定的、易于实现线性相位、允许设计多通带（或多阻带）滤波器，但与IIR相比，在满足同样阻带衰减的情况下需要的阶数较高，滤波器的阶数越高，占用的运算时间越多，因此在满足指标要求的情况下应尽量减少滤波器的阶数。

2．函数设计法

窗函数设计法是一种通过截短和计权的方法使无限长非因果序列成为有限长脉冲应响应序列的设计方法，通常在设计滤波器之前，应该先根据具体的工程应用确定滤波器的技术指标，在大多数实际应用中，数字滤波器常常被用来实现选频操作，所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。

用窗函数法设计FIR滤波器的步骤如下：

（1）根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N（或阶数M=N-1），窗函数类型可根据最小阻带衰减As独立选择，因为窗口长度N对最小阻带衰减As没有影响，在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N，设待求滤波器的过渡带宽为Δω，它与窗口长度N近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N，在N和窗函数类型确定后，即可调用MATLAB中的窗函数求出窗函数W（n）。

（2）根据待求滤波器的理想频率响应求出理想单位脉冲响应hd（n），如果给出待求滤波器频率应为Hd（ejw），则理想的单位脉冲响应可以用下面的傅里叶反变换式求出：

（13）

在一般情况下，hd（n）是不能用封闭公式表示的，需要采用数值方法表示；

从w=0到w=2π采样N点，采用离散傅里叶反变换（IDFT）即可求出。

（3）计算滤波器的单位脉冲响应h（n），它是理想单位脉冲响应和窗函数的乘积。

（4）算技术指标是否满足要求，为了计算数字滤波器在频域中的特性，可调用freqz子程序，如果不满足要求，可根据具体情况，调整窗函数类型或长度，直到满足要求为止。

使用窗函数法设计时要满足以下两个条件：

窗谱主瓣尽可能地窄，以获得较陡的过渡带；

尽量减少窗谱的最大旁瓣的相对幅度，也就是使能量尽量集中于主瓣，减小峰肩和纹波，进行增加阻带的衰减。

3、几种数字滤波器的窗函数设计

3.1数字低通滤波器的窗函数设计

假设理想低通数字滤波器的频率响应

，其幅频特性为

那么该滤波器的：

式中

表示截止频率（rad），

表示采样延迟。

由此可得理想数字低通滤波器的单位冲激响应

为：

（14）

根据式

（1），可知

是一个以α为中心的偶对称的无限长非因果序列，如果截取一段n=0～N-1的hd（n）作为h（n），则为要保证所得到是线性相位FIR滤波器，延时α应为h（n）长度N的一半,

α=（N-1）/2

（15）

其中WR（n）=RN（n）（为了表示一致）

设W（ejω）为该窗口函数的频谱:

用幅度函数和相函数来表示，则有

其线性相位部分

e-jωα则是表示延时一半长度α=（N-1）/2，对频响起作用的是它的幅度函数

图1矩形窗函数及其幅度函数

理想频响也可以写成幅度函数和相位函数的表示形式Hd（ejω）=Hd（ω）e-jωα

其中幅度函数为

二个信号时域乘积对应于频域卷积，所以有

如果也以幅度函数H（ω）和相位函数来表示H（ejω）,H（ejω）=H（ω）e-jw则实际FIR滤波器的幅度函数H（ω）为

正好是