现代控制理论 工程硕士 第四章李雅普诺夫稳定性理论
4李氏稳定性12
f
d
1 x1 x 2 x1 x2 x x
3 2
A 小范围稳定系统
令
1 0 x
xe1
2 0 x
0 xe2 1
0 xe3 1
0 0
平衡状态的进一步认识
对于线性定常系统 x’=Ax 平衡状态xe是满足下述方程的解。 Axe=0 当矩阵A为非奇异时,线性系统只有一个平衡状态 xe=0; 当A为奇异时,则存在无限多个平衡状态 对于非线性系统,通常可有一个或几个平衡状态,它们分别为 对应于式f(x,t)0的常值解。
当h 40， h 0.3时，变化K：
Tm s s K p K m K 1 K 0 0
3 2
无论怎样调整系统的参数，（如K、Tm）都不能使系统稳定
结构不稳定系统
s 3 2 0 .3 2 s sK 0 40 2 40
0 K 2 0.3 40
古典控制理论中的相关稳定性判据：
二、李雅普诺夫意义下的稳定 1.李氏意义下的稳定
如果对每个实数 0 ，都对应存在 另一个实数 ( , t0 ) 0 ，满足
x0 xe ( , t0 )
的任意初始态 x0 出发的运动轨迹
x(t ; x0 , t0 ) ，在t 都满足：
x(t; x0 , t0 ) xe , t t0
系统的传递函数：
根据上述稳定性的定义，可以用 (t ) 函数作 为扰动来讨论系统的稳定性。 设线性定常系统在初始条件为零时，输入一 个理想单位脉冲 (t ) ，即:系统在零平衡状态下， 受到一个扰动信号的作用，如果当t趋于无穷时， 系统的输出响应C(t)收敛到原来的零平衡状态， 即 lim C (t ) 0
李亚普诺夫稳定性分析
李亚普诺夫稳定性分析
可是在相当长的一段时间里,李雅普诺夫第二法并没有 引起研究动态系统稳定性的人们的重视,这是因为当时 讨论系统输入输出间关系的经典控制理论占有绝对地 位。 ➢ 随着状态空间分析法引入动态系统研究和现代控制 理论的诞生,李雅普诺夫第二法又重新引起控制领域 人们的注意,成为近40年来研究系统稳定性的最主要 方法,并得到了进一步研究和发展。 ➢ 本章节将详细介绍李雅普诺夫稳定性的定义,李雅普 诺夫第一法和第二法的理论及应用。
定理2 设定常系统的状态方程为 x f (x)
其中xe=0为其平衡状态。 ➢ 若存在一个有连续一阶偏导数的正定函数V(x),满足 下述条件: 1) 若 V ( x ) 为负定的; 2) 当||x||→,有V(x)→, 则该系统在原点处的平衡状态是大范围渐近稳 定的。
李亚普诺夫稳定性分析
对上述李雅普诺夫稳定性定理的使用有如下说明:
况，则 V ( x ) 为正半定或负半定。不属以上所有情况的V ( x ) 不定。
李亚普诺夫稳定性分析
2. 李雅普诺夫第二法的主要定理
下面分别介绍李雅普诺夫稳定性分析的如下3个定理: ➢ 渐近稳定性定理 ➢ 稳定性定理 ➢ 不稳定性定理
李亚普诺夫稳定性分析
2. 李雅普诺夫第二法的主要定理
(1) 定常系统大范围渐近稳定性定理1
✓ 但对于时变系统来说,则这两者的意义很可能不同。
对于李雅普诺夫渐近稳定性，还有如下说明： ➢ 稳定和渐近稳定，两者有很大的不同。 ✓ 对于稳定而言，只要求状态轨迹永远不会跑出球域 S(xe,)，至于在球域内如何变化不作任何规定。 ✓ 而对渐近稳定，不仅要求状态的运动轨迹不能跑出 球域，而且还要求最终收效或无限趋近平衡状态xe。
现代控制理论第版课后习题答案
现代控制理论第版课后习题答案Document number：WTWYT-WYWY-BTGTT-YTTYU-2018GT《现代控制理论参考答案》第一章答案1-1 试求图1-27系统的模拟结构图，并建立其状态空间表达式。
解：系统的模拟结构图如下： 系统的状态方程如下： 令y s =)(θ，则1x y =所以，系统的状态空间表达式及输出方程表达式为1-2有电路如图1-28所示。
以电压)(t u 为输入量，求以电感中的电流和电容上的电压作为状态变量的状态方程，和以电阻2R 上的电压作为输出量的输出方程。
解：由图，令32211,,x u x i x i c ===，输出量22x R y =有电路原理可知：•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为：1-4 两输入1u ，2u ，两输出1y ，2y 的系统，其模拟结构图如图1-30所示，试求其状态空间表达式和传递函数阵。
解：系统的状态空间表达式如下所示： 1-5系统的动态特性由下列微分方程描述列写其相应的状态空间表达式，并画出相应的模拟结构图。
解：令..3.21y x y x y x ===，，，则有相应的模拟结构图如下： 1-6 （2）已知系统传递函数2)3)(2()1(6)(+++=s s s s s W ,试求出系统的约旦标准型的实现，并画出相应的模拟结构图解：ss s s s s s s s W 31233310)3(4)3)(2()1(6)(22++++-++-=+++= 1-7 给定下列状态空间表达式[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321321321100210311032010x x x y u x x x x x x ‘（1） 画出其模拟结构图 （2） 求系统的传递函数 解：（2）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+-=-=31103201)()(s s s A sI s W 1-8 求下列矩阵的特征矢量（3）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=6712203010A 解：A 的特征方程 061166712230123=+++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+---=-λλλλλλλA I 解之得：3,2,1321-=-=-=λλλ当11-=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---3121113121116712203010p p p p p p 解得： 113121p p p -== 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P（或令111-=p ，得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P ） 当21-=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---32221232221226712203010p p p p p p 解得： 1232122221,2p p p p =-= 令212=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1423222122p p p P（或令112=p ，得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=21213222122p p p P ）当31-=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---33231333231336712203010p p p p p p 解得： 133313233,3p p p p =-= 令113=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3313323133p p p P1-9将下列状态空间表达式化成约旦标准型（并联分解）（2）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡32121321321110021357213311201214x x x y y u x x x x x x解：A 的特征方程 0)3)(1(311212142=--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=-λλλλλλA I 当31=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--3121113121113311201214p p p p p p 解之得 113121p p p == 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P当32=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--1113311201214312111312111p p p p p p 解之得 32222212,1p p p p =+= 令112=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0013222122p p p P当13=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--332313332313311201214p p p p p p 解之得3323132,0p p p == 令133=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1203323133p p p P约旦标准型1-10 已知两系统的传递函数分别为W 1(s)和W 2(s)试求两子系统串联联结和并联连接时，系统的传递函数阵，并讨论所得结果 解：（1）串联联结 （2）并联联结1-11 （第3版教材）已知如图1-22所示的系统，其中子系统1、2的传递函数阵分别为求系统的闭环传递函数解：1-11（第2版教材） 已知如图1-22所示的系统，其中子系统1、2的传递函数阵分别为求系统的闭环传递函数 解：1-12 已知差分方程为试将其用离散状态空间表达式表示，并使驱动函数u 的系数b(即控制列阵)为（1）⎥⎦⎤⎢⎣⎡=11b解法1： 解法2：求T,使得⎥⎦⎤⎢⎣⎡=-111B T 得⎥⎦⎤⎢⎣⎡=-10111T 所以 ⎥⎦⎤⎢⎣⎡-=1011T 所以，状态空间表达式为第二章习题答案2-4 用三种方法计算以下矩阵指数函数At e 。
现代控制理论-复习
离散时间状态空间模型：源自掌握基本原理，离散模型的导出，经典例子。
离散时间状态空间模型的分析。
能控、能观性 能控、能观性的定义、实际意义、判别条件、例子。 能控标准型是能控的； 一般的能控系统可以等价变换为能控标准型； 系统的离散化不能保持能控性； 输出能控性、和状态能控性的关系。 能控能观性的对偶原理 基于传递函数的能控、能观性条件：零极点对消 倒立摆的例子
06
分析：运动分析、能控性、能观性、稳定性。
07
设计：稳定化控制器、极点配置、观测器、基于观测器
08
的输出反馈控制器、线性二次型最优控制器。
09
要求：概念、方法、意义
状态空间模型 通过分析其内在变化规律列出相应动态方程； 通过输入输出数据建立传递函数模型，进而给出其状态空间实现； 掌握处理传递函数的状态实现方法，从特殊到一般的方法，掌握一些特殊状态空间实现的形式：能控标准型、能观标准型、对角型，它们的意义。 状态空间模型的状态变量图； 由状态空间模型确定传递函数； 状态空间模型的性质 等价模型的概念（可以简化结构），状态空间实现的不惟一性，等价模型具有相同传递函数、相同极点、相同能控、能观性
稳定性
李雅普诺夫稳定性、渐近稳定性的概念、意义；
具体例子的解释；
李雅普诺夫稳定性理论的实质：能量的变化
存在一个能量函数，沿系统轨线，能量函数衰减。
以上分析数学上的准备：函数的定号性！
重点：线性系统的稳定性分析
李雅普诺夫方程
线性系统李雅普诺夫稳定性定理的描述、举例应用。
李雅普诺夫稳定性定理的几何意义。
系统性能的分析（李雅普诺夫稳定性部分）。
线性二次型最优控制器的描述：
闭环性能指标：
最优闭环系统特性：稳定性。
现代控制理论总结
现代控制理论总结第一章：控制系统的状态空间表达式1、状态变量，状态空间与状态轨迹的概念：在描述系统运动的所有变量中，必定可以找到数目最少的一组变量，他们足以描述系统的全部运动，这组变量就称为系统的状态变量。
以状态变量X1,，X2，X3，……X n为坐标轴所构成的n维欧式空间（实数域上的向量空间）称为状态空间。
随着时间的推移，x（t）在状态空间中描绘出一条轨迹，称为状态轨迹。
2、状态空间表达式：状态方程和输出方程合起来构成对一个系统完整的动态描述，称为系统的状态空间表达式。
3、实现问题：由描述系统输入输出关系的运动方程或传递函数建立系统的状态空间表达式，这样的问题称为实现问题单入单出系统传函：W(s)=，实现存在的条件是系统必须满足m<=n，否则是物理不可实现系统最小实现是在所有的实现形式中，其维数最低的实现。
即无零，极点对消的传函的实现。
三种常用最小实现：能控标准型实现，能观标准型实现，并联型实现（约旦型）4、能控标准型实现，能观标准型实现，并联型实现（约旦型）传函无零点系统矩阵A的主对角线上方元素为1，最后一行元素是传函特征多项式系数的负值，其余元素为0，A为友矩阵。
控制矩阵b除最后一个元素是1，其他为0，矩阵A，b具有上述特点的状态空间表达式称为能控标准型。
将b与c矩阵元素互换，另输出矩阵c除第一个元素为1外其他为0，矩阵A，c具有上述特点的状态空间表达式称为能观标准型。
传函有零点见书p17页……..5、建立空间状态表达式的方法：①由结构图建立②有系统分析基里建立③由系统外部描述建立（传函）6、子系统在各种连接时的传函矩阵：设子系统1为子系统2为1）并联：另u1=u2=u，y=y1+y2的系统的状态空间表达式所以系统的传递函数矩阵为：2）串联：由u1=u，u2=y1，y=y2得系统的状态空间表达式为：W(S)=W2(S)W1(S)注意不能写反，应为矩阵乘法不满足交换律3）反馈：系统状态空间表达式：第二章：状态空间表达式的解：1、状态方程解的结构特征：线性系统的一个基本属性是满足叠加原理，把系统同时在初始状态和输入u作用下的状态运动x（t）分解为由初始状态和输入u分别单独作用所产生的运动和的叠加。
自动控制理论 第10章 李雅普诺夫稳定性分析
２）如果ｘｅ＝０为系统的平衡状态，则李氏函数应满足Ｖ（ｘｅ）＝ Ｖ（０）＝０。但当ｘ（ｔ）≠ ０
时， 不管其分量大于零或小于零，均能使Ｖ（ｘ）＞０。
基于上述的性质，人们常以状态矢量ｘ的二次型函数Ｖ（ｘ）作为李氏函数
的候选函数，即
式中，ｘ为实变数矢量。只要矩阵Ｐ是正定的，则上式所示的Ｖ（ｘ）就符 合对李氏函数性质的要求。
对于连续定常系统，李雅普诺夫第二方法是根据Ｖ（ｘ）和
的性
质去判别它的稳定性。因此需要研究以下两个问题：
１）具备什么条件的函数才是李雅普诺夫函数，简称李氏函数。
２）怎样利用李氏函数去判别系统平衡状态的稳定性？
由对图10-2所示系统的讨论，可知李氏函数必须要同时具有如下两个性质：
１）李氏函数是自变量为系统的状态矢量ｘ（ｔ）的标量函数。
态是不稳定的。
2021/6/18
第十章 李雅普诺夫稳定性分析
6
为了能更直观地理解上述平衡状态稳定性的概念，
下图在二维状态平面上分别画出了系统平衡状态的稳 定、渐近稳定和不稳定３种情况。
2021/6/18
第十章 李雅普诺夫稳定性分析
7
自动控制理论
第二节 李雅普诺夫第二方法
正定函数
2021/6/18
11
自动控制理论
由上式可见，除了ｘｅ＝０外，系统的能量Ｖ（ｘ）在运动过程中由于 受到了阻尼器的阻尼作用而不断地减小，最后使Ｖ（ｘ）＝０。这个例子很 容易把能量函数Ｖ（ｘ）与实际系统联系起来。然而，对一般的系统而言， 至今还没有一个普遍适用“能量函数” 的表达式。对此，李雅普诺夫提出了 一个虚拟的能量函数，人们称它为李雅普诺夫函数，用Ｖ（ｘ）表示。
则称系统的平衡状态ｘｅ是渐近稳定的。
基于MATLAB的李雅普诺夫第二法稳定性分析
基于MATLAB的李雅普诺夫第二法稳定性分析引言：对于一个给定的控制系统，稳定性是系统的一个重要特性。
稳定性是系统正常工作的前提，是系统的一个动态属性。
在控制理论工程中，无论是调节器理论、观测器理论还是滤波预测、自适应理，都不可避免地要遇到系统稳定性问题，而且稳定性分析的复杂程度也在急剧增长。
当已知一个系统的传递函数或状态空间表达式时, 可以对其系统的稳定性进行分析；当系统的阶次较高时，分析、计算的工作量很大, 给系统的分析带来很大困难。
运用MATLAB 软件，其强大的科学计算能力和可视化编程功能, 为控制系统稳定性分析提供了强有力的工具。
一．MATLAB 语言简介MATLAB 是MATrix LABoratory 的缩写, 它是MA TLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。
它具有强大的矩阵计算能力和良好的图形可视化功能, 为用户提供了非常直观和简洁的程序开发环境, 因此被称为第四代计算机语言。
MA TLAB 发展至今, 现已集成了许多工具箱, 一般来说, 它们都是由特定领域的专家开发的, 用户可以直接是用工具箱学习、应用和评估不同的方法而不需要自己编写代码，大大提高了分析运算的效率，为此MA TLAB 语言在控制工程领域已获得了广泛地应用。
二．控制系统稳定性的基本概念稳定性是控制系统的重要特性, 也是系统能够正常运行的首要条件。
如何分析系统的稳定性并提出保证系统稳定的措施, 是自动控制理论的基本任务之一。
1892年，俄国数学家李雅普诺夫（Lyaponov）提出了分析稳定性的两种方法。
第一种方法，通过对线性化系统特征方程的根的分析情况来判断稳定性，称为间接法。
此时，非线性系统必须先线性近似，而且只能使用于平衡状态附近。
第二种方法，从能量的观点对系统的稳定性进行研究，称为直接法，对线性、非线性系统都适用。
线性系统理论精简版 ——5.系统的稳定性
内部稳定性和外部稳定性在满足一定条件下是等 价的（后面讨论）。
经典理论判稳方法及局限性 间接判定：方程求解－(对非线性和时变通常很难)
直接判定：单入单出中，基于特征方程的根是否都
分布在复平面虚轴的左半部分；以及采用劳斯判据、 奈魁斯特频率判据等。局限性是仅适用于线性定常， 不适用于非线性和时变系统。
0 1
xe , 3
0 1
5.2.2 李雅普诺夫稳定 定义：若状态方程
f ( x, t ) x 所描述的系统，对于任意的>0和任意初始
x2
时刻t0，都对应存在一个实数(,t0)>0，
使得对于任意位于平衡态xe的球域S(xe,) 的初始状态x0，当从此初始状态x0出发的 状态方程的解x都位于球域S(xe,)内，则 称系统的平衡态xe是李雅普诺夫意义下稳
现代控制理论判稳方法： 李雅普诺夫稳定性理论是稳定性判定的通用方法，适用于 各种系统。 李亚普诺夫第一法：先求解系统微分方程，根据解的性质
判定稳定性－－间接法。
李亚普诺夫第二法：直接判定稳定性。思路：构造一个李
亚普诺夫函数V(x)，根据V(x)的性质判稳。－－对任何复
杂系统都适用。
5.2
V ( X ) 0 X 0 V ( X ) 0 X 0
例5-2
2 2 V ( X ) x1 2x2
当 x1 0, x2 0 时，V ( X ) 0； 当 x1 0, x2 0 时， V ( X ) 0。所以，V(X)是正定的。
(2) 正半定性(准正定) 如果对任意非零向量 X ( X 0) ，恒有 V ( X )≥0， 且当 X 0时V ( X ) 0 ，则称 V ( X ) 为正半定的。即
现代控制理论第版课后习题答案
现代控制理论参考答案第一章答案1-1 试求图1-27系统的模拟结构图,并建立其状态空间表达式; 解：系统的模拟结构图如下： 系统的状态方程如下： 令y s =)(θ,则1x y =所以,系统的状态空间表达式及输出方程表达式为1-2有电路如图1-28所示;以电压)(t u 为输入量,求以电感中的电流和电容上的电压作为状态变量的状态方程,和以电阻2R 上的电压作为输出量的输出方程; 解：由图,令32211,,x u x i x i c ===,输出量22x R y =有电路原理可知：•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为：1-4 两输入1u ,2u ,两输出1y ,2y 的系统,其模拟结构图如图1-30所示,试求其状态空间表达式和传递函数阵;解：系统的状态空间表达式如下所示： 1-5系统的动态特性由下列微分方程描述列写其相应的状态空间表达式,并画出相应的模拟结构图; 解：令..3.21y x y x y x ===，，,则有 相应的模拟结构图如下： 1-6 2已知系统传递函数2)3)(2()1(6)(+++=s s s s s W ,试求出系统的约旦标准型的实现,并画出相应的模拟结构图解：ss s s s s s s s W 31233310)3(4)3)(2()1(6)(22++++-++-=+++=1-7 给定下列状态空间表达式[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321321321100210311032010x x x y u x x x x x x ‘（1） 画出其模拟结构图（2） 求系统的传递函数 解：2⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+-=-=31103201)()(s s s A sI s W 1-8 求下列矩阵的特征矢量3⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=6712203010A 解：A 的特征方程 061166712230123=+++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+---=-λλλλλλλA I 解之得：3,2,1321-=-=-=λλλ当11-=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---3121113121116712203010p p p p p p 解得： 113121p p p -== 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P 或令111-=p ,得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P 当21-=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---32221232221226712203010p p p p p p解得： 1232122221,2p p p p =-= 令212=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1423222122p p p P或令112=p ,得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=21213222122p p p P 当31-=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---33231333231336712203010p p p p p p 解得： 133313233,3p p p p =-= 令113=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3313323133p p p P 1-9将下列状态空间表达式化成约旦标准型并联分解2⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡32121321321110021357213311201214x x x y y u x x x x x x解：A 的特征方程 0)3)(1(311212142=--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=-λλλλλλA I 当31=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--3121113121113311201214p p p p p p 解之得 113121p p p == 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P 当32=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--1113311201214312111312111p p p p p p解之得 32222212,1p p p p =+= 令112=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0013222122p p p P 当13=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--332313332313311201214p p p p p p 解之得3323132,0p p p == 令133=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1203323133p p p P约旦标准型1-10 已知两系统的传递函数分别为W 1s 和W 2s试求两子系统串联联结和并联连接时,系统的传递函数阵,并讨论所得结果 解：1串联联结 2并联联结1-11 第3版教材已知如图1-22所示的系统,其中子系统1、2的传递函数阵分别为 求系统的闭环传递函数 解：1-11第2版教材 已知如图1-22所示的系统,其中子系统1、2的传递函数阵分别为 求系统的闭环传递函数 解：1-12 已知差分方程为试将其用离散状态空间表达式表示,并使驱动函数u 的系数b 即控制列阵为 1⎥⎦⎤⎢⎣⎡=11b 解法1： 解法2：求T,使得⎥⎦⎤⎢⎣⎡=-111B T 得⎥⎦⎤⎢⎣⎡=-10111T 所以 ⎥⎦⎤⎢⎣⎡-=1011T 所以,状态空间表达式为第二章习题答案2-4 用三种方法计算以下矩阵指数函数At e ;2 A=1141⎛⎫⎪⎝⎭解：第一种方法： 令0I A λ-=则11041λλ--=-- ,即()2140λ--=;求解得到13λ=,21λ=- 当13λ=时,特征矢量11121p p p ⎡⎤=⎢⎥⎣⎦由 111Ap p λ=,得11112121311341p p p p ⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦即112111112121343p p p p p p +=⎧⎨+=⎩,可令112p ⎡⎤=⎢⎥⎣⎦当21λ=-时,特征矢量12222p p p ⎡⎤=⎢⎥⎣⎦由222Ap p λ=,得121222221141p p p p -⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦即1222121222224p p p p p p +=-⎧⎨+=-⎩ ,可令212p ⎡⎤=⎢⎥-⎣⎦则1122T ⎡⎤=⎢⎥-⎣⎦,111241124T -⎡⎤⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦第二种方法,即拉氏反变换法：第三种方法,即凯莱—哈密顿定理 由第一种方法可知13λ=,21λ=-2-5 下列矩阵是否满足状态转移矩阵的条件,如果满足,试求与之对应的A 阵;3()22222222t tt t ttt t e e e e t e e e e --------⎡⎤--Φ=⎢⎥--⎣⎦ 4()()()()3333112412t t t t t tt t e e e e t e e e e ----⎡⎤+-+⎢⎥Φ=⎢⎥⎢⎥-++⎢⎥⎣⎦解：3因为 ()10001I ⎡⎤Φ==⎢⎥⎣⎦,所以该矩阵满足状态转移矩阵的条件 4因为()10001I ⎡⎤Φ==⎢⎥⎣⎦,所以该矩阵满足状态转移矩阵的条件2-6 求下列状态空间表达式的解：初始状态()101x ⎡⎤=⎢⎥⎣⎦,输入()u t 时单位阶跃函数;解： 0100A ⎡⎤=⎢⎥⎣⎦ 因为 01B ⎡⎤=⎢⎥⎣⎦,()()u t I t =2-9 有系统如图所示,试求离散化的状态空间表达式;设采样周期分别为T=和1s,而1u 和2u 为分段常数; 图 系统结构图 解：将此图化成模拟结构图 列出状态方程则离散时间状态空间表达式为 由()At G T e =和()0TAt H T e dtB =⎰得：当T=1时 ()()()()11111001111k e e x k x k u k e ke ----⎡⎤-⎡⎤+=+⎢⎥⎢⎥--⎢⎥⎣⎦⎣⎦当T=时 ()()()()()0.10.10.10.11001110.90.1k e e x k x k u k e k e ----⎡⎤-⎡⎤⎢⎥+=+⎢⎥⎢⎥---⎣⎦⎣⎦第三章习题3-1判断下列系统的状态能控性和能观测性;系统中a,b,c,d 的取值对能控性和能观性是否有关,若有关,其取值条件如何 1系统如图所示： 解：由图可得： 状态空间表达式为：由于•2x 、•3x 、•4x 与u 无关,因而状态不能完全能控,为不能控系统;由于y 只与3x 有关,因而系统为不完全能观的,为不能观系统; 3系统如下式：解：如状态方程与输出方程所示,A 为约旦标准形;要使系统能控,控制矩阵b 中相对于约旦块的最后一行元素不能为0,故有0,0≠≠b a ;要使系统能观,则C 中对应于约旦块的第一列元素不全为0,故有0,0≠≠d c ; 3-2时不变系统试用两种方法判别其能控性和能观性; 解：方法一：方法二：将系统化为约旦标准形;⎥⎦⎤⎢⎣⎡=1-111T ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=21212121T 1- B T -1中有全为零的行,系统不可控;CT 中没有全为0的列,系统可观; 3-3确定使下列系统为状态完全能控和状态完全能观的待定常数i i βα和 解：构造能控阵：要使系统完全能控,则211αα≠+,即0121≠+-αα 构造能观阵：要使系统完全能观,则121αα-≠-,即0121≠+-αα 3-4设系统的传递函数是1当a 取何值时,系统将是不完全能控或不完全能观的 2当a 取上述值时,求使系统的完全能控的状态空间表达式; 3当a 取上述值时,求使系统的完全能观的状态空间表达式; 解：1 方法1 ：)6)(3)(1()()()(++++==s s s as s u s y s W 系统能控且能观的条件为Ws 没有零极点对消;因此当a=1,或a=3或a=6时,系统为不能控或不能观; 方法2：系统能控且能观的条件为矩阵C 不存在全为0的列;因此当a=1,或a=3或a=6时,系统为不能控或不能观;2当a=1, a=3或a=6时,系统可化为能控标准I 型3根据对偶原理,当a=1, a=2或a=4时,系统的能观标准II 型为 3-6已知系统的微分方程为：u y y y y 66116...=+++试写出其对偶系统的状态空间表达式及其传递函数; 解：63611603210=====b a a a a ，，，，系统的状态空间表达式为 传递函数为其对偶系统的状态空间表达式为： 传递函数为61166)(23+--=s s s s W 3-9已知系统的传递函数为 试求其能控标准型和能观标准型;解：345213486)(222++++=++++=s s s s s s s s W系统的能控标准I 型为 能观标准II 型为3-10给定下列状态空间方程,试判别其是否变换为能控和能观标准型;解：[]100210311032010=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=C b A ，， 3-11试将下列系统按能控性进行分解1[]111,100,340010121-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=C b A 解：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--==9310004102b A Ab bM rankM=2<3,系统不是完全能控的; 构造奇异变换阵c R ：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==010*********R Ab R b R ，，,其中3R 是任意的,只要满足c R 满秩;即⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=031100010c R 得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=-010*******c R 3-12 试将下列系统按能观性进行结构分解1 []111,100,340010121-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=C b A 解： 由已知得[]111,100,340010121-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=C b A 则有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=4742321112CA CA C Nrank N=2<3,该系统不能观构造非奇异变换矩阵10R -,有10111232001R --⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦ 则0311210001R --⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦3-13 试将下列系统按能控性和能观性进行结构分解1[]211,221,102322001=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=C b A 解：由已知得211121226202M A Ab Ab ⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥-⎣⎦rank M=3,则系统能控 rank N=3,则系统能观所以此系统为能控并且能观系统取211121226202c T ⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦,则1217344173215344c T -⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦则002105014A ⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦,12100c B T b -⎡⎤⎢⎥==⎢⎥⎢⎥⎣⎦,[]271323c c cT == 3-14求下列传递函数阵的最小实现; 1 ()111111w s s ⎡⎤=⎢⎥+⎣⎦解： 01α=,01111B ⎡⎤=⎢⎥⎣⎦,1001c A -⎡⎤=⎢⎥-⎣⎦ 1001c B ⎡⎤=⎢⎥⎣⎦,1111c C ⎡⎤=⎢⎥⎣⎦,0000c D ⎡⎤=⎢⎥⎣⎦ 系统能控不能观取101101R -⎡⎤=⎢⎥⎣⎦,则01101R -⎡⎤=⎢⎥⎣⎦所以10010ˆ01A R AR --⎡⎤==⎢⎥-⎣⎦,1011ˆ01c B R B -⎡⎤==⎢⎥⎣⎦ 010ˆ10c C C R ⎡⎤==⎢⎥⎣⎦,00ˆ00D ⎡⎤=⎢⎥⎣⎦所以最小实现为ˆ1m A =,[]ˆ11m B =,1ˆ1m C ⎡⎤=⎢⎥⎣⎦,00ˆ00m D ⎡⎤=⎢⎥⎣⎦ 验证：()()1111ˆˆˆ111m mm C sI A B w s s -⎡⎤-==⎢⎥+⎣⎦3-15设1∑和2∑是两个能控且能观的系统1试分析由1∑和2∑所组成的串联系统的能控性和能观性,并写出其传递函数； 2试分析由1∑和2∑所组成的并联系统的能控性和能观性,并写出其传递函数; 解： 11∑和2∑串联当1∑的输出1y 是2∑的输入2u 时,331222x x x x =-++010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=--+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]001y x = 则rank M=2<3,所以系统不完全能控; 当2∑得输出2y 是1∑的输入1u 时011034100021x x u ⎡⎤⎡⎤⎢⎥⎢⎥=--+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]210y x = 因为 2001016124M bAbA b ⎡⎤⎢⎥⎡⎤==-⎣⎦⎢⎥⎢⎥--⎣⎦rank M=3 则系统能控因为2210321654c N cA cA ⎡⎤⎡⎤⎢⎥⎢⎥==--⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦rank N=2<3 则系统不能观 21∑和2∑并联010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=--+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]211y x = 因为rank M=3,所以系统完全能控 因为rank N=3,所以系统完全能观现代控制理论第四章习题答案4-1判断下列二次型函数的符号性质：1222123122313()31122Q x x x x x x x x x x =---+-- 2222123122313()4262v x x x x x x x x x x =++--- 解：1由已知得110∆=-<,2112013-∆==>-,31111711302411112--∆=--=-<--- 因此()Q x 是负定的 2由已知得110∆=>,2113014-∆==>-,3111143160131--∆=--=-<--因此()Q x 不是正定的 4-2已知二阶系统的状态方程：试确定系统在平衡状态处大范围渐进稳定的条件;解：方法1：要使系统在平衡状态处大范围渐进稳定,则要求满足A 的特征值均具有负实部;即：有解,且解具有负实部; 即：1122112212210a a a a a a +<>且方法2：系统的原点平衡状态0e x =为大范围渐近稳定,等价于T A P PA Q +=-;取Q I =,令11121222PP P P P ⎡⎤=⎢⎥⎣⎦,则带入TA P PA Q +=-,得到 若 11211211222111221122122112222204()()0022a a a a a a a a a a a a a a +=+-≠,则此方程组有唯一解;即 其中11221221det A A a a a a ==- 要求P 正定,则要求 因此11220a a +<,且det 0A >4-3试用lyapunov 第二法确定下列系统原点的稳定性;11123x x -⎡⎤=⎢⎥-⎣⎦ 21111x x -⎡⎤=⎢⎥--⎣⎦解：1系统唯一的平衡状态是0e x =;选取Lyapunov 函数为2212()0V x x x =+>,则 ()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;2系统唯一的平衡状态是0e x =;选取Lyapunov 函数为2212()0V x x x =+>,则 ()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;4-6设非线性系统状态方程为： 试确定平衡状态的稳定性;解：若采用克拉索夫斯基法,则依题意有： 取P I =很明显,()Q x 的符号无法确定,故改用李雅普诺夫第二法;选取Lyapunov 函数为2212()0V x x x =+>,则()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;4-9设非线性方程：试用克拉索夫斯基法确定系统原点的稳定性; 解：1采用克拉索夫斯基法,依题意有：x →∞,有()V x →∞; 取P I =则2121013()132x Q x x ⎡⎤-+=⎢⎥-+⎣⎦,根据希尔维斯特判据,有： 2221121210310310132x x x -∆=∆==->-+，（）,()Q x 的符号无法判断; 2李雅普诺夫方法：选取Lyapunov 函数为421233()042V x x x =+>,则 ()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;4-12试用变量梯度法构造下列系统的李雅普诺夫函数 解：假设()V x 的梯度为： 计算()V x 的导数为：选择参数,试选112212211,0a a a a ====,于是得：12x V x ⎛⎫∇= ⎪⎝⎭,显然满足旋度方程12122121,0V V x xx x x x ∂∇∂∇∂∂===∂∂∂∂即,表明上述选择的参数是允许的;则有：如果121211202x x x x -><或,则()V x •是负定的,因此,1212x x <是12x x 和的约束条件; 计算得到()V x 为：()V x 是正定的,因此在121211202x x x x -><即范围内,0e x =是渐进稳定的;现代控制理论第五章习题答案5-1已知系统状态方程为：试设计一状态反馈阵使闭环系统极点配置为-1,-2,-3; 解：依题意有：2011012112M bAbA b ⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥⎣⎦3rankM =,系统能控; 系统0(,,)A b C =∑的特征多项式为：则将系统写成能控标准I 型,则有010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦; 引入状态反馈后,系统的状态方程为：()x A bK x bu =++,其中3K ⨯为1矩阵,设[]012K k k k =,则系统(,,)K A bK C =∑的特征多项式为：根据给定的极点值,得到期望特征多项式为：比较*()()f f λλ与各对应项系数,可解得：012599k k k =-=-=-，则有：[]-5-9-9K =;5-3有系统：（1） 画出模拟结构图;（2） 若动态性能不满足要求,可否任意配置极点 （3） 若指定极点为-3,-3,求状态反馈阵; 解1系统模拟结构图如下：2系统采用状态反馈任意配置极点的充要条件是系统0(,,)A b C =∑完全能控; 对于系统0(,,)A b C =∑有： []0111M bAb ⎡⎤==⎢⎥-⎣⎦2rankM =,系统能控,故若系统动态性能不满足要求,可任意配置极点;3系统0(,,)A b C =∑的特征多项式为：则将系统写成能控标准I 型,则有010231x x u ⎡⎤⎡⎤=+⎢⎥⎢⎥--⎣⎦⎣⎦; 引入状态反馈后,系统的状态方程为：()x A bK x bu =++,设[]01K k k =,则系统(,,)KA bK C =∑的特征多项式为：根据给定的极点值,得到期望特征多项式为：比较*()()f f λλ与各对应项系数,可解得：[]017373k k K =-=-=--，; 5-4设系统传递函数为试问能否利用状态反馈将传递函数变成 若有可能,试求出状态反馈K ,并画出系统结构图;解：6522)3)(2)(1()2)(1()(232--+-+=+-++-=s s s s s s s s s s s W由于传递函数无零极点对消,因此系统为能控且能观; 能控标准I 型为 令[] 210k k k K =为状态反馈阵,则闭环系统的特征多项式为由于状态反馈不改变系统的零点,根据题意,配置极点应为-2,-2,-3,得期望特征多项式为比较 )(λf 与 )(*λf 的对应项系数,可得 即[]52118---=K 系统结构图如下：5-5使判断下列系统通过状态反馈能否镇定;11222 A 011,01011b ---⎡⎤⎡⎤⎢⎥⎢⎥=-=⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦解：系统的能控阵为：2240010115M bAbA b -⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥-⎣⎦3rankM =,系统能控; 由定理 5.2.1可知,采用状态反馈对系统0(,,)A b C =∑任意配置极点的充要条件是(,,)A b C =∑完全能控;又由于3rankM =,系统0(,,)A b C =∑能控,可以采用状态反馈将系统的极点配置在根平面的左侧,使闭环系统镇定; 5-7设计一个前馈补偿器,使系统 解耦,且解耦后的极点为1,1,2,2----; 解：0()()() d W s W s W s = 5-10已知系统：试设计一个状态观测器,使观测器的极点为-r,-2rr>0;解：因为1001c N cA ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎣⎦满秩,系统能观,可构造观测器; 系统特征多项式为[]21det det 0I A λλλλ-⎡⎤-==⎢⎥⎣⎦,所以有10010,0,10a a L ⎡⎤===⎢⎥⎣⎦ 于是11001100x T ATx T bu x u --⎡⎤⎡⎤=+=+⎢⎥⎢⎥⎣⎦⎣⎦ 引入反馈阵12g G g ⎡⎤=⎢⎥⎣⎦,使得观测器特征多项式：根据期望极点得期望特征式：比较()f λ与()*f λ各项系数得：即223r G r ⎡⎤=⎢⎥⎣⎦,反变换到x 状态下2201321023r r G TG r r ⎡⎤⎡⎤⎡⎤===⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ 观测器方程为：。
现代控制理论课后习题及答案
《现代控制理论》课后习题及答案第一章控制系统的状态空间表达式1-1.试求图1-1系统的模拟结构图，并建立其状态空间表达式。
图1-27系统方块结构图图1-1 系统结构方块图解：系统的模拟结构图如下：图1-30双输入--双输出系统模拟结构图图1-2 双输入—双输出系统模拟结构图系统的状态方程如下：u K K x K K x K K x X K x K x x x x J K x J x J K x J K x x J K x x x pp p p n p b1611166131534615141313322211+--=+-==++--===••••••令y s =)(θ，则1x y =所以，系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡••••••6543211654321111111126543210000010000000000000010010000000000010x x x x x x y uK K x x x x x x K K K K K K J K J J K J K J K x x x x x x p p pp n p b1-2.有电路如图1-3所示。
以电压)(t u 为输入量，求以电感中的电流和电容上的电压作为状态变量的状态方程，和以电阻2R 上的电压作为输出量的输出方程。
U图1-28 电路图图1-3 电路图解：由图，令32211,,x u x i x i c===，输出量22x R y =有电路原理可知：•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000*********x x x R y u L x x x CCL L R L L R x x x 。