z变换分析法

合集下载

第五次实验心得体会

第五次实验心得体会

心得体会今天我们做的实验是离散信号与系统的Z 变换分析, Z 变换分析法是分析离散时间信号与系统的重要手段, 实验前我书上和资料上了解到Z 变换它是由拉氏变换而来的, 属于一种线性坐标变换, 它将差分方程化为代数方程, 是分析采样系统的主要数学工具。

在离散系统分析中为简化运算而建立的对函数序列的数学变换, 其作用与拉普拉斯变换在连续系统分析中的作用很相似。

在采样控制理论中,Z 变换是主要的数学工具。

Z 变换还在时间序列分析、数据平滑、数字滤波等领域有广泛的应用。

在MATLAB 语言中有专门对信号进行正反Z 变换的函数ztrans( ) 和itrans( )。

离散信号f(k)的Z 变换定义为:()()k k F z f k z ∞-=-∞=∑反Z 变换的定义为:11()()2k f k F z z dz j π-=⎰(1)求离散序列的Z 变换:1122()()cos()()k k f k k πε=程序:syms k zf=0.5^k*cos(k*pi./2);Fz=ztrans(f)运行结果: Fz =4*z^2/(4*z^2+1)(2)离散序列:3()()(5)f k k k εε=--程序: syms k z f=('Heaviside(k)-Heaviside(k-5)')Fz=ztrans(f)运行结果:f =Heaviside(k)-Heaviside(k-5)(3)但在离散序列:[]4()(1)()(5)f k k k k k εε=---程序: syms k z f=k*(k-1)*('Heaviside(k)-Heaviside(k-5)')Fz=ztrans(f)运行结果: Fz =2/z^4*(z^2+3*z+6)在两个离散序列出现了不同的结果, 前者直接输出原来的函数, 猜想是不是因为后者系数K (K-1)有关。

执行下列程序: syms k zf=k*(k-1)Fz=ztrans(f)运行结果: Fz =z*(1+z)/(z-1)^3-z/(z-1)^2(4)而3()()(5)f k k k εε=--的z 变换为: Fz=(z/z-1)-(z^(-5)*z/z-1)=(z-z^(-4))/z-1 和用MATLAB 仿真的f =Heaviside(k)-Heaviside(k-5)显然不符。

离散系统Z变换分析法02

离散系统Z变换分析法02

3.闭环 Z 传递函数的结构图1
闭环 Z 传递函数的结构图2
2.5.4 过渡过程特性
与连续系统用传递函数分析过渡过程类 似,可以用 Z传递函数来分析离散系统的过 渡过程特性。 • 分析离散系统的过渡过程特性的步骤: • • 1)Y(Z)=GC(Z)R(Z)


2)由Y(Z)求出y(kT)
例题12 例题12
2. 开环 Z 传递函数 • 线件离散系统的开环 Z传递函数 跟连续系统的开环传递函数具有类似 的特性。
串联环节的Z传递函数
例题9
z az , G2 ( z ) = , 设图2 − 10 a)中G1 ( z ) = ( − aT z −1 z −e 试求开环Z传递函数G ( z )。 z az 解:G ( z ) = G1 ( z )G2 ( z ) = z − 1 z − e − aT az 2 = ( z − 1)( z − e − aT )
(1)离散系统稳定的充要条件(时域) 设:系统差分方程
c(k ) + a1c(k − 1) + a2 c(k − 2) + L + an c(k − n) = b0 r (k ) + b1r (k − 1) + L + b0 r (k − m)
系统齐次方程
ቤተ መጻሕፍቲ ባይዱ
c(k ) + a1c(k − 1) + a2 c(k − 2) + L + an c(k − n) = 0
−1 −1 −Ts
1 1 )( − )] s s+a
= 1(t ) − 1(t − T ) − e − at + e − a ( t −T ) 对y (t )采样,离散化后,得 y (kT ) = 1(kT ) − 1(kT − T ) − e −akT + e − a ( kT −T ) 则 HG ( z ) = Z [ y (kT )] z 1 z 1 1 − e − aT = − − − = − aT − aT z −1 z −1 z − e z−e z − e −aT

基于Matlab语言的线性离散系统的Z变换分析法

基于Matlab语言的线性离散系统的Z变换分析法

基于Matlab语言的线性离散系统的Z变换分析法实验一基于Matlab语言的线性离散系统的Z变换分析法班级: 姓名: 学号: 日期:一、实验目的:1、学习并掌握Matlab语言离散时间系统模型建立方法;2.学习离散传递函数的留数分析与编程实现的方法;3.学习并掌握脉冲与阶跃的编程方法;4.理解与分析离散传递函数不同极点的时间响应特点。

二、实验工具:1MATLAB软件(6、5以上版本);2每人计算机一台。

