› 位相空間ポートレート 戻る
Nonlinear Dynamics Simulator

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

単振動・減衰・ファン・デル・ポール・ダフィング・ローレンツ系の位相空間軌跡をRK4で数値積分。アトラクター・リミットサイクル・カオスをリアルタイムに探索しよう。

使い方
• キャンバスをクリック → その点から軌跡を開始
• 複数の軌跡を同時表示可能
• パラメータ変更で系の挙動が変化
方程式
軌跡数: 0
積分法: RK4 / dt = 0.01

位相空間とは何か

力学系の状態は、位置 $x$ と速度 $\dot{x}$ の組 $(x, \dot{x})$ で完全に決まります。この2つを座標軸にとった平面を位相空間（phase space）といい、系の時間発展は位相空間上の曲線（軌跡）として描かれます。

各系の方程式と特徴

1. 単純調和振動子
$$\dot{x}=y,\quad \dot{y}=-\omega^2 x$$

エネルギー保存系。位相空間では閉じた楕円を描く。同心楕円族を形成し、軌跡は交わらない。

2. 減衰振動子
$$\dot{x}=y,\quad \dot{y}=-2\gamma y - \omega^2 x$$

エネルギー散逸系。軌跡は原点（安定固定点）に螺旋状に収束。$\gamma < \omega$ で振動的減衰、$\gamma > \omega$ で過減衰。

3. ファン・デル・ポール振動子
$$\dot{x}=y,\quad \dot{y}=\mu(1-x^2)y - x$$

$|x|<1$ では負の減衰（エネルギー注入）、$|x|>1$ では正の減衰（散逸）。すべての軌跡がリミットサイクルに収束。真空管・心拍のモデル。

4. ダフィング振動子（周期外力付き）
$$\dot{x}=y,\quad \dot{y}=-\delta y - \alpha x - \beta x^3 + \gamma\cos(\omega t)$$

非線形ばねと周期外力の組み合わせ。適切なパラメータでストレンジアトラクターとカオスが現れる。初期値鋭敏性（バタフライ効果）を視覚的に確認できる。

5. ローレンツアトラクター（x-z 投影）
$$\dot{X}=\sigma(Y-X),\quad \dot{Y}=X(\rho-Z)-Y,\quad \dot{Z}=XY-\beta Z$$

大気対流の簡略モデル。3次元の軌跡をx-z平面に投影表示。カオス的な挙動と「バタフライ型」の奇妙なアトラクターを観察できる。

Q & A

🧑‍🎓 位相空間ポートレートって、ふつうの時間グラフと何が違うんですか？

🎓 時間グラフは「横軸が時間、縦軸が変位」だけど、位相空間は「横軸が変位、縦軸が速度」なんだ。時間の情報を消した代わりに、系の定性的な構造が見えやすくなる。単振動なら時間グラフはサインカーブだけど、位相空間では閉じた楕円になる。楕円が閉じているということは「周期運動」を意味するわけ。

🧑‍🎓 なるほど！じゃあ減衰振動だと楕円が閉じないんですね？

🎓 そう、原点に向かって螺旋を巻きながら収束する。その収束先の点を「安定固定点」または「アトラクター」と呼ぶ。実際の機械振動では、このアトラクターの場所と収束速度が設計の核心で、「どれくらいの減衰係数で何秒後に振動が止まるか」がこのグラフから直感的にわかる。

🧑‍🎓 ファン・デル・ポール振動子のリミットサイクルって、なぜ「自然に」その軌跡に落ち着くんですか？

🎓 ミソは $\mu(1-x^2)y$ という項だよ。$|x|<1$ のとき $(1-x^2)>0$ だから速度 $y$ と同じ向きの力が加わる → エネルギーが注入される。$|x|>1$ のときは逆向き → エネルギーが散逸される。この2つのバランスが取れた閉曲線がリミットサイクルで、内側から来ても外側から来ても必ずそこに収束するんだ。

🧑‍🎓 心拍のモデルになるって書いてましたが、具体的にはどういうことですか？

🎓 心臓は「外から刺激がなくても一定のリズムで勝手に拍動し続ける」自律的な振動子なんだ。少し乱れても元のリズムに戻る性質、これがまさにリミットサイクルの特徴。$\mu$ を大きくしてみると波形がノコギリ状になる。あれが弛緩振動（relaxation oscillation）で、実際の心電図に近い形になる。

🧑‍🎓 カオスって「予測不可能」ということですか？ランダムとは違うの？

🎓 大事な区別だね。カオスは決定論的なんだ。方程式は完全に決まっていて、初期条件が同じなら全く同じ軌跡が再現される。でも初期条件がほんのわずかでも違うと、軌跡が指数関数的に離れていく。ランダムは「本質的に予測不可能」だけど、カオスは「原理的には予測可能、でも実用上は不可能」という違いがある。気象予報が「2週間先まで」しか信頼できないのはこのせいだよ。

