P.321の問題番号「8.6.3」 に対する解答

系$\dot{\theta}_1=\omega_1, \ \ \dot{\theta}_2=\omega_2$、$\omega_1/\omega_2$は無理数。
解）トーラス上の任意の点$(p_1,p_2)$、任意の初期条件$(q_1,q_2)$に対して、$p_1=0, \ \ p_2=q_2$としても一般性は変わらないので、このように設定する。
$\theta_2$方向に軌道が１周するときの時間は$2\pi/\omega_2$であるので、このとき、$\theta_1$方向には$2\pi \omega_1/\omega_2$だけ移動する。したがって、任意の初期条件$(q_1,q_2)$から初めて、$\theta_2$方向に軌道が$n$周したとすると、$\theta_1$方向には
\[ 2\pi \left[ \frac{\omega_1}{\omega_2} n+\delta_0 \right] \]に移動することになる。ここで、記号$[a]$を$a>0$に対してその整数部分を除いたものとして定義し、また、$q_1=2\pi \delta, \ \ 0<\delta_0<1$とおいた。
次に、自然数$n_0$を適切に選ぶと、
\[ \left[ \frac{\omega_1}{\omega_2}n_0+\delta_0 \right] < \min \left( \left[ \frac{\omega_1}{\omega_2} \right], \delta_0 \right) \]とすることができる。このときの$\theta_1$方向の位置を
\[ \delta_1=\left[ \frac{\omega_1}{\omega_2}n_0+\delta_0 \right] \]として、$2\pi \delta_1$とする。さらに、自然数$n_1$を適切に選ぶと、
\[ \left[ \frac{\omega_1}{\omega_2}n_1+\delta_1 \right] < \delta_1 \]とすることができる。このときの$\theta_1$方向の位置を
\[ \delta_2=\left[ \frac{\omega_1}{\omega_2}n_1+\delta_1 \right] \]として、$2\pi \delta_2$とする。
この操作を繰り返していくと、$N$回目の操作で、
\[ 2\pi \delta_N < \varepsilon \]とすることができる。
以上のことから、トーラス上の任意の点$q$から出発した軌道が$p$から距離$\varepsilon$以内を通るような
\[ t=\frac{2\pi}{\omega_2} \sum_{i=0}^N n_i < \infty \]が存在することがわかる。

ストロガッツ非線形ダイナミクスとカオス(9784621085806)

投稿者：goodbook 投稿日：2021-06-04 04:22:34

P.320の問題番号「8.6.2」 に対する解答

結合振動子系$\dot{\theta}_1 = \omega_1+K_1 \sin (\theta_2-\theta_1), \ \ \dot{\theta}_2 = \omega_2+K_2 \sin (\theta_1-\theta_2)$
(a) この系のヌルクラインは
\[ \sin (\theta_1-\theta_2)=\frac{\omega_1}{K_1}, \ \ \sin (\theta_1-\theta_2)=-\frac{\omega_2}{K_2} \]となるので、$\omega_1,\omega_2>0$および$K_1,K_2>0$ならば系は固定点をもたない。
(b) $\phi=\theta_1-\theta_2, \ \ \varphi=\theta_1+\theta_2$とおくと、
\[ \dot{\phi}=\omega_1-\omega_2-(K_1+K_2)\sin \phi, \ \ \dot{\varphi}=\omega_1+\omega_2-(K_1-K_2)\sin \phi \]の方程式が成り立つ。このとき、
\[ \frac{d \varphi}{d \phi}=\frac{\dot{\varphi}}{\dot{\phi}}=\frac{2\omega^*}{\omega_1-\omega_2-(K_1+K_2)\sin \phi}+\frac{K_1-K_2}{K_1+K_2} \]となる。ここで、
\[ \omega^*=\frac{K_1\omega_2+K_2\omega_1}{K_1+K_2} \]とおいた。この式より
\[ (K_1+K_2)d \varphi=\left[ (K_1-K_2)+\frac{2\omega^*}{s-\sin \phi} \right] d \phi \]と変形して両辺を積分した後、少し整理すると保存量を得ることができる。ここで、\[ s=\frac{\omega_1-\omega_2}{K_1+K_2} \]とおいた。
i) $|\omega_1-\omega_2|<K_1+K_2$のとき、
\[ E=\frac{s \tan(\phi/2)-1+\sqrt{1-s^2} \tanh [(K_1\theta_2+K_2\theta_1) \sqrt{1-s^2}/2\omega^*]}{\sqrt{1-s^2}+(s \tan(\phi/2)-1) \tanh [(K_1\theta_2+K_2\theta_1) \sqrt{1-s^2}/2\omega^*] } \] ii) $|\omega_1-\omega_2|>K_1+K_2$のとき、
\[ E=\frac{s \tan(\phi/2)-1-\sqrt{s^2-1} \tan [(K_1\theta_2+K_2\theta_1) \sqrt{s^2-1}/2\omega^*]}{\sqrt{s^2-1}+(s \tan(\phi/2)-1) \tan [(K_1\theta_2+K_2\theta_1) \sqrt{s^2-1}/2\omega^*] } \]となる。

