现代控制理论-3控制系统的状态方程求解.ppt

《现代控制理论-3控制系统的状态方程求解.ppt》由会员分享，可在线阅读，更多相关《现代控制理论-3控制系统的状态方程求解.ppt（134页珍藏版）》请在三一办公上搜索。

1、1,现代控制理论基础,2,第二章 控制系统的状态方程求解,2.1 线性定常系统状态方程的解2.2 线性定常系统状态转移矩阵的几种求法2.3 线性离散系统的状态空间表达式及连续系统 的离散化2.4 线性定常离散系统状态方程的求解,3,2.1.线性定常连续系统状态方程的解,可见，输出方程求解要依赖状态方程的解。关键是求解状态方程。本节重点来讨论这个问题。先讨论自由运动的规律，即求自由解。,前面我们详细讨论了状态空间表达式的建立及相互转换。在建立了新的数学模型之后，接着就是求解问题。,由于状态空间表达式由两部分组成,即,4,一、齐次状态方程的解所谓齐次状态方程，与齐次微分方程类似，即输入u(t)=0

2、的情况。故齐次方程为：设初始时刻 t0=0，初始状态为x0 1.采用拉氏变换法求解：对齐次方程两边取拉氏变换.,反变换即得齐次状态方程的解：,5,下面就来讨论：,-解的变化是按指数形式变化的。对于状态方程的解，是否也具有指数形式呢？,2.级数展开法：分析标量微分方程可知,6,7,逐项变换,即 x(t)=e-Atx0,当初始时刻为t00,初始状态为x(t0)时,所以齐次状态方程的解可写为,8,3.求齐次状态解的关键是求 转移矩阵 eAt，前面已给出了两种方法：,2.系统状态的变化实质上是从初始状态开始的状态转移，而转移规律取决于 eAt，eA(t-t0)故称其为状态转移矩阵.一般用,1.齐次状态

3、方程的解表示了系统在初始条件作用下的自由运动，又称为零输入解；,小结：,来表示。,9,a)拉氏变换法：,例：已知系统的状态方程为：试求在初始状态 时的状态解。,由于按幂级数计算不易写出闭式的结果，故通常用拉氏变换法。,b)幂级数法：,10,解：1.求eAt,11,所以,2.求x(t):,12,二.状态转移矩阵：,的解(t),定义为系统的状态转移矩阵。,1.定义：线性定常系统，初始时刻t0=0,满足以下矩阵微分方程和初始条件,在状态空间分析中状态转移矩阵是一个十分重要的概念。,13,讨论：（1）满足上述定义的解为(t)=eAt(t0=0),证明：,14,所以当(t)=eAt时，又因为(t)=eA

4、t(t=0时)eA0=I+A0+.=I 所以(0)=I故 eAt 是状态转移矩阵(t),（2）状态转移矩阵(t)是A阵同阶的方阵，其元素均为时间函数.,15,由于状态转移矩阵具有矩阵指数函数的形式，故可推出如下性质,2.性质：,（1）(t-t0)是非奇异阵.且,16,（2）,其中,17,（3）,（4）,18,由此关系 可用于从 eAt 反求 A.,例：已知,（5）,19,（6）若,则,20,21,块对角阵、约旦块矩阵见P87 2)、3),22,当系统输入u0 时，其S-E为.,直接用分离变量法积分求解方程与采用拉式变换法求解方程,其结果是一致的.第一种方法：直接求解法,三.非齐次状态方程的解：

5、,23,左乘 e-At：,移项：,即,在区间t0,t上积分,24,结论：非齐次状态方程的解由两部分组成：a).由初始状态产生的自由分量零输入解b).由输入引起的强迫分量零状态解,即,或：,25,例：已知系统,由前例得：,解：1.求 eAt：,试求：x(0)=0,u(t)=1(t)时的状态解。,26,2.求x(t),27,将该非齐次状态方程两边取拉氏变换,可得sX(s)-x0=AX(s)+BU(s)即 X(s)=(sI-A)-1x0+BU(s)其中X(s)和U(s)分别为x(t)和u(t)的拉氏变换。对上式两边取拉氏反变换,并利用卷积分公式,则有,上述求解的关键为等式右边第二项。,第二种方法：拉

