常微分方程式・相平面解析シミュレーター

相平面解析の対象は、時間$t$に陽に依存しない2次元自律系（Autonomous System）です。状態の時間変化は以下の連立常微分方程式で記述されます。

$$ \begin{aligned}\frac{dx}{dt}&= f(x, y) \\ \frac{dy}{dt}&= g(x, y) \end{aligned}$$

$x, y$: 系の状態変数（例：位置と速度、捕食者と被食者の個体数）。
$f, g$: 状態$(x,y)$によって決まる変化率を表す関数（多くの場合非線形）。
時間$t$でパラメータ化された解$(x(t), y(t))$を$(x,y)$平面上に描いた曲線を「軌跡」または「軌道」と呼びます。

平衡点$(x^ , y^ )$（$f(x^ , y^ )=0$かつ$g(x^ , y^ )=0$を満たす点）近傍の振る舞いを調べるために線形化を行い、その安定性はヤコビアン（Jacobi）行列$J$の固有値によって分類できます。

$$ J = \begin{pmatrix}\dfrac{\partial f}{\partial x}& \dfrac{\partial f}{\partial y}\\[8pt] \dfrac{\partial g}{\partial x}& \dfrac{\partial g}{\partial y}\end{pmatrix}_{(x^, y^)}$$

固有値$\lambda_1, \lambda_2$の実部$\mathrm{Re}(\lambda)$に注目します。
安定ノード/焦点: $\mathrm{Re}(\lambda_1), \mathrm{Re}(\lambda_2) \lt 0$ (全ての軌跡が平衡点に収束)。
不安定ノード/焦点: $\mathrm{Re}(\lambda_1), \mathrm{Re}(\lambda_2) \gt 0$ (平衡点から遠ざかる)。
鞍点: $\mathrm{Re}(\lambda_1) \cdot \mathrm{Re}(\lambda_2) \lt 0$ (一部の方向では収束、他では発散、常に不安定)。
中心: $\mathrm{Re}(\lambda_1) = \mathrm{Re}(\lambda_2) = 0$ (線形近似では中立安定、周囲を周期軌道)。

制御工学・ロボティクス：倒立振り子やドローンの姿勢制御など、非線形システムのコントローラ設計において、平衡点（目標状態）の安定性と収束領域を相平面で解析します。リミットサイクル（安定な周期軌道）の検出は、望ましい周期動作の設計に役立ちます。

生態学・人口動態：ロトカ・ボルテラ方程式は、被食者と捕食者の個体数変動をモデル化した古典的な例です。相平面で平衡点（共存平衡）の安定性を調べることで、個体数が周期的に変動し続ける条件や、一方が絶滅する条件を予測できます。

疫学・SIRモデル：感染症の流行過程を、感受性者(S)、感染者(I)、回復者(R)の3集団の動態としてモデル化します。相平面（例えばS-I平面）を用いることで、流行のピークや終息条件、基本再生産数$R_0$の影響を視覚的に理解できます。

構造工学・非線形振動：橋梁や建築物の風致振動、機械部材の大きな変形を伴う振動では、線形理論では説明できない現象（ジャンプ、サブハーモニック振動など）が起きます。これをダフィング振子やファン・デル・ポール振子などのモデルで表現し、相平面解析で安定性や分岐を調べます。

まず、「相平面は単なる軌跡のプロットではない」という点を押さえましょう。軌跡は初期値から始まる一本の線ですが、相平面解析の真価は、方向場（矢印）全体を見て「流れ」を把握することにあります。一つの軌跡だけを見て「この系は振動する」と早合点するのは危険。例えば減衰振動でも、初期エネルギーが大きければ一見複雑な動きをしますが、方向場全体を見れば、すべての矢印が最終的に一点（安定平衡点）に吸い込まれていく様子がわかります。

次に、パラメータ設定の「スケール」。例えば「振り子」モデルで、減衰係数bを0.1から1.0に変えると、挙動が劇的に変わります。しかし、表示範囲Rがデフォルトの±5のままだと、矢印が全て平衡点周辺に密集してしまい、流れが見えなくなることがあります。パラメータを変えたら、「表示範囲±R」スライダーで視野を広げたり狭めたりするのがコツです。実務でも、適切なスケーリングは必須です。

最後に、「リアルタイムシミュレーションは万能ではない」という落とし穴。このツールは数値計算で軌跡を描いています。非常に鋭敏な系（例えばサドル点の不安定多様体付近）では、ほんの少しの数値誤差で軌跡が全く異なる方向に行ってしまうことがあります。シミュレーション結果を盲信せず、「この結果は本当に信頼できるか？」と常に理論（ヌルクラインの形状、平衡点の線形化予測）と照らし合わせる習慣をつけましょう。