第2章-离散时间系统课件

2.1引言问题：

通过考察信号在采样时刻的行为，如何把一个连续时间系统转换为一个离散时间系统？注意：

1.采样数据系统是一个时变系统，本章回避这个问题，仅研究与计算机时钟相同步的那些时刻的信号。

2.面向计算机的数学模型仅仅给出在采样点上的特性，而物理过程本身仍是一个连续时间系统。2.1引言2.2连续时间信号的采样采样意味着“连续时间信号由一个数值序列所代替，这个序列代表了某些时刻的信号值”。过程变量——模拟变换器有关的采样——数值序列——处理后新的数值序列——转换为连续时间信号——作用于过程。(采样过程重构过程)计算机接收受控过程在离散时间上的测量值，并在离散时间上发送新的控制信号。描述信号在逐个样点上的变化，而不关心样点之间的特性。(差分方程)图2.1计算机控制系统的原理框图2.2连续时间信号的采样采样意味着“连续时间信号由一个数值2.3连续时间状态空间系统的采样问题：

1.利用A/D和D/A变换器可以把一个连续时间系统和计算机连接起来，如何描述这样的系统？

2.计算机里的信号序列{u(tk)}和{y(tk)}，如何寻求这两个序列的关系？把寻找一个与连续时间系统相对应的离散时间系统称为连续时间系统的采样，所得到的模型亦称为频闪模型。

图2.2连接有A/D和D/A转换器的连续时间系统框图

2.3连续时间状态空间系统的采样问题：图2.2连接有2.3连续时间状态空间系统的采样连续时间系统由下列状态空间方程表示：

(2.1)

系统具有r个输入，p个输出，且阶数为n。系统的零阶保持采样

在计算机控制中，普遍的把D/A变换器设计成这样，即在指定下一个变换之前，它一直保持模拟信号恒定不变。通常，这样称之为零阶保持电路。2.3连续时间状态空间系统的采样连续时间系统由下列状态空间采样时刻上系统变量之间的关系当给定系统在采样时刻tk时的状态，则通过求解方程(2.1)便能得到某个未来时刻t的系统状态。于是，当tk

t

tk+1时，t时刻的状态为：

(2.2)这样，t时刻的状态变量向量为x(tk)和u(tk)的线性函数。如果A/D和D/A变换器是完全同步的，并且变换时间可以忽略不计，那么我们可以把输入u和输出y看作是在同一瞬间进行的采样。

(2.3)采样时刻上系统变量之间的关系当给定系统在采样时刻tk时的状态采样系统在采样时刻上的系统方程为：

(2.3)式中，注意：方程(2.3)并不包含任何的近似。由于控制信号在两个采样时刻之间保持恒定，故方程(2.3)给出了状态变量和输出量在采样时刻上的准确值。把模型(2.3)式称之为系统(2.1)式的零阶保持采样。式(2.3)也称之为系统(2.1)的零阶保持等价系统。采样系统在采样时刻上的系统方程为：均匀采样的离散时间系统对于周期为h的周期采样，有：tk=k

h，这时，式(2.3)表示的模型便可以简化成时不变系统：

(2.4)式中：

(2.5)由式(2.5)得：

均匀采样的离散时间系统矩阵

和

满足方程：注意：单位矩阵I的维数等于输入信号的个数。采样周期h，矩阵

(h)和

(h)可以由下面的方阵得到：

(2.6)矩阵和满足方程：如何计算

和

下面给出可以采用的五种方法：Matlab中的数值计算；矩阵指数级数展开；Laplace变换；凯莱-哈密顿(Cayley-Hamiton)定理；变换成对角型或者约当型。如果n2，矩阵

和

：

