Web読書会サイト「RETHAC」

TOP

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

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

本の検索

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

最新情報

新着記事順

新着書籍順

投稿数順

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

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

	P.373の問題番号「9.1.3」に対する解答水車の方程式 \[ \begin{align} \dot{a}_1 &= \omega b_1-Ka_1 \\ \dot{b}_1 &= -\omega a_1+q_1-Kb_1 \\ \dot{\omega} &= -\frac{\nu}{I}\omega + \frac{\pi gr}{I}a_1 \end{align} \] 解）変数変換 \[ x=\frac{\omega}{K}, \ \ y=\frac{\pi gr}{K \nu}a_1, \ \ z=-\frac{\pi gr}{K \nu}\left(b_1-\frac{q_1}{K} \right), \ \ \tau=Kt \]を行うと、ローレンツ方程式 \[ \begin{align} \dot{x} &= \sigma(y-x) \\ \dot{y} &= rx-xz-y \\ \dot{z} &= xy-bz \end{align} \]が得られる。ここで、 \[ \sigma = \frac{\nu}{KI}, \ \ r=\frac{\pi grq_1}{K^2 \nu}, \ \ b=1 \]と表すことができる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-06 04:40:37
	P.372の問題番号「9.1.2」に対する解答 (a) $Q(\theta)=q_1 \cos \theta$のとき、$n\neq 1$のモードは \[ \dot{a}_n=n \omega b_n-Ka_n, \ \ \dot{b}_n=-n\omega a_n-Kb_n \]となる。ここで、$a_n^2+b_n^2$の時間微分を考えると、 \[ \frac{d}{dt}(a_n^2+b_n^2)=2a_n \dot{a}_n+2b_n \dot{b}_n=-2K(a_n^2+b_n^2) \]となるので、 \[ a_n^2+b_n^2=Ae^{-2Kt} \]となり、$t \to \infty$で$a_n,b_n \to 0$となる。 (b) $Q(\theta)$が \[ Q(\theta)=\sum_{n=0}^{\infty} q_n \cos(n \theta) \]の場合、$n\neq 1$のモードは \[ \dot{a}_n=n \omega b_n-Ka_n, \ \ \dot{b}_n=-n\omega a_n-Kb_n+q_n \]となる。ここで、これらの方程式を \[ \begin{align} \dot{a}_n &=n \omega \left(b_n-\frac{Kq_n}{n^2\omega^2+K^2} \right)-K \left( a_n-\frac{n\omega q_n}{n^2\omega^2+K^2} \right), \\ \dot{b}_n &=-n\omega \left( a_n-\frac{n\omega q_n}{n^2\omega^2+K^2} \right)-K \left(b_n-\frac{Kq_n}{n^2\omega^2+K^2} \right) \end{align} \]のように変形することができる。このとき、$\omega$が時間依存するので、 \[ \frac{d}{dt} \left( a_n-\frac{n\omega q_n}{n^2\omega^2+K^2} \right) = \dot{a}_n, \ \ \frac{d}{dt}\left(b_n-\frac{Kq_n}{n^2\omega^2+K^2} \right)=\dot{b}_n \]にはならないが、ある程度時間が経過したときに、$\omega$が適当な固定点$\omega^$に近づいていくとすると、$\omega$を定数とみなせると考えられる。このとき、$t \to \infty$で、 \[ a_n \to \frac{n\omega^ q_n}{n^2{\omega^}^2+K^2}, \ \ b_n \to \frac{K q_n}{n^2{\omega^}^2+K^2} \]になると考えられる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-05 05:10:46
	P.372の問題番号「9.1.1」に対する解答 (a) 角度$\theta$での水の質量分布$m(\theta,t)$に対して、慣性モーメントは$m(\theta,t)r^2$となるので、 \[ I_{\mathrm{water}}=\int_0^{2\pi}m(\theta,t)r^2d \theta=r^2\int_0^{2\pi}m(\theta,t)d\theta=Mr^2 \]と表示できる。 (b) まず、$\dot{M}$は、 \[ \dot{M}=\int_0^{2\pi}\frac{\partial m(\theta,t)}{\partial t} d\theta \]と表せる。次に、連続の式(9.1.2)を用いると、 \[ \begin{align} \dot{M} &= \int_0^{2\pi} \left[ Q(\theta)-Km(\theta,t)-\omega \frac{\partial m}{\partial \theta} \right] d\theta \\ &= \int_0^{2\pi}Q(\theta)d\theta-K \int_0^{2\pi}m(\theta,t)d\theta-\omega \int_0^{2\pi} \frac{\partial m}{\partial \theta} d\theta \\ &= Q_{\mathrm{total}}-KM-\omega(m(2\pi,t)-m(0,t)) \\ &= Q_{\mathrm{total}}-KM \end{align} \]となることがわかる。 (c) $\dot{M}=Q_{\mathrm{total}}-KM$より$M$は時間経過するにつれて、安定な固定点$M^*=Q_{\mathrm{total}}/K$に近づくことがわかる。したがって、慣性モーメントは$t \to \infty$で \[ I(t) \to I_{\mathrm{wheel}}+\frac{r^2 Q_{\mathrm{total}}}{K} \]になることがわかる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-03 04:49:01
	P.324の問題番号「8.7.12」に対する解答 $N$個の同一の振動子系 \[ \dot{\theta}_i=f(\theta_i)+\frac{K}{N}\sum_{j=1}^Nf(\theta_j) \ \ (i=1,\cdots,N) \]$f(\theta)$は滑らかで$2\pi$周期的、すべての$\theta$について$f(\theta)>0$。解）同相解を$\theta_1(t)=\cdots=\theta_N(t)=\theta^(t)$とおくと、$N$個すべての方程式は \[ \frac{d\theta^(t)}{dt}=(1+K)f(\theta^) \]となり、同相解は周期的となる。次に、同相解の安定性を決定するため、$\eta_i(t)$を無限小の摂動として、$\theta_i(t)=\theta^(t)+\eta_i(t)$とすると、 \[\dot{\eta}_i(t)=f'(\theta^)\eta_i+Kf'(\theta^)\frac{1}{N}\sum_{j=1}^N\eta_j \]が得られる。ここで、$\eta$の2次の項を落とした。変数を \[ \mu=\frac{1}{N}\sum_{j=1}^N\eta_j, \ \ \xi_i=\eta_{i+1}-\eta_i \ \ (i=1,\cdots,N-1) \]と変換することによって分解できて、 \[ \dot{\mu}(t)=(1+KN)f'(\theta^)\mu(t), \ \ \dot{\xi}_i(t)=f'(\theta^)\xi_i(t) \ \ (i=1,\cdots,N-1) \]となる。これらの方程式を変数分離すると、 \[ \frac{d\mu}{\mu}=\frac{(1+KN)f'(\theta^)}{(1+K)f(\theta^)}d\theta^, \ \ \frac{d\xi_i}{\xi_i}=\frac{f'(\theta^)}{(1+K)f(\theta^)}d\theta^ \]となるので、閉軌道$\theta^$を1周した後の摂動の変化は \[ \ln\frac{\mu(T)}{\mu(0)}=\frac{1+KN}{1+K}[ \ln f(\theta^) ]_0^{2\pi}=0, \ \ \frac{\xi_i(T)}{\xi_i(0)}=\frac{1}{1+K}[ \ln f(\theta^*) ]_0^{2\pi}=0 \]となり、$\mu(T)=\mu(0), \ \ \xi_i(T)=\xi_i(0)$が得られる。よって、すべての$i$に対して$\eta_i(T)=\eta_i(0)$である。つまり、閉軌道を1周しても摂動はすべて変化しない。したがって、特性乗数はすべて$\lambda_j=1$となる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-01 11:33:09
	P.324の問題番号「8.7.11」に対する解答 $N$次元の系 \[ \dot{\phi}_i=\Omega+a \sin \phi_i + \frac{1}{N} \sum_{j=1}^N \sin \phi_i \]に対し、変数変換$\phi_i \to \phi_i+\frac{\pi}{2}$とすると、系は \[ \dot{\phi}_i=\Omega+a \cos \phi_i + \frac{1}{N} \sum_{j=1}^N \cos \phi_i \]と書き直すことができる。このとき、系は \[ t \to -t, \ \ \phi_i \to -\phi_i + 2\pi n \]の変換に対して可逆となる。同相解$\phi^(t)$の方程式 \[ \frac{d \phi^}{dt}=\Omega+(a+1) \cos \phi^* \]も$t \to -t, \ \ \phi^* \to -\phi^* + 2\pi n$に対して可逆となる。これは、同相解が$\phi^(t)=-\phi^(-t)+2\pi n$を満たすことを示す。特に、$t=0$のとき、$\phi^(0)=-\phi^(0)+2\pi n$となるので、 \[ \phi^(0) = \pi n \]となる。次に、$T\neq T'$として、$\phi^(t)$が \[ \phi^(T)=\pi n + 2\pi, \ \ \phi^(-T') = \pi n - 2\pi \]を満たすと仮定する。これは、閉軌道を1周するごとに帰還時間が変化する、つまり、同相状態は吸引的であることを示す。しかし、これは$\phi^(T)=-\phi^(-T)+2\pi n$であることから、 \[ \phi^(-T)=\phi^(-T')=\pi n - 2\pi \]となり、$T \neq T'$であることに反する。したがって、同相状態は吸引的ではないことが示される。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-31 12:47:40
	P.324の問題番号「8.7.10」に対する解答閉軌道をもつ系$\dot{\boldsymbol{x}}=\boldsymbol{f}(\boldsymbol{x})$のフロケ乗数の導出方法断面$S$上の点を$\boldsymbol{x}_0$とし、$\boldsymbol{x}(0)=\boldsymbol{x}_0$とすると、この系の一般解は \[ \boldsymbol{x}(t)=\boldsymbol{\varphi}(t,\boldsymbol{x}_0) \]と表すことができる。このとき、帰還時間を$T$とすると、ポアンカレ写像$P$は \[ P(\boldsymbol{x}_0)=\boldsymbol{\varphi}(T, \boldsymbol{x}_0) \]となる。この式から固定点$\boldsymbol{x}^$における線形化したポアンカレ写像は \[ DP(\boldsymbol{x}^)=\left. \frac{\partial \boldsymbol{\varphi}(T, \boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} \right\|_{\boldsymbol{x}_0=\boldsymbol{x}^} \]と書くことができる。一方、 \[ \frac{d}{dt} \frac{\partial \boldsymbol{\varphi}(t,\boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} = \frac{\partial \boldsymbol{f}(\boldsymbol{x})}{\partial \boldsymbol{x}} \frac{\partial \boldsymbol{\varphi}(t,\boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} \]が成り立つので、 \[ \left. \frac{\partial \boldsymbol{\varphi}(t,\boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} \right\|_{t=0} = \boldsymbol{I} \]を初期値として、この式を$0$から$T$まで数値積分すると、$DP(x^)$を得ることができる。ここで、$\boldsymbol{I}$は単位行列。最後に、$DP(x^*)$の特性方程式を解くことにより、数値的にフロケ乗数を求めることができる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-27 04:30:26
	P.323の問題番号「8.7.9」に対する解答系$\dot{r}=r-r^2, \ \ \dot{\theta}=1$ (a) $r_0$を$S$上の初期点とする。$\dot{\theta}=1$なので、最初に$S$に戻るのは帰還時間$t=2\pi$の後となる。よって$r_1=P(r_0)$で、$r_1$は \[ \int_{r_0}^{r_1} \frac{dr}{r(1-r)}=\int_0^{2\pi}dt=2\pi \]を満たす。この式より \[ r_1=(1-(1-r_0^{-1})e^{-2\pi})^{-1} \]となる。ゆえにポアンカレ写像$P$は \[ P(r)=(1-(1-r^{-1})e^{-2\pi})^{-1} \]となる。 (b) 固定点は$P$のグラフが45°の直線と交わる点$r^=1$に生じる。また、ポアンカレ写像$P(r)$の$r^=1$での傾きは$e^{-2\pi}<1$となるので、クモの巣図法により、固定点$r^=1$が安定で唯一であることがわかる。 (c) $r^$における線形化したポアンカレ写像は \[ DP(r^*)=e^{-2\pi} \]となる。つまり、特性乗数は$e^{-2\pi}$となる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 16:26:58
	P.323の問題番号「8.7.8」に対する解答系$\dot{x}+x=F(t), \ \ F(t)$は滑らかな$T$周期関数 $t=0$のとき$x(0)=x_0$とすると、一般解は \[ x(t)=\left[x_0+\int_0^tF(t')e^{t'}dt' \right] e^{-t} \]となる。この一般解よりポアンカレ写像$P$は \[ P(x_0)=x(T)=\left[x_0+\int_0^T F(t)e^tdt \right] e^{-T} \]となる。この$P$のグラフは傾きが$e^{-T}<1$の直線となるので、これは対角線とただ1つの点で交わる。したがって、クモの巣図法により、固定点はただ１つで、大域的に安定である。つまり、系は安定な$T$周期解$x(t)$を必然的にもつと考えられる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 10:49:32
	P.323の問題番号「8.7.7」に対する解答系$\dot{t}=1, \ \ \dot{\theta}=\sin t-\sin \theta$ 相図は添付図のようになる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 07:11:40
	P.323の問題番号「8.7.6」に対する解答まず、系$\dot{\theta} + \sin \theta = 0$を考えると、$\theta=2n\pi, \ \ (2n+1)\pi$に固定点をもち、それぞれ$\theta=2n\pi$が安定な固定点、$\theta=(2n+1)\pi$が不安定な固定点となる。これを考慮して、系$\dot{\theta} + \sin \theta = \sin t$を考える。最初、不安定な固定点付近にある$\theta$は$\sin t$の項がなければ、すぐに安定な固定点に向かうが、$\sin t$の効果により、不安定な固定点の周りでしばらく振動する。この後、適当な時刻で、$\theta$は安定な固定点の方に動き、$\sin t$の効果により、安定な固定点の周りで振動を続ける。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 07:09:57

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

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

スポンサードリンク

ログイン

こんにちは！ゲストさん

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

お知らせ

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管理人による読書の新習慣

11月に紹介された本

この本の読書会サイト

	P.373の問題番号「9.1.3」に対する解答水車の方程式 \[ \begin{align} \dot{a}_1 &= \omega b_1-Ka_1 \\ \dot{b}_1 &= -\omega a_1+q_1-Kb_1 \\ \dot{\omega} &= -\frac{\nu}{I}\omega + \frac{\pi gr}{I}a_1 \end{align} \] 解）変数変換 \[ x=\frac{\omega}{K}, \ \ y=\frac{\pi gr}{K \nu}a_1, \ \ z=-\frac{\pi gr}{K \nu}\left(b_1-\frac{q_1}{K} \right), \ \ \tau=Kt \]を行うと、ローレンツ方程式 \[ \begin{align} \dot{x} &= \sigma(y-x) \\ \dot{y} &= rx-xz-y \\ \dot{z} &= xy-bz \end{align} \]が得られる。ここで、 \[ \sigma = \frac{\nu}{KI}, \ \ r=\frac{\pi grq_1}{K^2 \nu}, \ \ b=1 \]と表すことができる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-06 04:40:37
	P.372の問題番号「9.1.2」に対する解答 (a) $Q(\theta)=q_1 \cos \theta$のとき、$n\neq 1$のモードは \[ \dot{a}_n=n \omega b_n-Ka_n, \ \ \dot{b}_n=-n\omega a_n-Kb_n \]となる。ここで、$a_n^2+b_n^2$の時間微分を考えると、 \[ \frac{d}{dt}(a_n^2+b_n^2)=2a_n \dot{a}_n+2b_n \dot{b}_n=-2K(a_n^2+b_n^2) \]となるので、 \[ a_n^2+b_n^2=Ae^{-2Kt} \]となり、$t \to \infty$で$a_n,b_n \to 0$となる。 (b) $Q(\theta)$が \[ Q(\theta)=\sum_{n=0}^{\infty} q_n \cos(n \theta) \]の場合、$n\neq 1$のモードは \[ \dot{a}_n=n \omega b_n-Ka_n, \ \ \dot{b}_n=-n\omega a_n-Kb_n+q_n \]となる。ここで、これらの方程式を \[ \begin{align} \dot{a}_n &=n \omega \left(b_n-\frac{Kq_n}{n^2\omega^2+K^2} \right)-K \left( a_n-\frac{n\omega q_n}{n^2\omega^2+K^2} \right), \\ \dot{b}_n &=-n\omega \left( a_n-\frac{n\omega q_n}{n^2\omega^2+K^2} \right)-K \left(b_n-\frac{Kq_n}{n^2\omega^2+K^2} \right) \end{align} \]のように変形することができる。このとき、$\omega$が時間依存するので、 \[ \frac{d}{dt} \left( a_n-\frac{n\omega q_n}{n^2\omega^2+K^2} \right) = \dot{a}_n, \ \ \frac{d}{dt}\left(b_n-\frac{Kq_n}{n^2\omega^2+K^2} \right)=\dot{b}_n \]にはならないが、ある程度時間が経過したときに、$\omega$が適当な固定点$\omega^$に近づいていくとすると、$\omega$を定数とみなせると考えられる。このとき、$t \to \infty$で、 \[ a_n \to \frac{n\omega^ q_n}{n^2{\omega^}^2+K^2}, \ \ b_n \to \frac{K q_n}{n^2{\omega^}^2+K^2} \]になると考えられる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-05 05:10:46
	P.372の問題番号「9.1.1」に対する解答 (a) 角度$\theta$での水の質量分布$m(\theta,t)$に対して、慣性モーメントは$m(\theta,t)r^2$となるので、 \[ I_{\mathrm{water}}=\int_0^{2\pi}m(\theta,t)r^2d \theta=r^2\int_0^{2\pi}m(\theta,t)d\theta=Mr^2 \]と表示できる。 (b) まず、$\dot{M}$は、 \[ \dot{M}=\int_0^{2\pi}\frac{\partial m(\theta,t)}{\partial t} d\theta \]と表せる。次に、連続の式(9.1.2)を用いると、 \[ \begin{align} \dot{M} &= \int_0^{2\pi} \left[ Q(\theta)-Km(\theta,t)-\omega \frac{\partial m}{\partial \theta} \right] d\theta \\ &= \int_0^{2\pi}Q(\theta)d\theta-K \int_0^{2\pi}m(\theta,t)d\theta-\omega \int_0^{2\pi} \frac{\partial m}{\partial \theta} d\theta \\ &= Q_{\mathrm{total}}-KM-\omega(m(2\pi,t)-m(0,t)) \\ &= Q_{\mathrm{total}}-KM \end{align} \]となることがわかる。 (c) $\dot{M}=Q_{\mathrm{total}}-KM$より$M$は時間経過するにつれて、安定な固定点$M^*=Q_{\mathrm{total}}/K$に近づくことがわかる。したがって、慣性モーメントは$t \to \infty$で \[ I(t) \to I_{\mathrm{wheel}}+\frac{r^2 Q_{\mathrm{total}}}{K} \]になることがわかる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-03 04:49:01
	P.324の問題番号「8.7.12」に対する解答 $N$個の同一の振動子系 \[ \dot{\theta}_i=f(\theta_i)+\frac{K}{N}\sum_{j=1}^Nf(\theta_j) \ \ (i=1,\cdots,N) \]$f(\theta)$は滑らかで$2\pi$周期的、すべての$\theta$について$f(\theta)>0$。解）同相解を$\theta_1(t)=\cdots=\theta_N(t)=\theta^(t)$とおくと、$N$個すべての方程式は \[ \frac{d\theta^(t)}{dt}=(1+K)f(\theta^) \]となり、同相解は周期的となる。次に、同相解の安定性を決定するため、$\eta_i(t)$を無限小の摂動として、$\theta_i(t)=\theta^(t)+\eta_i(t)$とすると、 \[\dot{\eta}_i(t)=f'(\theta^)\eta_i+Kf'(\theta^)\frac{1}{N}\sum_{j=1}^N\eta_j \]が得られる。ここで、$\eta$の2次の項を落とした。変数を \[ \mu=\frac{1}{N}\sum_{j=1}^N\eta_j, \ \ \xi_i=\eta_{i+1}-\eta_i \ \ (i=1,\cdots,N-1) \]と変換することによって分解できて、 \[ \dot{\mu}(t)=(1+KN)f'(\theta^)\mu(t), \ \ \dot{\xi}_i(t)=f'(\theta^)\xi_i(t) \ \ (i=1,\cdots,N-1) \]となる。これらの方程式を変数分離すると、 \[ \frac{d\mu}{\mu}=\frac{(1+KN)f'(\theta^)}{(1+K)f(\theta^)}d\theta^, \ \ \frac{d\xi_i}{\xi_i}=\frac{f'(\theta^)}{(1+K)f(\theta^)}d\theta^ \]となるので、閉軌道$\theta^$を1周した後の摂動の変化は \[ \ln\frac{\mu(T)}{\mu(0)}=\frac{1+KN}{1+K}[ \ln f(\theta^) ]_0^{2\pi}=0, \ \ \frac{\xi_i(T)}{\xi_i(0)}=\frac{1}{1+K}[ \ln f(\theta^*) ]_0^{2\pi}=0 \]となり、$\mu(T)=\mu(0), \ \ \xi_i(T)=\xi_i(0)$が得られる。よって、すべての$i$に対して$\eta_i(T)=\eta_i(0)$である。つまり、閉軌道を1周しても摂動はすべて変化しない。したがって、特性乗数はすべて$\lambda_j=1$となる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-08-01 11:33:09
	P.324の問題番号「8.7.11」に対する解答 $N$次元の系 \[ \dot{\phi}_i=\Omega+a \sin \phi_i + \frac{1}{N} \sum_{j=1}^N \sin \phi_i \]に対し、変数変換$\phi_i \to \phi_i+\frac{\pi}{2}$とすると、系は \[ \dot{\phi}_i=\Omega+a \cos \phi_i + \frac{1}{N} \sum_{j=1}^N \cos \phi_i \]と書き直すことができる。このとき、系は \[ t \to -t, \ \ \phi_i \to -\phi_i + 2\pi n \]の変換に対して可逆となる。同相解$\phi^(t)$の方程式 \[ \frac{d \phi^}{dt}=\Omega+(a+1) \cos \phi^* \]も$t \to -t, \ \ \phi^* \to -\phi^* + 2\pi n$に対して可逆となる。これは、同相解が$\phi^(t)=-\phi^(-t)+2\pi n$を満たすことを示す。特に、$t=0$のとき、$\phi^(0)=-\phi^(0)+2\pi n$となるので、 \[ \phi^(0) = \pi n \]となる。次に、$T\neq T'$として、$\phi^(t)$が \[ \phi^(T)=\pi n + 2\pi, \ \ \phi^(-T') = \pi n - 2\pi \]を満たすと仮定する。これは、閉軌道を1周するごとに帰還時間が変化する、つまり、同相状態は吸引的であることを示す。しかし、これは$\phi^(T)=-\phi^(-T)+2\pi n$であることから、 \[ \phi^(-T)=\phi^(-T')=\pi n - 2\pi \]となり、$T \neq T'$であることに反する。したがって、同相状態は吸引的ではないことが示される。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-31 12:47:40
	P.324の問題番号「8.7.10」に対する解答閉軌道をもつ系$\dot{\boldsymbol{x}}=\boldsymbol{f}(\boldsymbol{x})$のフロケ乗数の導出方法断面$S$上の点を$\boldsymbol{x}_0$とし、$\boldsymbol{x}(0)=\boldsymbol{x}_0$とすると、この系の一般解は \[ \boldsymbol{x}(t)=\boldsymbol{\varphi}(t,\boldsymbol{x}_0) \]と表すことができる。このとき、帰還時間を$T$とすると、ポアンカレ写像$P$は \[ P(\boldsymbol{x}_0)=\boldsymbol{\varphi}(T, \boldsymbol{x}_0) \]となる。この式から固定点$\boldsymbol{x}^$における線形化したポアンカレ写像は \[ DP(\boldsymbol{x}^)=\left. \frac{\partial \boldsymbol{\varphi}(T, \boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} \right\|_{\boldsymbol{x}_0=\boldsymbol{x}^} \]と書くことができる。一方、 \[ \frac{d}{dt} \frac{\partial \boldsymbol{\varphi}(t,\boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} = \frac{\partial \boldsymbol{f}(\boldsymbol{x})}{\partial \boldsymbol{x}} \frac{\partial \boldsymbol{\varphi}(t,\boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} \]が成り立つので、 \[ \left. \frac{\partial \boldsymbol{\varphi}(t,\boldsymbol{x}_0)}{\partial \boldsymbol{x}_0} \right\|_{t=0} = \boldsymbol{I} \]を初期値として、この式を$0$から$T$まで数値積分すると、$DP(x^)$を得ることができる。ここで、$\boldsymbol{I}$は単位行列。最後に、$DP(x^*)$の特性方程式を解くことにより、数値的にフロケ乗数を求めることができる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-27 04:30:26
	P.323の問題番号「8.7.9」に対する解答系$\dot{r}=r-r^2, \ \ \dot{\theta}=1$ (a) $r_0$を$S$上の初期点とする。$\dot{\theta}=1$なので、最初に$S$に戻るのは帰還時間$t=2\pi$の後となる。よって$r_1=P(r_0)$で、$r_1$は \[ \int_{r_0}^{r_1} \frac{dr}{r(1-r)}=\int_0^{2\pi}dt=2\pi \]を満たす。この式より \[ r_1=(1-(1-r_0^{-1})e^{-2\pi})^{-1} \]となる。ゆえにポアンカレ写像$P$は \[ P(r)=(1-(1-r^{-1})e^{-2\pi})^{-1} \]となる。 (b) 固定点は$P$のグラフが45°の直線と交わる点$r^=1$に生じる。また、ポアンカレ写像$P(r)$の$r^=1$での傾きは$e^{-2\pi}<1$となるので、クモの巣図法により、固定点$r^=1$が安定で唯一であることがわかる。 (c) $r^$における線形化したポアンカレ写像は \[ DP(r^*)=e^{-2\pi} \]となる。つまり、特性乗数は$e^{-2\pi}$となる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 16:26:58
	P.323の問題番号「8.7.8」に対する解答系$\dot{x}+x=F(t), \ \ F(t)$は滑らかな$T$周期関数 $t=0$のとき$x(0)=x_0$とすると、一般解は \[ x(t)=\left[x_0+\int_0^tF(t')e^{t'}dt' \right] e^{-t} \]となる。この一般解よりポアンカレ写像$P$は \[ P(x_0)=x(T)=\left[x_0+\int_0^T F(t)e^tdt \right] e^{-T} \]となる。この$P$のグラフは傾きが$e^{-T}<1$の直線となるので、これは対角線とただ1つの点で交わる。したがって、クモの巣図法により、固定点はただ１つで、大域的に安定である。つまり、系は安定な$T$周期解$x(t)$を必然的にもつと考えられる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 10:49:32
	P.323の問題番号「8.7.7」に対する解答系$\dot{t}=1, \ \ \dot{\theta}=\sin t-\sin \theta$ 相図は添付図のようになる。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 07:11:40
	P.323の問題番号「8.7.6」に対する解答まず、系$\dot{\theta} + \sin \theta = 0$を考えると、$\theta=2n\pi, \ \ (2n+1)\pi$に固定点をもち、それぞれ$\theta=2n\pi$が安定な固定点、$\theta=(2n+1)\pi$が不安定な固定点となる。これを考慮して、系$\dot{\theta} + \sin \theta = \sin t$を考える。最初、不安定な固定点付近にある$\theta$は$\sin t$の項がなければ、すぐに安定な固定点に向かうが、$\sin t$の効果により、不安定な固定点の周りでしばらく振動する。この後、適当な時刻で、$\theta$は安定な固定点の方に動き、$\sin t$の効果により、安定な固定点の周りで振動を続ける。ストロガッツ　非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-07-24 07:09:57