力学系解析の手順

2次元連続力学系の定性的解析 - Step by Step

解析のワークフロー

  1. 系を定義: ẋ = f(x,y), ẏ = g(x,y) を書き下す
  2. 固定点を求める: f(x,y) = 0 かつ g(x,y) = 0 を解く
  3. ヤコビ行列を計算: J = [∂f/∂x, ∂f/∂y; ∂g/∂x, ∂g/∂y]
  4. 各固定点でJを評価: J|(x*,y*)
  5. 固有値を求める: τ = tr(J), Δ = det(J) から分類
  6. 安定性を判定: 固有値の実部の符号で決定
  7. ヌルクラインを描く: 相図の骨格を把握
  8. 位相図をスケッチ: 固定点・ヌルクライン・軌道を統合

例題を選択

0

系の定義

解析する2次元連続力学系を定義する:

ẋ = f(x, y) = x(3 - x - 2y)
ẏ = g(x, y) = y(2 - x - y)

物理的意味: 2種の競争系。x, y は種の個体数密度。

1

固定点を求める

固定点 (x*, y*) は ẋ = 0 かつ ẏ = 0 を満たす点:

手順:
  1. f(x, y) = 0 と g(x, y) = 0 を連立
  2. 代数的に解く(因数分解、代入法など)
  3. 全ての解を列挙(境界上の解も忘れずに)

計算過程

f(x,y) = x(3 - x - 2y) = 0 より: x = 0 または x + 2y = 3
g(x,y) = y(2 - x - y) = 0 より: y = 0 または x + y = 2

固定点の一覧

固定点 x* y* 物理的意味
P₁00両種絶滅
P₂30種xのみ生存
P₃02種yのみ生存
P₄11共存
2

ヤコビ行列を計算

系を線形化するために、ヤコビ行列(偏微分の行列)を計算:

J = [ ∂f/∂x∂f/∂y ∂g/∂x∂g/∂y ]

各偏微分の計算

f(x,y) = x(3 - x - 2y) = 3x - x² - 2xy

∂f/∂x = 3 - 2x - 2y
∂f/∂y = -2x

g(x,y) = y(2 - x - y) = 2y - xy - y²

∂g/∂x = -y
∂g/∂y = 2 - x - 2y

一般的なヤコビ行列

J(x,y) = [ 3 - 2x - 2y-2x -y2 - x - 2y ]
3

各固定点でヤコビ行列を評価

各固定点 (x*, y*) の座標をヤコビ行列に代入:

P₁ = (0, 0):
J(0,0) = [ 30 02 ]
P₂ = (3, 0):
J(3,0) = [ -3-6 0-1 ]
P₃ = (0, 2):
J(0,2) = [ -10 -2-2 ]
P₄ = (1, 1):
J(1,1) = [ -1-2 -1-1 ]
4

固有値を計算

2×2行列の固有値は、トレースと行列式から簡単に求まる:

τ = tr(J) = a + d(対角成分の和)
Δ = det(J) = ad - bc(行列式)

λ = (τ ± √(τ² - 4Δ)) / 2
判別式の意味:
  • τ² - 4Δ > 0: 実固有値(ノード)
  • τ² - 4Δ < 0: 複素固有値(スパイラル)
  • τ² - 4Δ = 0: 重根(退化ノード、星型ノード)

計算結果

固定点 τ (trace) Δ (det) τ² - 4Δ λ₁ λ₂
P₁(0,0) 5 6 1 > 0 3 2
P₂(3,0) -4 3 4 > 0 -3 -1
P₃(0,2) -3 2 1 > 0 -2 -1
P₄(1,1) -2 -1 8 > 0 -1+√2 ≈ 0.41 -1-√2 ≈ -2.41
5

固定点を分類・安定性判定

分類チャート(τ-Δ平面)

条件 分類
Δ < 0 サドル点
Δ > 0, τ < 0, τ²>4Δ 安定ノード
Δ > 0, τ > 0, τ²>4Δ 不安定ノード
Δ > 0, τ < 0, τ²<4Δ 安定スパイラル
Δ > 0, τ > 0, τ²<4Δ 不安定スパイラル
Δ > 0, τ = 0 中心点

安定性の判定基準

漸近安定: 全ての固有値の実部 < 0

⟺ τ < 0 かつ Δ > 0

不安定: 少なくとも1つの固有値の実部 > 0

⟺ τ > 0 または Δ < 0

各固定点の分類結果

P₁ = (0, 0): 不安定ノード

λ₁ = 3 > 0, λ₂ = 2 > 0(両方正)→ 全方向に発散

P₂ = (3, 0): 安定ノード

λ₁ = -3 < 0, λ₂ = -1 < 0(両方負)→ 全方向から収束

P₃ = (0, 2): 安定ノード

λ₁ = -2 < 0, λ₂ = -1 < 0(両方負)→ 全方向から収束

P₄ = (1, 1): サドル点

λ₁ ≈ 0.41 > 0, λ₂ ≈ -2.41 < 0(異符号)→ 不安定

6

ヌルクラインを描く

ヌルクラインは相図の「骨格」を与える:

x-ヌルクライン(ẋ = 0)

f(x,y) = x(3 - x - 2y) = 0
→ x = 0 または y = (3-x)/2

この曲線上では、軌道は垂直に横切る

y-ヌルクライン(ẏ = 0)

g(x,y) = y(2 - x - y) = 0
→ y = 0 または y = 2 - x

この曲線上では、軌道は水平に横切る

ヌルクラインの交点 = 固定点
両方のヌルクラインが交わる点では ẋ = ẏ = 0 となり、軌道は停止する。
7

位相図を描く

相平面(x-y)
x-ヌルクライン
y-ヌルクライン
安定固定点
不安定固定点
サドル点
8

解析のまとめ

系の大域的振る舞い

固定点のまとめ

固定点 分類 安定性 物理的解釈
(0, 0) 不安定ノード 不安定 両種絶滅状態は不安定
(3, 0) 安定ノード 安定 種xが種yを競争排除
(0, 2) 安定ノード 安定 種yが種xを競争排除
(1, 1) サドル点 不安定 共存は不安定(双安定系)

結論

この系は双安定(bistable)である。

初期条件により、種xが勝つ(3,0)か種yが勝つ(0,2)のどちらかに収束する。

共存点(1,1)はサドル点であり、その安定多様体が2つの吸引域の境界(分離線)となる。

解析のチェックリスト

必須項目

  • 全ての固定点を見つけたか?
  • 各固定点でヤコビ行列を評価したか?
  • 固有値を正しく計算したか?
  • 分類は τ-Δ チャートと一致するか?
  • ヌルクラインの交点と固定点は一致するか?

よくあるミス

  • 境界上の固定点(x=0, y=0)を見落とす
  • ヤコビ行列の偏微分を間違える
  • 行列式の符号を間違える
  • 中心点を「安定」と誤判定
  • 大域的な振る舞いを考慮しない