ロトカ・ヴォルテラの方程式のソースを表示

{{Otheruses|[[捕食-被食関係]]のロトカ・ヴォルテラの方程式|[[競争 (生物)|競争関係]]|ロトカ・ヴォルテラの競争方程式|[[KEYTALK]]の楽曲、シングル|ロトカ・ヴォルテラ}}
[[Image:Lotka-Volterra.png|right|290px|thumb|ロトカ・ヴォルテラ方程式の解の一例。縦軸は個体数、横軸は時間。捕食者（Predatori、青）と被食者（Prede、赤）の個体数変動の[[位相]]は一般にずれており、捕食者が増加すると、急速に被食者が減少し、さらに捕食者が減少する、という時間変化を示す。]]
'''ロトカ・ヴォルテラの方程式'''（ロトカ・ヴォルテラのほうていしき、{{lang-en|Lotka-Volterra equations}}）とは、生物の[[捕食-被食関係]]による個体数の変動を表現する[[数理モデル]]の一種。2種の[[個体群]]が存在し、片方が[[捕食者]]、もう片方が被食者のとき、それぞれの個体数増殖速度を二元連立非線形[[常微分方程式]]系で表現する。'''ロトカ・ヴォルテラの捕食式'''や'''ロトカ・ヴォルテラ捕食系'''、'''ロトカ-ヴォルテラの捕食者-被食者モデル'''などとも呼ばれる{{Sfn|日本生態学会（編）|2004|p=141}}{{Sfn|日本生態学会（編）|2015|p=44}}<ref>{{cite book|和書 |author=Steven H. Strogatz |translator=田中久陽・中尾裕也・千葉逸人 |title=ストロガッツ　非線形ダイナミクスとカオス―数学的基礎から物理・生物・化学・工学への応用まで |publisher=丸善出版 |year=2015 |isbn=978-4-621-08580-6 |page=208 }}</ref>。

具体的には以下の方程式で表される{{Sfn|巌佐|1990|p=35}}。

:<math>\begin{align}
\frac{dx}{dt} &= a x - b xy,\\
\frac{dy}{dt} &= c xy - d y.
\end{align}</math>

ここで ''x'' は被食者の個体数、 ''y'' は捕食者の個体数、''t'' は時間をあらわし、4つの係数 ''a'', ''b'', ''c'', ''d'' は正の[[実数]]のパラメータである。

被食者と捕食者の個体数変動パターンの一つの例として、被食者が自然増殖して増えていくとそれを餌とする捕食者も増殖し、捕食者が増殖したことによって被食頻度が増えて被食者が減少し、被食者が減少したことによってそれを餌とする捕食者も減少し、捕食者が減少したことによって被食者の自然増殖数が被食頻度を上回って被食者が増え、そして最初に戻り…、このような形で被食者と捕食者が交互に増減し続けることが考えられる{{Sfn|日本生態学会（編）|2004|p=141}}{{Sfn|日本生態学会（編）|2015|pp=40&ndash;41}}。ロトカ・ヴォルテラの方程式は、このような個体数の周期的な増減の様子を示すことができる簡素で基礎的なモデルとなっている{{Sfn|マレー|2014|p=71}}。

名称は、この方程式をそれぞれ独立発案したアメリカの数学者[[アルフレッド・ロトカ]]とイタリアの数学者[[ヴィト・ヴォルテラ]]に由来する{{Sfn|Berryman|1992|p=1531}}。ロトカは1910年に化学物質濃度の変動を説明するために<ref>Lotka, A.J., "Contribution to the Theory of Periodic  Reaction", ''Journal of Physical Chemistry A|J. Phys. Chem.'', '''14'''  (3), pp 271–274  (1910)</ref><ref name="Goelmany">Goel, N.S. et al., “On the Volterra and Other Non-Linear Models of  Interacting Populations”, ''Academic Press Inc.'',  (1971)</ref>、ヴォルテラは1926年に[[アドリア海]]の魚数の変動を説明するために発案した{{Sfn|マレー|2014|pp=65&ndash;66}}。

== 式の導出と前提条件 ==
===被食者の増殖速度===
[[File:Tiger chasing boar.jpg|thumb|[[トラ]]から逃げる[[イノシシ]]]]
モデルの連立方程式内の
:<math>\frac{dx}{dt} = a x - b xy</math>
は被食者の個体数増殖速度 ''dx''/''dt'' を表している。上記の式は、以下のような生態学的な前提条件から導出される。

まず、捕食者が存在しない場合を仮定すると、被食者の個体数 ''x'' は順調に自然増していくと考えられる。この自然増は、[[マルサスモデル]]のようにその個体数に比例して増殖速度が増え、制限なく指数関数的に増殖すると仮定する{{Sfn|マレー|2014|p=65}}。すなわち、被食者にとっての餌は不足することなく十分あるような環境にあると仮定する{{Sfn|ハーバーマン|1992|p=108}}。これを表しているのが、右辺第一項 ''ax'' である{{Sfn|日本生態学会（編）|2015|p=42}}。

