状态方程的解.ppt
《状态方程的解.ppt》由会员分享,可在线阅读,更多相关《状态方程的解.ppt(59页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、状态方程的解 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life, there is hope。有生命必有希望。有生命必有希望2 2 状态方程的解状态方程的解1 1、相关知识、相关知识032! 3! 21kkktktktttte1!nnsntaseat100)()()(xtxtAxtx )(tx齐00)()()()(xtxtButAxtx 特齐非齐xtxtx)()(2 2、齐次状态方程的解、齐次状态方程的解0)0()()(xxtAxtx 若为标量微分方程:若为标量微分方程:00)0()()(xxsaXxssX拉氏变换拉氏变换a
2、sxsX0)( 000!)(xkatxetxkkat拉氏反变换拉氏反变换标量微分方程为矩阵微分方程标量微分方程为矩阵微分方程n=1n=1的特例的特例000!)(xktAxetxkkkAt032! 3! 21kkktktktttte1!nnsntaseat1(1 1)说明)说明110!k kAtkA teLsIAk000)( 0 xetxtttA ,当0ttAe称状态转移矩阵,将状态由称状态转移矩阵,将状态由t t0 0时刻的时刻的x(tx(t0 0) )转移到转移到t t时刻时刻的的x(t)x(t)状态由状态由t t时刻转移到时刻转移到t t0 0时刻的转移矩阵时刻的转移矩阵00ttAttAe
3、e(2 2)矩阵指数的性质)矩阵指数的性质Ie 0A tAtAeee1AtAtee()A B tAtBtABBAeee设设P P是与是与A A同阶的非奇异矩阵,则有同阶的非奇异矩阵,则有PePeAtAPtP11求求eAt的一的一种方法种方法11PPeeAPtPAtkAtkAteeAtAtAtdeAee Adt矩阵指数的性质矩阵指数的性质201021()()()210 A ttA ttA tteeettt传递性传递性意义:解可以分段求意义:解可以分段求(3 3)几个特殊的状态转移矩阵)几个特殊的状态转移矩阵( (可直接使用!可直接使用!) )100nAA A为对角阵为对角阵100ntAtteee
4、2121 2!1 !0 1 (2)! 0 0 0 1mmAtttttmteetm几个特殊的状态转移矩阵几个特殊的状态转移矩阵A A为约当块为约当块100100m mA 几个特殊的状态转移矩阵几个特殊的状态转移矩阵A A为约当矩阵为约当矩阵100kJAJJ100kJ tAtJ teee几个特殊的状态转移矩阵几个特殊的状态转移矩阵A A通过非奇异阵通过非奇异阵P P化为对角形矩阵化为对角形矩阵APP1111PPePPeetAPtPAt AtAtettetcos sinsint cos (4)(4)矩阵指数的计算方法矩阵指数的计算方法: :直接法直接法( (一般不用一般不用) )0!kkkAtktA
5、e一般用于计算机求解一般用于计算机求解当当A A为特殊矩阵,如幂零阵,可用此法为特殊矩阵,如幂零阵,可用此法0 0 01 0 00 1 0A0 0 00 0 01 0 02A03A2 22220 0 1 1 0 00 0220 1 00 0 0 0 00 1 t0 0 10 0 00 0 00 0 1AtA teIAtttttt底友零阵,实质是特征根为底友零阵,实质是特征根为0 0的约当块的约当块例例3322! 31! 21tAtAAtIeAt232310000110102!3!00tttttt0110AAte已知已知 求求243535241cossin24!3!5!sincos13!5!2!
6、4!tttttttttttttt 矩阵指数的计算方法矩阵指数的计算方法: :拉氏变换法拉氏变换法( (常用常用) )11AteLsIAnnnnnnnnnasasasasBsBsBsBAsIAsIAsI1221112211*1nkIaABBnkABtrkaIBkkkkk, 3 , 2 , 2 , 1 1 111迭代公式迭代公式例例11220123xxxx321321000ssssAsI*12111311121222212121212ssIAsssssIAssIAssssss ttttttttAteeeeeeeeAsILe2222112222例例0 0 01 0 00 1 0A3221332211
7、asasasBsBsBAsI 1112112222322, 0,0 0 11110 0 002220 0 00 0 10 0 00 0 0BIatr ABtr ABABa IAatr ABtr AtrBABa IA 333433110,033atr ABtr ABABa I sssssssssAsI1 0 01 1 01 1 10 0 00 0 01 0 00 0 01 0 00 1 01 0 00 1 00 0 112322311 0 0 1 021 1211tttAsILeAt矩阵指数的计算方法矩阵指数的计算方法: :凯莱凯莱- -哈密尔顿法哈密尔顿法 0111nnnnaaaAIf本质:化
8、本质:化无穷级数无穷级数为为有限项之和有限项之和(1) (1) 凯莱凯莱- -哈密尔顿定理哈密尔顿定理12121nnnnnAa Aa AaAa I 的线性组合可表示成IAAAAnnn,21 0111IaAaAaAAfnnnnAaAaAaAaAAAnnnnnn211211AaAaAaIaAaAaAaannnnnnn2112122111IaaAaaaAaaaAaannnnn11123211221的线性组合可表示成IAAAAnnn,211线性表示都可由IAAAAAAnnnnn,212112121nnnnnAa Aa AaAa I 22! 21tAAtIeAt的线性组合可表示为IAAAennAt,21
9、111010)()()()(nnnkkkAtAtAtItAte(2) (2) 化化 为为A A的有限项的有限项Ate(3) (3) 的计算的计算)(tk1) A特征值互异时特征值互异时地位等同,则与由iA10( ) 1,2,intkkiketin112210121222210111212110)()()()()()()()()()()()(21nnnnntnntnntttttettttetttten112210121222210111212110)()()()()()()()()()()()(21nnnnntnntnntttttettttettttentttnnnnnnnneeettt2111
10、21222211211110 1 1 1)()()(求解步骤求解步骤n,0A-I ) 1 (21由10(2) ( )1,2,intkkiketin代入)( ) 3(tk求出10)( )4(nkkkAtAte写出例例AteA,求3- 2-1 0 212 -1(1) 320122 3IA ,1202221( )1 -12 -12(2) 1 -21 -1( )ttttttttteeeeteeee2012222222 0(3) ( )( ) 0 2 0 02 -2 -3-22 -ttAttttttttttttteeet It Aeeeeeeeeeeee2) A有重特征根A有一个重特征根有一个重特征根
11、为为m重重1nmm,21为单根,为单根,(n-m)个单根个单根mnnmmtmnntnntnnttmnntmtmettnnttettntttetttte111113111322211121111212110)(!1)(! 1!)(! 0!1)(21)(6)(2)(1)(2)()()()()(11111110121210111110)()()()()()()()()(21nnnntnmnmtnmnmttttetttetttenmm例例xx2- 1-1 02 -12101 2IA (两重)121)()()(111011ttettettttttetteet)()(1001 0 0 ( )( ) 0 2
12、 tttAttttttttttteteteet It Aeteteteeteteteete矩阵指数的计算方法矩阵指数的计算方法: :特征值与特征向量法特征值与特征向量法任一矩阵任一矩阵A A都可通过非奇异线性变换阵都可通过非奇异线性变换阵P P变成对角形或约当形。变成对角形或约当形。APP10 AIi0iiAI根据线性代数知识,有:根据线性代数知识,有:任意矩阵任意矩阵A A都可化为约当标准形(含对角形);都可化为约当标准形(含对角形);A A化对角形化对角形 A A有有n n个线性无关的特征向量;个线性无关的特征向量;A A化对角形化对角形 A A有有n n个不同的特征值(不同特征值对应特征
13、向量线性无个不同的特征值(不同特征值对应特征向量线性无关,这样就有关,这样就有n n个线性无关特征向量)个线性无关特征向量)111PPePPeetAPtPAti特征值特征值i对应于特征值对应于特征值 的特征向量的特征向量inP,21ttAPtPnneeeAPP00,001111110011PeePPPeettAPtPAtnA A有有n n个不同特征值个不同特征值n,21例例3- 2-1 0AAte求求 2 , 10233 21- 212AI(1) 求特征值求特征值(2)(2)求特征向量求特征向量002 2 1- 1-212111xxxxAI1121 xx令11 10201 2 1- 2-212
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 状态方程
限制150内