Numbers of snowshoe hare (yellow, background) and Canada lynx (black line, foreground) furs sold to the Hudson's Bay Company. Canada lynxes eat snowshoe hares.
← ホームへ戻る
ロトカ・ヴォルテラ方程式捕食者被食者個体群動態微分方程式

ロトカ・ヴォルテラ方程式捕食者と被食者のダイナミクスを解き明かす

🗓 2026年8月11日

自然界において、ある種の個体数が増えれば、それを餌とする別の種の個体数も増えるという相互作用が頻繁に見られます。このような捕食者(Predator)被食者(Prey)の複雑な関係を数学的に記述しようとしたのが「ロトカ・ヴォルテラ方程式」です。このモデルは、生物学のみならず経済学やマーケティングなど、競合する2つの要素が影響し合うシステムの解析に広く応用されています。

Key Facts

  • 相互依存のモデル: 被食者の増加が捕食者の増加を招き、その結果として被食者が減少するというサイクルを記述する。
  • 非線形微分方程式: 2つの1階非線形微分方程式で構成され、個体数の時間的変化を決定論的に示す。
  • 周期的な変動: 解は一般的に閉じた曲線となり、個体数が一定の周期で増減し続ける。
  • パラドックスの提示: 被食者の成長率を高めても、最終的に恩恵を受けるのは捕食者側であるという特性を持つ。

モデルの基礎概念と数学的構造

ロトカ・ヴォルテラモデルは、被食者の個体数密度を $x$、捕食者の個体数密度を $y$ とし、時間 $t$ による変化を以下の2つの式で表します。

  • 被食者の変化率: $\frac{dx}{dt} = \alpha x - \beta xy$
  • 捕食者の変化率: $\frac{dy}{dt} = -\gamma y + \delta xy$

ここで、$\alpha$ は被食者の最大成長率、$\beta$ は捕食者による被食者の死亡率への影響を示します。一方、$\gamma$ は捕食者の自然死亡率、$\delta$ は被食者の存在が捕食者の成長に与える影響を表しています。このモデルは、世代が絶えず重なり合っているという連続的な前提に基づいています。

この方程式系は、より広範な「コロモゴロフ個体群モデル」の一種であり、競争や共生、疾病などの多様な生態学的相互作用をモデル化する枠組みの一部となっています。

ภาพประกอบบทความ

生物学的な解釈と前提条件

このモデルを成立させるためには、いくつかの簡略化された仮定が置かれています。まず、被食者は常に十分な食物を得られること、そして捕食者の唯一の食物源は被食者であることです。現実の自然環境でこれらの条件が完全に満たされることは稀ですが、モデルが示す本質的な特性は多くの生物学的システムに共通しています。

特に注目すべきは、個体数の平衡状態における特性です。被食者の平衡密度は捕食者のパラメータに依存し、逆に捕食者の平衡密度は被食者のパラメータに依存します。例えば、被食者の成長率($\alpha$)を高めて環境を改善させても、結果として増加するのは捕食者の個体数であり、被食者の平衡密度は変わらないという現象が起こります。これは「豊饒化のパラドックス」とも関連しています。

歴史的な実例として、第一次世界大戦中のアドリア海での漁獲データが挙げられます。漁業活動の減少により被食魚の成長率が上がった結果、皮肉にも捕食魚の割合が増加したことが観察されており、ヴォルテラはこの現象を説明するためにモデルを構築しました。

Numbers of snowshoe hare (yellow, background) and Canada lynx (black line, foreground) furs sold to the Hudson's Bay Company. Canada lynxes eat snowshoe hares.

個体数変動のダイナミクスと解析

方程式を解析すると、捕食者の個体数変動は被食者の変動から90度遅れて追従する、単純調和振動に近い挙動を示すことがわかります。時間軸で個体数をプロットすると、両種が交互にピークを迎える周期的な波形が現れます。

Prey and predator populations over time

また、時間変数を除去して $x$ と $y$ の関係を直接的に表すと、解は相平面上で「閉じた曲線」を描きます。これは、ある一定の条件下で個体数が永遠にサイクルを繰り返すことを意味します。

Population dynamics for rabbit and fox problem mentioned aside.

しかし、数学的なモデルと現実には乖離もあります。例えば、計算上は個体数が極めて低い値まで減少しても回復しますが、現実には個体数が少なすぎると偶然の変動で絶滅に至ります。これは「アット・フォックス問題(atto-fox problem)」と呼ばれ、理論上の極小値が生物学的な絶滅閾値を下回るリスクを指摘しています。

Phase-space plot for the predator prey problem for various initial conditions of the predator population.

モデルの要約

ロトカ・ヴォルテラモデルの主要な構成要素と意味を以下の表にまとめます。

