動(dòng)力系統(tǒng):奇點(diǎn)類型與軌線分布

奇點(diǎn)類型

本篇僅討論平面自治系統(tǒng)的奇點(diǎn)類型傻挂,對(duì)于線性化系統(tǒng)求解其特征值蠢熄。特征值的實(shí)部和虛部決定了奇點(diǎn)的類型:

  • 實(shí)特征值
    如果所有特征值均為負(fù)實(shí)數(shù)逗噩,則該奇點(diǎn)是穩(wěn)定結(jié)點(diǎn)平窘;
    如果所有特征值均為正實(shí)數(shù)吓肋,則是不穩(wěn)定結(jié)點(diǎn);
    如果有正有負(fù)瑰艘,則是鞍點(diǎn)是鬼。
  • 復(fù)特征值:如果特征值是復(fù)數(shù),則需要看其實(shí)部紫新。
    如果實(shí)部為負(fù)均蜜,奇點(diǎn)是穩(wěn)定的螺旋點(diǎn);
    如果實(shí)部為正芒率,則是不穩(wěn)定的螺旋點(diǎn)囤耳。

例題:線性系統(tǒng)

系統(tǒng) 一:

\begin{cases} \frac{dx}{dt} = -4x + y, \\ \frac{dy}{dt} = -3x. \end{cases}

第一步:確定雅可比矩陣

這個(gè)系統(tǒng)的右側(cè)為:
f(x, y) = \begin{bmatrix} -4x + y \\ -3x \end{bmatrix}

對(duì) f(x, y) 求關(guān)于 xy 的偏導(dǎo)數(shù),得到雅可比矩陣 J

J = \begin{bmatrix} \frac{\partial f_1}{\partial x} & \frac{\partial f_1}{\partial y} \\ \frac{\partial f_2}{\partial x} & \frac{\partial f_2}{\partial y} \end{bmatrix} = \begin{bmatrix} -4 & 1 \\ -3 & 0 \end{bmatrix}

第二步:求特征值

我們接下來計(jì)算雅可比矩陣 J 的特征值。特征值 \lambda 滿足以下特征方程:

\det(J - \lambda I) = 0

J 和單位矩陣 I 帶入:

\begin{vmatrix} -4 - \lambda & 1 \\ -3 & -\lambda \end{vmatrix} = 0

計(jì)算行列式:

(-4 - \lambda)(-\lambda) - (1)(-3) = 0

展開并簡(jiǎn)化:

4\lambda + \lambda^2 + 3 = 0

這得到一個(gè)二次方程:

\lambda^2 + 4\lambda + 3 = 0

解這個(gè)二次方程紫皇,使用求根公式:

\lambda = \frac{-4 \pm \sqrt{(4)^2 - 4 \cdot 1 \cdot 3}}{2 \cdot 1} = \frac{-4 \pm \sqrt{16 - 12}}{2} = \frac{-4 \pm 2}{2}

所以得到兩個(gè)特征值:

\lambda_1 = -1, \quad \lambda_2 = -3

第三步:判斷奇點(diǎn)類型

  • 兩個(gè)特征值 \lambda_1 = -1\lambda_2 = -3 都是負(fù)實(shí)數(shù)慰安。
  • 這意味著奇點(diǎn)是一個(gè) 穩(wěn)定節(jié)點(diǎn),因?yàn)樗刑卣髦刀紴樨?fù)聪铺,軌線將向奇點(diǎn)收斂化焕。

因此,系統(tǒng) (1) 的奇點(diǎn)是一個(gè) 穩(wěn)定節(jié)點(diǎn)铃剔,軌線在奇點(diǎn)附近向奇點(diǎn)收斂撒桨。


系統(tǒng)二

\begin{cases} \frac{dx}{dt} = -x - 5 + y, \\ \frac{dy}{dt} = -3x. \end{cases}

第一步:確定雅可比矩陣

這個(gè)系統(tǒng)的右側(cè)為:
f(x, y) = \begin{bmatrix} -x - 5 + y \\ -3x \end{bmatrix}

對(duì) f(x, y) 求關(guān)于 xy 的偏導(dǎo)數(shù),得到雅可比矩陣 J
J = \begin{bmatrix} -1 & 1 \\ -3 & 0 \end{bmatrix}

第二步:求特征值

特征值 \lambda 滿足以下特征方程:
\det(J - \lambda I) = 0

計(jì)算行列式:
\begin{vmatrix} -1 - \lambda & 1 \\ -3 & -\lambda \end{vmatrix} = 0

展開行列式:
(-1 - \lambda)(-\lambda) - (1)(-3) = 0

簡(jiǎn)化得到:
\lambda^2 + \lambda - 3 = 0

使用求根公式解此二次方程:
\lambda = \frac{-1 \pm \sqrt{1 + 4 \cdot 3}}{2} = \frac{-1 \pm \sqrt{13}}{2}

因此键兜,特征值為:
\lambda_1 = \frac{-1 + \sqrt{13}}{2}, \quad \lambda_2 = \frac{-1 - \sqrt{13}}{2}

第三步:判斷奇點(diǎn)類型

  • 由于特征值的符號(hào)不同(一個(gè)為正凤类,一個(gè)為負(fù)),奇點(diǎn)是 鞍點(diǎn)普气。
  • 鞍點(diǎn)表示系統(tǒng)的奇點(diǎn)在某些方向上不穩(wěn)定谜疤。

系統(tǒng)三:

\begin{cases} \frac{dx}{dt} = x - \frac{3y}{4}, \\ \frac{dy}{dt} = 7x - 4y. \end{cases}

第一步:確定雅可比矩陣

J = \begin{bmatrix} 1 & -\frac{3}{4} \\ 7 & -4 \end{bmatrix}

