ロトカ・ヴォルテラの競争方程式

ロトカ・ヴォルテラの競争方程式は、競争関係にある生物の個体数の変動を表す数理モデルの一種である寺本。2つの種AとBが食料や生息場所などを巡って競争するとき、種Aの存在は、種Bの繁殖を阻害する。種Aにとっても、種Bの存在は繁殖を邪魔する。このような、お互いの存在がそれぞれの個体数増殖に相互にマイナスの影響を与える点を含め、2種の個体数の時間変化を表したモデルがロトカ・ヴォルテラの競争方程式である。

名称はアルフレッド・ロトカとヴィト・ヴォルテラに由来する巌佐。競争方程式以外に競争式や競争系、競争モデル、単にロトカ・ヴォルテラ方程式などとも呼ばれる日本生態学会巌佐寺本。英語では Lotka-Volterra competition model や Lotka-Volterra equations などと呼ばれる寺本Silvertown。

1. 導出
-----

モデルを導出する過程には様々なものがあるが、典型的にはロトカ・ヴォルテラの競争方程式はロジスティック方程式の拡張として導入できるハーバーマン。簡単な状況から始めるために、2種ではなく、他の種と生物間相互作用を持たない1種のみが環境内に存在する場合を考える。種が1つであっても、各個体は限られた資源を競い合う。そのため、個体数が増えるほど個体数増加は落ち着いていき、ある程度の個体数で飽和すると考えられる。それを簡潔に表したのがロジスティック方程式で

という常微分方程式で表される巌佐。ここで、//t// は時間、//N// は個体数、//dN// は //N// の時間微分で個体数増加率を示す。//r// は内的自然増加率と呼ばれる、//K// は環境収容力と呼ばれる正の定数である。右辺の  が個体数が増えるほど個体数増加率を抑える効果を与えているハーバーマン。

ここにもう一種が加わり、元からいた種と限られた資源を巡って競争する状況が考えられる。それぞれの種を「種1」「種2」と呼ぶとする。種1の存在は、種2の繁殖を阻害する。種1にとっても、種2の存在は繁殖を邪魔する。このように、競争関係にある2つの種が存在するとき、お互いの存在がそれぞれの個体数増殖に相互にマイナスの影響を与えるマレー。この点を含め、2種の個体数の時間変化を表したモデルがロトカ・ヴォルテラの競争方程式である。

ロジスティック方程式では、自らの個体数増加による個体数増加率への抑制効果を  で表していた。同じように、他方の種の個体数が多いほど、一方の種の個体数増加率は抑制されると考えられる。そのため2種の場合は、種1に対する抑制の効果を  で表すことが考えられる大串。ここで、//N//1, //N//2 は種1と2のそれぞれの個体数である。//α//12 が種2の個体数が種1の個体数増加率に与える影響の大きさを与える正の定数である北原。このようなモデリングによって、種1と種2についてのロトカ・ヴォルテラの競争方程式は、以下のような2変数の連立常微分方程式で表される日本生態学会。

ここで、//r//1, //r//2 は種1と種2のそれぞれの内的自然増加率、//K//1, //K//2 は種1と2それぞれの環境収容力である。//α//12 が種2の個体数が種1の個体数増加率に与える影響の大きさを、//α//21 が種1の個体数が種2の個体数増加率に与える影響の大きさを表している。//α//12 と //α//21 は競争係数と呼ばれる大串日本生態学会。

同じモデルは

という形式でも表され、使われる重定ハーバーマン。//β//1, //β//2 は種内競争係数と呼ばれ、//γ//1, //γ//2 は種間競争係数あるいは単に競争係数と呼ばれる日本数理生物学会。上のモデルの係数との関係は //K//1 = //r//1, //K//2 = //r//2, //α//12 = //γ//12 , //α//21 = //γ//21 となっている寺本。よって、//α// の値は種間競争の強さと種内競争の強さの相対値ともいえる巌佐。

2. 解の振る舞い
---------

2.1. アイソクライン
------------

ロトカ・ヴォルテラの競争方程式の陽な解は求まっていない重定。アイソクライン法によって解の大域的な振る舞いを知ることができる北原。//N//1-//N//2平面上のアイソクラインは

または

を満たす曲線である。この条件を満たす曲線は

直線1：直線2：直線3：直線4：という4つの直線である北原。1番目は //N//2 軸と一致する直線である。2番目は //N//1 軸と一致する直線である。3番目は//N//1切片が //K//1、//N//2切片が //K//1 の直線である。4番目は//N//1切片が //K//2、//N//2切片が //K//2 の直線である。これらの直線を境界にして、それぞれの個体数増加率の正負が切り替わる。現実の生物個体数は正の値であるから、特に関係するのは3番目と4番目の直線である。3番目の直線上では、//dN//1 = 0 であるから、この直線を通る解は平面上を上下方向にだけ動く。そのため、このアイソクライン直線を傾き無限大のアイソクラインと呼ぶ巌佐。一方、4番目の直線上では、//dN//2 = 0 であるから、この直線を通る解は平面上を左右方向にだけ動く。そのため、このアイソクライン直線を傾き0のアイソクラインと呼ぶ巌佐。

2.2. 定性的な挙動
-----------

