Web読書会サイト「RETHAC」

TOP

Web読書会サイト「RETHAC」とは

Web上で読書会を行うサイトです。
時間も場所も選ばずにいつでもどこからでも
読書感想や意見を交換することができます。

本の検索

これまで紹介された気になるフレーズ

最新情報

新着記事順

新着書籍順

投稿数順

キーワードを設定して記事を絞り込む

[1] << 前へ 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 次へ >> [48]

	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:24:19
	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.6.1」に対する解答系$\dot{\theta}_1 = \omega_1+\sin \theta_1 \cos \theta_2, \ \ \dot{\theta}_2 = \omega_2+\sin \theta_2 \cos \theta_1, \ \ \omega_1, \ \omega_2 \geq 0$ (a)(b) まず位相差$\phi=\theta_1-\theta_2$について調べる。このとき、 \[ \dot{\phi}=\dot{\theta}_1-\dot{\theta}_2=\omega_1-\omega_2+\sin \phi \]となる。したがって、 i) $\|\omega_1-\omega_2\|>1$ならば固定点はなし。 ii) $\|\omega_1-\omega_2\|=1$ならば1つの固定点、すなわち、半安定な周期解をもつ。 iii) $\|\omega_1-\omega_2\|<1$ならば2つの固定点、すなわち、安定な位相ロック解と不安定な位相ロック解をもつ。従って、系は$\|\omega_1-\omega_2\|=1$で周期軌道のサドルノード分岐を起こすことがわかる。一方、方程式は \[ \dot{\theta}_1= \omega_1+\frac{1}{2} \{ \sin (\theta_1+\theta_2)+\sin (\theta_1-\theta_2) \}, \\ \dot{\theta}_2= \omega_2+\frac{1}{2} \{ \sin (\theta_1+\theta_2)-\sin (\theta_1-\theta_2) \} \]と書ける。したがって、 \[ \sin ( \theta_1^+\theta_2^ ) = -\omega_1-\omega_2, \ \ \sin ( \theta_1^-\theta_2^ ) = -\omega_1+\omega_2 \]を満たす$(\theta_1^, \theta_2^)$を考えると、 i) $\omega_1+\omega_2>1$ならば固定点はなし。 ii) $\omega_1+\omega_2=1$ならば4つの固定点をもち、それぞれ半安定な固定点となる。 iii) $\omega_1+\omega_2<1$ならば$\|\omega_1+\omega_2\|=1$のときの4つの固定点がそれぞれサドルとノードに別れる。これは、$\omega_1+\omega_2=1$で無限周期分岐を起こすことを示している。以上の相図を添付(i)-(v)に示す。 (c) $(\omega_1, \omega_2)$平面上の安定性ダイアグラムを描くと添付図右下のようになる。 $\|\omega_1-\omega_2\|=1$を満たす直線上に周期軌道のサドルノード分岐をもち、$\omega_1+\omega_2=1$の直線上に無限周期分岐をもつ。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-30 13:07:44
	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:44:47
	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.5.1」に対する解答 $f(I)=(\ln I)^{-1}$とおくと、 \[ f^{(1)}=-\frac{1}{(\ln I)^2} \cdot \frac{1}{I} = -\frac{f^2}{I} \ \to \ If^{(1)}=-f^2 \]となる。この式の両辺を$I$で微分していくと、 \[ f^{(n+1)}=-I^{-1} \left( nf^{(n)}+\sum_{i=0}^n a_{i,n-i} f^{(i)} f^{(n-i)} \right) \]という関係式が得られる。ここで、係数$a_{i,n-i}$は \[ (1+x)^n = \sum_{i=0}^n a_{i,n-i} x^i \]を満たし、また$f^{(0)}=f$とした。 $\ln(I-I_c)$は$I_c$において発散するので、$f(I-I_c)$は$I_c$において$0$となる。次に、$(I-Ic)[ \ln (I-I_c)]^2$は$I_c$において$0$となるので、$f^{(1)}(I-I_c)$は$I_c$において発散する。したがって、上記で求めた$f^{(n+1)}$の関係式を用いると、$f^{(n+1)}(I-I_c)$は$I_c$において発散することがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-19 05:26:45
	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
	P.319の問題番号「8.4.11」に対する解答へのコメント \[ k=1, \ \ b=\frac{4}{3}, \ \ F=2 \](a) 平均化された系の相図は添付図のようになる。 (b) 演習問題8.4.10(a)で導入した$r_s, \ r_b$は数値的に求めることができ、これらの結果を使うと、$ a_s=2.723, \ \ a_b=4.063$と求めることができる。したがって、$a=2.8$のとき$a_s<2.8<a_b$を満たすので、演習問題8.4.11の結果から2つの安定な固定点をもつ。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-16 07:46:06
	P.319の問題番号「8.4.11」に対する解答 \[ k=1, \ \ b=\frac{4}{3}, \ \ F=2 \](a) 平均化された系の相図は添付図のようになる。 (b) 演習問題8.4.10(a)で導入した$r_s, \ r_b$は数値的に求めることができ、これらの結果を使うと、$ a_s=2.723, \ \ a_b=4.063$と求めることができる。したがって、$a=2.8$のとき$a_s<2.8<a_b$を満たすので、演習問題8.4.11の結果から2つの安定な固定点をもつ。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-16 07:44:11