6、氏变换法,28,下面先回顾卷积积分的拉氏变换法则。设W1(s)和W2(s)分别为原函数f1(t)和f2(t)的拉氏变换,则f1(t)和f2(t)的卷积的拉氏变换为,结果与直接求解法完全相同。,对上述状态方程的求解式利用卷积分公式,则有,29,所谓脉冲响应，即初始条件为零时，输入u为单位脉冲函数(t)，系统的输出称为脉冲响应。,四.系统的脉冲响应及脉冲相应矩阵：,根据这个定义，可求线性定常系统的脉冲响应。但是多变量系统的输入有r个，输出有m个。则脉冲响应显然与传递函数阵的维数不同,即系统地输出为Y(s)=G(s)U(s)是 m1维的列向量.而G(s)是mr维矩阵.,在单变量系统定义脉冲响应函数为

7、 h(t)=L-1G(s),30,即 h(t)=L-1G(s)mr,而 y(t)=L-1G(s)U(s)m1,为了将这一含义推广到多变量系统，我们按以下方式定义脉冲响应函数阵。以后将会知道，在多变量系统中,脉冲响应函数阵虽不等于真正的脉冲响应输出 y(t),但却等于传递矩阵的拉式反变换。,定义：mr 阶矩阵 h(t)=CeAtB 称为系统的脉冲响应矩阵。,31,状态解为：,初始时刻t0=0初始状态x(0)=0,设系统的状态空间表达式为,则输出解为：,32,讨论单变量系统的情况:当输入,-卷积,33,以上关系表明h(t)包含了G(s)的全部信息，也反映系统的基本传递特性。,反之,性质：1.h(t

8、)是传递矩阵的拉式变换,34,2.h(t)在线性变换下的不变性：,即,证明:令 线性变换后.,其中:,35,则状态转移矩阵满足以下性质：,一般有：,36,1.齐次状态方程的解：,小结：本节主要讨论了状态求解的问题：,2.非齐次状态方程的解：,37,4.脉冲响应矩阵：,定义：满足矩阵微分方程 的解（t）,3.状态转移矩阵：,38,2.2 线性定常连续系统(t)的算法,1.对低阶系统（三阶以下）计算较方便，写出的结果是解析式，在实际中最常用。,特点：,一.拉氏变换法：前面已在求状态解时推出,在线性定常系统状态方程的求解中，关键是求(t),本节介绍几种算法:,2.对于高阶系统，会遇到求逆的困难,如,

9、39,求逆阵可采用一些数值计算方法，用计算机计算。求逆变换关键是高阶分解因式，部分分式展开很麻烦。,二.幂级数法：此法是一种直接计算法，前面已介绍过。,40,特点：是一种无穷开式，很难写成闭式，一般采用近似计算，精度将取决于所取项数的多少，适合于计算机计算。,例：,已知系统的状态方程为：,试求其状态转移矩阵.,解：将A阵代入幂级数展开式,41,42,三.对角形法与约当标准形法：,1.矩阵A的特征值 12n 互不相同，其状态转移矩阵可由下式求得,其中：P是使A化成对角形的线性变换。,43,则,证明：12n 互异，必有非奇异矩阵P，将A化成对角形，即：,44,小结：利用对角线法 eAt的方法:1.