三、实验内容:1在Matlab语言平台上,通过给定的离散时间系统差分方程,理解课程中Z变换定义,掌握信号与线性系统模型之间Z传递函数的几种形式表示方法;2学习语言编程中的Z变换传递函数如何计算与显示相应的离散点序列的操作与实现的方法,深刻理解课程中Z变换的逆变换;3通过编程,掌握传递函数的极点与留数的计算方法,加深理解G(z)/z的分式方法实现过程;4通过系统的脉冲响应编程实现,理解输出响应的离散点序列的本质,即逆变换的实现过程;5通过编程分析,理解系统单位阶跃响应的Z变换就是系统的传递函数与单位阶跃函数Z变换,并完成响应的脉冲离散序列点的计算;6通过程序设计,理解课程中的不同的传递函数极点对系统动态行为的影响,如单独极点、复极点对响应的影响。

四、实验步骤:(一)传递函数的零极点程序: 结果:numg=[0、1 0、03 -0、07];deng=[1 -2、7 2、42 -0、72];g=tf(numg,deng,-1)get(g);[nn dd]=tfdata(g,'v')[zz,pp,kk]=zpkdata(g,'v')hold onpzmap(g), hold offaxis equal(二)留数法程序:numg=[2 -2、2 0、65];deng=[1 -0、6728 0、0463 0、4860];[rGoz, pGoz,other]=residue(numg,[deng 0])G=tf(numg,deng,-1)impulse(G)[y,k]=impulse(G);stem(k,y,'filled');impulse(G)结果:rGoz = 0、4905 + 0、0122i0、4905 - 0、0122i-2、31851、3374pGoz = 0、6364 + 0、6364i0、6364 - 0、6364i-0、6000other = []Transfer function:2 z^2 - 2、2 z + 0、65-----------------------------------z^3 - 0、6728 z^2 + 0、0463 z + 0、486Sampling time: unspecified(三)不同位置的根对系统的影响1)2个共轭极点(左圆内)+1实极点(圆内)P1 =0、6364 + 0、6364iP2=0、6364 - 0、6364iP3=-0、6000程序: 结果:zz3=[-0、2 0、4];pp3=[-0、6 0、6364+0、6364i 0、6364-0、6364i];kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,50);stem(k,y,'filled'),grid2)2个共轭极点(右圆内)+1实极点(圆内)P1= -0、8592 P2= -0、0932 + 0、4558i P3= -0、0932 - 0、4558i 程序: 结果:zz3=[-0、2 0、4];pp3=[-0、8592 -0、0932+0、4558i -0、0932-0、4558i]; kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,50);stem(k,y,'filled'),grid3)2个共轭极点(圆上)+1实极点(圆内)p1=0、6+0、8i p2=0、6-0、8i p3=-0、6程序: 结果:zz3=[-0、2 0、4];pp3=[-0、8592 -0、6+0、8i -0、6-0、8i];kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,100);stem(k,y,'filled'),grid4、2个共轭极点(虚轴上)+1实极点(圆内)p1=i p2= -i p3= -0、6程序: 结果:zz3=[-0、2 0、4];pp3=[-0、6 i -i];kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,100);stem(k,y,'filled'),grid5、2个实极点(圆内)+1个实极点(圆外)p1=2 p2=0、8 p3=-0、6程序: 结果:zz3=[-0、2 0、4];pp3=[2 0、8 -0、6];kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,100);stem(k,y,'filled'),grid6、2个实极点(圆内)+1个实极点(圆上)p1=1 p2=0、8 p3=-0、6程序: 结果:zz3=[-0、2 0、4];pp3=[1 0、8 -0、6];kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,100);stem(k,y,'filled'),gridp1=1 p2=-0、8 p3=-0、6程序: 结果:zz3=[-0、2 0、4];pp3=[1 0、8 -0、6];kk3=2;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);[y,k]=impulse(eg3,100);stem(k,y,'filled'),grid五、实验报告要求1、根据实验结果,分析离散传递函数不同极点的时间响应特点2、通过程序设计,分析不同的传递函数极点如:单极点、复极点、重根极点对系统动态行为的影响3、分析留数法的意义,根据系统的阶跃响应判别系统的稳定性4、对Z变换的进一步思考六、实验结果:1、根据实验结果,分析离散传递函数不同极点的时间响应特点。

信号分析第六章第一节z变换及收敛域

信号分析第六章第一节z变换及收敛域

X
15
4 斜变序列 x(k)k(k)


Z[k(k)]
kzk
z
k zk1
z [d (zk)]
k0
k0
k0 dz
z d[ d zk0
zk]z
z
z1
z (z1)2
k(k) z
(z 1)2
RO:C z1
kak1(k)
z (za)2
ROC: z a
X

16

k2(k) k2zk
z(z1)
离散系统的Z变换分析
连续系统的拉氏变换分析
X
4

第一节 Z 变 换