[1] << 前へ 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 次へ >> [48]

--------------------

スポンサードリンク

ログイン

こんにちは！ゲストさん

[ログイン] [新規登録]（登録無料）

お知らせ

2022年7月に紹介された本(詳細)

(2022-08-01 04:18:31)

2022年6月に紹介された本(詳細)

(2022-07-01 04:52:59)

2022年5月に紹介された本(詳細)

(2022-06-01 04:03:28)

2022年4月に紹介された本(詳細)

(2022-05-01 06:38:57)

2022年3月に紹介された本(詳細)

(2022-04-01 04:34:01)

2022年2月に紹介された本(詳細)

(2022-03-01 04:57:24)

2022年1月に紹介された本(詳細)

(2022-02-01 04:17:11)

2021年11月に紹介された本(詳細)

(2021-12-01 05:33:59)

2021年5月に紹介された本(詳細)

(2021-06-01 05:02:45)

2021年1月に紹介された本(詳細)

(2021-02-01 05:25:32)

2020年12月に紹介された本(詳細)

(2021-01-01 08:54:16)

2020年11月に紹介された本(詳細)

(2020-12-01 07:51:26)

2020年10月に紹介された本(詳細)

(2020-11-02 05:11:05)

2020年9月に紹介された本(詳細)

(2020-10-02 05:12:01)

2020年8月に紹介された本(詳細)

(2020-09-01 05:05:05)

画像を添付できるようになりました。(詳細)

(2020-08-29 07:22:17)

2020年7月に紹介された本(詳細)

(2020-08-01 04:34:39)

2020年6月に紹介された本(詳細)

(2020-07-02 03:48:44)

2020年4月に紹介された本(詳細)

(2020-05-01 04:19:16)

管理者ブログ始めました(詳細)

(2020-04-11 09:07:40)

2020年3月に紹介された本(詳細)

(2020-04-01 04:25:30)

2020年1月に紹介された本(詳細)

(2020-02-02 06:45:14)

2019年12月に紹介された本(詳細)

(2020-01-04 06:57:23)

2019年11月に紹介された本(詳細)

(2019-12-01 05:02:20)

2019年10月に紹介された本(詳細)

(2019-11-01 05:52:58)

「ユーザー情報」：投稿履歴を書籍順で表示(詳細)

(2019-08-08 06:23:46)

2019年7月に紹介された本(詳細)

(2019-08-01 04:40:37)

2019年6月に紹介された本(詳細)

(2019-07-20 07:25:59)

（新機能２）twitterアカウントによるログイン(詳細)

(2019-06-10 05:51:27)

（新機能１）記事検索(詳細)

(2019-06-10 05:44:30)

サイト名の変更　読書会サイト「RETHAC」(詳細)

(2018-04-30 09:35:15)

コメント入力画面の変更について(詳細)

(2018-03-20 05:46:19)

お知らせ等の送信について(詳細)

(2018-03-18 14:27:29)

数式の利用について(詳細)

(2018-01-27 06:53:35)

気になるフレーズに本の画像を表示(詳細)

(2017-06-21 05:54:20)

いくつかの機能を追加、修正しました(詳細)

(2017-06-17 05:00:03)

Web読書会ページ新設(詳細)

(2017-04-29 20:10:40)

ブログをはじめました(詳細)

(2017-03-07 06:05:49)

Twitterアプリの連携機能を追加(詳細)

(2016-09-09 04:58:53)

レスポンシブWebデザインに対応(詳細)

(2016-02-17 06:16:03)

自動ツイート機能を追加(詳細)

(2015-04-10 06:36:25)

読書支援サイトRETHAC（α版）リリース(詳細)

(2015-01-04 04:39:59)

リンク

RETHAC管理人による読書の新習慣

7月に紹介された本

