现代控制理论基础

2.3 线性离散系统的状态空间表达式及连续系统的离散化 2.3 线性离散系统的状态空间表达式及连续系统的离散化 2.3.1 离散系统的状态空间表达式 1.一般形式 式中：T是采样周期。方程中的矢量，各系数矩阵的名称和维数都与连续系统相同，为简单常省去T将方程写成如下形式

即： 2.结构图。上述方程可用结构图来表示

3.差分方程和脉冲传递函数与离散状态空间表达式之间的转换 在单变量离散系统中，数学模型分为差分方程和脉冲传递函数两类，它们与离散状态空间表达式之间的变换，和连续系统分析相类似。 连续 D.E T.F S.E 脉冲传函 状态空间表达式 离散 差分方程

例2-11 已知脉冲传递函数为 试求其状态空间表达式。 解：1) G(z) 差分方程 状态空间表达式 y(k+2) + 5y(k+1) + 6y(k) = u(k+2) + 2u(k+1) + u(k)

2) G(z) 部分分式法 状态空间表达式

4.状态空间表达式 G(z) G(z) = C(zI G ) 1H + D

2.3.2 定常连续系统的离散化 对连续系统，若常用数字计算机进行实时控制或求解，首先必须把连续系统转化成离散系统，这个过程称之为连续系统的离散化。 定常连续系统 离散化 离散系统

1、直接离散化： 离散化的实质就是用一个矩阵差分方程去代替一个矩阵微分方程。 取t0 = kt , t = (k+1)T 在kT    (k+1)T ，其输入向量u(t) = u(kT)，则状态方程的解为

对第二项积分作变量代换： 令t = (k+1)T   ； dt =  d  上限：  = (k+1)T，t = (k+1) T   = 0 下限：  = kT , t = (k+1)T   = T y(t) = Cx(t)+ Du(t) y(kT) = Cx(kT)+ Du(kT)

例2-12 求 的离散化方程。 解：先求eAt： (t) = eAt =L1[ (sI A ) 1 ]

x(k+1) = Gx(k) + Hu(k)

3) 按脉冲传函与标准型状态空间表达式的关系写出离散化的状态空间表达式 2、由脉冲传函实现离散化 步骤： 1) 首先求连续系统的传递函数 2) 按照离散系统的结构图求脉冲传函 3) 按脉冲传函与标准型状态空间表达式的关系写出离散化的状态空间表达式 U(s) G0(s) Y(s)

例2-13 已知连续系统的传递函数为 试求其离散化状态空间表达式 解：因为离散化后的系统结构图为： U(s)

2.3.3 定常连续系统的离散化的近似方法 近似方法出发点：用差商代替微商 x(k+1) = [I +TA]x(k) + TBu(k) G = I +TA H =TB

例2-13 求 的近似离散化方程。 解： H

当T = 0.1时

2.4 离散系统状态方程的解 定常离散系统状态方程的解：（两种方法） 2.4.1 迭代法 x(k+1) = Gx(k) + Hu(k) 2.4 离散系统状态方程的解 定常离散系统状态方程的解：（两种方法） 2.4.1 迭代法 x(k+1) = Gx(k) + Hu(k) 依次将采样时刻k=0，1，2，3,…代入上式即可。 k=0时,x(1)=Gx(0)+Hu(0) k=1时,x(2)=Gx(1)+Hu(1)=G2X(0)+GHu(0)+Hu(1) …

几点讨论： (1) 定常离散系统的状态解由两部分组成： 由初始状态引起的响应——反映系统的自由运动——零输入响应 由输入引起的响应——反映系统的强迫运动——零状态响应。 (2) 第k个采样时刻的状态，只与采样时刻0, 1, 2, …, k-1时的输入值有关系，而与第k个次采样时刻输入值无关，这是惯性系统的一个基本特征；

(3) G k 称为定常离散系统的状态转移矩阵，记为φ(k)= G k 初始条件： φ(0) = G 0 = I

2.4.2 z变换法 (k) = G k =Z1[ (zI G ) 1 z ]

1. z变换求出的解是一个完整解，其中解的结构可分为自由解和强迫解两部分，可分别求出，对分析运动过程有本质的帮助。解的形式是一个闭式，即解析式。 2. 迭代法求出的解是一个数值解。只能求出某一时刻的数值。但迭代公式本身就是状态方程，简单方便，而且不用求出状态转移矩阵Gk；如果已求出φ(k)=Gk，则可用解的迭代公式求出自由分量和强迫分量.

例2-15：求线性定常离散系统的解 已知 解： (1) 用迭代法求解

（2）用z变换法求解: (k) = G k =Z1[ (zI G ) 1 z ]

又知 u(k) = 1 X(z) = (zI G ) 1 [ z x(0)+HU(z) ]

x(k) = Z1[ X(z )] 令k = 0，1，2，3，… 代入上式，可得 以上两种方法计算结果完全一致，只是迭代法是一个数值解，而z变换法则得到了一个解析表达式。

2.4.3 离散系统的状态转移矩阵 1．直接法 根据离散系统递推迭代法中的定义 Φ(k) = G k 来计算。该方法简单，易于计算机来解，但不易得到Φ(k)的封闭式。 2．z变换法 根据z变换法求取离散系统状态方程解中的对应关系，状态转移矩阵Φ(k)为 Φ(k) = Z1[ (zI  G ) 1 z ]

3．化系统矩阵G为标准形法 （1）当离散系统矩阵G的特征值均为单根时 当离散系统矩阵G的特征根均为单根时，经过线性变换可将系统矩阵G化为对角线标准形，即 P 1 GP =  Φ(k) = G k = P k P 1

例2-16 齐次离散系统状态方程为 试求其状态转移矩阵Φ(k)。 解： 化系统矩阵G为对角线标准形的变换矩阵P为

Φ(k) = G k = P k P 1

(2) 当离散系统矩阵G的特征值有重根时 Φ(k) = G k = QJ kQ 1 4．化为G的有限项法 应用凯莱-哈密尔顿定理，系统矩阵G满足其自身的零化多项式。离散系统状态转移矩阵可化为G的有限项，即 式中αi(k) ( i = 0，1，… n −1)为待定系数，可仿照连续系统的方法来求取。

例2-17 线性定常离散系统的状态方程为 试求系统的状态转移矩阵Φ(k)。 解：离散系统特征方程为

结 束