自宅で参加できる読書会
NO. 00188748 DATE 2024 05 02

ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで の読書会ページ

ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806)

ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで

著者:Strogatz,StevenHenry,1959- 田中,久陽 中尾,裕也 千葉,逸人,1982-

出版社:丸善出版 (201501)

ISBN-10:4621085808

ISBN-13:9784621085806

意見、感想、コメントなど

投稿一覧に戻る

P.259の問題番号「7.6.13」への解答

ダフィン振動子$\ddot{x}+x+\varepsilon x^3=0, \ \ 0 < \varepsilon \ll 1, \ \ x(0)=a, \ \ \dot{x}(0)=0$

(a)ダフィン方程式をベクトル場として表すと、
$\dot{x}=y, \ \ \dot{y}=-x-\varepsilon x^3$
と書ける。また、この系のエネルギーは
\[ E= \frac{1}{2} \dot{x}^2 + \frac{1}{2} x^2 + \frac{1}{4} \varepsilon x^4 = \frac{1}{2} y^2 + \frac{1}{2} x^2 + \frac{1}{4} \varepsilon x^4 \]となり、初期条件$x(0)=a, \ \ \dot{x}(0)=0$を考慮すると、
\[ \frac{1}{2} y^2 + \frac{1}{2} x^2 + \frac{1}{4} \varepsilon x^4 = \frac{1}{2} a^2 + \frac{1}{4} \varepsilon a^4 \]が得られる。
次に、ベクトル場を極座標で
$x= r \cos \theta, \ \ y= r \sin \theta$
と表すと、エネルギーの式は
\[ \frac{1}{2} r^2 + \frac{1}{4} \varepsilon r^4 \cos^4 \theta = \frac{1}{2} a^2 + \frac{1}{4} \varepsilon a^4 \]となり、この式から、
\[ r^2 = \frac{-1 + \sqrt{ 1 + \varepsilon \cos^4 \theta ( 2a^2 + \varepsilon a^4 ) } }{\varepsilon \cos^4 \theta } \]を得る。また、
\[ \frac{d \theta}{d t} = \frac{x \dot{y} - y \dot{x} }{r^2} = -(1 + \varepsilon r^2 \cos^4 \theta ) \]となるので、周期$T(\varepsilon)$は
\[ \begin{align}
T(\varepsilon) &= \int^T dt = \int_0^{2 \pi} \frac{dt}{d \theta} d \theta \\
&= \int_0^{2 \pi} \frac{d \theta}{1 + \varepsilon r^2 \cos^4 \theta} \\
&= \int_0^{2 \pi} \frac{d \theta}{\sqrt{ 1 + \varepsilon \cos^4 \theta ( 2a^2 + \varepsilon a^4 ) }}
\end{align}
\]となる。

解答者:goodbook 解答日時:2021-03-02 04:30:37

コメントを書き込む

問題解答へのコメント

1

(b) 被積分関数を$\varepsilon$で展開すると、
\[ \begin{align}
(1 + \varepsilon \cos^4 \theta ( 2a^2 + \varepsilon a^4 ) )^{-\frac{1}{2}} &= 1 - \varepsilon a^2 \cos^4 \theta \\
&+ \varepsilon^2 a^4 \left[ - \frac{1}{2} \cos^4 \theta + \frac{3}{2} \cos^8 \theta \right] + \mathcal{O}(\varepsilon^3)
\end{align}
\]となるので、
\[ \int_0^{2 \pi} \cos^4 \theta d \theta = \frac{3}{8} \cdot 2 \pi, \ \ \int_0^{2 \pi} \cos^8 \theta d \theta = \frac{35}{128} \cdot 2 \pi \]を用いると、
\[ T(\varepsilon) = 2 \pi \left[ 1 - \frac{3}{8} \varepsilon a^2 + \frac{57}{256} \varepsilon^2 a^4 + \mathcal{O}(\varepsilon^3) \right] \]となる。すなわち、
\[ c_0 = 2 \pi, \ \ c_1 = -\frac{3 \pi a^2}{4}, \ \ c_2 = \frac{57 \pi a^4}{128} \]となる。最後に、
\[ \omega = \frac{2 \pi}{T} = \left[ 1 - \frac{3}{8} \varepsilon a^2 + \mathcal{O}(\varepsilon^2) \right]^{-1} = 1 + \frac{3}{8} \varepsilon a^2 + \mathcal{O}(\varepsilon^2) \]となり、これは(7.6.57)式と一致している。

投稿者:goodbook 投稿日時:2021-03-02 04:32:15