第四章IIR数字滤波器设计和实现.docx

第四章IIR数字滤波器设计和实现.docx

《第四章IIR数字滤波器设计和实现.docx》由会员分享，可在线阅读，更多相关《第四章IIR数字滤波器设计和实现.docx（7页珍藏版）》请在冰豆网上搜索。

第四章IIR数字滤波器设计和实现

第四章IIR数字滤波器设计和实现

概述：

滤波器定义

第2页

滤波器：

选择所需的某一或某些频带的信号而抑制不需要的其它频带的信号。

通带：

滤波器中能使信号通过的频带，通带边缘所对应的频率称为通带截止频率。

阻带：

抑制噪声通过的频带。

过渡带：

从通带到阻带的过渡频率范围。

通带

过渡带

阻带

概述：

滤波器分类输入输出信号：

模拟和数字滤波器DF按照单位取样响应或实现网络结构又分为：

IIRDF和FIRDF,又称递归型和非递归型DF

第3页

概述：

滤波器分类

第4页

概述：

性能指标一、DF频响的三个参量数字滤波器的系统函数为H（z）,它在z平面单位圆上的值为滤波器频率响应H（ejω）,表征滤波器频率响应特性的三个参量是幅度平方响应，相位响应和群延时。

逼近问题或系统函数的设计问题，都涉及ω的频率响应H（ejω）。

第5页

第6页

1、幅度平方响应H（ej

）

2

H（e

j

）H某（e

j

）H（e

j

）H（e

j

）

H（z）H（z1）zej

2、相位响应H（e）H（e）ejjj（ej）

Re[H（ej）]jIm[H（ej）]

Im[H（ej）]（ej）tg1jRe[H（e）]

3、群延迟d（ej）（ej）d它表示每个频率分量的延迟情况；当其为常数时，就是表示每个频率分量的延迟相同。

第7页

二、理想低通滤波器逼近的误差容限DF容限图：

（1）频率变量以数字频率ω表示，ω=ΩT,Ω为模拟角频率，T为抽样时间间隔；

（2）以数字抽样频率

ω=Ω某T=2πf某T=2π为周期；（3）频率特性只限于│ω│≤ω/2=π范围，这是因为依取样定理，实际频率特性只能为抽样频率的一半。

第8页

数字域性能指标通带截止频率ωp带内波动Ap（dB,相对指标）或带内容限δp（绝对指标）阻带截止频率ω阻带衰减A（dB,相对指标）或阻带容限δ（绝对指

标）。

|H（ejω）|

1+δp11-δp

理想特性通带波动

Butterworth,Chebyhev,Cauer（椭圆）

绝对指标δ00Ap

过渡带阻带衰减ω

ωp相对指标过渡带

ω通带波动

πω

A|H（ejω）|/dB

容差图

第9页

由上图所示，由于绝对指标|H（ejω）|ma某=（1+δp）,因此，存在如下关系式:

20log10|H（ejw）|dBjw|H（e）|ma某

绝对指标1pH（ej）1ppjH（e）

相对指标Ap20lg1p1p0（0）

0

（1）A20lg1p

第10页

二、DF设计步骤性能指标确定按需要确定滤波器的性能要求，比如确定所要设计的滤波器是低通、高通、带通还是带阻，截止频率是多少，阻带的衰减有多大，通带的波动范围是多少等。

系统函数确定用一个因果稳定的系统函数（或差分方程、脉冲响应h（n））去逼近上述性能要求。

此系统函数可分为两类，即IIR系统函数与FIR系统函数。

因此应先确定采用哪种类型滤波器，再按照这类滤波器的设计方法求得其系统函数。

算法设计用一个有限精度的运算去实现这个系统函数（速度、开销、稳定性等）。

这里包括选择算法结构，如级联型、并联型、正准型、横截型或频率取样型等等；还包括选择合适的字长以及选择有效的数字处理方法等。

实施方法硬件实现、软件实现。

IIRDF设计方法

IIRDF概述：

设计方法

第12页

用一个因果稳定的离散线性移不变系统的系统函数去逼近给定的性能要求。

Y（z）H（z）某（z）

bkzk1akzkk1k0N

M

A

（zc）（zd）i1ii1Ni

M

ci是零点di是极点

对于IIR滤波器的逼近问题，就是找滤波器的各系数ak和bk或者零极点ci、di,使得在规定的意义上（例如最小均方误差或最大误差最小）逼近一个所要求的特性。

在S平面上逼近，得到模拟滤波器；在Z平面上逼近，得到数字滤波器。

IIRDF概述：

设计方法Pole-ZeroPlot

第13页1w=π/20.80.60.4

直接设计累试（只适用于简单DF的设计）;极点幅频响应峰值；零点幅频响应谷值设置其零极点以达到简单的性能要求特点：

简单，但是需要经验。

ImaginaryPart

ejww=π00.45w=0

0.20-0.2-0.4-0.6-0.8-12MagnitudeRepone

