ダフィング方程式

カオス挙動を示す強制ダフィング方程式のポアンカレ切断面

数学におけるダフィング方程式（ダフィングほうていしき、テンプレート:Lang-en-short）あるいはダフィング振動子（Duffing oscillator）は、ある減衰的な駆動振動子をモデル化するために用いられる非線型の二階常微分方程式である。次の式で与えられる：

$\ddot{x}+\delta \dot{x}+\alpha x+\beta x^3=\gamma \cos (\omega t).$

ここで（未知）函数 テンプレート:Math は時間 テンプレート:Math での位置、$\dot{x}$ は テンプレート:Math の時間に関する一階導函数、すなわち速度で、$\ddot{x}$ は テンプレート:Math の時間に関する二階導函数、すなわち加速度である。数 $\delta ,\alpha ,\beta ,\gamma$ および $\omega$ は与えられた定数である。

この式は、（$\beta =\delta =0$ の場合に対応する）単振動よりも複雑なポテンシャルを持つ減衰振動子の動きを表す。例えば、物理学の言葉で言うと、ばねの剛性がフックの法則に従わないテンプレート:仮リンクのモデルと見なされる。

ダフィング方程式は、カオス的挙動を示す力学系の一例である。ジャパニーズ・アトラクタがダフィング方程式におけるカオスの例としてよく知られている。さらにダフィングシステムは、周波数ヒステリシスの挙動のような、跳躍共振現象を周波数反応において示すものである。

パラメータ

• $\delta$ は減衰の大きさを制御する。
• $\alpha$ は剛性の大きさを制御する。
• $\beta$ は復元力に含まれる非線型性の量を制御する。$\beta =0$ であるなら、ダフィング方程式は減衰かつ駆動（driven）な単振動子を表す。
• $\gamma$ は周期駆動力の振幅を制御する。$\gamma =0$ であるなら、駆動力の無いシステムとなる。
• $\omega$ は周期駆動力の周波数を制御する。

解の方法

$\gamma >0$ に対するダフィング振動子のリミットサイクル

$\gamma >0$ に対するダフィング振動子のリミットサイクルの相プロット

$\gamma <0$ に対するダフィング振動子のカオス振動

$\gamma <0$ に対するダフィング振動子のアトラクターのアニメーション

一般に、ダフィング方程式の厳密な記号解が得られるとは限らない。しかし、以下のような多くの近似手法が利用できる：

• フーリエ級数展開は、任意の正確さでの方程式の運動を与える。
• ダフィング項と呼ばれる $x^3$ の項は、小さいものとして近似でき、システムは摂動された単振動子として扱われる。
• テンプレート:仮リンクを利用すれば、複雑だが実行可能な解を得ることができる。
• オイラー法やルンゲ＝クッタ法のような様々な数値的手法が利用できる。

非減衰（$\delta =0$）かつ非駆動（$\gamma =0$）なダフィング方程式の特別な場合においては、テンプレート:仮リンクを利用することで厳密解を得ることができる。

非減衰かつ非強制振動子の解の有界性

非減衰かつ非強制（$\gamma =\delta =0$）なダフィング方程式に $\dot{x}$ を掛けると、次の式が得られる：

$\begin{array}{cc}& \dot{x}(\ddot{x}+\alpha x+\beta x^3)=0\\ & \Rightarrow \frac{d}{dt}[\frac{1}{2}{\left(\dot{x}\right)}^2+\frac{1}{2}\alpha x^2+\frac{1}{4}\beta x^4]=0\\ & \Rightarrow \frac{1}{2}{\left(\dot{x}\right)}^2+\frac{1}{2}\alpha x^2+\frac{1}{4}\beta x^4=H.\end{array}$

ここで テンプレート:Mvar は定数である。テンプレート:Mvar の値は初期条件 $x(0)$ および $\dot{x}(0)$ によって決まる。

テンプレート:Mvar に $y=\dot{x}$ を代入することで、システムはテンプレート:仮リンクであることが分かる：

$\dot{x}=+\frac{\partial H}{\partial y},$ テンプレート:Pad $\dot{y}=-\frac{\partial H}{\partial x}$ テンプレート:Pad with テンプレート:Pad $H=\frac{1}{2}{y}^2+\frac{1}{2}\alpha x^2+\frac{1}{4}\beta x^4.$

$\alpha$ と $\beta$ のいずれも正であるなら、解は有界である：

$|x|\le \sqrt{\frac{2H}{\alpha }}$ テンプレート:Pad and テンプレート:Pad $|\dot{x}|\le \sqrt{2H},$

ここでハミルトニアン テンプレート:Math は正である。

参考文献

その他

• テンプレート:Citation
• テンプレート:Citation

外部リンク

• テンプレート:Wayback - スカラーペディア百科事典「ダフィング方程式」の項目。
• テンプレート:MathWorld