1、第一章 控制系统的状态空间表达式 状态空间表达式n 阶 DuCxyBA1:r:mynA:rB:nmC:rD:A 称为系统矩阵，描述系统内部状态之间的联系；为输入（或控制）矩阵，表示输入对每个状态变量的作用情况；C 输出矩阵，表示输出与每个状态变量间的组成关系，D 直接传递矩阵，表示输入对输出的直接传递关系。 状态空间描述的特点考虑了“输入状态输出 ”这一过程，它揭示了问题的本质，即输入引起了状态的变化，而状态决定了输出。状态方程和输出方程都是运动方程。状态变量个数等于系统包含的独立贮能元件的个数，n 阶系统有 n 个状态变量可以选择。状态变量的选择不唯一。从便于控制系统的构成来说，把状态变量选
2、为可测量或可观察的量更为合适。建立状态空间描述的步骤：a 选择状态变量；b 列写微分方程并化为状态变量的一阶微分方程组；c 将一阶微分方程组化为向量矩阵形式，即为状态空间描述。状态空间分析法是时域内的一种矩阵运算方法，特别适合于用计算机计算。 模拟结构图（积分器 加法器 比例器）已知状态空间描述，绘制模拟结构图的步骤：积分器的数目应等于状态变量数，将他们画在适当的位置，每个积分器的输出表示相应的某个状态变量，然后根据状态空间表达式画出相应的加法器和比例器，最后用箭头将这些元件连接起来。 状态空间表达式的建立 由系统框图建立状态空间表达式：a 将各个环节（放大、积分、惯性等）变成相应的模拟结构图
3、；b 每个积分器的输出选作 ，输入则为 ；c 由模拟图写出状态方程和输出方程。ixix 由系统的机理出发建立状态空间表达式：如电路系统。通常选电容上的电压和电感上的电流作为状态变量。利用 KVL 和 KCL 列微分方程，整理。由描述系统的输入输出动态方程式（微分方程）或传递函数，建立系统的状态空间表达式，即实现问题。实现是非唯一的。方法：微分方程 系统函数 模拟结构图 状态空间表达式。注意：a 如果系统函数分子幂次等于分母幂次，首先化成真分式形式，然后再继续其他工作。b 模拟结构图的等效。如前馈点等效移到综合反馈点之前。c 对多输入多输出微分方程的实现，也可以先画出模拟结构图。5状态矢量的线性
4、变换。也说明了状态空间表达的非唯一性。不改变系统的特征值。特征多项式的系数也是系统的不变量。特征矢量 的求解：也就是求 的非零解。ip0)(xAIi状态空间表达式变换为约旦标准型（为任意矩阵）：主要是要先求出变换矩阵。a 互异根时，各特征矢量按列排。b 有重根时，设 3 阶系统， ， 为单根，对特征矢量 ， 求123 1p3法与前面相同， 称作 的广义特征矢量，应满足 。2p1 121)(pAI系统的并联实现：特征根互异；有重根。方法：系统函数 部分分式展开 模拟结构图 状态空间表达式。6由状态空间表达式求传递函数阵 )(sW的矩阵函数 表示第 j 个输入对第 i 个输出的传递关系。DBAsI
5、CW1)() rmijij状态空间表达式不唯一，但系统的传递函数阵 是不变的。)(s子系统的并联、串联、反馈连接时，对应的状态空间表达及传递函数阵 。方法：画出系统结)(sW构图，理清关系，用分块矩阵表示。第二章 控制系统状态空间表达式的解一线性定常系统齐次状态方程（ ）的解：Ax 0)(xetAt二矩阵指数函数状态转移矩阵1 表示 到 的转移。5 个基本性质。Ate)()0(xt2 的计算：ta 定义；b 变换为约旦标准型 ,ATJ1)(或 11TeeJttt或c 用拉氏反变换 记忆常用的拉氏变换对 11sILeAt 22212 cos;sin;)(1;!;)(1; sttatesntats
6、tt aatd 应用凯莱-哈密顿定理三线性定常系统非齐次方程（ ）的解： 。可由BuAx dButxt )()0()(拉氏变换法证明（当然给出拉氏变换法的求解思路） 。求解步骤：先求 ，然后将 B 和 u(t)Ate代入公式即可。特殊激励下的解。第三章 线性控制系统的能控性和能观性一能控性及能观性定义（线性连续定常）二线性定常系统的能控性判别（具有一般系统矩阵的多输入系统）判别方法（一）：通过线性变换 BuAx BuTAz111若 A 的特征值互异，线性变换（ ）为对角线标准型， ，能控性充要条件：Tz没有全为 0 的行。 变换矩阵 T 的求法。BT2若 A 的特征值有相同的，线性变换（ ）为
7、约当标准型， ，能控性充要条件：xJ1对应于相同特征值的部分，每个约当块对应的 中最后一行元素没有全为 0 的。 B1中对应于互异特征根部分，各行元素没有全为 0 的。变换矩阵 T 的求法。1这种方法能确定具体哪个状态不能控。但线性变换比较复杂，关键是求 、 、 。T1B判别方法（二）：直接从，判别能控的充要条件是 能控性判别矩阵 的秩为 n。BuAx ),(12BABMn在单输入系统中， 是一个 的方阵；Mn而多输入系统， 是一个 的矩阵，可通过r )(Trakn三线性定常系统的能观性判别判别方法（一）：通过线性变换 CxyATzyA1若 A 的特征值互异，线性变换（ ）为对角线标准型， ，
8、能观性充要条件：z AT1中没有全为 0 的列。 变换矩阵 T 的求法。TC若 A 的特征值有相同的，线性变换（ ）为约当标准型， ，能控性充要条件：xJ1对应于相同特征值的部分，每个约当块对应的 中第一列元素没有全为的。 对应于互C异特征根部分，对应的 中各列元素没有全为的。变换矩阵 T 的求法。C这种方法能确定具体哪个状态不能观。但线性变换比较复杂，关键是求 、 、 。1TC判别方法（二）：直接从，C 判别能观性的充要条件是 能观性判别矩阵 的秩为 n。1nCAN在单输入系统中， 是一个 的方阵；Nn而多输入系统， 是一个 的矩阵，可通过m)(TMrak六能控性与能观性的对偶原理若 ， ，
9、 ，则 与 对偶。TA12TCB12TB12),(11CA),(22BA对偶系统的传递函数阵是互为转置的。且他们的特征方程式是相同的。 与 对偶，则 能控性等价于 能观性， 能观性等价于 能控性。121212七能控标准型和能观标准型对于状态反馈，化为能控标准型比较方便；对于观测器的设计及系统辨识，能观标准型比较方便。 能控标准型（如果已知系统的状态空间表达式）判别系统的能控性。计算特征多项式 ，即可写出 。求011| aaAInn A变换矩阵 ， 。求 ，计算 ，111ncApT 1,1,0Bbn 1cT1bTc，也可以验证是否有 。1c 1cAT 能观标准型判别系统的能观性。计算特征多项式
10、，即可写出 。求011| aaAInn A变换矩阵 ， 。求 ，计算 ，112,TATno 101ncA02TbT102，也可以验证是否有 。002c 2o 如果已知传递函数阵，可直接写出能控标准型和能观标准型的状态空间表达。 01211)( assasWnnn 能控标准型： 12100naaA 0b110nc能观标准型： 12100naA 1210nb10c八线性系统的结构分解1按能控性分解（状态不完全能控，即 ） ，通过非奇异变换 完成。rkM1 xRc，前 个列矢量是 M 中 个线性无关的列，其他列矢量保证nnc RR 1211n非奇异的条件下是任意的。2按能观性分解（状态不完全能观，即
11、 ） ，通过非奇异变换 完成。rankN1 xRo，前 个行矢量是 N 中 个线性无关的行，其他行矢量保证 非奇异的条件下是noRR12111n 1o任意的。3按能控性和能观性分解（系统是不完全能控和不完全能观的） ，采用逐步分解法，虽然烦琐，但直观。步骤：首先按能控性分解（ 能控状态， 不能控状态） 。对不能控子系统按能观性分解（cxcx不能控能观状态， 不能控不能观状态） 。将能控子系统按能观性分解（ 能控能观状态，ocxocx cox能控不能观状态） 。综合各步变换结果，写出最后的表达式。另一种方法：化为约当标准型，判断各状态的能控性能观测性，最后按 4 种类型分类排列。九传递函数阵的实
12、现问题1实现的定义：由 写出状态空间表达式，甚至画出模拟结构图，称为传递函数阵的实现问)(sW题。条件：传递函数阵中每个元的分子分母多项式都是实常数； 元是 s 的真有理分式。注意：如果不是有理分式，首先求出直接传递矩阵 。)(limWDs2能控标准型和能观标准型实现单入单出系统， 是有理分式，可直接根据分子分母多项式系数写出能控标准 1 型和能观)(sW标准 2 型实现。3最小实现（维数最小的实现）为 最小实现的充要条件是 是完全能控能观的。CxyBuA)(s ),(CBA步骤：对给定的 ，初选一种实现（能控标准型或能观标准型） ，假设选能控标准型，判断是W否完全能观测，若完全能观测则就是最
13、小实现；否则进行能观性分解，进一步找出能控能观部分，即为最小实现。注意：传递函数阵 的实现不是唯一的，最小实现也不是唯一的。)(s十传递函数 中零极点对消与能控性和能观性之间的关系)(对单输入系统、单输出系统或者单输入单输出系统，系统能控能观的充要条件是传递函数没有零极点对消。而对多输入多输出系统，传递函数阵没有零极点对消只是最小实现的充分条件，也就是说，即使存在零极点对消，系统仍有可能是能控能观的（p147 例 3-19） 。对单输入单输出系统，若传递函数出现了零极点对消，还不能判断到底是不能控还是不能观，还是既不能控又不能观。第四章 稳定性与李雅普诺夫方法一 稳定性的定义李雅普诺夫给出了对
14、任何系统都普遍适用的稳定性定义。1平衡状态为齐次状态方程。满足对所有 t，都有 成立的状态矢量 称为系统的平衡状),(txf 0),(txfe ex态。稳定性问题都是相对于某个平衡状态而言的。通常只讨论坐标原点处的稳定性。2稳定性的几个定义李雅普诺夫意义下稳定（相当于自控里的临界稳定） ； 渐近稳定（相当于自控里的稳定） ；大范围渐近稳定，大范围渐近稳定的必要条件是整个状态空间只有一个平衡状态； 不稳定。二 李雅普诺夫第一法（间接法）1线性定常系统的稳定判据状态稳定性：平衡状态 渐近稳定的充要条件是 A 的所有特征值具有负实部。0ex输出稳定性：充要条件是传递函数的极点位于 s 的左半平面。2
15、非线性系统的稳定性线性化处理。 ； ，若 A 的所有特征值具有负实部，则原非线性系统在平xAexf衡状态 渐近稳定。若 A 的所有特征值至少有一个具有正实部，则原非线性系统在平衡状态 不ex ex稳定。若若 A 的所有特征值至少有实部为零，则稳定性不能有特征值的符号来确定。三李雅普诺夫第二法（直接法） 借助于一个李雅普诺夫函数来直接对平衡状态的稳定性做出判断。1预备知识是由 n 维矢量 x 定义的标量函数，且在 处，恒有 ，对任何非零矢量 x，如)(xV 0x0)(xV果 ，则称之为正定；0如果 ，则称之为负定；如果 则称之为半正定或非负定；如果 则称之)(x)(xV 0)(xV为半负定或非正
16、定；如果 或 ，则称之为不定。0)(xV为二次型标量函数， 为实对称阵。要判别 的符号只要判别 的符号即可。PxVT)( P)(xVP的定号判据（希尔维特斯判据）：首先求出 的各阶顺序主子式 ，若所有的 ，则i0i（ ）正定；若 的 ， 的 则 （ ）负定；)(x偶 数i0i奇 数i0iP)(x2李雅普诺夫函数对于一个给定系统，如果能找到一个正定的标量函数 ，而 是负定的，则这个系统是)(V)(渐近稳定的，这个标量函数 叫做李雅普诺夫函数。 )(xV李雅普诺夫第二法的关键问题就是寻找李雅普诺夫函数 的问题。)(x稳定性判据设 ，平衡状态为 ，如果存在标量函数 是正定的，即 时，有 ，)(xf
17、0ex)(V0x0)(xV时，有 ，且满足 ，则称原点平衡状态是渐近稳定的；如果当 时，0x0V)(V ，则系统是大范围渐近稳定的。)(设 ，平衡状态为 ，如果存在标量函数 是正定的，即 时，有 ，)(xf 0ex)(xV0x0)(xV时，有 ，且满足 ，但除 外，即 ， 不恒等于，则称原点平0x0V)(V0x0)(衡状态是渐近稳定的；如果当 时， ，则系统是大范围渐近稳定的。x)(设 ，平衡状态为 ，如果存在标量函数 是正定的，即 时，有 ，)(xf 0ex)(xV0x0)(xV时，有 ，且满足 ，但任意的 ， 恒等于，则称原点平衡状态是李0x0V)(V0雅普诺夫意义下稳定的。设 ，平衡状态
18、为 ，如果存在标量函数 是正定的，即 时，有 ，)(xf 0ex)(x0x0)(x时，有 ，且满足 ， ，则称原点平衡状态是不稳定的。0x0)(需要注意：这些判据定理知识充分条件，也就是说，没有找到合适的李雅普诺夫函数来证明原点的稳定性，不能说明原点一定是不稳定的。如果 是可找到的，那么通常是非唯一的，)(xV但不影响结论。 最简单的形式是二次型标量函数，但不一定都是简单的二次型。 构造)(xV需要较多技巧。)(xV四李雅普诺夫方法在线性系统中的应用线性定常连续系统渐近稳定判据定理： ，若 A 是非奇异的，原点 是唯一的平衡点。原点大范围渐近稳定的充要条件 0ex是对任意对称实正定矩阵 ，李雅
19、普诺夫方程 ，存在唯一的对称正定解 。QTAPQP该定理等价于的特征值具有负实部。但高阶系统求解特征值复杂。步骤：选定正定矩阵 ，通常为 ，代入李雅普诺夫方程，确定出 ，判断是否正定，进I而做出系统渐近稳定的结论。五非线性系统的李雅普诺夫稳定性分析雅可比矩阵法步骤： ，写出 ，计算雅可比矩阵 ，对给定正定矩阵 （通常 ） ，)(xf)(xf xfJ)( PI为正定的。并且 为系统的一个李雅普诺夫函数。)(PJxQT)(PffxVT第五章 线性定常系统的综合一线性反馈控制系统的基本结构及其特性1状态反馈 将系统的每一个状态变量乘以相应的反馈系数，然后反馈到输入端与参考输入相加，作为受控系统的控制
20、输入。K 称为状态反馈增益阵， 。 设原受控系统 ，nr),(0CBA=0。状态反馈闭环系统的状态空间表达式 简称CxyBvKAx)( ),(KK与原受控系统 比较，状态反馈增益阵的引入，并不增加系统的维数，但可以),(0CBA通过的选择改变闭环系统的特征值，从而使获得所要求的性能。2输出反馈 由输出端 y 引入输出反馈增益阵 H（ ） ，然后反馈到输入端与参考输入相加，mr作为受控系统的控制输入。状态空间表达式为 简称CxyBvAx)(),(CBHA通过的选择也可以改变闭环系统的特征值，从而改变性能，但可供选择的自由度远比小（通常 ） 。nm从输出到状态变量导数 的反馈 从输出 y 引入反馈
21、增益阵 G（ ）到状态变量的导x mn数 ，所得状态空间表达式为 简称x CxyBuGA)( ),(CBAH通过的选择也可以改变闭环系统的特征值，从而改变性能。以上三种反馈的共同点是，不增加新的状态变量，系统开环与闭环同维，其次，反馈增益阵都是常数矩阵，反馈为线性反馈。闭环系统的能控性与能观性a 状态反馈不改变受控系统 的能控性，但不保证系统的能观性不变。),(0CBAb 输出反馈不改变受控系统 的能观性，但不保证系统的能控性不变。二极点配置问题 就是通过选择反馈增益矩阵，将闭环系统的极点恰好配置在根平面所期望的位置，以获得所希望的动态性能。 只讨论单输入单输出系统采用状态反馈 对系统 任意配
22、置极点的充要条件是 完全能控。),(0cbA0给定 ，给定期望的极点，设计状态反馈控制器的方法：),(0cbA能控规范型法，适合于 。首先判断是否完全能控，是，则存在状态观测器。 通过线3n性变换 化为能控标准型，得到 。加入状态反馈增益矩阵xTc1 ),(cb,得到闭环系统 状态空间表达式，求出对应的闭环特征多项,0nkK KAK式 。由给定的期望极点，求出期望的闭环特征多项式 。|)(|)(bAIf )()*if将 与 比较，即可得到 。 把对应与 的 ，通过 *f ,110nk K1cT。进一步画出模拟结构图。,110nk当阶次较低时， ，可直接由反映物理系统的 A,b 矩阵求状态反馈增
23、益矩阵3，不通过非奇异变换，使设计工作简单。首先判断是否完全能控，是，则存,110nK在状态观测器。加入状态反馈增益矩阵 ,得到闭环系统 状态,110nkK ),(cbKA空间表达式，求出对应的闭环特征多项式 。由给定的期望极点，求出期|)(|)(bKAIf望的闭环特征多项式 。将 与 比较，即可得到 。()*if*f ,110nk进一步画出模拟结构图。注意，如果给定的是传递函数，则先画出其要求的模拟结构图，写出状态空间描述，然后做其他工作。2采用输出反馈不能任意极点配置，正是输出线性反馈的基本弱点。采用从输出到 的反馈 对系统 任意配置极点的充要条件是 完全能观。x ),(0cbA0设计 从
24、输出到 的反馈阵的问题就是其对偶系统 设计状态反馈阵的问题。00方法：（）能观标准型法，适合于 。首先判断是否完全能观，是，则存在输出反馈3n。通过线性变换 化为能观标准型，得到 。 加入输出反馈增益矩阵xTo2 ),(cbA,得到闭环系统 状态空间表达式，求出对应的闭环特征多项ngG,110 ),(cbGA式 。由给定的期望极点，求出期望的闭环特征多项式 。|)(|)(cGAIf )()*if将 与 比较，即可得到 。 把对应与 的 ，通过 *f Tng,110 GTO2。进一步画出模拟结构图。,110ng当阶次较低时， ，可直接由反映物理系统的 A,c 矩阵求状态反馈增益矩阵3，不通过非奇
25、异变换，使设计工作简单。首先判断是否完全能观，是，则,110nG存在输出反馈。加入从输出到 的反馈增益矩阵 ,得到闭环系统x ,110ngG状态空间表达式，求出对应的闭环特征多项式 。由给定),(cbAG |)(|)(GcAIf的期望极点，求出期望的闭环特征多项式 。将 与 比较，即可得到)()*if*f。进一步画出模拟结构图。,110ng三系统镇定问题 所谓系统镇定，是对受控系统 通过反馈使其极点均具有负实部，保证系统为渐),(0CBA近稳定。镇定问题是极点配置问题的一种特殊情况，它只要求把闭环极点配置在根平面的左侧，而并不要求将闭环极点严格地配置在期望极点上。状态反馈能镇定的充要条件是其不
26、能控子系统为渐近稳定。输出反馈能镇定的充要条件是结构分解中能控能观子系统是输出反馈能镇定的，其余子系统是渐近稳定的。输出到 的反馈实现镇定的充要条件是不能观子系统为渐近稳定。x五状态观测器作用：闭环极点的任意配置、系统解藕以及最优控制系统都离不开状态反馈。但状态变量并不是都能直接检测，有些根本无法检测，这就提出状态观测或状态重构问题。龙伯格提出的状态观测器理论，解决的状态重构问题，使状态反馈成为一种可实现的控制律。定义：动态系统 以 的输入 u 和输出 y 作为输入量，产生一组输出量 逼近于 ，即0 x，则称 为 的一个状态观测器。构造原则： 必须是完全能观或不能观子系统0|limxt 0是渐
27、近稳定的； 的输出 应以足够快的速度渐近于 ； 在结构上尽可能简单（具有尽可能低x x的维数） ，以便于物理实现。等价性指标 动态系统 原系统 xcyBuA 0cxyBuA得到)(x )(0xeAt只要系统是稳定的，即的特征值具有负实部，就可做到 与 是稳态等价的。重构状态方程原因：系统的状态是不能直接量测的，因此很难判断是否有 逼近于 ； 不一定能保证x的特征值均具有负实部。 克服这个困难，用对输出量的差值 的测量代替对状态误差y的测量，当 ，有 。 同时，引入反x0|limxt 0|)(|lim|li|li ccxyttt馈阵，使系统的特征值具有负实部。状态重构方框图为 p213 5.16
28、(a) 要求熟练记忆，这种状态观测器称为渐近观测器。状态观测器方程为 记为xCy BuGyxCAyGBuA )()( ),(GBCA这里的 G 称为输出误差反馈矩阵。可以证明，如果 的特征值具有负实部，那么状态误差 将逐渐衰减到，即估计状态 逼近于实际的状态 。逼近的速度取决于 G 的选择，即x的特征值的配置。CA观测器的存在性对于完全能观测的线性定常系统，其观测器总是存在的。观测器存在的充要条件是 不能观子系统是渐近稳定的。0观测器的极点配置定理：线性定常系统 ，其观测器 可以任意配置极点，即具有任意),(0CBA),(GBCA逼近速度的充要条件是 完全能观测。 极点配置方法：（1）能观标准
29、型法，适合于 。首先判断是否完全能观，是，存在观3n测器可以任意极点配置。 通过线性变换 化为能观标准型，得到 。 加入xT ),(cbA输出误差反馈阵 ,得到闭环系统状态空间表达式 ，求TngG,110 yGuBxx出对应的闭环特征多项式 。 由给定的期望极点，求出期望的闭环特征多|)(|)(cGAIf项式 。 将 与 比较，即可得到 。把对应与()*if*f TngG,110的 ，通过 。得观测器方程，GT ,g，进一步画出模拟结构图。)( yGBuxAyBuxcAx （或 当阶次较低时， ，可由特征值不变原理求状态反馈增益矩阵 ，不通过3n ,110ngG非奇异变换，使设计工作简单。首先判断是否完全能观，是，则存在观测器可以任意极点配置。引入输出误差反馈矩阵 ,得到观测器系统 状态空间表达式,110ng ),(BcA