しかし、捕食者が存在する場合、被食者の個体数は捕食によって減少し、捕食者の存在は被食者増殖速度を抑制する効果を持つ。よって、捕食者数 ''y'' に比例して被食者増殖速度 ''dx''/''dt'' が減少すると仮定できる{{Sfn|寺本|1997|p=25}}。またさらに、捕食者がランダムに被食者を探索しているとすれば、被食者個体数が多いほど出会う割合が高まると考えられる{{Sfn|巌佐|1990|p=35}}。よって、被食者増殖速度は被食者個体数にも比例して減少すると仮定できる{{Sfn|伊藤|1994|p=80}}。これを表しているのが、右辺第二項 &minus;''bxy'' である{{Sfn|日本生態学会（編）|2015|p=42}}。このような、それぞれの個体数の単純な積で個体数増殖速度への影響を表すことを、[[質量作用の法則]]や質量作用の仮定と呼ぶ<ref>{{cite book |和書 |others=瀬野裕美（責任編集）|editor=日本数理生物学会 |title=「数」の数理生物学 |publisher=共立出版 |series=シリーズ 数理生物学要論 巻1 |year=2008 |edition=初版 |isbn=978-4-320-05675-6 |page=9}}</ref>。ロトカ・ヴォルテラの方程式は、この原則を基礎としている{{Sfn|Berryman|1992|p=1534}}。

===捕食者の増殖速度===
[[File:Bobcat having caught a rabbit.jpg|thumb|[[ウサギ]]を捕食する[[ボブキャット]]]]
捕食者の個体数増殖速度 ''dy''/''dt'' は
:<math>\frac{dy}{dt} = c xy - d y</math>
と表される。上記の式は、以下のような生態学的な前提条件から導出される。

まず、被食者が存在しない場合を考える。被食者にとっての餌はこの方程式系に現れる変数とは別に常に十分あると仮定したが、捕食者にとっての餌は被食者のみとする{{Sfn|ハーバーマン|1992|pp=108&ndash;109}}。よって、被食者が存在しないことは食糧が尽きたことと同じであり、捕食者の死亡率は出産率を上回り、捕食者の個体数 ''y'' は減少の一途を辿ることになる。この減少の仕方も、被食者の自然増のように個体数が多ければ多いほど減少速度が大きくなる、すなわち個体数 ''y'' に減少速度 ''dy''/''dt'' が比例すると仮定する{{Sfn|マレー|2014|p=65}}。これを表しているのが、右辺第二項 &minus;''dy'' である{{Sfn|日本生態学会（編）|2015|p=42}}。

そして、捕食者が増える速度は、捕食に成功した回数に比例すると考えられる{{Sfn|日本生態学会（編）|2015|p=43}}。捕食による被食者減少速度が &minus;''bxy'' と仮定されたように、捕食による捕食者増殖速度も同じ理屈から被食者数 ''x'' と捕食者数 ''y'' に比例するといえる。これを表しているのが、右辺第一項 ''cxy'' である{{Sfn|伊藤|1994|p=80}}。

==個体数の振る舞い ==
このロトカ・ヴォルテラ方程式を解析的に解いて ''x'' と ''y'' の ''t'' に関する明示的な解を得ることはできない{{Sfn|大串|1994|p=71}}。しかし、以下のような解の挙動を分析し、それぞれの個体数がどのように振る舞うかを知ることができる。

===平衡点===
[[File:Lotka-Volterra equations equilibrium points.png|thumb|ロトカ・ヴォルテラ方程式における2つの平衡点]]
どのようなときに、個体数 ''x'', ''y'' が増えも減りもしない、つまり時間 ''t'' の経過によらず変化しない状態になるかについて考える。これは、方程式の ''dx''/''dt'' と ''dy''/''dt'' がともに 0 ということなので、次式が得られる。
:<math>x (a - b y) = 0</math>
:<math>y (c x - d) = 0</math>
この式を満たす ''x'', ''y'' の組合せは、次の2組である{{Sfn|ハーバーマン|1992|p=112}}。
:<math>(x=0,\quad y=0)</math> 
:<math>\left(x=\frac{d}{c},\quad y=\frac{a}{b}\right)</math>
''x'', ''y'' がこれら2組の値をとるとき、その ''x'', ''y'' の値は時間によらず一定となる。このような点を[[平衡点]]と呼ぶ{{Sfn|寺本|1997|p=77}}。平衡点 (0, 0) は、捕食者も被食者も全滅してしまった状態である{{Sfn|ハーバーマン|1992|p=112}}。一方、平衡点 (''d''/''c'', ''a''/''b'') では、捕食者・被食者ともにある個体数で共存する状態となっている{{Sfn|日本生態学会（編）|2015|p=43}}。

