一阶非齐次线性微分方程组的常数变易公式

内容来源

常微分方程(第四版) (王高雄,周之铭,朱思铭,王寿松) 高等教育出版社


  • 一阶非齐次线性微分方程组

x′=A(t)x+f(t)(1) x'=A(t)x+f(t)\tag{1} x=A(t)x+f(t)(1)

其中 A(t)A(t)A(t) 是区间 [a,b][a,b][a,b] 上连续的 n×nn\times nn×n 矩阵,f(t)f(t)f(t) 是同一区间上连续的 nnn 维向量

初值条件 x(t0)=ηx(t_0)=\etax(t0)=η

  • 对应的齐次线性微分方程组

x′=A(t)x(2) x'=A(t)x\tag{2} x=A(t)x(2)

Φ(t)\Phi(t)Φ(t) 是其基解矩阵

(2)(2)(2) 的通解为

φ(t)=Φ(t)c \varphi(t)=\Phi(t)c φ(t)=Φ(t)c

ccc 变易为 ttt 的向量函数

φ(t)=Φ(t)c(t) \varphi(t)=\Phi(t)c(t) φ(t)=Φ(t)c(t)

代入 (1)(1)(1)

Φ′(t)c(t)+Φ(t)c′(t)=A(t)Φ(t)c(t)+f(t) \Phi'(t)c(t)+\Phi(t)c'(t)=A(t)\Phi(t)c(t)+f(t) Φ(t)c(t)+Φ(t)c(t)=A(t)Φ(t)c(t)+f(t)

因为 Φ(t)\Phi(t)Φ(t)(2)(2)(2) 的基解矩阵,所以 Φ′(t)=A(t)Φ(t)\Phi'(t)=A(t)\Phi(t)Φ(t)=A(t)Φ(t),那么

Φ(t)c′(t)=f(t) \Phi(t)c'(t)=f(t) Φ(t)c(t)=f(t)

移项积分得

c(t)=∫t0tΦ−1(s)f(s)ds c(t)=\int^t_{t_0}\Phi^{-1}(s)f(s)\mathrm{d}s c(t)=t0tΦ1(s)f(s)ds

这样 c(t0)=0c(t_0)=0c(t0)=0

所以 (1)(1)(1) 的解为

φ(t)=Φ(t)∫t0tΦ−1(s)f(s)ds \varphi(t)=\Phi(t)\int^t_{t_0}\Phi^{-1}(s)f(s)\mathrm{d}s φ(t)=Φ(t)t0tΦ1(s)f(s)ds

且满足初值条件 φ(t0)=0\varphi(t_0)=0φ(t0)=0

所以满足初值条件 φ(t0)=η\varphi(t_0)=\etaφ(t0)=η 的解为

φ(t)=Φ(t)Φ−1(t0)η+Φ(t)∫t0tΦ−1(s)f(s)ds \varphi(t)=\Phi(t)\Phi^{-1}(t_0)\eta+\Phi(t)\int^t_{t_0}\Phi^{-1}(s)f(s)\mathrm{d}s φ(t)=Φ(t)Φ1(t0)η+Φ(t)t0tΦ1(s)f(s)ds

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值