微分方程-非齐次线性方程组的通解

非齐次线性方程组的通解

本文根据之前的结果进一步给出非齐次线性方程组

\dfrac{\text{d}\boldsymbol{x}}{\text{d}t}=\boldsymbol{A}(t)\boldsymbol{x}+\boldsymbol{f}(t)\quad(3.15)

的解的结构. 我们假设 \boldsymbol{A}(t) 是区间 [\alpha,\,\beta] 上的 n\times n 阶连续矩阵函数,\boldsymbol{f}(t) 是区间 [\alpha,\,\beta] 上的 n 维连续列向量函数. 当 \boldsymbol{f}(t)\equiv\boldsymbol{0} ,方程组(3.15)变为(3.9). 我们称方程组(3.9)时方程组(3.15)相应的齐次线性方程组.

容易验证,方程组(3.15)的解与方程组(3.9)的解之间有如下关系:如果 \boldsymbol{x}^{*}_{1}(t)\boldsymbol{x}^{*}_{2}(t) 是方程组(3.15)的两个解,则 \boldsymbol{x}_{1}^{*}-\boldsymbol{x}_{2}^{*} 是相应的齐次线性方程组(3.9)的解;反之,如果 \boldsymbol{x}^{*}(t)\boldsymbol{x}(t) 分别是方程组(3.15)和方程组(3.9)的解,则 \boldsymbol{x}^{*}(t)+\boldsymbol{x}(t) 也是当成组(3.15)的一组解. 一般地,非齐次线性方程组的通解结构如下:


定理 3.4

\boldsymbol{x}^{*}(t) 是方程组(3.15)的某个解,则方程组(3.15)的任一解 \boldsymbol{x}(t) 都可以表示为

\boldsymbol{x}(t)=\boldsymbol{X}(t)\boldsymbol{c}+\boldsymbol{x}^{*}(t)\quad(3.16)

其中 \boldsymbol{c}n 维常数列向量,\boldsymbol{X}(t) 是相应齐次方程组(3.9)的基解矩阵.

证明

由方程组(3.15)的解与方程组(3.9)的解之间的关系知,\boldsymbol{x}(t)-\boldsymbol{x}^{*}(t) 是齐次方程组(3.9)的解. 由之前的只是,存在一个 n 维常数列向量 \boldsymbol{c},使得 \boldsymbol{x}(t)-\boldsymbol{x}^{*}(t)=\boldsymbol{X}(t)\boldsymbol{c},从而得到(3.16).


定理 3.4 表明,要找出方程组(3.15)的全部解,只需找出它的一个特解以及它相应的齐次线性方程组(3.9)的基解矩阵即可. 下面我们运用常数变易法来证明:只要找到齐次方程组(3.9)的基解矩阵 \boldsymbol{X}(t),就能切丁非齐次线性方程组(3.15)的一个特解,从而给出其通解.

这就是下面的定理:

定理 3.5

若矩阵函数 \boldsymbol{X}(t),\;\alpha\leqslant t\leqslant\beta 是齐次线性方程组(3.9)的基解矩阵,则非齐次线性方程组(3.15)的通解为

\displaystyle\boldsymbol{x}(t)=\boldsymbol{X}(t)\boldsymbol{c}+\boldsymbol{X}(t)\int_{t_0}^t\boldsymbol{X}^{-1}(\tau)\boldsymbol{f}(\tau)\text{d}\tau\quad(3.17)

其中 \boldsymbol{c}n 维常数列向量. 方程组(3.15)满足初值条件 \boldsymbol{x}(t_0)=\boldsymbol{x}^0 的解为

\displaystyle\boldsymbol{x}(t)=\boldsymbol{X}(t,t_0)\boldsymbol{x}^0+\int_{t_0}^t\boldsymbol{X}(t,\tau)\boldsymbol{f}(\tau)\text{d}\tau\quad(3.18)

其中 t,t_0\in[\alpha,\,\beta],\boldsymbol{X}(t,t_0)=\boldsymbol{X}(t)\boldsymbol{X}^{-1}(t_0)

由于我们在证明中用到的方法i奥做常数变易法,因此,我们把(3.17)或(3.18)称为常数变易公式.

证明

由之前知识,齐次方程组(3.9)的通解为 \boldsymbol{X}(t)\boldsymbol{c},其中 \boldsymbol{c} 为任意 n 维常数列向量. 现在我们把常数向量 \boldsymbol{c} 变易为 t 的待定列向量函数 \boldsymbol{c}(t),以期寻找非齐次线性方程组(3.15)的形如

\boldsymbol{x}^{*}(t)=\boldsymbol{X}(t)\boldsymbol{c}(t)\quad(3.19)

的特解. 把(3.19)带入非齐次线性方程组(3.15),得到

\displaystyle\dfrac{\text{d}\boldsymbol{X}(t)}{\text{d}t}\boldsymbol{c}(t)+\boldsymbol{X}(t)\dfrac{\text{d}\boldsymbol{c}(t)}{\text{d}t}=\boldsymbol{A}(t)\boldsymbol{X}(t)\boldsymbol{c}(t)+\boldsymbol{f}(t)

由于 \boldsymbol{X}(t) 是齐次方程组(3.9)的解矩阵,所有

\dfrac{\text{d}\boldsymbol{X}(t)}{\text{d}t}=\boldsymbol{A}(t)\boldsymbol{X}(t).

因此向量函数 \boldsymbol{c}(t) 满足微分方程组

\boldsymbol{X}(t)\dfrac{\text{d}\boldsymbol{c}(t)}{\text{d}t}=\boldsymbol{f}(t).\quad(3.20)

因为 \boldsymbol{X}(t) 是方程组(3.9)的基解矩阵,所以它的逆 \boldsymbol{X}^{-1}(t) 存在. 从而由(3.20)得

\dfrac{\text{d}\boldsymbol{c}(t)}{\text{d}t}=\boldsymbol{X}^{-1}\boldsymbol{f}(t)

对上式两边积分得到

\displaystyle\boldsymbol{c}(t)=\boldsymbol{c}^0+\int_{t_0}^{t}\boldsymbol{X}^{-1}(\tau)\boldsymbol{f}(\tau)\text{d}\tau

其中 \boldsymbol{c}^0 为常向量. 特别地,取 \boldsymbol{c}^0=\boldsymbol{0},借的到方程组(3.15)得一个特解

\boldsymbol{x}^{*}(t)=\boldsymbol{X}(t)\int_{t_0}^{t}\boldsymbol{X}^{-1}(\tau)\boldsymbol{f}(\tau)\text{d}\tau

再由定理 3.4 即得(3.17)并相应得到(3.18). 定理得证.

尽管我们拥有了验证齐次线性方程组(3.9)得一个解矩阵是否为基解矩阵的方法,然而,要计算一个齐次线性微分方程组的基解矩阵仍然不是一件容易的事. 对常系数矩阵 \boldsymbol{A}(t)\equiv\boldsymbol{A} 的情形,之后会给出方法.

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,657评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,662评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,143评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,732评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,837评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,036评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,126评论 3 410
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,868评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,315评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,641评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,773评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,470评论 4 333
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,126评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,859评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,095评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,584评论 2 362
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,676评论 2 351

推荐阅读更多精彩内容