これらの平衡点から ''x'', ''y'' の状態点がわずかにずれて与えられるときに、状態点が時間にともなって平衡点に収束するのか、それとも離れていくのかを特徴づける[[安定性理論|安定性]]は、次のように判別できる。2次以上の項が無視できるほどズレが小さいとすれば、平衡点 (0, 0) 近傍で系は次のように表すことができる{{Sfn|マレー|2014|p=67}}。
:<math>\frac{dx}{dt}=ax</math>
:<math>\frac{dy}{dt}=-dy</math>
これを[[行列]]表記すると、
:<math>
\begin{pmatrix}
\frac{dx}{dt} \\
\frac{dy}{dt} 
\end{pmatrix}
=
\begin{pmatrix}
a & 0 \\
0 & -d 
\end{pmatrix}
\begin{pmatrix}
x \\
y
\end{pmatrix}
</math>
となる。
:<math>A=
\begin{pmatrix}
a & 0 \\
0 & -d 
\end{pmatrix}
</math>
と置いたとき、''A'' の[[固有値]]は ''a'' と &minus;''d'' となり、正と負の固有値をもつので、平衡点 (0, 0) は[[鞍点]]となっている{{Sfn|Hirsch et al.|2007|p=246}}。また、少なくとも1つの固有値は正なので、指数関数的にズレが増加する不安定な平衡点である{{Sfn|マレー|2014|p=67}}。

平衡点 (''d''/''c'', ''a''/''b'') についても同様に、 平衡点近傍で系を次のように表すことができる{{Sfn|ハーバーマン|1992|p=116}}。
:<math>
\begin{pmatrix}
\frac{dx}{dt} \\
\frac{dy}{dt} 
\end{pmatrix}
=
\begin{pmatrix}
0 & -\frac{bd}{c} \\
\frac{ac}{b} & 0 
\end{pmatrix}
\begin{pmatrix}
x \\
y
\end{pmatrix}
</math>
固有値は <math>\pm i \sqrt{ad}</math> となる{{Sfn|Hirsch et al.|2007|p=247}}（<math>i</math> は[[虚数単位]]）。固有値は[[複素共役]]の純虚数となっており、平衡点 (''d''/''c'', ''a''/''b'') は渦心点となっている{{Sfn|Hirsch et al.|2007|p=60}}。したがって、平衡点近傍においては、平衡点周りで状態点が近づきも離れもしない、中立安定な平衡点となる{{Sfn|ハーバーマン|1992|pp=116&ndash;117}}。

===アイソクライン法による概略===
''x'' と ''y'' を変数とする平面（相平面）上で、''dx''/''dt'' = 0 または ''dy''/''dt'' = 0 を満たす直線に注目することで、個体数がどのような振る舞いを起こしているかの概略を知ることができる。このような手法を'''アイソクライン法'''や'''等傾斜線法'''と呼ぶ{{Sfn|寺本|1997|p=21}}{{Sfn|ハーバーマン|1992|pp=71&ndash;73}}。

相平面で横軸を ''x''、縦軸を ''y'' とする。現実の生物では個体数は正の値であるので、''x'' と ''y'' の値が正である相平面の[[直交座標系|第一象限]]が興味の対象となる{{Sfn|ハーバーマン|1992|p=111}}。相平面上では、''dx''/''dt'' = 0 を満たす直線とは ''y'' = ''a''/''b'' と ''x'' = 0 の直線であり、''dy''/''dt'' = 0 を満たす直線とは ''x'' = ''d''/''c'' と ''y'' = 0 の直線である{{Sfn|日本生態学会（編）|2015|p=43}}。このような ''dx''/''dt'' = 0 または ''dy''/''dt'' = 0 を満たす直線を '''アイソクライン'''や'''等傾斜線'''と呼ぶ<ref>{{cite book ja-jp |author=重定南奈子 |chapter=第1章 数理生態学 |pages=8|editor=日本数理生物学会 |title=生命・生物科学の数理 |publisher=岩波書店 |series=岩波講座 応用数学 4 [対象 8] |year=1993 |isbn=4-00-010514-0 }}</ref>{{Sfn|ハーバーマン|1992|p=112}}。前者の直線上では ''dx''/''dt'' = 0 であるから、解曲線がこの直線を通るとき、''x'' の値は変化せず、''y'' の値のみが変化する。よって、解曲線は直線を上下方向（''y''軸方向）にだけ通過する。そのため、この直線を'''傾き無限大のアイソクライン'''と呼ぶ{{Sfn|日本生態学会（編）|2005|p=33}}。一方、後者の直線上では ''dy''/''dt'' = 0 であるから、同じ理屈から解曲線はこの直線を左右方向（''x''軸方向）にだけ通過する。そのため、この直線を'''傾きゼロのアイソクライン'''と呼ぶ{{Sfn|日本生態学会（編）|2005|p=33}}。

