ルンゲクッタ14次でJapanese Attractorを描いてみる(PARI/GP)。
Duffing方程式のストロボプロットでジャパニーズ・アトラクタを描いてみよう。35段14次のルンゲクッタ法で、PARI/GPを使う(グラフはExcelで)。60桁精度で。
計算している式は
dx/dt = y
dy/dt = -k*y - x^3 + B*z
z=cos(t)
k=0.05, B=7
初期値はx(0)=2.1, y(0)=0で。
こういうカオスを生じるものはどんなに精度上げて計算してもいつかは破綻するわけですが、4段4次のルンゲクッタで倍精度よりは(このくらいの時間なら)圧倒的に正確に描けているはず。とはいえ当時、上田さんがこれを出したときは精度がどうのこうの、とか言って日本の当時の権威のような老人方からは黙殺されたそうですが(で海外で認められて逆輸入ということで、情けない話)。
« はねるのトびらのしりたしゲームコーナーの中近東風のBGMはアバのヴーレ・ヴー | トップページ | 誕生寺→音声館あたりをぶらぶら »
「学問・資格」カテゴリの記事
- 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)
- Gemini CLIが使えるようになっていたので早速VSCodeのターミナルから使って、JavaScriptで連立一次方程式を計算するコードを書いてもらった。普通にガウスの消去法で計算するhtmlを作ってくれた。(2025.06.27)
« はねるのトびらのしりたしゲームコーナーの中近東風のBGMはアバのヴーレ・ヴー | トップページ | 誕生寺→音声館あたりをぶらぶら »
コメント