Модель Лотки — Вольтерры

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

История

Первые работы по математическому моделированию биологических систем относятся к началу XX века. Альфред Лотка, занимавшийся проблемами демографии и популяционной динамики, в 1925 году опубликовал книгу «Элементы физической биологии», в которой применил дифференциальные уравнения для описания взаимодействия видов. В том же году Вито Вольтерра, вдохновлённый наблюдениями своего зятя, биолога Умберто Д’Анконы, заинтересовался колебаниями численности рыб в Адриатическом море. Д’Анкона обратил внимание на то, что после Первой мировой войны, когда рыболовство сократилось, доля хищных рыб (акул и скатов) в уловах возросла. Вольтерра построил математическую модель, объясняющую этот феномен, и в 1926 году опубликовал статью «Variazioni e fluttuazioni del numero d’individui in specie animali conviventi» («Вариации и флуктуации численности особей в сосуществующих видах животных»). Впоследствии модель получила название по именам обоих учёных.

Математическая формулировка

Модель Лотки — Вольтерры описывается системой двух дифференциальных уравнений:

\[ \begin{cases} \frac{dx}{dt} = \alpha x - \beta x y \\ \frac{dy}{dt} = \delta x y - \gamma y \end{cases} \]

где:

• \(x(t)\) — численность популяции жертв (например, зайцев) в момент времени \(t\);
• \(y(t)\) — численность популяции хищников (например, волков) в момент времени \(t\);
• \(\alpha\) — коэффициент рождаемости жертв (скорость прироста популяции жертв в отсутствие хищников);
• \(\beta\) — коэффициент смертности жертв от хищников (скорость, с которой хищники потребляют жертв);
• \(\delta\) — коэффициент прироста хищников (эффективность потребления жертв, преобразование биомассы жертв в потомство хищников);
• \(\gamma\) — коэффициент смертности хищников (скорость естественной убыли популяции хищников в отсутствие жертв).

Все параметры \(\alpha, \beta, \gamma, \delta\) положительны.

Анализ поведения решений

Система имеет две особые точки (точки равновесия):

1. Тривиальная точка: \((x=0, y=0)\) — соответствует вымиранию обеих популяций. Эта точка является седловой, то есть неустойчивой.
2. Нетривиальная точка: \(\left(x = \frac{\gamma}{\delta}, y = \frac{\alpha}{\beta}\right)\) — соответствует сосуществованию видов. В этой точке численности обеих популяций остаются постоянными, если система находится в точном равновесии.

В окрестности нетривиальной точки решения системы имеют циклический характер: численности жертв и хищников совершают периодические колебания. Фазовые траектории представляют собой замкнутые кривые, окружающие точку равновесия. Период колебаний зависит от начальных условий и параметров системы. Амплитуда колебаний определяется удалённостью начальной точки от равновесия.

Ключевая особенность модели — закон сохранения (интеграл движения), который был найден Вольтеррой. Вдоль любой траектории выполняется соотношение:

\[ V(x, y) = \delta x + \beta y - \gamma \ln x - \alpha \ln y = \text{const} \]

Это означает, что система является консервативной (не диссипативной) — в ней нет затухания колебаний, и популяции могут бесконечно долго циклически изменять свою численность.

Интерпретация и биологический смысл

Модель Лотки — Вольтерры демонстрирует характерные циклы «хищник — жертва»:

• Когда численность жертв высока, хищники имеют abundant food, их популяция растёт.
• Рост числа хищников приводит к усилению давления на жертв, их численность начинает снижаться.
• Сокращение популяции жертв вызывает нехватку пищи для хищников, их численность падает.
• Уменьшение числа хищников позволяет популяции жертв восстановиться, и цикл повторяется.

Эти циклы не затухают, что является следствием математической структуры модели. В реальных экосистемах колебания обычно затухают из-за дополнительных факторов (ограниченность ресурсов, внутривидовая конкуренция, случайные воздействия), которые не учитываются в исходной модели.

