常微分方程考研讲义第五章 线性微分方程组(7)

2019-03-11 08:06

根据expAt的定义,一般来说,(expAt)?的分量是一个无穷级数,因而难于计算。这里的要点就是将初始向量?进行分解,从而使得(expAt)?的分量可以表示为t的指数函数与t的幂函数乘积的有限项的线性组合。

假设?1,?2,?,?k分别是矩阵A的n1,n2,?,nk重不同特征值。这时由(5.49),我们有

??v1?v???vn (5.50)

其中vj?Uj,(j?1,2,?k),因为子空间Uj是由方程组(5.48)产生的,vj一定是(5.48)的解。由此即得

(A??jE)lvj?0,l?nj,(j?1,2,?k)

(5.51)

注意到当矩阵是对角形时,由例1知道,expAt是很容易求得的,这时得到

?e?jt?0???t?t?ejexp(??jEt)?ej??????E

?0?e?jt???由此,并根据等式(5.51),既有

(expAt)vj?(expAt)e??exp(A??jE)t??vjn?1 t2tj?jt?nj?1?2?e?E?t(A??jE)?(A??jE)???(A??jE)?vj2!(nj?1)!?????jt??exp(??jEt)??vj?e?jt在根据等式(5.50),知微分方程组(5.33)的解?(t)?(expAt)?可表为

?(t)?(expAt)??(expAt)?vj??(expAt)vjj?1j?1kk??ej?1k?jt?ttnj?1?2E?t(A??E)?(A??E)???(A??E)??vjjjj2!(nj?1)!????2nj?1

所以,方程(5.33)满足?(0)??的解?(t)最后可以写成

?(t)??ej?1k?jt?nj?1ti?i??(A??jE)?vj (5.52) ?i?1i!?在特别情形,当A只有一个特征值时,无需将初始向量分解为(5.50),这时对于任何u,都有

(A??E)nu?0

这就是说,(A??E)n是一个零矩阵,这样一来,由expAt的定义,我们得到

expAt?e(A??E)t?e?t?tt(A??E)i (5.53) ?i?0i!n?1i为了要从(5.52)中得到expAt,只要注意到

expAt?(expAt)E??(expAt)e1,(expAt)e2,?,(expAt)en?

?1??0??0??0??1??0???????其中 e1????,e2??0?,?,en????

??????0??????0?????0???0???1??是单位向量,这就是说,依次令??e1,??e2,?,??en,求得n个解,以这n个解作为即可得到expAt。

例7 如果A是例4的矩阵,试解初值问题x??Ax,?(0)??,并求expAt。 解 从例4知道,?1?3是A的二重特征值,这时n1?2,只有一个子空间U1,将

??1?n1?2及????代入(5.52)即得

??2????11????1?3t??1?t(??1??1)? ?(t)?e?E?t(A?3E)???e?E?t?????e?????11????2???2?t(??1??1)??3t3t(5.54)

利用公式(5.53),即得

?t???11??3t?1?t expAt?e?E?t(A?3E)??e?E?t???e?????11????t1?t??3t3t或者,分别令

??e1???,??e2???

?0??1?然后代入(5.54),亦同样得到上面的结果

?1??0?t??1?t expAt?e3t????t1?t???41000??0?4100???0?400?,试求expAt 例8 如果A??0??000?40????0000?4??解 这里n?5,???4是A的5重特征值,直接计算可得(A?4E)3?0。因此,由公式(5.53)可得

expAt?e这样一来

?4t?t22? E?t(A?4E)?(A?4E)??2!????1??0???expAt?e?4t??0??0????0????1??0?0??0??0t10000000??0?01000???0100??t?0??0010??00001????01000??0?00100?2?t?0000???0?2!?0000??0?0000???00100???0000????0000???0000????0000???

?e?4t?t200?2!?t00?100??010??001???3x1?x2?x3?x1?3?11??????2x1?x3例9 考虑方程组?x2,这里系数矩阵A?201,试求满足初始条

???x??x?x?2x??1?12??123?3??1???件?(0)??2??的解?(t),并求expAt。 ?????3??解 A的特征方程为

det(?E?A)?(??1)(??2)2?0

?1?1,?2?2分别为n1?1,n2?2重特征值,为了确定三维欧几里得空间的子空间

U1和U2,根据(5.48),需要考虑下面方程组:

(A?E)u?0和(A?2E)2u?0

首先讨论

?2?11??u?0

(A?E)u??2?11????1?11??或

?2u1?u2?u3?0??2u1?u2?u3?0 ?u?u?u?0?123这个方程组的解为

?0??

u1??????????其中?为任意常数。子空间U1是由向量u1所张成的。其次讨论

?000??u?0

(A?2E)2u???110?????110??或

??u1?u2?0 ??u?u?0?12这个方程组的解为

????

u2??????????其中?,?是任意常数。子空间U2是由向量u2所张成的。

现在需要找出向量v1?U1,v2?U2使得能够将初始向量?写成(5.50)的形式。因为v1?U1,v2?U2,所以

?0?????,v????

v1???2??????????????其中?,?,?是某些常数,这样一来

??1??0??????????????? ?2????????3????????????因而???1,?????2,?????3,解之得到???2??1,???1,

???3??2??1,且

?1?0????,v???

v1??????21?21???????3??1????3??2??1??根据公式(5.52),我们得到满足初始条件?(0)??的解为

?(t)?etEv1?e2t(E?t(A?2E))v2??1?0??1?11??????e2t?E?t?2?21?????et?????211?????????????3??1???1?10?????3??2??1??????tt???1?0??1?t?

??e2t?2t1?2tt????et?????211?????????t1???3??1???t????3??2??1???0???1?t(?3??2??1)???e2t???t(?????)??et????321??21??1????3??1????3??2??1??为了得到expAt,依次令?等于

?1??0??0??0?,?1?,?0? ????????0????0????1??代入上式,我们得到三个线性无关的解。利用这三个解作为列,即得

?(1?t)e2t?expAt???et?(1?t)e2t??et?e2t??te2tet?te2tet?e2tte2t??te2t? e2t??应该指出,公式(5.52)是本节的主要结果。公式(5.52)告诉我们,常系数线性微分方程组(5.33)的任一解都可以通过有限次代数运算求出来。在常微分方程的


常微分方程考研讲义第五章 线性微分方程组(7).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:项目管理试题

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: