位相空間ポートレート — 非線形振動・カオスの可視化

力学系の時間発展は、一般に一階の連立常微分方程式（状態方程式）で記述されます。ここでは変位を$x$、速度を$y=\dot{x}$と置いています。

$\dot{x}, \dot{y}$はそれぞれ$x, y$の時間微分、$f$は系に応じた関数、$\boldsymbol{p}$はパラメータ（質量、減衰係数、非線形性の強さなど）を表します。このツールではこの方程式をRK4（4次ルンゲ・クッタ法）で数値的に解き、軌跡を描画しています。

具体例として、減衰振動（ダンパー付きバネ質量系）の方程式は以下の通りです。

ここで、$\gamma$は減衰係数（「パラメータ」の$d$）、$\omega$は固有角周波数です。$\gamma \gt 0$のとき、位相空間の軌跡はらせんを描きながら原点（平衡点）に収束します。これが「安定な固定点アトラクター」の振る舞いです。

機械・構造物の振動解析：自動車のサスペンションや建物の耐震設計では、位相ポートレートを用いて、外部からの衝撃（初期条件）後に振動がどのように減衰して安定点に戻るかを評価します。非線形性を考慮したダフィング振動子のモデルは、大きな変位で剛性が変化する部品の解析に使われます。

電子回路設計：発振回路（オシレータ）の設計では、ファン・デル・ポール型のリミットサイクルが安定した発振状態を保つかどうかを位相空間で確認します。パラメータを変えて望ましい発振波形と安定性を両立させる設計が行われます。

生体現象のモデリング：心筋細胞の電気的活動（心拍）や神経細胞の発火パターンは、非線形振動子モデルで記述されます。位相空間解析は、不整脈などの異常なリズムがどのようなメカニズムで発生するかを理解するのに役立ちます。

気象・流体のカオス解析：ローレンツ系は大気の対流の単純化モデルとして提唱されました。その位相空間に現れる「ストレンジアトラクター」は、わずかな初期値の差が巨大な結果の違いを生む「バタフライ効果」、すなわちカオスの可視化そのものです。長期予測が本質的に困難な現象の理解に貢献しています。

このツールを使い始めるとき、いくつかつまずきやすいポイントがあるよ。まず「初期条件を変えても軌跡が同じに見える」というケース。例えば「ダフィング振動子」で強制振動をオフにし、減衰をゼロにしてみて。初期条件を変えても、全ての軌跡が同じ閉曲線上に乗ってしまうんだ。これは「保存系」だから。エネルギーが保存されるので、初期の「押し加減」で決まるエネルギー一定の曲線（等高線みたいなもの）の上をずっと回り続ける。軌跡が重なっているわけじゃなく、異なるエネルギーを持つ独立した軌跡が、別々の閉曲線を描いているんだ。パラメータ設定で「減衰」を少し入れると、これらの軌跡がらせん状に収束する様子が観察できるから確認してみて。

次に数値計算の精度への過信。このツールはRK4という高精度な手法を使っているけど、特にカオス系（ローレンツ系など）では、ほんの少しの計算誤差が時間とともに爆発的に増幅される（バタフライ効果だね）。だから、長時間シミュレーションした結果の細かい軌跡を「絶対的真実」と思わないこと。実務では、アトラクターの「形」や「性質」（安定か不安定か）といった定性的な振る舞いを見るのに使うのが基本だ。

最後はパラメータ変更の「一気食い」。パラメータをいきなり大きく変えると、系の振る舞いが大きく変わって「何が起こったかわからない」状態になりがち。例えばファン・デル・ポール振動子の非線形性パラメータ$\mu$を、0.1から一気に5.0にすると形が激変するよね。変化を理解するには、0.1, 0.5, 1.0, 2.0…と少しずつ変えながら、リミットサイクルがどう連続的に変形するかを追跡するのがコツ。これで「パラメータの物理的意味」が体感できる。