ロトカ・ヴォルテラ(Lotka–Volterra)方程式の数値解法

はじめに

ロトカ・ヴォルテラ系は数理生物学で現れ，動物種の成長をモデル化する．ここにが捕食動物の数，が餌食の数を表す2つの種があるとする．以下のロトカ・ヴォルテラ系

で，点（ドット）は時間についての微分を表す．

ロトカ・ヴォルテラ系(1)には不変量があり，すべてのについて一定である．

不変量(2)の等高線は閉じているので，解は周期解である．(3)の数値解も周期解となるのが望ましいが，必ずしもそうなるとは限らない．(4)はポアソン(Poisson)方程式系

である．このは(5)で定義されている．

ポアソン方程式系とポアソン積分については，[HLW02]のChapter VII.2と[MQ02]に記述されている．

ここで定義済みの問題を含むパッケージをロードして，ロトカ・ヴォルテラ方程式系を選ぶ：

解を視覚化する効用関数を定義する：

系(6)を解くのに前進（陽的）オイラー法を使う：

RadauIIAの陰的ルンゲ・クッタ(Runge-Kutta)で後退（陰的）オイラー法を定義し，それを使って(7) を解く．結果の螺旋は，初期条件から固定点(2,1)に向かって時計回りに渦を巻く：

ロトカ・ヴォルテラ方程式の不変量(8)を使った前進オイラー法の投影により，周期解が得られる：

正しい質的挙動を得るための別のアプローチとして，(9)を2つの系に加算的に分割するというものがある．

(10)を適切に解くと，ポアソン(Poisson)積分を構築することができる．

ロトカ・ヴォルテラ方程式の分割の方程式を定義する：

(11)の各系に後退および前進オイラー法を使って，分割法についてシンプレクティックオイラー法を定義する：

シンプレクティックオイラー法を使った数値解は，周期解である：

ここで，ロトカ・ヴォルテラ方程式を分割し，(12)の各系の流れ（厳密解）を計算してみる．解は次のようにして見付けることができる．ここで定数は各ステップの初期条件に関連していなければならない．

これは明示的であり非常に効率的な積分法を生成するという利点がある．"LocallyExact"メソッドは初期化段階だけでDSolveを使って，各分割の流れを計算する一般的な方法を提供する．

記号と数値のハイブリッドの分割法を設定し，これを使ってロトカ・ヴォルテラ系を解く：

分割法を使った数値解は，周期解である：

トップへ