バタフライ・エフェクトをGIFアニメで(PARI/GPでのロジスティック写像の計算)
今朝、新聞のテレビ欄見ていたらNHKドラマ”書店員ミチルの身の上話”の紹介で、バタフライ効果が・・・などと書かれていた。
テレビ欄でバタフライ効果、とかいうのが面白いなと思ってなにか計算してみようと思った。
バタフライ効果、というのは蝶の羽ばたきのような小さな空気の震えが、遠く離れた場所で大規模な異常気象を引き起こすようなことを言うもので、カオスの説明によくつかわれる。というか長期予報が当たらないたとえによくつかわれると言うか、、、
いちばん簡単なのはロジスティック写像の初期値をちょっと変えて、何回も繰り返すとどうなるか見てみるということだが、そもそも計算の誤差もあるので、普通のExcelとかそういうのは×。
そこで高精度で1000桁でも簡単に計算できるPARI/GPを使ってみる。
計算は
xn+1 = a * xn * (1-xn)
でa=3.9としておいて、x0として0.1, 0.1+0.01, 0.1+0.001, 0.1+0.0001, ...
とちょっとずつ初期値が違うものと比較してみた。
結果をGIFアニメにしたのがこちら。
50回繰り返しても重なるには、x0=0.1+1E-13くらいのずれじゃないとだめで、x0=0.11くらいならもう1回目からずれていっているという。
« 「かのこちゃんとマドレーヌ夫人」を読んだ。 | トップページ | 千里セルシーに行ってきた。 »
「学問・資格」カテゴリの記事
- JavaScriptの数値計算ライブラリmathjsを使う(11)バーニングシップフラクタルを描いてみる。このココログでも計算できるようにした。(2023.01.31)
- JavaScriptの数値計算ライブラリmathjsを使う(10) リーマンゼータ関数(Riemann Zeta function)を計算、3次元化してPlotlyでぐりぐり動かす。(2023.01.27)
- JavaScriptの数値計算ライブラリmathjsを使う(9) 仏陀のお姿のフラクタル Buddhabrot(ブッダブロ)を描く。このココログでも計算できるようにしてみた。(2023.01.26)
- JavaScriptの数値計算ライブラリmathjsを使う(8) 4段4次のルンゲクッタ法でローレンツ方程式を計算。このココログでもPlotlyで3次元でぐりぐり動かせるようにしてみた。 (2023.01.25)
- JavaScriptの数値計算ライブラリmathjsを使う(7) マンデルブロ集合を描く。このココログでも計算して図示できるようにしてみた。(2023.01.16)
コメント