オレゴネータ(Oregonator)をドルマン・プリンス法で計算(Excel VBA)
化学反応のモデルとして有名なオレゴネータをルンゲクッタ8次であるDormand&Prince法のコード、DOP853のVBA移植版で計算しよう。式はこちら。
\[ \dot y_{1} = 77.27 (y_{2} + y_{1}(1-8.375\times10^{-6}y_{1}-y_{2})) \] \[ \dot y_{2} = \frac{1}{77.27} (y_{3} - (1+y_{1})y_{2}) \] \[ \dot y_{3} = 0.161 (y_{1} - y_{3}) \]
常微分方程式の数値解法I,IIに紹介されていたパラメータを使った。
本に載っていたようなグラフが再現できた。
« 「神様のパズル」を読んだ。 | トップページ | 「二の悲劇」を読んだ。 »
「学問・資格」カテゴリの記事
- 高周波(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)
コメント