控制系统的能控性和能观测性.ppt
第3章 控制系统的能控性和能观测性,在多变量控制系统中,能控性和能观测性是两个反映控制系统构造的基本特性,是现代控制理论中最重要的基本概念。,问题的提出:这是由于在经典控制理论中,只限于讨论控制作用(输入)对输出的控制。输入与输出这两个量的关系,唯一地由系统的传递函数所确定,只要系统是稳定的,系统就是能控的。另一方面,系统的输出量本身就是被控量,对于一个实际的物理系统来说,它当然是可以观测到的,所以在经典控制理论中没有必要涉及能控性和能观性。然而在现代控制理论中,是把反映系统内部运动状态的状态向量作为被控量,而且它们不一定是实际上可观测到的物理量,至于输出量则是状态向量的线性组合,这就产生了从输入量 到状态量 的能控性问题和从输出量 到状态量 的能观测性问题。,本章的内容为:,1.引言能控性、能观测性的基本概念2.能控性及其判据3.能观测性及其判据4.离散系统的能控性和能观测性5.对偶原理,6.能控标准形和能观测标准形,7.能控性、能观测性与传递函数的关系,8.系统的结构分解,9.实现问题,10.使用MATLAB判断系统的能控性和能观测性,3.1 引言,首先,通过例子介绍能控性、能观测性的基本概念。,例3-1 电路如下图所示。如果选取电容两端的电压 为状态变量,即:。电桥平衡时,不论输入电压 如何改变,不随着 的变化而改变,或者说状态变量不受 的控制。即:该电路的状态是不能控的。,显然,当电桥不平衡时,该电路的状态是能控的。,例3-2 电路如下图所示,如果选择电容C1、C2两端的电压为状态变量,即:,电路的输出 为C2上的电压,即,则电路的系统方程为,系统状态方程的解为,可见,不论加入什么样的输入信号,总是有,一般情况下,系统方程可以表示为,(1),状态能控与否,不仅取决于B 阵(直接关系),还取决于A 阵(间接关系)。,系统状态转移矩阵为,系统能观测问题是研究测量输出变量 y 去确定状态变量的问题。,例3-3 电路如下图所示。选取 为输入量,为输出量,两个电感上的电流分别作为状态变量,则系统方程为,系统状态方程的解为,从上式可知,不论初始状态为什么数值,输出 仅仅取决于其差值。当,则输出恒等于零。显然,无法通过对输出的观测去确定初始状态,称这样的系统是不能观测的。,对于不能观测的系统,其不能观测的状态分量与y 既无直接关系,又无间接关系。状态是否能观测不仅取决于C,还与A 有关。,一般情况下,系统方程如式(1)所示,状态能观测与否,不仅取决于C 阵(直接关系),还取决于A阵(间接关系)。,两个例子的分析结论是:能控与A,B阵有关;能观与A,C阵有关。,3.2 能控性及其判据,3.2.1 线性定常系统的能控性及其判据,1.能控性定义,2)如果在有限时间区间 内,存在容许控制,使系统从状态空间坐标原点推向预先指定的状态,则称系统是状态能达的;由于连续系统的状态转移矩阵是非奇异的,因此系统的能控性和能达性是等价的容许控制(对于一个实际的控制问题,输入控制的u(t)的取值必定要受一定条件的约束。满足约束条件的控制作用u(t)的一个取值对应于r维空间的一个点,所有满足条件的控制作用u(t)的取值构成r维空间的一个集合,记为,称之为容许控制集。凡是属于容许控制集的控制,都是容许控制。)。,3)只有整个状态空间中所有的有限点都是能控的,系统才是能控的。,2.能控性判据,定理3-1(2)式的线性定常系统为状态能控的充分必要条件是下面的nn维格拉姆矩阵满秩,(5),(这个定理为能控性的一般判据,所谓满秩就是每个状态能控。但是,由于要计算状态转移矩阵,比较繁琐。实际上,常用下面介绍的判据。),定理3-2(2)式的线性定常系统为状态能控的充分必要条件是下面的nnr 维能控性矩阵满秩。,(6),(7),于是,(9),如果系统能控,必能够从(9)式中解得,。这样就要求,(本判据本身很简单,因此是最为常用的方法。),定理3-3(PBH判别法)(2)式的线性定常系统为状态能控的充分必要条件是,对A 的所有特征值,都有,(10),(证明略),例3-6 有如下两个线性定常系统,判断其能控性。,(1),(2),解,根据定理3-4,系统(1)不能控;系统(2)能控。,则系统能控的充分必要条件是矩阵 中与每一个约当子块最下面一行对应行的元素不全为零。,(12),例3-7 有如下两个线性定常系统,判断其能控性。,(1),(2),(定理(3-4)、定理(3-5)不仅可以判断系统能控性,而且对于不能控的系统,可以知道哪个状态分量不能控。)说明:1.上面通过几个定理给出判断系统能控性的判据。虽然它们的表达形式、方法不同,但是,在判断线性定常系统能控性时是等价的。2.在线性连续定常系统中,由于能达性和能控性是等价的,因此,能控性判据同样可以判断能达性。,3.2.2 线性时变系统的能控性判据,(13),线性时变系统的状态方程为,定理3-7 状态在时刻 能控的充分必要条件是存在一个有限时间,使得以下格拉姆矩阵非奇异。,(14),(16),当,定理3-8 如果线性时变系统的 和 的元是(n1)阶连续可微的。如果存在一个有限的,使得,(17),则系统在 是能控的。,解,3.3 能观测性判据,3.3.1 线性定常系统能观测性及其判据,1.能观测性定义,(18),线性定常系统方程为,如果在有限时间区间()内,通过观测,能够惟一地确定系统的初始状态,称系统状态在 是能观测的。如果对任意的初始状态都能观测,则称系统是状态完全能观测的。,2)如果根据 内的输出 能够惟一地确定任意指定状态,则称系统是可检测的。连续系统的能观测性和能检测性等价。,),(,x,1,t,3)状态空间中所有有限点都是能观测的,则系统才是能观测的。,4)系统的输入 以及确定性的干扰信号 均不改变系统的能观测性。,(这个定理为能观测性的一般判据。但是,由于要计算状态转移矩阵,比较繁琐。实际上,常用下面介绍的判据。),定理3-10(18)式所描述的系统为能观测的充分必要条件是以下能观性矩阵满秩,即,(21),(22),证明 设,系统的齐次状态方程的解为,(23),应用凯-哈定理,有,则,或者写成,由于 是已知函数,因此,根据有限时间 内的 能够唯一地确定初始状态 的充分必要条件为 满秩。,例3-9 系统方程如下,试判断系统的能控性,(由于以上判据很简单,因此最为常用),定理3-12 如果(18)式描述的系统的A 阵特征值 互异,经过非奇异线性变换成为对角阵,则系统为能观测的充分必要条件是 C 矩阵中不包含元素全为零的列。,例3-10 有如下两个线性定常系统,判断它们的能观测性。,(1),(2),解 根据定理3-12可以判断,系统(1)是不能观测的。系统(2)是能观测的。,则系统能观测的充分必要条件是矩阵 中与每一个约当子块第一列对应的列,其元素不全为零。,例3-11 如下线性定常系统,试判别系统的能观测性。,解 应用定理3-13可知,系统能观测。,(定理(3-12)、定理(3-13)不仅可以判断系统能观测性,而且对于不能观测的系统,可以知道哪个状态分量不能观测。)说明:1.上面通过几个定理给出判断系统能观测性的判据。虽然它们的表达形式、方法不同,但是,在判断线性定常系统能观测性时是等价的。2.在线性连续定常系统中,由于能检测性和能观测性是等价的,因此,能观测性判据同样可以判断能检测性。,3.3.2 线性时变系统的能观测性判据,线性时变系统方程为,(25),定理3-14 状态在时刻 能观测的充分必要条件是存在一个有限时刻,使得函数矩阵 的n个列在 上线性无关。,定理3-15 状态在时刻 能观测的充分必要条件是存在一个有限时间,使得以下能观性格拉姆矩阵非奇异。,定义,(26),(27),小结,1、关于控制系统的能控能观性:现代控制理论的重要概念,核心内容(最优控制和最优估计)的基础。分别分析输入对状态变量的控制能力和输出对状态变量的反应能力。2、控制系统是否能控能观,如何判断?定义:能控性-在非容许控制作用下,状态变量能在有限时间域内从初始位置到达指定的终端位置;能观性-在有限时间域内,输出能唯一确定系统在初始时刻的状态。判别准则有:一是根据系统模拟结构图,依据定义判别;二是根据A,B阵或是A,C阵,通过矩阵运算来判别。3、已介绍线性连续定常、时变系统的判别方法,继续介绍线性离散系统的判别方法以及其它相关知识。,输入量,状态变量,输出量,能控性,输入量,输出量,G(s),能观性,状态方程,输出方程,3.4 离散系统的能控性和能观测性,线性定常离散系统方程为,(29),3.4.2 能控性判据,(30),定理3-17 系统(29)能控的充分必要条件是能控性矩阵 的秩为n,即,当然也可以直接依据定义来求解,如有系统,利用递推法有,3.4.3 能观测性定义,对于(29)式所描述的系统,根据有限个采样周期的,可以惟一地确定系统的任一初始状态,则称系统是状态完全能观测的。,例3-13 线性定常离散系统方程为,试判断系统的能观测性。,第8次学生经典部分回顾及MATLAB实践讲解题目,1、介绍线性时变系统的运动分析。3、介绍线性连续系统方程的离散化。4、运用SAIMULINK实现离散时间系统和混合系统的集成。,第7次课堂交流给分,1、NYQUIST稳定判据及应用2、控制系统的相对稳定性3、SIMULINK软件包的介绍(1)4、SIMULINK软件包的介绍(2)屈辉:试图讲清楚,但仍然不清楚。给C+田茂:感觉没有责任感,交流不清楚。C;任勇勇:有自信,但交待不太清楚,B张翼:交待清楚,且流畅,不易。对于老师或是同学们或是自己都是负责的。A+曾冰:还可以。B,3.5 对偶原理,线性定常系统方程为,(34),构造一个系统,(35),系统(34)和(35)互为对偶系统。,(上面介绍了系统能控性和能观测性。从概念上和形式上都很相似。它给人们一个启示,即能控性和能观测性之间存在某种内在的联系。这个联系就是系统的对偶原理),(式(35)的系数矩阵为,输入矩阵为,输出矩阵为),对偶系统具有两个基本特征,1.对偶的两个系统传递函数矩阵互为转置,2.对偶的两个系统特征值相同,对偶原理:系统(34)的能控性等价于系统(35)的能观测性;系统(34)的能观测性等价于系统(35)的能控性。,例3-15 线性定常系统如下,判断其能观测性。,该对偶系统的能控性矩阵,对偶系统能控,根据对偶原理,原系统能观测。,有了对偶原理,一个系统的能控性问题可以通过它的对偶系统的能观测性问题的解决而解决;而系统的能观测性问题可以通过它的对偶系统的能控性问题的解决而解决。这在控制理论的研究上有重要意义。,3.6 能控标准形和能观测标准形,(36),3.6.1 能控标准形,线性定常系统,设A的特征多项式,能控性矩阵,推论:具有能控标准形的系统一定能控。,(证明参见参考教材1的109页),例3-16 已知能控的线性定常系统,(2)A 的特征多项式,(3)计算变换矩阵 P,(4)计算,(5)能控标准形,例3-16+试将下列状态空间表达式变换成能控标准I型,3.6.2 能观测标准形,系统(36)的能观测性矩阵为,则系统能观测,推论:具有能观标准形的系统一定能观。,变换矩阵可取为,(39),例:试将例题列状态空间表达式变换为能观标准I型,3.7 能控性、能观性与传递函数的关系,考察SISO线性定常系统,(40),其传递函数为,(41),传递函数的分子、分母分别为,可以看出,在没有零极点对消的情况下,传递函数的特征根和系统矩阵A 的特征值相同。,定理3-24 SISO系统(40)能控又能观的充分必要条件是 不存在零、极点对消。,例3-17 线性定常系统方程如下,求系统传递函数,并且判断系统能控性与能观性。,解 传递函数为,能控性,能观性,例3-19 MIMO线性定常系统方程为,传递函数矩阵,能控性,能观性,可见,传递函数矩阵虽然有零极点对消,但是系统既能控又能观。这是因为极点(s-1)还剩一个,并未消失,只是降低系统重极点的重数。,(42),MIMO线性定常系统,定理3-25 若系统(42)的状态向量和输入向量之间的传递函数矩阵 的各行线性无关,则系统能控。,定理3-26 若系统(42)的输出向量和状态向量之间的传递函数矩阵 的各列线性无关,则系统能观。,3.8 系统的结构分解,一个不能控、不能观测的系统,从结构上来说,必定包括能控、不能控以及能观测、不能观测的子系统。如何按照能控性或能观测性进行分解呢?我们知道,线性变换不改变系统的能控性和能观测性。因此,可采用线性变换方法将其分解。这里必须解决3个问题:1、如何分解?2、分解后系统方程的形式为何?3、变换矩阵如何确定?下面介绍结构分解问题。,线性定常系统,(43),3.8.1 按能控性分解,并且 维子系统为,(45),变换矩阵 的确定方法:因为即矩阵 中有n1个线性无关的列向量,再补充 个列向量,从而构成非奇异的矩阵,例3-20 系统方程如下,要求按能控性进行结构分解。,经过线性变换后,3.8.2 按能观性分解,定理3-28 若系统(43)不能观,且状态 有 个状态分量能观,则存在线性变换,使其变换成下面形式,(47),能观性矩阵 中有 个线性无关的行向量,在它们的基础上,再补充 个行向量,构成变换矩阵。,例3-21 系统方程如下,要求按能观性进行结构分解。,解,从 中任选两个行向量,例如,再补充一个与之线性无关的行向量。,3.8.3 同时按能控性和能观性进行结构分解,定理3-29 若系统(43)不能控,不能观,且存在线性变换,使其变换成下面形式,系统传递函数矩阵,(50),(51),3.9 实现问题,(52),如果给定一个传递函数,求得一个系统方程,注:当传递函数分子的阶次小于分母的阶次时,有(52)式形式;当传递函数分子的阶次等于分母的阶次时,有(53)式形式。,在基于状态空间方法分析和设计控制系统时,要知道系统的状态空间表达式。然而在有的情况下,只知道系统的传递函数(矩阵),这时就要将给定的传递函数(矩阵)描述变成与之输入输出特性等价的状态空间表达式描述。这个问题称为系统实现问题。这里只讨论SISO系统的实现问题。举例则是一个MIMO系统。,3.9.1 能控标准形实现,系统传递函数为,即:,进行拉普拉斯反变换,选择系统的状态变量,于是有,,写成矩阵形式,其中,2.含零点,,,引入中间变量x1=u/D(s),详细推导见刘豹第二版P22,写成矩阵形式,其中,3.9.2 能观标准形实现,系统传递函数为,如果令,于是,写成矩阵形式,例3-22:试求W(s)的能控和能观标准型实现(见P129),最小实现,传递函数阵只能反映系统中能控且能观的子系统动力学行为。将系统中不能控不能观的状态分量消去,不影响系统的传递函数阵。寻求最小实现的步骤,第9次学生经典部分回顾及MATLAB实践讲解题目,1、控制系统的校正及综合设计(概念),2、定常线性系统的相位超前校正3、SIMULINK的分析工具4、SIMULINK的运行机理(S函数模块),3.9.3 并联形实现,为简单起见,以两阶系统传递函数为例,进行介绍。,1)传递函数极点互异,选取,有,则,2)传递函数有重极点,矩阵形式,3.9.4 串联形实现,设,3.10 MATLAB的应用,3.10.1 判断线性系统的能控性和能观测性,用MATLAB可以很方便地求出线性控制系统的能控性矩阵和能观测性矩阵,并且求出它们的秩。从而判断系统的能控性和能观测性。函数ctrb()和obsv()分别计算系统的能控性矩阵和能观测性矩阵。格式为:Qc=ctrb(A,B),Qo=obsv(A,C)。,例 3-23 判断下面的线性系统是否能控?是否能观测?,其中,解 先分别计算系统的能控性矩阵和能观测性矩阵。然后,再用rank()函数计算这两个矩阵的秩。,输入以下语句,这些语句的执行结果为,从计算结果可以看出,系统能控性矩阵和能观测性矩阵的秩都是3,为满秩,因此该系统是能控的,也是能观测的。,注:当系统的模型用sys=ss(A,B,C,D)输入以后,也就是当系统模型用状态空间的形式表示时,我们也可以用Qc=ctrb(sys),Qo=obsv(sys)的形式求出该系统的能控性矩阵和能观测性矩阵。,3.10.2 线性系统按能控性或者能观测性分解,在用MATLAB进行结构分解时,不能控(不能观)的系统,其结构分解的系统方程形式与本章3.8节不同。,当系统能控性矩阵的秩 时,我们可以使用函数命令ctrbf()可以对线性系统进行能控性分解。其调用格式为。其中,T为相似变换矩阵。,输出为一个向量,sum(K)可以求出能控的状态分量的个数。,输出为一个向量,sum(K)可以求出能观测的状态分量的个数。,解 输入下列语句,语句执行结果为,从输出的向量可以看出有两个状态分量是能控的。可以验证,输入语句,得到的结果为,可见,A1=Abar,所得到的结果是正确的。,3.10.3 线性系统转换成能控标准形和能观标准形,下面通过两个例子来说明将系统变换成能控标准形和能观标准形的方法。,运行结果为,表明系统为能控,因此可以变换成能控标准形。而且求出A 的特征多项式为,(即:,),2)计算变换矩阵,输入以下语句,计算结果为,计算结果为,表明经过变换以后的系统方程为,例3-26 系统方程为,其中,求线性变换,将其变换成能观标准形。,解 1)判断系统是否为能观测,并且求出A阵的特征多项式,输入下面语句,运行结果为,表明系统为能观测,因此可以变换成能观标准形。而且求出的特征多项式为,(即:,),2)计算变换矩阵,输入以下语句,计算结果为,3)计算出能观标准形,输入以下语句,计算结果为,表明经过变换以后的系统方程为,