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

2019-03-11 08:06

的基解矩阵。

解首先,我们证明?(t)是解矩阵。令?1(t)表示?(t)的第一列,这时

?et??11??et??11??1?(t)??????????01??1(t)

01??0????0??这表示?1(t)是一个解。同样,如果以?2(t)表示?(t)的第二列,我们有

?(t?1)et??11??tet??11??(t)???2?2(t) ?t????????t?e??01??e??01?这表示?2(t)也是一个解。因此,?(t)???1(t),?2(t)?是解矩阵。

其次,根据定理2*,因为det?(t)?e2t?0,所以?(t)是基解矩阵。

推论1* 如果?(t)是(5.15)在区间a?t?b上的基解矩阵,C是非奇异n?n常数矩阵,那么,?(t)C也是(5.15)在区间a?t?b上的基解矩阵。

证明 首先,根据解矩阵的定义易知,方程(5.15)的任一解矩阵X(t)必满足关系

X?(t)?A(t)X(t),(a?t?b)

反之亦然。现令

?(t)??(t)C,(a?t?b)

微分上式,并注意到?(t)为方程的基解矩阵,C为常数矩阵,得到

??(t)???(t)C?A(t)?(t)C?A(t)?(t)

即?(t)是(5.15)的解矩阵。又由C的非奇异性,我们有

det?(t)?det?(t)?detC?0(a?t?b)

因此由定理2*知,?(t)即?(t)C是(5.15)的基解矩阵。

推论2* 如果?(t),?(t)在区间a?t?b上是x??A(t)x的两个基解矩阵,那么,存在一个非奇异n?n常数矩阵C,使得在区间a?t?b上?(t)??(t)C。 证明 因为?(t)为基解矩阵,故其逆矩阵??1(t)一定存在。现令

??1(t)?(t)?X(t) (a?t?b)

?(t)??(t)X(t) (a?t?b)

易知X(t)是n?n可微矩阵,且

detX(t)?0 (a?t?b)

于是

A(t)?(t)???(t)???(t)X(t)??(t)X?(t) (a?t?b)

???A(t)?(t)X(t)??(t)X(t)?A(t)?(t)??(t)X(t)由此推知?(t)X?(t)?0,或X?(t)?0(a?t?b),即X(t)为常数矩阵,记为C。因此我们有

?(t)??(t)C (a?t?b)

其中C???1(a)?(a)为非奇异的n?n常数矩阵推论2*得证。 5.2.2 非齐线性微分方程组 本段讨论非齐线性微分方程组

x??A(t)x?f(t) (5.14)

的解的结构问题,这里A(t)是区间a?t?b上的已知n?n连续矩阵,f(t)是区间

a?t?b上的已知n维连续列向量,向量f(t)通常称为强迫项,因为如果(5.14)描

述一个力学系统,f(t)就代表外力。 容易验证(5.14)的两个简单性质:

性质1 如果?(t)是(5.14)的解,?(t)是(5.14)对应的齐线性方程组(5.15)的解,则?(t)??(t)是(5.14)的解。

?(t)和?(t)是(5.14)的两个解,则??(t)??(t)是(5.15)的解。 性质2 如果? 下面的定理7给出(5.14)的解的结构。

定理7 设?(t)是(5.15)的基解矩阵,?(t)是(5.14)的某一解,则(5.14)的任一解?(t)都可表为

?(t)??(t)c??(t) (5.23)

这里c是确定的常数列向量。

?(t)??(t)是(5.15)的解,再由5.2.1的定理1*,得到 证明 由性质2我们知道??(t)??(t)??(t)c

这里c是确定的常数列向量,由此即得

?(t)??(t)c??(t)

定理证毕。

定理7告诉我们,为了寻求(5.15)的任一解,只要知道(5.14)的一个解和它对应的齐线性方程组(5.15)的基解矩阵。在知道(5.15)的基解矩阵?(t)的情况下,寻求(5.14)的解?(t)的简单的方法 常数变易法。

