离散时间系统状态方程的解

现代控制理论知识点汇总Revised at 2 pm on December 25, 2020.第一章 控制系统的状态空间表达式１． 状态空间表达式 n 阶DuCx y Bu Ax x+=+= 1:⨯r u 1:⨯m y n n A ⨯: r n B ⨯: n m C ⨯:r m D ⨯:A 称为系统矩阵，描述系统内部状态之间的联系；Ｂ为输入（或控制）矩阵，表示输入对每个状态变量的作用情况；C 输出矩阵，表示输出与每个状态变量间的组成关系，Ｄ直接传递矩阵，表示输入对输出的直接传递关系。

２． 状态空间描述的特点①考虑了“输入－状态－输出”这一过程，它揭示了问题的本质，即输入引起了状态的变化，而状态决定了输出。

②状态方程和输出方程都是运动方程。

③状态变量个数等于系统包含的独立贮能元件的个数，n 阶系统有n 个状态变量可以选择。

④状态变量的选择不唯一。

⑤从便于控制系统的构成来说，把状态变量选为可测量或可观察的量更为合适。

⑥建立状态空间描述的步骤：a 选择状态变量；b 列写微分方程并化为状态变量的一阶微分方程组；c 将一阶微分方程组化为向量矩阵形式，即为状态空间描述。

⑦状态空间分析法是时域内的一种矩阵运算方法，特别适合于用计算机计算。

３． 模拟结构图（积分器 加法器 比例器）已知状态空间描述，绘制模拟结构图的步骤：积分器的数目应等于状态变量数，将他们画在适当的位置，每个积分器的输出表示相应的某个状态变量，然后根据状态空间表达式画出相应的加法器和比例器，最后用箭头将这些元件连接起来。

４． 状态空间表达式的建立① 由系统框图建立状态空间表达式：a 将各个环节（放大、积分、惯性等）变成相应的模拟结构图；b 每个积分器的输出选作i x ，输入则为i x；c 由模拟图写出状态方程和输出方程。

② 由系统的机理出发建立状态空间表达式：如电路系统。

通常选电容上的电压和电感上的电流作为状态变量。

利用KVL 和KCL 列微分方程，整理。

第三章控制系统状态方程求解3－1 线性连续定常齐次方程求解所谓齐次方程解，也就是系统的自由解，是系统在没有控制输入的情况下，由系统的初始状态引起的自由运动，其状态方程为：………………………………………………………(3－1)上式中，X是n×1维的状态向量，A是n×n的常数矩阵。

我们知道，标量定常微分方程的解为： (3)2〕与〔3－2〕式类似，我们假设〔3－1〕的解X(t)为时间t的幂级数形式，即：………………………………(3－3) 其中为与X〔t〕同维的矢量。

将〔3－3〕两边对t求导，并代入〔3－1〕式，得：上式对任意时间t都应该成立，所以变量t的各阶幂的系数都应该相等，即：即:……………………………………………〔3－4〕将系统初始条件代入〔3－3〕，可得。

代入〔3－4〕式可得：…………………………………………………………………〔3－5〕代入〔3－3〕式可得〔3－1〕式的解为： (3)6)我们记： (3)7〕其中为一矩阵指数函数，它是一个n×n的方阵。

所以〔3－6〕变为：……………………………………………………………………〔3－8〕当〔3－1〕式给定的是时刻的状态值时，不难证明：………………………………………………………………〔3－9〕从〔3－9〕可看出，形式上是一个矩阵指数函数，且也是一个各元素随时间t变化的n×n矩阵。

但本质上，它的作用是将时刻的系统状态矢量转移到t时刻的状态矢量，也就是说它起到了系统状态转移的作用，所以我们称之为状态转移矩阵(The State Transition Matrix)，并记：……………………………………………………………〔3－10〕所以：【例3－1】，求解：根据〔3－7〕式，3－2 的性质及其求法性质1：【证】根据的定义式〔3－7〕，【证毕】性质2：①②③【证】：①：根据(3－7)式，即有：②：由性质1及其关系①，③:由②式两边同时左乘，注意本身是一个n×n的方阵，，所以：即：从上式可知，矩阵指数函数的逆矩阵始终存在，且等于。

§7-1 概述一、 离散时间信号与离散时间系统离散时间信号：只在某些离散的时间点上有值的信号。

离散时间系统：处理离散时间信号的系统。

混合时间系统：既处理离散时间信号，又处理连续时间信号的系统。

二、 连续信号与离散信号连续信号可以转换成离散信号，从而可以用离散时间系统（或数字信号处理系统）进行处理：三、 离散信号的表示方法：1、 时间函数：f(k)<——f(kT)，其中k 为序号，相当于时间。

例如：)1.0sin()(k k f =2、 （有序）数列：将离散信号的数值按顺序排列起来。

例如：f(k)={1,0.5,0.25,0.125,……,}时间函数可以表达任意长（可能是无限长）的离散信号，可以表达单边或双边信号，但是在很多情况下难于得到；数列的方法表示比较简单，直观，但是只能表示有始、有限长度的信号。

四、 典型的离散时间信号1、 单位样值函数：⎩⎨⎧==其它001)(k k δ下图表示了)(n k -δ的波形。

连续信号离散信号 数字信号 取样量化这个函数与连续时间信号中的冲激函数)(t δ相似，也有着与其相似的性质。

例如：)()0()()(k f k k f δδ=， )()()()(000k k k f k k k f -=-δδ。

2、 单位阶跃函数：⎩⎨⎧≥=其它001)(k k ε这个函数与连续时间信号中的阶跃函数)(t ε相似。

用它可以产生（或表示）单边信号（这里称为单边序列）。

3、 单边指数序列：)(k a k ε比较：单边连续指数信号：)()()(t e t e t a at εε=，其底一定大于零，不会出现负数。

4、 单边正弦序列：)()cos(0k k A εφω+(a) 0.9a = (d) 0.9a =-(b) 1a = (e) 1a =-(c) 1.1a = (f) 1.1a =-双边正弦序列：)cos(0φω+k A五、 离散信号的运算1、 加法：)()()(21k f k f k f +=<—相同的k 对应的数相加。

离散系统状态方程一、概述离散系统是指系统的输入和输出都是离散信号的系统。

离散系统的状态方程描述了系统的状态随时间变化的规律。

常见的离散系统状态方程的求解方法有前向欧拉法和后向欧拉法。

二、离散系统状态方程的基本形式离散系统状态方程通常采用以下一般形式表示：x(k+1) = Ax(k) + Bu(k)y(k) = Cx(k) + Du(k)其中，x(k)是系统的状态向量，A是状态转移矩阵，B是输入矩阵，u(k)是系统的输入向量，y(k)是系统的输出向量，C是输出矩阵，D是直达矩阵。

三、前向欧拉法1. 基本原理前向欧拉法是一种离散时间下的常微分方程近似数值解法。

对于给定的状态方程，前向欧拉法通过欧拉逼近求解系统状态在下一个离散时间点的值。

2. 求解步骤给定系统状态方程 x(k+1) = Ax(k) + Bu(k)，可以用前向欧拉法求解下一个时刻的状态值：x(k+1) = (I + hA)x(k) + hBu(k)其中，h为离散时间间隔。

3. 优缺点前向欧拉法求解简单且易于实现，但由于近似计算的原因，在一些特定系统下可能存在数值不稳定的情况。

四、后向欧拉法1. 基本原理后向欧拉法是另一种常微分方程的数值解法。

对于给定的状态方程，后向欧拉法通过在欧拉逼近的基础上，利用下一个时刻的状态值反推当前时刻的状态值。

2. 求解步骤给定系统状态方程 x(k+1) = Ax(k) + Bu(k)，可以用后向欧拉法求解当前时刻的状态值：x(k) = (I - hA)^(-1)x(k+1) - h(I - hA)^(-1)Bu(k+1)其中，h为离散时间间隔。

3. 优缺点后向欧拉法相比于前向欧拉法，更加稳定且精确，但求解过程较为复杂，通常需要使用数值计算方法求解线性方程组。

五、总结离散系统状态方程描述了离散系统的状态随时间变化的规律，对于工程控制和信号处理等领域有着重要的应用。

前向欧拉法和后向欧拉法都是常见的离散系统状态方程求解方法，各有优缺点，应根据实际情况选择合适的方法进行求解。