如何计算和下面给出可以采用的五种方法：例1：一阶连续时间系统：应用式(2.5)，得到：因此，采样系统变为：例1：一阶连续时间系统：例2双重积分器双重积分器的微分方程为：令y=x1，=x2，其状态空间表达式为：得到：双重计分器的离散时间模型为：(2.7)例2双重积分器(2.7)采样之逆问题：从离散时间描述中获取相应的连续时间系统是否可能？需要满足什么条件下才使可能的？采样之逆问题：考察一阶差分方程：从例子1中发现，相应的连续时间系统是从：得到：表明：当a>0时，才能得到一个具有实系数的连续时间系统。一般情况下，从式(2.6)可以得到：此处的ln(

)为矩阵对数函数。表明：连续时间系统可由对一个方阵取它的矩阵对数得到。当矩阵

在负实轴上没有特征值时。对数才唯一存在。考察一阶差分方程：具有时延的系统的采样在工业过程的数学模型中，时间延迟是很常见的。由于具有时间延时的连续时间系统是无限维系统，所以这种系统的理论十分复杂。设系统描述为：(2.8)第1种情况：

h图2.3

u(t)、u(t-

)和h之间的关系具有时延的系统的采样在工业过程的数学模型中，时间延迟是很常见先计算系统(2.8)的零阶保持采样。式(2.8)在一个采样周期上的积分为：(2.9)信号u(t)在整个采样间隔上是分段恒定的，故，延迟信号u(t-

)也是分段恒定的；2.延迟信号在各个采样时刻之间会有变化。要计算式(2.9)的积分项，方便的办法是：把积分区间分成两部分，使得u(t-

)在每一部分中都是恒定的，即：先计算系统(2.8)的零阶保持采样。(2.9)信号u(t)在连续时间系统(2.8)的采样系统为：

(2.10)式中：(2.11)

