数字滤波器的结构.
数字滤波器的基本结构(3)-sw_OK
5.2 无限长单位冲激响应(IIR)滤波器 的基本结构
一、直接I型
表述一个IIR滤波器的系统函数和差分方程分别 由(5-1)和(5-2)式表述,
M
N
y(n) bk x(n k) ak y(n k)
k 0
k 1
(5-2)
根据(5-2)式可以看出,y(n)可以分为两部分之和
M
第一部分为 bk x(n k) 对应输入x(n)及其各延迟 k 0
(2)将输入x(n)和输出y(n)互换位置。
18
5.2 无限长单位冲激响应(IIR)滤波器 的基本结构
x(n)
b0
b1
z 1 a1
b2
z 1 a2
y(n)
bM 1
bM
z 1
aN 1
z 1
aN
图8 直接 II 型的转置型
19
5.2 无限长单位冲激响应(IIR)滤波器 的基本结构
[例 1]设IIR数字滤波器的系统函数为
图6可以看作是图5的极点网络和零点网络互换级联 位置而成的。
观察图6
∵w1=w2 ∴前后两部分对应的延迟支路输出节点变量 也相等,即图中的w1(n-1)=w2(n-1),w1(n-i)=w2(n-i),
故可将前后两部分对应的延迟支路合并,合并后的信 号流图为
15
5.2 无限长单位冲激响应(IIR)滤波器 的基本结构
H (z) 8z3 4z2 11z 2
(z 1)(z2 z 1)
4
2
试画出该IIR数字滤波器的直接II型及其转置型的结构。
8 4z1 11z2 2z3 解: H (z) 1 5 z1 3 z2 1 z3
448
20
5.2 无限长单位冲激响应(IIR)滤波器 的基本结构
数字滤波器的基本结构
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 +
数字滤波器的基本结构
未来研究方向
新型算法研究
针对实际应用中的挑战,未来研究将进一步探索新型的数字滤波器 算法,以提高其性能、稳定性和适应性。
高性能硬件实现
随着集成电路和计算机工程的发展,未来研究将进一步探索高性能 、低功耗的数字滤波器硬件实现方法。
跨领域应用
数字滤波器在许多领域都有广泛的应用前景,如医疗、航空航天、环 保等,未来研究将进一步拓展数字滤波器的应用领域。
梯度下降法
通过迭代地更新滤波器的 系数,使得误差的梯度下 降最快,从而逐渐逼近最 优解。
牛顿法
利用牛顿定理,通过迭代 来寻找最优解,具有较高 的收敛速度和精度。
最优滤波器设计
最小均方误差(MMSE)滤波器
以最小化输出信号与期望信号之间的均方误差为优化目标,设计最优的滤波器 。
卡尔曼滤波器
一种递归滤波器,通过预测和更新来估计系统的状态,具有较高的稳定性和精 度。
控制系统
数字滤波器可以用于控制系统 的处理,如伺服控制、PID控制
、卡尔曼滤波等。
02
CHAPTER
数字滤波器的基本结构
数字滤波器的基本结构 直接形式
直接形式是数字滤波器的基本结构之 一。它是一种直观的形式,由一个输 入和一个输出组成,输入信号经过一 个或多个线性时不变系统后得到输出 信号。直接形式的结构简单,易于理 解和实现。
硬件优化
随着集成电路和计算机工程的发展,数字滤波器的硬件实 现越来越高效,低功耗、高速度和小型化成为主要趋势。
软件算法改进
数字滤波器的算法不断优化,以适应更复杂和多变的应用 场景,如神经网络、深度学习等算法的引入使得滤波效果 更加精确。
嵌入式应用
随着嵌入式系统的发展,数字滤波器在嵌入式设备上的应 用越来越广泛,这要求数字滤波器具有更强的稳定性和适 应性。
数字滤波器的基本结构ya
三、转置定理 如果将原网络中所有支路方向加以倒转,且将输入 和输出交换其系统函数仍不改变。
x(n)
a1
bb Z1 0 1
a2
b Z1 2
y(n)
bM1
a Z1
N1 bM
aN
Z1
(原网络)
y(n)
b0
b a1
Z1
1
a2
b Z1 2
bM1
x(n)
aN1
aN
b Z1 M
Z1
(转置后的网络)
5.3 FIR滤波器的基本结构
x(n) h(n) y(n)
y(n) x(n) h(n)
进行傅氏变换得:
Y (e j) X (e j) H (e j)
这种关系可用差分方程、单位冲激响应及系统函数进行描述。
X (e j )
0
H (e j )
0
Y(e j )
0
H(ejω)为矩形窗时 的情形
ωc
πω
ωc
πω
ω
πω
二、数字滤波器的系统函数与差分方程
a)输入节点或源节点x,(n) b)输出节点或阱 节点y(,n)
所处的节点; 所处的节点;
c)分支节点,一个输入,一个或一个以上输 出的节点;将值分配到每一支路;
d)相加器(节点)或和点,有两个或两个以 上输入的节点。
支路不标传输系数时,就认为其传输系数为1;
任何一节点值等于所有输入支路的信号之和。
y(n) a1 y(n 1) a2 y(n 2) b0 x(n)
x(n) b0
y(n)
b0x(n) a2 y(n 2)
Z 1
a1
a1y(n 1) Z 1
y(n 1)
数字信号处理第四章-数字滤波器的结构
3).H (z)
Y (z) X (z)
(1 bz1) (1 az1)
y(n) ay(n 1) x(n) bx(n 1)
9
10
11
w w
12
转置流图:
w(n) y(n)
原流图:
w(n) ay(n 1) x(n) bx(n 1) 两边作Z变换:
w(n) x(n) aw(n 1) y(n) w(n) bw(n 1) 两边作Z变换:
乘法系数为复数,运算量增加; 系统的稳定性依赖于零、极点相互抵消,对实
现的精度要求很高。在存在有限字长效应的情 况下,有可能造成系统不稳定。
54
确保所有零点、极点在单位圆内。 55
(h(n)为实数)
第k对 极点, 即第k 个与第 N-k个 谐振器 合并
56
谐振频 率不变
还有两点需要注意:(存在实根) 57
1
前言
线性时不变系统用单位冲击响应来表示 系统函数实际上单位冲击响应的Z变换 系统函数反映线性时不变系统的特性 大多数的信号处理可看成是对信号的滤波操作 数字滤波器实际上就是线性时不变系统
因此数字滤波器可以表示为:
2
前言
M
bk zk
H(z) Y(z) / X (z)
k 0 N
1 ak zk
从信号流图中:
可以清楚地看到系统中的运算步骤和运 算结构。FFT时用到了该特点。
运算结构可以直观反映所需的存储单元 和运算次数。由于是数字实现,必然存 在系统误差,运算结构同时也可以反映 系统误差的累积问题。 下面讨论的IIR和FIR滤波器结构将涉及 上述问题。
14
1
15
无限冲击响应滤波器的特点
82
第四章-数字滤波器的基本结构
将(4-7)式关系代入上式,得
H ( z)
N 11 2
h(n)
[zn
z(N 1n) ]
h(
N
1)
N 1
z2
(4-9)
n0
2
(4-8)(4-9)式中+号代表偶对称,-号代表奇对称。
当h(n)奇对称时,由于
h(n)
h(
N
1
n), 故h(
N 1) 2
0
下面的图19、图20分别画出N为偶数和N为奇数时 的线性相位FIR滤波器的结构。
W k N
WN( N k )
各并联支路的极点为
r
j 2 k
e N
,k
0,1, 2,
, N 1
为使系数为实数,可将共轭根合并,在z平面上 这些共轭根在半径为r的圆周上以实轴为轴成对 称分布,即 zN k zk
也就是 W (N k )
j 2 ( N k )
e N
(e
j
2 k N
)
WNk
27
4.3 有限长单位冲激响应(FIR)滤波器 的基本结构
级联型的每级对应一组由 (0i , 1i , 2i ) 参数决定的零点
6
4.3 有限长单位冲激响应(FIR)滤波器 的基本结构
三、线性相位的FIR滤波器结构: 在许多实际应用,如图像处理中,要求数字滤波器具
有线性相位 具有线性相位特性的滤波器传输函数H(ej)为
H(e j ) H() e j ()
则(4-12)式可写成:
1
N 1
H (z)
N
HC (z)
k 0
HK (z)
(4-13)
N 1
上式表明H(z)可看成是由 HC (z)和 HK (z) 两部分级 k 0
数字滤波器的基本结构
x[k ]
e0f [k]
e1f [k] K1
z 1 e0b [k ]
K1 z1 e1b [k ]
e
f p
1[k
]
ebp1[k ]
z 1
e2f [k] K2 K2 e2b [k ]
Kp Kp
e
f p
[k
]
y[k ]
Kp
z 1
Kp
ebp [k] yb[k]
e
f p
[k
]
反射系数
e
b p
[k
]
AZ系统的基本格形单元
第28页/共35页
三种滤波器的系统函数
全零点(AZ)滤波器 p A(z) 1 a p (n)z n n1
全极点(AP)滤波器
H(z) 1
1
A(z)
p
1 a p (n)z n
AZAP滤波器
n1
p
H (z) m0 bm z m B(z) A(z) A(z)
第29页/共35页
一、全零点(AZ)滤波器的格型结构
第30页/共35页
反射系数Kp的确定
根据系统函数,由高阶系数递推各低阶反射系数Kp
K p a p ( p)
a p1(i)
ap
(i)
K 1
pap
K
2 p
(
p
i)
(i 1,2,, p 1)
K
p1
a p1 (
p
1)
ap(
p
1) K
1
K
2 p
p a p (1)
第31页/共35页
二、全极点(AP)滤波器的格型结构
e
f p
[k
第4章 数字滤波器的算法结构
第4章思考题参考解答1、答:直接型结构是由差分方程或系统函数直接得到,容易获得该结构是其优点。
这种结构最大的缺点是对参数的量化非常敏感。
当N较大时,由于参数量化会导致系统的零极点位置有较大的位移,如果极点移到单位圆上或单位圆外,将导致系统不稳定。
因此,经常使用的是一阶、二阶IIR直接结构。
2、答:级联型结构是由一阶系统或若干二阶系统级联而成的系统。
结构中的每个一阶系统决定一个零点和一个极点,每个二阶系统决定一对零点和一对极点,通过分别调节零点系数和极点系数,可以分别改变一对零点或一对极点的位置。
相对于直接型结构,调整方便是级联型的优点。
此外,级联结构中后面网络输出不会再流入到前面,运算误差的积累相对直接型较小。
3、答:并联型结构通常由一个一阶网络与若干分子阶次较分目的阶次低的二阶网络并联而成的系统。
结构中的每个一阶网络决定一个实数极点,每个二阶网络决定一对共轭极点,因此调整极点位置方便,但零点位置调整不如级联型方便。
这种结构中,由于基本网络是并联的,产生的运算误差互不影响,相对于直接型和级联型,并联型运算误差最小,并且由于基本网络并联,可同时对输入信号进行运算,相对于直接型和级联型,并联型运算速度最高。
4、答:并联型结构的FIR滤波器优点有两个:一是其在频率采样点ωk处H(e jω)=H(k)为乘法器的系数,调整该系数就能直接调整频率特性;二是对于任意滤波器,只要h(n)的长度N相等,除乘法器的系数H(k)不同外,其他结构包括梳状滤波器和并联一阶或二阶网络结构均一样,便于实现标准化、模块化。
其缺点是当N较大时,二阶网络很多,使乘法器和延时器很多,造成结构复杂。
5、略。
参考教材6.1.1节。
6、略。
参考教材189页。
7、答:“转置定理”适用于所有介绍过的网络的转置变换,而“求逆准则”只适用于零、极点格型结构的相互转换。
8.答:数字谐振器的零点有两种放置法,可以放置在原点,或将两个零点放置在z=±1处。
实验四 数字滤波器的结构
制直接型和并联型的信号流图。
六、实验报告要求
简述实验目的和实验原理。
列写练习题的代码并绘制程序产生的图
形。
总结实验中你的收获和体会。
a=[1];
[sos,g]=tf2sos(b,a)
程序运行结果:
sos =1.0000 g =2 级联型的表达式
0.9500
0 1.2500
1.0000 1.0000
0 0
0 0
1.0000 -0.5000
H ( z ) 2(1 0.95 z 1 )(1 0.5 z 1 1.25 z 2 )
一个离散LSI系统可以用系统函数表示:
Y ( z ) b( z ) H ( z) X ( z ) a( z )
m0 N
b
M
m
z m
1 ak z k
k 1
b0 b1 z 1 b2 z 2 ... bm z m 1 a1 z 1 a2 z 2 ... ak z k
实验四 数字滤波器的结构
一、实验目的
加深对数字滤波器分类与结构的了解。 掌握数字滤波器的基本结构及其相互间
的转换方法。
学习利用MATLAB语言进行数字滤波器各
种结构相互间转换。
二、实验原理及方法
数字滤波器的分类
数字滤波器从滤波功能上可以分为低通、高通、带通、
带阻以及全通滤波器;根据系统的单位冲激响应的特 性,又可以分为有限长(FIR)和无限长(IIR)冲激 响应滤波器。
也可用差分方程表示:
y (n) ak y (n k ) bm x(n m)
k 1 m 0 N M
以上两个公式中,当ak至少有一个不为0时,则在有限
5.1 数字滤波器结构的表示方法
b z−k ∑k a 1−∑ k z−k
k= 1 k=0 N
M
——电子信息工程 电子信息工程
三、基本运算单元及两种方法表示
方框图表示 单位延时 信号流图法
x(n)
乘常数
a
ax(n)
x(n)
a
ax(n)
x1(n)
相加
+
x1(n)+x2(n) x1(n)
x1(n)+x2(n)
x2(n)
x2(n)
——电子信息工程 电子信息工程
流图结构: 流图结构:
x(n) b0 1 5 a1 a2
2 z 3 -1 z 4
-1
y(n)
输出结点 或阱结点
输入结点 或源结点
结点1、 结点1、5 相当于相 加器 五、画信号流图的规则
2、3、4相当于分 、 、 相当于分 支结点(一个输入, 支结点(一个输入, 一个或多个输出) 一个或多个输出)
节点的值=所有输入支路的值之和 节点的值 所有输入支路的值之和 支路的值=支路起点处的节点值 支路的值 支路起点处的节点值 ×传输系数
四、数字滤波器结构的表示方法: 数字滤波器结构的表示方法: 方框图法、 方框图法、信号流图法
例:二阶数字滤波器
y(n) = a y(n −1 + a2 y(n −2) +b x(n) ) 1 0
方框图结构:x(n) 方框图结构:
b0
+
z-1
y(n)
+
a2
a1 y(n-1) z-1
——电子信息工程 电子信息工程
一、数字滤波器的功能
把输入序列变换(运算)成为所要求的输出序列。 把输入序列变换(运算)成为所要求的输出序列。
数字滤波器的基本结构IV
k1
k0
17
特点:
第一个网络实现零点,即实现x(n)加权延时:
N
bkx(n k)
k0
第二个网络实现极点,即实现y(n)加权延时:
N
ak y(n k)
k 1
可见,第二网络是输出延时,即反馈网络。
*共需(M+N)个存储延时单元。
18
二、直接II(典范)型结构
直接型结构是由两个网络级联组成:
52
这个设备是由输入输出 延时部分、系数ai、bi 存储器、运算器及控制 器组成。
每一部分都可以用数字 硬件来构成。
6
5.1 数字滤波器的基本概念
数字滤波器的描述 数字滤波器的分类
7
5.1 数字滤波器结构的表示方法
一个数字滤波器可以用差分方程来描述:
N
M
y(n)aky(nk)bkx(nk)
46
N为奇数时
N1
H(z) h(n)zn n0 N n 2 1 0 1h(n)znh N 2 1 zN 2 1nN N 1 11h(n)zn 2 令 nN1m N n 2 1 0 1h(n) znz(N 1 n) h N 2 1 zN 2 1
47
h(n)偶对称,取“+”
Hz=
m0 N
1 anzn
Yz X z
n1
将系统函数整理为:
Hz=11 .5 02 .3 .1 zz1 1 00 .2 .4 zz2 2
1.52.1z10.4z2 1 0.3z10.2z2
23
H(z)1.52.1z10.4z2 10.3z10.2z2 得 a1 0.3,a2 0.2 b0 1.5 b1 2.1,b2 0.4 直接I型结构:
数字滤波器的结构
1)
数字滤波器实现的方法有两种
硬件实现——利用专用数字信号处理器或通用数字信号处理器来实
现,通常称为DSP芯片,这类芯片是解决实时处理要求的单片可编程
处理器。
2)
计算机软件实现——利用计算机和通用软件编程实现。把滤波器要 完成的运算编成程序,通过计算机来执行,称为计算机软件实现。 无论是利用硬件实现还是利用计算机软件实现,由于使用的设 备不会是无限精度的。因此用计算机处理数字信号的过程中,都是 用有限字长的算术运算来逼近无限精度的运算,所以运算过程中必 然会引入种种误差,这些误差来源主要有三个方面:
种情况:或者是实根,或者是共轭复根。如果是共轭复根,
则零极点必须共轭成对,也就是说 ci 必有 ci 同时存在。为 了结构上的一致性,可用实系数二阶多项式表示H(Z) 。
H ( z) A
(1
i 1 M i 1
M
1i
z 1 2i z 2 ) z 1 2i z 2 )
H ( z)
1 bi Z i
i 1 N
有影响,分母系数对极点有影响。
要调整滤波器的零、极点,只能间接调整 ai 、bi 的系数, 因为系数 ai 、bi 与系统函数零、极点的关系不是直接的而是 间接关系,所以系数对滤波器的性能控制不明显。
②系统的变化过于灵敏(是指对有限精度运算过于灵敏) 因为系数 ai 、bi 和所有的零、极点有关,所以对 ai 、bi 的精度要求比较高。当阶数越高时,要求的精度就越高 。由于有限字长效应的影响,造成这种结构容易产生较 大误差,甚至系统会出现不稳定的情况。
1、采样信号的量化误差:
模拟信号在数字化过程中,经采样后的采样信号的幅度是一个连续 量,需要把它量化为有限个电平值。如下图:
《第七章 数字滤波器的结构》
2、由基本信号流图求系统函数H(z) 由基本信号流图求系统函数H(z) 对给定的信号流图,设置中间节点变量; 对给定的信号流图,设置中间节点变量; 节点变量w(n)等于该节点的所有输入支路变量之和; 节点变量w(n)等于该节点的所有输入支路变量之和; w(n)等于该节点的所有输入支路变量之和 确定流图的输入与输出关系,求出系统函数H(z)。 确定流图的输入与输出关系,求出系统函数H(z)。 H(z) 已知基本信号流图如下,求其系统函数H(z) [例]:已知基本信号流图如下,求其系统函数H(z) b0
第七章 数字滤波器的结构 本章主要内容
数字滤波器结构的表示方法; 数字滤波器结构的表示方法; 无限长脉冲响应(IIR)数字滤波器的基本结构 无限长脉冲响应( ) 有限长脉冲响应( 有限长脉冲响应(FIR)数字滤波器的基本结构 )
数字滤波器的表示方法 差分方程:描述时域离散系统输入输出之间的关系。 差分方程:描述时域离散系统输入输出之间的关系。 单位脉冲响应:系统对δ 的零状态响应 的零状态响应。 单位脉冲响应:系统对δ(n)的零状态响应。 系统函数:系统单位脉冲响应 系统函数:系统单位脉冲响应h(n)的Z变换 的 变换 如果系统输入输出服从N阶差分方程, 如果系统输入输出服从 阶差分方程,即: 阶差分方程
对于某一确定的差分方程,可对应多种不同的算法,例如: 对于某一确定的差分方程,可对应多种不同的算法,例如:
1 H1(z) = 1− 0.8z−1 + 0.15z−2 −1.5 2.5 H2 (z) = + −1 1− 0.3z 1− 0.5z−1 1 1 H3(z) = ⋅ −1 1− 0.3z 1− 0.5z−1
写出差分方程如下: 写出差分方程如下 解:由H(z)写出差分方程如下:
数字信号处理-第五章数字滤波器的基本结构(new)
H ( z) A
将两个一阶因子组合成二阶因子,则
数字信号处理-第五章 数字滤波器络结构及 FIR数字滤波器的基本网络结构
数字信号处理-第五章 数字滤波器的基本结构
滤波器表示方式
(1)系统函数
k b z k M
Y ( z) H ( z) X ( z)
1 ak z k
k 1
k 0 N
1 ak z k
k 1
k 0 N
N2 M N Ak Bk (1 g k z 1 ) k G z k 1 1 * 1 1 c z ( 1 d z )( 1 d z ) k 1 k 1 k 0 k k k N1
一般IIR滤波器满足
N1
数字信号处理-第五章 数字滤波器的基本结构
5.2 无限长单位冲激响应(IIR)滤波器的基本结构)
IIR滤波器有以下几个特点: (1)系统的单位冲激响应 (2)系统函数
h( n)
是无限长的
H ( z)
在有限z平面(
0 z
)上有极点存在
(3)结构上存在输出到输入的反馈,也就是结构是递归的 1、直接Ⅰ型 一个IIR滤波器的有理系统函数为:
x n
3 1.5 -1.5 0.5
z 1 z 1 z 1
-3.5 2.5
y n
数字信号处理-第五章 数字滤波器的基本结构 级联型:
3z 3 3.5z 2 2.5z 3 3.5z 1 2.5z 2 1 H ( z) 2 2 z z 1 z 0.5 1 z z 1 0.5z 1
第九章数字滤波器的分类及结构
画出该滤波器的直接型结构。
解答:如右图所示。 直接型结构的特点:
x ( n) 5/4 -3/4 1/8 z-1 z-1 8 -4 11
y(n)
所需要的延迟单元最少;
系统调整不方便; 受有限字长影响较大。
z-1
-2
3. IIR 滤波器的结构
对系统函数 H(z) 进行因式分解:
H ( z)
M
1.
2. 3. 4. 5.
数字滤波器的分类
数字滤波器结构的表示方法 IIR 滤波器的结构 FIR 滤波器的结构 离散时间系统 的 Lattice 结构
1. 滤波器的分类
根据单位冲激响应 h(n) 的时间特性分类:
无限冲激响应数字滤波器(IIR) 有限冲激响应数字滤波器(FIR)
根据实现方法和形式分类:
r 0
M
3. IIR 滤波器的结构
表示为两个系统级联的形式:
x ( n)
H 1 ( z) b0 z-1
H 1 ( z ) b( r ) z
r 0 M r
y′ (n)
H 2 ( z)
y(n)
x ( n)
y′(n)
-a1 -a2 z-1 z- 1
y(n)
1 1 a(k ) z k
k 1 N
b1
z-1
b2
bM-1
H 2 ( z)
z-1
bM
-aN-1 -aN z-1
称为直接Ⅰ型结构。
3. IIR 滤波器的结构
直接Ⅰ型的变型:
x ( n)
H 2 ( z) y′ (n) H 1 ( z) b0 b1 b2 bM-1 -aN-1
H 1 ( z ) b( r ) z r
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
二、数字滤波器的工作原理
设x(n)是系统的输入,X (e )是其付氏变换。 y (n)是系统的输出,Y (e jw )是其付氏变换。 则:
x(n) h(n) y(n)
jw
则LTI系统的输出为:
y ( n)
看出:输入序列的频谱X (e jw )经过滤波器 jw jw jw (其系统性能用H (e )表示)后变成X (e ) H (e ) jw jw jw 选取H (e ), 使滤波器输出X (e ) H (e )符合我们的要求, 这就是数字滤波器的工作原理。
m
1 jw jw h ( n m ) x ( m ) F [ X ( e ) H ( e )]
三、数字滤波器表示方法
• 有两 种表示方法:方框图表示法;流图 表示法. • 数字滤波器中,信号只有延时,乘以常数 和相加三种运算。 • 所以DF结构中有三个基本运算单元:加 法器,单位延时,乘常数的乘法器。
六、本章介绍主要的内容
1.分别介绍FIR、IIR滤波器实现的基本结构。 2.介绍一种特殊的滤波器结构实现形式:格型 滤波器结构.
第二节 IIR DF的基本结构
一、IIR DF特点
1.单位冲激响应h(n)是无限长的n→∞
2.系统函数H(z)在有限长Z平面(0<|Z|<∞)有 极点存在。 3.结构上存在输出到输入的反馈,也即结构 上是递归型的。 4.因果稳定的IIR滤波器其全部极点一定在单 位园内。
第 6章
数字滤波器的结构 DF (Digital Filter)
第一节 引言
一、什么是数字滤波器
• 顾名思义:其作用是对输入信号起到滤波 的作用;即DF是由差分方程描述的一类特 殊的离散时间系统。 • 它的功能:把输入序列通过一定的运算变 换成输出序列。不同的运算处理方法决定 了滤波器的实现结构的不同。
4.模拟滤波器的理想幅频特性
H ( j)
LPAF
c
H ( j)
c
HPAF
c c
H ( j) c
BPAF
H ( j)
c
BSAF
c 2 c1 c1 c 2
5.数字滤波器的理想幅频特性
H (e jw )
…….
H (e )
jw
2
c
…….
1、经典滤波器
• 假定输入信号x(n)中的有用成分和希望去除 的成分,各自占有不同的频带。当x(n)经 过一个线性系统(即滤波器)后即可将欲去 除的成分有效地去除。但如果信号和噪声的 频谱相互重叠,那么经典滤波器将无能为力。
|X(ejw)| 有用 无用 |H(ejw)| |Y(ejw)|
wc
w
wc
LPDF
H (e )
jw
2
3
…….
HPDF
H (e )
jw
2
…….
BPDF
2
BSDF
五、研究DF实现结构意义
1.滤波器的基本特性(如有限长冲激响应FIR与无 限长冲激响应IIR)决定了结构上有不同的特点。
2.不同结构所需的存储单元及乘法次数不同,前 者影响复杂性,后者影响运算速度。 3.有限精度(有限字长)实现情况下,不同运算 结构的误差及稳定性不同。 4.好的滤波器结构应该易于控制滤波器性能,适 合于模块化实现,便于时分复用。
1、方框图、流图表示法
方框图表示法: 单位延时
Z-1
信号流图表示法:
Z-1
系数乘
a
a
相加
把上述三个基本单元互联,可构成不同数字网络或运 算结构,也有方框图表示法和流图表示法。
2.例子
例:二阶数字滤波器:
y(n) a1 y(n 1) a2 y(n 2) b0 x(n)
其方框图及流图结构如下:
x(n) b0
a1 Z-1
y(n)
x(n)
b0
a1 Z-1 Z-1
y(n)
a2
Z-1
a2
看出:可通过流图或方框图看出系统的运算步骤和运算结构。
以后我们用流图来分析数字滤波器结构。DF网络结构或DF运 算结构二个术语有微小的差别,但大抵一样,可以混用。
四、数字滤波器的分类
• • • • 滤波器的种类很多,分类方法也不同。 1.从功能上分;低、带、高、带阻。 2.从实现方法上分:FIR、IIR 3.从设计方法上来分:Chebyshev(切比雪 夫),Butterworth(巴特沃斯) • 4.从处理信号分:经典滤波器、现代滤波器 • 等等。
wc
w
2.现代滤波器
它主要研究内容是从含有噪声的数据记录(又称 时间序列)中估计出信号的某些特征或信号本身。一 旦信号被估计出,那么估计出的信号将比原信号会有 高的信噪比。 现代滤波器把信号和噪声都视为随机信号,利用 它们的统计特征(如自相关函数、功率谱等)导出一 套最佳估值算法,然后用硬件或软件予以实现。
以下我们讨论M<=N情况。
则这一系统差分方程为:
y (n) ai y (n i ) bi x(n i )
i 0 i 0
N
M
2、直接I型 (1)直接I型流图
• IIRDF的差分方程就代表了一种最直接的计算公式, 用流图表现出来的实现结构即为直接I型结构(即由 差分方程直接实现。) 方程看出:y(n)由两部分组成:
二、IIR DF基本结构
IIR DF类型有:直接型、级联型、并联型。 直接型结构:直接I型、直接II型(正准型、 典范型)。
1、 IIR DF系统函数及:
i b Z i M
H(z )
1
i 0 N
i a Z i i 1
Y ( z) X ( z)
现代滤波器理论源于维纳在40年代及其以后的工 作,这一类滤波器的代表为:维纳滤波器,此外,还 有卡尔曼滤波器、线性预测器、自适应滤波器。 本课程主要讲经典滤波器,外带一点自适应滤波器
3.模拟滤波器和数字滤波器
• 经典滤波器从功能上分又可分为: 低通滤波器(LPAF/LPDF):Low pass analog filter 带通滤波器(BPAF/BPDF):Bandpass analog filter 高通滤波器(HPAF/HPDF):High pass analog filter 带阻滤波器(BSAF/BSDF):Bandstop analog filter • 即它们每一种又可分为:数字(Digital)和模拟 (Analog)滤波器。