Классические примеры

Наиболее известный пример, иллюстрирующий циклы Лотки — Вольтерры, — данные о динамике численности канадских рысей и зайцев-беляков, собранные компанией Hudson’s Bay Company на основе закупок пушнины в XIX — начале XX века. Наблюдаемые колебания с периодом около 10 лет хорошо согласуются с предсказаниями модели, хотя современные исследования показывают, что в реальности эти циклы обусловлены также влиянием климатических факторов и доступности растительной пищи.

Другой пример — колебания численности популяций в лабораторных экспериментах с простейшими (например, инфузориями и хищными коловратками), проведённых Г. Ф. Гаузе в 1930-х годах. В этих экспериментах наблюдались затухающие колебания, что указывало на ограниченность применимости простой модели.

Модификации и обобщения

Исходная модель Лотки — Вольтерры имеет ряд ограничений, которые привели к созданию её модификаций:

Модель с внутривидовой конкуренцией

Добавление логистического роста для жертв (ограничение по ёмкости среды) превращает систему в более реалистичную:

\[ \frac{dx}{dt} = \alpha x \left(1 - \frac{x}{K}\right) - \beta x y \]

где \(K\) — ёмкость среды для жертв. В этом случае колебания становятся затухающими, и система стремится к устойчивому равновесию.

Функциональный отклик

Вместо линейного члена \(\beta x y\) используются более сложные зависимости, например, функция Холлинга (тип II или III), учитывающая насыщение хищника:

\[ \frac{dx}{dt} = \alpha x - \frac{\beta x y}{1 + h x} \]

где \(h\) — время обработки жертвы. Это позволяет моделировать ситуации, когда хищник не может бесконечно увеличивать потребление жертв.

Модели с несколькими видами

Обобщения на случай нескольких видов жертв и хищников (системы Лотки — Вольтерры для \(n\) видов) используются в теоретической экологии для изучения устойчивости сообществ.

Применение

Модель Лотки — Вольтерры и её модификации находят применение в различных областях:

• Экология: анализ динамики популяций, управление промыслом (например, рыболовством), прогнозирование последствий интродукции видов.
• Эпидемиология: моделирование распространения инфекций (аналогия «жертва — здоровые, хищник — возбудитель»).
• Экономика: описание конкуренции фирм или взаимодействия отраслей (например, модели «инноватор — консерватор»).
• Химия: моделирование автоколебательных реакций (например, реакция Белоусова — Жаботинского).
• Социология: анализ динамики социальных групп, распространения мнений.

Критика и ограничения

Основные недостатки классической модели Лотки — Вольтерры:

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

Несмотря на эти ограничения, модель Лотки — Вольтерры остаётся важным дидактическим инструментом и отправной точкой для построения более сложных моделей в математической экологии.

Интересные факты

• Вито Вольтерра разработал модель, чтобы объяснить наблюдения своего зятя, но сам не был биологом — он был математиком, физиком и одним из основателей теории функционального анализа.
• Система Лотки — Вольтерры является одной из немногих нелинейных систем, для которой существует точный интеграл движения (закон сохранения).
• В 1930-х годах советский биолог Георгий Францевич Гаузе экспериментально проверил модель на простейших и показал, что в реальности колебания затухают, что привело к созданию «принципа Гаузе» (конкурентного исключения).
• Модель используется не только в естественных науках: в экономике она применялась для описания циклов «занятость — безработица» (модель Гудвина, 1967 год).

Источники

• Лотка А. Дж. «Элементы физической биологии» (1925)
• Вольтерра В. «Вариации и флуктуации численности особей в сосуществующих видах животных» (1926)
• Базыкин А. Д. «Математическая биофизика взаимодействующих популяций» (1985)
• Мюррей Дж. Д. «Математическая биология» (2002)
• Гаузе Г. Ф. «Борьба за существование» (1934)