式(2.10)的状态空间模型为：注意：1.引入了r个状态变量u(kh-h)，它表示为控制信号的先前值。2.连续时间系统(2.8)是无限维的，但对应的采样系统却是有限维的。3.要想规定系统的状态，就必须在整个等于时间延迟的时间区间上存储输入信号。4.应用零阶保持重构法，输入信号总可以用有限数目的采样之来表示连续时间系统(2.8)的采样系统为：(2.10)式中：(2长时间的延迟第2种情况：

>h设：其中，d为常数，于是可以导出如下方程：式中，

0和

1由式(2.11)确定，但需要用

代替原来式中的

对应的状态空间描述为：

(2.12)

注意：如果

>0，那么，额外的引入了d

r个状态变量来描述时间延迟，这里r是输入的个数。状态空间描述的特征多项式为

drA(

)，其中A(

)为矩阵

的特征多项式。长时间的延迟第2种情况：>h其中，d为常数，于是可以例3简单的造纸机模型

模型的状态方程为：采样间隔h=1。d=3，

=0.6。根据(2.12)，可得：其中：

例3简单的造纸机模型模型的状态方程为：采样间隔h=1。具有内部时延的系统设系统由下列方程描述：

(2.13)

图2.4具有内部时间延迟的系统

具有内部时延的系统设系统由下列方程描述：(2.13)图2假设，u(t)在采样区间h上分段恒定。现在，试图找出x1(kh)和x2(kh)的递归方程。对系统(2.13)进行采样，当

=0，且采样区间h，得到分块系统：定理：内部时间延迟对系统(2.13)以采样区间为h，且0<

h进行的周期性采样，得到的采样数据表达式为：(2.14)

假设，u(t)在采样区间h上分段恒定。定理：(2.14)(2.14)

式中，(2.15)注意：具有时间延迟

的采样数据系统(2.14)是由对系统(2.13)采样而得到的，其中采样间隔h、

h、

之间无时间延迟。这样，相对于所需的采样时间间隔给出了

1、

2、

21、

1、

2。这意味着应用采样系统的标准软件就可得到式(2.14)。(2.14)式中，(2.15)注意：2.4

离散时间系统集中讨论差分方程的特性。时不变离散时间系统可以用差分方程来表示：(2.16)

为了简单起见，采样时间取作单位时间，即h=1。系统方程的解为了分析离散时间系统，需要解系统方程(2.16)。假设初始条件x(k0)，输入信号u(k0),

u(k0+1)，…均为已知。那么，状态如何演变呢？2.4离散时间系统集中讨论差分方程的特性。(2.16)为通过下列简单的迭代法可以求解式(2.16)：(2.17)

上述的解包括有两个部分：一个依赖于初始条件；另一个为输入信号的加权和。通过下列简单的迭代法可以求解式(2.16)：(2.17)上例4差分方程的解考察离散时间系统：已知x(0)=[11]T，很容易证得：如果

i

<1，i=1，2，则x(k)将收敛于原点。如果

的一个特征值的绝对值大于1，则会有一个或两个状态均发散。例4差分方程的解考察离散时间系统：已知x(0)=[12.5状态空间模型的坐标系变换考察离散时间系统(2.16)，讨论如何引进新的坐标系。(2.16)

假设T是一个非奇异矩阵，且定义新的状态向量z(k)=Tx(k)，则：且：状态空间表示取决于表示状态所选择的坐标系。

2.5状态空间模型的坐标系变换考察离散时间系统(2.16)定理2.2

特征方程的不变量当通过非奇异变换矩阵T来引入新的状态变量时，特征方程：保持不变。证明：寻找一个变化矩阵，就如同从线性方程组：中解出T的n2个元素。可以选择合适的坐标变换，以便获得形式简单的系统方程。定理2.2保持不变。证明：寻找一个变化矩阵，就如同从线性

对角型假设

的具有相异的特征值，那么存在矩阵T，使得：成立。式中，

i为矩阵

的特征值，矩阵T的计算以后讨论。得到一组解耦的一阶差分方程：上述方程组的解就很简单，每一种振型都具有如下形式的解：(2.18)对角型假设的具有相异的特征值，那么存在矩阵T，使得：成立

约当型如果

具有多重特征值，则通常

无法对角化。设

是一个n

n的矩阵，则引入表示式：式中，Lk为k

k的矩阵，于是存在矩阵T使得：(2.19)

式中，k1+k2+

+kr=n。

i矩阵

的特征值，它们未必是相异的。方程(2.19)称作为约当型。约当型如果具有多重特征值，则通常无法对角化。式中，Lk2.6输入－输出模型1.脉冲响应函数在单输入－单输出的离散时间系统中，在一个有限区间上的输入和输出信号，可用下列有限维的向量来表示：把向量Y和U关联起来的一般线型模型便可写成：其中是N

N矩阵，Yp是考虑初始条件的结果。2.6输入－输出模型1.脉冲响应函数把向量Y和U关联起如果Y与U有因果关系，那么，一定是下三角矩阵，它的元素在m>k时为零。于是普通线性系统的输入和输出关系便可以写成：函数称为系统的脉冲响应函数其中yp一项表示考虑了系统的初始条件的结果。说明：对于零初始条件，脉冲响应的取值表示在m时刻上的单位脉冲在k时刻产生的输出。对于多输入和多输出系统，脉冲响应就直接是一个简单的矩阵值函数。对于是不变系统，脉冲响应只是k－m的函数，即：如果Y与U有因果关系，那么，一定是下三角矩阵，它的元素其因此，可得：离散时间系统的脉冲响应函数为：性质（脉冲响应的不变量）离散时间系统的脉冲响应与状态空间模型的坐标变换无关。因此，可得：离散时间系统的脉冲响应函数为：性质（脉冲响应的不脉冲传递算子位移算子运算法把系统看成是输入信号映射到输出信号的一种算子。该算子运算中，所有的信号都是看成是双向无限序列，采用周期为单位采样。正向平移算子用符号q表示：后向平移算子用符号q－1表示：性质：位移算子具有单位范数。脉冲传递算子应用平移算子简化高阶差分方程的运算，考虑方程：其中，d＝na－nb，称为系统的极点盈数，

na

nb。应用平移算子可得：如果引入多项式：和上述差分方程可以表示为：应用平移算子简化高阶差分方程的运算，考虑方程：其中，d＝n同样，如果采用后向平移算子，上述差分方程可以写成：把多项式A的系数顺序反过来就获得下列多项式：该多项式称之为互反多项式，引入互反多项式后，系统可以写成：平移算子何时可以进行乘法，除法，加法和减法的运算呢？同样，如果采用后向平移算子，上述差分方程可以写成：把多项式A如果成立，那么，是否亦成立？需要条件吗？举例说明：考虑差分方程：采用平移算子表示为：可得解为：如果采用迭代结果，如果y(k0)=y0，得到方程的解为：由于q－1具有单位范数，等式右边可以表示为收敛级数：显然，两种结果不一致。采用迭代结果，如果y(k0)=y0，得到方程的解为：由于q注意：如果假设存在某个k0，使得k

k0时的序列全部为零，那么就有可能建立一个算子代数学，它允许除以q的任意多项式。也就是说，当差分方程的初始条件为零时，可以进行代数学的运算。（2）脉冲传递算子采用算子演算可以很方便的把输入－输出关系表示成正向平移算子或者后向平移算子的有理函数，这种有理函数称为脉冲传递算子。根据状态空间模型，可得：因此，注意：（2）脉冲传递算子因此，从而导出：于是系统的脉冲传递算子为：同样，系统的脉冲传递算子也可以用后向平移算则表示为：性质（脉冲传递算子的不变量）离散时间系统的脉冲传递算子H(q)与状态空间模型的坐标变换无关。从而导出：于是系统的脉冲传递算子为：同样，系统的脉冲传递算子第2章--离散时间系统课件第2章--离散时间系统课件第2章--离散时间系统课件3脉冲传递函数（1）z变换它是一种研究带初始条件和不带初始条件的线性差分方程的便利工具。将一个半无限的时间序列映射为一个复变量的函数。说明平移算子和z变换变量范围。算子运算考虑的是双无限时间序列，z变换考虑的是半无限时间序列。初始条件。z变换考虑初始条件，变量z为复变量。形式。平移算子演算与z变换计算形式上非常紧密，处理差分方程可以任选其一。q是一个作用与序列的算子，z是一个复变量。从纯数学的观点来说，可以明显的看出两者的不同，在举例2中说明。3脉冲传递函数（1）z变换定义z变换：考虑离散时间信号{f(kh)：k=0,1,…}，z变换定义为：其中，

z是复变量，f的z变换记作为Zf或者F。z反变换定义为：式中，积分周线包含F(z)的全部奇点。定义z变换：考虑离散时间信号{f(kh)：k=0,1,…z变换的性质线性性质：时移性质：卷积性质：z变换的性质线性性质：时移性质：卷积性质：初值定理：终值定理：如果(1－z－1)F(z)在单位圆上或者在单位圆外没有任何极点的话，那么初值定理：终值定理：（2）脉冲传递函数差分方程：如果对等式两边取z变换，可得：因此：且（2）脉冲传递函数差分方程：如果对等式两边取z变换，可得：离散系统的脉冲响应与脉冲传递函数是z变换对，即：Z{h(k)}=H(z))系统的脉冲传递函数离散系统的脉冲传递算子与脉冲传递函数形式相同，但是含义不同。离散系统的脉冲响应与脉冲传递函数是z变换对，即：系统的脉冲传离散系统脉冲传递函数的计算根据连续时间传递函数来直接确定脉冲传递函数。设系统的传递函数为G(s)，它前面接一个零阶保持器。脉冲传递函数由给定的信号的响应所惟一确定。图连续时间系统的采样步骤如下：1）确定传递函数G(s)的系统的阶跃响应。2）确定相应的阶跃响应的z变换。3）用阶跃函数的z变换除。离散系统脉冲传递函数的计算根据连续时间传递函数来直接确定脉冲通过使用上述方法，就可以推导下面的表达式：如果对于一个大数值的s来说，传递函数G(s)以至少等于

s

－1的递减速率趋于零，并有相异的均不在原点的极点，则得到：式中，si为G(s)的极点，Res表示为残数。通过使用上述方法，就可以推导下面的表达式：如果对于一个大数值举例1（脉冲传递函数的计算）：考察一个单位阶跃输入，这样序列{u(kh)}是1的