ロトカ・ヴォルテラモデルの構成要素まとめ
変数・パラメータ 意味 生物学的影響
$x$ 被食者の個体数密度 捕食者の餌となる資源量
$y$ 捕食者の個体数密度 被食者を抑制する圧力
$\alpha$ 被食者の成長率 餌が豊富にある時の増加速度
$\beta$ 捕食効率 捕食者が被食者を捕食する割合
$\gamma$ 捕食者の死亡率 餌がない時の減少速度
$\delta$ 増殖効率 被食者の摂取が捕食者の出生に寄与する度合い

歴史的背景と応用

このモデルは、1910年にアルフレート・ロトカが化学反応の理論として提案したことに始まります。その後、1920年代に生物学的システムへと拡張されました。ほぼ同時期に、イタリアの数学者ヴィト・ヴォルテラが独立して同様の方程式を導き出したため、現在は両者の名を冠して呼ばれています。

また、この数学的構造は生物学以外にも応用されています。1960年代にはリチャード・グッドウィンらによって経済理論に導入され、景気循環や市場のダイナミクスを分析するためのツールとして活用されました。

Frequently Asked Questions

ロトカ・ヴォルテラ方程式で予測される個体数の動きはどうなりますか?

被食者が増えると、それを餌とする捕食者も増えます。しかし、捕食者が増えすぎると被食者が激減し、それに伴い餌を失った捕食者も減少します。捕食者が減ると再び被食者が増加し、このサイクルが周期的に繰り返されます。

このモデルの最大の弱点は何ですか?

現実の生態系にある「環境収容力(餌や空間の限界)」が考慮されていない点です。また、個体数が極端に少なくなった際に、数学的には回復可能でも現実には絶滅してしまうという「離散的な個体数」の問題(アット・フォックス問題)があります。

「豊饒化のパラドックス」とは具体的にどういうことですか?

被食者の餌を増やして被食者が生きやすい環境を作っても、その結果として捕食者の数だけが増え、最終的な被食者の個体数は増えないという現象です。環境改善のメリットが、食物連鎖の上位者に吸収されてしまうことを意味します。

生物学以外にどのような分野で使われていますか?

主に経済学やマーケティングで利用されています。例えば、2つの競合製品の市場シェアの奪い合いや、資本と雇用のサイクルなど、相互に影響し合う2つの変数の変動を分析する際に用いられます。

References

  1. Cooke, D.; Hiorns, R. W.; et al. (1981). The Mathematical Theory of the Dynamics of Biological Populations. Vol. II. Academic Press.
  2. Freedman, H. I. (1980). Deterministic Mathematical Models in Population Ecology. .
  3. Brauer, F.; Castillo-Chavez, C. (2000). Mathematical Models in Population Biology and Epidemiology. .
  4. Hoppensteadt, F. (2006). "Predator-prey model". . 1 (10): 1563. :2006SchpJ...1.1563H. :10.4249/scholarpedia.1563.
  5. Beals, M.; Gross, L.; Townsend, C.R. (1999). "Predator-Prey Dynamics". www.tiem.utk.edu. Archived from the original on 2012-12-15. Retrieved 2018-01-09.
  6. Gilpin, M. E. (1973). "Do hares eat lynx?". American Naturalist. 107 (957): 727–730. :1973ANat..107..727G. :10.1086/282870.  84794121.
  7. Jost, C.; Devulder, G.; Vucetich, J.A.; Peterson, R.; Arditi, R. (2005). "The wolves of Isle Royale display scale-invariant satiation and density dependent predation on moose". J. Anim. Ecol. 74 (5): 809–816. :10.1111/j.1365-2656.2005.00977.x.
  8. Pan, A.; Pourziaei, B.; Huang, H. (2015-06-03). "Effect of Ocean Iron Fertilization on the Phytoplankton Biological Carbon Pump". Advances in Applied Mathematics and Mechanics. 3 (1): 52–64. :10.4208/aamm.10-m1023.  124606355.
  9. Prasolov, Alexander V. (2016). Some quantitative methods and models in economic theory. Economic issues, problems and perspectives. New York: Nova Publishers.  .
  10. Hung, Hui-Chih; Chiu, Yu-Chih; Wu, Muh-Cherng (2017). "A Modified Lotka–Volterra Model for Diffusion and Substitution of Multigeneration DRAM Processing Technologies". Mathematical Problems in Engineering. 2017: 1–12. :10.1155/2017/3038203.  1024-123X.

📸 フォトギャラリー

Numbers of snowshoe hare (yellow, background) and Canada lynx (black line, foreground) furs sold to the Hudson's Bay Company. Canada lynxes eat snowshoe hares.
Prey and predator populations over time
Population dynamics for rabbit and fox problem mentioned aside.
Phase-space plot for the predator prey problem for various initial conditions of the predator population.
ภาพประกอบบทความ