Web読書会サイト「RETHAC」

TOP

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

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

本の検索

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

最新情報

新着記事順

新着書籍順

投稿数順

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

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

	P.316の問題番号「8.2.12」に対する解答 (a) 系$\dot{x}=-y+xy^2, \ \ \dot{y}=x-x^2$ 解）この系では$\omega=1, \ \ f(x,y)=xy^2, \ \ g(x,y)=-x^2$となる。原点では \[ f_x = y^2, \ \ f_y=2xy, \ \ g_x=-2x, \ \ g_y=0, \\ f_{xx}=0, \ \ f_{xy}=2y=0, \ \ f_{yy}=2x=0, \ \ g_{xx}=-2, \ \ g_{xy}=0, \ \ g_{yy}=0, \\ f_{xxx}=0, \ \ f_{xyy}=2, \ \ g_{xxy}=0, \ \ g_{yyy}=0 \]となるので、$a=1/8$が得られる。 (b) 系$\dot{x}=-y+ \mu x +xy^2, \ \ \dot{y}=x+ \mu y-x^2$は演習問題8.2.2より原点でホップ分岐を起こしている。また、(a)の結果からこの分岐は亜臨界であることがわかる。これは演習問題8.2.2-8.2.4の結果と一致する。ストロガッツ非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-04-26 06:25:54
	P.316の問題番号「8.2.11」に対する解答減衰を受けるダフィン振動子$\ddot{x}+\mu \dot{x} + x-x^3=0$ (a) $\dot{x}=y$とおくと、方程式は \[ \dot{x}=y, \ \ \dot{y}=-\mu y -x+x^3 \]と書くことができる。この系の固定点は$(0,0), \ \ (\pm 1, 0)$。以下、各固定点について分類。 \[ (\pm 1, 0) \ : \ J = \begin{pmatrix} 0 & 1 \\ 1 & -\mu \end{pmatrix}, \ \ \tau=-\mu, \ \ \Delta = -1<0 \]となるので、サドルとなる。 \[ (0, 0) \ : \ J = \begin{pmatrix} 0 & 1 \\ -1 & -\mu \end{pmatrix}, \ \ \tau=-\mu, \ \ \Delta = 1>0 \]となるので、$\mu>0$のとき、安定スパイラルとなり、$\mu$が$0$未満まで減少すると、不安定スパイラルに変わることがわかる。 (b) この系の相図を添付図に示す。相図より、$\mu=0$において、原点を囲む閉軌道の連続的な帯をもっているので、$\mu=0$における分岐がホップ分岐の退化版であることがわかる。ストロガッツ非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-04-25 17:30:18
	P.315の問題番号「8.2.10」に対する解答 \[ \dot{x}=B-x-\frac{xy}{1+qx^2}, \ \ \dot{y}=A-\frac{xy}{1+qx^2} \] 解）この系の固定点は \[ \left( B-A, \frac{A(1+q(B-A)^2)}{B-A} \right) \]となる。この固定点でのヤコビ行列は \[ J = \begin{pmatrix} -1+\frac{A(q(B-A)^2-1)}{(B-A)(1+q(B-A)^2)} & -\frac{B-A}{1+q(B-A)^2} \\ \frac{A(q(B-A)^2-1)}{(B-A)(1+q(B-A)^2)} & -\frac{B-A}{1+q(B-A)^2} \end{pmatrix}, \\ \tau=\frac{q(2A-B)(B-A)^2-B-(B-A)^2}{(B-A)(1+q(B-A)^2)}, \ \ \Delta = \frac{B-A}{1+q(B-A)^2} \]となる。ホップ分岐が起こる条件は$\tau=0, \ \ \Delta>0$であることから、 \[ q = \frac{B+(B-A)^2}{(2A-B)(B-A)^2}, \ \ B>A \]のとき、ホップ分岐が起こることがわかる。また、$q>0$であるので、$A,B$に対する条件は、 \[ A<B<2A \]となる。これらの式を満たす$A,B,q$をそれぞれ$A_c, B_c, q_c$とおく。添付図に$A_c=2, \ \ B_c=3, \ \ q_c=4$としたとき、$B<B_c$と$B>B_c$の場合の相図を示す。$B<B_c$のとき、$\tau<0$となり、固定点は安定スパイラルになる。一方、$B>B_c$のとき$\tau>0$となり、固定点は不安定スパイラルとなり、その周りに安定なリミットサイクルが現れる。したがって、この系は超臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-04-25 07:56:32
	P.315の問題番号「8.2.9」に対する解答 \[ \dot{x}=x(b-x-\frac{y}{1+x}), \ \ \dot{y}=y(\frac{x}{1+x}-ay) \] (a) この系のヌルクラインは$\dot{x}=0, \ \dot{y}=0$より、$x=0, \ \ y=(1+x)(b-x), \ \ y=0, \ \ y=\frac{x}{a(1+x)}$となる。これらを描くと、添付図の青いラインのようになる。 (b) すべての$a,b>0$において、ヌルクライン$y=(1+x)(b-x)$は点$(b,0)$と点$(0,b)$を通り、一方、$y=\frac{x}{a(1+x)}$は原点を通り、$y=1/a$に$x$が正の向きに漸近する関数になっているので、この2つのヌルクラインは必ず第１象限で交わる。つまり、正の固定点$x^>0, \ \ y^>0$が存在することがわかる。 (c) この系の固定点は$(0,0), \ \ (b,0), \ \ (x^,y^)$の3点。このうち、$(0,0)$は孤立していない固定点と予想され、$(b,0)$はサドルとなる。 \[ (x^,y^) \ : \ J = \frac{1}{1+x^} \begin{pmatrix} x^((b-1)-2x^) & -x^ \\ b-x^* & -x^* \end{pmatrix}, \\ \tau=\frac{x^((b-2)-2x^)}{1+x^}, \ \ \Delta = \frac{x^(b-bx^+2x^{2})}{(1+x^)^2} \]となる。ここで、$\tau=0$より$2x^=b-2$が得られ、このとき、 \[ \Delta=\frac{4(b-2)}{b^2}, \ \ \tau^2-4\Delta=-\frac{16(b-2)}{b^2} \]となるので、$b>2$のとき、ホップ分岐が起こることがわかる。また、このとき、 \[ a_c=\frac{x^}{(b-x^)(1+x^)^2}=\frac{4(b-2)}{b^2(b+2)} \]となる。 (d) まず、$a$と$x^$の関係を考える。$b$の値が固定されているとすると、ヌルクライン$y=(1+x)(b-x)$も固定される。このとき、$a$を大きくしていくと、ヌルクライン$y=\frac{x}{a(1+x)}$は$y=0$に近づいていくので、$x^$は大きくなりながら、$b$に近づいていくことがわかる。したがって、$a>a_c$のとき$\tau<0$、$a<a_c$のとき$\tau>0$となる。これを踏まえて、$b=3, \ \ a_c=4/45$のときの相図を描くと添付図のようになり、この系は$a>a_c$のとき$(x^,y^)$は安定スパイラルとなり、$a<a_c$のとき$(x^,y^*)$は不安定スパイラルで、その周りに安定なリミットサイクルを持つので、超臨界ホップ分岐を起こすことわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-23 05:31:43
	P.315の問題番号「8.2.8」に対する解答へのコメント系$\dot{x}=x[x(1-x)-y], \ \ \dot{y}=y(x-a)$ (a) ヌルクラインは添付図の青いライン。 (b) 固定点は$(0,0), \ \ (1,0), \ \ (a, a-a^2)$の3点。以下、各固定点ごとに分類。 \[ (0,0) \ : \ J = \begin{pmatrix} 0 & 0 \\ 0 & -a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = 0, \ \ \tau^2 - 4 \Delta = a^2 \geq 0 \]となるので、孤立していない固定点と予想される。 \[ (1,0) \ : \ J = \begin{pmatrix} -1 & -1 \\ 0 & 1-a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = a-1, \ \ \tau^2 - 4 \Delta = (a-2)^2 \geq 0 \]となるので、$0 \leq a<1$のときサドル、$a>1$のとき安定ノードとなる。 \[ (a,a-a^2) \ : \ J = \begin{pmatrix} a-2a^2 & -a \\ a-a^2 & 0 \end{pmatrix}, \ \ \tau=a-2a^2, \ \ \Delta = a^2(1-a), \ \ \tau^2 - 4 \Delta = a^2(4a^2-3) \]となるので、$0 \leq a < \frac{1}{2}$のとき不安定スパイラル、$\frac{1}{2}<a<1$のとき安定な固定点となり、$a>1$のときサドルとなる。また、固有値は \[ \lambda = \frac{a(1-2a) \pm a \sqrt{4a^2-3}}{2} \]となる。 (c) $a>1$における相図は添付図のようになる。ほぼすべての軌道は安定固定点$(1,0)$に向かう。つまり、捕食者は絶滅する。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-21 05:23:34
	P.315の問題番号「8.2.8」に対する解答系$\dot{x}=x[x(1-x)-y], \ \ \dot{y}=y(x-a)$ (a) ヌルクラインは添付図の青いライン。 (b) 固定点は$(0,0), \ \ (1,0), \ \ (a, a-a^2)$の3点。以下、各固定点ごとに分類。 \[ (0,0) \ : \ J = \begin{pmatrix} 0 & 0 \\ 0 & -a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = 0, \ \ \tau^2 - 4 \Delta = a^2 \geq 0 \]となるので、孤立していない固定点と予想される。 \[ (1,0) \ : \ J = \begin{pmatrix} -1 & -1 \\ 0 & 1-a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = a-1, \ \ \tau^2 - 4 \Delta = (a-2)^2 \geq 0 \]となるので、$0 \leq a<1$のときサドル、$a>1$のとき安定ノードとなる。 \[ (a,a-a^2) \ : \ J = \begin{pmatrix} a-2a^2 & -a \\ a-a^2 & 0 \end{pmatrix}, \ \ \tau=a-2a^2, \ \ \Delta = a^2(1-a), \ \ \tau^2 - 4 \Delta = a^2(4a^2-3) \]となるので、$0 \leq a < \frac{1}{2}$のとき不安定スパイラル、$\frac{1}{2}<a<1$のとき安定な固定点となり、$a>1$のときサドルとなる。また、固有値は \[ \lambda = \frac{a(1-2a) \pm a \sqrt{4a^2-3}}{2} \]となる。 (c) $a>1$における相図は添付図のようになる。ほぼすべての軌道は安定固定点$(1,0)$に向かう。つまり、捕食者は絶滅する。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-21 05:20:16
	P.315の問題番号「8.2.7」に対する解答系$\dot{x} = \mu x + y -x^2, \ \ \dot{y} = -x + \mu y + 2 x^2$ この系の原点でのヤコビ行列は \[ J = \begin{pmatrix} \mu & 1 \\ -1 & \mu \end{pmatrix}, \ \ \tau=2 \mu, \ \ \Delta = \mu^2+1>0, \ \ \tau^2 - 4 \Delta = -4 \]となり、その固有値は \[ \lambda = \mu \pm i \]となるので、この系は$\mu=0$のとき原点でホップ分岐を起こすことがわかる。この系の相図を描くと、添付図のようになる。 $\mu<0$のとき、原点は安定スパイラルで、$\mu>0$のとき原点は不安定スパイラルとなり、また原点周りに安定なリミットサイクルが現れるので、この系は超臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-18 16:12:33
	P.315の問題番号「8.2.6」に対する解答系$\dot{x} = \mu x + y -x^3, \ \ \dot{y} = -x + \mu y + 2 y^3$ この系の原点でのヤコビ行列は \[ J = \begin{pmatrix} \mu & 1 \\ -1 & \mu \end{pmatrix}, \ \ \tau=2 \mu, \ \ \Delta = \mu^2+1>0, \ \ \tau^2 - 4 \Delta = -4 \]となり、その固有値は \[ \lambda = \mu \pm i \]となるので、この系は$\mu=0$のとき原点でホップ分岐を起こすことがわかる。この系の相図を描くと、添付図のようになる。 $\mu<0$のとき、原点は安定スパイラルで、また原点周りに不安定なリミットサイクルが現れ、$\mu>0$のとき原点は不安定スパイラルとなり、不安定なリミットサイクルは消えるので、この系は亜臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-18 14:58:34
	P.314の問題番号「8.2.5」に対する解答系$\dot{x} = y+ \mu x, \ \ \dot{y} = -x + \mu y -x^2 y$ この系の原点でのヤコビ行列は \[ J = \begin{pmatrix} \mu & 1 \\ -1 & \mu \end{pmatrix}, \ \ \tau=2 \mu, \ \ \Delta = \mu^2+1>0, \ \ \tau^2 - 4 \Delta = -4 \]となり、その固有値は \[ \lambda = \mu \pm i \]となるので、この系は$\mu=0$のとき原点でホップ分岐を起こすことがわかる。この系の相図を描くと、添付図のようになる。 $\mu<0$のとき、原点は安定スパイラルで、$\mu>0$のとき原点は不安定スパイラルとなり、また原点周りに安定なリミットサイクルが現れるので、この系は超臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-18 08:35:29
	P.314の問題番号「8.2.4」に対する解答系$\dot{x}=-y+ \mu x + x y^2, \ \ \dot{y}=x+ \mu y - x^2$ (a) この系を極座標で表すと、 \[ \dot{r} = \mu r - \frac{r^2}{4} (\sin \theta + \sin 3 \theta) + \frac{r^3}{8} (1- \cos 4 \theta) \\ \dot{\theta} = 1 - \frac{r}{4} (3 \cos \theta + \cos 3 \theta) - \frac{r^2}{8} (2 \sin 2 \theta - \sin 4 \theta) \]となる。 (b) (a)の結果から、$r \ll 1$で振動的な項を無視すると、 \[ \dot{\theta} \approx 1, \ \ \dot{r} \approx \mu r + \frac{r^3}{8} + \cdots \]が得られる。 (c) $\mu=-0.005$のときの原点付近の軌道を描くと添付図のようになる。点$(0.3,0)$からはじまる軌道は原点から遠ざかる方向に向かっていることがわかる。一方、点$(0.15,0)$からはじまる軌道は原点に近づく方向に向かっていることがわかる。したがって、$\mu=-0.005<0$で半径$r \approx \sqrt{-8 \mu} = 0.2$の不安定リミットサイクルが存在することが示唆される。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-17 08:40:52

