现代控制理论讲义整理ppt课件.ppt

现代控制理论,控制系统的状态空间分析与综合,2,引 论,经典控制理论: 数学模型:线性定常高阶微分方程和传递函数; 分析方法: 时域法(低阶13阶) 根轨迹法 频域法 适应领域:单输入单输出（SISO）线性定常系统 缺 点:只能反映输入输出间的外部特性，难以揭示系统内部的结构和运行状态。现代控制理论： 数学模型:以一阶微分方程组成差分方程组表示的动态方程 分析方法:精准的时域分析法 适应领域：（1）多输入多输出系统（MIMO、SISO、MISO、SIMO） （2）非线性系统 （3）时变系统 优越性：（1）能描述系统内部的运行状态 （2）便于考虑初始条件（与传递函数比较） （3）适用于多变量、非线性、时变等复杂大型控制系统 （4）便于计算机分析与计算 （5）便于性能的最优化设计与控制 内容：线性系统理论、最优控制、最优估计、系统辨识、自适应控制,近似分析,3,第一章 控制系统的状态空间描述,第二章 线性系统的运动分析,第三章 控制系统的李雅普诺夫稳定性分析,第四章 线性系统的可控性和可观测性,第五章 线性系统非奇异线性变换及系统的规范分解,第六章 线性定常控制系统的综合分析,4,1.1 系统数学描述的两种基本方法 1.2 状态空间描述常用的基本概念1.3 系统的传递函数矩阵 1.4 线性定常系统动态方程的建立,第一章 控制系统的状态空间,5,典 型 控 制 系 统 方 框 图,1.1 系统数学描述的两种基本方法,6,典型控制系统由被控对象、传感器、执行器和控制器组成。 被控过程具有若干输入端和输出端。 数学描述方法： 输入输出描述（外部描述）:高阶微分方程、传递函数矩阵。 状态空间描述（内部描述）：基于系统内部结构，是对系统的一种完整的描述。,7,输入：外部对系统的作用（激励）； 控制：人为施加的激励； 输入分控制与干扰。输出：系统的被控量或从外部测量到的系统信息 。若输出是由传感器测量得到的，又称为观测。状态、状态变量和状态向量 ：能完整描述和唯一确定系统时域行为或运行过程的一组独立（数目最小）的变量称为系统的状态；其中的各个变量称为状态变量。当状态表示成以各状态变量为分量组成的向量时，称为状态向量。状态空间：以状态向量的各个分量作为坐标轴所组成的n维空间称为状态空间。状态轨线：系统在某个时刻的状态，在状态空间可以看作是一个点。随着时间的推移，系统状态不断变化，并在状态空间中描述出一条轨迹，这种轨迹称为状态轨线或状态轨迹。 状态方程：描述系统状态变量与输入变量之间关系的一阶向量微分或差分方程称为系统的状态方程，它不含输入的微积分项。一般情况下，状态方程既是非线性的，又是时变的，可以表示为 输出方程：描述系统输出变量与系统状态变量和输入变量之间函数关系的代数方程称为输出方程，当输出由传感器得到时，又称为观测方程。输出方程的一般形式为动态方程：状态方程与输出方程的组合称为动态方程，又称为状态空间表达式 。一般形式为,1.2 状态空间描述常用的基本概念,8,或离散形式,线性系统：线性系统的状态方程是一阶向量线性微分或差分方程，输出方程是向量代数方程。线性连续时间系统动态方程的一般形式为线性定常系统：线性系统的A，B，C，D或G，H，C，D中的各元素全部是常数。即,或离散形式,若有,9,分别写出状态矩阵 A、控制矩阵 B、输出矩阵 C、前馈矩阵 D :,已知：,为书写方便，常把连续系统和离散系统分别简记为S(A,B,C,D)和S(G,H,C,D)。,线性系统的结构图 ：线性系统的动态方程常用结构图表示。,图中，I为( )单位矩阵，s是拉普拉斯算子，z为单位延时算子。,10,讨论： 1、状态变量的独立性。 2、由于状态变量的选取不是唯一的，因此状态方程、输出方程、动态方程也都不是唯一的。但是，用独立变量所描述的系统的维数应该是唯一的，与状态变量的选取方法无关。 3、动态方程对于系统的描述是充分的和完整的，即系统中的任何一个变量均可用状态方程和输出方程来描述。 例11 试确定图8-5中（a）、（b）所示电路的独立状态变量。图中u、i分别是是输入电压和输入电流，y为输出电压，xi为电容器电压或电感器电流。,解 并非所有电路中的电容器电压和电感器电流都是独立变量。对图8-5（a），不失一般性，假定电容器初始电压值均为0，有,11,因此，只有一个变量是独立的，状态变量只能选其中一个，即用其中的任意一个变量作为状态变量便可以确定该电路的行为。实际上，三个串并联的电容可以等效为一个电容。 对图（b） x1 = x2，因此两者相关，电路只有两个变量是独立的，即（x1和x3）或(x2和x3)，可以任用其中一组变量如（x2，x3）作为状态变量。,12,令初始条件为零，对线性定常系统的动态方程进行拉氏变换，可以得到,系统的传递函数矩阵（简称传递矩阵）定义为,例1-2 已知系统动态方程为,试求系统的传递函数矩阵。解 已知,故,1.3 系统的传递函数矩阵,13,1.4 .1 由物理模型建动态方程根据系统物理模型建立动态方程,1.4 线性定常系统动态方程的建立,RLC 电路,例1-3 试列写如图所示RLC的电路方程，选择几组状态变量并建立相应的动态方程，并就所选状态变量间的关系进行讨论。解 有明确物理意义的常用变量主要有：电流、电阻器电压、电容器的电压与电荷、电感器的电压与磁通。根据独立性要求，电阻器的电压与电流、电容器的电压与电荷、电感器的电流与磁通这三组变量不能选作为系统的状态。 根据回路电压定律,电路输出量 y 为,1)设状态变量为电感器电流和电容器电压，即 则状态方程为,输出方程为,14,其向量-矩阵形式为,简记为,式中，,2)设状态变量为电容器电流和电荷，即 则有,3)设状态变量 （ 无明确意义的物理量），可以推出,15,其向量-矩阵形式为,可见对同一系统，状态变量的选择不具有唯一性，动态方程也不是唯一的。,例1-4 由质量块、弹簧、阻尼器组成的双输入三输出机械位移系统如图所示，具有力F和阻尼器气缸速度V 两种外作用，输出量为质量块的位移，速度和加速度。试列写该系统的动态方程。 分别为质量、弹簧刚度、阻尼系数；x为质量块位移。,双输入三输出机械位移系统,解 根据牛顿力学可知，系统所受外力F与惯性力m 、阻尼力f( V )和弹簧恢复力 构成平衡关系，系统微分方程如下: 这是一个二阶系统，若已知质量块的初始位移和初始速度，系统在输入作用下的解便可唯一确定，故选择质量块的位移和速度作为状态变量。设 。由题意知系统有三个输出量，设,16,于是由系统微分方程可以导出系统状态方程,其向量-矩阵形式为,1.4.2 由高阶微分方程建动态方程1) 微分方程不含输入量的导数项 ：,选n个状态变量为 有,得到动态方程,17,式中,系统的状态变量图,2) 微分方程输入量中含有导数项 ：,一般输入导数项的次数小于或等于系统的阶数n。首先研究情况，为了避免在状态方程中出现输入导数项，可按如下规则选择一组状态变量，设,例15,18,其展开式为,式中， 是n个待定常数。是n个。,由上式的第一个方程可得输出方程是n个。,其余（n）个状态方程如下 n个。,对式求导，有 ：,19,由展开式将 均以 及 u 的各阶导数表示，经整理可得,令上式中 u 的各阶导数的系数为零，可确定各 h 值,记,故,则系统的动态方程为,式中,20,若输入量中仅含次导数且 ，可将高于次导数项的系数置0，仍可应用上述公式。,1.4 .3 由系统传递函数建立动态方程,应用综合除法有,式中， 是直接联系输入、输出量的前馈系数，当G(s)的分母次数大于分子次数时， ， 是严格有理真分式，其分子各次项的系数分别为,下面介绍由 导出几种标准型动态方程的方法：1） 串联分解 如图，取z为中间变量，将 分解为相串联的两部分，有,选取状态变量,21,则状态方程为,输出方程为,其向量-矩阵形式,式中，,当 具有以上形状时， 阵称为友矩阵，相应的状态方程则称为可控标准型。,时， 的形式不变，,22,当 时， 不变，,当 时，若按下式选取状态变量,式中，T为转置符号，则有,注意 的形状特征。若动态方程中的 具有这种形式，则称为可观测标准型。自行证明：可控标准型和可观测标准型是同一传递函数的不同实现。可控标准型和可观测标准型的状态变量图如图 ：,（对偶关系 ）,可控标准型状态变量图,可观测标准型状态变量图,23,例1-6 设二阶系统微分方程为 ，试列写可控标准型、可观测标准型动态方程，并分别确定状态变量与输入，输出量的关系。 解 系统的传递函数为,于是，可控标准型动态方程的各矩阵为,由G(s)串联分解并引入中间变量z有,对y求导并考虑上述关系式，则有,令 可导出状态变量与输入，输出量的关系；,可观测标准型动态方程中各矩阵为,24,状态变量与输入，输出量的关系为,该系统的可控标准型与可观测标准型的状态变量图 ：,（a）可控标准型实现 （b）可观测标准型实现,2） 只含单实极点时的情况 当 只含单实极点时，动态方程除了可化为可控标准型或可观测标准型以外，还可化为对角型动态方程，其A阵是一个对角阵。设D(s)可分解为 D(s)= 式中， 为系统的单实极点，则传递函数可展成部分分式之和,25,而 ，为 在极点 处的留数，且有Y(s)= U(s),若令状态变量 其反变换结果为,展开得,其向量-矩阵形式为 （其状态变量如图(a)所示 ）,26,若令状态变量则 Y(s)=,进行反变换并展开有,其向量-矩阵形式为,其状态变量图如图(b)所示 ，两者存在对偶关系 对角型动态方程状态变量图 如下：,27,（a） （b）,对角型动态方程状态变量图,3） 含重实极点时的情况 当传递函数除含单实极点之外还含有重实极点时，不仅可化为可控标准型或可观测标准型，还可化为约当标准型动态方程，其A阵是一个含约当块的矩阵。设D(s)可分解为 D(s)= 式中 为三重实极点， 为单实极点，则传递函数可展成为下列部分分式之和：,28,其状态变量的选取方法与之含单实极点时相同，可分别得出向量-矩阵形式的动态方程：,29,其对应的状态变量图如图（a），（b）所示。上面两式也存在对偶关系。,约当型动态方程状态变量图,30,1.4 .4 由差分方程和脉冲传递函数建立离散动态方程单输入-单输出线性定常离散系统差分方程的一般形式为：,两端取z变换并整理得,G(z)称为脉冲传递函数 ，利用z变换关系 和 ，可以得到动态方程为：,简记为,31,1.4 .5 由传递函数矩阵建动态方程 （传递函数矩阵的实现） 给定一传递函数矩阵G(s)，若有一系统（A，B，C，D）能使 成立，则称系统（A，B，C，D）是G(s)的一个实现。这里仅限于单输入-多输出和多输入-单输出系统。SIMO系统的实现：,单输入多输出系统结构图,1）系统可看作由q个独立子系统组成，传递矩阵为：,32,式中，d为常数向量； 为不可约分的严格有理真分式（即分母阶次大于分子阶次）函数。通常 , , 的特性并不相同，具有不同的分母，设最小公分母为：,的一般形式为,将 作串联分解并引入中间变量Z，令若将A阵写为友矩阵，便可得到可控标准型实现的状态方程：,每个子系统的输出方程 ：,33,每个子系统的输出方程 ：,可以看到，单输入，q维输出系统的输入矩阵为q维列向量，输出矩阵为（q n）矩阵，故不存在其对偶形式，即不存在可观测标准型实现。MISO系统的实现:,多输入单输出系统结构图,系统由p个独立子系统组成，系统输出由子系统输出合成为 ：,34,式中,同理设 , , 的最小公分母为D（s），则,若将A阵写成友矩阵的转置形式，便可得到可观测标准型实现的动态方程 ：,35,可见，p维输入，单输入系统的输入矩阵为（n p）矩阵输出矩阵为一行矩阵，故不存在其对偶形式，即不存在可控标准型实现。,例1-7 已知单输入-多输出系统的传递函数矩阵为 ，求其传递 矩阵的可控标准型实现及对角型实现。,例1-7 已知单输入-多输出系统的传递函数矩阵为 ，求其传递 矩阵的可控标准型实现及对角型实现。,解 由于系统是单输入，多输出的，故输入矩阵只有一列，输出矩阵有两行。将 化为 严格有理真分式,各元素的最小公分母D（s）为,故,则可控标准型动态方程为 ：,36,由 可确定系统极点为-1，-2，它们构成对角形状态矩阵的元素。鉴于输入矩阵只有一列，这里不能选取极点的留数来构成输入矩阵，而只能取元素全为1的输入矩阵。于是，对角型实现的状态方程为 ：,其输出矩阵由极点对应的留数组成， 在-1，-2处的留数分别为：,故其输出方程为,37,本章作业：83，84，85，87,38,第二章 线性系统的运动分析,2.1 线性定常连续系统的自由运动2.2 状态转移矩阵的性质2.3 线性定常连续系统的受控运动2.4 线性定常离散系统的分析2.5 连续系统的离散化,39,在控制u=0情况下，线性定常系统由初始条件引起的运动称为线性定常系统的自由运动，可由齐次状态方程描述 ：齐次状态方程求解方法：幂级数法、拉普拉斯变换法和凯莱哈密顿定理法。 幂级数法:设齐次方程的解是t的向量幂级数式中， 都是n维向量，且 ，求导并考虑状态方程，得,2.1 线性定常连续系统的自由运动,等号两边对应的系数相等，有,40,故,定义,则,称为矩阵指数函数，简称矩阵指数 ,又称为状态转移矩阵，记为 ： 求解齐次状态方程的问题，核心就是计算状态转移矩阵的问题 。拉普拉斯变换法 ：,对 进行拉氏变换，有： 进行拉氏反变换，有： 与 相比有：它是 的闭合形式。,例 2-1 设系统状态方程为 ，试用拉氏变换求解。,解,41,状态方程的解为 :,凯莱哈密顿定理 矩阵A满足它自己的特征方程。即若设n阶矩阵A的特征多项式为,则有 :,42,从该定理还可导出以下两个推论：推论1 矩阵A的 次幂，可表为A的（n-1）阶多项式 ：,推论2 矩阵指数 可表为A的（n-1）阶多项式，即： 且各作为时间的函数是线性无关的。,在式推论1中用A的特征值替代A后等式仍能满足：,利用上式和k个就可以确定待定系数 ：若 互不相等 ：可写出各所构成的n元一次方程组为 ：,43,求解上式，可求得系数 ， ， ，它们都是时间t的函数，将其代入推论2式后即可得出 。,例2-2 已知 ，求 。,解 首先求A的特征值：,将其代入 ，有：,44,若矩阵 A 的特征值是 m 阶的： 则求解各系数的方程组的前m个方程可以写成：,其它由 组成的（k - m）个方程仍与第一种情况相同，它们上式联立即可解出各待定系数。,45,例2-3 已知 ，求 。,解 先求矩阵 A 的特征值，由得：,46,2.2 状态转移矩阵的性质,状态转移矩阵 具有如下运算性质：,1),2),3),4),表明 与 可交换，且,在式 3）中，令 便可证明；表明 可分解为 的乘积，且 是可交换的。,证明:由性质3）有,根据 的这一性质，,对于线性定常系统，显然有,5),证明 ：由于,则,即由,转移至,的状态转移矩阵为,47,6),证明：由,和,得到,7),8) 若,，则,证明：,例2-4 已知状态转移矩阵为,，试求,。,解：根据状态转移矩阵的运算性质有,9) 若,，则,48,2.3 线性定常连续系统的受控运动,线性定常系统的受控运动: 线性定常系统在控制作用下的运动，数学描述为：,主要有如下两种解法：,1) 积分法 由上式,由于,积分后有,即,式中，第一项为零输入响应；第二项是零状态响应。通过变量代换，上式又可表示为：,若取 作为初始时刻，则有,49,2) 拉普拉斯变换法 将 式两端取拉氏变换，有,进行拉氏反变换有,例2-5 设系统状态方程为,且,试求在,作用下状态方程的解。,解 由于,前面已求得,50,51,2.4 线性定常离散系统的分析,1)递推法(线性定常系统) 重写系统的动态方程如下：,令状态方程中的k=0，1，k-1，可得到T，2T，kT 时刻的状态，即：,k=0:,k=2:,k=1:,k=k-1:,于是，系统解为:,52,2.5 连续系统的离散化,2.5.1 线性定常连续系统的离散化,已知线性定常连续系统状态方程,在,及 作用下的解为:,令,，则,；令,则,并假定在 区间内， ，于是其解化为,若记,变量代换得到,故离散化状态方程为,式中，,与连续状态转移矩阵,的关系为,53,2.5.2 非线性时变系统的离散化及分析方法,对于非线性时变系统，常采用近似的离散化处理方法。当采样周期T足够小时，按导数定义有,代入（8-5a）得到离散化状态方程,对于非线性时变系统，一般都是先离散化，然后再用递推计算求数值解的方法进行系统的运动分析。,本章作业：88，89，811,54,3.1 李雅普诺夫稳定性概念 3.2 李雅普诺夫稳定性间接判别法 3.3 李雅普诺夫稳定性直接判别法 3.4 线性定常系统的李雅普诺夫稳定性分析,第三章 控制系统的李雅普诺夫稳定性分析,55,如果对于所有t，满足 的状态 称为平衡状态（平衡点）。,1) 平衡状态:,3.1 李雅普诺夫稳定性概念,平衡状态的各分量不再随时间变化；若已知状态方程，令 所求得的解 x ,便是平衡状态。,（1）只有状态稳定，输出必然稳定； （2）稳定性与输入无关。,2) 李雅普诺夫稳定性定义：,如果对于任意小的 0，均存在一个 ，当初始状态满足 时，系统运动轨迹满足lim ，则称该平衡状态xe 是李雅普诺夫意义下稳定的，简称是稳定的。 表示状态空间中x0点至xe点之间的距离，其数学表达式为：,3) 一致稳定性：,通常与、t0 都有关。如果与t0 无关，则称平衡状态是一致稳定的。定常系统的与t0 无关，因此定常系统如果稳定，则一定是一致稳定的。,56,4）渐近稳定性：,系统的平衡状态不仅具有李雅普若夫意义下的稳定性，且有：,称此平衡状态是渐近稳定的。,5）大范围稳定性：,当初始条件扩展至整个状态空间，且具有稳定性时，称此平衡状态是大范围稳定的，或全局稳定的。 此时 。,6）不稳定性 ：,不论取得得多么小，只要在 内有一条从x0 出发的轨迹跨出 ，则称此平衡状态是不稳定的。,注意：按李雅普诺夫意义下的稳定性定义，当系统作不衰减的振荡运动时则认为是稳定的，同经典控制理论中的稳定性定义是有差异的。经典控制理论的稳定是李雅普诺夫意义下的一致渐近稳定。,57,稳定性定义的平面几何表示,设系统初始状态 x0 位于平衡状态 xe 为球心、半径为的闭球域内，如果系统稳定，则状态方程的解在的过程中，都位于以 xe 为球心，半径为的闭球域内。,（a）李雅普诺夫意义下的稳定性 （b）渐近稳定性 （c）不稳定性,58,3.2 李雅普诺夫稳定性间接判别法,李雅普诺夫第一法（间接法） 是利用状态方程解的特性来判断系统稳定性的方法，它适用于线性定常、线性时变及可线性化的非线性系统。 线性定常系统的特征值判据 系统 渐近稳定的充要条件是：系统矩阵A的全部特征值位于复平面左半部，即证明：(略),59,李雅普诺夫第二法（直接法）基本原理 ：根据物理学原理，若系统贮存的能量（含动能与位能）随时间推移而衰减，系统迟早会到达平衡状态。 实际系统的能量函数表达式相当难找，因此李雅普诺夫引入了广义能量函数，称之为李雅普诺夫函数。它与 及t 有关，是一个标量函数，记以 ；若不显含t ，则记以 。 考虑到能量总大于零，故为正定函数。能量衰减特性用 或 表示。 实践表明，对于大多数系统，可先尝试用二次型函数 作为李雅普诺夫函数。,3.3 李雅普诺夫稳定性直接判别法,60,3.3.1 标量函数定号性,正定性：标量函数 在域S中对所有非零状态 有 且 ，则称 均在域S内正定。如 是正定的。负定性：标量函数 在域S中对所有非零x有 且 ，则称 在域S内负定。如 是负定的。 如果 是负定的，则 一定是正定的。负（正）半定性： ，且 在域S内某些状态处有 ，而其它状态处均有 （ ），则称 在域S内负（正）半定。 设 为负半定，则 为正半定。如 为正半定不定性: 在域S内可正可负，则称 不定。如 是不定的。,二次型函数 是一类重要的标量函数，记,其中，P 为对称矩阵，有 。,61,当的各顺序主子行列式均大于零时 ，即,则 正定，且称 P为正定矩阵。当 P的各顺序主子行列式负、正相间时，即,则 负定，且称 P为负定矩阵。若主子行列式含有等于零的情况，则 为正半定或负半定。不属以上所有情况的 不定。,62,设系统状态方程为 ，其平衡状态满足 ,不失一般性地把状态空间原点作为平衡状态,并设在原点邻域存在 对 x 的连续一阶偏导数。,3.3.2 李雅普诺夫第二法诸稳定性定理,定理1 若（1）,正定，（2）,负定；则原点是渐近稳定的。,负定表示能量随时间连续单调地衰减，故与渐近稳定性定义叙述一致。,定理2 若（1）,正定；（2）,负半定，且在非零状态不恒为零；则原点是渐近,稳定的。,负半定表示在非零状态存在,，但在从初态出发的轨迹,上，不存在,的情况，于是系统将继续运行至原点。状态轨迹仅是经历能量不变的状态，而不,会维持在该状态。,定理3 若（1）,正定；（2）,负半定，且在非零状态恒为零；则原点是李雅普,，表示系统能维持等能量水平运行，使系统维持在非零状态,沿状态轨迹能维持,诺夫意义下稳定的。,而不运行至原点。,定理4 若（1）,正定；（2）,正定；则原点是不稳定的。,正定表示能量函数随时间增大，故状态轨迹在原点邻域发散。,正定，当,正半定，且在非零状态不恒为零时，则原点不稳,参考定理2可推论：,定。,63,注意：李雅普诺夫第二法诸稳定性定理所述条件都是充分条件。,具体分析时，先构造一个李雅普诺夫函数,，通常选二次型函数，求其导数,再将状态方程代入，最后根据,是否有恒为零：令,将状态方程代入，若能导出非零解，表示对,，,若导出的是全零解，表示只有原点满足,的条件。,的定号性判别稳定性。,的条件是成立的；,例3-1 试用李雅普诺夫第二法判断下列非线性系统的稳定性。,解 令,及,，可以解得原点（,）是系统的唯一平衡状态。,，则,将状态方程代入有,显然,负定，根据定理1，原点是渐近稳定的。鉴于只有一个平衡状态，该非线性,与t 无关，系统大范围一致渐近稳定。,取李雅普诺夫函数为,系统是大范围渐近稳定的。因,判断在非零状态下,64,例3-2 试判断下列线性系统平衡状态的稳定性。,，,解 令,得知原点是唯一的平衡状态。选,则,当,时，,；当,时，,故,不定，不能对稳定性作出判断，应重选,选,，则考虑状态方程后得,对于非零状态（如,）存在,，对于其余非零状态，,，故,根据定理2，原点是渐近稳定的，且是大范围一致渐近稳定。,负半定。,例3-3 试判断下列线性系统平衡状态的稳定性。,，,解 由,可知原点是唯一平衡状态。选,，考虑状态方程则有,对所有状态，,，故系统是李雅普诺夫意义下稳定的。,65,例3-4 试判断下列线性系统平衡状态的稳定性。,解 原点是唯一平衡状态。选,，则,，,与,故存在非零状态（如,使,而对其余任意状态有,，故,根据定理4的推论，系统不稳定。,无关，,),正半定。,解,是系统的唯一平衡状态，方程中的常数项可以看作是阶跃输入作用的,，,得到,原状态方程在,状态空间（1，1）处稳定性判别问题就变成变换后状态方程在 X,对其求导考虑状态方程得到,系统原点是大范围一致渐近稳定的，因而原系统在平衡状态（1，1）处是大,结果。作坐标变换,选,状态空间原点处稳定性的判别问题。,围一致渐近稳定的。,注意：一般不能用李雅普诺夫函数去直接判别非原点的平衡状态稳定性。,例3-5 试判断下列线性系统平衡状态的稳定性。,66,例3-6 试判断下列非线性系统平衡状态的稳定性。,解 这实际上是一个可线性化的非线性系统的典型例子。令,得知系统有两个平衡状态，,和,对位于原点的平衡状态，选,于是，当,时，系统在原点处的平衡状态是局部,根据定理4，当,时原点显然是不稳定的,时原点也是不稳定的,从状态方程直接看出。,，作坐标变换,，得到新的状态方程,因此，通过与原状态方程对比可以断定：对于原系统在状态空间,处的平衡状态，当,时是局部一致渐近稳定的；当,时是不稳定的，,时也是不稳定的。,一致渐近稳定的。,或系统发散，,也可以,当,对于平衡状态,当,有,67,3.4 线性定常系统的李雅普诺夫稳定性分析,3.4.1 连续系统渐近稳定的判别,设系统状态方程为,，,平衡状态。可以取下列正定二次型函数,作为李雅普诺夫函数,根据定理1，只要,正定（即,负定）则系统是大范围一致渐近稳定的。于是线性,，存在满足 式的,为非奇异矩阵，故原点是唯一,求导并考虑状态方程,令,得到,定常连续系统渐近稳定的判定条件可表示为：给定一正定矩阵,正定矩阵,。,(#),(#),先指定正定的,阵，然后验证,阵是否正定。,注：,(),68,定理5 （证明从略）系统,渐近稳定的充要条件为：给定正定实对称矩阵,正定实对称矩阵,使 式成立。,，存在,该定理为系统的渐近稳定性判断带来实用上的极大方便。,(),例3-7 试用李雅普诺夫方程确定使图所示系统渐近稳定的,值范围。,例3-7 系统框图,解 由图示状态变量列写状态方程,稳定性与输入无关，可令,。由于,，,非奇异，原点为唯一的平衡状,为正半定矩阵,态。取,则,，,负半定。令,，有,，考虑状态方程中,，解得,；考虑到,，解得,，表明唯有原点存在,69,令,展开的代数方程为6个，即,，,，,，,，,解得,使,正定的条件为：,及,。故,时，系统渐近稳定。由于是线性定,常系统，系统大范围一致渐近稳定。,70,3.4.2 离散系统渐近稳定的判别,设系统状态方程为,，式中,以,代替,，有,阵非奇异，原点,考虑状态方程，有,是系统的一个李雅普诺夫函数，于是,式称为李雅普诺夫代数方程。,定理7 系统,渐近稳定的充要条件是：给定任一正定实对称矩阵,（常,），存在正定对称矩阵,，使 式成立。,令,取正定二次型函数,是平衡状态。,(#),(#),(#),71,本章作业：814，815,72,可控性和可观测性的概念,线性系统的可控性和可观测性,线性定常系统的可控性,线性定常系统的可观测性,可控性,可观测性与传递函数矩阵的关系,返回,连续系统离散化后的可控性与可观测性,73,4.1 可控性和可观测性的概念,可控性,如果系统所有状态变量的运动都可以通过有限点的控制输入来使其由任意的初态达到任意设定的终态，则称系统是可控的，更确切的说是状态可控的；否则，就称系统是不完全可控的，简称为系统不可控。,可观性,如果系统所有的状态变量任意形式的运动均可由有限点的输出测量完全确定出来，则称系统是可观测的，简称为系统可观测；反之，则称系统是不完全可观测的，简称为系统不可观测。,可控性与可观测性的概念，是用状态空间描述系统引伸出来的新概念，在现代控制理 论中起着重要的作用。可控性、可观测性与稳定性是现代控制系统的三大基本特性。,第四章 线性系统的可控性和可观测性,74,下面举几个例子直观地说明系统的可控性和可观测性。,上图所示的结构图，其中左图显见,受,的控制，但,与,无关，故系统不可,，但,是受,影响的，,能间接获得,中图中的,、,均受,的控制，故系统可控，但,与,中的,、,均受u的控制，且在,中均能观测到,、,故系统是可控可观测的。,控。系统输出量,的信息，故系统是可观测的。,无关，故系统不可观测。,又图,只有少数简单的系统可以从结构图或信号流图直接判别系统的可控性与可观测性，如果系统结构、参数复杂，就只能借助于数学方法进行分析与研究，才能得到正确的结论。,75,4.2 线性定常系统的可控性,可控性分为状态可控性和输出可控性，若不特别指明，一般指状态可控性。状态可控性只与状态方程有关，与输出方程无关。, 4.2.1离散系统的可控性,(1) 单输入离散系统,为导出系统可控性的条件，设单输入系统状态方程为,定义,其解为,由于,和,取值都可以是任意的，因此,的取值也可以是任意的。,76,记,称,为可控性矩阵。,个方程中有,个未知数,称为可控性判据。此为充要条件。,当rank S1n时，系统不可控，表示不存在能使任意,转移至任意,的控制。,(4-1),或,则,(4-2),式(4-1)是一个非齐次线性方程组，,77,从以上推导看出，状态可控性取决于,和 ，当,不受约束时，可控系统的状态转移,个采样周期便可以完成，有时状态转移过程还可能少于,上述过程不仅导出了单输入离散系统可控性条件，而且 还给出了求取控制指,过程至多以,个采样周期。,令的具体方法。, 4.2.1 多输入离散系统,设系统状态方程为,可控性矩阵为,多输入线性定常散离系统状态可控的充分必要条件是,或,(4-1),78,的行数总小于列数，在列写,时，若能知道,的秩为,，便不必把,和列写出来。,阶行列式,多输入线性定常离散系统,转移过程一般可少于,个采样,周期。,例8-30 设单输入线性定常散离系统状态方程为,试判断可控性；若初始状态,，确定使,的控制序列,；研究使,的可能性。,解 由题意知,故该系统可控。,技巧：,便可确定可控性。,（2）利用,计算一次,（1）,的其余列,都计算,79,可按式（8-90）求出,令k=0，1，2，可得状态序列,。为了避免矩阵求逆，下面用递推法来求。,令,，即解下列方程组,80,其系数矩阵即可控性矩阵S1，它的非奇异性可给出如下的解,若令,，即解下列方程组,容易看出其系数矩阵的秩为2，但增广矩阵,两个秩不等，方程组无解，意为不能在第二个采样周期内使给定初态转移至原点。若,的秩为3，,该两个秩相等时，便意味着可用两步完成状态转移。,81,例8-31 输入线性定常离散系统的状态方程为,试判断可控性，设初始状态为,，研究使,的可能性。,解：,由前三列组成的矩阵的行列式不为零，故该系统可控，一定能求得控制使,给出,系统从任意初态在三步内转移到原点。由,设初始状态为,82,由于,可求得,，在一步内使该初态转移到原点。当初始状态为,亦然，只是,。但本例不能对任意初态，使之在一步内转移到原点。,时，, 4.2.1 连续系统的可控性（1）单输入系统,，如果存在无约束的分段连续控制函数,从任意初态,转移至任意终态,，则称该系统是状态完全可控的，简称是可控的。,间间隔内,设状态方程为,定义,终态解为,显然，,的取值也是任意的。于是有,，能使系统,定义：在有限时,83,利用凯莱-哈密顿定理的推论,有,令,则有,记,其状态可控的充分必要条件是,（2）多输入系统,记可控性矩阵,状态可控的充要条件为,或,84,例8-32 试用可控性判据判断图8-20所示桥式电路的可控性。解 选取状态变量：,。电路的状态方程如下：,可控性矩阵为,当,时，,，系统可控；反之当,，即电桥处于平衡状态时，,，系统不可控，显然， 不能控制 。,85,图8-20 电桥电路 图8-21并联电路,例8-33 试用可控性判断图8-21并联网络的可控性。,解 网络的微分方程为,式中，,状态方程为,于是,当,时，系统可控。当,，,，有,，,系统不可控；实际上，设初始状态,， 只能使,，而不能将,与,分别转移到不同的数值，即不能同时控制住两个状态。,，,86,例8-34 判断下列状态方程的可控性,解,显见S4矩阵的第二、三行元素绝对值相同，,（3）A为对角阵或约当阵时的可控性判据,，系统不可控。,设二阶系数A、b矩阵为,其可控性矩阵S3的行列式为,由此可知：A阵对角化且有相异元素时，只需根据输入矩阵没有全零行即可判断系统可控。,时，则不能这样判断，这时,，系统总是不可控的。,若,87,又设二阶系数A、b矩阵为,其可控性矩阵S3的行列式为,矩阵中与约当块最后一行所对应的行不是全零,由此可知：当A阵约当化且相同,矩阵中的其它行是否为零行是无关的。,以上判断方法可推广到A阵对角化、约当化的n阶系统。设系统状态方程为,A为对角阵时的可控性判据可表为：A为对角阵且元素各异时，,输入矩阵不存在全零行。,特征值分布在一个约当快时，只需根据输入,行，即可判断系统可控，与输入,88,当A为对角阵且含有相同元素时，上述判据不适用，应根据可控性矩阵的秩来判断。设系统状态方程为,全零行（与约当块其它行所对应的行允许是全零行）；输入矩阵中与相异特征值所对应的行不存在全零行。,A阵约当化时的可控性判据可表为：输入矩阵中与约,当A阵的相同特征值分布在两个或更多个约当块时，例如,适用，也应根据可控性矩阵的秩来判断。,，以上判据不,当块最后一行所对应的行不存在,89,例8-35 下列系统是可控。1),2),3),例8-36 下列系统不可控1),2),3),90,（4）可控标准型问题,其可控性矩阵为,与该状态方程对应的可控性矩阵,一定是可控的，这就是式(4-3)称为可控标准型的由来。,是一个右下三角阵，且其主对角线元素均为1，系统,(4-3),91,4.3 线性定常系统的可观测性,4.3.1离散系统的状态可观测性,及,，则称系,因为是讨论可观性，可假设输入为0，其解为,将,写成展开式,定义：已知输入向量序列,输出向量序列,，能唯一确,确定任意初始状态向量,统是完全可观测的。,92,其向量-矩阵形式为,令,为线性定常离散系统可观测性矩阵。,可观测的充分必要条件为,93,例8-37 判断下列线性定常离散系统的可观测性，并讨论可观测性的物理解释。其输出矩阵取了两种情况。,解 计算可观测性矩阵V1 （1）,故系统可观测。由输出方程,由于,可见，在第k步便可由输出确定状态变量 .,故在第（k+1）步便可确定,。由于,故在第（k+2）步便可确定,该系统为三阶系统，可观测意味着至多以三步便能由y（k），y（k+1），y（k+2）的输出测量值来确定三个状态变量。,。,94,（2）,故系统不可观测。由输出方程,可看出三步的输出测量值中始终不含 ，故 是不可观测状态变量。只要有一个状态变量不可观测，称系统不完全可观测，简称不可观测。,连续系统的状态可观测性 其定义为：已知输入,及有限时间间隔,到的输出,，能唯一确定初始状态,，则称系统是完全可观测的，简称系统可观测。,内测量,95,4.3.2,连续系统的可观测性,定义,已知输入u(t)及有限时间间隔,对于多输入系统,状态可观测的充分必要条件是,或,均称为可观测性矩阵。,96, 4.3.3 A为对角阵或约当阵时的可观测性判据,（1）单输入对角二阶系统,可观测矩阵,的行列式为,判据：A阵对角化且有相异特征值时，只需根据输出矩阵中没有全零列即可判断系统,时，则不能这样判断，这时,，系统总是不可观测的。,可观测。若,（2）单输入约当二阶系统,则,97,有时A阵的相同特征值分布在两个或更多个约当块内时，例如,，以上判断方法不适用。,以下 推广到A阵对角化、约当化的n阶系统。设系统动态方程为（令u=0）,式中,为系统相异特征值，状态变量间解耦，输出解为,判据：输出矩阵中与约当块最前一列所对应的列不是全零列。,98,A为对角阵时可观测判据：可表为： A为对角阵且元素各异时，输出矩阵不,存在全零列。,当A为对角阵但含有相同元素时，上述判据不适用，应根据可观测矩阵的秩来判断。,设系统动态方程为,为二重特征值且构成一个约当块，,，为相异特征值。动态方程解为,99,输出矩阵中与约当块最前一列对应的列不存在全零列（与约当块其它列所对应的列允许是全零列）；输出矩阵中与相异特征值所对应的列不存在全零列。 对于相同特征值分布在两个或更多个约当块内的情况，以上判据不适用，仍应用可观测矩阵来判断。,故A为约当,例8-38 下列系统可观测，试自行说明。,1）,2）,阵且相同特征值分布在一个约当块内时，可观测判据：,100,例8-39 下列系统不可观测，试自行说明。（1）,（2）,4.3.4 可观测标准型问题 动态方程中的A、c矩阵具有下列形式,101,其可观测性矩阵,V2是一个右下三角阵，,，系统一定可观测，这就是形如（8125）所示的A、C,矩阵称为可观测标准型名称的由来。一个可观测系统，当A、C阵不具有可观测标准型时，也可选择适当的变换化为可观测标准型。,102,4.4 可控性、可观测性与传递函数矩阵的关系, 4.4.1 SISO系统,当A阵具有相异特征值,时，通过线性变换定可是A对角化为,利用A阵对角化的可控、可观测性判据可知：当,时，,不可控；当,时，,测。试看传递函数,所具有的相应特点。由于,不可观,其中,（令初始条件