基于Matlab的IIR数字滤波器设计毕业论文Word下载.docx

基于Matlab的IIR数字滤波器设计毕业论文Word下载.docx

《基于Matlab的IIR数字滤波器设计毕业论文Word下载.docx》由会员分享，可在线阅读，更多相关《基于Matlab的IIR数字滤波器设计毕业论文Word下载.docx（37页珍藏版）》请在冰豆网上搜索。

第二章IIR数字滤波器设计方法 ５

第一节IIR数字滤波器的设计步骤 ５

第二节用脉冲相应不变法设计IIR数字滤波器 ６

一、设计原理 ６

二、脉冲响应不变法优缺点 ８

第三节双线性变换法设计IIR数字滤波器 ９

一、设计原理 ９

二、双线性变换法优缺点 １１

第三章IIR滤波器的MATLAB设计 １３

第一节IIR数字滤波器的典型设计法 １４

第二节IIR数字滤波器的直接设计法 １８

第三节FDATool介绍和界面设计 ２３

第四节FDATOOL设计IIR数字滤波器 ２４

第五节SIMULINK 仿真IIR滤波器 ２６

总结 ２９

致谢.............................................错误！

未定义书签。

参考文献 ３０

结束语............................................错误！

２

前言

随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。

目前数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。

在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，DigitalFilter），根据其单位冲激响应函数的时域特性可分为两类：

无限冲激响应IIR（InfiniteImpulse

Response）滤波器和有限冲激响应FIR（FiniteImpulseResponse）滤波器。

与FIR滤波器相比，IIR的实现采用的是递归结构，极点须在单位圆内，在相同设计指标下，实现IIR滤波器的阶次较低，即所用的存储单元少，从而经济效率高。

MATLAB是英文MATrixLABoratory（矩阵实验室）的缩写。

它是美国的

MathWorks公司推出的一套用于科学计算和图形处理可视化、高性能语言与软件环境。

MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。

工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。

只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常方便。

­

１­

第一章数字滤波器

第一节数字滤波器的概念

滤波器是指用来对输入信号进行滤波的硬件和软件。

数字滤波器是对数字信号实现滤波的线性时不变系统。

数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。

数字滤波实质上是一种运算过程，实现对信号的运算处理。

数字滤波器和模拟滤波器相比，因为信号的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。

输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。

描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。

时域离散系统的频域特性:

Y（ejw）=X（ejw）H（ejw）



（式1-1）

其中Y（ejw）、X（ejw）分别是数字滤波器的输出序列和输入序列的频域特性（或

称为频谱特性）,

H（ejw）是数字滤波器的单位取样响应的频谱，又称为数字滤

波器的频域响应。

输入序列的频谱X（ejw）经过滤波后X（ejw）H（ejw）,因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H（ejw），使得滤波后的X（ejw）H（ejw）满足设计的要求，这就是数字滤波器的滤波原理。

第二节数字滤波器的分类

按照不同的分类方法，数字滤波器有许多种类，但总起来可以分成两大类：

经典滤波器和现代滤波器。

经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的频率成分占有不同的频带，通过一个合适的选频滤波器滤除干扰，得到纯净信号，达到滤波的目的。

但是，如果信号和干扰的频谱相互重叠，则经典滤波器不能有效地滤除干扰，最大限度地恢复信号，这时就需要现代滤波器，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。

现代滤波器是根据

２­

随机信号的一些统计特性，在某种最佳准则下，最大限度地抑制干扰，同时最大限度地回复信号，从而达到最佳滤波的目的。

经典数字滤波器从滤波特性上分类，可以分为：

低通滤波器、高通滤波器、带通滤波器、带阻滤波器。

Ha（jΩ）

高通

带阻

低通

带通

0 Ω 0 Ω

c Ω 0 Ω

图1­

1各种理想滤波器的幅频特性

IIR数字滤波器的特征是，具有无限持续时间冲激响应，需要用递归模型来实现，

其差分方程为：

N N

y（n）=å

aix（n-i）+å

biy（n-i）

（式1-2）

系统函数为：

i=0

i=1

H（z）=

M

-r

å

brZ

r=0

N

（式1-3）

k

1+å

aiZ-

k=1

设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H（z），使其频率响应H（z）满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。

３­

第三节数字滤波器的设计要求

滤波器的指标常常在频域给出。

数字滤波器的频响特性函数H（ejw）一般为复函数，所以通常表示为：

z=ejw

H（ejw）=H（z）| =

H（ejw）ejF（w）（式1­

4）

其中，|H（ejw）|称为幅频特性函数，Φ（w）称为相频特性函数。

幅频特性表示信号通过该滤波器后各频率成分的衰减情况，而相频特性反映各频率通过滤波器后在时间上的延时情况。

一般IIR数字滤波器，通常只用幅频响应函数

|H（ejw）|来描述设计指标，相频特性一般不作要求。

IIR滤波器指标参数如下图所示。

图中，ωp和ωs分别为通带边界频率和阻带边界频率；

δ1和δ2分别为通带波纹和阻带波纹；

允许的衰减一般用dB数表示，通带内所允许的最大衰减（dB）和阻带内允许的最小衰减（dB）分别为α

p和αs表示：

p

a=-20lg1-d1

1+d1

=20lg1+d1（式1­

5）

1-d1

一般要求：

as=-20lgd2

（式1­

6）

当0£

w£

wp时，-20lgH（e

jw）

£

ap；

jw

当w£

p时，as£

H（e 。

图1­

2低通滤波器的技术要求

４­

第二章IIR数字滤波器设计方法

IIR数字滤波器是一种离散时间系统，其系统函数为

H（Z）=

bkz-

k=0

1-å

akz-

=r（z）

x（z）

（式2-1）

假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个（M-N）的FIR子系统的级联。

IIR数字滤波器的设计实际上是求解滤波器的系数ak和bk，它是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近

系统的特性。

如果在S平面上去逼近，就得到模拟滤波器；

如果在z平面上去逼近，就得到数字滤波器。

第一节IIR数字滤波器的设计步骤

IIR数字滤波器的设计一般有两种方法：

一个是借助模拟滤波器的设计方法进行。

其设计步骤是，先设计模拟滤波器，再按照某种方法转换成数字滤波器。

这种方法比较容易一些，因为模拟滤波器的设计方法已经非常成熟，不仅有完整的设计公式，还有完善的图表供查阅；

另外一种直接在频率或者时域内进行，由于需要解联立方程，设计时需要计算机做辅助设计。

其设计步骤是：

先设计过渡模拟滤波器得到系统函数Ha（s），然后将Ha（s）按某种方法转换成数字滤波器的系统函数H（z）。

这是因为模拟滤波器的设计方法已经很成熟，不仅有完整设计公式，还有完善的图表和曲线供查阅；

另外，还有一些典型的优良滤波器类型可供我们使用。

为了保证转换后的H（z）稳定且满足技术指标要求，对转换关系提出两点要

求:

（1）因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。

（2）数字滤波器的频率相应模仿模拟滤波器的频响特性，s平面的虚轴映射为

z平面的单位圆，相应的频率之间呈线性关系。

利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：

５­

（1）确定数字低通滤波器的技术指标：

通带边界频率wp、通带最大衰减ap、阻带截止频率ws、阻带最小衰减as。

（2）将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。

（3）按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。

（4）用所选的转换方法，将模拟滤波器Ha（s）转换成数字低通滤波器系统函数H（z）。

IIR数字滤波器的设计流程图如下：

数字滤波器

技术指标

变换