[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.316の問題番号「8.2.12」に対する解答 (a) 系$\dot{x}=-y+xy^2, \ \ \dot{y}=x-x^2$ 解）この系では$\omega=1, \ \ f(x,y)=xy^2, \ \ g(x,y)=-x^2$となる。原点では \[ f_x = y^2, \ \ f_y=2xy, \ \ g_x=-2x, \ \ g_y=0, \\ f_{xx}=0, \ \ f_{xy}=2y=0, \ \ f_{yy}=2x=0, \ \ g_{xx}=-2, \ \ g_{xy}=0, \ \ g_{yy}=0, \\ f_{xxx}=0, \ \ f_{xyy}=2, \ \ g_{xxy}=0, \ \ g_{yyy}=0 \]となるので、$a=1/8$が得られる。 (b) 系$\dot{x}=-y+ \mu x +xy^2, \ \ \dot{y}=x+ \mu y-x^2$は演習問題8.2.2より原点でホップ分岐を起こしている。また、(a)の結果からこの分岐は亜臨界であることがわかる。これは演習問題8.2.2-8.2.4の結果と一致する。ストロガッツ非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-04-26 06:25:54
	P.316の問題番号「8.2.11」に対する解答減衰を受けるダフィン振動子$\ddot{x}+\mu \dot{x} + x-x^3=0$ (a) $\dot{x}=y$とおくと、方程式は \[ \dot{x}=y, \ \ \dot{y}=-\mu y -x+x^3 \]と書くことができる。この系の固定点は$(0,0), \ \ (\pm 1, 0)$。以下、各固定点について分類。 \[ (\pm 1, 0) \ : \ J = \begin{pmatrix} 0 & 1 \\ 1 & -\mu \end{pmatrix}, \ \ \tau=-\mu, \ \ \Delta = -1<0 \]となるので、サドルとなる。 \[ (0, 0) \ : \ J = \begin{pmatrix} 0 & 1 \\ -1 & -\mu \end{pmatrix}, \ \ \tau=-\mu, \ \ \Delta = 1>0 \]となるので、$\mu>0$のとき、安定スパイラルとなり、$\mu$が$0$未満まで減少すると、不安定スパイラルに変わることがわかる。 (b) この系の相図を添付図に示す。相図より、$\mu=0$において、原点を囲む閉軌道の連続的な帯をもっているので、$\mu=0$における分岐がホップ分岐の退化版であることがわかる。ストロガッツ非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-04-25 17:30:18
	P.315の問題番号「8.2.10」に対する解答 \[ \dot{x}=B-x-\frac{xy}{1+qx^2}, \ \ \dot{y}=A-\frac{xy}{1+qx^2} \] 解）この系の固定点は \[ \left( B-A, \frac{A(1+q(B-A)^2)}{B-A} \right) \]となる。この固定点でのヤコビ行列は \[ J = \begin{pmatrix} -1+\frac{A(q(B-A)^2-1)}{(B-A)(1+q(B-A)^2)} & -\frac{B-A}{1+q(B-A)^2} \\ \frac{A(q(B-A)^2-1)}{(B-A)(1+q(B-A)^2)} & -\frac{B-A}{1+q(B-A)^2} \end{pmatrix}, \\ \tau=\frac{q(2A-B)(B-A)^2-B-(B-A)^2}{(B-A)(1+q(B-A)^2)}, \ \ \Delta = \frac{B-A}{1+q(B-A)^2} \]となる。ホップ分岐が起こる条件は$\tau=0, \ \ \Delta>0$であることから、 \[ q = \frac{B+(B-A)^2}{(2A-B)(B-A)^2}, \ \ B>A \]のとき、ホップ分岐が起こることがわかる。また、$q>0$であるので、$A,B$に対する条件は、 \[ A<B<2A \]となる。これらの式を満たす$A,B,q$をそれぞれ$A_c, B_c, q_c$とおく。添付図に$A_c=2, \ \ B_c=3, \ \ q_c=4$としたとき、$B<B_c$と$B>B_c$の場合の相図を示す。$B<B_c$のとき、$\tau<0$となり、固定点は安定スパイラルになる。一方、$B>B_c$のとき$\tau>0$となり、固定点は不安定スパイラルとなり、その周りに安定なリミットサイクルが現れる。したがって、この系は超臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス(9784621085806) 投稿者：goodbook 投稿日：2021-04-25 07:56:32
	P.315の問題番号「8.2.9」に対する解答 \[ \dot{x}=x(b-x-\frac{y}{1+x}), \ \ \dot{y}=y(\frac{x}{1+x}-ay) \] (a) この系のヌルクラインは$\dot{x}=0, \ \dot{y}=0$より、$x=0, \ \ y=(1+x)(b-x), \ \ y=0, \ \ y=\frac{x}{a(1+x)}$となる。これらを描くと、添付図の青いラインのようになる。 (b) すべての$a,b>0$において、ヌルクライン$y=(1+x)(b-x)$は点$(b,0)$と点$(0,b)$を通り、一方、$y=\frac{x}{a(1+x)}$は原点を通り、$y=1/a$に$x$が正の向きに漸近する関数になっているので、この2つのヌルクラインは必ず第１象限で交わる。つまり、正の固定点$x^>0, \ \ y^>0$が存在することがわかる。 (c) この系の固定点は$(0,0), \ \ (b,0), \ \ (x^,y^)$の3点。このうち、$(0,0)$は孤立していない固定点と予想され、$(b,0)$はサドルとなる。 \[ (x^,y^) \ : \ J = \frac{1}{1+x^} \begin{pmatrix} x^((b-1)-2x^) & -x^ \\ b-x^* & -x^* \end{pmatrix}, \\ \tau=\frac{x^((b-2)-2x^)}{1+x^}, \ \ \Delta = \frac{x^(b-bx^+2x^{2})}{(1+x^)^2} \]となる。ここで、$\tau=0$より$2x^=b-2$が得られ、このとき、 \[ \Delta=\frac{4(b-2)}{b^2}, \ \ \tau^2-4\Delta=-\frac{16(b-2)}{b^2} \]となるので、$b>2$のとき、ホップ分岐が起こることがわかる。また、このとき、 \[ a_c=\frac{x^}{(b-x^)(1+x^)^2}=\frac{4(b-2)}{b^2(b+2)} \]となる。 (d) まず、$a$と$x^$の関係を考える。$b$の値が固定されているとすると、ヌルクライン$y=(1+x)(b-x)$も固定される。このとき、$a$を大きくしていくと、ヌルクライン$y=\frac{x}{a(1+x)}$は$y=0$に近づいていくので、$x^$は大きくなりながら、$b$に近づいていくことがわかる。したがって、$a>a_c$のとき$\tau<0$、$a<a_c$のとき$\tau>0$となる。これを踏まえて、$b=3, \ \ a_c=4/45$のときの相図を描くと添付図のようになり、この系は$a>a_c$のとき$(x^,y^)$は安定スパイラルとなり、$a<a_c$のとき$(x^,y^*)$は不安定スパイラルで、その周りに安定なリミットサイクルを持つので、超臨界ホップ分岐を起こすことわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-23 05:31:43
	P.315の問題番号「8.2.8」に対する解答へのコメント系$\dot{x}=x[x(1-x)-y], \ \ \dot{y}=y(x-a)$ (a) ヌルクラインは添付図の青いライン。 (b) 固定点は$(0,0), \ \ (1,0), \ \ (a, a-a^2)$の3点。以下、各固定点ごとに分類。 \[ (0,0) \ : \ J = \begin{pmatrix} 0 & 0 \\ 0 & -a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = 0, \ \ \tau^2 - 4 \Delta = a^2 \geq 0 \]となるので、孤立していない固定点と予想される。 \[ (1,0) \ : \ J = \begin{pmatrix} -1 & -1 \\ 0 & 1-a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = a-1, \ \ \tau^2 - 4 \Delta = (a-2)^2 \geq 0 \]となるので、$0 \leq a<1$のときサドル、$a>1$のとき安定ノードとなる。 \[ (a,a-a^2) \ : \ J = \begin{pmatrix} a-2a^2 & -a \\ a-a^2 & 0 \end{pmatrix}, \ \ \tau=a-2a^2, \ \ \Delta = a^2(1-a), \ \ \tau^2 - 4 \Delta = a^2(4a^2-3) \]となるので、$0 \leq a < \frac{1}{2}$のとき不安定スパイラル、$\frac{1}{2}<a<1$のとき安定な固定点となり、$a>1$のときサドルとなる。また、固有値は \[ \lambda = \frac{a(1-2a) \pm a \sqrt{4a^2-3}}{2} \]となる。 (c) $a>1$における相図は添付図のようになる。ほぼすべての軌道は安定固定点$(1,0)$に向かう。つまり、捕食者は絶滅する。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-21 05:23:34
	P.315の問題番号「8.2.8」に対する解答系$\dot{x}=x[x(1-x)-y], \ \ \dot{y}=y(x-a)$ (a) ヌルクラインは添付図の青いライン。 (b) 固定点は$(0,0), \ \ (1,0), \ \ (a, a-a^2)$の3点。以下、各固定点ごとに分類。 \[ (0,0) \ : \ J = \begin{pmatrix} 0 & 0 \\ 0 & -a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = 0, \ \ \tau^2 - 4 \Delta = a^2 \geq 0 \]となるので、孤立していない固定点と予想される。 \[ (1,0) \ : \ J = \begin{pmatrix} -1 & -1 \\ 0 & 1-a \end{pmatrix}, \ \ \tau=-a, \ \ \Delta = a-1, \ \ \tau^2 - 4 \Delta = (a-2)^2 \geq 0 \]となるので、$0 \leq a<1$のときサドル、$a>1$のとき安定ノードとなる。 \[ (a,a-a^2) \ : \ J = \begin{pmatrix} a-2a^2 & -a \\ a-a^2 & 0 \end{pmatrix}, \ \ \tau=a-2a^2, \ \ \Delta = a^2(1-a), \ \ \tau^2 - 4 \Delta = a^2(4a^2-3) \]となるので、$0 \leq a < \frac{1}{2}$のとき不安定スパイラル、$\frac{1}{2}<a<1$のとき安定な固定点となり、$a>1$のときサドルとなる。また、固有値は \[ \lambda = \frac{a(1-2a) \pm a \sqrt{4a^2-3}}{2} \]となる。 (c) $a>1$における相図は添付図のようになる。ほぼすべての軌道は安定固定点$(1,0)$に向かう。つまり、捕食者は絶滅する。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-21 05:20:16
	P.315の問題番号「8.2.7」に対する解答系$\dot{x} = \mu x + y -x^2, \ \ \dot{y} = -x + \mu y + 2 x^2$ この系の原点でのヤコビ行列は \[ J = \begin{pmatrix} \mu & 1 \\ -1 & \mu \end{pmatrix}, \ \ \tau=2 \mu, \ \ \Delta = \mu^2+1>0, \ \ \tau^2 - 4 \Delta = -4 \]となり、その固有値は \[ \lambda = \mu \pm i \]となるので、この系は$\mu=0$のとき原点でホップ分岐を起こすことがわかる。この系の相図を描くと、添付図のようになる。 $\mu<0$のとき、原点は安定スパイラルで、$\mu>0$のとき原点は不安定スパイラルとなり、また原点周りに安定なリミットサイクルが現れるので、この系は超臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-18 16:12:33
	P.315の問題番号「8.2.6」に対する解答系$\dot{x} = \mu x + y -x^3, \ \ \dot{y} = -x + \mu y + 2 y^3$ この系の原点でのヤコビ行列は \[ J = \begin{pmatrix} \mu & 1 \\ -1 & \mu \end{pmatrix}, \ \ \tau=2 \mu, \ \ \Delta = \mu^2+1>0, \ \ \tau^2 - 4 \Delta = -4 \]となり、その固有値は \[ \lambda = \mu \pm i \]となるので、この系は$\mu=0$のとき原点でホップ分岐を起こすことがわかる。この系の相図を描くと、添付図のようになる。 $\mu<0$のとき、原点は安定スパイラルで、また原点周りに不安定なリミットサイクルが現れ、$\mu>0$のとき原点は不安定スパイラルとなり、不安定なリミットサイクルは消えるので、この系は亜臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-18 14:58:34
	P.314の問題番号「8.2.5」に対する解答系$\dot{x} = y+ \mu x, \ \ \dot{y} = -x + \mu y -x^2 y$ この系の原点でのヤコビ行列は \[ J = \begin{pmatrix} \mu & 1 \\ -1 & \mu \end{pmatrix}, \ \ \tau=2 \mu, \ \ \Delta = \mu^2+1>0, \ \ \tau^2 - 4 \Delta = -4 \]となり、その固有値は \[ \lambda = \mu \pm i \]となるので、この系は$\mu=0$のとき原点でホップ分岐を起こすことがわかる。この系の相図を描くと、添付図のようになる。 $\mu<0$のとき、原点は安定スパイラルで、$\mu>0$のとき原点は不安定スパイラルとなり、また原点周りに安定なリミットサイクルが現れるので、この系は超臨界ホップ分岐を起こしていることがわかる。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-18 08:35:29
	P.314の問題番号「8.2.4」に対する解答系$\dot{x}=-y+ \mu x + x y^2, \ \ \dot{y}=x+ \mu y - x^2$ (a) この系を極座標で表すと、 \[ \dot{r} = \mu r - \frac{r^2}{4} (\sin \theta + \sin 3 \theta) + \frac{r^3}{8} (1- \cos 4 \theta) \\ \dot{\theta} = 1 - \frac{r}{4} (3 \cos \theta + \cos 3 \theta) - \frac{r^2}{8} (2 \sin 2 \theta - \sin 4 \theta) \]となる。 (b) (a)の結果から、$r \ll 1$で振動的な項を無視すると、 \[ \dot{\theta} \approx 1, \ \ \dot{r} \approx \mu r + \frac{r^3}{8} + \cdots \]が得られる。 (c) $\mu=-0.005$のときの原点付近の軌道を描くと添付図のようになる。点$(0.3,0)$からはじまる軌道は原点から遠ざかる方向に向かっていることがわかる。一方、点$(0.15,0)$からはじまる軌道は原点に近づく方向に向かっていることがわかる。したがって、$\mu=-0.005<0$で半径$r \approx \sqrt{-8 \mu} = 0.2$の不安定リミットサイクルが存在することが示唆される。ストロガッツ非線形ダイナミクスとカオス : 数学的基礎から物理・生物・化学・工学への応用まで(9784621085806) 投稿者：goodbook 投稿日：2021-04-17 08:40:52