アレンストーフ軌道を35段14次ルンゲクッタ法(適応刻み幅機能付き)で計算(PARI/GP)
高次ルンゲクッタを導入したのはいいのだが、何も考えず刻み幅を決めるととんでもなく時間がかかってしまう。もともとの公式は埋め込み型で、刻み幅が調整できるのにサボっていてその機能使っていなかった。そこでようやく今頃になってその機能入れてみた。解いたのは、
y1" = y1 + 2* y1' - μ'* (y1 + μ)/D1 - μ*(y1 - μ') / D2
y2" = y2 - 2* y2' - μ' *y2/D1 - μ*y2 / D2
D1 = ((y1+μ)^2+y2^2)^(3/2)
D2 = ((y1-μ')^2+y2^2)^(3/2)
μ=0.012277471
μ' = 1-μ
y1(0) = 0.994, y1'(0)=0,
y2(0)=0,y2'(0)=-2.00158510637908252240537862224
計算結果はこちら。同じ計算結果を得るのに前は9000ポイント必要だったけど、1800ポイントまで減少した。図は両方を重ねているけど、全く重なっている。
ソースはこちら。↓
« 「バガボンド33巻」を読んだ。「あきらめろ」というセリフに衝撃! | トップページ | van der Pol方程式を35段14次ルンゲクッタ法で計算(PARI/GP) »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス)を生むThomas-Rössler方程式のパラメータbを色々変えて、Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
« 「バガボンド33巻」を読んだ。「あきらめろ」というセリフに衝撃! | トップページ | van der Pol方程式を35段14次ルンゲクッタ法で計算(PARI/GP) »
コメント