数字信号处理第五章(经典)
合集下载
第五章 数字信号处理课后答案刘顺兰版
H α ( s ) = H αN (
ω c = 2πf c T = 2π × 400 HZ / 6000 HZ =
Ωc = 2 ωc 2 π 2 tg = tg ( ) = 0.2 × T T 2 T 15
2π 15
s=
2 1 − z −1 , T 1 + z −1
s = Ωc
1 − z −1 −1 π tg ( ) 1 + z 15 1 1
=
1 + 3z −1 + 3z −2 + z −3 0.005376(1 + 3z −1 + 3z −2 + z −3 ) = 186 − 412 z −1 + 318 z − 2 − 84 z −3 1 − 2.215 z −1 + 1.71z − 2 − 0.4516 z −3
5.24 用双线性变换设计一个三阶巴特沃思数字高通滤波器，采样频率为 f s = 6 KHZ ，截止 频率为 f c = 1.5 KHZ （不计 3KHZ 以上的频率分量） 。 解法 1：三阶巴特沃思低通模拟滤波器的原型函数：
按照冲激不变条件，可以写出
因此系统函数为
H ( z ) = ∑ h(n) z − n
n =0
∞
1 1 2 2 = + 1 − e − aT e − jbT z −1 1 − e −aT e jbT z −1 = 1 − (e − aT cos bT ) z −1 (1 − e − aT e − jbT z −1 )(1 − e −aT e jbT z −1 )
所以
ω1 + ω 2
H BP ( z ) = H αN ( s )
s=
1 1+ z − 2 3 1− z − 2
数字信号处理-第五章
系 统 函 数
：
H (z) n M 0h (n )z n Y X ( (z z ) ) b 0 b 1 z 1 b 2 z 2 b M z M
单位脉冲响应的值等于差分方程系数：
h
h(n)=bn
n=0,1,·····,M
33
FIR数字滤波器的特点：
系统函数：
N1
H(z) h(n)zn n0
有N-1个零点分布于z平面 z=0处 是N-1阶极点
h
26
还可以如下式这样进行分解： H (z)1 1 0 0..6 4z z 1 11 1 0 0..5 3z z 1 1H 3(z)H 4(z)
h
27
级联型结构的特点：
调整某一路的分子系数能单独调整滤波器的一组 零点，而不影响其它零极点； 调整某一路的分母系数能单独调整滤波器的一组 极点，而不影响其它零极点；便于调整滤波器频 率响应性能
直接型
h
38
将H(z)进行因式分解，得到： H(z)=(0.6+0.5z-1)(1.6+2z-1+3z-2) 按照上式画出它的级联型结构如图所示。
级联型
h
39
5.5 线性相位网络结构
FIR滤波器单位抽样响应h(n)为实数，0 n N 1 且满足：
偶对称： h (n ) h (N 1 n ) 或奇对称 h (n ) h (N 1 n ) ： 即对称中心在 (N-1) / 2处 则这种FIR滤波器具有严格线性相位。
bi zi
i0 N 1 ak zk
k 1
基本运算：加法，乘法（乘以常数），移位（时延）
h
3
信号流图由基本支路构成
1.基本支路箭头表示信号流向，两个圆点表示输入输出节点，箭头旁边的 符号表示增益（缺省为1）
第五章 时域离散系统的基本网络结构
数字滤波器和FFT一样，是数字信号处理 的重要内容。
本章的主要内容就是描述数字滤波器的基 本网络结构。（IIR、FIR）
引言
时域离散系统或网络可以用差分方程、单 位脉冲响应以及系统函数进行描述。
M
N
y(n) bi x(n i) ai y(n i)
i0
i 1
系统函数H(z)为
M
H (z)
(2) 流图环路中必须存在延时支路；
(3) 节点和支路的数目是有限的。
信号流图表达的系统含义
每个节点连接的有输入支路和输出支路，节点变 量等于所有输入支路的输出之和.
根据信号流图可以求出系统函数（节点法、梅逊 公式法）。
1(n) 2 (n 1) 2 (n) 2 (n 1) 2 (n) x(n) a12 (n) a21n y(n) b21(n) b12 (n) b02(n)
画出H(z)的直接型结构和级联型结构。
级联型
解： 将H(z)进行因式分解，得到： H(z)=(0.6+0.5z-1)(1.6+2z-1+3z-2)
其直接型结构和级联型结构如图所示。
x(n)
0.6
z－ 1 0.5
1.6 z－ 1
2 z－ 1
3
y(n) x(n)
z－ 1
z－ 1
z－ 1
0.96 2
2.8 1.5 y(n)
0 j
y(n)
1 j
z－ 1 1j
1 j
z－ 11 j
(a)
2 j
z－
1
2
j
(b)
一阶和二阶直接型网络结构 (a)直接型一阶网络结构；(b)直接型二阶网络结构
IIR的级联型例题
本章的主要内容就是描述数字滤波器的基 本网络结构。（IIR、FIR）
引言
时域离散系统或网络可以用差分方程、单 位脉冲响应以及系统函数进行描述。
M
N
y(n) bi x(n i) ai y(n i)
i0
i 1
系统函数H(z)为
M
H (z)
(2) 流图环路中必须存在延时支路；
(3) 节点和支路的数目是有限的。
信号流图表达的系统含义
每个节点连接的有输入支路和输出支路，节点变 量等于所有输入支路的输出之和.
根据信号流图可以求出系统函数（节点法、梅逊 公式法）。
1(n) 2 (n 1) 2 (n) 2 (n 1) 2 (n) x(n) a12 (n) a21n y(n) b21(n) b12 (n) b02(n)
画出H(z)的直接型结构和级联型结构。
级联型
解： 将H(z)进行因式分解，得到： H(z)=(0.6+0.5z-1)(1.6+2z-1+3z-2)
其直接型结构和级联型结构如图所示。
x(n)
0.6
z－ 1 0.5
1.6 z－ 1
2 z－ 1
3
y(n) x(n)
z－ 1
z－ 1
z－ 1
0.96 2
2.8 1.5 y(n)
0 j
y(n)
1 j
z－ 1 1j
1 j
z－ 11 j
(a)
2 j
z－
1
2
j
(b)
一阶和二阶直接型网络结构 (a)直接型一阶网络结构；(b)直接型二阶网络结构
IIR的级联型例题
第五章 数字信号处理- 微弱信号处理
第五章微弱信号处理5.1 微弱信号检测技术中气体浓度检测仪中的应用微弱信号不仅意味着信号的幅度小，而且主要指被噪声淹没中的信号。
为了检测被背景噪音淹没的信号，就需要分析噪声产生的原因和规律，研究被测信号的特点、相关性以及噪声的统计特性，以寻找出从背景噪声中检测有用信号的方法。
因此，微弱信号检测技术的首要任务是提高信噪比。
它不同于一般的检测技术，它注重的不是传感器的物理模型和传感原理、相应的信号转换电路和仪表实现方法，而是如何抑制噪声和提高信噪比。
由于被测量的信号微弱，传感器的固有噪声、放大电路及测量仪器的固有噪声以及外界的干扰噪声往往比有用信号的幅度大得多，放大被测信号的过程同时也放大了噪声，而且必然还会附加一些额外的噪声，因此只靠放大是不能把微弱信号检测出来的。
只有在有效地抑制噪声的条件下增大微弱信号的幅度，才能提取出有用的信号。
为了表征噪声对信号的淹没程度，引入信噪比SNR来表示，它指的是信号的有效值S与噪音的有效值N之比。
而评价一种微弱信号检测方法的优劣，经常采用两种指标：一种是信噪改善比SNIR，它等于系统输出端的信噪比SNR和系统输入段oSNR之比。
SNIR越大，表明系统抑制噪声的能力越强。
i另一个指标是检测分辨率，指的是检测仪器指示值可以响应与分辨的最小输入值的变化值。
检测分辨率不同于检测灵敏度，后者表示的是检测系统标定曲线的斜率，定义为输出变化量y∆之比。
一般情况下，∆的输入变化量x∆与引起y灵敏度越高，分辨率越好。
但提高系统的放大倍数虽可提高灵敏度，但却不一定能提高分辨率，因为分辨率要受噪声和误差额制约。
5.1.1 本检测系统的噪声源广义的噪声是扣除被测信号真实值以后的各种测量值，可以分为两类：一是干扰；另一被称为电子噪声（狭义）。
干扰是指被非被测信号或非测量系统所引起的噪声。
从理论上讲，干扰是属于理想上可排除的噪声。
不少干扰源发出的干扰是有规律的，有些具有周期性，有些只是瞬时值。
《信号、系统与数字信号处理》第五章 Z变换与离散系统的频域分析
同理
sinh0nun
1 2
e0n
e0n
un
1 z
2
z
e0
z z e0
z2
z sinh0 2z cosh0
1
z max e0 , e0
2、双边z变换的移位 n0 0
若 xn X z
RX
z
R X
则 x n n0 z n0 X z
RX
z
R X
证明： Z x n n0
n
xT t nT estdt
n
xnT esnT
n
令 z esT 引入新的复变量， 将上式写为
X s s xnT zn
n
此式是复变量 z 的函数（T 是常数），记为
X z xnzn
n
x 2z2 x 1z x0 x1z1 x2z2
Z xn 2un z2 X z z1x1 x 2
3） 若 xn 为因果序列 xnun X z
则 xn mun zm X z
m0
xn
mun
zm
X
z
m1 k 0
xk
z
k
例5-9 求周期序列的单边z变换
解： 周期序列 xn xn rN
m0
令 n 0 ~ N 1 的主值区序列为 x1 n ，
( z 1)
4、指数序列加权
若 xn X z RX z RX
则 an xn X a1z
RX a 1z RX
证：Z an xn an xnzn
n
xn a1z n X z / a
n
RX a 1z RX
a
R X
z
a
R X
利用
数字信号处理 第五章
+ 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 +
数字信号处理第五章-IIR数字滤波器的设计
24
2、由模平方函数确定系统函数
模拟滤波器幅度响应常用幅度平方函数表示：
| H ( j) |2 H ( j)H *( j)
由于冲击响应h(t)为实函数，H ( j) H *( j)
| H ( j) |2 H ( j)H ( j) H (s)H (s) |s j
H (s)是模拟滤波器的系统函数，是s的有理分式；
分别对应：通带波纹和阻带衰减(阻带波纹)
（4种函数）
只介绍前两种
31
32
33
无论N多大，所 有特性曲线均通 过该点
特性曲线单调减小，N越大，减小越慢 阻
特性曲线单调减小，N越大，减小越快
34
20Nlog2：频率增加一倍，衰减6NdB
35
另外：
36
无论N多大，所 有特性曲线均通 过Ωc点: 衰减3dB, Ωc 为 3dB带宽
8
根据
(线性相位滤波器)
非线性相位滤波器
9
问题：
理想滤波器的幅度特性中，频带之间存 在突变，单位冲击响应是非因果的；
只能用逼近的方法来尽量接近实际的要 求。
滤波器的性能要求以频率响应的幅度特 性的允许误差来表征，如下图：
10
p
11
低通滤波器的频率响应包括：
通带：在通带内，以幅度响应的误差δp逼近 于1；
20
3、数字滤波器设计的基本方法
利用模拟理论进行设计 先按照给定的技术指标设计出模拟滤波 器的系统函数H(s)，然后经过一定的变 换得到数字滤波器的系统函数H(z),这实 际上是S平面到Z平面的映射过程： 从时域出发，脉冲响应不变法 从频域出发，双线性变换法 适合于设计幅度特性较规则的滤波器， 如低通、高通等。
由于系统稳定， H(s)的极点一定落在s的左半 平面，所以左半平面的极点一定属于H(s)，右 半平面的极点一定属于H(-s)。
2、由模平方函数确定系统函数
模拟滤波器幅度响应常用幅度平方函数表示：
| H ( j) |2 H ( j)H *( j)
由于冲击响应h(t)为实函数，H ( j) H *( j)
| H ( j) |2 H ( j)H ( j) H (s)H (s) |s j
H (s)是模拟滤波器的系统函数，是s的有理分式；
分别对应：通带波纹和阻带衰减(阻带波纹)
（4种函数）
只介绍前两种
31
32
33
无论N多大，所 有特性曲线均通 过该点
特性曲线单调减小，N越大，减小越慢 阻
特性曲线单调减小，N越大，减小越快
34
20Nlog2：频率增加一倍，衰减6NdB
35
另外：
36
无论N多大，所 有特性曲线均通 过Ωc点: 衰减3dB, Ωc 为 3dB带宽
8
根据
(线性相位滤波器)
非线性相位滤波器
9
问题：
理想滤波器的幅度特性中，频带之间存 在突变，单位冲击响应是非因果的；
只能用逼近的方法来尽量接近实际的要 求。
滤波器的性能要求以频率响应的幅度特 性的允许误差来表征，如下图：
10
p
11
低通滤波器的频率响应包括：
通带：在通带内，以幅度响应的误差δp逼近 于1；
20
3、数字滤波器设计的基本方法
利用模拟理论进行设计 先按照给定的技术指标设计出模拟滤波 器的系统函数H(s)，然后经过一定的变 换得到数字滤波器的系统函数H(z),这实 际上是S平面到Z平面的映射过程： 从时域出发，脉冲响应不变法 从频域出发，双线性变换法 适合于设计幅度特性较规则的滤波器， 如低通、高通等。
由于系统稳定， H(s)的极点一定落在s的左半 平面，所以左半平面的极点一定属于H(s)，右 半平面的极点一定属于H(-s)。
数字信号处理第三版(姚天任、江太辉) 答案 第五章-unprotected
E[x(i)x( j)] −
2
N −1 N −1
E[x(i)x( j)]
N n=0
N2 i=0 j=0
N2 i=0 j=0
∑ ∑∑ =
1 N
N −1
E[x2 (n)] −
n=0
1 N2
N −1 N −1
E[x(i)x( j)]
i=0 j=0
∑ ∑ ∑∑ =
1
N −1
E[x2 (n)] −
1
N −1
N −1 N −1
∫ = 1
q
0 −q
xdx
=
1 2q
x2
|0−q =
−
q 2
∞
∫ mx2 = E[x2 ] = −∞ xpx2 (x)dx
∫ = 1
q
q/2 −q/2
xdx
=
1 2q
x2
|−q
/2 q/
2
=
0
∞
∫ mx3 = E[x3 ] = −∞ xpx3 (x)dx
∫ = 1
2π
2π 0
xdx =
1 4π
x2
|02π = π
∞ −∞
(x
−
mx2
)2
px2
( x)dx
∫ = 1 q
q/2 −q / 2
x2dx
=
1 3q
x3
|q / 2
−q/
2
=
q2 12
∫ σ 2 x3
=
E[( x3
− mx3 )2 ] =
∞ −∞
(x
−
mx3
)2
px3
( x)dx
∫ = 1
2π
2π 0
2
N −1 N −1
E[x(i)x( j)]
N n=0
N2 i=0 j=0
N2 i=0 j=0
∑ ∑∑ =
1 N
N −1
E[x2 (n)] −
n=0
1 N2
N −1 N −1
E[x(i)x( j)]
i=0 j=0
∑ ∑ ∑∑ =
1
N −1
E[x2 (n)] −
1
N −1
N −1 N −1
∫ = 1
q
0 −q
xdx
=
1 2q
x2
|0−q =
−
q 2
∞
∫ mx2 = E[x2 ] = −∞ xpx2 (x)dx
∫ = 1
q
q/2 −q/2
xdx
=
1 2q
x2
|−q
/2 q/
2
=
0
∞
∫ mx3 = E[x3 ] = −∞ xpx3 (x)dx
∫ = 1
2π
2π 0
xdx =
1 4π
x2
|02π = π
∞ −∞
(x
−
mx2
)2
px2
( x)dx
∫ = 1 q
q/2 −q / 2
x2dx
=
1 3q
x3
|q / 2
−q/
2
=
q2 12
∫ σ 2 x3
=
E[( x3
− mx3 )2 ] =
∞ −∞
(x
−
mx3
)2
px3
( x)dx
∫ = 1
2π
2π 0
精品文档-数字信号处理(吴瑛)-第5章
第5章 数字滤波器概论
5.3 实际滤波器的设计指标
5.3.1 图5.3.1是理想低通滤波器的幅频响应，该理想低通滤波
器具有截止频率ωd。可以看出，理想滤波器在通带内幅度为常 数(非零)，在阻带内幅度为零。另外，一般理想滤波器 要求具有线性相位(在第8章讨论)，这里假设相频响应 θ(ω)=0
h(n) sin(nd )
第5章 数字滤波器概论
1. 根据H(ejω) 一般数字滤波器从滤波功能上分类，和模拟滤波器一样， 可以分成低通、高通、带通和带阻等滤波器。它们的理想幅频 响应如图5.2.2
第5章 数字滤波器概论
图5.2.2 (a) 低通; (b) 高通; (c) 带通; (d) 带阻
第5章 数字滤波器概论
需要注意的是，数字滤波器的频率响应H(ejω)都是以2π 为周期的，滤波器的低通频带处于2π的整数倍处，而高通频 带处于π
5.3.2 当滤波器形状为非理想时，要用一些参数指标来描述其关
键特性。图5.3.5表示低通滤波器的幅频响应。滤波器的通带 定义了滤波器允许通过的频率范围。在阻带内，滤波器对 信号严重衰减。ωp和ωs分别称为通带截止频率(或通带上限频 率)和阻带截止频率(或阻带下限频率)。参数δ1定义了通带波 纹(Pass Band Ripple)，即滤波器通带内偏离单位增 益的最大值。参数δ2定义了阻带波纹(Stop Band Ripple)，即 滤波器阻带内偏离零增益的最大值。
截短脉冲响应自然会对频率响应产生影响。截短后，滤波 器幅频响应曲线不再是理想矩形，通带不再平坦，有过渡带， 同时阻带衰减不再为零。图5.3.4给出了因果脉冲响应 的幅频响应。当然，脉冲响应保留的采样点越多，即滤波器阶
第5章 数字滤波器概论 图5.3.4 非理想低通滤波器因果脉冲响应的幅频响应
数字信号处理 第五章
G0
G0
x(n )
•
• • •
•
11
•
•
z
01
1
11
• •
•
y (n )
• •
或
x (n )
21
•
z 1
•
• 01 • • 11 z 11 • • • z y (n ) • 21 • •
1 1
并联型网络结构的优缺点： 优点：运算速度快，调整极点方便，乘法运算量化误差在输出端的 噪声功率最小。 缺点：调整零点不方便，当H(z)有多阶极点时，部分分式展开较麻烦。
b z
k 0 N k k 1
M
k
1 ak z k
w 令： 1 ( n) bk x( n k )
k 0
M
有： y( n) w1 ( n) ak y(n k )
k 1
N
直接Ⅰ型
x (n)
z
• b1 z 1 b2 •
z 1 •
1
•
b0
w1 ( n)
• • • •
1
• •
1
• • •
y(n)
z 1
a1
a2
•
•
z 1
bM
aN
直接Ⅰ型 之1
•
z 1
x (n)
1
w2 (n)
• • •
a1
a2
• 1 z •
z
1
z 1
• 1 z • b1 z 1 • b2
z 1
w2 (n) b0
• • •
•
y(n)
x (n )
• • a1 z 1 • • 1 a z • 2 •
G0
x(n )
•
• • •
•
11
•
•
z
01
1
11
• •
•
y (n )
• •
或
x (n )
21
•
z 1
•
• 01 • • 11 z 11 • • • z y (n ) • 21 • •
1 1
并联型网络结构的优缺点： 优点：运算速度快，调整极点方便，乘法运算量化误差在输出端的 噪声功率最小。 缺点：调整零点不方便，当H(z)有多阶极点时，部分分式展开较麻烦。
b z
k 0 N k k 1
M
k
1 ak z k
w 令： 1 ( n) bk x( n k )
k 0
M
有： y( n) w1 ( n) ak y(n k )
k 1
N
直接Ⅰ型
x (n)
z
• b1 z 1 b2 •
z 1 •
1
•
b0
w1 ( n)
• • • •
1
• •
1
• • •
y(n)
z 1
a1
a2
•
•
z 1
bM
aN
直接Ⅰ型 之1
•
z 1
x (n)
1
w2 (n)
• • •
a1
a2
• 1 z •
z
1
z 1
• 1 z • b1 z 1 • b2
z 1
w2 (n) b0
• • •
•
y(n)
x (n )
• • a1 z 1 • • 1 a z • 2 •
数字信号处理_吴镇扬_第二版_第五章习题答案
5.7 （1）由于h2(n)是h1(n)圆周移位的序列，根据DFT的 2π 性质有： −j 4k − jπ k
H 2 (k ) = e
8
H 1 (k ) = e
H 1 (k )
~ ~ H1 ( k ) = H 2 ( k ) 成立 所以
(2)由于h1 (n ) 和h2 (n ) 均为偶对称序列，以其构成的低通滤波器
(3)若采用海明窗设计，则
⎡ ⎛ 2πn ⎞⎤ wHam ( n) = ⎢0.54 − 0.46 cos ⎜ ⎟ ⎥ RN ( n ) ⎝ N − 1 ⎠⎦ ⎣ 2 h( n) = sin[(n − α )ωc ]cos[(n − α )ω0 ]wHam (n) N 为奇数时， (n − α )π
h( n N 为偶数时， ) =
0 −ωc
e − jωα e jω nd ω
可见h(n)关于(N-1)/2偶对称，即 h( n) = h( N − 1 − n)
(1)当 N 为奇数时，为第一类滤波器。 (2)当N为偶数时，为第二类滤波器
⎧hd ( n) h( n) = hd ( n) ⋅ R(n ) = ⎨ ⎩0 0 ≤ n ≤ N −1
解：由经验公式可知若 不小于 At 40dB , 则
β = 0.5842 At - 21)0.4 + 0.07886(At - 21) ≈ 3.3953 （ At − 8 40 − 8 N= = ≈ 22.28 2.286∆ω 2.286× 0.2π ωc + ωr ωc′ = = 0.2π 2 ′ ⎧ωc ′ ⎪ π Sa[ωc (n − α )] n ≠ α ′ 1 ωc − jωα jωn ⎪ hd (n) = ∫ ′ e e dω = ⎨ ′ 2π −ωc ωc ⎪ n =α ⎪ ⎩ π
数字信号处理 第五章
第五章离散时间信号的数字处理Q5.1运行程序P5.1，产生连续时间序号及其抽样形式，并显示它们。
clf;t = 0:0.0005:1;f = 13;xa = cos(2*pi*f*t);subplot(2,1,1)plot(t,xa);gridxlabel('时间, msec');ylabel('振幅');title('连续时间序号 x_{a}(t)');axis([0 1 -1.2 1.2])subplot(2,1,2);T = 0.1;n = 0:T:1;xs = cos(2*pi*f*n);k = 0:length(n)-1;stem(k,xs);grid;xlabel('时间 n');ylabel('振幅');title('离散事件序号 x[n]');axis([0 (length(n)-1) -1.2 1.2])Q5.2 正弦信号的频率是多少赫兹？抽样周期是多少秒？正弦信号的频率f=13Hz，抽样周期T=0.1s。
Q5.3 解释两个axis命令的效果。
给x,y轴标刻度。
Q5.4 以比在程序P5.1中列出的抽样周期低的两个抽样周期和高的两个抽样周期的四个其他值，运行程序P5.1.评论你的结果。
T=0.04s T=0.08sT=0.15s T=0.3s由上图可以发现：当取的T越小时，得到的图形越接近原图形。
Q5.5 通过将正弦信号的频率分别变为3HZ和7HZ，重做习题Q5.1。
相应的等效离散时间信号与习题Q5.1中产生的离散时间信号之间有差别么？若没有，为什么没有？f=3Hz f=7Hz由图可以看出，变换频率得到的两个图没有区别，因为他们的抽样周期一样。
Q5.6 运行程序P5.2，产生离散时间信号x[n]及其连续时间等效ya[t]，并显示它们。
clf;T = 0.1;f = 13;n = (0:T:1)';xs = cos(2*pi*f*n);t = linspace(-0.5,1.5,500)';ya = sinc((1/T)*t(:,ones(size(n))) - (1/T)*n(:,ones(size(t)))')*xs;plot(n,xs,'o',t,ya);grid;xlabel('时间, msec');ylabel('振幅');title('重构的连续时间序号 y_{a}(t)');axis([0 1 -1.2 1.2]);图1 图2Q5.7 在程序P5.2中，t的范围和时间增量的值是什么？在图中，t的范围是什么？改变t的范围，显示上述程序所计算的全范围ya[t]并再次运行程序P5.2,。
中国石油大学《数字信号处理》第五章 有限离散傅立叶变换与快速傅立叶变换
X ( f ) x(n)e
n 0
N 1
j2nf
第一节 离散傅立叶变换
对于一个周期信号，只需研究其在一个周期内的 变化过程和规律。
因此，对XΔ(f )只要分析[0,1/Δ]范围内的频谱即可。
频谱XΔ(f )在一个周期内[0,1/Δ]是连续的，即所 对应的采样点数是无穷多个。实际上在计算有限离 散序列的频谱时，只能计算有限个点上的值，因此 1 把XΔ(f )在区间[0,1/Δ]分成N等份，每份的间隔为 N 则得到有限离散频谱：
对于非周期连续时间信号x(t)，其频谱X(f )为一 个无限长的非周期连续频谱函数 。
x(t ) X ( f )
F
若用采样间隔Δ对信号x(t)进行离散采样，则得到离 散序列x(nΔ) ，其频谱XΔ(f )是周期连续，可以表示 为：以1/Δ为周期将X(f )向两边延拓而得到的。则有
m X( f ) X ( f ) m
第一节 离散傅立叶变换 二、离散傅立叶变换(DFT)
上面讨论的离散傅立叶级数中xp(n)和Xp(k)都是周期性 的序列，是无限长的。
1、有限长离散序列和周期离散序列的关系
有限长序列x(t)
周期延拓 T T N
周期性的序列xp(n)
离散采样
x(t ) x p (t ) x p (n)
第一节 离散傅立叶变换
数字信号处理
Digital Signal Processing（D S P ）
中国石油大学
数字信号处理
Digital Signal Processing（D S P ）
第五章 有限离散傅立叶变换与快速傅立叶变换
第五章 有限离散傅立叶变换与快速傅立叶变换
数字信号处理第五章
b1 w1 b2
y(n)
w1(n) w2 (n 1) w2 (n) w'2 (n 1) w'2 (n) x(n) a1w2 (n) a2 w1(n)
-a2 (a)
y (n) b2 w1(n) b1w2 (n) b0 w'2 (n)
W1 ( z ) W2 ( z ) z 1 W2 ( z ) W '2 ( z ) z 1
x(n)
z-1
h(0)
N 1
z-1
h(1)
z-1
h(2)
z-1
h(N-2) h(N-1)
y(n)
2、级联型结构
将H(z)系统函数进行因式分解，并将共轭成对的零点放在一 起，形成一个系数为实数的二阶形式，这样级联型网络结构 就是由一阶或二阶因子构成的级联结构，其中每一个因式都 用直接型实现。
H（z ) h(n) Z n ( 0i 1i z 1 2i z 2 )
§5.3 无限长脉冲响应基本网络结构
IIR网络的特点:信号流图中含有反馈支路，即存在 环路，其单位脉冲响应是无限长的。有三种结构： 直接型、级联型、并联型 1、直接型
y ( n)
M
i b z i M
H1(z)
H(z)=H1(z)H 2(z)
N
H (z)
1 ai z i
i 1
b0 x ( n) w’2 z-1 -a1 w2 z-1 b1 w1 b2 (b) y(n)
H(z)
-a2 (a)
基本信号流图
非基本信号流图
根据信号流图可以求出网络的系统函数，方法是列出各个节 点变量方程，求解该方程，推导出输出与输入之间的关系。
程佩青_数字信号处理_经典版(第四版)_第5章_5.1-5.3
N阶FIR 数字滤波器
N 1
N 1
H (z) h(k)zk bk zk
k 0
k 0
M
yn bk xn k k 0
h(n)* x(n) (5.3.1)
M
N
特点： yn bk xn k ak yn k
k 0
k 1
（1）系统单位抽样响应h(n)是有限长的。
(5.1.1)
（2）必须全部的ak = 0（k=1,2,…,N），即结构上没有输出 到输入的反馈（没有递归结构，只有非递归结构）。
5.2 IIR滤波器的基本结构
5.2.1 IIR滤波器的特点
h(n)是无限长的， h(n)、 ak 、bk 均为实数 有某一个ak 0（k=1,2,3，…,N），即一定是递归结构 至少有一个ak 0（k=1,2,3，…,N），故H(z)在有限z平 面上一定有极点存在，即不一定稳定
若只有b0 0，其它bk = 0（k=1,2,3，…,N），称为全极 点型IIR滤波器，或自回归（AR）系统
试画出其直接型、
H(z)
33
级联型和并联型结构。
(1 1 z 1)(1 1 z 1 1 z 2 )
3
22
解： 并联型
将系统函数H(z)表达为 H (z) 2
1 z 1
部分分式之和的形式
1 1 z 1 1 1 z 1 1 z 2
1
3
22
2
一阶节
1/ 3
z 1
yy([nk)]
x(xn[k) ]
x[k] A
11 21
z 1 11 z 1 21
1L
z
1
1L
2L
z
1
2L
数字信号处理-第五章数字滤波器的基本结构(new)
1 2 ( 1 p z ) ( 1 z z ) 1k k 2k 1 k 1 N1 1 1 2 ( 1 c z ) ( 1 a z a z ) k 1k 2k k 1 k 1 k 1 N2 M1 M2
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
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
数字信号处理第5章
第5章 数字滤波器的基本结构5.1 学习要求1 掌握IIR 数字滤波器的基本网络结构，包括直接型、级联型和并联型；2 掌握FIR 数字滤波器的基本网络结构，包括直接型、级联型和频率抽样型；3 了解数字信号处理中的量化效应和数字信号处理的实现。
5.2 学习要点5.2.1 数字滤波器的结构特点与表示方法一个数字滤波器可以用系数函数表示为：01()()()1Mkk k N kk k b zY z H z X z a z -=-===-∑∑ (5-1) 直接由此式可得出表示输入输出关系的常系数线性差分方程为：1()()()N Mk k k k y n a y n k b x n k ===-+-∑∑ (5-2)由式(5-2)看出，实现一个数字滤波器需要几种基本的运算单元—加法器、单位延时和常数乘法器。
这些基本的单元可以有两种表示法：方框图法和信号流图法，如图5-1所示。
用方框图表示较明显直观，用流图表示则更加简单方便。
z ⊕aa单位延时乘常数相加方框图表示法信号流图表示法图5-1 基本运算过程的表示5.2.2 无限长单位脉冲响应(IIR)滤波器的基本结构无限长单位脉冲响应(IIR)滤波器有以下几个特点：(1) 系统的单位脉冲响应()h n 是无限长的；(2) 系统函数()H z 在有限z 平面(0z <<∞)上有极点存在； (3) 结构上存在着输出到输入的反馈，也就是结构上是递归型。
同一种系统函数()H z 的基本网络结构有直接I 型、直接Ⅱ型、级联型和并联型四种。
1直接I 型直接型按式(5-2)差分方程式将输入采样值（序列）)(n x 延迟并乘以系数k b ，将输出采样（序列）)(n y 延迟并乘以系数k a ，再把它们加起来，这种结构称为直接I 型，结构流图如图5-2所示。
由图可看出，总的网络)(z H 由Mkk k b z-=∑和11Nkk k a z-=-∑两部分网络级联组成，第一个网络实现零点，第二个网络实现极点，从图中又可看出，直接I 型结构需要N M +级延时单元。
数字信号处理第五章chhy
1 1 H1 ( z ) 1 H1 ( z ) 1 0.8 z 1 0.15z 2 1 2 1 1 0.8 z 0.15z 2 1.5 2.5 1.5 2.5 H 2 ( z) H 2 ( z ) 1 0.3z 1 1 0.5 z 1 1 1 2 1 0.3z 1 1 0.5 z 1 1 1 1 1 H3 ( z) H 3 ( z ) 1 0.3z 1 1 0.5 z 1 1 1 3 1 0.3z 1 1 0.5 z 1
5.3 IIR系统的基本网络结构
例：已知IIR数字滤波器的系统函数，画出该滤波器的直接型结构。
8 4 z 1 11z 2 2 z 3 H ( z) 5 3 1 1 z 1 z 2 z 3 4 4 8
5 3 1 y (n 1) y (n 2) y (n 3) 4 4 8 5 3 1 y (n) y (n 1) y (n 2) y (n 3) 8 x(n) 4 x(n 1) 11x(n 2) 2 x(n 3) y (n)
i 0 i 1
M
N
对输入信号的直接算法，
其系统函数H(z)为： M
刻以前的y(n-i)，可以 y (n ) bi x (n i ) ai y (n
i 0 递推出y(n) i 1
已知x(n)、ai、bi和n时 M
N
N i i 0 YY z z Y i Y( 1) ii aNz NaaYY ( ()z )zjiz,Hj((,zH z )( z( ))) ( z)0 i bi z i j i i i0 i ) ) X X z i ( z )z Y ( z )Y ( zY (bi bi( z()z (z z z ja j Y z , H z ) ( z bi X ) i 0 H ( z ) X ( )) X ( z )NN j jN i 1 j 1 1 X (z z N i 0 i 0 i 0 j j i=11 a j j z i X ( z ) 11 a1z a 1 j1 ai z j 1 j 1
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
LPAF
H ( j )
c
HPAF
c
H ( j )
c
BPAF
c
H ( j )
c1
BSAF
c2
H ( j )
c2 c1 c1
c2
5、数字滤波器的理想幅频特性
LPAF
H ( e j )
……… ……… ………
2 c c
IIR DF类型有：
直接型、级联型、并联型。 直接型结构：
直接I型、直接II型（正准型、典范型）
1、 IIR DF系统函数及差分方程
一个N阶IIR DF有理的系统函数可能表示为：
H (z)
b Z
i 0 N i i 1
M
i
1 ai Z
Байду номын сангаас
i
Y (z) X (z)
注：以下我们讨论M<=N情况。
H (z)
bi Z i 1 ai Z i
i 1 i 0 N
M
A
(1 ci z 1 ) (1 d i z 1 )
i 1 i 1 N
M
零、极点ci 和d i 只有两种情况： (a )或者是实根 (b)或者是共轭复根
系统函数系数分析
H (z) A (1 ci z 1 ) (1 d i z 1 )
3、直接II型（正准型/典范型）
直接II型原理 一个线性时不变系统，若交换其级联子系统的次 序，系统函数不变。把此原理应用于直接I型结构。 即： （1）交换两个级联网络的次序。
（2）合并两个具有相同输入的延时支路。
直接II型的结构流图：过程1--对调
x(n)
b0 a1 a2
y(n) z-1
β1M
a21
则这一系统差分方程为：
y( n) a i y( n i ) bi x( n i )
i 1 i 0
N
M
2、直接I型
直接I型流图
y( n) a i y( n i ) bi x( n i )
i 1 i 0
N
M
IIR DF的差分方程就代表了一种最直接的计算公式，用流图表
其方框图及流图结构如下：
x(n) b0 y(n)
x(n) b0
a1 a2
z-1 z-1
y(n)
a1
a2
z-1 z-1
说明：可通过流图或方框图看出系统的运算步骤和运算结构。 以后我们用流图来分析数字滤波器结构。
四、数字滤波器的分类
滤波器的种类很多，分类方法也不同。 1、从功能上分；低通、带通、高通、带阻。 2、从实现方法上分：FIR、IIR 3、从设计方法上来分： Butterworth（巴特沃斯）、 Chebyshev(切比雪夫）、 Ellips（椭圆）等。 4、从处理信号分：经典滤波器、现代滤波器
注：本课程主要讲经典滤波器
3、模拟滤波器和数字滤波器
经典滤波器从功能上分又可分为： 1、低通滤波器（LPAF/LPDF)
(Low pass analog filter / Low pass digital filter)
2、高通滤波器（HPAF/HPDF)
(High pass analog filter / High pass digital filter)
五、研究数字滤波器结构意义
滤波器的基本特性（如有限长冲激响应FIR与无限 长冲激响应IIR）决定了结构上有不同的特点。 不同结构所需的存储单元及乘法次数不同，前者影 响复杂性，后者影响运算速度。
有限精度（有限字长）实现情况下，不同运算结构 的误差及稳定性不同。
好的滤波器结构应该易于控制滤波器性能，适合于 模块化实现，便于时分复用。
8 z 4 z 11z 2 5 2 3 1 3 z z z 4 4 8
3 2
8 4 z 1 11z 2 2 z 3 5 1 3 2 1 3 1 z z z 4 4 8
必须将H(z)代为z-1的有理式
8 4 z 11z 2 z H (z) 5 1 3 2 1 3 1 z z z 4 4 8
第二节 IIR DF的基本结构
一、IIR DF特点 1、单位冲激响应h(n)是无限长的：n→∞ 2、系统函数H(z)在有限长z平面（0<|z|<∞) 有极点存在。 3、结构上存在输出到输入的反馈，也即结构上 是递归型的。 4、因果稳定的IIR滤波器其全部极点一定在单位圆内。
二、IIR DF基本结构
1、经典滤波器
假定输入信号x(n)中的有用成分和希望去除的成分，各自 占有不同的频带。 当x(n)经过一个线性系统（即滤波器）后即可将欲去除的 成分有效地去除。
|X(ej)|
无用 有用 c
|H(ej)|
|Y(ej)|
c
c
如果信号和噪声的频谱相互重叠，那么经典滤波器将 无能为力，此时可以设计现代滤波器来解决。
x(n) a1 对调 a2 z-1 z-1
b0 z-1 b1 z-1 b2
y(n)
z-1 b1 z-1 z-1 b2 bM
z-1
a N-1 z-1 aN z-1
-1 a N-1 z-1 z bM
aN
z-1
第一部分 对调
第二部分
直接II型的结构流图：过程2—合并
x(n)
a1 z-1
b0
z-1 b1
z-1
z-1
b2
bM
a2
a N-1 z-1
aN
z-1
直接I型DF结构的特点
1、两个网络级联：第一个横向结构M节延时网络实 现零点，第二个有反馈的N节延时网络实现极点。
2、共需(N+M)级延时单元。
3、系数ai、bi不是直接决定单个零极点，因而不能很 好地进行滤波器性能控制。 4、极点对系数的变化过于灵敏，从而使系统频率响应 对系统变化过于灵敏，也就是对有限精度（有限字 长）运算过于灵敏，容易出现不稳定或产生较大误 差。 注：直接型结构多用于低阶（2~3阶）滤波器。
二、数字滤波器的工作原理
设：x(n)是系统的输入，X(ej)是其傅立叶变换； y(n)是系统的输出，Y(ej)是其傅立叶变换； 则：
x(n)
LTI系统的输出为：
h(n)
y(n)
y( n)
m
x( m )h( n m ) F 1[ X (e j ) H (e j )]
3、带通滤波器（BPAF/BPDF)
(Bandpass analog filter / Bandpass digital filter)
4、带阻滤波器（BSAF/BSDF)
(Bandstop analog filter / Bandstop digital filter)
4、模拟滤波器的理想幅频特性
若将每一对共轭因子合并起来构成一个实系数的二阶因子，则：
H (z) A
(1 g i z 1 ) (1 1i z 1 2 i z 2 ) (1 pi z 1 ) (1 1i z 1 2 i z 2 )
i 1 i 1 i 1 N1 i 1 N2
M1
M2
基本二阶节的级联结构
H (z) A (1 g i z 1 ) (1 1i z 1 2 i z 2 ) (1 pi z 1 ) (1 1i z 1 2 i z 2 )
i 1 i 1
1 N1
M1 i 1 N1
M2
直接I型
x(n) 8 y(n) 5/4 -3/4 1/8 x(n) Z-1 5/4 Z-1 11 -2 Z-1 -3/4 1/8 11
1
2
3
直接II型
8 -4 y(n)
Z-1
Z-1 Z-1
-4
Z-1
Z-1
Z-1
-2
注意反馈部分系数符号
4、级联型结构
系统函数因式分解
一个N阶系统函数可用它的零、极点来表示, 即系统函数的分子、分母进行因式分解：
i 1 i 1 N M
A
(1 g i z 1 ) (1 hi z 1 )(1 hi* z 1 ) (1 pi z 1 ) (1 qi z 1 )(1 qi* z 1 )
i 1 i 1 i 1 N1 i 1 N2
M1
M2
式中：gi , pi为实根；hi , qi为复根。 其中 : N 1 2 N 2 N；M 1 2 M 2 M
第五章 数字滤波器结构 Digital Filter Structures
第一节 引言
一、什么是数字滤波器
顾名思义：其作用是对输入信号起到滤波的作用； 即DF是由差分方程描述的一类特殊的 离散时间系统。 功能： 把输入序列通过一定的运算变换成输出 序列。不同的运算处理方法决定了滤波 器的实现结构的不同。
三、数字滤波器表示方法
表示方法：方框图表示法、流图表示法
三种运算：相加、乘以常数、延时 基本运算单元：加法器、单位延时、乘常数的 乘法器。
1、方框图、流图表示法 方框图表示法 信号流图表示法
相加
a
乘常数
a
延时
z-1
z-1
例：二阶数字滤波器：
y(n) a1 y(n 1) a2 y(n 2) b0 x(n)
2、现代滤波器
它主要研究内容是从含有噪声的数据记录（又称时间序列） 中估计出信号的某些特征或信号本身。一旦信号被估计出， 那么估计出的信号将比原信号会有高的信噪比。
现代滤波器把信号和噪声都视为随机信号，利用它们的统计 特征（如自相关函数、功率谱等）导出一套最佳估值算法， 然后用硬件或软件予以实现。 现代滤波器理论源于维纳在40年代及其以后的工作，这一类 滤波器的代表为维纳滤波器，此外，还有卡尔曼滤波器、线 性预测器、自适应滤波器。
H ( j )
c
HPAF
c
H ( j )
c
BPAF
c
H ( j )
c1
BSAF
c2
H ( j )
c2 c1 c1
c2
5、数字滤波器的理想幅频特性
LPAF
H ( e j )
……… ……… ………
2 c c
IIR DF类型有：
直接型、级联型、并联型。 直接型结构：
直接I型、直接II型（正准型、典范型）
1、 IIR DF系统函数及差分方程
一个N阶IIR DF有理的系统函数可能表示为：
H (z)
b Z
i 0 N i i 1
M
i
1 ai Z
Байду номын сангаас
i
Y (z) X (z)
注：以下我们讨论M<=N情况。
H (z)
bi Z i 1 ai Z i
i 1 i 0 N
M
A
(1 ci z 1 ) (1 d i z 1 )
i 1 i 1 N
M
零、极点ci 和d i 只有两种情况： (a )或者是实根 (b)或者是共轭复根
系统函数系数分析
H (z) A (1 ci z 1 ) (1 d i z 1 )
3、直接II型（正准型/典范型）
直接II型原理 一个线性时不变系统，若交换其级联子系统的次 序，系统函数不变。把此原理应用于直接I型结构。 即： （1）交换两个级联网络的次序。
（2）合并两个具有相同输入的延时支路。
直接II型的结构流图：过程1--对调
x(n)
b0 a1 a2
y(n) z-1
β1M
a21
则这一系统差分方程为：
y( n) a i y( n i ) bi x( n i )
i 1 i 0
N
M
2、直接I型
直接I型流图
y( n) a i y( n i ) bi x( n i )
i 1 i 0
N
M
IIR DF的差分方程就代表了一种最直接的计算公式，用流图表
其方框图及流图结构如下：
x(n) b0 y(n)
x(n) b0
a1 a2
z-1 z-1
y(n)
a1
a2
z-1 z-1
说明：可通过流图或方框图看出系统的运算步骤和运算结构。 以后我们用流图来分析数字滤波器结构。
四、数字滤波器的分类
滤波器的种类很多，分类方法也不同。 1、从功能上分；低通、带通、高通、带阻。 2、从实现方法上分：FIR、IIR 3、从设计方法上来分： Butterworth（巴特沃斯）、 Chebyshev(切比雪夫）、 Ellips（椭圆）等。 4、从处理信号分：经典滤波器、现代滤波器
注：本课程主要讲经典滤波器
3、模拟滤波器和数字滤波器
经典滤波器从功能上分又可分为： 1、低通滤波器（LPAF/LPDF)
(Low pass analog filter / Low pass digital filter)
2、高通滤波器（HPAF/HPDF)
(High pass analog filter / High pass digital filter)
五、研究数字滤波器结构意义
滤波器的基本特性（如有限长冲激响应FIR与无限 长冲激响应IIR）决定了结构上有不同的特点。 不同结构所需的存储单元及乘法次数不同，前者影 响复杂性，后者影响运算速度。
有限精度（有限字长）实现情况下，不同运算结构 的误差及稳定性不同。
好的滤波器结构应该易于控制滤波器性能，适合于 模块化实现，便于时分复用。
8 z 4 z 11z 2 5 2 3 1 3 z z z 4 4 8
3 2
8 4 z 1 11z 2 2 z 3 5 1 3 2 1 3 1 z z z 4 4 8
必须将H(z)代为z-1的有理式
8 4 z 11z 2 z H (z) 5 1 3 2 1 3 1 z z z 4 4 8
第二节 IIR DF的基本结构
一、IIR DF特点 1、单位冲激响应h(n)是无限长的：n→∞ 2、系统函数H(z)在有限长z平面（0<|z|<∞) 有极点存在。 3、结构上存在输出到输入的反馈，也即结构上 是递归型的。 4、因果稳定的IIR滤波器其全部极点一定在单位圆内。
二、IIR DF基本结构
1、经典滤波器
假定输入信号x(n)中的有用成分和希望去除的成分，各自 占有不同的频带。 当x(n)经过一个线性系统（即滤波器）后即可将欲去除的 成分有效地去除。
|X(ej)|
无用 有用 c
|H(ej)|
|Y(ej)|
c
c
如果信号和噪声的频谱相互重叠，那么经典滤波器将 无能为力，此时可以设计现代滤波器来解决。
x(n) a1 对调 a2 z-1 z-1
b0 z-1 b1 z-1 b2
y(n)
z-1 b1 z-1 z-1 b2 bM
z-1
a N-1 z-1 aN z-1
-1 a N-1 z-1 z bM
aN
z-1
第一部分 对调
第二部分
直接II型的结构流图：过程2—合并
x(n)
a1 z-1
b0
z-1 b1
z-1
z-1
b2
bM
a2
a N-1 z-1
aN
z-1
直接I型DF结构的特点
1、两个网络级联：第一个横向结构M节延时网络实 现零点，第二个有反馈的N节延时网络实现极点。
2、共需(N+M)级延时单元。
3、系数ai、bi不是直接决定单个零极点，因而不能很 好地进行滤波器性能控制。 4、极点对系数的变化过于灵敏，从而使系统频率响应 对系统变化过于灵敏，也就是对有限精度（有限字 长）运算过于灵敏，容易出现不稳定或产生较大误 差。 注：直接型结构多用于低阶（2~3阶）滤波器。
二、数字滤波器的工作原理
设：x(n)是系统的输入，X(ej)是其傅立叶变换； y(n)是系统的输出，Y(ej)是其傅立叶变换； 则：
x(n)
LTI系统的输出为：
h(n)
y(n)
y( n)
m
x( m )h( n m ) F 1[ X (e j ) H (e j )]
3、带通滤波器（BPAF/BPDF)
(Bandpass analog filter / Bandpass digital filter)
4、带阻滤波器（BSAF/BSDF)
(Bandstop analog filter / Bandstop digital filter)
4、模拟滤波器的理想幅频特性
若将每一对共轭因子合并起来构成一个实系数的二阶因子，则：
H (z) A
(1 g i z 1 ) (1 1i z 1 2 i z 2 ) (1 pi z 1 ) (1 1i z 1 2 i z 2 )
i 1 i 1 i 1 N1 i 1 N2
M1
M2
基本二阶节的级联结构
H (z) A (1 g i z 1 ) (1 1i z 1 2 i z 2 ) (1 pi z 1 ) (1 1i z 1 2 i z 2 )
i 1 i 1
1 N1
M1 i 1 N1
M2
直接I型
x(n) 8 y(n) 5/4 -3/4 1/8 x(n) Z-1 5/4 Z-1 11 -2 Z-1 -3/4 1/8 11
1
2
3
直接II型
8 -4 y(n)
Z-1
Z-1 Z-1
-4
Z-1
Z-1
Z-1
-2
注意反馈部分系数符号
4、级联型结构
系统函数因式分解
一个N阶系统函数可用它的零、极点来表示, 即系统函数的分子、分母进行因式分解：
i 1 i 1 N M
A
(1 g i z 1 ) (1 hi z 1 )(1 hi* z 1 ) (1 pi z 1 ) (1 qi z 1 )(1 qi* z 1 )
i 1 i 1 i 1 N1 i 1 N2
M1
M2
式中：gi , pi为实根；hi , qi为复根。 其中 : N 1 2 N 2 N；M 1 2 M 2 M
第五章 数字滤波器结构 Digital Filter Structures
第一节 引言
一、什么是数字滤波器
顾名思义：其作用是对输入信号起到滤波的作用； 即DF是由差分方程描述的一类特殊的 离散时间系统。 功能： 把输入序列通过一定的运算变换成输出 序列。不同的运算处理方法决定了滤波 器的实现结构的不同。
三、数字滤波器表示方法
表示方法：方框图表示法、流图表示法
三种运算：相加、乘以常数、延时 基本运算单元：加法器、单位延时、乘常数的 乘法器。
1、方框图、流图表示法 方框图表示法 信号流图表示法
相加
a
乘常数
a
延时
z-1
z-1
例：二阶数字滤波器：
y(n) a1 y(n 1) a2 y(n 2) b0 x(n)
2、现代滤波器
它主要研究内容是从含有噪声的数据记录（又称时间序列） 中估计出信号的某些特征或信号本身。一旦信号被估计出， 那么估计出的信号将比原信号会有高的信噪比。
现代滤波器把信号和噪声都视为随机信号，利用它们的统计 特征（如自相关函数、功率谱等）导出一套最佳估值算法， 然后用硬件或软件予以实现。 现代滤波器理论源于维纳在40年代及其以后的工作，这一类 滤波器的代表为维纳滤波器，此外，还有卡尔曼滤波器、线 性预测器、自适应滤波器。