相平面に ''y'' = ''a''/''b'' の水平線と ''x'' = ''d''/''c'' の鉛直線を描くと、平衡点 (''d''/''c'', ''a''/''b'') で2つの直線は交わり、相平面は4つの領域に分類される。''y'' = ''a''/''b'' の直線より上側の領域では、''dx''/''dt'' の値は常に負となっている。一方、下側の領域は ''dx''/''dt'' の値は常に正となる{{Sfn|ハーバーマン|1992|p=114}}。ここで、''dx''/''dt'' の値が正ということは ''x'' の値が増加している状態であり、負ということは ''x'' の値が減少している状態である{{Sfn|大串|1994|pp=71&ndash;72}}。よって、方程式の解の曲線は、''y'' = ''a''/''b'' の直線より上側の領域では左向きに進み、下側の領域では右向きに進むことが予測できる{{Sfn|日本生態学会（編）|2004|p=144}}。

また同様に、''x'' = ''d''/''c'' の直線より左側の領域では ''dy''/''dt'' の値は常に負で、右側の領域は ''dy''/''dt'' の値は常に正となる{{Sfn|ハーバーマン|1992|p=114}}。これによって上記と同じように、方程式の解の曲線は、''x'' = ''d''/''c'' の直線より左側の領域では下向きに進み、右側の領域では上向きに進むことが予測できる{{Sfn|日本生態学会（編）|2004|p=144}}。これらを組み合わせると、解の曲線は、平衡点 (''d''/''c'', ''a''/''b'') を中心にして反時計回りに回転する[[軌道 (力学系)|軌道]]となっていることが明らかになる{{Sfn|ハーバーマン|1992|p=114}}。
{{multiple image
 | align     = center
 | direction = horizontal
 | width     = 200
 | image1    = Lotka-Volterra equations isocline method 2.png
 | caption1  = 1) 解曲線は、''y'' = ''a''/''b'' の直線より上側領域では左向きに進み、下側領域では右向きに進む
 | image2    = Lotka-Volterra equations isocline method 1.png
 | caption2  = 2) 解曲線は、''x'' = ''d''/''c'' の直線より左側領域では下向きに進み、右側領域では上向きに進む
 | image3    = Lotka-Volterra equations isocline method 3.png
 | caption3  = 3) 解曲線は、平衡点 (''d''/''c'', ''a''/''b'') を中心にして反時計回りに回転する軌道となる
}}