第二步:求特征值

特征值 \lambda 滿足以下特征方程:
\det(J - \lambda I) = 0

計(jì)算行列式:
\begin{vmatrix} 1 - \lambda & -\frac{3}{4} \\ 7 & -4 - \lambda \end{vmatrix} = 0

展開行列式:
(1 - \lambda)(-4 - \lambda) + \frac{21}{4} = 0

整理得到:
\lambda^2 + 3\lambda - \frac{7}{4} = 0

解得特征值為復(fù)數(shù),說明系統(tǒng)具有螺旋奇點(diǎn)现诀。


系統(tǒng)四

\begin{cases} \frac{dx}{dt} = -x - y + 1, \\ \frac{dy}{dt} = x - y - 5. \end{cases}

第一步:確定雅可比矩陣

系統(tǒng)的右側(cè)為:
f(x, y) = \begin{bmatrix} -x - y + 1 \\ x - y - 5 \end{bmatrix}

對(duì) f(x, y) 求關(guān)于 xy 的偏導(dǎo)數(shù)夷磕,得到雅可比矩陣 J
J = \begin{bmatrix} -1 & -1 \\ 1 & -1 \end{bmatrix}

第二步:求特征值

特征值 \lambda 滿足以下特征方程:
\det(J - \lambda I) = 0

計(jì)算行列式:
\begin{vmatrix} -1 - \lambda & -1 \\ 1 & -1 - \lambda \end{vmatrix} = 0

展開行列式:
(-1 - \lambda)(-1 - \lambda) - (-1)(1) = 0

整理得到:
\lambda^2 + 2\lambda = 0

將方程分解為:
\lambda(\lambda + 2) = 0

因此,特征值為:
\lambda_1 = 0, \quad \lambda_2 = -2

第三步:判斷奇點(diǎn)類型

  • 特征值中包含 0 和一個(gè)負(fù)特征值 -2仔沿,因此系統(tǒng)的奇點(diǎn)是一個(gè) 鞍結(jié)點(diǎn)(即鞍點(diǎn)和中心點(diǎn)的混合類型)坐桩。
  • 鞍結(jié)點(diǎn)表示系統(tǒng)的奇點(diǎn)在某些方向上可能穩(wěn)定,而在其他方向上不穩(wěn)定封锉。

最后绵跷,我們用Python來繪制了上面四個(gè)系統(tǒng)的軌線分布圖:


image.png

例題:非線性系統(tǒng)

對(duì)于非線性系統(tǒng),我們一般采取線性近似法來近似其在奇點(diǎn)附近的狀態(tài)成福。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末碾局,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子闷叉,更是在濱河造成了極大的恐慌擦俐,老刑警劉巖,帶你破解...
    沈念sama閱讀 221,548評(píng)論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件握侧,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡嘿期,警方通過查閱死者的電腦和手機(jī)品擎,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,497評(píng)論 3 399
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來备徐,“玉大人萄传,你說我怎么就攤上這事。” “怎么了秀菱?”我有些...
    開封第一講書人閱讀 167,990評(píng)論 0 360
  • 文/不壞的土叔 我叫張陵振诬,是天一觀的道長(zhǎng)。 經(jīng)常有香客問我衍菱,道長(zhǎng)赶么,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 59,618評(píng)論 1 296
  • 正文 為了忘掉前任脊串,我火速辦了婚禮辫呻,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘琼锋。我一直安慰自己放闺,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,618評(píng)論 6 397
  • 文/花漫 我一把揭開白布缕坎。 她就那樣靜靜地躺著怖侦,像睡著了一般。 火紅的嫁衣襯著肌膚如雪谜叹。 梳的紋絲不亂的頭發(fā)上础钠,一...
    開封第一講書人閱讀 52,246評(píng)論 1 308
  • 那天,我揣著相機(jī)與錄音叉谜,去河邊找鬼旗吁。 笑死,一個(gè)胖子當(dāng)著我的面吹牛停局,可吹牛的內(nèi)容都是我干的很钓。 我是一名探鬼主播,決...
    沈念sama閱讀 40,819評(píng)論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼董栽,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼码倦!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起锭碳,我...
    開封第一講書人閱讀 39,725評(píng)論 0 276
  • 序言:老撾萬榮一對(duì)情侶失蹤袁稽,失蹤者是張志新(化名)和其女友劉穎,沒想到半個(gè)月后擒抛,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體推汽,經(jīng)...
    沈念sama閱讀 46,268評(píng)論 1 320
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,356評(píng)論 3 340
  • 正文 我和宋清朗相戀三年歧沪,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了歹撒。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 40,488評(píng)論 1 352
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡诊胞,死狀恐怖暖夭,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤迈着,帶...
    沈念sama閱讀 36,181評(píng)論 5 350
  • 正文 年R本政府宣布竭望,位于F島的核電站,受9級(jí)特大地震影響裕菠,放射性物質(zhì)發(fā)生泄漏咬清。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,862評(píng)論 3 333
  • 文/蒙蒙 一糕韧、第九天 我趴在偏房一處隱蔽的房頂上張望枫振。 院中可真熱鬧,春花似錦萤彩、人聲如沸粪滤。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,331評(píng)論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽杖小。三九已至,卻和暖如春愚墓,著一層夾襖步出監(jiān)牢的瞬間予权,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,445評(píng)論 1 272
  • 我被黑心中介騙來泰國(guó)打工浪册, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留扫腺,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,897評(píng)論 3 376
  • 正文 我出身青樓村象,卻偏偏與公主長(zhǎng)得像笆环,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子厚者,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,500評(píng)論 2 359

推薦閱讀更多精彩內(nèi)容