ロトカ・ヴォルテラ方程式を
dx/dt = (a-qy)x
dy/dt = (-b+rx)y
※xは被食者の個体数、yは捕食者の個体数であり、共に0より大きい数とする
(a,b,q,rは正定数)
の2式としたときに、不動点は(x,y) = (b/r,a/q)であり、アイソクラインによってxy平面を4つに分割(不動点(b/r,a/q)を中心として4つの領域にわかれる)すると、それぞれの領域でdx/dt、dy/dtの正負が反転するというところまではわかったのですが、肝心の解の挙動を求めることができません。
xy平面で表すと、不動点を中心とした楕円軌道を描くのだと予想しているのですが、根拠がありません。
可能性としてはx=b/rとy=a/qを漸近線とした曲線が描けるような気もします。
ここはやはり2式を微分方程式として解いてしまうのがてっとり早いのでしょうか?
よろしければご教授よろしくお願い致します。
A 回答 (5件)
- 最新から表示
- 回答順に表示
No.5
- 回答日時:
>肝心の解の挙動を求めることができません
軌道の式を求めることができたとしても、手元に計算機がないと軌道の様子は分かりません。解が安定かどうか、周期的かどうかを調べるには、リヤプノフ関数を見つけることができればいいんですが、そのリヤプノフ関数をどのように見つけるかが、これまた難しい問題です。
話のどこかで変数分離という言葉が出てきたようですから、試みに
L=P(x)+Q(y)
としてみてはどうですか。
No.3
- 回答日時:
で、
xとyの関係を求めて各定数に数値を入れてグラフを書いてみたの?
予想より実際の軌道が知りたいな
この回答への補足
すいません、先ほども記述しましたが、私の力量ではxとyの関係式を導き出すことはできません。変数分離を用いても困難でした。
どうやらguuman様の回答を見てるとやはり微分方程式を解いてみないことには始まらないみたいですね。
出直してきます。
No.2
- 回答日時:
辺々割ってdtを消去したら変数分離法で解けるのでは?
軌道はyとxの関係が分かればいいので軌道だけは求まるのでは?
時間変化は式から解析するしかないかも知れんが
補足にどんな軌道だったか書け
この回答への補足
回答ありがとうございます。
軌道はあくまで予想でしかないのですが、xy平面の第1象限をアイソクライン(直線x=b/r,y=a/q)によって4つの領域にわけ、直線x=b/rと直線y=a/qが交わる不動点を中心とした、反時計周りの楕円軌道になるのではないかと考えています。(前のお礼で書いたベクトル(dx/dt,dy/dt)の各四領域における正負から考えて)
No.1
- 回答日時:
軌道だけ分かればいいのならば微分法定式を解いてみたら
a,b,q,rを色々変えて図を書けいいんじゃないの
回答ありがとうございます。
説明が足りませんでした。
微分方程式を解くのが一番てっとり早いとは思っていたのですが、肝心の解法がわからないので、微分方程式を解かないで解の挙動を調べる方法はないものかと模索しておりました。
xy平面上で(dx/dt,dy/dt)のベクトルを考えたときに、4分割した右下の領域ではベクトルが右上を向くこと、右上の領域では左上を向くこと、左上の領域では左下と向くこと、左下の領域では右下を向くことのは確かなんですが、それが直線x=b/rとy=a/qにどのように接近、接触するのかが導き出せないのです。
感覚的な解釈になってしまいますが、捕食者と被食者の関係を考えると直線x=b/rとy=a/qは漸近線には成りえないとは思うのですが、それを数学的に解釈することができません。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・14歳の自分に衝撃の事実を告げてください
- ・架空の映画のネタバレレビュー
- ・「お昼の放送」の思い出
- ・昨日見た夢を教えて下さい
- ・ちょっと先の未来クイズ第4問
- ・【大喜利】【投稿~10/21(月)】買ったばかりの自転車を分解してひと言
- ・メモのコツを教えてください!
- ・CDの保有枚数を教えてください
- ・ホテルを選ぶとき、これだけは譲れない条件TOP3は?
- ・家・車以外で、人生で一番奮発した買い物
- ・人生最悪の忘れ物
- ・【コナン30周年】嘘でしょ!?と思った○○周年を教えて【ハルヒ20周年】
- ・ハマっている「お菓子」を教えて!
- ・最近、いつ泣きましたか?
- ・夏が終わったと感じる瞬間って、どんな時?
- ・10秒目をつむったら…
- ・人生のプチ美学を教えてください!!
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・都道府県穴埋めゲーム
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
積分で1/x^2 はどうなるのでし...
-
e^-2xの積分
-
フーリエ級数の問題で、f(x)は...
-
∫1/(x^2+1)^2 の不定積分がわ...
-
∫1/√x dx 積分せよ 教えて下さい
-
∫x^2√(4-x^2)dxの積分
-
項の右端につく縦棒の意味を教...
-
【数3 微分法の応用】 α=dv/dt ...
-
写真の赤丸のようになぜ、(d²y/...
-
積分 Xの-2乗を積分するとどう...
-
∫e^cos(x) dx の計算
-
x^2 * exp(x^2) dxの不定積分
-
フーリエ変換の問題について
-
微積分 dの意味
-
異常積分
-
媒介変数の積分ってなぜあのよ...
-
【数学Ⅱ・Ⅲ】微分の問題
-
x−1分の2の微分の仕方を教えて...
-
『偏微分』3年前から気になって...
-
ニュートン法(2次)よりも高次収...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
積分で1/x^2 はどうなるのでし...
-
e^-2xの積分
-
∫1/(x^2+1)^2 の不定積分がわ...
-
∫1/√x dx 積分せよ 教えて下さい
-
積分 Xの-2乗を積分するとどう...
-
1/X^2の積分ってlogX^2ですか?
-
項の右端につく縦棒の意味を教...
-
微積分 dの意味
-
∫e^cos(x) dx の計算
-
x/(a^2+x^2)の積分について
-
フーリエ級数の問題で、f(x)は...
-
数学 定積分の問題です
-
2次微分の変数変換
-
【数学Ⅱ・Ⅲ】微分の問題
-
確率密度関数をf(x)=1-|x-1|と...
-
定積分∫[3→0]|x^2-4|dxの答え...
-
∫x^2√(4-x^2)dxの積分
-
∮a^xdxこれを公式的に導いてほ...
-
媒介変数の積分ってなぜあのよ...
-
微分の記号
おすすめ情報