第5章 数字滤波器基本结构
数字滤波器的基本结构 ppt课件

算子zw-11(表n) 示b0,x(n它) 表w5示(n)单 b位0x延(n)时 a。1y(n 1) a2 y(n 2)
y(n) w2 (n) w1(n)
y(n) a1 y(n 1pp)t课件a2 y(n 2) b0x(n)
6
第5章 数字滤波器的基本结构
5.2 IIR滤波器的基本结构
入的数字序列通过一定的运算后转变为另一组输出的
数字序列,因此它本身就是一台数字式的处理设备。
数字滤波器一般可用两种方法实现:1)根据描述数字滤
波器的数学模型或信号流图,用数字硬件装配成一台
专门的设备,构成专用的信号处理机;2)直接利用通用
计算机,将所需要的运算编成程序让计算机来执行,
即用软件来实现数字滤波器。
M
N
M
ak y(n k) bk x(n k)
bk x(nk1k) k 0
N
k 点 共(M+N)个延时单元
实现系统函数极点
图5-4 实现N阶差p分pt课方件 程的直接I型结构
9
第5章 数字滤波器的基本结构
二、直接Ⅱ型(典范型、正准型)结构
方框图表示法
信号流图表示法
图 5-1 基本运算的方框图表示及信号流图表示
ppt课件
5
第5章 数字滤波器的基本结构
二阶数字滤波器: y(n) a1y(n 1) a2 y(n 2) b0x(n)
源节点或 输入节点
阱节点或 输出节点
加法器
●
分支节点
输入支w2(路n) 的 y信(n)号值等于这一支路起点处节点信号值 乘值以,支www则354(((路认nnn))) 上为信 来 方aww1的其23w向号 代((3nn传(传,流 表n)11输有图一))输a向是条系2系yyw线((一支4数nn数(段n种路。)12为上)有,) 如a标1向箭1y,果注(图头n而出支的,1支延)方它路路向用迟a上2的代箭y支不(传n表头路标输信的2)值则传号有。用输流向动线延系的段数迟
第5章数字滤波器的基本结构

1、横截型(卷积型、直接型)
差分方程:
2、级联型
将H(z)分解成实系数二阶因式的乘积形式:
级联型的特点
• 每个基本节控制一对零点,便于控制滤波器的 传输零点
• 系数比直接型多,所需的乘法运算多
3、频率抽样型
N个频率抽样H(k)恢复H(z)的内插公式:
子系统: 是梳状滤波器
在单位圆上有N个等间隔角度的零点:
5.3 FIR数字滤波器的基本结构
• FIR数字滤波器的特点: 系统函数:
有N-1个零点分布于z平面 z=0处 是N-1阶极点
1)系统的单位抽样响应 h(n)有限长,设长度为N
2)系统函数H(z)在
处收敛,有限z平面只
有零点,全部极点在 z = 0 处(因果稳定系统)
3)无输出到输入的反馈,一般为非递归型结构
• 原网络中所有支路方向倒转,并将输入x(n)和 输出y(n)相互交换,则其系统函数H(z)不改变。
例:设IIR数字滤波器差分方程为:
试用四种基本结构实现此差分方程。 解:对差分方程两边取z变换,得系统函数:
得直接Ⅰ型结构:
典范型结构:
将H(z)因式分解: 得级联型结构:
将H(z)部分分式分解: 得并联型结构:
频率响应:
子系统:
单位圆上有一个极点:
与第k个零点相抵消,使该频率 率响应等于H(k)
Hale Waihona Puke 处的频频率抽样型结构的优缺点
• 调整H(k)就可以有效地调整频响特性
• 若h(n)长度相同,则网络结构完全相同,除了 各支路增益H(k),便于标准化、模块化
• 有限字长效应可能导致零极点不能完全对消, 导致系统不稳定
对其进行傅氏变换得:
数字滤波器的基本结构

