日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

第二十八讲 解非齐次线性方程组

發(fā)布時間:2025/4/16 编程问答 19 豆豆
生活随笔 收集整理的這篇文章主要介紹了 第二十八讲 解非齐次线性方程组 小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.

一,關(guān)于二階方程組x?′=Ax?{\vec{x}}'=A\vec{x}x=Ax的理論(對n階方程也成立):
(假設(shè)A是常數(shù)矩陣)

  • 定理Ax?′=Ax?{\vec{x}}'=A\vec{x}x=Ax的通解是x?=c1x1?+c2x2?\vec{x}=c_{1}\vec{x_{1}}+c_{2}\vec{x_{2}}x=c1?x1??+c2?x2??x1?\vec{x_{1}}x1??x2?\vec{x_{2}}x2??線性無關(guān))
  • 證明可以用線性疊加原理,這里不做詳細(xì)說明了。
  • 定理B:朗斯基行列式W(x1?,x2?):=∣x1?x2?∣W(\vec{x_{1}},\vec{x_{2}}):=|\vec{x_{1}} \vec{x_{2}}|W(x1??,x2??):=x1??x2??,(符號:=:=:=表示定義或等價),這里x1?\vec{x_{1}}x1??x2?\vec{x_{2}}x2??不一定線性無關(guān)。如果二階方程有兩個解,那么朗斯基行列式是自變量t的函數(shù),并且有兩種可能性,要么W(t)?≡0W(t)\not\equiv 0W(t)??0(當(dāng)x1?\vec{x_{1}}x1??x2?\vec{x_{2}}x2??線性無關(guān)時),要么W≡0W\equiv 0W0(當(dāng)x1?\vec{x_{1}}x1??x2?\vec{x_{2}}x2??線性相關(guān)時)。
  • 方程組的基本矩陣x?′=Ax?{\vec{x}}'=A\vec{x}x=Ax,特征向量矩陣X:=[x1?x2?]X:=\begin{bmatrix}\vec{x_{1}} & \vec{x_{2}}\end{bmatrix}X:=[x1???x2???]x1?\vec{x_{1}}x1??x2?\vec{x_{2}}x2??線性無關(guān)
  • X的性質(zhì)1:|X|對于任意自變量t都≠0,因為x1?\vec{x_{1}}x1??x2?\vec{x_{2}}x2??線性無關(guān)
  • X的性質(zhì)2X′=[x1?′x2?′]=[Ax1?Ax2?]=A[x1?x2?]=AX{X}'=\begin{bmatrix}{\vec{x_{1}}}' & {\vec{x_{2}}}'\end{bmatrix}=\begin{bmatrix}A\vec{x_{1}} & A\vec{x_{2}}\end{bmatrix}=A\begin{bmatrix}\vec{x_{1}} & \vec{x_{2}}\end{bmatrix}=AXX=[x1???x2???]=[Ax1???Ax2???]=A[x1???x2???]=AX
  • 二,解非齊次線性方程組:
    一般形式:{x′=ax+by+r1(t)y′=cx+dy+r2(t)\left\{\begin{matrix}{x}'=ax+by+{\color{Red} r_{1}(t)}\\ {y}'=cx+dy+{\color{Red} r_{2}(t)}\end{matrix}\right.{x=ax+by+r1?(t)y=cx+dy+r2?(t)?
    簡化形式:x?′=Ax?+r?(t){\vec{x}}'=A\vec{x}+\vec{r}(t)x=Ax+r(t)
    定理C:微分方程組的通解xg?=xc?+xp?\vec{x_{g}}=\vec{x_{c}}+\vec{x_{p}}xg??=xc??+xp??,其中xc?\vec{x_{c}}xc??x?′=Ax?{\vec{x}}'=A\vec{x}x=Ax的通解,xp?\vec{x_{p}}xp??是微分方程組的一個特解。可以用線性疊加原理證明。
    找到特解xp?\vec{x_{p}}xp??是求解的關(guān)鍵。

    三,例題:

    圖中,箭頭表示流向,數(shù)字表示流速,單位是L/h,x表示左邊容器中鹽的含量,y表示右邊容器中鹽的含量,兩個容器的容量都是1L。
    假設(shè)輸入項為:外部流入左邊容器的液體濃度是5e?t5e^{-t}5e?t,外部流入右邊容器的液體濃度是0。輸入項不全為0,決定了方程組是非齊次方程組。

  • 建立微分方程組:
    x′=?3x+2y+5e?t{x}'=-3x+2y+5e^{-t}x=?3x+2y+5e?t
    含義:左邊容器x的變化率=-流出速度X左容器中鹽的濃度+內(nèi)部流入速度X右容器中鹽的濃度+從外部流入速度X外部液體的濃度。
    y′=3x?4y+0{y}'=3x-4y+0y=3x?4y+0
    含義:右邊容器y的變化率=內(nèi)部流入速度X左容器中鹽的濃度+流出速度X右容器中鹽的濃度+從外部流入速度X外部液體的濃度。
    {x′=?3x+2y+5e?ty′=3x?4y+0\left\{\begin{matrix}{x}'=-3x+2y+5e^{-t}\\ {y}'=3x-4y+0\end{matrix}\right.{x=?3x+2y+5e?ty=3x?4y+0?

  • 矩陣化:
    [x′y′]=[?323?4][xy]+[5e?t0]\begin{bmatrix}{x}'\\ {y}'\end{bmatrix}=\begin{bmatrix}-3 & 2\\ 3 & -4\end{bmatrix}\begin{bmatrix}x\\ y\end{bmatrix}+\begin{bmatrix}5e^{-t}\\ 0\end{bmatrix}[xy?]=[?33?2?4?][xy?]+[5e?t0?]

  • 參數(shù)變分法求特解xp?\vec{x_{p}}xp??
    xp?=v1(t)x1?+v2(t)x2?\vec{x_{p}}=v_{1}(t)\vec{x_{1}}+v_{2}(t)\vec{x_{2}}xp??=v1?(t)x1??+v2?(t)x2??
    和定理A類似,只不過把常數(shù)c改成了參數(shù)v
    化為基本矩陣:xp?=[x1?x2?][v1(t)v2(t)]=Xv?\vec{x_{p}}=\begin{bmatrix} \vec{x_{1}}& \vec{x_{2}}\end{bmatrix}\begin{bmatrix}v_{1}(t)\\ v_{2}(t)\end{bmatrix}=X\vec{v}xp??=[x1???x2???][v1?(t)v2?(t)?]=Xv

  • 將特解xp?\vec{x_{p}}xp??代入方程組x?′=Ax?+r?(t){\vec{x}}'=A\vec{x}+\vec{r}(t)x=Ax+r(t),求出v?\vec{v}v
    代入得:xp?′=Axp?+r?(t){\vec{x_{p}}}'=A\vec{x_{p}}+\vec{r}(t)xp??=Axp??+r(t)
    等式左邊:xp?′=(Xv?)′=X′v?+Xv?′{\vec{x_{p}}}'={(X\vec{v})}'={X}'\vec{v}+X{\vec{v}}'xp??=(Xv)=Xv+Xv,(乘積的求導(dǎo)公式)
    等式右邊:Ax?+r?(t)=AXv?+r?(t)A\vec{x}+\vec{r}(t)=AX\vec{v}+\vec{r}(t)Ax+r(t)=AXv+r(t)
    因為X是方程組的基本矩陣,所以根據(jù)X的性質(zhì)2:AX=X′AX={X}'AX=X
    等式右邊:Ax?+r?(t)=AXv?+r?(t)=X′v?+r?(t)A\vec{x}+\vec{r}(t)=AX\vec{v}+\vec{r}(t)={X}'\vec{v}+\vec{r}(t)Ax+r(t)=AXv+r(t)=Xv+r(t)
    左邊=右邊:X′v?+Xv?′=X′v?+r?(t)?Xv?′=r?(t){X}'\vec{v}+X{\vec{v}}'={X}'\vec{v}+\vec{r}(t)\Rightarrow X{\vec{v}}'=\vec{r}(t)Xv+Xv=Xv+r(t)?Xv=r(t)
    v?′=X?1r?(t){\vec{v}}'=X^{-1}\vec{r}(t)v=X?1r(t),根據(jù)X的性質(zhì)1,X存在逆矩陣
    v?=∫X?1r?(t)dt\vec{v}=\int X^{-1}\vec{r}(t)dtv=X?1r(t)dtX?1r?(t)X^{-1}\vec{r}(t)X?1r(t)是一個列向量,元素都是t的函數(shù),只要逐個積分就算出來了。
    結(jié)果:xp?=Xv?=X∫X?1r?(t)dt\vec{x_{p}}=X\vec{v}=X\int X^{-1}\vec{r}(t)dtxp??=Xv=XX?1r(t)dt,只要找到一個特解就行,因此不用在公式后加積分常數(shù)。

  • x?′=Ax?{\vec{x}}'=A\vec{x}x=Ax的通解xc?\vec{x_{c}}xc??的部分省略了

  • 總結(jié)

    以上是生活随笔為你收集整理的第二十八讲 解非齐次线性方程组的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網(wǎng)站內(nèi)容還不錯,歡迎將生活随笔推薦給好友。

    主站蜘蛛池模板: 男生女生搞鸡视频 | 久久艹久久 | 精品伦精品一区二区三区视频密桃 | 高清av免费| 一本加勒比波多野结衣 | 成人黄色一区二区三区 | 亚洲最新在线观看 | 国产一区二区在线免费观看视频 | 实拍女处破www免费看 | 国产女人18毛片水18精品 | 黑人操亚洲人 | 91久久精品国产91性色tv | 欧美顶级毛片在线播放 | www九九九 | 91这里只有精品 | 中文字幕在线网站 | 热久久影院 | 一区二区三区视频在线免费观看 | 精品国产av一区二区三区 | 北条麻妃av在线播放 | 成人免费网站视频 | 伊人一二三 | 污污视频网站在线免费观看 | 国产1级片| 国产成人无码精品久久久久久 | 法国空姐在线观看视频 | 国产三级国产精品国产国在线观看 | 韩国黄色网址 | 干欧美少妇 | 麻豆视频二区 | 久久久久黄 | 日韩av资源站 | 成年人免费av | 亚洲h动漫 | 91福利在线播放 | 97狠狠操| 国产无遮挡又黄又爽在线观看 | 51热门大瓜今日大瓜 | 一区二区三区午夜 | 欧美在线一二 | 欧美在线观看免费高清 | 日本一区二区三区久久 | 五月激情综合婷婷 | 香蕉成人在线视频 | 国产又粗又黄视频 | 久热网站| 尤物视频免费观看 | 亚洲综合99 | 少妇性l交大片免潘金莲 | 国产女主播视频 | 97国产在线播放 | 精品国产乱码一区二区 | 久久精品亚洲精品国产欧美 | 操穴网站 | 九一国产精品 | 日韩成人精品在线 | www.香蕉.com | 精品久久国产字幕高潮 | 午夜欧美成人 | 国产特级av| 97超碰成人| 中文字幕亚洲综合 | 欧美xxxx精品| 国产乱码av | 国产一区日本 | 秋霞影院午夜丰满少妇在线视频 | 人人爱国产 | 制服丝袜中文字幕在线 | 美女扒开下面让男人捅 | 成人激情五月天 | 亚洲AV午夜精品 | 久久久精品动漫 | 东京热无码av一区二区 | av不卡在线 | 毛茸茸日本熟妇高潮 | 男人的天堂狠狠干 | 欧美乱大交xxxxx潮喷l头像 | 欧洲一区二区三区 | 亚洲一区二区三区激情 | 免费成人在线看 | 国产suv精品一区二区四 | 欧美黄色一级网站 | 欧美日韩性生活视频 | 一级片a级片| 国产福利专区 | 午夜片在线观看 | 99热这里都是精品 | 丝袜美腿中文字幕 | 日韩精品人妻一区二区中文字幕 | 日韩乱论 | av视 | 色综合成人 | 久久久精品久 | 欧美一级免费在线 | 污污视频网站在线免费观看 | 成人动漫一区二区三区 | 欧美拍拍 | 少妇人妻邻居 | 欧美五月婷婷 |