ストロガッツ非線形ダイナミクスとカオス(9784621085806)

投稿者：goodbook 投稿日：2021-06-03 06:04:42

P.320の問題番号「8.5.3」 に対する解答

ロジスティック方程式$\dot{N}=rN(1-N/K(t))$ （$K(t)$は正で、tについて滑らかで$T$周期的）
(a) 曲線$N(t)=K(t)$を考える。このとき、$N(t)$の変化はこの曲線より上側では$\dot{N}<0$となり下向き、この曲線より下側では$\dot{N}>0$となり上向きとなる。特に、任意の$N>0$から始まる$N(t)$の値はやがて$K_{\mathrm{min}} \leq N \leq K_{\mathrm{max}} $の細い帯に入り、そこに永久に留まる。
時刻$t=0$で$N$であったものが、$T$時間後に$P(N)$となるポアンカレ写像を考える。$t=0, \ \ N=K_{\mathrm{min}}$で始まる$N$の時間変化を考えると、$P(K_{\mathrm{min}})>K_{\mathrm{min}}$となる。一方、$t=0, \ \ N=K_{\mathrm{max}}$で始まる$N$の時間変化を考えると、$P(K_{\mathrm{max}})<K_{\mathrm{max}}$となる。ポアンカレ写像$P(N)$が単調な連続関数であることを考慮すると、中間値の定理より$P(N)$のグラフは45°の対角線をどこかで横切らなくてはならない。この交点を$N^*$とすると、$P(N^*)=N^*$となり、これは周期$T$のリミットサイクルが存在することを示す。
(b) ポアンカレ写像$P(N)$がある区間で$P(N) \equiv N$となり、無限に多数の閉軌道の帯ができるとする。これらの閉軌道のうち、2つの閉軌道$N_1(t) > N_2(t)$を選らんでくる。この場合、ある時刻$t$において、$N_1(t) > K(t) > N_2(t)$となり、$\dot{N}_1<0. \ \ \dot{N}_2>0$となるはずである。このような$N_1, \ \ N_2$として互いに近接したもので、非常に短い時間間隔$\Delta t$後に、$N_1(t+\Delta t)=N_2(t+\Delta t)=K(t+\Delta t)$となるようなものを選ぶことができるはずである。しかしこれは軌道が互いに交わらないということに反する。したがって、リミットサイクルは唯一であると考えられる。

ストロガッツ　非線形ダイナミクスとカオス(9784621085806)

投稿者：goodbook 投稿日：2021-05-23 08:48:32

P.320の問題番号「8.5.2」 に対する解答

駆動された振動子$\phi''+\alpha \phi'+\sin \phi = I$
$y=\phi'$とすると、系は
\[ \phi'=y, \ \ y'=I-\sin \phi - \alpha y \]と表すことができる。

相図１）$\alpha$が一定で十分に小さければ、$I$を減少させることにより、系の安定なリミットサイクル(i)はサドルノード分岐を起こし（ii,iii)、双安定な状態（iii)になる。さらに$I$を減少させると、ホモクリニック分岐を起こし、リミットサイクルは破壊される（iv,v)。

ストロガッツ　非線形ダイナミクスとカオス(9784621085806)

投稿者：goodbook 投稿日：2021-05-22 08:43:23

P.320の問題番号「8.4.12」 に対する解答

点$(\mu,1)$を通過する軌道は安定多様体（ここでは$y$軸）から$\mathcal{O}(\mu)$の距離だけずれた軌道を進み、サドル（ここでは原点）に近づいていく。その後、軌道はサドルから逃げ出し、不安定多様体（ここでは$x$軸）から$\mathcal{O}(\mu)$の距離だけずれた軌道をサドルから離れる方向に進んでいくと考えられる。したがって、軌道がサドルから逃げ出して、$x(t) \approx 1$に達するまでにかかる時間は
\[ T = \int dt \approx \int_{\mu}^1 \frac{dt}{dx} dx = \int_{\mu}^1 \frac{1}{\lambda_u x} dx = \frac{1}{\lambda_u} \ln \left( \frac{1}{\mu} \right) \]と見積もられる。

ストロガッツ　非線形ダイナミクスとカオス(9784621085806)

投稿者：goodbook 投稿日：2021-05-17 05:43:44