动力系统:奇点类型与轨线分布

奇点类型

本篇仅讨论平面自治系统的奇点类型,对于线性化系统求解其特征值。特征值的实部和虚部决定了奇点的类型:

  • 实特征值
    如果所有特征值均为负实数,则该奇点是稳定结点;
    如果所有特征值均为正实数,则是不稳定结点;
    如果有正有负,则是鞍点。
  • 复特征值:如果特征值是复数,则需要看其实部。
    如果实部为负,奇点是稳定的螺旋点;
    如果实部为正,则是不稳定的螺旋点。

例题:线性系统

系统 一:

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

第一步:确定雅可比矩阵

这个系统的右侧为:
f(x, y) = \begin{bmatrix} -4x + y \\ -3x \end{bmatrix}

f(x, y) 求关于 xy 的偏导数,得到雅可比矩阵 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 的特征值。特征值 \lambda 满足以下特征方程:

\det(J - \lambda I) = 0

J 和单位矩阵 I 带入:

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

计算行列式:

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

展开并简化:

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

这得到一个二次方程:

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

解这个二次方程,使用求根公式:

\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}

所以得到两个特征值:

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

第三步:判断奇点类型

  • 两个特征值 \lambda_1 = -1\lambda_2 = -3 都是负实数。
  • 这意味着奇点是一个 稳定节点,因为所有特征值都为负,轨线将向奇点收敛。

因此,系统 (1) 的奇点是一个 稳定节点,轨线在奇点附近向奇点收敛。


系统二

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

第一步:确定雅可比矩阵

这个系统的右侧为:
f(x, y) = \begin{bmatrix} -x - 5 + y \\ -3x \end{bmatrix}

f(x, y) 求关于 xy 的偏导数,得到雅可比矩阵 J
J = \begin{bmatrix} -1 & 1 \\ -3 & 0 \end{bmatrix}

第二步:求特征值

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

计算行列式:
\begin{vmatrix} -1 - \lambda & 1 \\ -3 & -\lambda \end{vmatrix} = 0

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

简化得到:
\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}

第三步:判断奇点类型

  • 由于特征值的符号不同(一个为正,一个为负),奇点是 鞍点
  • 鞍点表示系统的奇点在某些方向上不稳定。

系统三:

\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

计算行列式:
\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

解得特征值为复数,说明系统具有螺旋奇点。


系统四

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

第一步:确定雅可比矩阵

系统的右侧为:
f(x, y) = \begin{bmatrix} -x - y + 1 \\ x - y - 5 \end{bmatrix}

f(x, y) 求关于 xy 的偏导数,得到雅可比矩阵 J
J = \begin{bmatrix} -1 & -1 \\ 1 & -1 \end{bmatrix}

第二步:求特征值

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

计算行列式:
\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

第三步:判断奇点类型

  • 特征值中包含 0 和一个负特征值 -2,因此系统的奇点是一个 鞍结点(即鞍点和中心点的混合类型)。
  • 鞍结点表示系统的奇点在某些方向上可能稳定,而在其他方向上不稳定。

最后,我们用Python来绘制了上面四个系统的轨线分布图:


image.png

例题:非线性系统

对于非线性系统,我们一般采取线性近似法来近似其在奇点附近的状态。

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

推荐阅读更多精彩内容