10、求 12n（条件：12n 互异）;2.求特征矢量：P1P2Pn;3.写出变换阵 P=P1P2Pn,求出P-1 4.求 eAt：特点：求P阵比较麻烦，常用于理论推导。,45,例：已知用对角形求（t）解：1.求特征值：,46,2.求特征矢量：,即,解出：,47,48,49,3.求P，P-1：4.求 eAt：,50,51,2.矩阵A有相重特征值:定理：若矩阵A有相重特征值，其状态转移矩阵可由下式求得,52,eAt=QeJtQ-1 其中：Q是使A化为约当标准形J的线性变换阵。证明:若A阵具有重特征值，且每个互异特征值对应一个独立的特征矢量，则必存在一个非奇异阵Q，使A阵化为约当标准形J。即 其中 则

11、J=QAQ-1,53,其中：若Ji为J的约当块，则eJit为（t）中对应的约当块。,54,证明:以Ji有三重特征值为例证明。,此时,55,56,步骤：求 eAt 的方法同对角形求法相一致 1.求i；2.求Qi；3.求eAt=QeJtQ-1,57,四.化 eAt 为A的有限项法：,由于 eAt 可展开无穷级数，但计算时只取有限项，计算结果是不准确的，若能把无穷项级数化成有限项，则计算会简便准确。,1.化有限项的有关理论：,凯哈定理及最小多项式的概念在现代理论中经常用到.下面简要介绍一下有关内容：,1）矩阵A的零化多项式：定理：设有变量s的多项式，矩阵A是nn阶方阵，若满足：,58,则称 为矩阵的

12、零化多项式。2）凯哈密顿定理定理：矩阵的特征多项式是的零化多项式。即：,证明：,59,又因为中各元为(n-1)次多项式，故可一般表示为：代入上式有：用A代替s将上式展开 得,60,3）矩阵A的最小多项式：定义：A的零化多项式中，次数最低的零化多项式称为A的最小多项式。用 表示。的求法：定理：设A的伴随矩阵 全部元素的最大公因子为d(s)则.,61,注：1.该定理证明要用到矩阵多项式的概念.2.计算 要先求。将 各元变为因子相乘的多项式。从中找出各元的最大公因子，且 取首1多项式的形式.,例：已知：试求A的最小多项式并验证凯哈定理。,62,解：1.,63,所以最大因子：故A的最小多项式为：进一步

13、可验证上式是以A为根的零化多项式,64,2.验证凯哈定理：,65,则An可表示成低于n阶幂矩阵的线性组合。,2.eAt 能化成有限项的依据：,由凯哈定理知：矩阵A的特征多项式是A的零化多项式，即,66,由此可推得：,上式表明：对于kn，Ak均可用 An-1，A,I这n个独立项的线性组合来表示。所以可将eAt无穷项化成有限项。,67,故可令：,设n个根为1 2.n，按上式对每个根都有以下结果即,特征方程,第一种情况：A的特征值互异,2.待定系数 的求法,式中，n个待定系数，是t的标量函数。,68,于是对于,其中系数与前面 eAt 的系数相同。,当kn 时，ik的各项均可用 的线性组合表示，得出下

14、列方程组：,69,解此方程组，得系数,例：已知,试用化 为A的有限项法求,70,解：1.求特征值,2.求系数,71,72,3.求,73,74,第二种情况：A有相重特征值 设A有n重特征值1，则按以上方法必有下式,但由于是n重根，不能按同样形式写出n个方程，对上式依次对1求导，直至（n-1）次，可得到（n-1）个导数方程。然后联立这n个方程解出n个待定系数。即,75,解方程组即可求得系数。,76,第三种情况：系统有单根，也有重特征根 设系统矩阵A的特征值中，1为m重特征值，m+1，n为互异的单特征值，根据情况二列写m个方程，根据情况一列写（n m）个方程，解上述n个方程，即可得出系数 的计算公式

