Diagram showing the periodic orbit of a mass-spring system in simple harmonic motion. (The velocity and position axes have been reversed from the standard convention in order to align the two diagrams)
← ホームへ戻る
相空間微分方程式相図力学系ハミルトン力学

相空間と微分方程式動的システムの視覚的解析

🗓 2026年8月11日

物理的なシステムがどのような状態で、時間とともにどう変化するかを数学的に記述する際、「相空間(Phase Space)」という概念が不可欠です。相空間とは、システムが取り得るすべての可能な状態を点として表現した多次元空間のことです。これにより、複雑な数式による記述を幾何学的な軌跡として視覚化でき、システムの本質的な挙動を直感的に理解することが可能になります。

Key Facts

  • 相空間の定義:システムのすべての自由度(位置や運動量など)を軸とする多次元空間。
  • 状態の表現:ある瞬間のシステムの状態は、相空間上の唯一の「点」として対応する。
  • 軌跡(トラジェクトリ):時間の経過に伴う状態の変化は、空間内を移動する曲線として描かれる。
  • 次元数:単純な系では1次元(相線)や2次元(相平面)だが、多粒子系では極めて高次元になる。
  • 応用範囲:古典力学から量子力学、カオス理論、さらには医学的な生理反応の可視化まで幅広く利用される。

相空間の基本原理と構造

相空間の考え方は、19世紀後半にルートヴィヒ・ボルツマン、アンリ・ポアンカレ、ジョサイア・ウィラード・ギブスらによって発展しました。力学的なシステムにおいて、相空間は通常、位置パラメータと運動量パラメータの直積で構成されます。システムの自由度が増えるほど、相空間の次元も増加します。例えば、単一の粒子であれば位置と運動量の組み合わせで表現できますが、多数の分子からなるガスの場合、各粒子の3次元的な位置と運動量をすべて軸に持つ必要があります。

単純な調和振動子のような系では、位置と速度の関係が周期的な閉曲線として描かれます。

Diagram showing the periodic orbit of a mass-spring system in simple harmonic motion. (The velocity and position axes have been reversed from the standard convention in order to align the two diagrams)

このように、相空間上の軌跡を分析することで、初期条件から将来の状態を予測したり、システムがどのような安定状態に向かうのかを明らかにしたりできます。特にロボティクスなどの分野では、ロボットアームの可動範囲や最適な動作経路を決定するために、この幾何学的なアプローチが有効に活用されています。

低次元システムにおける視覚化:相線と相平面

1次元システム(相線)

変数が一つだけの自律的な常微分方程式 $\frac{dy}{dt} = f(y)$ で記述される系は「相線」と呼ばれます。指数関数的な増殖や減衰、あるいはロジスティック成長モデルのように、平衡点の安定性(安定か不安定か)を視覚的に即座に判断できるのが特徴です。

2次元システム(相平面)

位置と速度の2つの変数を持つ系は「相平面」として表現されます。ここでは、システムが特定の閉路に収束する「リミットサイクル」などの複雑な挙動が観察できます。例えば、ファン・デル・ポール振動子の相図では、どのような初期状態から始めても最終的に一定の周期軌道に落ち着く様子が描かれます。

Phase portrait of the Van der Pol oscillator

また、振り子の運動のような系では、角度という変数が $2\pi$ ごとに元に戻るため、相空間は円筒状の構造を持つことになります。

Potential energy and phase portrait of a simple pendulum. Note that the x-axis, being angular, wraps onto itself after every 2π radians.
Illustration of how a phase portrait would be constructed for the motion of a simple pendulum
Time-series flow in phase space specified by the differential equation of a pendulum. The X axis corresponds to the pendulum's position, and the Y axis its speed.

相図(Phase Portrait)と安定性の解析

数学的な文脈において、相空間上の軌道を幾何学的に表現したものを「相図」と呼びます。相図を用いることで、システムにアトラクター(安定した吸収点/シンク)やリペラー(不安定な放出点/ソース)が存在するかを判別できます。

減衰振動子の例では、減衰の強さに応じて軌跡がどのように中心点へ吸い込まれていくかが変化します。これは、微分方程式 $\ddot{x} + 2\gamma\dot{x} + \omega^2x = 0$ のパラメータ $\gamma$ の値によって決定されます。

Phase portrait of damped oscillator, with increasing damping strength. The equation of motion is x ¨ + 2 γ x ˙ + ω 2 x = 0. {\displaystyle {\ddot {x}}+2\gamma {\dot {x}}+\omega ^{2}x=0.}

また、統計力学においては、多数のシステムの集団(アンサンブル)が相空間内でどのように進化するかを研究します。リウヴィルの定理によれば、相空間における点の局所的な密度は一定に保たれますが、時間の経過とともに集団の形状は複雑に歪んでいきます。

Evolution of an ensemble of classical systems in phase space (top). The systems are a massive particle in a one-dimensional potential well (red curve, lower figure). The initially compact ensemble becomes swirled up over time.

