matlab解微分方程_常微分方程:(第五章)线性微分方程组

参考《常微分方程》第三版(王高雄)

常微分方程 王高雄 第五章 线性微分方程组_哔哩哔哩 (゜-゜)つロ 干杯~-bilibili​www.bilibili.com
a8742e2dbcc1bdbd6b5b2ca877ff5ef1.png

前面介绍了线性微分方程的求解,现在考虑线性微分方程组的求解。

5.1 存在唯一性定理

记号和定义:

考虑如下方程组:(p186)

9874446c3363fb54abef2808e59c984b.png

作以下变换:

6f5ec4b4901f04f71906844ebce8349f.png

则(5.1)等价为以下形式:

bf0166902b390c20745c7a5aed000991.png

齐次线性与非齐次线性微分方程组:(p202)

fcb053dd719e920f31086f166c69593d.png

参考第四章,我们也先引入一些用到的定义。

一.常用定义:

1.朗斯基行列式(p203)

08f04b19e17e1e074c6ad34f8d735360.png

2.向量函数的线性相关&线性无关(p202)

6b74152137b1c95647512d45851c550f.png

a3f6752c2e8d9047aad39e500896cac6.png

3.基本解组(p205)

cdcb91e6dc6bbef378950366184e356a.png

4.矩阵或向量在闭区间上连续、可微、可积:如果它的每一个元都在该闭区间上连续、可微、可积

5.(5.4)的解的定义(p189)

633f0476544d4f8cd6126d6c4ff7584e.png

6.矩阵可微的性质(p188)

38899971c8b612897c0c80014284f680.png

7.矩阵和向量的范数定义及其性质(p195)

c59d8531f3b9848e8db80af0ff523de6.png

8.解矩阵、基解矩阵(p208)

f6b4c2aebe1f1aa32b28794c899e8a01.png

9.基解矩阵的性质(p208)

ebff35ccb14a627f369434f74275c285.png

3d17f47c2a83bfeab51f1194c76c0d3b.png

1417f60517f79a0e803adbb681aa92e3.png

dd116efef992cf3af7c00cdb22cb0bf2.png

二.齐次线性微分方程组基本性质

a)存在唯一性定理(p196)

2e9cdc17f9f94ea11fc09fe5acf91aeb.png

b)叠加原理(p202)

f23b70fd2ba64b4e3644ff05275e2bc7.png

c) 向量函数线性相关性与朗斯基行列式的关系(p203)

4f85c984e9793b11d9ec3acfbc17c6d5.png

d)方程解(函数)线性无关,则朗斯基行列式不为0(p204)

5982b5941676ab67fb919778411ebfdb.png

e)通解结构定理(p205)

e36016e74af0fec266be7ecea98983f0.png

6e6a013466f937c06efff4673ce1d912.png

三.非齐次线性微分方程组的基本性质及常数变易法

a)齐次线性微分方程组与非齐次线性微分方程组解之间的关系(p211)

c6aba17d0a9354baef68f3cc7d1d3d8e.png

b)常数变易法(p211)

37e242e0db8287c6af4073b154b95a6f.png

定理

7bd5996c0eb8438b12585e90903a8a3b.png

此外:注意到5.1.1关于n阶线性微分方程的初值问题(5.6)与线性微分方程组的初值问题(5.7)等价性的论述,我们可以得到关于n阶非齐次线性微分方程的常数变易公式.

ae7986c74504514ec412c8b8260aaaf0.png

158596fae24b28cf60cfe9d28fbc4d73.png
  • 综上,对于n阶非齐次线性微分方程(5.28),其满足初值条件的解由下面公式给出:

537d288cbdd24cc89590a6b7675e6001.png

(5.28)的任一解

都具有形式:

496c6c35876909866429010f66f835b9.png
  • 特别的,当n=2时,公式(5.29)就是:

9689bbaf2060658e7a1cf930288d1856.png

通解就是

79fc33fa2b65699f26049c2652356bbd.png

5.3 常系数线性微分方程(p219)

,A是n*n的常数矩阵

5.3.1 矩阵指数expA的定义和性质

定义:

c9a0dea873649242f9a84c2eb2670a3e.png

我们规定

其中

29e9e6186ddf95782ab7632c95245a69.png

expA的性质

9770d2fbe91b9b9a108b4f8a6ddcf155.png

48fae7320646f8b84627c6669af3d8d7.png

29ce7c5956984f636a47073aebe747c6.png

常系数齐次线性微分方程组(5.33)的解:

75a544302be0d47fdb44c234e66f43fc.png

c6f4486d4a7ee86ec35c12a0c28023e9.png

5.3.2 基解矩阵的计算公式

a)矩阵的特征值、特征向量、特征多项式、特征方程

dbb31a35121f2e7c03f19f7756e025eb.png

b)常系数微分方程的基解矩阵定理(p227)

76c71743c4b7b0a27ee80712b0033d00.png

c)当A是任意的n*n矩阵时,(5.33)的基解矩阵的计算方法(p228)

一般地,有

bfaeda22e9a27951678d86bbb46bad77.png

9e60556254ee21cc005fd72aa9e02b15.png

其中

d4ac74c2315f31f41ad8aad1068e4584.png

(5.52)的应用

9e78123c3208ee695d806ee8aa757114.png

matlab解微分方程

y
  • p213例2:

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

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

  • 3
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值