可解カオスの例として、楕円関数（ヤコビのSn)が出てくるシュレーダー写像(Schröder map)を厳密解と精度保証付き102桁までの計算で比較（カシオの高精度計算サイト利用）
先日、ロジスティック写像のa=4の厳密解と精度保証付き桁数の多い計算を比較することをやって、結構盛り上がった。
カオスを生じるロジスティック写像を厳密解と精度保証付き桁数の多い計算（カシオの高精度計算サイトkeisan.casio.jp)で桁数102桁まで変えて比較してみる。
今回は別の厳密解がある例を計算してみよう。Schroder mapです。
こちらを参考に。
こんな写像で、かつヤコビの楕円関数Snを使った厳密解がある。
これはa=4のロジスティック写像の厳密解が三角関数でかけたのの拡張版で楕円関数が出てくる。
ではいつものようにカシオの高精度計算サイトkeisan.casio.jpでやってみる。
このサイトは精度保証付きで102桁まで計算できるのだ。
しかももちろん楕円関数も標準装備！なので簡単に試せる。ただ、、、初期値だけはどうしょうもないので桁数上げてニュートン法でξ0を計算してスタート。
結果はこちら。k=0.1の場合。
ロジスティック写像のときと同じく、いくら精度保証していても桁数が少ないともう100回も行かずにどんどんずれていく。
« 丸で三角で四角なもの、なーんだ？をGoogle検索の3次元グラフ機能で描いてみる（再掲） | トップページ | #仮面ライダービルド 第30話の話数を表す数式は、自分以下の数で互いに素なものが全部素数、で最大なものが30ということ(totative of nの最大が30） »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(11/21) IEEE Microwave Magazineにこれまでの無線通信についてのレビュー、Microwave Journalは6Gに向けたパワーアンプ、Bluetooth SIGが6GHz帯を検討、電波のバンドがめちゃくちゃややこしい件の解説。(2022.11.21)
- Raspberry PI 3を数年ぶりに触る。Mathematicaをインストールすると最新版の13.1が入った！新機能で分数階微分とその微分方程式が解けるようになってるので例題通りやってみた。(2022.11.10)
- 先延ばしの普遍的法則（課題提出、夏休みの宿題、進捗どうですか、などなど）N(t)=M/(D-t+C)-M/(D+C)を導出、カシオの高精度計算サイトkeisan.casio.jpにUP!締め切りギリギリじゃないとみんなやらない感覚にあってるな。(2022.11.05)
- Visual C# (C_sharp)の数値計算ライブラリ MathNET Numericsを使う(10) 数値積分としてガウス・クロンロッド積分公式と二重指数関数型積分公式を試す。 (2022.11.01)
- ハロウィンということで昔Scratchで作ったハロウィンのかぼちゃランタンの三体問題の8の字解を再掲。シンプレクティックオイラー法で計算してます。(2022.10.30)
「日記・コラム・つぶやき」カテゴリの記事
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（11/27更新) 日本も急増しているが、それ以上に中国の増え方がすごくなってる。韓国もじわじわ増えてきてる。(2022.11.27)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(11/21) IEEE Microwave Magazineにこれまでの無線通信についてのレビュー、Microwave Journalは6Gに向けたパワーアンプ、Bluetooth SIGが6GHz帯を検討、電波のバンドがめちゃくちゃややこしい件の解説。(2022.11.21)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（11/13更新) 日本は完全に増加に転じた。人口比18.4%まで到達。中国と韓国はもう公式サイトで陽性者数発表しなくなった。。。(2022.11.13)
- Raspberry PI 3を数年ぶりに触る。Mathematicaをインストールすると最新版の13.1が入った！新機能で分数階微分とその微分方程式が解けるようになってるので例題通りやってみた。(2022.11.10)
- 2022/11/8の皆既月食をiPhone 12 Pro Maxで撮影。だが手持ち固定なし設定一切合切デフォルトではこれが限界か…(2022.11.09)
トラックバック
この記事へのトラックバック一覧です： 可解カオスの例として、楕円関数（ヤコビのSn)が出てくるシュレーダー写像(Schröder map)を厳密解と精度保証付き102桁までの計算で比較（カシオの高精度計算サイト利用）:
« 丸で三角で四角なもの、なーんだ？をGoogle検索の3次元グラフ機能で描いてみる（再掲） | トップページ | #仮面ライダービルド 第30話の話数を表す数式は、自分以下の数で互いに素なものが全部素数、で最大なものが30ということ(totative of nの最大が30） »
コメント