数字信号处理 第五章
数字信号处理-第五章
系 统 函 数
：
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）
第五章 时域离散系统的基本网络结构
本章的主要内容就是描述数字滤波器的基 本网络结构。（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的级联型例题
数字信号处理-原理实现及应用(高西全-第3版)第5章 信号的相关函数及应用
rxy (m) ryx (m)
性质2 rxy (m) rx (0)ry (0) ExEy
性质3
lim
m
rxy (m)
0
因为一般能量信号都是有限非零时宽的，所以，当 m 时，二者的非零区不重叠， 所以，该性质成立。
信息与通信工程系—数字信号处理
2.自相关函数性质
性质1
若 x(n) 是实信号，则 rx (m)是实偶函数，即
[h(m) h(m)][x(m) x(m)]
rh (m) rx (m)
ry (0) rh (m) rx (m) m0
= rh (n)rx (n m) = rh (n)rx (n)
n
m0 n
系统稳定，则h(n)为能量信号
rh (m) 存在；
如果 rx (m) 存在，则 ry (m) 存在。
观测信号 y(n) x(n) w(n)；y(n) 的自相关函 ry (m)
(a) 2
w(n)
0
-2 10 20 30 40 50 60 70 80 90 100 n
(b) 2
y(n)
0
-2
10
20
30
40
50 n
60
70
80
90 100
噪声自相关
(c)
函数导致
1
ry(m)
0
-1
-50 -40 -30 -20 -10
h(m) [x(m) x(m)]
h(m) rx (m)
所以，ryx (m)可以看成线性时不变系统对输入序列的响应输出。
rx (m)
LTI系统 h(n)
ryx (m)
信息与通信工程系—数字信号处理
系统输出信号的自相关函数:
第五章 数字信号处理- 微弱信号处理
第五章微弱信号处理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 +
《数字信号处理导论_第5章》
1 H ( ) b(n )cos n 2 n 1
N /2
1 时 cos n 0 2
则 H ( ) 0 z 1是零点
H ( )对 0, 2 呈偶对称 H ( )对 呈奇对称
0
10
20
y2 ( n )
0
-2 -10
0
10
20
定义：
d ( ) g ( ) d
为系统的群延迟 (Group Delay, GD)
显然，若系统具有线性相位，则其GD为
常数。
GD可作为相频响应是否线性的一种度量，同 时，它也表示了系统输出的延迟。
若： 则：
x(n) xa (n) cos( 0 n), c 0 x(n) : Narrowband Signal
z e j
j N21 N 1 N 1 " " h(n)cos 2 n e n 0 N 1 N 1 j je 2 h ( n )sin N 1 n " " 2 n 0
为第一类线性相位
N 1 2
2）h(n)奇对称
h( n) h( N 1 n)
j N 1 N 1 2
频率响应：
j
N 1 H (e ) H ( z ) z e j je h(n)sin 2 n n 0 N 1 j j N 1 N 1 2 2 e h(n)sin 2 n n 0
z 1为零点 故不能设计成高通、带阻滤波器
3）h(n)奇对称，N为奇数
数字信号处理第五章-IIR数字滤波器的设计
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)。
精品文档-数字信号处理(吴瑛)-第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 非理想低通滤波器因果脉冲响应的幅频响应
数字信号处理 第五章07
an example
2
2
1.5
h(n)
1.5
a(n)
1
1
0.5
0.5
0
0
-0.5
-0.5
-1
-1
-1.5
-1.5
-2
0
2
4
6
-2
0
2
4
6
4 2 0 -2 -4
H()
0 -5 -10 -15 -20
()
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.2
0.4
0.6
0.8
h(n)
2 2
d(n)
1
1
0
0
-1
-1
-2
-2
-3
0
2
4
6
-3
0
2
4
6
5
0
H()
0
-5 -10 -15
()
-5
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
-20
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
5
4 2
arg[H(ej)]
0
0
|H(ej)|
-5 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
h(n ) e
j n
e
j (N n 1)
数字信号处理 第五章
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 •
数字信号处理 第五章
第五章离散时间信号的数字处理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 ）
第五章 有限离散傅立叶变换与快速傅立叶变换
第五章 有限离散傅立叶变换与快速傅立叶变换
数字信号处理第五章 FIR滤波器、窗函数法
(5.23)
0
h(n)偶对称N为偶数时的频 率响应特性
不适用高通或带阻滤波器
2）H(w)以
为奇对称，以
偶对称
以
为奇对称，以
为偶对称
22
1. 成立条件
5.2.3 恒群延时单独成立
解线性微分方程
＝
的取值条件？
0
的取值条件？
有什么特性？
(5.24)
23
离散傅立叶变换： 由相位定义： 化简得：
(5.25)
n
28
相位函数：
(5.30) 线性相位
幅度函数：
(5.31)
其中
1）H(w)在
关于
处奇对称
奇对称
2）传输函数H(z)在
处有零点
0
h(n)奇对称N为奇数时的频 率响应特性
不适用于低通、高通、带阻滤波器
29
3. h(n)奇对称N为偶数（线性相位IV型）时的频率响应
所以
456 7
0123
n
(5.36)
数字信号处理 Digital Signal Processing
第五章 FIR滤波器设计和实现
Finite Impulse Response Filter Design and Implementation
2
FIR数字滤波器
• Finite Impulse Response (FIR)：有限冲激响应，非递归 滤波器，输出只与当前和有限个过去输入有关
偶对称
0 h(n)偶对称N为奇数时的 频率响应特性
20
3. h(n)偶对称N为偶数（线性相位II型）时的频率响应
分拆
0123456 7 n
变量代换及
[工学]第五章 数字信号处理1AD_DA_采样定律、能量泄漏
![[工学]第五章 数字信号处理1AD_DA_采样定律、能量泄漏](https://img.taocdn.com/s3/m/b3343a3efad6195f312ba6e1.png)
西安工业大学机电学院
（2)采样定理
上述情况表明，如果fs＞2fm，就不发生频混现象，因此对 采样脉冲序列的间隔Ts须加以限制，即采样频率fs(1／Ts） 必须大于信号x(t)中的最高频率fm的两倍，即
f s ＞ 2 fm 或 ω s ＞ 2 ω m
为了保证采样后的信号能真实地保留原始模拟信号 的信息，采样信号的频率必须大于原信号中最高频 率成分的2倍。这是采样的基本法则，称为采样定 理。
西安工业大学机电学院
3.3 采样定理
采样是将采样脉冲序列p(t)与连续时间信 号x(t)相乘，取离散点x(nt)的值的过程。
西安工业大学机电学院
X(0), X(1), X(2), ……, X(n)
西安工业大学机电学院
1)采样信号的频谱 采样：等间隔的脉冲序列p(t)与连续时间信号x(t)相乘。 理想脉冲采样过程如下：
当采样信号的频率低于被采样信号的最高频率时，采样所得的信 号中混入了虚假的低频分量，这种现象叫做频率混叠。
西安工业大学机电学院
(1) 频混现象 频混现象又称频谱混叠效应：因采样频率过低造成信号失 真的现象。
f
频带范围为-fm～fm
f s f
其周期为fs＝1/Ts
f s f
f
西安工业大学机电学院
•当采样周期Ts较小时，fs＞2fm，周期谱图相互分离, 如图b； •当采样周期Ts较大时，fs<2fm，周期谱图相互重叠， 即谱图之间高频与低频部分发生重叠，如图 c，此即 频混现象，这将使信号复原时丢失原始信号中的高频 信息。
X ( f kf )
s
•一个连续信号经过理想采样以后，它的频谱将沿着 频率轴每隔一个采样频率fs ,重复出现一次，即其 频谱产生了周期延拓。
数字信号处理 第5章
【例5.2.1】 求图5.2.2（a）信号流图决定的系统函数H(z)。
解 图5.2.2（a）信号流图的节点变量方程为(5.2.1)式， 对 （5.2.1）式进行Z变换，得到：
w1 (n) w2 (n 1) w (n) w (n－1) 2 2 －a －a w2 (n) x(n) 1 w2 (n) 2 w1 (n) y (n) b2 w1 (n) b1 w2 (n) b0 w2 (n)
23
第4章 FFT
Hj(z)如下式：
H j ( z)
0j 1j z 1 2 j z 2
1 a1 j z a2 j z
1 2
（5.3.2）
式中，β0j、β1j、β2j、α1j和α2j均为实数。这样H(z)就分解
成一些一阶或二阶的子系统函数的相乘形式：
H ( z ) H1 ( z ) H 2 ( z ) H k ( z )
的。例如，一个简单的一阶IIR网络的差分方程为
y(n) ay(n 1) x(n)
其单位脉冲响应h(n)=anu(n)。这两类不同的网络结构各 有不同的特点，下面分类叙述其网络结构。
14
第4章 FFT
5.3 无限长脉冲响应基本网络结构
IIR网络的基本网络结构有三种：直接型、级联型、并 联型。 1． 直接型 将N阶差分方程重写如下：
第5章 时域离散系统的网络结构
5.1 引言 5.2 用信号流图表示网络结构 5.3 无限长脉冲响应基本网络结构
5.4 有限长脉冲响应基本网络结构
5.5 线性相位结构
5.6 频率采样结构
5.7 格型网络结构
1
第4章 FFT
5.1 引方程、单位脉冲
数字信号处理-第五章数字滤波器的基本结构(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
数字信号处理第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)。
第五章离散时间信号的数字处理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,。
评论这种改变后产生的曲线。
t的范围是[-0.5 1.5]，增量是0.004.在图中，t的范围是[0 1]。
T=[0 0.5]如图2，改变后的图跟原来的图显示的波形差了点。
时间T=0.5后边的只抽样。
Q5.8 恢复原始显示范围并通过分别改变正弦信号的频率3Hz和7Hz，重复程序P5.2。
相应的等效离散时间信号与习题Q5.6中产生的离散时间信号有差别么？若没有，为什么没有？f=3Hz f=7Hz从上图可以看出两个频率得到的结果没有区别；因为选择的抽样周期相同，频率的不同不会影响离散时间信号。
Q5.9 在程序P5.3中，连续时间函数xa[t]是什么？xa[t]的连续时间傅里叶变换是如何计算的？连续时间函数：xa = 2*t.*exp(-t),计算：wa = 0:10/511:10;ha = freqs(2,[1 2 1],wa);Q5.10运行程序P5.3，产生并显示离散时间信号及其连续时间等效，以及它们各自的傅里叶变换。
有任何明显的的混叠影响吗？clf;t = 0:0.005:10;xa = 2*t.*exp(-t);subplot(2,2,1)plot(t,xa);gridxlabel('时间, msec');ylabel('振幅');title('连续时间序号 x_{a}(t)');subplot(2,2,2)wa = 0:10/511:10;ha = freqs(2,[1 2 1],wa);plot(wa/(2*pi),abs(ha));grid;xlabel('频率, kHz');ylabel('振幅');title('|X_{a}(j\Omega)|');axis([0 5/pi 0 2]);subplot(2,2,3)T = 1;n = 0:T:10;xs = 2*n.*exp(-n);k = 0:length(n)-1;stem(k,xs);grid;xlabel('时间序号 n');ylabel('振幅');title('离散时间信号 x[n]');subplot(2,2,4)wd = 0:pi/255:pi;hd = freqz(xs,1,wd);plot(wd/(T*pi), T*abs(hd));grid;xlabel('频率, kHz');ylabel('振幅');title('|X(e^{j\omega})|');axis([0 1/T 0 2])图1 图2从图1可以看出，没有混叠影响。
Q5.11 将抽样周期增加到1.5，再运行程序P5.3.有任何明显的混叠影响吗？从图2看出，没有明显的混叠影响。
Q5.12 对于的影响，修改程序P5.3，并重做系统Q5.10和Q5.11将程序部分修改xa = exp(-pi.*t.*t)，得到如下两图。
T=1 T=1.5Q5.13 在程序P5.4中，通带波纹RP和最小阻带衰减RS是多少dB?通带和阻带边界频率是多少Hz？通带波纹RP=0.5dB，最小阻带衰减RS=30dB；通带边界频率3500Hz，阻带边界频率是4500Hz。
Q5.14 运行程序P5.4并显示增益响应。
所设计的滤波器满足给定的指标么？所涉及的滤波器阶数N和单位为Hz的3dB截止频率是多少？clf;Fp = 3500;Fs = 4500;Wp = 2*pi*Fp; Ws = 2*pi*Fs;[N, Wn] = buttord(Wp, Ws, 0.5, 30,'s');[b,a] = butter(N, Wn, 's');disp(‘N=’);disp(N);disp(‘Wn=’);disp(Wn)wa = 0:(3*Ws)/511:3*Ws;h = freqs(b,a,wa);plot(wa/(2*pi), 20*log10(abs(h)));gridxlabel('频率, Hz');ylabel('增益, dB');title('增益响应');axis([0 3*Fs -60 5]);满足给定的指标。
所涉及的滤波器阶数N是18，单位为Hz的3dB截止频率是2333.8Hz。
Q5.15 用cheb1ord和cheby1修改程序P5.4，以设计与程序P5.4有着相同指标的一个切比雪夫1型低通滤波器。
运行修改的程序并显示增益响应，所设计的滤波器满足给定的指标么？所涉及的滤波器阶数N和单位为Hz的通带边界频率是多少？[N, Wn] = cheb1ord(Wp, Ws, 0.5, 30,'s');[b,a] = cheby1(N, 30,Wn, 's');满足给定的指标Q5.16 用cheb2ord和cheby2修改程序P5.4，以设计与程序P5.4有着相同指标的一个切比雪夫2型低通滤波器。
运行修改的程序并显示增益响应，所设计的滤波器满足给定的指标么？所涉及的滤波器阶数N和单位为Hz的阻带边界频率是多少？[N, Wn] = cheb2ord(Wp, Ws, 0.5, 30,'s');[b,a] = cheby2(N, 30,Wn, 's');满足给定的指标Q5.17 用ellipord和ellip修改程序P5.4，以设计与程序P5.4有着相同指标的一个椭圆低通滤波器。
运行修改的程序并显示增益响应，所设计的滤波器满足给定的指标么？所涉及的滤波器阶数N和单位为Hz的通带边界频率是多少？[N, Wn] = ellipord(Wp,Ws,0.5, 30,'s');[b,a] = ellip(N,0.5,30,Wn, 's');满足给定的指标Q5.18 在程序P5.5中，运算符==的功能是什么？==的功能是判断后边字符是否恒等于前边字符。
Q5.19 用程序P5.5对于字长值6和8，产生如下十进制小数的源码形式的二进制等效（a）0.80165，（b）-0.80165，（c）0.64333，（d）-0.9125。
通过手算验证结果。
d = input('输入十进制小数 = ');b = input('输入所需的字长 = ');d1 = abs(d);beq = [zeros(1,b)];for k = 1:bint = fix(2*d1);beq(k) = int;d1 = 2*d1 - int;endif sign(d) == -1;bin = [1 beq];elsebin = [0 beq];enddisp('二进制等效是');disp(bin)Q5.20 使用程序P5.6确定习题Q5.19中产生的二进制小数的十进制等效。
你的答案和原十进制小数有多接近？bin = input('输入二进制小数 = ');b = length(bin) - 1; d = 0;for k = 1:bd = d + bin(k+1)*2^(-k);endif sign(bin(1)) == 0;dec = d;elsedec = - d;enddisp('十进制等效是 ');disp(dec);所得结果均为0，跟原来的不接近。
Q5.21 在程序P5.7中，运算符~的目的是什么？~目的是取反。
Q5.22 使用程序P5.7,确定并验证习题Q5.19中展开的二进制的反码表示。
bin = input('输入二进制数 = ');if sign(bin(1)) == 0;onescomp = bin;elsebin(1) = 0;onescomp = ~bin;enddisp('反码等效是');disp(onescomp);Q5.23 在程序P5.8中，运算符|和&的目的是什么？|代表或，&代表与。