位相空間ポートレートとは

🧑‍🎓
位相空間って何ですか？普通のグラフと何が違うんですか？
🎓
ざっくり言うと、運動の「状態」を一発で見られる便利なグラフだよ。普通のグラフは「時間」と「変位」をプロットするけど、位相空間は「変位」と「速度」を軸にとるんだ。例えばこのシミュレーターで「単振動」を選んでみて。軌跡が綺麗な楕円になるでしょ？これが「状態が時間とともにどう動くか」を一つの図で表しているんだ。
🧑‍🎓
え、じゃあその楕円の大きさは何で決まるんですか？上の「初期条件」のスライダーを動かすと変わりますね。
🎓
その通り！初期条件、つまり最初の「押し加減」で決まるんだ。スライダーで初期位置$x_0$や初速度$y_0$を変えると、楕円の大きさや形が変わる。これが「初期条件の違いがそのまま軌跡の違いになる」ってこと。実務では、設計した機械がどんな初期乱れからも安全な動きに戻るか、この軌跡の収束の仕方を見て評価するんだ。
🧑‍🎓
「ファン・デル・ポール」に変えると、ぐるぐる巻きついていく軌跡が出てきました！これがリミットサイクル？
🎓
鋭いね！それが「リミットサイクル」というアトラクターだ。パラメータ$\mu$をいじってみて。$\mu$が小さいとほぼ円に近いけど、大きくなると歪な形になるだろ？これは非線形性の強さを表している。電子回路の発振器や心臓の拍動モデルなど、一定のリズムを生み出す現象の本質を、この閉じた軌道が捉えているんだよ。

物理モデルと主要な数式

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

$$ \dot{x}= y, \quad \dot{y}= f(x, y, t; \boldsymbol{p}) $$

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

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

$$ \dot{x}= y, \quad \dot{y} = -2\gamma y - \omega^2 x $$

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

実世界での応用

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

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

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

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

よくある誤解と注意点

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

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

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

関連する工学分野

この位相空間の考え方は、一見関係なさそうな分野にも広く応用されているんだ。例えば制御工学では「状態空間」と呼ばれ、まさに位相空間そのもの。ドローンやロボットアームの姿勢（変位）と角速度（速度）を状態ベクトルとして捉え、最適な制御入力を設計するLQR（線形二次レギュレータ）などの基礎になっている。このシミュレーターで減衰振動の軌跡が原点に収束する様子を見るのは、制御系の「安定性」を視覚的に学ぶ第一歩だ。

電気電子工学では、先ほど出た発振回路だけでなく、PLL（位相同期回路）の引き込み現象や、スイッチング電源の安定性解析にも使われる。電圧と電流を位相平面の軸にとり、スイッチのON/OFFでジャンプする軌跡を追う「不連続系」の解析だ。また、メカトロニクスでは、モーターの回転角と角速度の関係を位相平面で分析し、位置決め制御の応答を評価する。

さらに化学工学の連続 stirred-tank reactor (CSTR) では、反応物濃度と温度がカップリングした非線形振動（化学振動）を起こすことがあり、その安定性を位相ポートレートで調べる。生物学以外にも、経済物理学で市場の変動を非線形振動子モデルで表現する試みなど、応用範囲はどんどん広がっているんだ。

発展的な学習のために

このツールに慣れて「もっと知りたい」と思ったら、次のステップに進んでみよう。まず数学的背景としては、「線形代数」と「常微分方程式」の基礎知識が必須。特に、平衡点（固定点）近傍での振る舞いを調べる線形化（ヤコビ行列の計算）と固有値解析を学べば、「なぜ原点に引き込まれるのか」「なぜぐるぐる回るのか」が数式ではっきり理解できる。例えば、減衰振動の線形化行列の固有値が複素数で、かつ実部が負だから「減衰振動しながら収束する」と説明できるようになる。

次に挑戦したいのは「分岐解析」だ。このツールでパラメータを連続的に変えると、ある臨界点を境に系の振る舞いが突然変わることがある（例：安定な固定点が不安定化し、代わりにリミットサイクルが生まれる）。これがホップ分岐と呼ばれる現象で、ファン・デル・ポール振動子で$\mu$を0から正にしていくとまさにこれが起きる。分岐図を描くことで、系の振る舞いの地図を作成できるようになる。

最後に、カオスへの入門として「ストレンジアトラクター」の性質を深掘りしよう。ローレンツアトラクターはフラクタル構造（自己相似性）を持ち、初期値敏感性を示す。これを理解するには、リアプノフ指数（軌跡の発散/収束の平均的な率を表す）やポアンカレ写像（連続的な軌跡から離散的な点の集合を得るテクニック）といった概念に触れてみることをお勧めする。これらの道具を使うことで、一見ランダムな振る舞いの中に潜む秩序を読み解く力が身につくんだ。