超高精度でロジスティック写像を計算(PARI/GP)
カオスを生じる写像として非常によく知られているロジスティック写像
xn+1 = a xn ( 1- xn)
ですが、計算の桁数が少ないとすぐ信用されない結果しか出さないのでも有名。
で、PARI/GPで100000桁!!!精度で計算して
Excelと50桁計算できるサイトkeisan.casio.jpと比較してみた。
a=3.9, 初期値は0.1としておく。
誤差の結果がこれ。
Excelはもはや60くらいの繰り返しでアウト(倍精度実数のCでも同じことっす)。keisan.casio.jpは善戦するも。200超えたくらいでアウト。
ということで、こういうのの計算にPARI/GPはぴったりです。
« 御影近くのNTTのアンテナ | トップページ | 「アトポス」を読んだ。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年7月13日 Pythonの高周波ライブラリscikit-rfがv1.8.0に、SamsungがGalaxy Z Fold7など発表→QualcommがSnapdragon 8 Eliteが使われていると発表、NGMNが基地局アンテナの推奨事項をまとめる、STMicroとMetalenzがメタサーフェス光学のライセンス締結(2025.07.14)
- Interface2025年8月号Pythonで体験!はじめての暗号を買った。上杉暗号からRSA、AES、DHなど、特に楕円曲線暗号についてはコードも実際に動かすところまで詳しくかかれていた。耐量子暗号や聞いたことなかったY-00暗号や関数型暗号も記載。(2025.07.10)
- 高周波・RFニュース 2025年7月8日 NordicとSercommのセルラーIoTモジュール、iFixitがFairphone 6を分解、スコアは10/10、RCR wireless newsのウェビナー2件(6GとIndustry4.0)、SEMCOが高耐圧C0G MLCCを車載急速充電に提案(2025.07.09)
- 高周波・RFニュース 2025年7月2日 5G Americasが6Gに向けセンシングと通信ホワイトペーパー発行、KYOCERA AVXが3dBハイブリッドカプラ発表、TDKが車載薄膜インダクタ発表、Nordicが1次電池向けPMIC発表、ローデ・シュワルツの6GとAI/ML解説記事(2025.07.02)
- 高周波・RFニュース 2025年6月30日 QualcommがAIを用いた6Rxアンテナ解説、Next G Allianceと日本のXGMFが5G,6Gで協力、5G Americasが25Q1で5G加入者増加と発表、TechInsigtsがHuawei Pura 80 Pro+分解、Qorvoが5-7GHzのWi-Fi 7 FEM発表(2025.06.30)
コメント