一.Z变换的提出—由拉氏变换引出
连续信号 等间隔采样 抽样信号
x s(t) x (t)T (t) x (t) (t k)T x (k)T ( t k)T
k 0
k 0
单边拉氏变换
X s (s)
0
x(kT) (t kT)est dt
★反因果序列的ROC为 z R的2 圆内区域;
即X(z) 最小的模值极点为半径的圆内区域 注意:收敛域是否包含z=0需判断. ★双边序列的因果和反因果序列的收敛域存在公共域,
ROC为R1 z R2圆环状,不存在公共区域z变换不存在.
★ ROC内不包含任何极点(以极点为边界);
★有限长序列的ROC为整个 z 平面 0 z
k 0
x(kT) (t kT)est dt 0 k 0
x(kT)eskT 引入连续复变z 量 esT
k 0
取 T1 X S(s) x(k)Z kX (Z ) k 0
X
5
说明:
第 页

Z域变换分析方法

Z域变换分析方法
[1 0.7 z 0.1z ]Y ( z) 0.7 y(1) 0.1z y(1) 0.1y(2)
1 2 1
第8章 Z变换
(2 z 2.6)z 代入初始条件,整理得 : Y ( z ) 2 z 0.7 z 0.1 Y ( z) (2 z 2.6) 12 10 z ( z 0.2)(z 0.5) ( z 0.5) ( z 0.2)
例8-10: 已知某离散LTI系统的单位阶跃响应为:
s[n] (2 3 5 10)u[n]
n n
(1)求系统单位抽样响应 (2)求此二阶差分方程
解: ( 1)
h[n] s[n] s[n 1] 1 n 12 n ( 2 5 )u[n] 11.1 [n] 2 5稳定系统全部极点就一定是位于单位圆内的呢?
第8章 Z变换
三、由极点分布决定系统稳定性 系统稳定的充要条件是单位样值响应绝对可和。即:
n
h( n )

因果稳定系统的充要条件为 :h(n)是单边的而且是有 界的。即: 因果
稳定
h(n) h(n)u (n) 非因果也 可以稳定 h( n) a<1 n
一、系统函数的求取 定义一:系统单位样值响应h[n]的Z变换
激励与单位样值响应的卷积为系统零状态响应
y[n] x[n] h[n]
由卷积定理
Y ( z) X ( z)H ( z)
Y ( z) H ( z) X ( z)
H ( z ) h[n]z
n 0

n
第8章 Z变换
定义二:系统零状态响应的Z变换与输入的Z变换之比 若x(n)是因果序列, 则在系统零状态下:

第三章 Z变换

第三章 Z变换