15、。,例：已知系统矩阵,试用化eAt为A的有限项法求eAt。,77,2.求系数 i(t)：,解：1.求特征值：,78,即,79,3.求 eAt：,80,可见，以上几例求出的 eAt 中各元都是 的线性组合。,81,2-3线性离散系统的状态空间表达式及连续系统的离散化,一.离散系统的状态空间表达式,1.一般形式。由离散状态方程和离散输出方程组成。,式中：T是采样周期。方程中的矢量，各系数矩阵的名称和维数都与连续系统相同，为简单起见常省去T将方程写成如下形式,82,即：,2.结构图。上述方程可用结构图来表示,83,3.差分方程和脉冲传递函数与离散状态空间表达式之间的转换 在单变量离散系统中，数学模型

16、分为差分方程和脉冲传递函数两类，它们与离散状态空间表达式之间的变换，和连续系统分析相类似。,84,解：1,G(z)差分方程 状态空间表达式,例：已知脉冲传递函数为,试求其状态空间表达式,差分方程为,85,所以,86,2.G（z）部分分式法 状态空间表达式,则,87,3.状态空间表达式 G(z),88,对连续系统，若常用数字计算机进行实时控制或求解，首先必须把连续系统转化成离散系统，这个过程称之为连续系统的离散化。,二.定常连续系统的离散化,离散方程,设定常连续系统,89,连续系统其状态解为：,即取t0=kT,t=(k+1)T,1、直接离散化(精确离散化方法)：离散化的实质就是用一个矩阵差分方程

17、去代替一个矩阵微分方程。,90,在 其输入向量u(t)=u(kT)，初始时刻t0=kT，则状态方程的解为,对第二项积分作变量代换：令t=(k+1)T-；dt=-d上限：=(k+1)T，t=(k+1)T-=0下限：=kT,t=(k+1)T-=T,91,与离散系统的状态空间模型比较可得：,92,例：求 的离散化方程,解：先求eAt：,由拉氏变换法得,93,94,95,2、由脉冲传函实现离散化,步骤：,1首先求连续系统的传递函数,2按照离散系统的结构图求脉冲传函,3按脉冲传函与标准型状态空间表达式的关系写出离散化的状态空间表达式,96,U(s),解：因为离散化后的系统结构图为：,上图的传递函数为：,

18、例1-26 已知连续系统的传递函数为，试求其离散化状态空间表达式,97,对上式取z变换：,最后由G（z）写出其能控标准型的状态空间表达式,98,2-4 离散系统状态方程的解,k=0时,x(1)=Gx(0)+Hu(0)k=1时,x(2)=Gx(1)+Hu(1)=G2X(0)+GHu(0)+Hu(1),一.定常离散系统状态方程的解：（两种方法）,1迭代法：状态方程本身就是一个基本的迭代方程,依次将采样时刻k=0，1，2，3代入上式即可。,已知：初始时刻KT=0，初始状态为x(0),99,几点讨论：,2).第k个采样时刻的状态，只与采样时刻0,1,2k-1时的输入值有关系，而与第k个次采样时刻输入值

19、无关，这是惯性系统的一个基本特征；,由初始状态引起的响应反映系统的自由运动零输入响应由输入引起的响应反映系统的强迫运动零状态响应。,1).定常离散系统的状态解由两部分组成：,100,(k)也满足状态转移矩阵的两个定义条件：矩阵差分方程：(k+1)=G(k)初始条件：(0)=G0=I证明：,3).与连续系统的状态解比较上式中的Gk称为定常离散系统的状态转移矩阵，记为(k)=Gk,101,4).(k)的基本性质,102,序列：,或,5).引入(k)后，状态解又可表示为：,序列：,103,2.z变换法：对方程两边取z变换,与第一种方法比较可知：,求反变换：,104,所以,105,2.迭代法求出的解是

20、一个数值解。只能求出某一时刻的数值。但迭代公式本身就是状态方程，简单方便，而且不用求出状态转移矩阵Gk；如果已求出(k)=Gk，则可用解的迭代公式求出自由分量和强迫分量.,1.z变换求出的解是一个完整解，其中解的结构可分为自由解和强迫解两部分，可分别求出，对分析运动过程有本质的帮助。解的形式是一个闭式，即解析式。,注：,106,例：求线性定常离散系统的解,解：(1)用迭代法求解,已知,107,直至,108,(2)用标准型求Gk,再代入解的迭代公式,也可先求出,109,又知 u(k)=1,于是,（3）用z变换法求解:,先计算(zI G)-1,110,111,令k=0，1，2，3，代入上式，可得,

