用z变换解差分方程.
§8.7 用z变换解差分方程
n
( n ≥ 0)
第
例8-7-2
已知系统框图 列出系统的差分方程。 列出系统的差分方程。
n
6 页
x(n)
1 E
+ + +
−3
1 E 1 E
y(n)
(− 2) n ≥ 0 x(n) = , y(0) = y(1) = 0, 0 n<0
求系统的响应 y(n)。 。 解： （1） 列差分方程，从加法器入手 ） 列差分方程，
第
一．应用z变换求解差分方程步骤
一．步骤
(1)对差分方程进行单边 变换（移位性质）； 对差分方程进行单边 变换（移位性质） 对差分方程进行单边z变换 (2)由z变换方程求出响应 由 变换方程求出响应 变换方程求出响应Y(z) ； (3) 求Y(z) 的反变换，得到 的反变换，得到y(n) 。
3 页
0.9y ( −1) z 0.05z2 Y ( z) = + ( z −1)( z − 0.9) z − 0.9
z −1
Y ( z) A A2 1 = + z z −1 z − 0.9
第 5 页
Y ( z) A A2 1 = + z z −1 z − 0.9
A = 0.5 1
A2 = 0.45
z z Y ( z) = 0.5 + 0.45 z −1 z − 0.9
§8.7 用z变换解差分方程
第
序言
2 页
描述离散时间系统的数学模型为差分方程。 描述离散时间系统的数学模型为差分方程。求解差分 方程是我们分析离散时间系统的一个重要途径。 方程是我们分析离散时间系统的一个重要途径。 求解线性时不变离散系统的差分方程有两种方法： 求解线性时不变离散系统的差分方程有两种方法： •时域方法 时域方法 •z变换方法 变换方法 •差分方程经 变换→代数方程； 差分方程经z变换 代数方程； 差分方程经 变换→ •可以将时域卷积→频域（z域）乘积； 可以将时域卷积→ 可以将时域卷积 频域（ 域 乘积； •部分分式分解后将求解过程变为查表； 部分分式分解后将求解过程变为查表； 部分分式分解后将求解过程变为查表 •求解过程自动包含了初始状态（相当于0-的 求解过程自动包含了初始状态（ 求解过程自动包含了初始状态 相当于0 条件）。 条件）。
z变换求解差分方程步骤
z变换求解差分方程步骤嘿，咱今儿就来讲讲这用 z 变换求解差分方程的步骤哈。
这可就像是解开一道神秘的谜题呢！你想想，差分方程就像是一个调皮的小精灵，藏着好多秘密等我们去发现。
而 z 变换呢，就是那把神奇的钥匙啦。
首先呢，得把差分方程给它表示清楚咯，可不能模模糊糊的。
就像你要找东西，总得先知道要找啥样的不是？然后对这个差分方程进行 z 变换，这就好比给它施了个魔法，一下子就变得不一样啦。
在这个过程中啊，你得细心点儿，可别弄错啦。
这就跟走迷宫似的，一步错步步错呀。
接着呢，就会得到一个关于 z 的表达式，这可就是我们前进的线索呢。
然后呢，咱得把这个表达式给它化简化简，把那些复杂的东西都去掉，就像给苹果削皮一样，让它露出最精华的部分。
这时候可就考验咱的本事啦，得有耐心，还得有那么点儿小技巧。
再接下来呀，就得求解啦！这就像是终于找到了宝藏的位置，要把它挖出来一样。
把 z 的值求出来，这可不容易呢，但咱不能怕呀，要勇往直前！等求出了 z 的值，可别以为就大功告成咯。
还得把它变回原来的世界，也就是反变换回去。
这就像是把变了形的东西再变回来，可神奇啦。
哎呀，你说这过程是不是挺有意思的？就好像是一场冒险，每一步都充满了挑战和惊喜。
你要是能熟练掌握这 z 变换求解差分方程的步骤，那可就厉害咯，就像是拥有了超能力一样！你想想，以后遇到那些复杂的差分方程，别人都抓耳挠腮不知道咋办的时候，你就能轻松搞定，那多牛呀！这就好比别人还在走路，你都开上小汽车啦，一下子就把他们甩在后面啦。
所以呀，可得好好学这 z 变换求解差分方程的步骤哦，别偷懒，多练练，肯定能掌握得牢牢的。
到时候，不管啥样的难题都难不倒你啦！这多棒呀，是不是？。
Z变换及差分方程的求解
Z变换及差分⽅程的求解第⼆讲离散时间动态经济系统运动分析及稳定性分析2.1离散时间函数与Z变换⽬的要求：通过本节的学习使学⽣掌握离散时间函数及Z变换的概念，会使⽤Z变换的性质解决问题，掌握差分⽅程及离散时间系统的运动分析⽅法。
教学内容：我们经常会遇到利⽤离散时间函数表⽰的差分⽅程或差分⽅程组，这在经济管理中经常遇到。
现介绍离散时间函数，差分⽅程后⾯介绍。
⼀、离散时间函数例1 ⼈⼝离散时间函数设全国⼈⼝普查每年进⾏⼀次。
每年的7⽉1⽇凌晨零点的⼈⼝数代表该年的⼈⼝数。
我们以t=0 代表1990年7⽉1⽇凌晨的这个时刻，那么t=1,2,3,……分别表⽰1991年、1992年、1993年等各年度7⽉1⽇凌晨零点。
各年度普查的实际⼈⼝数如下表所⽰中国实际⼈⼝数据（亿⼈）x(0)=11.4333, x(1)=11.5823, x(2)=11.7171,x(3)=11.8517, x(4)=11.9850, x(5)=12.1121,x(6)=12.2389, x(7)=12.3626,……由于在离散时间离取值，故称之为离散时间函数例2 国民⽣产总值GNP(gross national product)离散时间函数。
则，GNP(t)表⽰第t年的GNP数值。
GNP(O)=33560.5, GNP(1)=46670.0, GNP(2)=57494.9,……例3 企业⽉产量离散时间函数。
表为电视机⼯⼚⽣产⽉报表（万台）则，Y(0)=1.5, Y(1)=2, Y(2)=1.8,……可以看出，经济管理实践中基本上采⽤离散时间函数来表达各种变量的变化，并该函数没有解析表达式，只有图象、列表表达式。
其⾃变量为离散时间。
⼆、Z 变换及其逆变换导⾔：Z 变换是怎么发明出来的？⽜顿、莱布尼兹等发明了微积分，之后发明了常系数微分⽅程及⽅程组。
在求解⽅程时总结经验，简化计算，如⽤符号s 表⽰微分运算s=d/dt,即s 〃f(t)=df(t)/dt 。
差分方程及其Z变换法求解
= b0 r[(k + m)T ] + b1r[(k + m − 1)T ] + .......bm−1r[(k + 1)T ] + bm r (kT )
zX 1 ( z ) − zx1 (0) = X 2 ( z )
x2(kT)
z −1
x1(kT) z −1 x2(z) y[(k+1)T] KT
-
x1(0) 1 x1(z)
例2：画出例2所示离散系统的模拟图
y[(k + 1)T ] = -( KT -1) y (kT ) + KTr (kT ) r(kT)+ 1)T ] + ( KT -1) y (kT ) = KTr (kT ) y (k + 1) + ( K -1) y (k ) = Kr (k )
KT-1
三、差分方程的解
差分方程的求解：迭代法、z变换法。 迭代法：将原系统的差分方程化为如下形式：
y[(k + n)T ] = −a1 y[(k + n − 1)T ] − ...... − an −1 y[(k + 1)T ] − an y[kT ] + b0 r[( k + m)T ] + b1r[(k + m − 1)T ] + .......bm −1r[( k + 1)T ] + bm r (kT )
y (kT ) = 0.446 + 1.429(-0.4) k -1.875(-0.6) k
差分方程的Z变换解
F=ztrans(f) 式中，f为时间函数的符号表达式，F为Z变换式，也是符号表达式。
反变换的调用格式为 f=iztrans(F)
式中，F为Z变换式的符号表达式，f为时间函数，是符号形式。
为了改善公式的可读性，MATLAB提供了pretty函数，调用格式为 Pretty(f)
全响应
n
n
-6 (-1) + 7/3 (-2) + 2/3
10
实验内容 1
求下列序列的变 换，并注明收敛域。
(a) (b) (c) (d)
11
实验内容 2
求下列
(a)
的逆变换 。
(b)
(c)
(d容 3
用单边 变换解下列各差分方程。
(a)
，
(b)
，
，
13
实验内容 4*
已知
，收敛域
，求 。
解 部分分式展开式为
或
8
计算示例 3
描述某离散系统的差分方程为
激励信号
，若初始条件
，
试分别求其零输入响应
、 零状态响应
和
全响应。
9
计算示例 3
程序运行后在命令窗口显示的结果：
>> 零状态响应
n
n
-(-1) + 1/3 (-2) + 2/3
零输入响应
n
n
-5 (-1) + 2 (-2)
式中，f为符号表达式。
3
实验原理与说明
2、求反Z变换的部分分式法
若
为有理式，则可表达为
MATLAB提供了一个对
信号与系统5-2差分方程的Z变换解课件
电信学院
1
前向差分方程
查公式
考虑二阶系统：
y(k 2) a1y(k 1) a0 y(k) b2 f (k 2) b1 f (k 1) b0 f (k)
初始值：yzi (0), yzi (1)
两边取Z变换有：
(z2 a1z a0 )Y (z) yzi (0)z2 yzi (1)z a1yzi (0)z (b2z2 b1z b0 )F(z)
1
(z
1)( z
2)
z
1
3
z
1
z
1
3
z
2
全响应
yzs (k )
[2 3
(1)k
1 3
(2)k
] (k)
y(k)
yzi
(k)
yzs (k )
[
2 3
6(1)k
2 3
(2)k
]
(k)
电信学院 返回
8
例 5.12 解 法 二
y(k 2) 3y(k 1) 2y(k) f (k 1) 3 f (k) yzi(1)=1, yzi(2)=3
F(z)
Y (z) Yzi (z) Yzs (z) 零输入响应
零状态响应
电信学院
3
系统函数
定义
H
(z)
零状态响应的z变换 激励信号的z变换
Yzs (z) F(z)
二阶系统零状态响应
Yzs (z)
b2z2 b1z b0 z2 a1z a0
F(z)
H (z)F (z)
对n阶LTI系统的系统函数
(b2z2 b1z b0 )F(z) b2 f (0)z2 b2 f (1)z b1 f (0)z
令：M (z) [ y(0) b2 f (0)]z2 [ y(1) a1y(0) b2 f (1) b1 f (0)]z
利用z变换解差分方程
于是
Y(z) =
br z−r ∑ ak z−k ∑
k＝ 0 M r＝ 0 N
M
X(z)
令
H(z) =
∑b z
r r＝ 0 N k＝ 0
−r
ak z−k ∑
则
Y(z) = X (z)H(z)
−1
此时对应的序列为 y(n) = F [X(z)H(z)]
差分方程为 例：若描述离散系统的 1 1 y(n) + y(n −1) − y(n − 2) = x(n) 2 2 x(n) = 2n u(n) , y( 已知激励 初始状态 −1) =1, y(−2) = 0, 求系统的零输入响应、 零状态响应和全响应。 求系统的零输入响应、 零状态响应和全响应。
ak z−k [Y(z) = ∑br z−r [X(z) + ∑x(m)z−m] ∑
k＝ 0 r＝ 0 m=−r N M −1
如果激励x(n)为因果序列， 如果激励x(n)为因果序列，上式可以写成 x(n)为因果序列
ak z−k [Y(z) = ∑br z−r X(z) ∑
k＝ 0 r＝ 0 N M
8.5节已经给出利用 节已经给出利用z 在8.5节已经给出利用z变换解差分方程的简 单实例，本节给出一般规律。 单实例，本节给出一般规律。这种方法的原 理是基于z变换的线性和位移性， 理是基于z变换的线性和位移性，把差分方程 转化为代数方程，从而使求解过程简化。 转化为代数方程，从而使求解过程简化。
k＝ 0 l =−k r＝ 0 m=−r −1
若激励x(n)=0,即系统处于零输入状态，此时 若激励x(n)=0,即系统处于零输入状态， x(n)=0,即系统处于零输入状态 差分方程（ 差分方程（1）成为齐次方程∑a y(n −源自) =0k=0 kN
差分方程及其Z变换法求解
例1：右图所示的一阶系统描述它的微分方程为
y(t ) Ke(t ) K (r (t ) y(t ))
y(t ) Ky(t ) Kr (t )
用一阶前向差分方程近似：
(1)
r( t ) e( t ) -
K
1/s
y( t )
y (k 1)T y (kT ) dy y (t ) lim dt T 0 T
由图：x1 (k 1)T x2 (kT )
zX 1 ( z ) zx1 (0) X 2 ( z )
x2(kT)
z
1
x1(kT)
z 1
x1(0) 1
x1 ( z)
x2(z) y[(k+1)T]
例2：画出例2所示离散系统的模拟图
y[(k 1)T ] -( KT -1) y(kT ) + KTr (kT ) r(kT)
y (k 1)T y (kT ) T
(T 很小)
(2)
式中：T为采样周期，（2）代入（1）得：
y (k 1)T (KT 1) y(kT ) KTr(kT )
y(k 1) ( K 1) y(k ) Kr (k )
(3)
二、离散系统差分方程的模拟图
连续系统采用积分器s-1作为模拟连续系统微分方程的主要器件； 与此相对应，在离散系统中，采用单位延迟器z-1。 单位延迟器：把输入信号延迟一个采样周期T秒或延迟1拍。
再利用初始条件，逐次迭代得到各采样时刻的值。
特点：适用于计算机处理求解。 例3：用迭代法解二阶差分方程 y(k+2) +3y(k+1)+2y(k)=1(k)
利用初始条件 y(0)=0, y(1)=1，则有： y(k+2) =-3y(k+1) -2y(k)+1(k) y(2) =-3y(1) -2y(0)+1(0)= -3*1-2*0+1= -2
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
分母高于分子的次数是响应的起点
Y
z
z
2z
1z
22
从n 2开始yn有不为零的值。
三．差分方程解的验证
原方程迭代出 y0, y1, y2两种迭代结果相同, 解的表达式迭代出y0, y1, y2解答是正确的
一．应用z变换求解差分方程步骤
一．步骤
(1)对差分方程进行单边z变换（移位性质）； (2)由z变换方程求出响应Y(z) ； (3) 求Y(z) 的反变换，得到y(n) 。
二．差分方程响应y(n)的起始点确定
全响应y(n)根据输入信号加上的时刻定
对因果系统y(n)不可能出现在x(n)之前
观察Y(z)分子分母的幂次
§8.7 用z变换解差分方程
序言
描述离散时间系统的数学模型为差分方程。求解差分方 程是我们分析离散时间系统的一个重要途径。
求解线性时不变离散系统的差分方程有两种方法： •时域方法——第七章中介绍，烦琐 •z变换方法 •差分方程经z变换→代数方程； •可以将时域卷积→频域（z域）乘积； •部分分式分解后将求解过程变为查表； •求解过程自动包含了初始状态（相当于0-的 条件）。