矩阵函数的性质及其应用
- 格式:doc
- 大小:360.50 KB
- 文档页数:9
§7 矩阵函数的性质及其应用一、矩阵函数的性质: 设n n C B A ⨯∈. 1.A e Ae e dtd At At At⋅== proof : 由 ()∑∑⋅==∞=m m m m AtA t m At m e!1!1对任何t 收敛。
因而可以逐项求导。
()∑∞=--=∴01!11m mm At A t m e dt d ()()⎪⎪⎭⎫ ⎝⎛-⋅=∑∞=-11!11m m At m A ()⎪⎪⎭⎫ ⎝⎛⋅=∑k At k A !1AteA ⋅=()()()A e A At m A A t m Atm m m m m ⋅=⋅⎪⎪⎭⎫ ⎝⎛-=⋅-=∑∑∞=∞=---01111!11!11 可见,A 与At e 使可以交换的,由此可得到如下n 个性质 2.设BA AB =,则 ①.At At Be B e =⋅ ②.B A A B B A e e e e e +=⋅=⋅③.()()AA A AA AB A B A B A BA B A B A BA cos sin 22sin sin cos 2cos sin cos cos sin sin sin sin cos cos cos 22=-=⇒+=+-=+= proof :①,由m m BA B A BA AB =⇒=而∑∑∞=∞==⎪⎭⎫⎝⎛=00!1!1m m m m m m AtB A t m B t A m B e()∑∑∞=∞=⋅==00!1!1m mm m m At m B BA t m At e B ⋅=②令()()A B t At Bt C t e e e +--=⋅⋅ 由于()0=t C dtd)(t C ∴为常数矩阵 因而E e e e C C t C =-⋅===000)0()1()(当1=t 时,E e e e B A B A =⋅⋅--+ …………………. (@) 特别地 A B -= 有E e e e A A =⋅⋅-0∴ 有 ()A Ae e --=1∴同理有()B Be e --=1代入(@)式 因而有B A B A e e e ⋅=+ 3.利用绝对收敛级数的性质,可得 ①A i A e iA sin cos +=()()iA iAiA iAe e iA e e A ---=+=⇒21sin 21cos ②()()A A A A sin sin cos cos -=-=-4.E A A =+22cos sin()()A E A AE A cos 2cos sin 2sin ππ+=+A E i A e e =+π2二、矩阵函数在微分方程组中的应用—常用于线性监测系统中 1. 一阶线性常系数齐次方程组的通解 AX dtdX = 其中()Tn n n x x x X C A ,,,21 =∈⨯ 则有()K e t X At ⋅=,其中()Tn k k k K ,,,21 =1eg解方程:⎪⎪⎪⎩⎪⎪⎪⎨⎧+=+-=+-=313212211234xx dtdx x x dtdxx x dt dx解:原方程变为矩阵形式AX dt dX = ⎪⎪⎪⎭⎫⎝⎛--=201034011A ()Tx x x X 321,,= 由()()212--=-λλλA E 得⎪⎪⎪⎭⎫ ⎝⎛=→100110002J A1200000-⎪⎪⎪⎭⎫⎝⎛=∴P e e e e P e t tt tAt⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎭⎫ ⎝⎛=∴-3211200000)(k k k P e e e e P t X t tt t2. 一阶线性常系数微分方程组的定解问题:1Th :一阶线性常数微分方程组的定解问题:()()⎪⎩⎪⎨⎧==Tn x x x X AXdt dX)0(,),0(),0(210 有唯一解)0(X e X At ⋅= proof :实际上,由AX dtdX=的通解为K e t X At ⋅=)( 将初值)0(X 代入,得)0(X k =)0(X e X At =∴由1Th 可的定解问题()⎪⎩⎪⎨⎧==Tn t x t x t x t X AX dt dX)(,),(),()(002010的唯一解为()()00)(t X e t X t t A ⋅=-2eg 求定解问题:()()⎪⎩⎪⎨⎧==Tx Axdt dx1,00,⎪⎪⎭⎫ ⎝⎛--=1221A 的解解:由0=-A E λ 得i x 32,1±=对应的特征向量记为:Ti ⎪⎪⎭⎫ ⎝⎛+=231,1α ⎪⎪⎭⎫⎝⎛-=231,1i β 则,于是矩阵:⎪⎪⎭⎫ ⎝⎛-+=23123111i iP13300--⋅⎪⎪⎭⎫ ⎝⎛⋅=∴P e e P eit itAt⎪⎪⎪⎪⎭⎫⎝⎛+=⎪⎪⎭⎫ ⎝⎛=t t t e t X At 3sin 313cos 3sin 3210)(练习:求微分方程组1132123313383625dx x x dt dx x x x dt dx x x dt ⎧=+⎪⎪⎪=-+⎨⎪⎪=--⎪⎩满足初始条件123(0)1,(0)1,(0)1x x x ===的解。
解:令1233081316,(),(0)12051x A x t x x x ⎛⎫⎡⎤⎡⎤ ⎪⎢⎥⎢⎥=-== ⎪⎢⎥⎢⎥ ⎪⎢⎥⎢⎥--⎝⎭⎣⎦⎣⎦ 可求得3d e t ()(1)I A λλ-=+,而A 的最小多项式2()(1)A m λλ=+。
可设01()g c c λλ=+,由()()010111(1)1tt t tg c c ec t e g c te c te ----⎧⎧-=-==+⎪⇒⎨⎨'-===⎪⎩⎩01At e c I c A ∴=+14083162014t t t e t t t t -+⎛⎫ ⎪= ⎪ ⎪--⎝⎭;)0(X e X At =∴112916t t e t t -+⎛⎫⎪= ⎪ ⎪-⎝⎭3.一阶常系数非齐次方程组的定解问题:()()⎪⎩⎪⎨⎧=+==00t x x t F Ax dt dxt t 其中()()()()()Tnt F t F t F t F ,,,21 = ()t F Ax dtdx=- 两边同乘以At e -得:()()t F e x e dtd At At--= 从0t 到t 上积分得:()()ττd F e t x et x ett AE At At⎰---=-00)(()()()()τττd F e t x et x tt t A t t A ⎰--+=∴000)(3eg .求:非齐次微分方程组()()⎪⎩⎪⎨⎧=+=T x t F AX dtdx1,0)0(的解: 其中⎪⎪⎭⎫⎝⎛-=3553A ()⎪⎪⎭⎫ ⎝⎛=-0t e t F 解:由i A E 5302,1±=⇒=-λλ对应特征向量为:⎪⎪⎭⎫⎝⎛=i 1α ⎪⎪⎭⎫ ⎝⎛=1i β 得可逆矩阵⎪⎪⎭⎫ ⎝⎛=11i i P ⎪⎪⎭⎫ ⎝⎛--=-11211i i P()()ti i Ate t t t t P e e P e3153535cos 5sin 5sin 5cos 00⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛=∴--+ ()⎰⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=--t t A Atd e e e t x 0010)(τττ τd e t t t t t t t t e t t e tt t t 40335sin 5cos 5cos 5sin 5sin 5sin 5cos 5cos 5cos 5sin -⎰⎪⎪⎭⎫ ⎝⎛+-++⎪⎪⎭⎫ ⎝⎛= 练习:求微分方程组112212313214221t dx x x dt dx x x dt dx x x e dt ⎧=-++⎪⎪⎪=-++⎨⎪⎪=++-⎪⎩满足初始条件123(0)1,(0)1,(0)1x x x ===-的解。
解:令12321011420,()2,(),(0)110111t x A F t X t x X e x -⎛⎫⎡⎤⎡⎤⎡⎤ ⎪⎢⎥⎢⎥⎢⎥=-=== ⎪⎢⎥⎢⎥⎢⎥⎪⎢⎥⎢⎥⎢⎥--⎝⎭⎣⎦⎣⎦⎣⎦ 可求得2det()(1)I A λλλ-=-。
可设2012()g c c c λλλ=++,由()()()001120120110111t g c c c t g c t c e t g c c c ⎧===⎧⎪⎪'===⇒⎨⎨⎪⎪=--=++=⎩⎩2012At e c I c A c A ∴=++1204120121t t t tt tt t e e t e ⎛⎫-⎪=-+ ⎪ ⎪+---⎝⎭;(0)At e X ∴112t t t t e -⎛⎫ ⎪=- ⎪ ⎪-⎝⎭,()A e F ττ-1201141202210121e e e e ττττττττττ-⎛⎫--⎡⎤⎡⎤⎪⎢⎥⎢⎥=-=⎪⎢⎥⎢⎥⎪⎢⎥⎢⎥---+-⎣⎦⎣⎦⎝⎭, 0()220tAt A At tt t e e F d e t t te t τττ-⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎰。
故:()()()01()01(1)tA t At t X t e X e F d t e τττ-⎡⎤⎢⎥∴=+=⎢⎥⎢⎥-⎣⎦⎰ 三、矩阵分析在求方程组最小二乘解等问题的应用。
例4 设,(1,2,,)m n m i i A R b R i k ⨯∈∈= ,证明:(0)n x R ∈为函数:221()ki ii f x A x b ==-∑的极小值点的充要条件是(0)n x R ∈为方程组11()k kTT i i i i i i A A x A b ===∑∑的解或方程组1122k k A b A b x A b ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦* 的最小二乘解。
证明:必要性:由于1()()()kT i i i i i f x A x b A x b ==--∑1()k T T T T T T i i i i i i i i i x A A x x A b b A x b b ==--+∑由于(0)x 为()f x 的极小点,则应有(0)(0)1(22)0kT T i i i i x i df A A x A b dx==-=∑即 (0)11k k T Ti i i i i i A A x A b ==⎛⎫= ⎪⎝⎭∑∑这就是说(0)x 是方程组*的代数方程组11()kkTT ii i i i i A A x A b ===∑∑的解,也就是方程组*的最小二乘解。
充分性:(0)x 是方程组*的最小二乘解,根据定义,它应该是函数()f x 的极小点。
练习:设,,,m n n k n k A R b R B R d R ⨯⨯∈∈∈∈,且rankA n =,方程Bx d =有解,试求约束极小化问题22min Bx dAx b =-的解,也就是求函数22()i i f x A x b =-在约束Bx d =下的极小点和极小值。