由定理1*可知,如果c是常数列向量,则?(t)??(t)c是(5.15)的解,它不可能是(5.14)的解。因此,将c变易为t的向量函数,而试图寻求(5.14)的形如

?(t)??(t)c(t) (5.24)

的解。这里c(t)是待定的向量函数。

假设(5.14)存在形如(5.24)的解,这时,将(5.24)代入(5.14)得到

??(t)c(t)??(t)c?(t)?A(t)?(t)c(t)?f(t)

因为?(t)是(5.15)的基解矩阵,所以??(t)?A(t)?(t),由此上式中含有

A(t)?(t)c(t)的项消去了。因而c(t)必须满足关系式

?(t)c?(t)?f(t) (5.25)

因为在区间a?t?b上?(t)是非奇异的,所以??1(t)存在。用??1(t)左乘(5.25)两边,得到

c(t)????1(s)f(s)ds,t0,t??a,b?

tt0其中c(t0)?0。这样,(5.24)变为

?(t)??(t)???1(s)f(s)ds,t0,t??a,b? (5.26)

tt0因此,如果(5.14)有一个形如(5.24)的解?(t),则?(t)由公式(5.26)决定。

反之,用公式(5.26)决定的向量函数?(t)必定是(5.14)的解。事实上,微分(5.26)得到

??(t)???(t)???1(s)f(s)ds??(t)??1(t)f(t)t0t?A(t)?(t)??(s)f(s)ds?f(t)t0t

?1再利用公式(5.26),即得

??(t)?A(t)?(t)?f(t)

显然,还有?(t0)?0,这样一来,我们就得到了下面的定理8。 定理8 如果?(t)是(5.15)的基解矩阵,则向量函数

?(t)??(t)???1(s)f(s)ds

t0t是(5.14)的解,且满足初始条件

?(t0)?0

由定理7和定理8容易看出(5.14)的满足初始条件

?(t0)??

的解?(t)由下面公式给出

?(t)??(t)?(t0)???(t)???1(s)f(s)ds (5.27)

?1t0t这里?h(t)??(t)??1(t0)?是(5.15)的满足初始条件

?h(t0)??

的解。公式(5.26)或公式(5.27)称为非齐线性微分方程组(5.14)的常数变易公式。 第五章

?e?t??x1??11???1?x?x?例2 x??? x(0)?????????01??x2??0??1??et解 在例1中我们已经知道?(t)???0tet?? et?是对应的齐线性方程组的基解矩阵。取矩阵?(t)的逆,我们得到:

?es?0??1?(t)?ses??es??1?s??s??e ?2se?01?这样,由定理8,满足初始条件

?0??(0)???

?0?的解就是

?et?(t)???0?e???0t?ettet?t?s?1?s??e?s?e????ds??t??001e????0??0?1??1?tet??(1?e?2t)??(et?e?t)??22t?????e?00????tet?t?e?2s??dst??0?e??0?

因为?(0)?E,对应的齐线性方程组满足初始条件

?h(0)???

1??的解就是

??1???1??(t?1)et??h(t)??(t)????t?

?1??e?由公式(5.27),所求解就是

1?1????(t?1)et??(et?e?t)??tet?(et?e?t)??(t)??h(t)??(t)??t??2? 2????e??0et????注意到5.1.1关于n阶线性微分方程的初值问题(5.6)与线性微分方程组的初值问题(5.7)等价性的讨论,我们可以得到关于n阶非齐线性微分方程的常数变易公式。 推论3 如果a1(t),a2(t),?,an(t),f(t)是区间a?t?b上的连续函数,

x1(t),x2(t),?,xn(t)是区间a?t?b上齐线性方程

x(n)?a1(t)x(n?1)???an(t)x?0 (5.21)

的基本解组,那么,非齐线性方程

x(n)?a1(t)x(n?1)???an(t)x?f(t) (5.28)


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

下一篇:项目管理试题

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

马上注册会员

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