H (z)
A
m1 N1
m1 N2
(1 ck z1) (11k z1 2k z2 )
k 1
k 1
将单实根因子看作二阶因子的特例:
46
M 1 2
(1 1m z1 2m z2 )
H (z) A m1 N 1 2 (1 1k z1 2k z2 ) k 1
:表示取整。
其中
Hi
(z)
1 1i z1 11i z1
2i 2i
z 2 z 2
,
级联结构:
i 0,1,..., m
X(n) H1(Z)
H2(Z)
。。。
Y(n) Hm(Z)
48
H(Z)的实现结构即可表示为基本二阶节 的级联形式。每个二阶节用典范型实现:
Z-1
Z 1 a1
y(n 1)
Z 1
a2
y(n 2)
Z 1 bM
x(n M )
Z 1
aN 1
y(n N 1)
Z 1
aN
y(n N)
实现N阶差分方程的直接I型结构
36
M=N
37
1)可直接差分方程或系统函数的标准形式画 出。两个网络级联:第一个横向结构M节延 时网络实现零点(分子,输入),第二个有 反馈的N节延时网络实现极点(分母,输 出) 。需要N+M级延时单元。
32
◦ 系统函数 ◦ 差分方程
M
bk z k
H(z)
k 0 N
1 ak zk
Y (z) X (z)
k 1
N
M
y(n) ak y(n k) bk x(n k)
数字信号处理 第五章

+ a2 z-1
数字信号处理—第五章
6
举例:二阶数字滤波器
y ( n ) a 1 y ( n 1) a 2 y ( n 2 ) b 0 x ( n )
x(n) b0 +
-1 a1 z
y(n)
+ a2 z-1
数字信号处理—第五章
7
举例:二阶数字滤波器
y ( n ) a 1 y ( n 1) a 2 y ( n 2 ) b 0 x ( n )
z z
2 2
H (z)
1 1k z 1 1k z
1 1
x(n)
H 1(z)
y (n )
H 2(z)
H k (z)
数字信号处理—第五章
22
数字信号处理—第五章
23
IIR数字滤波器的级联型结构优点
1) 每个二阶或一阶子系统单独控制零、极点。 2)级联顺序可交换,零、极点对搭配任意,因此级联 结构不唯一。有限字长对各结构的影响是不一样的, 可通过计算机仿真确定子系统的组合及排序。 3)级联各节之间要有电平的放大和缩小,以使变量值 不会太大或太小。太大可能导致运算溢出;太小可 能导致信噪比太小。 4)级联系统也属于最少延时单元实现,需要最少的存 储器,但乘法次数明显比直接型要多。 4)级联结构中后面的网络输出不会再流到前面,运算 误差积累比直接型小。
数字信号处理—第五章
4
基本单元(数字滤波器结构)有两种表 示方法
数字信号处理—第五章
5
举例:二阶数字滤波器
y ( n ) a 1 y ( n 1) a 2 y ( n 2 ) b 0 x ( n )
x(n) b0 +
数字信号处理课件(第五章 数字滤波器的基本结构)

F Ai 0i 1i z 1 H ( z) A0 1 1 pi z 1 1i z 1 2i z 2 (5-8) i 1 i 1 E
M 1
b
M
a
a
N 1
N
图 5-4 直接Ⅰ型结构
y (n N 1) 1 z y (n N )
第5章 数字滤波器的基本结构
5.2.2 直接Ⅱ型
直接Ⅱ型结构又称为典范型结构。由图5-4,直接Ⅰ型结构
的系统函数H(z)也可以看成是两个独立的系统函数的乘积。输入
信号x(n)先通过系统H1(z),得到中间输出变量y1(n),然后再把 y1(n)通过系统H2(z)得到输出信号y(n)。 即
若系统函数H(z)的分子阶数和分母阶数相等,即M=N时,其结构 如图5-5示。
输入信号x(n)先经过反馈网络H2(z),得到中间输出变量
y2 (n) ai y2 (n i ) x(n )
i 1
N
然后,将y2(n)通过系统H1(z),得到系统的输出y(n)
y (n) bi y2 (n i )
z-1 x(n-N)
z-1 y(n-N)
图 5-4 直接Ⅰ型结构
…
第5章 数字滤波器的基本结构
x(n) a1 a2 z-1 z-1
y2 (n) y2 (n-1) y2 (n-2) z-1 z-1
b0 b1 b2
y(n)
… … …
b N-1 bN z-1
…
a N-1
…
-1 aN z
图 5-5 直接Ⅰ型的变形结构
对应的差分方程为
第五章数字滤波器的基本结构

zk
WNk
j2k
e N
也就是说,此一阶网络在频率为
2
k
处响应为无穷大。
N
缺点:极点在单位圆上,由系数WN-k决定,当系数 量化时,极点会移动,若极点移到Z平面的单位圆
外,系统就不稳定了。
频率抽样结构修正:即将所有零极点都移动到单位 圆内某一靠近单位圆、半径为r的圆上。
四、快速卷积结构
利用“时域序列的圆周卷积等效于频域的离散频谱
k 1
当M=N时,H(z)表示为
H (z) G 0 k N 1 11 A c k kz 1 k N 2 11 1 0 k kz 1 1 kz 2 1 kz 2
1 1 az 1
r0 r1z1 1a1z1 a2z2
并联型特点: 各并联基本节间的误差相互没有影响,比级联型 的误差稍小。
转置定理: 如果将原网络中 所有去路方向倒 转,并将输入 x(n)和输出y(n) 相互交换,则其 系统函数H(z)不 变。
例:已知一个因果线性移不变滤波器的系统函数为
10.87 z15 H (z)(10.2z10.9z2)1(0.7z1)
画出该系统以下形式的流图。 (1)直接Ⅰ型 (2)直接Ⅱ型 (3)级联型 (4)并联型
(A
C)
( B 0.2C 0.7 A ) z 1 (0.9C (1 0 .2 z 1 0 .9 z 2 )(1 0 .7 z 1 )
0.7 B ) z 2
AC 1
B
0 .2 C
0.7 A
0 .875
0 .9 C 0 .7 B 0
A 0 .2794
B 0 .9265
C 0 .7206
mm,2m N
to
2(m1)
数字滤波器的基本结构
1 2 cos( 2 )z 1 z 2
N 1
实系数频率取样型结构流图
x[k] zN
1/N y[k]
1
z1
1
2 cos( 2 ) z1
N
1
z1 2 cos( 2 )
N
优点:1. H[m]零点较多时,实现较为简单。
2. 可以构成滤波器组,实现信号的频谱分析。
k 0
x[k]
z 1
z 1
z 1
1
1
1
z 1
z 1
1 z 1
z 1 1
h[0]
h[1]
h[2]
y[k]
h[ M 3] h[ M 1]
2
2
相同系数的共用乘法器,只需(M+1) /2个乘法器
三、 FIR 数字滤波器的级联型结构
将H(z)分解为若干个实系数一阶二阶因子相乘
第5章 数字滤波器的基本结构
IIR数字滤波器的基本结构 FIR数字滤波器的基本结构 格型结构
IIR数字滤波器的基本结构
直接型结构 级联型结构 并联型结构
一、IIR数字滤波器的直接型结构
Y (z)
M
bi z i
H2(z) W (z)
H(z)
i0 N
1 a j z j
第 p-1阶
e bp2 [k ]
e1f [k]
e0f [k]
第1阶
e1b [k ]
e0b [k ]
cp
c p1
c p2
c1
图中的方框是如下基本格型单元
c0 y[k]
e
f p
[k
]
e
数字滤波器的结构
y n
z 1
b2
b1
x n
z
1
z
1
z
a1 a0
1
aN
z
1
z
1
z
1
x n 1
b1 b2
bN 1
y n 1
z
1
aN 1
y n
x n
系数aibi 存储器
y n
z 1
z 1
运算器
z 1 z 1
输 出 寄 存 器
z 1
控制器 图5.1 硬件结构数字滤波器
z 1
4
第五章 数字滤波器的结构
同样这个运算也可以在 通用计算机上实现。 以一阶数字滤波器为例:
x 1 0 y 1 0
输入a0 b1 {x(n)}
8
第五章 数字滤波器的结构
x n (1) 直接型 一个N阶IIR滤波器的传递函数可 x n 1 以表达为
a0
z 1 z 1
y n
a1
a2
b1
b2
z 1 z 1
y n 1
y n 2
H ( z)
a z
i 0 i N i 1
N
1
x n 2
② a ⑥ b ⑤ 1 1
图5.5一阶数字滤波器的信号流图表达
可以看到,用信号流图表达数字网络的结构可以更简洁,我们在下面 将普遍采用信号流图的办法来分析数字滤波器的结构。
7
第五章 数字滤波器的结构
运算结构的不同将会影响系统的精度、误差、稳定性、经济性 以及运算速度等许多重要的性能。对于无限长单位脉冲响应(IIR)滤 波器与有限长单位脉冲响应(FIR)滤波器,它们在结构上各自有自己 不同的特点,下面将对它们分别加以讨论。
第5章数字滤波器结构
二阶因式表示存在复共轭根。为了简化级联形式, 特别是在时分多路复用时,采用相同形式的子网络 结构就更有意义,因而将实系数的两个一阶因子组 合成二阶因子,则整个H(z)就可以完全分解成实系 数的二阶因子的形式:
H (z)A 1 1 1 1k kz z 1 1 2 2k kz z 2 2A H k(z)
第5章 数字滤波器的结构
内容与重点:
➢ 内容:
数字滤波器结构的表示方法 IIR数字滤波器结构 FIR数字滤波器结构 数字滤波器的格型结构
➢ 重点
表示方法 IIR滤波器结构
• 直接型(I型、II型)、级联型、并联型
FIR滤波器结构
• 直接型、级联型、线性相位型
5.1 数字滤波器结构的表示方法
数字滤波器是一个数字信号处理系统。假 设数字滤波器系统将输入序列x(n),经过处理 后,得到输出序列y(n),用T[.]表示数字滤波 器的处理,则处理过程可以表示为:
y(n)T[x(n)]
若数字滤波器为线性系统,则可以用常系数差分方程表示 该处理过程:
N
M
y(n) aky(nk) bkx(nk)
k1
a22
z-1 β12
β22
y(n)
a13 z-1
a23
z-1 β13
β23
➢ 图5-16 六阶IIR滤波器的级联结构
26
➢ 例5-2:设IIR数字滤波器的系统函数为:
822z127z218z3 H(z)18z115z256z3
➢ 画出该滤波器的级联型结构。
27
➢ 解: 将H(z)的分子、分母进行因式分解,
将H(Z)展成部分分式形式:
H (Z ) k N 1 1 1 c A k k Z 1 k N 2 1(1 d B k k Z (1 1 )g 1 k ( Z d 1 k * ) Z 1 ) M k 0 N G k Z k
第5章_数字滤波器的基本结构
1.系统函数因式分解
一个N阶系统函数可用它的零、极点来表示。
M M
i 1 H (z) i0 A N N i 1 1 a z ( 1 d z i ) i i 1 i 1
bz
i
i
( 1Cz
i
1
)
将系统函数进行进一步分解,使分子、分母中 每个因式的次数不高于2,这样可以使各项系数都 是实数。
0 H ( z ) i H ( z ) H ( z ) H ( z ) H ( z ) 1 2 2 1 N i 1 a iz i 1
M
i b z i
M
其中: H 1 ( z )
i0
bi z i 1
N
H 2(z)
1
i 1
a i z i
x(n) z-1 z-1 z-1 z-1
b0 b1 b2 b M+1 bM 第一部分 对调 a1 a2 a N-1 aN
y(n) z-1 z-1 z-1 z-1
x(n) a1 a2 对调 a N-1 aN z-1 z-1 z-1 z
b0 z-1 b1 z-1 b2
y(n)
z-1 b M+1 bM
-1 z -1
第二部分
由于对调后前后两路都有一条内容完全相同的 延时链,可以合并为一条即可。
§5.3 无限长脉冲响应基本网络结构
主要特点:
①系统的单位冲激响应h(n)是无限长的(n→∞); ②系统函数H(z)在有限z平面上(0<|z|<∞)有极点存在。 ③结构上存在着输出到输入的反馈,即结构是递归的。 ④因果稳定的IIR滤波器其全部极点一定在单位园内。 同一种系统函数H(z)可以由多种不同结构,它的基本 结构有:直接型、级联型、并联型
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第5章数字滤波器的基本结构5.1 数字滤波器的结构特点与表示方法5.2 IIR滤波器的结构5.3 FIR滤波器的结构5.2 IIR滤波器的结构5.2.1 直接型直接型((Ⅰ型)一个N阶的IIR滤波器的输入输出关系可以用如式(5-1)所示的N阶的差分方程来描述。
把式(5-1)重写如下:∑∑= =−+−=Ni iMiiinyanxbny1) ()1()(从这个差分方程表达式直接可以看出,系统的输出y (n )由两部分构成:第一部分是一个对输入x (n )的M 阶延时链结构,每阶延时抽头后加权相加,构成一个横向结构网络。
第二部分是一个对输出y (n )的N 阶延时链的横向结构网络,是由输出到输入的反馈网络。
由这两部分相加构成输出,取M =N (当然M 可不等于N )可得结构图如图5-2。
从图上可以看出,直接Ⅰ型结构需要2N 个延时器和2N +1个乘法器(M=N 情形下)。
)(0i n x b Mi i −∑=)(1i n y a Ni i −∑=图5-2 直接Ⅰ型结构z -1z -1z -1…b N -1bNb2b1b 0x (n )x (n -1)x (n -2)x (n -N )z -1z -1z -1…a N -1a N a 2a 1y (n )y (n -1)y (n -2)y (n -N )…………5.2.2 直接直接ⅡⅡ型直接Ⅱ型结构又称为典范型结构。
由图5-2,直接Ⅰ型结构的系统函数H(z)也可以看成是两个独立的系统函数的乘积。
输入信号x(n)先通过系统H1(z),得到中间输出变量y 1(n),然后再把y1(n)通过系统H2(z)得到输出信号y(n)。
即∑∑=−=−−==NiiiMiiizazbzHzHzH1211)()()(式中,∑=−=Mi ii zb z H 01)(对应的差分方程为:∑∑=−=−=−=Ni ii Mi i za z H i n xb n y 120111)()()(对应的差分方程为)()()(11n y i n y a n y Ni i +−=∑=假设所讨论的IIR 数字滤波器是线性非时变系统,显然交换H 1(z )和H 2(z )的级联次序不会影响系统的传输效果,即)()()()()(1221z H z H z H z H z H ==若系统函数H (z )的分子阶数和分母阶数相等,即M=N (当然二者同样可以不相等)时,其结构如图5-3所示。
输入信号x (n )先经过反馈网络H 2(z ),得到中间输出变量)()()(122n x i n y a n y Ni i +−=∑=然后,将y 2(n )通过系统H 1(z ),得到系统的输出y (n ))()(02i n y b n y Mi i −=∑=结构图5-3中有两条完全相同的对中间变量y(n)进行延2迟的延时链,我们可以合并这两条延时链,得到如图5-4所示的直接Ⅱ型结构(图中取M=N)。
比较图5-2和图5-4可知: 直接Ⅱ型比直接Ⅰ型结构延时单元少,用硬件实现可以节省寄存器,比直接Ⅰ型经济;若用软件实现则可节省存储单元。
但对于高阶系统直接型结构都存在调整零、极点困难,对系数量化效应敏感度高等缺点。
图5-3 直接Ⅰ型的变形结构z -1z -1z -1…a N -1a Na 2a 1z -1z -1z -1…b N -1bNb2b10y 2(n )y 2(n -1)y 2(n -2)y 2(n -N )…………图5-4 直接Ⅱ型结构z -1z -1z -1…a N -1a Na 2a 1…b N -1b Nb 2b 10………5.2.3 级联型若把式(5-2)描述的N 阶IIR 滤波器的系统函数H (z )的分子和分母分别进行因式分解,得到多个因式连乘积的形式∏∏∑∑=−=−=−=−−−=−=Ni iM i iN i ii Mi ii zd zc Aza zb z H 11111)1()1(1)((5-4)式中:A 为常数,c i 和d i 分别表示H (z )的零点和极点。
若H (z )的分子和分母都是实系数多项式,而实系数多项式的根只有实根和共轭复根两种情况。
将每一对共轭零点(极点)合并起来构成一个实系数的二阶因子,并把单个的实根因子看成是二次项系数等于零的二阶因子,或者将两个实系数的一阶因子组合成一个二阶因子,则可以把H (z )表示成多个实系数的二阶数字网络H j (z )的连乘积形式,如式(5-5)所示:∏==Kj j z H A z H 1)()((5-5)式中:2211221101)(−−−−−−++=zz zz z H j j j j j j ααβββ若每一个实系数的二阶数字网络的系统函数H j (z )的网络结构均采用前面介绍的直接Ⅱ型结构,则可以得到系统函数H (z )的级联型结构,如图5-5所示。
图5-5 级联型结构x (n )y (n )z -1z -1α11α21β11β21β01……z -1z -1α1Kα2Kβ1K β2Kβ0KA在级联型结构中,每一个一阶网络只关系到滤波器的一个零点、一个极点;每个二阶网络只关系到滤波器的一对共轭零点和一对共轭极点。
调整系数β0j 、β1j和β2j只会影响滤波器的第j对零点,对其他零点并无影响;同样, 调整分母多项式的系数α1j 和α2j也只单独调整了第j对极点。
因此,与直接型结构相比,级联型结构便于准确地实现滤波器零、极点的调整。
此外,因为在级联结构中,后面的网络的输出不会流到前面,所以其运算误差也比直接型小。
5.2.4 并联型把传递函数H (z )展开成部分分式之和的形式,就可以得到滤波器的并联型结构。
∑∑∑∑∑−=−=−−−=−=−=−+−−−+−=−=N M k kkN k k k k k N k k k N k kk Mk kkzGz d z d z g B z c A za zbz H 011*11111021)1)(1()1(11)(这一公式是最一般的表达式。
公式中N=N 1+2N 2。
当M 〈N 时,上述公式中不包含最后一个求和项;如果M=N ,则上述公式中的最后一个求和项变为G 0一项。
一般IIR 滤波器都满足M 小于或等于N 的条件,则在该条件下,系统是由N 1个一阶子系统、N 2个二阶子系统等并联组成的,而这些子系统都可以采用典范性结构来实现。
当M=N 时,其表达式可进一步写成公式(5-6),其中A 0=G 0,p i =c i ,E=N 1,F=N 2。
∑∑=−−−=−−−++−+=Fi i i i i Ei i izz z z p A A z H 1221111011011)(ααγγ(5-6)由式(5-6)知,滤波器可由N1个一阶网络、N2个二阶网络和一个常数支路并联构成,其结构如图5-6所示。
并联型结构也可以单独调整极点位置,但对于零点的调整却不如级联型方便,而且当滤波器的阶数较高时,部分分式展开比较麻烦。
在运算误差方面,由于各基本网络间的误差互不影响,没有误差积累,因此比直接型和级联型误差稍小一点。
图5-6 并联型结构…x (n )α11z -1α21z-1γ01γ11α1F z -1α2Fz -1γ0F γ1F…z -1A 1p 1y (n除了以上一些基本结构外,还有一些其他结构,这取决于线性信号流图理论中的多种运算处理方法。
当然各种流图都保持输入到输出的传输关系不变,即H(Z)不变。
其中有一种方法称为流图的转置,这种方法是利用了流图的如下转置定理:转置定理:若将线性移不变网络中的所有支路方向倒转,转置定理并将输入x(n)和输出y(n)相互交换,则其系统函数H(Z)不改变。
5.3 FIR 滤波器的结构5.3.1 直接型直接型 设FIR 数字滤波器的单位脉冲响应h (n )的长度为N ,其传递函数和差分方程分别为:∑−=−=10)()(N n nzn h z H (5-7)∑−=−=1)()()(N m m n x m h n y (5-8)根据式(5-7)或式(5-8)可直接画出如图5-7所示的 FIR 滤波器的直接型结构。
由于该结构利用输入信号x (n )和滤波器单位脉冲响应h (n )的线性卷积来描述输出信号y (n ),所以FIR 滤波器的直接型结构又称为卷积型结构,有时也称为横截型结构。
当然也可以利用转置定理来获得转置型直接结构。
图5-7 FIR 的直接型结构z -1x (n )h (0)h (1)z -1h (2)……h (N -3)z -1h (N -2)z -1h (N -1)y (n )5.3.2 级联型当需要控制系统传输零点时,将传递函数H (z )分解成二阶实系数因子的形式:∑∏−==−−−++==10122110)()()(N n Mi i i i nz a z a a z n h z H 图5-8 FIR 的级联型结构x (n )y (n )z -1z -1a 11a 21a 01……z -1z -1a 12a 22a 02z -1z -1a 1Ma 2Ma 0M。