Scratch(プログラム言語)でルンゲクッタ法(4段4次)を実装した。まずはローレンツ方程式を計算。
Scratchで数値計算シリーズもようやくルンゲクッタまで来たよ。
うちのブログではおなじみ、8次(Dormand&Prince DOP853)や14次のルンゲクッタ法を実装しようかと思ったが、あまりにめんどくさいので断念。かわりにおなじみの4段4次のふつーーーのRunge-Kuttaを実装した。
こちら。
http://scratch.mit.edu/projects/11698142/
今回、初めてMake a blockという自分でブロックを作れる機能で、呼び出す関数を実装してみた。
最初に解くのは、、、もちろん、ローレンツ方程式。
dx/dt = -σ*(x-y)
dy/dt=-y-x*z+r*x
dz/dt=x*y-b*z
これだけいろんな人がプロジェクト作っていて、誰もルンゲクッタを実装していないというのは微分方程式ファンの私には驚き、、、
さて、これを応用していろいろやってみるか(続く)。
« 天神社でお参り。 | トップページ | 仮面ライダーのプロデューサー、平山亨さんがお亡くなりに、、、ご冥福をお祈りします。 »
「パソコン・インターネット」カテゴリの記事
- ユニクロとAkamaiのコラボTシャツに書かれているコードを解読してみる。base64でデコードするとbashのシェルスクリプトが出てきて実行すると♥PEACE♥FOR♥ALL♥FOR♥ALL♥PEACE♥FOR♥ALL♥という文字が色付きで正弦波で流れた。(2025.05.01)
- Google Gemini 2.5 Pro experimentalに高周波で使われるSパラメータのTouchstoneファイルを読み込んでプロットするC#コードを書いてもらうと570行のコードができて動いた。ファイルの拡張子snpのnでポート数を判別するが人間を信じないのでデータ数えて確認するのに笑った。(2025.04.21)
- Google ColabのJulia言語で搭載されているGeminiを使って一行もコードを書かずに2次元拡散方程式を差分法で計算してGIFアニメにする。次に同じように2次元波動方程式もやってもらう。(2025.04.09)
- Google ColabのJulia言語で主成分分析(PCA)をやってみる。データはおなじみアヤメ(iris)で、標準で特異値分解(SVD)が入っているのですぐできた。(2025.04.08)
- Google ColabのJulia言語でマンデルブロ集合、仏様のようなブッダブロ、燃える船・バーニングシップフラクタルを描いてみる。どれも計算が速い。(2025.04.04)
« 天神社でお参り。 | トップページ | 仮面ライダーのプロデューサー、平山亨さんがお亡くなりに、、、ご冥福をお祈りします。 »
コメント