Хищник-жертва: популяционная динамика Лотки-Вольтерры

simulator intermediate ~10 min
Загрузка симуляции...
Период осцилляций ≈ 15,7 ед. времени, пик жертв ≈ 44, пик хищников ≈ 18

При параметрах по умолчанию (α=1,1, β=0,4, γ=0,4, δ=0,1) система Лотки-Вольтерры создаёт устойчивые осцилляции с периодом около 15,7 единиц времени. Популяции жертв и хищников циклируют в противофазе — пик хищников отстаёт от пика жертв на четверть периода, что является отличительной чертой связанной осциллирующей динамики.

Формула

dx/dt = αx − βxy
dy/dt = δxy − γy
Equilibrium: x* = γ/δ, y* = α/β
Conserved quantity: V = δx − γ·ln(x) + βy − α·ln(y)

Модель хищник-жертва Лотки-Вольтерры

В 1925 году Альфред Лотка предложил математическую модель для осциллирующих химических реакций, которую независимо переоткрыл Вито Вольтерра в 1926 году для объяснения колебаний уловов рыбы в Адриатике. Получившиеся уравнения Лотки-Вольтерры стали одной из влиятельнейших моделей математической экологии, обеспечив теоретический фундамент для понимания динамики хищник-жертва.

Уравнения

Модель состоит из двух связанных дифференциальных уравнений:

dx/dt = αx − βxy — жертвы растут экспоненциально со скоростью α при отсутствии хищников. Член βxy описывает хищничество: встречи хищников и жертв пропорциональны произведению их численностей (принцип массового действия).

dy/dt = δxy − γy — хищники вымирают экспоненциально со скоростью γ без пищи. Член δxy описывает размножение хищников за счёт успешного хищничества, где δ — эффективность конвертации съеденных жертв в новых хищников.

Бесконечные осцилляции

Самое поразительное свойство модели — предсказание бесконечных незатухающих осцилляций. Система обладает сохраняющейся величиной (аналогом энергии в физике): V = δx − γ·ln(x) + βy − α·ln(y), которая остаётся постоянной вдоль любой траектории. Этот закон сохранения гарантирует, что орбиты на фазовой плоскости — замкнутые кривые: популяции циклируют вечно без затухания и усиления.

Точка равновесия (x* = γ/δ, y* = α/β) находится в центре этих орбит. Траектории ближе к равновесию имеют меньшую амплитуду осцилляций; дальние — большую. Начальные условия определяют, по какой орбите движется система.

Экологические прозрения

Запаздывание пика хищников на четверть периода от пика жертв — ключевое предсказание, подтверждённое во многих реальных экосистемах. Классический пример — записи Компании Гудзонова залива, показывающие связанные осцилляции зайца-беляка и канадской рыси на протяжении более века. Хотя реальные экосистемы включают сложности, отсутствующие в базовой модели — ёмкость среды, функциональные ответы, пространственную структуру — рамка Лотки-Вольтерры остаётся отправной точкой для понимания всех взаимодействий хищник-жертва.

Исследуйте фазовый портрет в правой части визуализации: обратите внимание на замкнутую орбиту, прочерчивающую вечный цикл. Форма и размер орбиты меняются с параметрами — увеличение эффективности хищничества (δ) или скорости рождения жертв (α) кардинально меняет динамику.

Частые вопросы

Что такое модель Лотки-Вольтерры?

Уравнения Лотки-Вольтерры — пара связанных нелинейных ОДУ первого порядка, моделирующих динамику двух взаимодействующих видов — хищника и жертвы. Независимо предложены Альфредом Лоткой (1925) и Вито Вольтеррой (1926), модель предсказывает бесконечные осцилляции обеих популяций.

Каковы уравнения Лотки-Вольтерры?

dx/dt = αx − βxy (рост жертв минус хищничество) и dy/dt = δxy − γy (рост хищников от хищничества минус естественная гибель), где x — популяция жертв, y — хищников, α, β, γ, δ — положительные параметры.

Почему популяции хищника и жертвы осциллируют?

Осцилляции возникают из петли обратной связи: обилие жертв подпитывает рост хищников, но рост хищников сокращает жертв, что затем вызывает спад хищников от голода, позволяя жертвам восстановиться. Эта отрицательная обратная связь с задержкой создаёт устойчивые циклы.

Что такое фазовый портрет в экологии?

Фазовый портрет строит популяцию хищников (y) против популяции жертв (x), устраняя время как явную переменную. В модели Лотки-Вольтерры траектории образуют замкнутые орбиты вокруг точки равновесия (x* = γ/δ, y* = α/β), форма орбиты зависит от начальных условий.

Источники

Встроить

<iframe src="https://homo-deus.com/lab/evolution/predator-prey/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub