《状态空间模型的线性变换和约旦规范形.ppt》由会员分享，可在线阅读，更多相关《状态空间模型的线性变换和约旦规范形.ppt（89页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、Ch.2 Ch.2 控制系统的状态空控制系统的状态空间模型间模型目录目录(1/1)目目 录录q概述概述q2.1 状态和状态空间模型状态和状态空间模型q2.2 根据系统机理建立状态空间模型根据系统机理建立状态空间模型q2.3 根据系统的输入输出关系建立状态空间模型根据系统的输入输出关系建立状态空间模型 q2.4 状态空间模型的线性变换和约旦规范型状态空间模型的线性变换和约旦规范型q2.5 传递函数阵传递函数阵q2.6 线性离散系统的状态空间描述线性离散系统的状态空间描述 q2.7 Matlab问题问题 q本章小结本章小结状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(1(1
2、/8)/8)2.4 状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形q从上一节的讨论可知,同一个系统的状态空间模型,即使其维数相同,但其具体结构和系数矩阵也是多种多样的,如系统矩阵A可以为对角线矩阵的或者约旦矩阵的,也可以为其他形式的。即,状态空间模型不具有唯一性状态空间模型不具有唯一性。状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(2/8)2/8)为何同一个系统具有不同的状态空间模型?原因原因:状态变量的不同选择这就产生了一个问题:各种不同选择的状态变量之间,以及它们所对应的状态空间模型之间的关系如何?状态空间模型的线性变换和约旦规范形状态空间模型
3、的线性变换和约旦规范形(3/8)3/8)q此外,在控制系统的分析和设计中,某些特殊的系统数学模型对讨论问题相对简单得多,如前面建立的对角线规范形的和约旦规范形。于是自然会提出如下问题:如何把一般形式的状态空间模型变换成特定形式的状态空间模型,以降低系统的分析问题和设计问题的难度。解决上述两个问题,就需引入状态空间的线性变换。什么是状态空间的线性变换什么是状态空间的线性变换?如何理解?本章关键喔!状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(4/8)4/8)q状态变量是一组实变量,它们所组成的状态空间为一个实线性空间。由线性代数知识可知,线性空间中,随着表征空间坐标的基底
4、的基底的选取的不同选取的不同,空间中的点关于各种基底的坐标亦不同坐标亦不同。这些基底之间的关系为进行了一次坐标变换,而空间中的点的 坐标则相当于作了一次相似变换。如,在如右图所示的平面直角坐标系中,A点在两个坐标系下的坐标存在如下变化关系(其中P为非可逆的变换矩阵)状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(5/8)5/8)n维空间中的旋转变换、极坐标变换,线性空间中的相似变换,都属于空间变换。其中旋转变换和相似变换还属于线性变换。状态空间中由于状态变量的不同选择类似于线性空间中的坐标架的不同选择，同一个系统不同选择状态变量组之间存在类似于线性空间不同坐标架之间的线性
5、变换,因此我们将在状态空间中坐标变换称为状态空间的线性变换。状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(6/8)6/8)q引入坐标变换和状态空间线性变换等概念,实际上就回答了上述两个问题:1.不同选取状态变量之间存在一个坐标变换坐标变换,其相应的状态空间模型之间也存在一个相应的相似变换相似变换。2.既然可以对状态变量和状态空间模型进行线性变换,则在一定条件下应可以将一般形式的状态空间模型变换成某种特殊的状态空间模型。状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(7/8)7/8)q本节主要讨论状态空间的线性变换,以及如何系统状态空间描述的约旦规范
6、形。本章关键问题:1.线性变换的几何及空间意义,建立空间想象力2.如何作系统线性变换3.系统的对角规范形和约旦规范形描述4.代数重数、几何重数与约旦矩阵5.如何求矩阵的广义特征向量建立空间概念,可是学好控制理论的关键喔状态空间模型的线性变换和约旦规范形状态空间模型的线性变换和约旦规范形(8/8)8/8)主要内容为:状态空间的线性变换状态空间的线性变换系统特征值的不变性与系统的不变量系统特征值的不变性与系统的不变量化状态方程为对角线规范形化状态方程为对角线规范形化状态方程为约旦规范形化状态方程为约旦规范形状态空间的线性变换状态空间的线性变换(1/(1/2)2)2.4.1 状态空间的线性变换状态空
7、间的线性变换q对于一个n阶动态系统,可通过选择适当的n个状态变量以建立状态空间模型来描述它。但是,这n个状态变量的选择却不是唯一的。这一点可利用线性代数中的基底不唯一来理解。一个n维线性独立的状态变量向量,在n维状态空间中构成一个坐标系,即相当于空间中的一个基底。根据线性代数知识,在这个空间中还存在另外的坐标系,且与原坐标系存在一个线性变换关系。状态空间的线性变换状态空间的线性变换(2/2)2/2)下面分别讨论:状态空间的线性变换状态空间的线性变换状态空间模型的线性变换状态空间模型的线性变换上述状态变量向量x与 间的变换,称为状态的线性变换。由线性代数知识可知,它们之间必有如下变换关系状态空间
8、的线性变换状态空间的线性变换(1/1)(1/1)1.状态空间的线性变换q设描述同一个线性线性状态空间的两个n维的状态变量向量分别为其中P为nn维的非奇异变换矩阵。值得指出的是:变换矩阵P只有为非奇异的,才能使x和 间的变换关系是等价的、唯一的和可逆的。q两种表达式式之间存在什么关系两种表达式式之间存在什么关系?状态空间的线性变换状态空间的线性变换(1(1/14)/14)2.状态空间模型的线性变换q设在状态变量x和 下,系统状态空间模型分别为将变换关系x=P 代入(A,B,C,D)的状态方程中有状态空间的线性变换状态空间的线性变换(2/14)2/14)由于变换矩阵P非奇异,因此有则有应该注意的是
9、,系统的初始条件也必须作相应的变换,即将上式与状态空间模型 比较,则线性系统(A,B,C,D)在线性变换矩阵P下的各矩阵具有如下对应关系其中t0为系统运动的初始时刻。q在进行状态空间的线性变换中,需要计算矩阵的逆,下面先简要复习一下逆矩阵的计算。补充资料补充资料:逆矩阵的计算逆矩阵的计算。q常用的逆矩阵计算方法有如下2种:计算伴随矩阵法计算伴随矩阵法。初等变换法初等变换法(三角矩阵变换法三角矩阵变换法)。下面分别介绍。状态空间的线性变换状态空间的线性变换(3/14)逆矩阵的计算逆矩阵的计算1.计算伴随矩阵法计算伴随矩阵法该方法计算式如下:P-1=adj(P)/|P|其中adj(P)和|P|分别
10、为矩阵P的伴随矩阵和行列式。伴随矩阵伴随矩阵的定义与计算如下:设有矩阵P为状态空间的线性变换状态空间的线性变换(4/14)计算伴随矩阵法计算伴随矩阵法状态空间的线性变换状态空间的线性变换(5/14)计算伴随矩阵法计算伴随矩阵法则其伴随矩阵为:其中pij*则为矩阵P的元素pij的代数余子式代数余子式。代数余子式代数余子式pij*为nn矩阵P去掉第i行第j列余下的n-1行n-1列的行列式值乘以符号。状态空间的线性变换状态空间的线性变换(6/14)计算伴随矩阵法计算伴随矩阵法例例 计算下述矩阵的逆矩阵。计算过程如下:(1)先计算代数余子式先计算代数余子式。状态空间的线性变换状态空间的线性变换(7/1
11、4)计算伴随矩阵法计算伴随矩阵法(2)计算伴随矩阵计算伴随矩阵。(3)计算行列式值计算行列式值。状态空间的线性变换状态空间的线性变换(8/14)计算伴随矩阵法计算伴随矩阵法(4)计算逆矩阵计算逆矩阵。状态空间的线性变换状态空间的线性变换(9/14)初等变换法初等变换法2.初等变换法初等变换法初等变换法(即三角矩阵变换法)求逆矩阵P-1的思想为:将矩阵P与单位矩阵I组成增广矩阵P I。通过对增广矩阵P I作初等行变换初等行变换,将P I变换成I T,则P-1=T。初等行变换初等行变换有如下三种:v将某两行交换v将某行乘以一非零常数v将某行乘以一常数加到另一行对前述的矩阵P求逆,计算过程例示如下:
12、状态空间的线性变换状态空间的线性变换(10/14)初等变换法初等变换法状态空间的线性变换状态空间的线性变换(11/14)-初等变换法初等变换法下面简单总结2种计算逆矩阵方法的计算时间复杂性(计算时间估计)。由逆矩阵与伴随矩阵的计算式可知,基于伴随矩阵计算的逆矩阵的计算时间复杂性(计算时间估计)为v1次nn的行列式值vn2次(n-1)(n-1)的行列式值合计为O(n5)而初等变换法的计算时间复杂性为O(n3)状态空间的线性变换状态空间的线性变换(12/14)例例2-5q基于上述逆矩阵计算,下面介绍状态空间模型变换的算例。q例例2-5 试将以下状态空间模型作变换矩阵为下式所示的线性变换状态空间的线
13、性变换状态空间的线性变换(13/14)q解 线性变换P的逆矩阵为因此,有状态空间的线性变换状态空间的线性变换(14/14)故系统在新的状态变量下的状态空间模型为q值得指出的是,状态空间的线性变换只是对状态变量作变换,对系统的输入和输出未作变换,因此系统的输入输出间的动态和静态关系对状态变换保持不系统的输入输出间的动态和静态关系对状态变换保持不变变。系统特征值的不变性与系统的不变量系统特征值的不变性与系统的不变量(1/2)2.4.2 系统特征值的不变性与系统的不变量系统特征值的不变性与系统的不变量q由前面的讨论可知,当选择不同的状态变量,则获得不同的状态空间模型描述。实际上,状态空间模型只是系统
14、在不同的状态变量选择下对系统的一种描述,它随状态变量选择的不同而不同,并不具有唯一性和不变性。那么,到底系统在状态空间中有哪些描述,哪些性质是不变的,是不随状态变量的选取不同而变化的?q线性定常系统的特征结构由特征值和特征向量所表征。系统的特征结构对系统运动的特性和行为具有重要的影响,决定了系统的基本特性。系统特征值的不变性与系统的不变量系统特征值的不变性与系统的不变量(2/2)下面我们将讨论系统经状态线性变换后,其特征值不变,亦即状态线性变换不改变系统的基本特性。系统矩阵的特征值是一种描述系统本质特征的,并具有唯一性的不变量,即不随状态变量的选取不同而变化的不变量,它在系统分析和综合上起着重
15、要的作用。下面将分别讨论:系统的特征值和特征向量系统的特征值和特征向量系统系统特征值的不变性特征值的不变性特征向量的计算特征向量的计算广义特征向量和特征向量链广义特征向量和特征向量链难点喔!重点喔系统的特征值和特征向量系统的特征值和特征向量(1/4)1.系统的特征值和特征向量q状态空间的线性变换,只是改变了描述系统的角度(或说坐标系),系统的本质特征应保持不变。对于线性定常系统来说,系统的特征值(极点)决定了系统的基本特性。特征值应是系统不变的本质特征之一。系统经状态线性变换后,其本质特征之一的特征值应保持不变,亦即状态线性变换不改变系统的基本特性。下面先讨论矩阵特征值和特征向量的定义。系统的
16、特征值和特征向量系统的特征值和特征向量(2/4)特征值和特征向量定义特征值和特征向量定义q定义定义2-2 设v是n维非零向量,A是nn矩阵。若方程组Av=v成立,则称为矩阵A的特征值特征值,非零向量v为所对应的矩阵A的特征向量特征向量。q将上述特征值的定义式写为(I-A)v=0 其中I为nn的单位矩阵。因此,由代数方程论可知,上式有非零特征向量v的解的充要条件为|I-A|=0 并称上式为矩阵A的特征方程特征方程,而|I-A|为A的特征多项式特征多项式。系统的特征值和特征向量系统的特征值和特征向量(3/4)特征值和特征向量定义特征值和特征向量定义q将|I-A|展开，可得|I-A|=n+a1n-1
17、+an-1+an=0其中ai(i=1,2,n)称为特征多项式的系数。因此,nn维的矩阵A的特征多项式为n阶多项式。若矩阵A为实矩阵,则对应的特征方程为一实系数代数方程,共有n个根。这n个根或为实数,或为成对出现的共轭复数。求解矩阵特征值的方法即为求解矩阵A的特征方程。n阶的特征方程的n个根1,2,n即为矩阵A的n个特征值。在得到特征值i后,由式(2-46)或式(2-47)可求得矩阵对应于i的特征向量vi。系统的特征值和特征向量系统的特征值和特征向量(4/4)q如下定义所示,矩阵特征值的概念可推广至线性定常系统(A,B,C,D)。q定义 对于线性定常系统(A,B,C,D),系统的特征值即为系统矩
18、阵A的特征值。q关于系统特征值,几点注记:A.一个n维线性定常系统必然有n个特征值与之对应。B.对于物理上可实现的系统,其系统矩阵必为实矩阵。因此,线性定常系统的特征多项式必为实系数多项式,即系统的特征值或为实数,或为成对出现的共轭复数。系统系统特征值的不变性特征值的不变性(1/2)(1/2)2.系统特征值的不变性q系统的特征值表征了系统本质的特征。而线性变换只是相当于对系统从另外一个角度来描述而已,并未改变系统的本质。刻划了系统本质特征的系统特征值应不随线性变换而改变,即有如下结论:线性定常系统特征值对线性变换具有不变性。系统系统特征值的不变性特征值的不变性(2/2)(2/2)q对于这个结论
19、,亦可证明如下:设系统原状态空间模型中的系统矩阵为A,经线性变换后,系统矩阵为可见,系统经线性变换后,其特征值不变。矩阵 的特征多项式为即证明了A的特征多项式等于的 特征多项式。特征向量的计算(1(1/9)/9)3.特征向量的计算q如何求解特征值i对应的特征向量?求解特征向量,即求如下齐次矩阵代数方程的非零解(iI-A)vi=0q由于i为A的特征值,故iI-A不可逆。因此,由代数方程理论可知,该方程组的解并不唯一。由特征向量的定义可知,我们需求解的是线性独立的特征向量。实际上,具体求特征向量时,可假定其特征向量的某个或几个元素的值,然后再求得该特征向量其他元素的值。特征向量的计算(2(2/9)
20、/9)q当特征方程存在重根时,线性独立的特征向量可能不唯一。因此,就产生如下问题:问题问题:对应于特征值i究竟有几个独立的特征向量?答案答案:矩阵的重特征值i所对应的线性独立的特征向量可能不止一个。v它的独立特征向量的数目等价于系统的维数与线性方程组(2-47)的线性独立的方程数之差,即为n-rank(iI-A)其中rank为矩阵的秩。特征向量的计算(3(3/9)/9)q因此,r重的特征值可能存在1至r个线性独立的特征向量。由此,导出如下问题:独立的特征向量数到底具有什么意义?它与特征值的重数之间有何关系?下面引入代数重数与几何重数两个概念。不要混淆喔!特征向量的计算(4(4/9)/9)q两个
21、基本概念:代数重数代数重数。由特征方程求得的特征值i的重数称为特征值i的代数重数。几何重数几何重数。特征值i线性独立的特征向量数称为特征值i的几何重数。代数重数和几何重数是两个不同的概念。几何重数具有几何上空间表征的意义,它代表在空间分解上不变的几何子空间的数目。而代数重数仅具有代数意义,它代表特征值在特征方程的重数。特征向量的计算(5(5/9)/9)例例2-62-6q例2-6 求如下矩阵的特征向量q解解 1.由特征方程|I-A|=0求得系统的特征值。特征向量的计算(6/9)6/9)例例2-62-6解该特征方程,可求得系统的特征值为1=1 2=3=2即2为系统的二重特征值,其代数重数为22.计
22、算1=1的特征向量。按定义有(1I-A)v1=0即特征向量的计算(7/9)7/9)例例2-62-6解之得特征向量v1的通解为v1=v11 v11 2v11令v11=1,解之得v1=v11 v12 v13=1 1 2特征向量的计算(8/9)8/9)例例2-62-63.计算重特征值2=3=2的特征向量。按定义有(2I-A)v2=0即特征向量的计算(9/9)-9/9)-例例2-62-6由于n-rank(2I-A)=2因此,特征值应有2个独立特征向量,故该重特征值的几何重数亦为2。解之得特征向量v2的通解为v2=v21 v22 v21令v21=1,v22=0和1,解之得v2=1 0 1 和 v3=1
23、1 1即重特征值2有两个线性独立的特征向量。广义特征向量和特征向量链(1(1/12)/12)4.广义特征向量和特征向量链q某些重特征值的线性独立特征向量数(几何重数)小于其代数重数,从而使得矩阵所有特征值所对应的线性独立特征向量数之和小于矩阵维数。为此,为能进行空间的结构分解和分析,下面引入一组辅助的空间变换基向量-广义特征向量和特征向量链。q定义 广义特征向量是重特征值i所对应的某个线性独立的特征向量vj满足如下方程组的向量vj,k:广义特征向量和特征向量链(2(2/12)/12)q解上述方程组一直到无解为止,就可求得特征值i的特征向量vj所对应的所有广义特征向量vj,k。q重特征值i的所有
24、线性独立特征向量vj及其对应的广义特征向量vj,k的个数等于其代数重数,否则就还存在其他特征向量或广义特征向量。值得指出的是,并不是重特征值i的任何一组线性独立的特征向量,都能求出所有的广义特征向量。若i的某一组特征向量vj及其相应广义特征向量vj,k的个数小于该特征值的代数重数,则应重新选取其他一组线性独立的特征向量并求取相应的广义特征向量。广义特征向量和特征向量链(3(3/12)/12)q重特征值i的特征向量vj的广义特征向量vj,1,vj,2,组成的向量链称为i的特征向量vj对应的特征向量链。广义特征向量并不是矩阵的特征向量,它只是与对应的特征向量组成该矩阵在n维线性空间中的一个不变子空
25、间。矩阵的所有特征向量和广义特征向量线性独立,并且构成n维线性空间的一组基底。这在矩阵分析中是相当重要的。广义特征向量和特征向量链(4(4/12)/12)q下面通过一个例子来简单介绍线性空间的特征子空间分解。例,某5维线性空间,存在一个3重特征值和一个2重特征值。3重特征值有2个独立特征向量,2重特征值有1个独立特征向量。则该线性空间可分解为如下3个独立的不变特征子空间。广义特征向量和特征向量链(5(5/12)/12)广义特征向量和特征向量链(6(6/12)/12)若该5维线性空间,3重特征值有1个独立特征向量,2重特征值有2个独立特征向量。则该线性空间可分解为如下3个独立的不变特征子空间。广
26、义特征向量和特征向量链(7(7/12)/12)例例2-72-7q例例2-7 求如下矩阵的特征向量和特征向量链q解解 1.由特征方程|I-A|=0可求得系统的特征值为1=2=3=-1即-1为系统的三重特征值,其代数重数为3。2.计算对应于三重特征值-1的特征向量。按定义有(1I-A)v1=0广义特征向量和特征向量链(8(8/12)/12)例例2-72-7即由于n-rank(1I-A)=2因此,该特征值应有2个独立特征向量,故该重特征值的几何重数亦为2。由于该重特征值的几何重数小于代数重数，因此存在广义特征向量。解之得如下特征向量的通解式 v1=v11 v12 -(v11+v12)/2广义特征向量
27、和特征向量链(9(9/12)/12)例例2-72-7分别令两组独立的v11 v12即可求得三重特征值1的两个线性独立的特征向量。三重特征值-1只有两个线性独立特征向量,其几何重数为2。因此,重特征值-1的两个独立特征向量中有一个一定存在广义特征向量。下面通过求广义特征向量来辅助决定选取合适的v11和v12。广义特征向量和特征向量链(10/12)10/12)例例2-72-73.计算对应于特征向量的广义特征向量和特征向量链。按定义式(2-51),特征向量v1的广义特征向量v1,2满足(1I-A)v1,2=-v1即因此,根据方程的可解性,存在广义特征向量的特征向量v1中的v11和v12满足v11=-
28、3v123倍关系广义特征向量和特征向量链(1(11/12)1/12)例例2-72-7此时的广义特征向量的解为v1,2=r1 r2 -(r1+r2-v12)/2其中r1和r2为任意数。因此存在广义特征向量的特征向量v1为和其对应的广义特征向量可以分别取为v1=v11 v12 -(v11+v12)/2 =-3v12 v12 v12 =1 -1/3 -1/3v1,2=r1 r2 -(r1+r2-v12)/2 =1 2/3 -1广义特征向量和特征向量链(1(12/12)2/12)例例2-72-7另外一个不存在广义特征向量的三重特征值1的特征向量为v2=v11 v12 -(v11+v12)/2=1 0
29、-1/2q本例共求得3个特征向量和广义特征向量。由于矩阵A的维数为33,因此对应于上述特征向量和广义特征向量,已不存在其他广义特征向量。故特征值1对应于特征向量v1的特征向量链为v1和v1,2。化状态方程为对角线规范形化状态方程为对角线规范形(1(1/12)/12)2.4.3 化状态方程为对角线规范形化状态方程为对角线规范形q对角线规范形是指系统矩阵A为对角线矩阵的一类状态空间模型。对于该类状态空间模型,由于在系统分析和综合时,清晰直观,使问题得以简化该类系统可简化成n个一阶惯性环节的并联故在状态空间分析法中是较重要的一类特殊状态空间模型。任何具有任何具有n个线性独立特征向量的状态空间模型一定
30、能个线性独立特征向量的状态空间模型一定能经状态变换变换成对角线规范形。经状态变换变换成对角线规范形。该结论可详细地并构造性地证明如下。化状态方程为对角线规范形化状态方程为对角线规范形(2(2/12)/12)q结论 已知线性定常系统的状态方程为其中系统矩阵若A的n个特征值1,2,n所对应的特征向量线性独立,则必存在变换矩阵P,使其进行状态变换x=P 后为对角线规范形,即系统的状态方程为为对角线矩阵,并且变换矩阵P可取为P=p1 p2 pn其中pi为矩阵A对应于特征值i的特征向量。三、化状态方程为对角线规范形三、化状态方程为对角线规范形(3(3/12)/12)q证明 若pi为对应与特征值i的独立特
31、征向量,则必有Api=ipi因此有Ap1 Ap2 Apn=1p1 2p2 npn对上式两边分别有Ap1 Ap2 Apn=Ap1 p2 pn=AP化状态方程为对角线规范形化状态方程为对角线规范形(4(4/12)/12)故AP=Pdiag1 2 n即P-1AP=diag1 2 n即证明了结论。对原状态方程进行线性变换 的后,可得化状态方程为对角线规范形化状态方程为对角线规范形(5/12)-5/12)-例例2-82-8q例2-8 试将下列状态空间模型变换为对角线规范形化状态方程为对角线规范形化状态方程为对角线规范形(6/12)-6/12)-例例2-82-8q解解 1.先求A的特征值。由特征方程可求得
32、特征值为1=-1 2=-2 3=-32.求特征值所对应的特征向量。由前述的方法可求得特征值1,2和3所对应的特征向量分别为p1=1 0 1 p2=1 2 4 p3=1 6 93.取A的特征向量组成变换矩阵P并求逆阵P-1,即有化状态方程为对角线规范形化状态方程为对角线规范形(7/12)7/12)例例2-82-84.计算各矩阵5.系统在新的状态变量下的状态空间模型为化状态方程为对角线规范形化状态方程为对角线规范形(8/12)8/12)q下面给出快速计算矩阵特征向量及对角线规范形的一个特例:在第三节讨论的状态空间模型中,其系统矩阵为其特征多项式为|I-A|=n+a1n-1+an-1+an即该类矩阵
33、的最后一行与特征多项式的系数一一对应。该类特殊系统矩阵A称为友矩阵友矩阵。单位矩阵化状态方程为对角线规范形化状态方程为对角线规范形(9/12)9/12)友矩阵的特征向量的特点:当特征值为i时,其对应的特征向量为该结论可由下式证明。即pi为友矩阵的特征值i对应的特征向量。化状态方程为对角线规范形化状态方程为对角线规范形(10/12)-10/12)-例例2-92-9因此,当友矩阵的特征值互异时,将友矩阵变换成对角线矩阵的变换矩阵恰为下述范德蒙矩阵范德蒙矩阵q例2-9 试将下列状态空间模型变换为对角线规范形三、化状态方程为对角线规范形三、化状态方程为对角线规范形(1(11/12)-1/12)-例例2
34、-92-9q解解 1.先求A的特征值。由特征方程可求得特征值为1=0 2=-1 3=-22.由于A为友矩阵,故将A变换成对角线矩阵的变换矩阵P及其逆阵P-1分别为三、化状态方程为对角线规范形三、化状态方程为对角线规范形(1(12/12)2/12)例例2-92-93.计算各矩阵4.系统在新的状态变量下的状态空间模型为化状态方程为约旦规范形化状态方程为约旦规范形(1/1)(1/1)2.4.4 化状态方程为约旦规范形化状态方程为约旦规范形q若系统存在重特征值且线性独立特征向量数小于该特征值的重数时,则系统矩阵A不能变换成对角线矩阵。在此种情况下,A可变换成约旦矩阵,系统表达式可变换成约旦规范形。下面
35、将分别讨论约旦块和约旦矩阵约旦块和约旦矩阵约旦规范形及其计算约旦规范形及其计算约旦块和约旦矩阵(1/(1/3)3)1.约旦块和约旦矩阵q矩阵的约旦块的定义为q由l个约旦块Ji组成的块对角的矩阵称为约旦矩阵,如J=block-diagJ1 J2 Jl约旦块和约旦矩阵(2/(2/3)3)q下述矩阵均为约旦矩阵上述第一个约旦矩阵有两个约旦块,分别为11维的特征值2的约旦块和33维的特征值-1的约旦块;第二个约旦矩阵有三个约旦块,分别为11维的特征值3的约旦块以及11维和22维的特征值-1的两个约旦块。约旦块和约旦矩阵(3/3)3/3)q由约旦块和约旦矩阵的定义可知,对角线矩阵可视为约旦矩阵的特例对角
36、线矩阵可视为约旦矩阵的特例,其每个约旦块的维数为11。在本课程中,若未加以特别指出的话若未加以特别指出的话,则所有对约旦矩阵则所有对约旦矩阵有关的结论都同样适用于对角线矩阵。有关的结论都同样适用于对角线矩阵。约旦规范形及其计算(1/16)(1/16)2.约旦规范形及其计算q定义 系统矩阵A为约旦矩阵的状态空间模型称为约旦规范形。q与对角线规范形一样,约旦规范形也是线性定常系统的状态空间分析中一种重要的状态空间模型。下面讨论一般状态空间模型与约旦规范形之间的线性变换的计算问题。q对于任何有重特征值且其线性独立特征向量数小于其维数的矩阵,虽然不能通过相似变换化成对角线矩阵,但可经相似变换化为约旦矩
37、阵。可经相似变换化为约旦矩阵。约旦规范形及其计算(2/16)(2/16)q状态空间模型变换与对角线规范形、约旦矩阵规范形的关状态空间模型变换与对角线规范形、约旦矩阵规范形的关系系?一般状态空间表达式对角线规范形约旦规范形n个独立特征向量代数重数=几何重数代数重数几何重数n个独立特征向量与广义特征向量特例线性变换Understand?约旦规范形及其计算(3/16)(3/16)q若将对角线矩阵视为约旦矩阵的特例的话,则任何矩阵皆可经相似变换化为约旦矩阵。相应地,任何状态空间模型都可经状态变换变换成约旦规范形。任何矩阵都可变换成约旦矩阵,但能变换成有几个约旦块的约旦矩阵,则与系统的特征向量有关。对此
38、有如下结论:矩阵所变换成的约旦矩阵的约旦块数等于该矩阵的线性独立特征向量数(即几何重数)。约旦规范形及其计算(4/16)(4/16)q由前面讨论可知:任何状态空间模型一定能经状态变换变换成约旦规范形。该结论可详细地并构造性地叙述并证明如下。约旦规范形及其计算(4/16)(4/16)q结论 已知线性定常系统的状态方程为x=Ax+Bu若A的共有p(pn)个互异的特征值,l(pln)个线性独立特征向量pi及相应地广义特征向量pi,j(i=1,2,l;j=1,2,mi),则必存在变换矩阵P,使其进行状态变换x=P 后为约旦规范形,即系统的状态方程为其中系统矩阵为约旦矩阵,并且变换矩阵P可取为P=P1
39、P2 Pl约旦规范形及其计算(5/16)(5/16)变换矩阵P P=P1 P2 Pl中的Pi为矩阵A对应于线性独立特征向量pi的特征向量链组成的如下分块矩阵q证明 若pi和pi,j为对应与特征值i的独立特征向量和广义特征向量,则必有约旦规范形及其计算(6/16)(6/16)q因此有其中Ji为相应的约旦块。Api=ipi约旦规范形及其计算(7/16)(7/16)即P-1AP=block-diagJ1 J1 Jl故APi=PiJi约旦规范形及其计算(8/16)(8/16)例例2-102-10即对原状态方程进行线性变换 的后,可得=P-1AP=block-diagJ1 J2 Jl即证明了结论。q例2
40、-10 试将下列状态空间模型变换为约旦规范形约旦规范形及其计算(9/16)(9/16)例例2-102-10q解解 1.先求A的特征值。由特征方程可求得特征值为1=2=3=2 4=-12.求特征值所对应的特征向量。由前述的方法可求得特征值2由如下两个线性独立特征向量P1,1=1 1 -1 1/3 P2,1=1 0 0 -1其中p1,1无广义特征向量,而p2,1的广义特征向量为P2,2=1 1 0 -1特征值-1的特征向量为P3,1=0 0 0 1约旦规范形及其计算(10/16)(10/16)例例2-102-103.取A的特征向量和广义特征向量组成变换矩阵P并求逆阵P-1,即有约旦规范形及其计算(
41、11/16)(11/16)例例2-102-104.计算各矩阵约旦规范形及其计算(12/16)(12/16)例例2-102-105.系统在新的状态变量下的状态空间模型为约旦规范形及其计算(13/16)(13/16)q对前面讨论的特殊矩阵-友矩阵,它的广义特征向量的快速计算方法为:当特征值为i时,其对应的特征向量和广义特征向量分别为约旦规范形及其计算(14/16)(14/16)例例2-112-11q解解 1.先求A的特征值。由特征方程可求得特征值为1=-1 2=3=-2其中mi为该特征值的代数重数。该结论可由广义特征向量和友矩阵的定义证明。q例2-11 试将下列状态空间模型变换为约旦规范形约旦规范形及其计算(15/16)(15/16)例例2-112-113.计算各矩阵2.由于A为友矩阵,故将A变换成对角线矩阵的变换矩阵P及其逆阵P-1分别为约旦规范形及其计算(16/16)-(16/16)-例例2-112-114.系统在新的状态变量下的状态空间模型为