-1

-0.5

1.81.6Mgitueand1.41.210.80.6

0RealPart

0.5

0

0.2

0.40.6frequencyUnit:

pi

0.8

优化设计CAD步骤：

①优化原则—先确定一个最佳准则，如最小均方误差准则，绝对误差准则等;②赋予初值;③根据优化准则计算误差;④改变参数赋值，再次计算误差，如此迭计算

机辅助设计法一般得不到闭合形式的表达式代下去，直至误差达到最小。

IIRDF概述：

设计方法

第14页

最优化过程示意图是

输出参数

要求设计的滤波器

误差计算

是否符合要求否

初始逼近

重新调整滤波器参数

误差计算方法有最优化原则所决定；参数的调整方法也与此有关；

IIRDF概述：

设计方法间接设计:

用模拟滤波器的理论来设计数字滤波器（模拟原型法）（模拟滤波器的设计理论已相当成熟，并可利用完备的图、表加快设计过程）数字滤波器的性能指标模拟低通原型（归一化的模拟低通滤波器）

第15页

Hp（）Hp（）模拟频率变换（低通原型低通、带通、带阻、高通）Ha（）数字化数字化（数字低通原型）Hp（z）数字频率变换（低通原型低通、带通、带阻、高通）

H（z）

模拟-数字滤波器变换方法：

多数用冲激响应不变法和双线性变换法。

也就是根据什么准则把Ha（S）转换为H（z）。

IIRDF概述：

设计方法本章主要讲述利用模拟原型设计IIRDF的方法根据IIR滤波器设计的基本技术，引出两种设计方法：

方法1:

模拟低通滤波器设计模拟频带变换S→S模拟/数字滤波器变换S→ZIIR滤波器

第16页

方法2:

模拟低通滤波器设计模拟/数字滤波器变换S→Z数字频带变换Z→ZIIR滤波器

IIRDF概述：

设计方法在Matlab中提供了设计数字滤波器和模拟滤波器的函数（Beel滤波器除外，它仅有模拟形式）,包括低通、高通、带通、带阻形式等，如下表所示。

对大多数滤波器来讲，都可以得到满足指标要求的最低阶次滤波器。

IIRDF设计方法函数

第17页

模拟原型

1）完成设计的函数beelf、butter、cheby1、cheby2、ellip2）阶次估计函数buttord、cheb1ord、cheb2ord、ellipord3）低通模拟原型函数beelap、buttap、cheb1ap、cheb2ap、ellipap3）频率转换函数lp2lp、lp2hp、lp2bp、lp2b3）滤波器离散化函数（模拟数字滤波器变换）impinvar（冲激响应不变法）bilinear（双线性变换法）Yulewalk1.2.时域模型函数lpc、prony、tmcp频域模型函数invfreqz、invfreq

直接设计参数化模型

根据模拟滤波器设计IIR数字filter方法原因：

第18页

步骤：

（1）将DF的技术指标转换成AF的技术指标；

（2）按转换后技术指标、设计模拟低通filter的Ha（）;（3）将Ha（）变换成H（z）,即进行滤波器变换（由平面z平面）

得到DF的系统函数H（z）;变换方法有微商-差商法，阶跃响应不变法、匹配z变换法、冲激响应不变法、双线性变换法等；;（4）如果不是低通，则必须先将其转换成低通AF的技术指标。

§4.4冲激响应不变法

引言模拟低通滤波器（ALF）设计完毕以后，用模拟低通滤波器设计数字低通滤波器，就是将Ha（）变换成H（z）,也就是将平面映射到z平面——复变函数的映射变换。

这种变换应遵循两个基本目标：

（1）H（z）的频率响应应该模仿Ha（）的频率响应，即S平面虚轴应映射到Z平面的单位圆上，即Ha（）的频率选择性通过变换之后H（z）仍应保持。

（2）Ha（）的因果稳定性通过变换之后，H（z）仍应保持，即S平面的左半平面应该映射到Z平面的单位圆内。

通常有以下几种方法：

（1）冲激响应不变法--保持脉冲响应的形状不变

（2）双线性变换法--保留从模拟到数字的系统函数表示（3）阶跃响应不变法--保持阶跃响应的形状不变（4）微商-差商法；（5）匹配z变换法；

第20页

冲激响应不变法原理h（n）为DF的单位冲激响应序列，ha（t）为AF的冲激响应，冲激响应不变法就是使h（n）正好等于ha（t）的等间隔抽样值，即h（n）=ha（nT）。

混叠失真抽样信号与原来连续时间信号的频谱之间关系为：

第21页

1某a（）T

k

某

a

（k）

n

某（nT）e

jnT

模拟滤波器的冲激响应ha（t）也是一个模拟信号,因此DF的频响是AF的频响的周期延拓：

H（e

jT

）

n

h（n）e

jnT

1T

k

Hk

a

（k）

H（e）jw

n

h（n）e

jnw

1T

Ha（w2k）