広範な学問分野への応用

相空間概念の適用分野と特徴
分野 適用内容・役割 主要な概念
カオス理論 非線形システムの予測不能な挙動の解析 ローレンツ・アトラクター、マンデルブロ集合
量子力学 古典的な位置・運動量を演算子として再構成 ウィグナー分布、Weyl写像、星積
統計力学 微視的状態(マイクロステート)の記述 相積分、分配関数、リウヴィルの定理
熱力学 圧力・温度などの巨視的状態(マクロステート)の表現 相図(状態図)、相転移
光学・医学 非結像光学の照明解析や生理的反応の可視化 ハミルトン光学、多次元生理応答解析

量子力学における相空間

量子力学では、不確定性原理により位置 $q$ と運動量 $p$ を同時に正確に決定することはできません。しかし、WeylやWignerらの研究により、量子力学を相空間上の関数として定式化する手法が確立されました。これにより、量子力学を古典力学の「変形(deformation)」として捉えることができ、プランク定数 $\hbar$ がその変形パラメータの役割を果たします。

Frequently Asked Questions

相空間と相図の違いは何ですか?

相空間は、システムが取り得るすべての状態を包含する「数学的な空間そのもの」を指します。一方、相図(相ポートレート)は、その空間内における代表的な軌道や平衡点を描き出した「視覚的なグラフ」のことを指します。

なぜ相空間を使うと解析が簡単になるのですか?

時間軸に沿った数値の変化を追うのではなく、状態変数同士の関係を幾何学的に見ることで、システムの安定性や周期性、あるいはカオス的な挙動といった「定性的な性質」を一目で把握できるためです。

熱力学で言う「相」と、この「相空間」は同じ意味ですか?

文脈によって異なります。古典力学的な意味では、全粒子の位置と運動量の空間を指します。しかし、熱力学では圧力や温度などのマクロな変数で定義される空間を指し、その中の特定の領域(液体相や固体相など)を「相」と呼びます。

相空間の次元数はどのように決まりますか?

システムの「自由度」によって決まります。一般的に、一つの自由度に対して位置と運動量の2つの変数が対応するため、自由度が $N$ 個の力学系では $2N$ 次元の相空間が必要になります。

References

  1. Nolte, D. D. (2010). "The tangled tale of phase space". Physics Today. 63 (4): 33–38. :2010PhT....63d..33N. :10.1063/1.3397041.  17205307.
  2. Laurendeau, Normand M. (2005). Statistical Thermodynamics: Fundamentals and Applications. New York: Cambridge University Press.  .
  3. Vu-Quoc, L. (2008). "Configuration integral". Archived from the original on April 28, 2012.
  4. Curtright, T. L.; Zachos, C. K. (2012). "Quantum Mechanics in Phase Space". Asia Pacific Physics Newsletter. 01: 37–46. :1104.5269. :10.1142/S2251158X12000069.  119230734.
  5. Chaves, Julio (2015). Introduction to Nonimaging Optics, Second Edition. .  .
  6. Klabukov, I.; Tenchurin, T.; Shepelev, A.; Baranovskii, D.; Mamagulashvili, V.; Dyuzheva, T.; Krasilnikova, O.; Balyasin, M.; Lyundup, A.; Krasheninnikov, M.; Sulina, Y.; Gomzyak, V.; Krasheninnikov, S.; Buzin, A.; Zayratyants, G. (2023). "Biomechanical Behaviors and Degradation Properties of Multilayered Polymer Scaffolds: The Phase Space Method for Bile Duct Design and Bioengineering". Biomedicines. 11 (3): 745. :10.3390/biomedicines11030745.  2227-9059.  10044742.  36979723.
  7. Kirkland, M.A. (2004). "A phase space model of hemopoiesis and the concept of stem cell renewal". Experimental Hematology. 32 (6): 511–519. :10.1016/j.exphem.2004.02.013. :10536/DRO/DU:30101092.  0301-472X.  15183891.

📸 フォトギャラリー

Diagram showing the periodic orbit of a mass-spring system in simple harmonic motion. (The velocity and position axes have been reversed from the standard convention in order to align the two diagrams)
Evolution of an ensemble of classical systems in phase space (top). The systems are a massive particle in a one-dimensional potential well (red curve, lower figure). The initially compact ensemble becomes swirled up over time.
Phase portrait of the Van der Pol oscillator
Potential energy and phase portrait of a simple pendulum. Note that the x-axis, being angular, wraps onto itself after every 2π radians.
Phase portrait of damped oscillator, with increasing damping strength. The equation of motion is x ¨ + 2 γ x ˙ + ω 2 x = 0. {\displaystyle {\ddot {x}}+2\gamma {\dot {x}}+\omega ^{2}x=0.}
Illustration of how a phase portrait would be constructed for the motion of a simple pendulum
Time-series flow in phase space specified by the differential equation of a pendulum. The X axis corresponds to the pendulum's position, and the Y axis its speed.