21、以上两种方法计算结果完全一致，只是迭代法是一个数值解，而z变换法则得到了一个解析表达式。,112,二、离散系统的状态转移矩阵,离散系统状态转移矩阵(k)的求取与连续系统转移矩阵(t)极为类似。,2z变换法 根据z变换法求取离散系统状态方程解中的对应关系，状态转移矩阵(k)为,113,那么，离散系统的状态转移矩阵(k)为,式中，为对角线标准形，若特征方程I G=0的特征根为1，2，n，则有,3化系统矩阵G为标准形法（1）当离散系统矩阵G的特征值均为单根时 当离散系统矩阵G的特征根均为单根时，经过线性变换可将系统矩阵G化为对角线标准形，即,114,式中，P为化系统矩阵G为对角线标准形的变换矩阵。,

22、115,例:齐次离散系统状态方程为,试求其状态转移矩阵(k)。解：,其特征值 1=0.2 2=0.8化系统矩阵G为对角线标准形的变换矩阵P为,116,则,117,(2)当离散系统矩阵G的特征值有重根时,式中J 约当标准形；Q 化系统矩阵G为约当标准形的变换矩阵。,4化为G的有限项法 应用凯莱-哈密尔顿定理，系统矩阵G满足其自身的零化多项式。离散系统状态转移矩阵可化为G的有限项，即,118,式中i(k)(i=0，1，n 1)为待定系数，可仿照连续系统的方法来求取。,例:线性定常离散系统的状态方程为,试求系统的状态转移矩阵(k)。,解：离散系统特征方程为,119,其特征值 1=1 2=2待定系数可

23、按下式求取,解之得,则离散系统状态转移矩阵为,120,121,第二章总结,一、线性定常连续系统1.线性定常系统状态转移矩阵,它包含了系统运动的全部信息，可以完全表征系统的动态特征。（i）定义条件,（ii）求法(1)幂级数法,122,(2)拉氏变换法,(3)对角形法或约当形法,对角形:,123,约当形:,化eAt为A的有限项法 凯莱-哈密尔顿定理：矩阵A满足其本身的零化特征多项式。,124,则有,i(t)的计算按A的特征值互异或有重根时分别计算。,(5)最小多项式法最小多项式为,式中d(s)为伴随矩阵(sI A)各元的最大公因子。则A也要满足其零化的最小多项式，即(A)=0。求eAt的方法与化e

24、At为的有限项法完全相似。,125,2.线性定常系统齐次方程的解可表示为,3.线性定常连续系统非齐次方程的解分为零输入的状态转移和零状态的状态转移，即,二、线性定常离散系统1.线性定常离散系统状态空间表达式,126,2.线性定常连续系统状态方程的离散化1）直接法:采样周期为T，离散化后系统矩阵和输入矩阵分别为,二、线性定常离散系统1.线性定常离散系统状态空间表达式,127,2)脉冲传递函数法:,状态空间表达式,其中离散化过程是通过求脉冲传函来完成的.,128,(1)状态转移矩阵(k),c.性质:,a.定义:,b.求法:,3.线性定常离散系统状态方程的解,或,129,(2).状态方程的解1)迭代法 a.状态方程法 把初始条件和输入函数直接代入状态方程表达式即可。,b.状态解法 如果已求出状态转移矩阵Gk,则把初始条件和输入函数直接代入状态解表达式即可。,130,2)z变换法,例:线性定常连续系统的状态空间表达式为,设采样周期T=1s，求离散化后系统的离散状态空间表达式。,131,解：先求连续系统状态转移矩阵,132,133,离散化后系统的离散状态空间表达式为,134,结 束,