2018/04/10僕はロトカヴォルテラ系が大好きなんですが、今日本を読んでいたら演習問題で出てきたので「ちょっくら復習するかw」とか思ってたら全然理解していないことに気づいて焦っていろいろ計算したのでそれをしたためておきます。ロトカヴォルテラ系は平衡点に関してヤコビ行列の固有値が純虚数になるということで微分方程式系の安定性解析で引き合いにだされることが多いと思います。大体の流れがじゃあ実際に2次元系で安定性解析をやってみよう→平衡点についてヤコビ行列の固有値が純虚数になってしまいましたね→これだけでは安定性についてはなにも言えません。いう感じで、もしかするとそのあとにリアプノフ関数が登場したりしてるかもしれません。定数項に関する記述が不完全(なにいってるかは以下読んでもらえればわかると思います)だったり、計算過程が省かれていたりして、リアプノフ関数の要件を満たしていることをきちんと確認している書籍やウェブサイトはあまり見ない気がします。いうことで計算してみました。興味の対象のロトカヴォルテラ系は以下のもの。dxy\end{align}平衡点を求めると次の2つ。今回興味があるのはもちろん後者です。いつも通り線形安定性解析をやってみることにします。\sqrt{ac}i\end{align}となっていつも通り純虚数が出てきます。天下り的にリアプノフ関数の候補として次の\(E(x,y)\)を考えることにします。\(x\)方向と\(y\)方向それぞれについて関数の形を調べると次のような増減表が得られます。0\end{align}つまり関数\(E(x,y)\)は平衡点で値0をとり、それ以外の点では正の値をとることになります。ここまででリアプノフ関数の一つ目の要件が満たせました。次に\(E(x,y)\)の軌道に沿った時間微分を計算します。0\end{align}これで2つ目の要件も満たせました。\right)\)は中立安定であるということがわかりました。定数項が適当というのは、\(E(x,y)\)の後ろについているあたかもなにかしらの帳尻合わせが働いているであろう項について単に\(K\)とかされていたりすることで、計算の過程が省かれているというのは、リアプノフ関数の一つ目の要件の確認がなされていないことを冒頭では指し示していました。リアプノフ関数そのものの導出や定数項の帳尻合わせについては後日追記したいと思います。その2力学系とそこからの本いろいろRメモ(行列の乗算、行列の生成、乱数)ロトカヴォルテラ系の安定性解析(2)コメントを残すコメントをキャンセルPREV:クレローの微分方程式NEXT:ロトカヴォルテラ系の安定性解析(2)プロフィール管理人の軟骨魚です。力学系とコーヒー牛乳が大好きな大学院生です。