===保存量===
ロトカ・ヴォルテラの方程式は[[力学系]]における[[保存系]]に該当し、保存量と呼ばれる量を持つ{{Sfn|マレー|2014|p=67}}。式から微分 ''dx''/''dy'' を求めると、
:<math>\frac{dx}{dy}=\frac{{dx}/{dt}}{{dy}/{dt}}=\frac{a x - b xy}{c xy - d y}</math>
となる。この変数分離形は
:<math>\frac{cx-d}{x}dx=\frac{a-by}{y}dy</math>
となり、両辺を積分して
:<math>H=cx + by - d \log x - a \log y</math>
が得られる{{Sfn|ハーバーマン|1992|p=119}}。ここで、log は[[自然対数]]である。右辺の ''H'' は一定の値を取る[[定数]]である。この式の意味は、時間経過に従って ''x'' と ''y'' が色々な値に変化しても、上式で与えられる ''H'' の値は常に同じに保たれるということである{{Sfn|寺本|1997|p=99}}。このような量は保存量や積分不変量と呼ばれ、保存量を持つ系は保存系と呼ばれる{{Sfn|寺本|1997|p=99}}。実際に ''H'' を ''t'' で微分すると、''dH''/''dt'' = 0 となり、''H'' が定数であることが確認できる{{Sfn|巌佐|1990|p=36}}{{efn|<math>\begin{align}
\frac{dH}{dt} &=\frac{\partial H}{\partial x}\frac{dx}{dt}+\frac{\partial H}{\partial y}\frac{dy}{dt} \\
&=\left( c-\frac{d}{x} \right) (ax-bxy) + \left( b-\frac{a}{y} \right) (cxy-dy) \\
&=acx - bcxy -ad +bdy +bcxy -bdy -acx +ad \\
&=0 
\end{align}
</math>
}}。平衡点 (''d''/''c'', ''a''/''b'') で ''H'' は最小値を取り、その値は
:<math>H_{min} = a+d-a \log \left( \frac{a}{b} \right)-d \log \left( \frac{d}{c} \right)</math>
となる<ref>{{cite journal |author=Shagi-Di Shih |year=1997 |month=December |title=THE PERIOD OF A LOTKA-VOLTERRA SYSTEM |url=http://journal.tms.org.tw/index.php/TJM/article/view/1423 |journal=Taiwanese Journal of Mathematics |publisher=The Mathematical Society of the Republic of China |volume=1 |issue=4 |page=453 |issn=2224-6851 |accessdate=2016-03-02 }}</ref>。''H'' &minus; ''H<sub>min</sub>'' はこの系における[[リアプノフ関数]]でもある{{Sfn|Hirsch et al.|2007|p=248}}。

===解曲線と個体数振動===
[[File:VolterraLotka.PNG|thumb|230px|解曲線は平衡点 (''d''/''c'', ''a''/''b'') を周回する[[閉曲線]]となっており、1つの閉曲線が一意の保存量を持つ。初期値によってどの閉曲線となるかが決定される]]
[[File:Lotka–Volterra equations, conserved quantity 3D plot.svg|thumb|230px|''x''-''y''相平面に高さ軸 ''H'' を加え、保存量 ''H'' と各閉曲線の関係を3次元的に示した図]]
上記のアイソクライン法による解析だけでは、解曲線の形状は確定しない。解曲線は、平衡点 (''d''/''c'', ''a''/''b'') を中心に反時計回りに回転していることは分かったが、平衡点を中心としてそこから離れていく[[渦巻]]形状なのか、逆に平衡点へ近づいていく渦巻形状なのか、あるいは円や楕円のように一周して元の点に戻る[[閉曲線]]なのか、などの可能性がある{{Sfn|日本生態学会（編）|2015|p=45}}。ロトカ・ヴォルテラの方程式の解は、これらの中の閉曲線に該当し、相平面の第一象限上で解曲線は平衡点 (''d''/''c'', ''a''/''b'') を中心にして一周する閉じた軌道を描く。これは、前述の保存量 ''H'' の存在などから証明される{{Sfn|Hirsch et al.|2007|p=248}}。

解曲線の形状は、純粋な円や楕円というよりは卵のような形となっている{{Sfn|寺本|1997|p=100}}。どの大きさの軌道を取るかは、被食者 ''x'' と捕食者 ''y'' の初期値 ''x''<sub>0</sub>, ''y''<sub>0</sub> によって決まる{{Sfn|日本生態学会（編）|2015|p=45}}。保存量 ''H'' の値は初期値 ''x''<sub>0</sub>, ''y''<sub>0</sub> によって決まり、''H'' の各値に1つの閉曲線が対応する{{Sfn|ハーバーマン|1992|p=119}}。さらに、''x'' と ''y'' の１周期中の平均量を計算すると、それらの値は、それぞれの平衡点 ''d''/''c'' と ''a''/''b'' に一致する{{Sfn|ハーバーマン|1992|pp=125&ndash;126}}。

[[File:Lotka Volterra equation Maple plot.png|right|330px|thumb|縦軸は個体数、横軸は時間で、捕食者（青）と被食者（赤）の個体数変動の時間変化を示している]]
解曲線が閉じた曲線であることは、被食者と捕食者の個体数は一定周期で[[振動]]していることも意味する{{Sfn|マレー|2014|p=66}}。個体数の時間発展波形は複雑な形状となる{{Sfn|ハーバーマン|1992|pp=122}}。捕食者と被食者の個体数変動の[[位相]]は1/4周期ほどずれており、
#被食者増加後に、捕食者増加
#捕食者増加後に、被食者減少
#被食者減少後に、捕食者減少
#捕食者減少後に、被食者増加
という変動の繰り返しを示す{{Sfn|日本生態学会（編）|2004|p=144}}。

個体数の範囲を平衡点近傍に限り、線形安定解析によって近似的な解析を行えば、それぞれの個体数変動の振動数を得ることもできる{{Sfn|ハーバーマン|1992|pp=115&ndash;117}}。このときの ''x'' と ''y'' は、上記の保存量 ''H'' と同じように、次のような関係で表される{{Sfn|ハーバーマン|1992|p=117}}。
:<math>C=\frac{a^2c^2}{b^2}x^2  +ad y^2.</math>
ここで、''C'' は一定値である。また、それぞれの個体数変動の[[振動数]] ''&omega;'' あるいは[[周期]] ''T'' は
:<math>\omega=\sqrt{ad},\quad T=\frac{2\pi}{\sqrt{ad}}</math>
で与えられる{{Sfn|ハーバーマン|1992|p=116}}{{Sfn|大串|1994|p=73}}。

===安定性===
前述のとおり、点 (''d''/''c'', ''a''/''b'') は中立安定な平衡点となっている。その周りに存在し得る軌道も初期値によって一つに決定され、一定の閉曲線を保ち続ける。すなわち、平衡点以外の軌道も、そこから離れも近づきもしない状態となっている。被食者も捕食者も絶滅することはなく、一方で、どちらの個体数も際限なく増え続けるということもない{{Sfn|Hirsch et al.|2007|p=249}}。

これは、系の外部から小さな乱れが加わった場合には、元の軌道から離れ、元に戻らないことも意味している。このような性質を「[[構造安定|構造的に不安定]]」などという{{Sfn|寺本|1997|p=100}}。現実にある多くの系を考えると、構造的に不安定であることは非現実的であることも多い{{Sfn|ハーバーマン|1992|pp=128&ndash;129}}{{Sfn|日本生態学会（編）|2015|pp=45&ndash;46}}。そのためより現実に合うようにモデルの改善が模索され、例えば、大域的に安定な[[リミットサイクル]]となるようにモデルの修正がされる{{Sfn|マレー|2014|pp=71&ndash;73}}{{Sfn|日本生態学会（編）|2015|pp=46&ndash;49}}。

==実際の生物における例==
===ダンコナとヴォルテラの研究===
[[File:Vito Volterra.jpg|thumb|150px|[[ヴィト・ヴォルテラ]] (Vito Volterra)]]
イタリアの生態学者ウンベルト・ダンコナ([[:en:Umberto D'Ancona|Umberto D'Ancona]]) は、漁業操業が低下した第一次世界大戦中に食用魚よりもサメなどの[[軟骨魚]]の年間漁獲率が増加したことに疑問を持った{{Sfn|ブラウン|2012|pp=224&ndash;225}}。これについて[[ヴィト・ヴォルテラ]]に相談を持ち掛け、ヴォルテラがこの現象を説明するためのモデル作成に取り組んだことが、ヴォルテラがロトカ・ヴォルテラの方程式を発案したきっかけである{{Sfn|Whittaker|1941|p=707}}{{Sfn|ブラウン|2012|pp=224&ndash;225}}。

ヴォルテラは、食用魚が被食者、軟骨魚が捕食者としてモデル（ロトカ・ヴォルテラの方程式）を作成した{{Sfn|ブラウン|2012|p=225}}。上記で説明したように、被食者の平均個体数は ''d''/''c'' で、捕食者の平均個体数は ''a''/''b'' である。漁業操業が行われており食用魚も軟骨魚も漁獲されているとすると、その効果は食用魚自然増加率の ''a'' を小さくして、軟骨魚自然減少率の ''d'' を大きくするように働くと考えることができる{{Sfn|Whittaker|1941|p=710}}。通常の操業量からある時期から操業量が低下したとする。これによって、通常の操業状態と相対的にみると、 ''a'' が大きくなり、 ''d'' が小さくなったということになる。したがって操業量低下により、被食者の平均個体数は減少し、捕食者の平均個体数が増加するということになる。これがダンコナの疑問に対するヴォルテラの説明である{{Sfn|ブラウン|2012|pp=229&ndash;230}}。

===周期的変動の例===
[[File:Milliers fourrures vendues en environ 90 ans odum 1953 en.jpg|thumb|250px|[[カンジキウサギ]]と[[カナダオオヤマネコ]]の捕獲頭数記録 (1845年-1935年)]]
[[File:Snowshoe Hare- Alert (6990916044).jpg|thumb|150px|[[カンジキウサギ]]]]
[[File:Canadian lynx by Keith Williams.jpg|thumb|150px|[[カナダオオヤマネコ]]]]
ロトカ・ヴォルテラの方程式で示された、被食者と捕食者の個体数が位相差を持ちながら一定振動を続ける振る舞いに近いといえる例は、実際の生物においていくつか確認されている。

野外環境における例としては、[[カナダ]]において、[[カンジキウサギ]]{{efn|[[カワリウサギ]]と記す文献もある{{Sfn|日本生態学会（編）|2015|p=40}}。}}とその捕食者である[[カナダオオヤマネコ]]の個体数が長期間にわたって振動していたデータがよく挙げられる{{Sfn|日本生態学会（編）|2004|pp=141&ndash;142}}{{Sfn|ハーバーマン|1992|pp=107&ndash;108}}。2つの個体数振動は、周期はほぼ同じで、位相は少しずれている{{Sfn|日本生態学会（編）|2015|p=40}}。ただし、このデータは個体数を直接観測したものではなく、毛皮取引を行っていた[[ハドソン湾会社]]による1845年から1935年までのカンジキウサギとカナダオオヤマネコの毛皮捕獲記録から、間接的に生息個体数を推定したものである{{Sfn|マレー|2014|pp=68&ndash;69}}。また、1973年のギルピン(M. E. Gilpin) による解析によれば、これらの個体数変動を相平面上にプロットすると軌道が時計回りとなっており、カンジキウサギがカナダオオヤマネコを捕食していると解釈できる奇妙な結果となっている{{Sfn|マレー|2014|pp=68&ndash;69}}。

環境を制御した飼育実験における例としては、ハフェイカー(C. B. Huffaker) による[[コウノシロハダニ]]とその捕食者である[[カブリダニ]]による飼育実験、[[内田俊郎]]による[[アズキゾウムシ]]とその寄生者である[[コマユバチ]]による飼育実験のデータが挙げられる{{Sfn|伊藤|1994|pp=80&ndash;81}}{{Sfn|日本生態学会（編）|2004|pp=141&ndash;142}}。ハフェイカーの実験では、単純な環境だと捕食が早すぎてどちらかの絶滅が起きてしまった。そのため、橋を設けたり扇風機を回したり環境を複雑にすることで、長期間にわたってそれぞれの個体数が振動しながら共存するデータを得ている{{Sfn|伊藤|1994|pp=80&ndash;81}}。

==モデルの改良==
現実にある多くの系を考えると、ロトカ・ヴォルテラの方程式
:<math>\begin{align}
\frac{dx}{dt} &= a x - b xy,\\
\frac{dy}{dt} &= c xy - d y\\
\end{align}</math>
は単純過ぎる部分がある。そのため、ロトカ・ヴォルテラの方程式を基礎としつつ、色々なモデルの研究がされてきた{{Sfn|マレー|2014|p=73}}。以下はその一例である。

問題点としてまず挙げられるのは、捕食者がいないときの被食者の増殖速度が ''ax'' となっており、青天井で増加し続ける点である。実際の系では、[[ロジスティック方程式]]のように、ある程度以上増加したら資源不足などが発生し、その増殖速度にブレーキがかかると考えるのが合理的である{{Sfn|マレー|2014|p=72}}。これを考慮に入れて、例えば、第1式の右辺第1項 ''ax'' をロジスティック型の ''ax''(1 &minus; ''x''/''K'') に置き換えたモデルが考えられる。ここで ''K'' は正の定数で、ロジスティックモデルにおける[[環境収容力]]である{{Sfn|Berryman|1992|p=1531}}。

また、被食者数に比例して無制限に捕食者増殖速度が増加する点も不自然である。これもある程度以上で飽和すると考えられる{{Sfn|マレー|2014|p=72}}。そのため、第1式の右辺第2項 &minus;''bxy'' を &minus;''bxy''/(1 + ''hx'') などと変形することが考えられる。ここで ''h'' は正の定数で、''x'' が増加してもこの項による捕食者1個体当たり増殖速度は ''b''/''h'' で飽和する{{Sfn|日本生態学会（編）|2015|p=46}}。

== ロトカ・ヴォルテラの競争モデル ==
{{main|ロトカ・ヴォルテラの競争方程式|競争 (生物)}}
類似の'''ロトカ・ヴォルテラの競争モデル'''

{{Indent|<math>
\begin{align}
\frac{dx}{dt} &= r_1x\,\frac{K_1-x-a_{21}y}{K_1},\\
\frac{dy}{dt} &= r_2y\,\frac{K_2-y-a_{12}x}{K_2}
\end{align}
</math>}}に関しては、[[ロトカ・ヴォルテラの競争方程式]]を参照。このモデルは、2種の[[個体群]]が捕食-被食関係というよりも[[競争 (生物)|競争関係]]にある場合を表している。このモデルも単に'''ロトカ‐ヴォルテラの式'''などと呼ばれることもある<ref>{{Cite web|和書|url=https://kotobank.jp/word/%E3%83%AD%E3%83%88%E3%82%AB%E2%80%90%E3%83%B4%E3%82%A9%E3%83%AB%E3%83%86%E3%83%A9%E3%81%AE%E5%BC%8F-791883 |title=法則の辞典の解説 ロトカ‐ヴォルテラの式【Lotka-Volterra equation】 |work=コトバンク |publisher=朝倉書店 |accessdate=2016-06-11}}</ref>。

ロトカ・ヴォルテラの競争モデルの解は捕食者-被食者モデルの場合と様相が異なり、それぞれの個体数 ''x'', ''y'' が周期変動しながら共存する解は存在しない。係数の値が ''K''<sub>1</sub> < ''K''<sub>2</sub>/''a''<sub>21</sub> かつ ''K''<sub>2</sub> < ''K''<sub>1</sub>/''a''<sub>12</sub> を満たすとき、''x'' と ''y'' は平衡点に収束し、それぞれの種が個体数一定で共存する。それ以外の場合にはどちらかの種が絶滅し、残った種の個体数は[[環境収容力]] ''K''<sub>1</sub> または ''K''<sub>2</sub> に落ち着く{{Sfn|巌佐|1990|p=15}}。

== 注釈 ==
{{notelist}}

==出典==
===脚注===
{{reflist|2}}

===文献リスト===
※文献内の複数個所に亘って参照したものを示す。
*{{cite book ja-jp
|author=R. ハーバーマン
|translator =稲垣宣生
|title=生態系の微分方程式
|year=1992
|edition=初版
|publisher=現代数学社
|isbn =4-7687-0307-0
|ref ={{Sfnref|ハーバーマン|1992}}
}}
*{{cite book ja-jp
|author=寺本英
|editor =川崎廣吉・重定南奈子・中島久男・東正彦・山村則男
|title=数理生態学
|year=1997
|edition=初版
|publisher=朝倉書店
|isbn =4-254-17100-5
|ref ={{Sfnref|寺本|1997}}
}}
*{{cite book ja-jp
|author=巌佐庸
|title=数理生物学入門―生物社会のダイナミックスを探る
|year=1990
|edition=初版
|publisher=HBJ出版局
|isbn =4-8337-6011-8
|ref ={{Sfnref|巌佐|1990}}
}}
*{{cite book ja-jp
|author=伊藤嘉昭
|title=生態学と社会―経済・社会系学生のための生態学入門
|year=1994
|edition=初版
|publisher=東海大学出版会
|isbn =4-486-01272-0
|ref ={{Sfnref|伊藤|1994}}
}}
*{{cite book ja-jp
|author=大串隆之
|chapter  = 3章 昆虫の個体群と群集
|pages=49&ndash;98
|title=昆虫生態学
|year=2014
|edition=初版
|publisher=朝倉書店
|isbn =978-4-254-42039-5
|ref ={{Sfnref|大串|1994}}
}}
*{{cite book ja-jp
|editor= 日本生態学会
|others=巌佐庸・舘田英典（担当編集委員）
|title=集団生物学
|series=シリーズ 現代の生態学 1
|year=2015
|edition=初版
|publisher=共立出版
|isbn =978-4-320-05744-9
|ref ={{Sfnref|日本生態学会（編）|2015}}
}}
*{{cite book ja-jp 
|editor= 日本生態学会 
|title=生態学入門 
|year=2004 
|edition=初版 
|publisher=東京化学同人 
|isbn =4-8079-0598-8 
|ref ={{Sfnref|日本生態学会（編）|2004}}
}}
*{{cite book ja-jp
|others =三村昌泰（総監修）、瀬野裕美・河内一樹・中口悦史・三浦岳（監修）
|author=ジェームス・D・マレー
|translator =勝瀬一登・吉田雄紀・青木修一郎・宮嶋望・半田剛久・山下博司
|title=マレー数理生物学入門
|year=2014
|edition=初版
|publisher=丸善出版
|isbn =978-4-621-08674-2
|ref ={{Sfnref|マレー|2014}}
}}
*{{cite book ja-jp 
|author=Morris W. Hirsch; Stephen Smale; Robert L. Devaney
|translator=桐木紳・三波篤朗・谷川清隆・辻井正人
|title=力学系入門　原著第2版―微分方程式からカオスまで
|publisher=共立出版
|edition=初版
|year=2007
|isbn=978-4-320-01847-1
|ref={{Sfnref|Hirsch et al.|2007}}
}}
*{{cite book ja-jp 
|author=M. ブラウン
|title=微分方程式 下―その数学と応用
|translator = 一樂重雄・河原正治・河原雅子・一樂祥子
|editor=シュプリンガー・ジャパン
|publisher=丸善出版
|year=2012
|isbn=978-4-621-06234-0
|ref={{Sfnref|ブラウン|2012}}
}}
*{{cite journal 
|author=Alan A. Berryman
|year=1992 
|month =Oct. 
|title=The Orgins and Evolution of Predator-Prey Theory
|url=http://www.jstor.org/stable/1940005  
|journal=Ecology 
|publisher=Ecological Society of America 
|volume=73 
|issue=5 
|pages=1530&ndash;1535 
|doi=10.2307/1940005
|ref ={{Sfnref|Berryman|1992}}
}}
*{{cite journal
|author=E. T. Whittaker
|title=Vito Volterra. 1860-1940
|year=1941 
|month=December
|journal=Obituary Notices of Fellows of the Royal Society
|volume=3
|issue=10
|publisher=Royal Society
|pages=690&ndash;729 
|jstor =769174 
|ref ={{Sfnref|Whittaker|1941}}
}}

==外部リンク==
{{Commonscat}}
*{{MathWorld|title=Lotka-Volterra Equations|urlname=Lotka-VolterraEquations}}
*{{Wayback|url=http://www.scholarpedia.org/article/Predator-prey_model |title=Predator-prey model |date=20061229085458}} - [[スカラーペディア]]百科事典「捕食者-被食者モデル」の項目。Lotka-Volterra Modelについての説明も含む。
*{{Kotobank|ロトカ‐ヴォルテラの式|2=法則の辞典}}

{{DEFAULTSORT:ろとかうおるてらのほうていしき}}
[[Category:常微分方程式]]
[[Category:力学系]]
[[Category:群集生態学]]
[[Category:個体群生態学]]
[[Category:数理生物学]]
[[Category:数学に関する記事]]