状态空间分析方法精选PPT.ppt
《状态空间分析方法精选PPT.ppt》由会员分享,可在线阅读,更多相关《状态空间分析方法精选PPT.ppt(216页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、关于状态空间分析方法第1页,讲稿共216张,创作于星期二第第9 9章章 状态空间分 析方法基本要求9-1 状态空间方法基础9-2 线性系统的可控性和可观性9-3 状态反馈和状态观测器9-4 有界输入、有界输出的稳定性9-5 李雅普诺夫第二方法返回主目录第2页,讲稿共216张,创作于星期二引言引言:前面几章所学的内容称为经典控制理论;下面要学的内容称为现代控制理论。两者作一简单比较。经典控制理论经典控制理论(50年代前年代前)现代控制理论现代控制理论(50年代后年代后)研究对象研究对象单输入单输出的线单输入单输出的线性定常系统性定常系统可以比较复杂可以比较复杂数学模型数学模型传递函数传递函数(输
2、入、输出描述输入、输出描述)状态方程状态方程(可描述内部行为可描述内部行为)数学基础数学基础运算微积、复变函运算微积、复变函数数线性代数、矩阵理论线性代数、矩阵理论设计方法的设计方法的特点特点非唯一性、试凑成非唯一性、试凑成份多份多,经验起很大作经验起很大作用。主要在复数域用。主要在复数域进行。进行。设计的解析性,与计设计的解析性,与计算机结合,主要在时算机结合,主要在时间域进行。间域进行。第3页,讲稿共216张,创作于星期二基本要求基本要求 掌握由系统输入输出的微分方程式、系统动态结构图、及简单物理模型图建立系统状态空间模型的方法。熟练掌握矩阵指数的计算方法,熟练掌握由时域和复数域求解状态方
3、程的方法。熟练掌握由动态方程计算传递函数的公式。正确理解可逆线性变换,熟练掌握可逆线性变换前、后动态方程各矩阵的关系。正确理解可控性和可观测性的概念,熟练掌握和运用可控性判据和可观性判据。返回子目录返回子目录第4页,讲稿共216张,创作于星期二熟练掌握可逆线性变换矩阵的构成方法,能将可控系统 化为可控标准形。能将不可控系统进行可控性分解。正确理解对偶原理,会将原系统的有关可观测性的问题转化为对偶系统的可控性问题来研究。正确理解单变量系统零、极点对消与动态方程可控、可观测的关系。熟练掌握传递函数的可控性标准形实现、可观性标准形实现的构成方法。正确理解状态反馈对可控性,可观性的影响,正确理解状态反
4、馈可任意配置闭环极点的充要条件。第5页,讲稿共216张,创作于星期二熟练掌握全维状态观测器的公式和设计方法,熟练掌握由观测器得到的状态估计值代替状态值构成的状态反馈系统,可进行闭环极点配置和观测器极点配置。正确理解系统齐次方程渐近稳定和系统BIBO稳定的概念,熟练掌握判别渐近稳定的方法和判别系统BIBO稳定的方法。k正确理解李雅普诺夫方程正定对称解存在的条件和解法,能通过解李雅普诺夫方程进行稳定性分析。第6页,讲稿共216张,创作于星期二9-1 状态空间方法基础在经典控制理论中,用传递函数来设计和分析单输入、单输出系统。在现代控制理论中,用状态变量来描述系统。采用矩阵表示法可以使系统的数学表达
5、式简洁明了,为系统的分析研究提供了有力的工具。返回子目录返回子目录第7页,讲稿共216张,创作于星期二状态:状态:动力学系统的状态可以定义为信息的集合。动力学系统的状态可以定义为信息的集合。一、状态空间的基本概念已知已知 时状态,时状态,时的输入,可确定时的输入,可确定 时任一变量的运动状况。时任一变量的运动状况。状态变量状态变量:确定动力学系统状态的最小一组变量确定动力学系统状态的最小一组变量 。第8页,讲稿共216张,创作于星期二状态空间:由 张成的n维向量空间。状态向量状态向量:如果完全描述一个给定系统的动如果完全描述一个给定系统的动态行为需要态行为需要n n个状态变量,那么状态个状态变
6、量,那么状态向量定义为向量定义为X(t)X(t)对于确定的某个时刻,状态表示为状态空间中一个点,状态随时间的变化过程,构成了状态空间中的一条轨迹。第9页,讲稿共216张,创作于星期二例9-2设一RLC网络如图所示。回路方程为图9-2 RLC网络第10页,讲稿共216张,创作于星期二选择状态变量则有写成输出第11页,讲稿共216张,创作于星期二写成若选另一组状态变量则有第12页,讲稿共216张,创作于星期二 若给出(t=0)时的初值 、和 时就可确定系统的行为。单输入单输入-单输出线性定常系统单输出线性定常系统选取状态变量二、系统的状态空间表达式第13页,讲稿共216张,创作于星期二(9-17)
7、第14页,讲稿共216张,创作于星期二或写成(9-19)第15页,讲稿共216张,创作于星期二系统结构图如图所示图9-3第16页,讲稿共216张,创作于星期二例9-3输入为输入为 u u,输出为,输出为y y。试求系统的状态方程和输出方程。试求系统的状态方程和输出方程。考虑用下列常微分方程描述的系统考虑用下列常微分方程描述的系统第17页,讲稿共216张,创作于星期二解:状态方程为写成取状态变量第18页,讲稿共216张,创作于星期二输出图9-4 例9-3系统的结构图第19页,讲稿共216张,创作于星期二多输入-多输出系统图9-6 多变量系统第20页,讲稿共216张,创作于星期二 为状态变量;为输
8、入量;为输出变量。第21页,讲稿共216张,创作于星期二矩阵形式:式中第22页,讲稿共216张,创作于星期二.输出变量方程第23页,讲稿共216张,创作于星期二式中式中第24页,讲稿共216张,创作于星期二图9-7 系统结构图第25页,讲稿共216张,创作于星期二三、线性定常系统状态方程的解式中式中 均为列向量。均为列向量。(9-28)齐次向量微分方程齐次向量微分方程(9-29)方程的解为方程的解为1、齐次状态方程的解第26页,讲稿共216张,创作于星期二可得代入方程 将方程两边系数必相等方程两边系数必相等,即即第27页,讲稿共216张,创作于星期二我们定义(9-31)(9-32)因此,齐次状
9、态方程的解为将 t=0 代入(9-29)中得第28页,讲稿共216张,创作于星期二(9-33)(9-34)(9-35)为nn矩阵,称矩阵指数。于是齐次状态方程的解为于是齐次状态方程的解为用拉氏变换法求解用拉氏变换法求解第29页,讲稿共216张,创作于星期二拉氏反变换后得到(9-37)(9-38)第30页,讲稿共216张,创作于星期二最终得到与前一种解法所得结果一致。式中(9-41)第31页,讲稿共216张,创作于星期二状态转移矩阵具有以下性质:状态转移矩阵具有以下性质:第32页,讲稿共216张,创作于星期二图9-8 状态转移特性性质性质3第33页,讲稿共216张,创作于星期二例9-5设系统的状
10、态方程为设系统的状态方程为试求状态转移矩阵。试求状态转移矩阵。第34页,讲稿共216张,创作于星期二解:求状态转移矩阵为其中可以写出方程解为第35页,讲稿共216张,创作于星期二例9-6设系统状态方程为设系统状态方程为试求状态方程的解。试求状态方程的解。第36页,讲稿共216张,创作于星期二解:用拉氏变换求解。先求出矩阵指数用拉氏变换求解。先求出矩阵指数 第37页,讲稿共216张,创作于星期二状态方程之解为 将上式进行拉氏反变换将上式进行拉氏反变换第38页,讲稿共216张,创作于星期二图9-9 系统的瞬态解(a)与相轨迹(b)第39页,讲稿共216张,创作于星期二改写为 用 左乘等式两边 2
11、2 非齐次状态方程的解非齐次状态方程的解非齐次方程(9-53)(9-54)第40页,讲稿共216张,创作于星期二用 左乘上式两边(9-54)则式(9-54)可以写成(9-55)积分上式得第41页,讲稿共216张,创作于星期二讨论非齐次状态方程的拉氏变换解法讨论非齐次状态方程的拉氏变换解法拉氏反变换得拉氏反变换得由于由于由卷积定理有由卷积定理有第42页,讲稿共216张,创作于星期二因此由于由于最后得到第43页,讲稿共216张,创作于星期二例9-7求下述系统状态的时间响应求下述系统状态的时间响应控制量控制量u u为单位阶跃函数。为单位阶跃函数。第44页,讲稿共216张,创作于星期二解:由状态转移矩
12、阵第45页,讲稿共216张,创作于星期二若初始状态为零状态,则若初始状态为零状态,则第46页,讲稿共216张,创作于星期二四、传递函数矩阵(9-58)系统状态方程系统状态方程(9-59)输出方程输出方程拉氏变换为拉氏变换为第47页,讲稿共216张,创作于星期二解出解出定义传递函数矩阵为(9-63)第48页,讲稿共216张,创作于星期二所以所以特征方程为第49页,讲稿共216张,创作于星期二例9-8设系统的动态方程为试求该系统的传递函数矩阵。第50页,讲稿共216张,创作于星期二解:已知已知故故第51页,讲稿共216张,创作于星期二第52页,讲稿共216张,创作于星期二例9-9设系统的状态方程为
13、设系统的状态方程为试求系统的特征方程和特征值。试求系统的特征方程和特征值。第53页,讲稿共216张,创作于星期二解:系统的特征方程为系统的特征方程为特征方程的根为-1、-2和-3。矩阵A的特征值也为-1、-2和-3。两者是一样的。第54页,讲稿共216张,创作于星期二五、动态方程的可逆线性变换五、动态方程的可逆线性变换其中 P 是nn 矩阵第55页,讲稿共216张,创作于星期二特征多项式特征多项式没有改变。第56页,讲稿共216张,创作于星期二传递函数阵传递函数阵传递函数阵没有改变传递函数阵没有改变第57页,讲稿共216张,创作于星期二例9-10对例9-9之系统进行坐标变换,其变换关系为试求变
14、换后系统的特征方程和特征值。第58页,讲稿共216张,创作于星期二解:根据题意求变换矩阵代入第59页,讲稿共216张,创作于星期二特征方程为特征值为-1,-2,-3,与例9-9结果相同。可得第60页,讲稿共216张,创作于星期二9-2 9-2 线性系统的可控性和可观测性线性系统的可控性和可观测性在状态空间法中,对系统的描述可由状态方程和输出方程来表示。状态方程是描述由输入和初始状态所引起的状态的变化;输出方程则是描述由于状态变化而引起输出的变化可控性和可观测性的概念,就是回答可控性和可观测性的概念,就是回答“系统的输入是系统的输入是否能控制状态的变化否能控制状态的变化和和“状态的变化能否由输出
15、状态的变化能否由输出反映出来反映出来这样两个问题。这样两个问题。返回子目录返回子目录第61页,讲稿共216张,创作于星期二一、准备知识一、准备知识设设A A 是 nn 矩阵,x x 是 n1 向量,齐次方程组若|A|=0,(9-70)式存在非零解;若|A|0,(9-70)式只有零解。Ax=0(9-70)1 1、齐次方程组的非零解、齐次方程组的非零解第62页,讲稿共216张,创作于星期二2、Cayley-Hamilton定理 Cayley-Hamilton定理指出,矩阵A满足自己的特征多项式。则A满足(9-71)(9-72)A的特征多项式第63页,讲稿共216张,创作于星期二应用Cayley-H
16、amilton 定理(9-78)对于矩阵指数 可以用来表示。第64页,讲稿共216张,创作于星期二例9-11解:矩阵A的特征多项式要求计算矩阵 的第65页,讲稿共216张,创作于星期二矩阵A满足自己的特征多项式,有本题中n=100,故有第66页,讲稿共216张,创作于星期二3 引理的充分必要条件是:的充分必要条件是:存在存在 使使(9-80)非奇异。这里非奇异。这里A:nn,b:n1.A:nn,b:n1.第67页,讲稿共216张,创作于星期二若对任意状态若对任意状态 ,存在一个有限时刻,存在一个有限时刻 和控制量和控制量 ,能在能在 时刻将状态时刻将状态 转移到转移到0 0,则称此系统的状态完
17、全可控。,则称此系统的状态完全可控。二、线性系统的可控性二、线性系统的可控性1 定义对于任意时刻对于任意时刻 和和 ,若存在控制向量,若存在控制向量 ,能将,能将 的的每个初始状态每个初始状态 转移到转移到 时刻的另一任意状态时刻的另一任意状态 ,则称此系统的状态完全可控。则称此系统的状态完全可控。等价的定义第68页,讲稿共216张,创作于星期二例如图9-10 二维系统状态转移过程如图所示二维系统状态转移过程如图所示系统可控。系统可控。第69页,讲稿共216张,创作于星期二2 可控性判据其中 A(nn),b(n1),c(1n),d(11)系统可控的充分必要条件是(9-84)(9-85)(9-8
18、6)单变量线性定常系统第70页,讲稿共216张,创作于星期二证明:将u(t)代入式(9-54),可得(9-87)若式若式(9-86)(9-86)成立,由前面准备知识的引理,存在成立,由前面准备知识的引理,存在t t1 100,使得,使得(1-(1-30)30)式定义的式定义的W(0,tW(0,t1 1)矩阵非奇异,取矩阵非奇异,取t t1 1为可控性定义中的为可控性定义中的t tf f ,且在且在0,t0,tf f 上定义上定义第71页,讲稿共216张,创作于星期二由定义可知式(9-86)成立时,系统可控。第72页,讲稿共216张,创作于星期二再证明若系统可控,则式(9-86)成立 根据凯莱哈
19、密尔顿定理(9-88)(9-89)假定系统由任意初始状态被控制到零状态,即 x(tf)=0。根据(9-54)式,则有第73页,讲稿共216张,创作于星期二把(9-89)式代入(9-88)式,得记这时(9-90)第74页,讲稿共216张,创作于星期二由于x(0)是任意的n维向量,(9-90)式要有解,一定有(9-86)式成立,即由上述可控性判据可知,系统的可控性只取决于由上述可控性判据可知,系统的可控性只取决于(9-84)(9-84)式中的式中的A A阵和阵和b b阵。今后为了方便起见,将可控性矩阵记为阵。今后为了方便起见,将可控性矩阵记为S S,这样,这样,可控的充要条件就写成:可控的充要条件
20、就写成:rankS=n rankS=n 或或 detS0detS0。第75页,讲稿共216张,创作于星期二图9-11 不可控系统第76页,讲稿共216张,创作于星期二例子系统可控系统第77页,讲稿共216张,创作于星期二3 约当型方程的可控性判据 约当块的一般形式为由前面讨论可知,等价变换不改变可控性。第78页,讲稿共216张,创作于星期二可控的充分必要条件为同一特征值对应的约当块只有一块,即各约当块的特征值不同。每一约当块最后一行,所对应的b中的元素不为零。这一充分必要条件又称为单输入系统约当形方程的可控性判据。第79页,讲稿共216张,创作于星期二例9-12系统状态方程为系统状态方程为试确
21、定系统可控时,试确定系统可控时,应满足的条件。应满足的条件。第80页,讲稿共216张,创作于星期二解:如果用直接计算可控性矩阵的方法也可得到同样结果.因为因为A A阵有两个若当块,根据判据的阵有两个若当块,根据判据的(1)(1)应有应有 ,由判据的,由判据的(2)(2),A A的第二行所对应的的第二行所对应的b b中的元中的元素素b2 2,b4 4均不为零,因此系统可控的充要条件均不为零,因此系统可控的充要条件为为第81页,讲稿共216张,创作于星期二4、可控标准形(9-92)则系统一定可控。一个单输入系统,如果具有如下形式第82页,讲稿共216张,创作于星期二(9-92)式的形式被称为单输入
22、系统的可控标准形可控标准形。对于一般的单输入n维动态方程 (9-93)其中A,b分别为nn,n1的矩阵。成立以下定理:若n维单输入系统可控,则存在可逆线性变换,将其变换成可控标准形。第83页,讲稿共216张,创作于星期二下面给出变换矩阵P的构成方法 计算可控性矩阵S;计算 ,并记 的最后一行为h。构造矩阵 P令 即可求出变换后的系统状态方程。即可求出变换后的系统状态方程。第84页,讲稿共216张,创作于星期二例9-13设系统状态方程为 试将系统状态方程化为可控标准形。第85页,讲稿共216张,创作于星期二解:先判断可控性,再计算变换矩阵,将状态方程化为可控标准形。故系统可控。一定可将它化为可控
23、标准形。第86页,讲稿共216张,创作于星期二此时标准形中的系统矩阵的最后一行系数就是A阵特征式的系数,但符号相反。则变换矩阵为第87页,讲稿共216张,创作于星期二可求出第88页,讲稿共216张,创作于星期二5 系统按可控性进行分解系统按可控性进行分解 系统可控时,可通过可逆线性变换变换为可控标准形,现在研究不可控的情况,这时应有下面的结果被称为系统按可控性进行分解的定理 第89页,讲稿共216张,创作于星期二若单变量系统(9-84,85)式的可控性矩阵满足(9-103)式,则存在可逆线性变换矩阵P,使得变换后的系统方程具有以下形式 式中 是n1维向量,是n2维向量,并且(9-106)(9-
24、107)第90页,讲稿共216张,创作于星期二(9-106)式表明下面的动态方程是可控的:(9-107)式表明的动态方程式(9-108,109)和原来的n维动态方程式(9-84,85)具有相同的传递函数。或者说传递函数中未能反映系统中不可控的部分。(9-108)(9-109)第91页,讲稿共216张,创作于星期二证明:证明:(9-110)考察考察(9-103)(9-103)式,并将它重新写出如下式,并将它重新写出如下进而可以证明进而可以证明补充选取线性无关的向量补充选取线性无关的向量并使得向量组并使得向量组 线性无关。线性无关。第92页,讲稿共216张,创作于星期二令若将若将(9-104(9-
25、104,105)105)式所表示的系统用方框图表示,可控式所表示的系统用方框图表示,可控性分解的意义就能更直观地体现出来,性分解的意义就能更直观地体现出来,(9-104(9-104,105)105)式式的系统方块图如图的系统方块图如图9-129-12所示。所示。即可证明 具有定理所要求的(9-104)的形式。第93页,讲稿共216张,创作于星期二图9-12 系统按可控性分解第94页,讲稿共216张,创作于星期二从图9-12中可见,控制输入不能直接改变 也不能通过影响 间接改变 ,故 这一部分状态分量是不受输入影响的,它是系统中的不可控部分。由图上还可看出系统的传递函数完全由图中虚线以上的部分所
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 状态 空间 分析 方法 精选 PPT
限制150内