Mackey-Glass方程式(遅延微分方程式)をExcel VBAで計算。
久しぶりに遅延微分方程式を計算してみる。
まずは、マッキー・グラス方程式。私はこれを”フラクタル科学”で初めて見た。
http://www.scholarpedia.org/article/Mackey-Glass_equation
dx(t)/dt = β*x(t-τ) / (1+[x(t-τ)]^n) -γ* x(t)
という形。
これは血液中の血球濃度x(t)の変化を表したものということでした。
γの項は老化で失われること、またβの項は濃度が小さいと増え方は濃度に比例しているが、だんだん大きくなると急に増え方が少なくなる、ということを表している。過去τだけ前ってことはτかかって血球が作られる、ということ。
で計算方法は”常微分方程式の数値解法”に載っていたものを使う。
基本、ルンゲクッタであるが、過去のxの変わりにgをつかうというもの。
複雑そうだけど、実際にプログラムするのはルンゲクッタ4次なら簡単。
で計算結果(n=9.65、γ=1, β=2, τ=2)
« プーチンの引き分け、というのを聞いて”わいの人生は殺すか殺されるかだ。半殺しはない。”by 岩鬼正美を思い出した。 | トップページ | ”空を飛ぶイカ”をGIFアニメで。 »
「学問・資格」カテゴリの記事
- 高周波(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)
« プーチンの引き分け、というのを聞いて”わいの人生は殺すか殺されるかだ。半殺しはない。”by 岩鬼正美を思い出した。 | トップページ | ”空を飛ぶイカ”をGIFアニメで。 »
コメント