参考《常微分方程》第三版(王高雄)
常微分方程 王高雄 第五章 线性微分方程组_哔哩哔哩 (゜-゜)つロ 干杯~-bilibiliwww.bilibili.com![a8742e2dbcc1bdbd6b5b2ca877ff5ef1.png](https://img-blog.csdnimg.cn/img_convert/a8742e2dbcc1bdbd6b5b2ca877ff5ef1.png)
前面介绍了线性微分方程的求解,现在考虑线性微分方程组的求解。
5.1 存在唯一性定理
记号和定义:
考虑如下方程组:(p186)
![9874446c3363fb54abef2808e59c984b.png](https://img-blog.csdnimg.cn/img_convert/9874446c3363fb54abef2808e59c984b.png)
作以下变换:
![6f5ec4b4901f04f71906844ebce8349f.png](https://img-blog.csdnimg.cn/img_convert/6f5ec4b4901f04f71906844ebce8349f.png)
则(5.1)等价为以下形式:
![bf0166902b390c20745c7a5aed000991.png](https://img-blog.csdnimg.cn/img_convert/bf0166902b390c20745c7a5aed000991.png)
齐次线性与非齐次线性微分方程组:(p202)
![fcb053dd719e920f31086f166c69593d.png](https://img-blog.csdnimg.cn/img_convert/fcb053dd719e920f31086f166c69593d.png)
参考第四章,我们也先引入一些用到的定义。
一.常用定义:
1.朗斯基行列式(p203)
![08f04b19e17e1e074c6ad34f8d735360.png](https://img-blog.csdnimg.cn/img_convert/08f04b19e17e1e074c6ad34f8d735360.png)
2.向量函数的线性相关&线性无关(p202)
![6b74152137b1c95647512d45851c550f.png](https://img-blog.csdnimg.cn/img_convert/6b74152137b1c95647512d45851c550f.png)
![a3f6752c2e8d9047aad39e500896cac6.png](https://img-blog.csdnimg.cn/img_convert/a3f6752c2e8d9047aad39e500896cac6.png)
3.基本解组(p205)
![cdcb91e6dc6bbef378950366184e356a.png](https://img-blog.csdnimg.cn/img_convert/cdcb91e6dc6bbef378950366184e356a.png)
4.矩阵或向量在闭区间上连续、可微、可积:如果它的每一个元都在该闭区间上连续、可微、可积
5.(5.4)的解的定义(p189)
![633f0476544d4f8cd6126d6c4ff7584e.png](https://img-blog.csdnimg.cn/img_convert/633f0476544d4f8cd6126d6c4ff7584e.png)
6.矩阵可微的性质(p188)
![38899971c8b612897c0c80014284f680.png](https://img-blog.csdnimg.cn/img_convert/38899971c8b612897c0c80014284f680.png)
7.矩阵和向量的范数定义及其性质(p195)
![c59d8531f3b9848e8db80af0ff523de6.png](https://img-blog.csdnimg.cn/img_convert/c59d8531f3b9848e8db80af0ff523de6.png)
8.解矩阵、基解矩阵(p208)
![f6b4c2aebe1f1aa32b28794c899e8a01.png](https://img-blog.csdnimg.cn/img_convert/f6b4c2aebe1f1aa32b28794c899e8a01.png)
9.基解矩阵的性质(p208)
![ebff35ccb14a627f369434f74275c285.png](https://img-blog.csdnimg.cn/img_convert/ebff35ccb14a627f369434f74275c285.png)
![3d17f47c2a83bfeab51f1194c76c0d3b.png](https://img-blog.csdnimg.cn/img_convert/3d17f47c2a83bfeab51f1194c76c0d3b.png)
![1417f60517f79a0e803adbb681aa92e3.png](https://img-blog.csdnimg.cn/img_convert/1417f60517f79a0e803adbb681aa92e3.png)
![dd116efef992cf3af7c00cdb22cb0bf2.png](https://img-blog.csdnimg.cn/img_convert/dd116efef992cf3af7c00cdb22cb0bf2.png)
二.齐次线性微分方程组基本性质
a)存在唯一性定理(p196)
![2e9cdc17f9f94ea11fc09fe5acf91aeb.png](https://img-blog.csdnimg.cn/img_convert/2e9cdc17f9f94ea11fc09fe5acf91aeb.png)
b)叠加原理(p202)
![f23b70fd2ba64b4e3644ff05275e2bc7.png](https://img-blog.csdnimg.cn/img_convert/f23b70fd2ba64b4e3644ff05275e2bc7.png)
c) 向量函数线性相关性与朗斯基行列式的关系(p203)
![4f85c984e9793b11d9ec3acfbc17c6d5.png](https://img-blog.csdnimg.cn/img_convert/4f85c984e9793b11d9ec3acfbc17c6d5.png)
d)方程解(函数)线性无关,则朗斯基行列式不为0(p204)
![5982b5941676ab67fb919778411ebfdb.png](https://img-blog.csdnimg.cn/img_convert/5982b5941676ab67fb919778411ebfdb.png)
e)通解结构定理(p205)
![e36016e74af0fec266be7ecea98983f0.png](https://img-blog.csdnimg.cn/img_convert/e36016e74af0fec266be7ecea98983f0.png)
![6e6a013466f937c06efff4673ce1d912.png](https://img-blog.csdnimg.cn/img_convert/6e6a013466f937c06efff4673ce1d912.png)
三.非齐次线性微分方程组的基本性质及常数变易法
a)齐次线性微分方程组与非齐次线性微分方程组解之间的关系(p211)
![c6aba17d0a9354baef68f3cc7d1d3d8e.png](https://img-blog.csdnimg.cn/img_convert/c6aba17d0a9354baef68f3cc7d1d3d8e.png)
b)常数变易法(p211)
![37e242e0db8287c6af4073b154b95a6f.png](https://img-blog.csdnimg.cn/img_convert/37e242e0db8287c6af4073b154b95a6f.png)
定理
![7bd5996c0eb8438b12585e90903a8a3b.png](https://img-blog.csdnimg.cn/img_convert/7bd5996c0eb8438b12585e90903a8a3b.png)
此外:注意到5.1.1关于n阶线性微分方程的初值问题(5.6)与线性微分方程组的初值问题(5.7)等价性的论述,我们可以得到关于n阶非齐次线性微分方程的常数变易公式.
![ae7986c74504514ec412c8b8260aaaf0.png](https://img-blog.csdnimg.cn/img_convert/ae7986c74504514ec412c8b8260aaaf0.png)
![158596fae24b28cf60cfe9d28fbc4d73.png](https://img-blog.csdnimg.cn/img_convert/158596fae24b28cf60cfe9d28fbc4d73.png)
- 综上,对于n阶非齐次线性微分方程(5.28),其满足初值条件的解由下面公式给出:
![537d288cbdd24cc89590a6b7675e6001.png](https://img-blog.csdnimg.cn/img_convert/537d288cbdd24cc89590a6b7675e6001.png)
(5.28)的任一解
![496c6c35876909866429010f66f835b9.png](https://img-blog.csdnimg.cn/img_convert/496c6c35876909866429010f66f835b9.png)
- 特别的,当n=2时,公式(5.29)就是:
![9689bbaf2060658e7a1cf930288d1856.png](https://img-blog.csdnimg.cn/img_convert/9689bbaf2060658e7a1cf930288d1856.png)
通解就是
![79fc33fa2b65699f26049c2652356bbd.png](https://img-blog.csdnimg.cn/img_convert/79fc33fa2b65699f26049c2652356bbd.png)
5.3 常系数线性微分方程(p219)
5.3.1 矩阵指数expA的定义和性质
定义:
![c9a0dea873649242f9a84c2eb2670a3e.png](https://img-blog.csdnimg.cn/img_convert/c9a0dea873649242f9a84c2eb2670a3e.png)
我们规定
其中
![29e9e6186ddf95782ab7632c95245a69.png](https://img-blog.csdnimg.cn/img_convert/29e9e6186ddf95782ab7632c95245a69.png)
expA的性质
![9770d2fbe91b9b9a108b4f8a6ddcf155.png](https://img-blog.csdnimg.cn/img_convert/9770d2fbe91b9b9a108b4f8a6ddcf155.png)
![48fae7320646f8b84627c6669af3d8d7.png](https://img-blog.csdnimg.cn/img_convert/48fae7320646f8b84627c6669af3d8d7.png)
![29ce7c5956984f636a47073aebe747c6.png](https://img-blog.csdnimg.cn/img_convert/29ce7c5956984f636a47073aebe747c6.png)
常系数齐次线性微分方程组(5.33)的解:
![75a544302be0d47fdb44c234e66f43fc.png](https://img-blog.csdnimg.cn/img_convert/75a544302be0d47fdb44c234e66f43fc.png)
![c6f4486d4a7ee86ec35c12a0c28023e9.png](https://img-blog.csdnimg.cn/img_convert/c6f4486d4a7ee86ec35c12a0c28023e9.png)
5.3.2 基解矩阵的计算公式
a)矩阵的特征值、特征向量、特征多项式、特征方程
![dbb31a35121f2e7c03f19f7756e025eb.png](https://img-blog.csdnimg.cn/img_convert/dbb31a35121f2e7c03f19f7756e025eb.png)
b)常系数微分方程的基解矩阵定理(p227)
![76c71743c4b7b0a27ee80712b0033d00.png](https://img-blog.csdnimg.cn/img_convert/76c71743c4b7b0a27ee80712b0033d00.png)
c)当A是任意的n*n矩阵时,(5.33)的基解矩阵的计算方法(p228)
一般地,有
![bfaeda22e9a27951678d86bbb46bad77.png](https://img-blog.csdnimg.cn/img_convert/bfaeda22e9a27951678d86bbb46bad77.png)
![9e60556254ee21cc005fd72aa9e02b15.png](https://img-blog.csdnimg.cn/img_convert/9e60556254ee21cc005fd72aa9e02b15.png)
其中
![d4ac74c2315f31f41ad8aad1068e4584.png](https://img-blog.csdnimg.cn/img_convert/d4ac74c2315f31f41ad8aad1068e4584.png)
(5.52)的应用
![9e78123c3208ee695d806ee8aa757114.png](https://img-blog.csdnimg.cn/img_convert/9e78123c3208ee695d806ee8aa757114.png)
matlab解微分方程:
y
- p213例2:
![99a1958265a43ae1f9b538866241c95a.png](https://img-blog.csdnimg.cn/img_convert/99a1958265a43ae1f9b538866241c95a.png)
matlab实现:
[x1,x2]=dsolve('Dx1=x1+x2+exp(-t)','Dx2=x2','x1(0)=-1','x2(0)=1')
x1 =
t*exp(t) - exp(t)*(exp(-2*t)/2 + 1/2)
x2 =
exp(t)
可见结构与书上结果相同。
- p215例3:
![dcc36842eb8aaf46196bcd0f34a1d23b.png](https://img-blog.csdnimg.cn/img_convert/dcc36842eb8aaf46196bcd0f34a1d23b.png)
matlab实现:
x=dsolve('D2x+x=tan(t)')
x =
C7*cos(t) - sin(2*t)/2 + C8*sin(t) - cos(t)*(2*atanh(tan(t/2)) - sin(t))
虽然与书中结果不一致,但是也是正确的,注意:tan(t),对于自变量t,要加括号。
2020.11.20