1
n
第一项为右边序列(因果)其收敛域为: z
第二项为左边序列,其收敛域为: 当Rx-<Rx+时,其收敛域为
Rx
0 z Rx
Rx z Rx
j Im[ z ]
Re[ z ]
Rx Rx
例:求序列
解:这相当
x(n) (n) 的Z变换及收敛域。
n1 n2 0 时的有限长序列,
留数的求法:
1、当Zr为一阶极点时的留数:
Re s[ X ( z ) z n1 ]Z Z r [( z zr ) X ( z) z n1 ]z zr
2、当Zr为l阶(多重)极点时的留数:
Re s[ X ( z ) z
l 1
n 1
] z zr
1 d l n 1 [( z z r ) X ( z ) z ] z zr l 1 (l 1)! dz
X ( z) 4 A ]z 2 1 [( z 2) z 3 X ( z) 1 A2 [( z 0.5) ] z 0.5 z 3 4 z 1 z X ( z) 3 z2 3 z 0.5
又 z 2, 得 4 n 1 n 2 (0.5) , n 0 x ( n) 3 3 ,n 0 0
Ak Re s[( z zk ) z ]z zk r k 1 d r x( z ) Ck r k [( z zi ) (r k )! dz z z z , k 1,2 i
r
分别求出各部分分式的z反变换(可查 P43 表3-1),然后相加即得X(z)的z反变换。
j Im[ z ]
Re[ z ]
收敛域

正态分布z变换

正态分布z变换

正态分布z变换
正态分布z变换是由数学家斯特凡黑格尔于1809年发明的一种
统计分析方法,正态分布z变换通常用来描述数据的分布,它能够帮助研究者估计个体特性,以及解释结果的实际意义。

正态分布z变换被广泛地应用于各种领域,包括金融、商业、政策分析等。

例如,在金融领域,该变换可用来确定给定变量的相关性,以及分析投资市场中的投资组合。

在商业领域,它可以用来识别客户群体中的潜在多样性,以及识别不同产品组合之间的关系。

在政策分析领域,它可以用来分析政策的影响,以及研究新政策的有效性。

正态分布z变换的基本原理是使用均值为零和标准差为一的正
态分布标准化随机变量,从而解决了处理非正态分布数据的问题。

与传统的标准化方法(例如均值除以标准偏差)相比,正态分布z变换更具有效率。

当处理的数据服从正态分布时,正态分布z变换的精度较高。

正态分布z变换通过对原始数据做几何均值和平方和变换,并将所得结果舍入到最接近的整数,来将原始数据变换为标准的正太分布数据。

标准正太分布数据可用来计算相关性和分析数据。

正态分布z变换并不是可以应用于所有数据的情况。

有些情况下,数据分布不满足正态分布,因此必须采用更加复杂的变换。

正态分布z变换也不能应用于有明显偏离正态分布的数据。

例如,在某些情况下,采用双峰分布或者无峰分布的数据,就不能采用正态分布z变换。

正态分布z变换具有可靠性,可以用来提供更准确的分析结果和
更有效的数据分析。

因此,可以说正态分布z变换在当今许多领域中起着至关重要的作用。

研究者可以利用正态分布z变换来更好地分析复杂的数据,从而获得更准确的结果。

第七章离散系统的Z变换分析方法

第七章离散系统的Z变换分析方法

第七章离散系统的Z变换分析⽅法第七章线性离散系统与Z 变换第⼀节概述离散系统(采样数字系统),与连续系统的根本区别在于所处理的信号是离散型的。

在离散控制系统中,认为系统变量仅是在离散的时刻上才发⽣变化,⽽在两个相邻时刻之间是不发⽣变化的。

离散信号的时间函数如图7-1所⽰。

图7-1 离散的时间函数在离散控制系统中最常⽤的计算机控制系统,其原理图7-2如所⽰。

图7-2 计算机控制系统原理图◆线性连续系统的动态特性可以由微分⽅程描述,分析线性定常连续系统采⽤拉⽒变换;◆线性离散系统的动态特性可以⽤线性差分⽅程描述,分析线性定常离散系统采⽤Z 变换法。

Z 变换是分析单输⼊单输出、线性定常离散系统的有⼒⼯具。

第⼆节 Z 变换Z 变换是由拉⽒变换引出的,可以把Z 变换看成拉⽒变换的⼀种变形。

⼀、采样函数的拉⽒变换设连续时间函数()x t 可以进⾏拉普拉斯变换,其拉⽒变换为()X s 。

连续时间函数 ()x t 经采样周期为 0T 的采样器采样后,变成离散信号*()x t+-++-+-+=)()()2()2()()()()0()(000000*nT t nT x T t T x T t T x t x t x δδδδ=()()n x nT t nT δ∞=-∑ (7-1)对上式进⾏拉普拉斯变换,⼜snT e nT t L 0]([0-=-δ可得0**0000000()[()][()()]()[()]()n nT sn n X s L x t L x nT t nT x nT L t nT x nT e δδ∞=∞∞-====-=-=∑∑∑ (7-2)⼆、采样函数的Z 变换在式(7-2)中,由于s 在指数⾥,给运算带来许多困难。

为此引进新的变量0T s z e =,则式7-2变形为∑∞=-=00)()(n nz nT x z X (7-3)称()X z 为离散时间函数 *()x t 的Z 变换,记为 *[()]()Z x t X z =或者[()]()Z x n T X z =。

实验一 基于matlab语言的线性离散系统的z变换分析法1(1)

实验一 基于matlab语言的线性离散系统的z变换分析法1(1)

实验一基于MATLAB语言的线性离散系统的Z变换分析法一、实验目的1. 学习并掌握 Matlab 语言离散时间系统模型建立方法;2.学习离散传递函数的留数分析与编程实现的方法;3.学习并掌握脉冲和阶跃响应的编程方法;4.理解与分析离散传递函数不同极点的时间响应特点。

二、实验工具1. MATLAB 软件(6.5 以上版本);2. 每人计算机一台。

三、实验内容1. 在Matlab语言平台上,通过给定的离散时间系统差分方程,理解课程中Z变换定义,掌握信号与线性系统模型之间Z传递函数的几种形式表示方法;2. 学习语言编程中的Z变换传递函数如何计算与显示相应的离散点序列的操作与实现的方法,深刻理解课程中Z变换的逆变换;3. 通过编程,掌握传递函数的极点与留数的计算方法,加深理解G(z)/z 的分式方法实现过程;4. 通过系统的脉冲响应编程实现,理解输出响应的离散点序列的本质,即逆变换的实现过程;5. 通过编程分析,理解系统的Z传递函数等于单位脉冲响应的Z变换,并完成响应的脉冲离散序列点的计算;6. 通过程序设计,理解课程中脉冲传递函数极点对系统动态行为的影响,如单独极点、复极点对响应的影响。

四、实验步骤1.创建系统How to create digital system g Four examples are as follows:numg=[0.1 0.03 -0.07];deng=[1 -2.7 2.42 -0.72];g=tf(numg,deng,-1)get(g);[nn dd]=tfdata(g,'v')[zz,pp,kk]=zpkdata(g,'v')Unite circle region with distrbuting zeros points and poles points hold onpzmap(g), hold offaxis equal运行结果:2.转换为零极点标准形式Convert from tf(z-function) to zpk(z-function) Part C exercise form gg=zpk(g)[zz,pp,kk tts]=zpkdata(gg,'v')[z,p k,ts]=zpkdata(g,'v')运行结果:3.四个例子Four examples are as follows:Part A exerciseeg1mun=[1.25 -1.25,0.30];eg1den=[1 -1.05 0.80 -0.10];eg1=tf(eg1mun,eg1den,-1);eg1zpk=zpk(eg1);[zz1,pp1,kk1,tts1]=zpkdata(eg1zpk,'v');Part B exerciseeg2mun=[0.84 -0.062 -0.156 0.058];eg2den=[1 -1.03 0.22 0.094 0.05];eg2=tf(eg2mun,eg2den,-1);eg2zpk=zpk(eg2);[zz2,pp2,kk2,tts2]=zpkdata(eg2zpk,'v');Part C exercisezz3=[-0.2 0.4];pp3=[0.6 0.5+0.75i 0.5-0.75i 0.3];kk3=150;tts3=-1;eg3zpk=zpk(zz3,pp3,kk3,tts3);eg3=tf(eg3zpk);Part D exercisezz4=[-0.3 0.4+0.2i 0.4-0.2i];pp4=[-0.6 -0.3,0.5 0.6];kk4=5;tts4=-1;eg4zpk=zpk(zz4,pp4,kk4,tts4);eg4=tf(eg4zpk);4.留数法Residue method and impluse response numg=[2 -2.2 0.65];deng=[1 -0.6728 0.0463 0.4860]; [rGoz, pGoz,other]=residue(numg,[deng 0]) [mag_pGoz,theta_pGoz] =xy2p(pGoz)[mag-rGoz,theta-rGoz]=xy2p(rGoz)G=tf(numg,deng,-1)impulse(G)[y,k]=impulse(G);stem(k,y,'filled');impulse(G)运行结果:5.复杂极点响应When transfer function is G(Z) with complex ,t=t*ts;pole of z=e^(+-j*30*pi/3) and z=-0.5,as well as its gain value is unit step signal,its collecting cycle is 0.5 second,how to analyze its response.gcfts=0.3;num=[1 0.5];den=conv([1 -exp(i*pi/3)],[1 -exp(-i*pi/3)]);g1=tf(num,den,ts)[y,k]=impulse(g1,20);stem(k,y,'filled'),grid运行结果:6.重极点响应How to analyze response with repeating poles dtime=[0:90];y(k+2)-1.8y(k+1)+0.81y(k)=3u(k+1)-1.2u(k) yi=impulse(gstep,dtime)gcfnum=[3 -1.2];den=[1 -1.8 0.81];[rGoz, pGoz,other]=residue(num,[den 0])t=0:60;y=rGoz(2,1).*(t.*(pGoz(2,1).^(t-1)))+rGoz(1,1).*(pGoz(1,1).^(t)) y1=zeros(1,61);y1(1,1)=rGoz(3,1);y=y+y1;t=ts*t;stem(t,y,'filled'),gridSpecial example about difference real pole tosystem response[rGoz,pGoz,other]=residue(num,[den,0])num1=[rGoz(1) 0];den1=[1 -pGoz(1)]gg1=tf(num1,den1,ts)[y,t]=impulse(gg1,50)stem(t,y,'filled'),grid运行结果:7.阶跃响应numg=[2 -2.2 0.56];deng=[1 -0.6728 0.0463 0.4860];g=tf(numg,deng,1);numgstep=[numg 0];dengstep=conv(deng,[1 -1]);gstep=tf(numgstep,dengstep,1)dtime=[0:90];yi=impulse(gstep,dtime)subplot(2,1,1)stem(dtime,yi,'filled')ys=step(g,dtime);subplot(2,1,2)stem(dtime,ys,'filled')dcgain(g)ys_ss=ys(end)ys_ss=ys(max(dtime))运行结果:Example 1: Analysis of subsection input function subplot(1,1,1)num=[2 -2.2 0.56];den=[1 -0.6728 0.0463 0.4860];ts=0.2;g=tf(num,den,ts);dtime=[0:ts:8]';u=2.0*ones(size(dtime));ii=find(dtime>=2.0); u(ii)=0.5;y=lsim(g,u,dtime);stem(dtime,y,'filled'),gridhold onplot(dtime,u,'o')hold offtext(2.3,-1.8,'output')text(1.6,2.3,'input')运行结果:五、实验思考1、根据实验结果,分析离散传递函数不同极点的时间响应特点。

第五章 Z域分析

第五章 Z域分析
m

m 1
x(k ) z
max( R 11 , R 21 ) z
2. 位移性
a. 双边Z变换

x(n) X ( z )
m
x(n m ) z
X ( z ), m 为整数
收敛域:1)不包括
0,
处,收敛域不变 处,需重新判断
2) 包括 0,
证明: z [ x ( n m )] 令k=n+m
z [ x ( n m )] z
x ( n )u ( n ) x ( z )
m

x ( n m )u ( n ) z [ X ( z )

m 1
x(k ) z
k
]
k 0
证明:
z [ x ( n m ) u ( n )]


x(n m ) z
n
n0
令k=n+m 则:
z [ x ( n m ) u ( n )] z [ X ( z )

a. X(z)/z 有N个单极点
则:
Z1 Z N
X z

N
Ai z z zi
i0
Ai
X (z) z
( z zi )
z zi
b X(z)有一个r阶重极点
X z A0
Z1
j

d
r
Ajz ( z z1 )
(r j) (r j)

j 1
z
1<|z|<2
k2
k 1
( k 1)
z ( z 2)
2
|z|<2
2 ( k 1)

第8章 z变换离散时间系统的z变换分析

第8章 z变换离散时间系统的z变换分析
1 z Z[u( n)] u( n)z z -1 1 z z 1 n 0 n 0
-n -n
收敛域 为 z >1
3. 斜变序列
间接求 解方法 已知 两边对(z -1)求导
两边乘(z -1)

同理,两边再求导,得

4. 指数序列
x(n) a n u(n)
运用留数定理来进行运算。又称为留数法,即
f (n) Res[F ( z )z n1 ]z pm
m
略!
二、幂级数展开法(长除法)
F ( z ) f (n)z n f (0) f (1)z 1 f ( 2)z -2
n 0


一般为变量z的有理分式,可用长除法,

s = 2,
例题 解
求x(n) = ?


见P60~61,表8-2、8-3、8-4(逆z变换表) 作业:P103,8-5 (1)(2)
8.5 z变换的基本性质
一、线性 若 x(n) ←→ X(z) y(n) ←→ Y(z)

Rx1 < |z| < Rx2 Ry1 < |z| < Ry2
ax(n) + by(n) ←→ aX(z) + bY(z)
F ( z ) f (0) f (1) z 1 f (2) z 2
所以
f (0) 0, f (1) 1, f (2) 0, f (3) 3, f (4) 4,
重点!
三、部分分式展开法
一般Z变换式是有理函数
以下研究因果序列的逆变换,即
X(z) (|z|>R) ← Z → x(n)
对于N阶LTI离散系统的差分方程:

第七章 z变换、离散时间系统的z域分析 PPT课件

第七章 z变换、离散时间系统的z域分析 PPT课件

1
n
u(n)的z变换,
2
3
并标明收敛域,绘出零极点图。
解:Zx(n)
x(n)zn
1
n
z

n
1
n
z
n
1
n

1
n
n-
n0 2
n0 3
n0 2z n0 3z
当 1 2z
1即 z
1时,
1
n
2 n0 2z
1 1-1/(2z)
z z1
2
当1 3z
1即 z
1时,
1
n
X (z) k A
m
z
m0 z z
m
其中,z 是 X (z)的极点,z 0。
m
z
0
A m
z
z m
X (z) z
zzm
k
X (z)
Az m
m0 z z
m
k
m0
A m
z m
n
u
(
n),
(右边Fra bibliotek序列

x(n)
Z
X 1
(z)
Z
1
k
m0
A m
z
z z
m
k
m0
A m
z m
n
u(n
1),(左边序列)
级数的系数就是序列x(n)。
• 右边序列,N(z)、D(z)按z的降幂(或z-1的升幂)排列
X (z) x(n)zn x(0)z0 x(1)z1 x(2)z2 n0
• 左边序列,N(z)、D(z)按z的升幂(或z-1的降幂)排列
1
X (z) x(n)zn x(1)z1 x(2)z2 x(3)z3 n

传递函数

传递函数

第2章 线性离散系统的Z变换分析法
—计算机控制技术—
2-4-2 用Z变换求解差分方程
平移定理
差分方程
代数方程
用Z变换求解差分方程的步骤:
(1)对差分方程作Z变换; (2)利用已知初始条件代入Z变换式,求出Y(z)表达式; (3)对Y(z)求Z反变换,求出差分方程的解:
y(kT )
1
[Y ( z )]
x(0)=0, x(1)=1,求其时间响应式。
解: 根据超前定理,其差分方程的Z
z X ( z) - z x(0) - zx(1) 3zX ( z) - 3zx(0) 2 X ( z) 0
2 2
整理后得
( z 2 3z ) x(0) z x(1) X ( z) 2 z 3z 2
Z传递函数的推导:
设n阶定常离散系统的差分方程为:
y(k ) a1 y(k 1) an y(k n) b0 u(k ) b1u(k 1) bm u(k m)
在零初始条件下,取Z变换
(1 a1 z 1 an z n )Y ( z ) (b0 b1 z 1 bm z m )U ( z )
i 0
k 1
通解 + 特解
第2章 线性离散系统的Z变换分析法
—计算机控制技术—
2.古典解法(解析法)
通解求法:
与式(2-1)对应的齐次方程为
y(kT ) a1 y(kT T ) a2 y(kT 2T ) an y(kT nT ) 0
(2-3)
k A 通解具有 的形式,代入式(2-3),有
第2章 线性离散系统的Z变换分析法
—计算机控制技术—

第八章 Z变换与Z域分析

第八章 Z变换与Z域分析

z (k ) z 1 z k 3 ( k 1) z 3
由线性性质得
|z|>1 |z|<3
z z 2z 4z F ( z) z 1 z 3 ( z 1)( z 3)
2
1<|z|<3
2、移位特性 (1)双边z变换 若f (k )是双边序列,其双边z变换为 f (k ) F ( z )
3<|z|<∞
根据时域乘ak性质,得
1 k F ( z ) Z [ f(k) Z f1 (k ) F1 (2 z ) ] 2 3 (2 z )2 4z2 2z 3 2z 3
2 k 0 1 k 2
z
k
z
2
z a 2 a 1 z 1 za
或者
a 2 z za
|z|>|a|
a 2 z F ( z ) Z [a k 2 ] Z [a a
例 8.2-3 已知f(k)=3k[ε(k+1)-ε(k-2)],求f(k)的双边Z变换

ZT [ e
n
z e j 0 z n j 0 n ZT [ e ] z e j 0 ZT [ cos0 n] ZT [ (e
n n j 0 n
j 0 n
]
z
e
j 0 n
) / 2]
z ( )/2 j 0 j 0 z e z e z ( z cos0 ) 2 z 2 z cos0 2 ( z )
Rx 2
6、双边序列 F ( z)
k
f (k ) z

k
f (k ) z
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
z
二 双边z 变换
例1:f ( k ) 3( 2)k ( k ) ( 5)k ( k 1)
分析:若存在双边z变换,必须存在公共收敛域
z 左边序列: ( 5) ( k 1) z5 2 z 5 (4 z 13) z F (z) 2 z 3 z 10
rzs (k ) 0 k 0
r(k-1) r(k-2) r(z)/z r(z)/z2
(1 a1 z 1 a0 z 2 )Rzs ( z ) (b1 z 1 b0 z 2 ) E( z )
R zs ( z ) b1 z 1 b0 z 2 b1 z b0 系统函数 H(z)= 2 1 2 E ( z ) 1 a1 z a0 z z a1 z a0
y(k+2)+5y(k+1)+6y(k)=e(k+1)+e(k) 2) y (0)=-1,y(1)=1
初始条件:1)yzi(0)=-1,yzi(1)=1
当 e(k)=ε(k),求在两种初始条件下系统全响应y (k) 方法二: y (0)=-1,y(1)=1 全响应的初始值 设全响应y(k)的单边z变换为Y(z),直接对等式两边做z变换 注意:要 考虑激励 函数的初 始值 y(k+1) y(k+2) e(k+1) z[Y(z)-y(0)] z2[Y(z)-y(0)-z-1y(1)] z[E(z)-e(0)]
z变换分析法
重点:
1:零状态
2:直接求全响应
例:双边序列f(k)的双边z变换为F(z),那f(-k)的双边z变 换
f(-2) f(-1)
f(k)
f(1)
f(0) f(1) f(2) f(2)
f(-1) f(-2) f(0)
f(k)
0
2
k
1
0
0
k
1 2
f (k ) f (2)z f (1)z f (0)z f (1)z f (2)z f (k ) f (2)z2 f (1)z1 f (0)z0 f (1)z1 f (2)z2 f (2)z f (1)z f (0)z f (1)z f (2)z
rzs (k ) (k 1) (k 2)
第四节 双边Z变换
重点:双边z变换及反变换
一 左边序列的z变换
f ( k ) b ( k 1)
k
0 ( k 1) 1
k0 k 1
F ( z ) b 1 z b 2 z 2 z
初始条件:1)yzi(0)=-1,yzi(1)=1 解: 1)yzi(0)=-1,yzi(1)=1
当 e(k)=ε(k),求在两种初始条件下系统全响应y (k) 全响应=零输入响应+零状态
2):先求零状态响应 1 1 1 k k y zs ( k ) [ ( 2) ( 3) ] ( k ) 3 2 6
2:当2<|z|<5
-5
-5 -2
Re[ z ]z ] Re[
3z 右边序列: 3 ( 2) k ( k ) z2 z ( 5) k ( k 1) 左边序列: z5 f ( k ) 3( 2)k ( k ) ( 5)k ( k 1)
3:当|z|<2
零状态响应
Rzs ( z ) E ( z ) H ( z )
1 求H(z) H ( z ) H ( S ) | S z 2 Z域相乘 Rzs(z)=H(z)E(z) 3 反Z变换 rzs(k)=Z-1{H(z)E(z)}
三 系统的全响应求解 例3:一线性因果系统 y(k+2)+5y(k+1)+6y(k)=e(k+1)+e(k) 2) y (0)=-1,y(1)=1
0 1 2 2 1
一 零输入响应 例1:一线性因果系统 y(k)+5y(k-1)+6y(k-2)=0
初始条件:yzi(-1)= 2/3 ,yzi(-2)= -7/18 ,求零输入响应。
二 零状态响应
rzs (k ) a1rzs (k 1) a0rzs (k 2) b1e(k 1) b0e(k 2)
1:当|z|>5
j Im[z]
j Im[z]
f (k ) [3(2)k (5)k ] (k )
2:当2<|z|<5
-5
-5 -2
Re[ z ]z ] Re[
3z 右边序列: 3 ( 2) k ( k ) z2 z ( 5) k ( k 1) 左边序列: z5 f ( k ) 3( 2)k ( k ) ( 5)k ( k 1)
Y (z) 8 / 3 3 / 2 1 / 6 z z 2 z 3 z 1
8 3 1 k k y( k ) [ ( 2) ( 3) ] ( k ) 3 2 6
1 已知零输入的初始条件:全响应=零输入响应+零状态 2 已知全响应的初始条件:直接对等式两边做z变换
例4:一线性因果系统 单位函数响应h(k)={0,1,-2} (k=0,1,2), 当 e(k)= ε(k),求在系统零状态响应响应y (k) 解: H ( z ) z 1 2z 2
Rzs ( z ) H ( z ).E( z )
z2 z . 2 z z 1
Rzs ( z ) 1 1 2 (1 ) z z z 1 1 z Rzs ( z ) 2 z z ( z 1)
例3:一线性因果系统
y(k+2)+5y(k+1)+6y(k)=e(k+1)+e(k)
初始条件:y (0)=-1,y(1)=1 当 e(k)= ε(k),求在两种初始条件下系统全响应y (k)
z ( z 4) 2z Y (z) ( z 2)( z 3) ( z 2)( z 3)( z 1)
收敛条件:
f(k)
zb
b 1 z 1 z b
j Im[z]
b-1
b-2 b-3
0
Rx2
Re[z ]
k
左边序列的收敛域是最里面极点的圆内区域
f (k ) F ( z )
1 f (k ) F ( ) z
左边序列的z变换:
1 f(k)反褶为右边g(k)
2 g(k)变换G(z)
3 G( z ) z 1 F ( z )
例4:一线性因果系统 单位函数响应h(k)={0,1,-2} (k=0,1,2), 当 e(k)= ε(k),求在系统零状态响应响应y (k) 解:
h(k ) (k 1) 2 (k 2) rzs (k ) e(k ) h(k ) (k 1) 2 (k 2)
m b z m n z
a n
1
b z m 1 n 1 z
m
1

b z 1 1 a z 1
b 0 a0 z
注意:
1 H ( z) Z h(k )
2 H ( z ) H ( S ) |S z N (z) D( z )
例2:已知一线性系统
对H(z)做反变换得到h(n)
k
3z 右边序列: 3 ( 2) ( k ) z2
k
z 2 z 5
j Im[z]
-5
-2
Re[z ]
双边z 变换的收敛域是一环状区域,其中 收敛域里的极点对应的是右边序列; 收敛域外的极点对应的是左边序列;
(4 z 13) z 3z z 例2:F ( z ) 2 = z 3 z 10 z 2 z 5
y(k+2)+5y(k+1)+6y(k)=e(k+1)+e(k) 且 e(k)=ε(k), 求系统函数H(z),单位函数响应h(k)和系统零状态响应yzs(k) 解:
(1 / 2) z ( 2 / 3) z z 1 1/ 6 H (z) 2 z2 z3 z 5z 6
1 1 2 h(k ) (k ) [ ( 2) k ( 3) k ] ( k ) 6 2 3
rzs ( k ) e( k ) h( k )
根据卷积性质:
Rzs ( z ) E ( z ) H ( z )
注意:1
H ( z ) Z h(k )
rzs (k ) a1rzs (k 1) a0rzs (k 2) b1e(k 1) b0e(k 2)
自然响应
Yzs ( z ) H ( z ) E ( z )
z( z 1) ( z 2)( z 3)( z 1)
1/ 3 1/ 2 1/ 6 ( )z z 2 z 3 z 1
1 1 1 k k yzs (k ) [ (2) (3) ] ( k ) 3 2 6
3:当|z|<2
f (k ) [3(2) (5) ] (k 1)
k k
(4 z 13) z 3z z 例2:F ( z ) 2 = z 3 z 10 z 2 z 5
1:当|z|>5
j Im[z]
j Im[z]
f (k ) [3(2)k (5)k ] (k )
y zs (0) 0, y zs (1) 1
求零输入响应:
y zi (0) 1, y zi (1) 0
y(k)=[C 1(-2 )k+ C 2(-3 )k]ε(k)
8 3 1 k k y( k ) [ ( 2) ( 3) ] ( k ) 3 2 6
例3:一线性因果系统
f (k ) [3(2) (5) ] (k 1)
相关文档
最新文档