この本の読書会サイト

	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:24:19
	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.6.1」に対する解答系$\dot{\theta}_1 = \omega_1+\sin \theta_1 \cos \theta_2, \ \ \dot{\theta}_2 = \omega_2+\sin \theta_2 \cos \theta_1, \ \ \omega_1, \ \omega_2 \geq 0$ (a)(b) まず位相差$\phi=\theta_1-\theta_2$について調べる。このとき、 \[ \dot{\phi}=\dot{\theta}_1-\dot{\theta}_2=\omega_1-\omega_2+\sin \phi \]となる。したがって、 i) $\|\omega_1-\omega_2\|>1$ならば固定点はなし。 ii) $\|\omega_1-\omega_2\|=1$ならば1つの固定点、すなわち、半安定な周期解をもつ。 iii) $\|\omega_1-\omega_2\|<1$ならば2つの固定点、すなわち、安定な位相ロック解と不安定な位相ロック解をもつ。従って、系は$\|\omega_1-\omega_2\|=1$で周期軌道のサドルノード分岐を起こすことがわかる。一方、方程式は \[ \dot{\theta}_1= \omega_1+\frac{1}{2} \{ \sin (\theta_1+\theta_2)+\sin (\theta_1-\theta_2) \}, \\ \dot{\theta}_2= \omega_2+\frac{1}{2} \{ \sin (\theta_1+\theta_2)-\sin (\theta_1-\theta_2) \} \]と書ける。したがって、 \[ \sin ( \theta_1^+\theta_2^ ) = -\omega_1-\omega_2, \ \ \sin ( \theta_1^-\theta_2^ ) = -\omega_1+\omega_2 \]を満たす$(\theta_1^, \theta_2^)$を考えると、 i) $\omega_1+\omega_2>1$ならば固定点はなし。 ii) $\omega_1+\omega_2=1$ならば4つの固定点をもち、それぞれ半安定な固定点となる。 iii) $\omega_1+\omega_2<1$ならば$\|\omega_1+\omega_2\|=1$のときの4つの固定点がそれぞれサドルとノードに別れる。これは、$\omega_1+\omega_2=1$で無限周期分岐を起こすことを示している。以上の相図を添付(i)-(v)に示す。 (c) $(\omega_1, \omega_2)$平面上の安定性ダイアグラムを描くと添付図右下のようになる。 $\|\omega_1-\omega_2\|=1$を満たす直線上に周期軌道のサドルノード分岐をもち、$\omega_1+\omega_2=1$の直線上に無限周期分岐をもつ。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-30 13:07:44
	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:44:47
	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.5.1」に対する解答 $f(I)=(\ln I)^{-1}$とおくと、 \[ f^{(1)}=-\frac{1}{(\ln I)^2} \cdot \frac{1}{I} = -\frac{f^2}{I} \ \to \ If^{(1)}=-f^2 \]となる。この式の両辺を$I$で微分していくと、 \[ f^{(n+1)}=-I^{-1} \left( nf^{(n)}+\sum_{i=0}^n a_{i,n-i} f^{(i)} f^{(n-i)} \right) \]という関係式が得られる。ここで、係数$a_{i,n-i}$は \[ (1+x)^n = \sum_{i=0}^n a_{i,n-i} x^i \]を満たし、また$f^{(0)}=f$とした。 $\ln(I-I_c)$は$I_c$において発散するので、$f(I-I_c)$は$I_c$において$0$となる。次に、$(I-Ic)[ \ln (I-I_c)]^2$は$I_c$において$0$となるので、$f^{(1)}(I-I_c)$は$I_c$において発散する。したがって、上記で求めた$f^{(n+1)}$の関係式を用いると、$f^{(n+1)}(I-I_c)$は$I_c$において発散することがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-19 05:26:45
	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
	P.319の問題番号「8.4.11」に対する解答へのコメント \[ k=1, \ \ b=\frac{4}{3}, \ \ F=2 \](a) 平均化された系の相図は添付図のようになる。 (b) 演習問題8.4.10(a)で導入した$r_s, \ r_b$は数値的に求めることができ、これらの結果を使うと、$ a_s=2.723, \ \ a_b=4.063$と求めることができる。したがって、$a=2.8$のとき$a_s<2.8<a_b$を満たすので、演習問題8.4.11の結果から2つの安定な固定点をもつ。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-16 07:46:06
	P.319の問題番号「8.4.11」に対する解答 \[ k=1, \ \ b=\frac{4}{3}, \ \ F=2 \](a) 平均化された系の相図は添付図のようになる。 (b) 演習問題8.4.10(a)で導入した$r_s, \ r_b$は数値的に求めることができ、これらの結果を使うと、$ a_s=2.723, \ \ a_b=4.063$と求めることができる。したがって、$a=2.8$のとき$a_s<2.8<a_b$を満たすので、演習問題8.4.11の結果から2つの安定な固定点をもつ。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-05-16 07:44:11