关于常系数齐次线性微分方程组通解的两种方法.docxVIP

关于常系数齐次线性微分方程组通解的两种方法 在一般方程教材中，关于线性差分矩阵的同解解释方法通常有两种方法。其中之一是直接计算矩阵指数函数的eat。另一种方法是从矩阵的资源方程出发，计算相应的资源向量，并根据资源值的不同重数求出同解。在这项工作中，我们给出了上述两种方法的不同解法，以便计算简单、清晰、离散。 首先说明如何利用凯莱-哈密顿定理直接计算. 矩阵指数函数eAt的定义: eAt=∞Σn=01n!Antn. 上式对所有的方阵A都有意义. 引理1 (凯莱-哈密顿定理)所有方阵A都满足它的特征方程.即如果 | (λt-A) |=bnλn+…+b1λ+b0, 则 bnAn+…+b1A+b0I=0, 定理1 若A是n阶方阵, 则 eAt=an-1An-1tn-1+…+a1At+a0I, 其中a0,a1, …,an-1是与A相对应的关于t的函数. 证明因为At是n阶方阵, 由引理1有 bnAntn+…+b1At+b0I=0, 所以 An+1tn+1=Antn·At= -bn-1bnAn-1tn-1+?+-b1bnAt+-b0bnΙ At=-bn-1bn -bn-1bnAn-1tn-1+?+-b1bnAt+-b0bnΙ +-bn-2bnAn-1tn-1+?+-b1bnA2t2+-b0bnAt=b2n-1-bnbn-2b2nAn-1tn-1+?+bn-1b1-bnb0b2nAt+bn-1b0b2nΙ, 以此类推, 当mn时,Amtm都可以表示为 An-1tn-1,An-2tn-2, …,At,I 的多项式. 定理2 若A是定理1中的矩阵, 定义 r(λ) =an-1λn-1+an-2λn-2+…+a1λ+a0, 如果λi是At的特征值 (A的特征值的t倍) , 则 eλi=r(λi) . 证明因为λi是At的特征值, 由引理1有 bnλni+…+b1λi+b0=0, 再由定理1即可证eλi=r(λi) . 进一步, 如果λi是k(k1) 重特征值, 那么下列等式成立: eλi=ddλr(λ)λ=λi, eλi=d2dλ2r(λ)λ=λi, ??????,eλi=dk-1dλk-1r(λ)|λ=λi, 该法适用于特征根的任何情形, 尤其当方程有不同的重特征根时, 利用该法计算基解矩阵特别简单. 例1 求eAt, 其中 A=(0100011-12). 解此处n=3, 有 eAt=a2A2t2+a1At+a0Ι=(a0a1ta2t20-a2t2+a02a2t2+a1t0-2a2t2-a1t3a2t2+2a1t+a0),r(λ)=a2λ2+a1λ+a0,r′(λ)=2a2λ+a1, 而At的特征值为 λ1=0,λ2=λ3=t. 于是有 et=r(t) , et=r′ (t) , e0=r(0) , 即得方程组 {et=a2t2+a1t+a0,et=2a2t+a1,e0=a0, 解之得 a0=e0=1,a1=-tet+2et-2t,a2=tet-et+1t2, 代入eAt的表达式得: eAt=(1-tet+2et-2tet-et+10-tet+ettet0-tettet+et). 以下说明如何利用循环公式求方程组的通解. 先来推导循环公式.对于方程 dxdt=Ax, 一般地,A必相似于一个约当形矩阵, 即必存在可逆矩阵T, 使得 T-1AT=diag (J1,J2, …,Js) , 其中 Jk=(λk1λk1?1λk)mk×mkk=1,2,?,s,sΣk=1mk=n. 令 x=Ty, 代入方程 dxdt=Ax, 得 dydt=Τ-1AΤy, 记 y=(?y1?y1??ys),?yk=(yk1yk2?ykmk)(k=1,2,?,s). 则 ddt(?y1?y1??ys)=(J1J2?Js)(?y1?y1??ys), 于是 ddt(yk1yk2?ykmk)=(λk1λk1?1λk)(yk1yk2?ykmk) 即对每一个k=1, 2, …,s都有 dykjdt=λkykj+yk(j+1)(j=1,2,?,mk-1),dykmkdt=λkykmk. 由此解出 ykmk=Ckmkeλkt,yk(mk-1)=(Ck(mk-1)+Ckmkt)eλkt,??????,yk1=(Ck1+Ck2t+?+Ckmktmk-1(mk-1)!)eλkt, 记 T= (H1,H2, …,Hs) , Hk= (hk1,hk2, …,hkmk) , 于是 x=Τy=(Η1,Η2,?,Ηs)(?y1?y1??ys)=sΣk=1Ηk?yk=sΣk=1(hk1,hk2,?,hkmk)(yk1yk2?ykmk)=sΣk=1(hk1yk1+hk2yk2+?+hkmkykmk)=sΣk=1[Ck1hk1+Ck2(thk1+hk2)+?+Ckmk(tmk-1(mk-1)!hk1+?+hkm