アイソクラインを使って、解の定性的な挙動を知ることができる。すなわち、//N//1-//N//2相平面上の点がアイソクライン直線の内側、外側、あるいは直線上にあるかどうかに注目すれば、個体数増加率の値そのものは不明でも、個体数増加率の正負は知ることができる日本生態学会。直線3の場合は

直線内側のとき： //dN//1 > 0
直線外側のとき： //dN//1 < 0
直線上のとき：//dN//1 = 0なので、内側では //N//1 は増加する方向、外側では //N//1 は減少する方向、直線上では //N//1 は増減しないことがわかる。直線4の場合も同様に、内側では //N//2 は増加する方向、外側では //N//2 は減少する方向、直線上では //N//2 は増減しないことがわかる日本生態学会巌佐。

直線3と直線4を重ね合わせることで、各個体数が変化する方向が判明する。//N//1-//N//2相平面上で直線3と直線4の相対的な位置関係は、係数の値によって次のように4種類ある重定。

　 //K//2 < //K//1 かつ //K//1 > //K//2 のとき　 //K//2 > //K//1 かつ //K//1 < //K//2 のとき　 //K//2 > //K//1 かつ //K//1 > //K//2 のとき　 //K//2 < //K//1 かつ //K//1 < //K//2 のときしたがって、それぞれの場合ごとに解の大局的な挙動が異なり、下図のようになる寺本大串。



2.3. 平衡点と最終結果
-------------

解の平衡点は //dN//1 = 0 かつ //dN//2 = 0 を満たす点であるマレー。それらは、

平衡点1：平衡点2：平衡点3：平衡点4：という4点である巌佐。これらの平衡点の安定性は、競争係数と環境収容力の値によって次のように決まる寺本。

　 //K//2 < //K//1 かつ //K//1 > //K//2 のとき：平衡点4は第1象限に存在せず、平衡点1, 3は不安定、平衡点2が大域的に安定　 //K//2 > //K//1 かつ //K//1 < //K//2 のとき：平衡点4は第1象限に存在せず、平衡点1, 2は不安定、平衡点3が大域的に安定　 //K//2 > //K//1 かつ //K//1 > //K//2 のとき：平衡点1, 4は不安定、平衡点2, 3が局所的に安定　 //K//2 < //K//1 かつ //K//1 < //K//2 のとき：平衡点1, 2, 3は不安定、平衡点4が大域的に安定これらの結果を具体的に言い換えると、十分な時間経過後、それぞれの個体数は次のような結果になるといえる。

　種2は絶滅して種1が残る。種1の個体数は //K//1 の値に収束する。　種1は絶滅して種2が残る。種2の個体数は //K//2 の値に収束する。　種2は絶滅して種1が残るか、種1は絶滅して種2が残るか、初期の個体数によってどちらかの結果となる。個体数は //K//1 または //K//2 の値に収束する。　種1と種2が共存して残る。それぞれの個体数は //K//1 − //α//12 //K//2 と //K//2 − //α//21 //K//1 の値に収束する。

3. 競争排除則
--------

モデルから得られた結果は、, ,  のケースでは片方の種はもう一方の種によって排除されるということであったマレー。この結果は多くの場合に共存は不可能ということを示唆している寺本。ロシアの生態学者ゲオルギー・ガウゼは、ロトカ・ヴォルテラの競争方程式が示す種間の排他性の存在をゾウリムシやイースト菌を用いた実験で確認した寺本。この結果をもとにガウゼは、同じニッチにある複数の種は平衡状態で長期的に共存できないという原則を提唱し、今日では競争排除則と呼ばれる日本生態学会寺本。

4. 名称
-----

名称については、競争方程式以外に、競争式や競争系、競争モデルなどとも呼ばれる日本生態学会巌佐寺本。表記揺れとしてヴォルテラではなくボルテラとも記される巌佐。他にロトカとヴォルテラの名を冠した生物の個体数のモデルとして、ロトカ・ヴォルテラの方程式がある。こちらは競争関係ではなく捕食・被食関係にある2種の個体数変動を表すモデルである。ただし、本記事で説明してきた競争関係を表すモデルについても単に「ロトカ・ヴォルテラの方程式」と呼ぶことも多い。

ロトカ・ヴォルテラの競争方程式はアメリカの数学者アルフレッド・ロトカとイタリアの数学者ヴィト・ヴォルテラに由来する。方程式の初出として、ロトカの1925年の著作、ヴォルテラの1926年あるいは1931年の論文・著作がしばし引用される大串。一方で、科学技術史研究者のシャロン・キングスランドは、ロトカは捕食・被食モデルの方には重点的に取り組んでいたものの、競争モデルへの関心は比較的薄く、競争モデル自体の論文を提出したのは1932年の1つだけであることを指摘している。それに対してヴォルテラは競争モデルに積極的に取り組み、ゲオルギー・ガウゼはヴォルテラの結論を前述のように実験的に検証した。そのためキングスランドは、ジョージ・イヴリン・ハッチンソンが呼んでいたように、Volterra-Gause equationsと呼ぶ方がより正確だろうと述べている。

5. 脚注
-----

6. 参照文献
-------

7. 外部リンク
--------