离散系统的Z域分析
- 格式：doc
- 大小：181.50 KB
- 文档页数：8
下载文档原格式
合集下载
离散系统的Z域分析法
X(z)
D z-1
X(k-1)
z-1X(z)
X
注意:z域框图只能求系统零状态响应 注意 域框图只能求系统零状态响应
第
例题
1.求如图系统的单位响应 求如图系统的单位响应h(k)和单位阶跃响应 和单位阶跃响应g(k) 求如图系统的单位响应 和单位阶跃响应
8 页
2. 知 阶 散 统 初 条 为 zi (0) = 2, yzi (1) =1 已 二 离 系 的 始 件 y , 当 入 (k) = ε (k)时 输 f , 1 5 k k 输 全 应 (k) =[ + 4⋅ 2 − ⋅ 3 ]ε (k), 出 响 y 2 2 求 差 方 .， 画 系 框 。 此 分 程 并 出 统 图
求解线性时不变离散系统的差分方程有两种方法： 求解线性时不变离散系统的差分方程有两种方法： •时域方法: y(k) =y zi (k) + yzs (k) yzs (k) = h(k) ∗ x(k) 时域方法: 时域方法 •z变换方法: 变换方法: 变换方法 Yzs (z) = H(z) ⋅ X(z)
X
第
二.系统框图的z域分析法
基本思路: 基本思路 时域框图 z域框图 域框图 z域代数方程 域代数方程 Yzs(z)
7 页
yzs(k)
x(k) ⇒ X (z) yzs (k) ⇒ Yzs (z) 延迟单元 x(k)
x(k)ε (k) ↔ X(z) x(k −1)ε (k) ↔ z−1X(z) + x(−1)
y(k)
X
第 5 页
优点: 优点:
•差分方程经 变换→代数方程； 差分方程经z变换 代数方程； 差分方程经 变换→ •将时域卷积→z域乘积； 将时域卷积→ 域乘积； 将时域卷积 域乘积 •部分分式展开后求解z逆变换较容易； 部分分式展开后求解z 部分分式展开后求解 逆变换较容易； •z变换过程自动引入了系统初始状态（相当于0变换过程自动引入了系统初始状态（相当于0 变换过程自动引入了系统初始状态 的条件） 可同时求出零输入和零状态响应。 的条件）,可同时求出零输入和零状态响应。 ， 注意：z域求解系统只需 －状态[y(-1)，y(-2), …,] 注意： 域求解系统只需0 状态 域求解系统只需 时域求解系统要递推出0 状态确定待定系数。 时域求解系统要递推出 ＋状态确定待定系数。
D z-1
X(k-1)
z-1X(z)
X
注意:z域框图只能求系统零状态响应 注意 域框图只能求系统零状态响应
第
例题
1.求如图系统的单位响应 求如图系统的单位响应h(k)和单位阶跃响应 和单位阶跃响应g(k) 求如图系统的单位响应 和单位阶跃响应
8 页
2. 知 阶 散 统 初 条 为 zi (0) = 2, yzi (1) =1 已 二 离 系 的 始 件 y , 当 入 (k) = ε (k)时 输 f , 1 5 k k 输 全 应 (k) =[ + 4⋅ 2 − ⋅ 3 ]ε (k), 出 响 y 2 2 求 差 方 .， 画 系 框 。 此 分 程 并 出 统 图
求解线性时不变离散系统的差分方程有两种方法： 求解线性时不变离散系统的差分方程有两种方法： •时域方法: y(k) =y zi (k) + yzs (k) yzs (k) = h(k) ∗ x(k) 时域方法: 时域方法 •z变换方法: 变换方法: 变换方法 Yzs (z) = H(z) ⋅ X(z)
X
第
二.系统框图的z域分析法
基本思路: 基本思路 时域框图 z域框图 域框图 z域代数方程 域代数方程 Yzs(z)
7 页
yzs(k)
x(k) ⇒ X (z) yzs (k) ⇒ Yzs (z) 延迟单元 x(k)
x(k)ε (k) ↔ X(z) x(k −1)ε (k) ↔ z−1X(z) + x(−1)
y(k)
X
第 5 页
优点: 优点:
•差分方程经 变换→代数方程； 差分方程经z变换 代数方程； 差分方程经 变换→ •将时域卷积→z域乘积； 将时域卷积→ 域乘积； 将时域卷积 域乘积 •部分分式展开后求解z逆变换较容易； 部分分式展开后求解z 部分分式展开后求解 逆变换较容易； •z变换过程自动引入了系统初始状态（相当于0变换过程自动引入了系统初始状态（相当于0 变换过程自动引入了系统初始状态 的条件） 可同时求出零输入和零状态响应。 的条件）,可同时求出零输入和零状态响应。 ， 注意：z域求解系统只需 －状态[y(-1)，y(-2), …,] 注意： 域求解系统只需0 状态 域求解系统只需 时域求解系统要递推出0 状态确定待定系数。 时域求解系统要递推出 ＋状态确定待定系数。
信号与系统第8章 离散时间系统的z域分析
零状态响应为
Yf
(z)
(1 z 1 z 2 ) 2 3z 1 z 2
1 1 z 1
1/ 6 0.5 5 / 6 1 z1 1 z1 1 0.5z1
yf [k] Z 1{Yf (z)}{1/ 6 0.5(1)k (5/ 6)(0.5)k}u[k]
y[k] yx[k] yf [k] {1/ 6 3.5(1)k (4 / 3)(0.5)k}u[k]
离散时间信号与系统的Z域分析
• 离散时间信号的Z域分析 • 离散时间系统的Z域分析 • 离散时间系统函数与系统特
性
离散时间信号的Z域分析
• 理想取样信号的拉普拉斯变换 • 单边Z变换定义 • 单边Z变换的收敛域 • 常用序列的Z变换 • 单边Z变换的性质 • Z反变换
理想取样信号的拉普拉斯变换
fs (t) f (t) (t kT) f (kT) (t kT)
Re(z)
三、常用序列的Z变换
1) Z{ (k)} 1, z 0
2) 3)
Z{u(k)} 1 1 z
Z{aku(k)}
1 , 1
1 a
z
z
1
1 z
a
4)
Z{e
j0k
u(k
)}
1
e
1
j0
z
1
z z e j0
5)
Z{e-
j0k u (k
)}
1
1 e- j0
z
1
z z e- j0
z e j0 z e j0
解代数方程
二阶系统响应的z域求解
y[k] a1 y[k 1] a2 y[k 2] b0 f [k] b1 f [k 1] k 0
初始状态为y[1], y[2] 对差分方程两边做Z变换，利用
信号与系统实验(MATLAB 西电版)实验17 离散系统的Z域分析
(4)
实验17 离散系统的Z域分析
离散系统的分析方法可分为时域解法和变换域解法两大 类。其中离散系统变换域解法只有一种，即Z变换域解法。Z 变换域没有物理性质，它只是一种数学手段，之所以在离散 系统的分析中引入Z变换的概念，就是要像在连续系统分析 时引入拉氏变换一样，简化分析方法和过程，为系统的分析 研究提供一条新的途径。
F=ztrans(f)： 实现函数f(n)的Z变换，默认返回函数F是 关于z
F=ztrans(f，w)：实现函数f(n)的Z变换，返回函数F是关 于w
F=ztrans(f，k，w)：实现函数f(k)的Z变换，返回函数F是 关于w的函数。
实验17 离散系统的Z域分析
2. 单边逆Z变换函数iztrans 功能：iztrans可以实现信号F(z)的逆Z
实验17 离散系统的Z域分析
3) 一个离散LTI系统，差分方程为y(k)-0.81y(k-2)=f(k)-f(k-2)，
(1) 系统函数H(z); (2) 单位序列响应h(k)的数学表达式，并画出波形; (3) 单位阶跃响应的波形g(k); (4) 绘出频率响应函数H(ejθ)
实验17 离散系统的Z域分析
实验17 离散系统的Z域分析
MATLAB %确定信号的Z syms n z% f1=3^n f1_z=ztrans(f1) f2=cos(2*n) f2_z=ztrans(f2)；
实验17 离散系统的Z域分析
f1 = 3^n f1_z = 1/3*z/(1/3*z-1) f2 = cos(2*n) f2_z = (z+1-2*cos(1)^2)*z/(1+2*z+z^2-4*z*cos(1)^2)
极点图见图17.3
实验17 离散系统的Z域分析
实验17 离散系统的Z域分析
离散系统的分析方法可分为时域解法和变换域解法两大 类。其中离散系统变换域解法只有一种，即Z变换域解法。Z 变换域没有物理性质，它只是一种数学手段，之所以在离散 系统的分析中引入Z变换的概念，就是要像在连续系统分析 时引入拉氏变换一样，简化分析方法和过程，为系统的分析 研究提供一条新的途径。
F=ztrans(f)： 实现函数f(n)的Z变换，默认返回函数F是 关于z
F=ztrans(f，w)：实现函数f(n)的Z变换，返回函数F是关 于w
F=ztrans(f，k，w)：实现函数f(k)的Z变换，返回函数F是 关于w的函数。
实验17 离散系统的Z域分析
2. 单边逆Z变换函数iztrans 功能：iztrans可以实现信号F(z)的逆Z
实验17 离散系统的Z域分析
3) 一个离散LTI系统，差分方程为y(k)-0.81y(k-2)=f(k)-f(k-2)，
(1) 系统函数H(z); (2) 单位序列响应h(k)的数学表达式，并画出波形; (3) 单位阶跃响应的波形g(k); (4) 绘出频率响应函数H(ejθ)
实验17 离散系统的Z域分析
实验17 离散系统的Z域分析
MATLAB %确定信号的Z syms n z% f1=3^n f1_z=ztrans(f1) f2=cos(2*n) f2_z=ztrans(f2)；
实验17 离散系统的Z域分析
f1 = 3^n f1_z = 1/3*z/(1/3*z-1) f2 = cos(2*n) f2_z = (z+1-2*cos(1)^2)*z/(1+2*z+z^2-4*z*cos(1)^2)
极点图见图17.3
实验17 离散系统的Z域分析
第6章 离散时间系统的z域分析
1 | z | 1 2 | z | 2
例 求序列f (k ) cosh (2k ) (k )的z变换。
1 2k 由于 cosh ( k ) (e e 2 k ) 2 2 在单边指数序列a k ( k )的z变换中令a e 2 , 可得 z e (k ) , | z || e 2 | z e2 根据z变换的线性性质可得
f (k )
3
f ( k ) ( k ) 3
2
2
1
1 o 1 2
f ( k 1) 3 2
k
1 o 1 2
f ( k 1) ( k ) 3 2
1
k
1
1 o 1 2
f ( k 1)
k
1 o 1 2
f ( k 1) ( k )
3
k
3
2 1
1 o 1 2
k
1 o 1 2
k
(1)双边Z变换的移位 若 f (k ) F ( z )
k 0
该式称为单边Z变换。
将f ( k )的Z变换简记为Z [ f ( k )] ，象函数F ( z )的逆z变换 简记为Z
1
[ F ( z )] f ( k )与F ( z )两者间的关系简记为 ，
f (k ) F ( z )
在拉普拉斯变换分析中重点讨论了单边拉普拉斯 变换，这是由于在连续时间系统中，非因果信号 的应用较少。 对于离散系统，非因果信号也有一定的应用范围， 因此对单、双边z变换都进行讨论。
a
b
O
Re(z )
6.1.3 常见序列的Z变换
(k )
1
O
k
(k ) 1
离散时间信号与系统的Z域分析
《信号与系统》课程实验报告变换。
zz z z z z F 2112)(232+++-=一、实验原理的验证 1、离散系统零极点图实验原理如下：离散系统可以用差分方程描述：∑∑==-=-Mm m Ni i m k f b i k y a 0)()(Z 变换后可得系统函数：NN MM z a z a a z b z b b z F z Y z H ----++++++==......)()()(110110 可以用root 函数可分别求零点和极点。
例7-4 求系统函数零极点图131)(45+-+=z z z z H实验结果如下：2、离散系统的频率特性实验原理如下：离散系统的频率特性可由系统函数求出，既令ωj e z =,函数freqz 可计算频率特性，调用格式是：[H ，W]=freqz(b,a,n)，b 和a 是系统函数分子分母系数，n 是π-0范围内n 个等份点，默认值为512，H 是频率响应函数值，W 是相应频率点； 例7-5 系统函数z z z H 5.0)(-=10个频率点的计算结果为幅频特性曲线相频特性曲线freqz语句直接画图例7-7已知系统函数114/11)1(4/5)(----=z z z H ，画频率响应和零极点图。
零极点图幅频特性曲线相频特性曲线二、已知离散系统的系统函数如下所示：1422)(232+-++=z z z z z H试用MATLAB 实现下列分析过程： （1）求出系统的零极点位置；（2）绘出系统的零极点图，根据零极点图判断系统的稳定性； （3）绘出系统单位响应的时域波形，并分析系统稳定性与系统单位响应时域特性的关系。
（1）由计算结果可知：系统的极点为p0=-3.3028、p1=1、p2=0.3028。
由计算结果可知：系统的零点为z0=1.4142i 、z1=-1.4142i 。
（2）系统的零极点图如下：程序清单如下： a=[1 2 -4 1]; b=[1 0 2]; ljdt(a,b)p=roots(a)q=roots(b)pa=abs(p)由图可知：第一个极点（p0）在单位圆外部，第二个极点（p1）在单位圆上，第三个极点（p2）在单位圆内部，因为有一个极点在单位圆外部，故该系统是不稳定的系统（稳定系统要求极点全部在单位圆内）。
离散系统的Z域分析
z
k
cos(
0
k
)
k
z
z2 z2 z cos 2z2 2z cos 0
0
1
2
..........
k
sin 0k
k
z
2z2
z 2
sin 0 z cos 0
1 2
.........
k k
k
z (z )2k kk Nhomakorabea1
五、ZT & DTFT
求和收敛
设f(k)
为因果序列、则
F (e j ) f k e jk
Z eS Ts e e Ts jTs e j
k
F (z) f (k)zk k 0
e Ts
Ts
2 s
S 域中的一点→ → Z 域中的一点；Z 域中的一点→ → S 域中的无穷个点。
S 1 Ln z 1 Ln(e j ) 1 Ln j
Ts
Ts
Ts
Ts
三、收敛域: F (z) f k zk
ak (k) bk (k 1) z z ∣a∣< |z|< |b|
za zb
jIm[z]
|b|
|a|
o
Re[z]
四、常用 z 变换
(k+1) ←→z; (k-1) ←→z-1;……
(k) ←→1 (k) ←→z/(z-1) ←→ - (- k-1)
零、极 点分布
k k z k k 1
F(z)
K1 e j z
z e j
K1 e j z
z e j
若z> , f(k)=2K1kcos(k+)(k)，… …
(3) F(z)有重极点 推导记忆：
k
cos(
0
k
)
k
z
z2 z2 z cos 2z2 2z cos 0
0
1
2
..........
k
sin 0k
k
z
2z2
z 2
sin 0 z cos 0
1 2
.........
k k
k
z (z )2k kk Nhomakorabea1
五、ZT & DTFT
求和收敛
设f(k)
为因果序列、则
F (e j ) f k e jk
Z eS Ts e e Ts jTs e j
k
F (z) f (k)zk k 0
e Ts
Ts
2 s
S 域中的一点→ → Z 域中的一点；Z 域中的一点→ → S 域中的无穷个点。
S 1 Ln z 1 Ln(e j ) 1 Ln j
Ts
Ts
Ts
Ts
三、收敛域: F (z) f k zk
ak (k) bk (k 1) z z ∣a∣< |z|< |b|
za zb
jIm[z]
|b|
|a|
o
Re[z]
四、常用 z 变换
(k+1) ←→z; (k-1) ←→z-1;……
(k) ←→1 (k) ←→z/(z-1) ←→ - (- k-1)
零、极 点分布
k k z k k 1
F(z)
K1 e j z
z e j
K1 e j z
z e j
若z> , f(k)=2K1kcos(k+)(k)，… …
(3) F(z)有重极点 推导记忆：
第六章 离散系统的z域分析
3z 例： 2δ(k)+ 3ε(k) ←→ 2 + δ ε z −1
第1-12页 12页
z > 1
青岛科技大学信息科学技术学院
信号与系统 电子教案
6.2
z变换的性质 z变换的性质
二、移位特性
双边z 双边z变换
若： f (k) ←→F (z) ， α<z<β，且有整数 β 且有整数m>0， ， 则： f(k±m) ←→ z±mF(z)， α<z<β ± ， β
2 2
z > a
青岛科技大学信息科学技术学院
信号与系统 电子教案
6.2
z变换的性质 z变换的性质
四、卷积定理
若： f1 (k) ←→F1(z) ， α1<z<β1 β f2 (k) ←→F2(z) ， α2<z<β2 β 则： f1(k) * f2(k) ←→ F1(z)F2(z)， ， 例 收敛域至少为 相交部分 求单边序列 (k+1)akε(k)的z变换，(0<a<1)。 的 变换， 。 变换
三、z域尺度变换(序列乘ak) 域尺度变换(序列乘a
若： f (k) ←→F (z) ， α<z<β，且对整数m>0， β 且对整数 ， 则： ak f(k) ←→ F(z/a)， αa<z<βa ， β 变换。 例：求指数衰减正弦序列 aksin(βk)ε(k) 的z变换。 β 解:
6.1 z 变 换
b k , k < 0 f 2 (k ) = b k ε (−k − 1) = 0, k ≥ 0
解: 反因果序列的 变换为： 反因果序列的z变换为 变换为：
第1-12页 12页
z > 1
青岛科技大学信息科学技术学院
信号与系统 电子教案
6.2
z变换的性质 z变换的性质
二、移位特性
双边z 双边z变换
若： f (k) ←→F (z) ， α<z<β，且有整数 β 且有整数m>0， ， 则： f(k±m) ←→ z±mF(z)， α<z<β ± ， β
2 2
z > a
青岛科技大学信息科学技术学院
信号与系统 电子教案
6.2
z变换的性质 z变换的性质
四、卷积定理
若： f1 (k) ←→F1(z) ， α1<z<β1 β f2 (k) ←→F2(z) ， α2<z<β2 β 则： f1(k) * f2(k) ←→ F1(z)F2(z)， ， 例 收敛域至少为 相交部分 求单边序列 (k+1)akε(k)的z变换，(0<a<1)。 的 变换， 。 变换
三、z域尺度变换(序列乘ak) 域尺度变换(序列乘a
若： f (k) ←→F (z) ， α<z<β，且对整数m>0， β 且对整数 ， 则： ak f(k) ←→ F(z/a)， αa<z<βa ， β 变换。 例：求指数衰减正弦序列 aksin(βk)ε(k) 的z变换。 β 解:
6.1 z 变 换
b k , k < 0 f 2 (k ) = b k ε (−k − 1) = 0, k ≥ 0
解: 反因果序列的 变换为： 反因果序列的z变换为 变换为：
差分方程离散系统的z域分析法稳定性
A/D：模拟信号→数字信号，图中还包括 连续信号→离散信号的采样过程
D/A：数字信号→模拟信号，图中还包括 离散信号→连续信号的保持过程 计算机
r
e
数字
控制器
u(t) 执行
D/A
机构
受控 y(t)对象A/D Nhomakorabea测量
计算机控制系统原理图
4
计算机控制系统的主要特点
修改控制器结构及参数很方便（改变控制程序）； 便于实现各种先进控制，能完成复杂的控制任务； 控制精度高，抗干扰能力强，能有效抑制噪声； 有显示、报警等多种功能。 有利于实现“智能化”、“网络化”、“管控一体 化”、多级分布式控制等；
表示为为便于数学处理将2t2t实际上被保持器抵消了该系数在有保持器的系统中率特性系统的传递函不影响离散信号的频系数相差一个变换的角度看两者只脉冲为时刻的单位幅值时刻的单位幅值脉冲表示ntsst二采样信号的数学表达式变换即得到z变换是离散信号拉氏变换的有理式表达形式个采样时刻的取值的系数为信号在第n10仿真实验
z-n的系数为信号在第n个采样时刻的取值
Z变换只表达了连续函数在采样时刻的特性，不包含采样 时刻之间的信息。
对f(t) 采样后的 f (t) 是唯一的，但 f (t) 所对应的 f(t) 不 唯一； f (t) 与 F(z) 之间的变换是唯一的。
19
S平面与Z平面的对应关系：
根据Z变换定义，有 z eTs
而 uh( t ) 1( t ) - 1( t - T )
U( s ) ,
Uh(
s
)
1
- e-Ts s
零阶保持器实际的传递函数为
Gh (
s
)
离散信号与系统的Z域分析
序列相加减（线性加权）后，所得序列z变换的ROC，有 可能比原序列z变换的ROC大。位移特性常用来分析单边 周期信号，单边周期信号总具有相似的形式。
8 离散信号与系统的 Z 域分析 p 16
例： F(z) = 1/(za) |z| a 求f [k]。 解：
1 F ( z) z 1 1 az
z 例： (3) u[k ] , z 3 z 3
k
类似于傅氏、拉氏变换的尺度变换特性。
1 1 s L f (at ) F ( j ) f (at ) F ( ), a a a a
F
8 离散信号与系统的 Z 域分析 p 18
a 0, a 0
例*：求aksin(0k) u[k] 的z变换及收敛域
1 cos 0 z 1 1 2 z 1 cos 0 z 2 sin 0 z 1 1 2 z 1 cos 0 z 2
五、单边z变换的主要性质
f [k ] F ( z), z R f
f1[k ] F1 ( z), z R f 1
1 2
sin 0 z 1 za 2 2 z 1 cos 0 z 2
8 离散信号与系统的 Z 域分析 p 19
五、单边z变换的主要性质
4. z域微分特性（时域线性加权）
dF ( z ) kf [k ] z dz
Z
Z Rf
m d m d F ( z) Z m m 或写成 : ( z ) F ( z ) k f [k ] ( z ) m dz dz
2 2
8 离散信号与系统的 Z 域分析 p 13
五、单边z变换的主要性质
2. 位移特性(记忆)
因果序列的位移
8 离散信号与系统的 Z 域分析 p 16
例： F(z) = 1/(za) |z| a 求f [k]。 解：
1 F ( z) z 1 1 az
z 例： (3) u[k ] , z 3 z 3
k
类似于傅氏、拉氏变换的尺度变换特性。
1 1 s L f (at ) F ( j ) f (at ) F ( ), a a a a
F
8 离散信号与系统的 Z 域分析 p 18
a 0, a 0
例*：求aksin(0k) u[k] 的z变换及收敛域
1 cos 0 z 1 1 2 z 1 cos 0 z 2 sin 0 z 1 1 2 z 1 cos 0 z 2
五、单边z变换的主要性质
f [k ] F ( z), z R f
f1[k ] F1 ( z), z R f 1
1 2
sin 0 z 1 za 2 2 z 1 cos 0 z 2
8 离散信号与系统的 Z 域分析 p 19
五、单边z变换的主要性质
4. z域微分特性（时域线性加权）
dF ( z ) kf [k ] z dz
Z
Z Rf
m d m d F ( z) Z m m 或写成 : ( z ) F ( z ) k f [k ] ( z ) m dz dz
2 2
8 离散信号与系统的 Z 域分析 p 13
五、单边z变换的主要性质
2. 位移特性(记忆)
因果序列的位移
离散信号与系统的 Z 域分析
第 六 章 离散信号与系统的 Z 域分析引言与线性连续系统的频域分析和复频域分析类似，线性离散系统的频域分析是输入信号分解为基本信号e jΩk 之和,则系统的响应为基本信号的响应之和。
这种方法的数学描述是离散时间傅里叶变换和逆变换。
如果把复指数信号e jΩk 扩展为复指数信号Z k ,Z=re jΩ ,并以Zk 为基本信号， 把输入信号分解为基本信号Z k 之和， 则响应为基本信号Z k 的响应之和。
这种方法的数学描述为Z 变换及其逆变换,这种方法称为离散信号与系统的Z 域分析法.如果把离散信号看成连续时间信号的 抽样值序列，则Z 变换可由拉普拉斯变换引入.因此离散信号与系统的Z 域分析 和连续时间信号与系统的复频域分析有许多相似之处.通过Z 变换，离散时间信 号的卷积运算变成代算,离散时间系统的差分方程变成Z 域的代数方程,因此可 以比较方便的分析系统的响应。
Z 变换从拉普拉斯变换到Z 变换对连续信号f(t)进行理想抽样，即f（t）乘以单位冲击序列δT （t），T 为 抽样间隔，得到抽样信号为f s （t）=f（t）δT （t）= =对fs（t）取双边拉普拉斯变换，得F s （s）=￡[fs（t）]=令z=e sT , 则Fs（s）=F(z) ,得F（z）=因为T为常数,所以通常用f(k)表示f(kT),于是变为F(z)=称为f(k)的双边Z变换,z为复变量。
z和s的关系为:z=e sTs=（1/T）㏑z由复变函数理论,可以得到f(k)= ∮cF(z)z k-1 dz式(7.1-5)称为F(z)的双边Z逆变换(后面讨论).双边Z变换的定义和收敛域§双边 Z 变换的定义对于离散序列f(k)(k=0,±1,±2,┄),函数(z的幂级数)F(z)=称为f(k)的双边Z变换,记为F(z)=Z[f(k)].F(z)又称为f(k)的象函数,f(k)又 称为F(z)的原函数.为了表示方便,f(k)与F(z)之间的对应关系可表示为 f(k) F(z)§双边 Z 变换的收敛域f(k)的双边Z变换为一无穷级数,因此存在级数是否收敛的问题.只有当 (7.1-6)式的级数收敛,F(z)才存在.F(z)存在或级数收敛的充分条件是 ∞在f(k)给定的条件下,式(7.1-6)级数是否收敛取决于z的取值.在z复平面上, 使级数收敛的z取值区域称为F(Z)的收敛域。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验名：离散系统的Z 域分析一、实验目的1、掌握离散序列z 变换的计算方法。
2、掌握离散系统系统函数零极点的计算方法和零极点图的绘制方法，并能根据零极点图分析系统的因果性和稳定性。
3、掌握利用MATLAB 进行z 反变换的计算方法。
二、实验原理与计算方法1、z 变换 离散序列x (n )的z 变换定义为：∑∞-∞=-=n nzn x Z X )()(。
在MATLAB 中可以利用符号表达式计算一个因果序列的z 变换。
其命令格式为：syms n; f=(1/2)^n+(1/3)^n; ztrans(f)2、离散系统的系统函数及因果稳定的系统应满足的条件一个线性移不变离散系统可以用它的单位抽样响应h (n )来表示其输入与输出关系，即y (n )= x (n )* h (n )对该式两边取z 变换，得： Y (z )= X (z )· H (z )则： )()()(z X z Y z H =将H (z )定义为系统函数，它是单位抽样响应h (n )的z 变换，即∑∞-∞=-==n nzn h n h Z z H )()]([)(对于线性移不变系统，若n <0时，h (n )=0，则系统为因果系统；若∞<∑∞-∞=n n h |)(|，则系统稳定。
由于h (n )为因果序列，所以H (z )的收敛域为收敛圆外部区域，因此H (z )的收敛域为收敛圆外部区域时，系统为因果系统。
因为∑∞-∞=-=n nzn h z H )()(，若z =1时H (z )收敛，即∞<=∑∞-∞==n z n h z H |)(||)(1，则系统稳定，即H(z)的收敛域包括单位圆时，系统稳定。
因此因果稳定系统应满足的条件为：1,||<∞≤<ααz ，即系统函数H (z )的所有极点全部落在z 平面的单位圆之内。
3、MATLAB 中系统函数零极点的求法及零极点图的绘制方法MATLAB 中系统函数的零点和极点可以用多项式求根函数roots ()来实现，调用该函数的命令格式为：p=roots(A)。
其中A 为待求根多项式的系数构成的行向量，返回向量p 是包含该多项式所有根位置的列向量。
如：求多项式8143)(2++=z z z A 的根的MATLAB 命令为：A=[1 3/4 1/8]; p=roots(A) 运行结果为： p=也可以用[z,p,k]=tf2zp(B,A)函数求得。
其中z 为由系统的零点构成的向量，p 为由系统的极点构成的向量，k 表示系统的增益；B 、A 分别为系统函数中分子分母多项式的系数向量。
离散系统的系统函数可能有两种形式，一种是分子和分母多项式均按z 的正次幂降幂排列，如12232)(23431+++++=z z z z zz z H ，另一种是分子分母多项式均按z 的负次幂升幂排列，如2112412111)(---+++=zz z z H ，在构造多项式系数向量时，分子和分母多项式系数向量的维数一定要相同，缺项用0补齐。
对于H 1(z )其分子多项式的系数向量应为：B=[0 1 0 2 0]；分母多项式的系数向量应为：A=[1 3 2 2 1]。
对于H 2(z )其分子多项式的系数向量应为：B=[1 1 0]；分母多项式的系数向量应为：A=[1 1/2 1/4]。
绘制系统函数的零极点图可由MATLAB 中的zplane 函数实现。
该函数的调用方法为：zplane(B,A)或者zplane(z,p,k)，其中B ，A ，z ，p ，k 的含义与tf2zp 函数相同。
若调用zplane(B,A)绘图，则首先将系统函数中分子分母多项式变换成按z 的正次幂降幂排列的系数向量，再求零极点。
4、z 反变换的计算方法z 反变换可由部分分式展开法求得。
由于指数序列a n u (n )的z 变换为az z-，因此求反变换时，通常对zz X )(进行展开： kk z z A z z A z z A z z X -+-+-=2211)( 其中),2,1()()(k i z z X z z A i z z i i =-==称为有理函数zz X )(的留数。
分两种情况进行讨论：（1）X (z )的所有极点均为单实极点此时kk z z z A z z zA z z z A z X -+-+-= 2211)(，则X (z )的z 反变换为：∑=⋅+=ki n i i z A A n x 10)()(（2）X (z )有共轭极点设X (z )有一对共轭极点βαj e p ±=2,1，则k k z z z A z z zA p z z r p z z r z X -+-+--=112211)(，其中留数的计算方法与单极点相同，即θj p z e r zz X p z r ||)()(1111=-==，r 2=r 1 *因此，只要求出zz X )(部分分式展开的系数（留数），就可以直接求出X (z )的z 反变换x (n )。
在MATLAB 中可利用函数residue()求解。
令B 和A 分别是zz X )(的分子和分母多项式构成的系数向量，则函数[r,p,k]=residue (B,A)将产生三个向量r 、p 、k ，其中r 为包含z z X )(部分分式展开系数r i (i =1,2,…,N )的列向量，p 为包含z z X )(所有极点的行向量，k为包含z z X )(部分分式展开的多项式项的系数c j (j =1,2,…,M -N )的列向量，若M ≤N ，则k为空阵。
用residue()函数求出zz X )(部分分式展开的系数后，便可根据其极点位置分布情况直接求出X (z )的反变换x (n )。
如：已知23)(22++=z z z z X ，求其z 反变换x (n )。
首先利用residue()函数求出23)(2++=z z zz z X 的部分分式展开的系数和极点，相应的MATLAB 命令为： B=[0 1 0]; A=[1 3 2];[r,p,k]=residue (B,A) 运行结果为： r = 2 -1 p = -2 -1 k = [ ]由以上结果可得：1122)(+-++=z z z X ；即X (z )只有两个单极点，其z 反变换为：[])()1()2(2)(n u n x n n ---=。
已知122)(232-+-+=z z z zz z X ，求其z 反变换x (n )。
利用residue()函数求出1221)(23-+-+=z z z z z z X 的部分分式展开的系数和极点，可得： B=[0 0 1 1]; A=[1 -2 2 -1]; [r,p,k]=residue (B,A)r =+ - p =+ - k =[ ]可见，z z X )(包含一对共轭极点，用abs()和angle()函数即可求出共轭极点的模和相位，相应命令为： p1=abs(p') p1 = a1=angle(p')/pi a1 =0 即共轭极点为：32,1πjep ±=，则12)(33-+--+--=-z z ez z ez z z X jjππ，其z 反变换为：)(23cos 2)(n u n n x ⎥⎦⎤⎢⎣⎡+⎪⎭⎫ ⎝⎛-=π三、实验内容（1）求下列序列的z 变换：2-nu (n )；-(1/2)n u (n )；(1/2)n +(1/3)nu (n )①2-nu (n )的Z 变换程序如下：syms n ; f=(1./2)^(n) ztrans(f)结果为：f = (1/2)^n ans = 2*z/(2*z-1)②-(1/2)nu (n )的Z 变换程序为：syms n ; f=-(1./2)^(n) ztrans(f)结果为：f = -(1/2)^n ans = -2*z/(2*z-1)③(1/2)n+(1/3)nu (n )因为是非因果系统，所以 Z[(1/2)n+(1/3)nu (n )]=()()23nn n n u n z ∞---=-∞+∑=331zz - （2）已知两个离散因果系统的系统函数分别为：142)(2321+-++=z z z z z z H ；2132122112)(-----+++-=z z z z z z H分别求出各系统的零极点，绘制零极点图，分析系统的稳定性；求出各系统单位抽样响应。
① 程序为：A=[1 2 -4 1] B=[0 1 1 0][z,p,k]=tf2zp(B,A)%求零极点[r,p,k]=residue (B,A)%部分分式展开的系数和极点 zplane(B,A)零极点为 ：z = 0 -1 p =零极图为：稳定性：由上图看出收敛域不包括单位圆，即不是稳定系统 系统抽样响应： 由下得r = p = k = []3028.01569.016667.03028.30.4902)(1--+-++=z z z z H 系统抽样响应为：h （n ）=()()[])(3028.01569.06667.0)3028.3(4902.0)(n u n h nn-+-=②程序为：A=[1 1 1/2 0]B=[0 2 -1 1][z,p,k]=tf2zp(B,A)%求零极点[r,p,k]=residue (B,A)%部分分式展开的系数和极点p1=abs(p')a1=angle(p')/pizplane(B,A)零极点为：z =+-p =+-k =2零极图为：稳定性：由上图看出收敛域包括单位圆，所以系统稳定 系统抽样响应： 由r = + - p = + -含一对共轭极点，用abs()和angle()函数即可求出共轭极点的模和相位 p1 =0 a1 =0 即共轭极点为：432,1πj ep ±=，则12)(4343-+--+--=-z z ez z ez z z X j j ππ，其z 反变换为：)(243cos 2)(n u n n x ⎥⎦⎤⎢⎣⎡+⎪⎭⎫ ⎝⎛-=π及为系统的抽样响应。