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

Система уравнений Лотки-Вольтерры - одна из самых известных математических моделей в экологии. Она описывает взаимодействие двух популяций: жертвы (например, зайцев) и хищника (например, рысей). Модель предсказывает устойчивые периодические колебания численности обоих видов - феномен, наблюдаемый в природе повсеместно. Разберём структуру уравнений, их биологический смысл и ограничения модели.

Система уравнений и её параметры

Модель Лотки-Вольтерры записывается как система двух обыкновенных дифференциальных уравнений. Пусть $x(t)$ - численность жертв, $y(t)$ - численность хищников:

$\frac{dx}{dt} = \alpha x - \beta x y$

$\frac{dy}{dt} = \delta x y - \gamma y$

Здесь параметры имеют чёткий биологический смысл:

• $\alpha$ - скорость естественного прироста жертвы (рождаемость минус смертность без хищника);
• $\beta$ - коэффициент хищничества (частота встреч и успешных нападений);
• $\delta$ - эффективность конверсии жертвы в прирост хищника;
• $\gamma$ - скорость естественной убыли хищника (смертность без пищи).

Первое уравнение читается так: популяция жертв растёт экспоненциально с темпом $\alpha$, но каждая встреча с хищником убирает $\beta y$ единиц из прироста. Второе уравнение: хищники прибавляются пропорционально числу съеденных жертв $\delta xy$, но гибнут с темпом $\gamma y$.

Рассчитать для своей задачи

Равновесные состояния и траектории фазового портрета

Система имеет два равновесия. Тривиальное - исчезновение обоих видов $(0, 0)$ - неустойчиво. Нетривиальное (сосуществование) находится из условий $dx/dt = 0$ и $dy/dt = 0$:

$x^* = \frac{\gamma}{\delta}, \quad y^* = \frac{\alpha}{\beta}$

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

На фазовом портрете $(x, y)$ траектории представляют собой замкнутые кривые вокруг точки равновесия. Это означает, что орбита зависит от начальных условий: при малом начальном числе хищников амплитуда колебаний меньше, при большом - больше. Система является нейтрально устойчивой - возмущение не затухает и не нарастает, а переводит систему на другую замкнутую орбиту.

Биологические допущения и их ограничения

Классическая модель строится на нескольких упрощениях, важных для понимания её применимости:

1. Неограниченный рост жертвы в отсутствие хищника. Реальная популяция ограничена ресурсами среды (ёмкостью биотопа). Модель Лотки-Вольтерры игнорирует это, что ведёт к нереалистичным взрывам численности.

2. Линейная функциональная реакция. Модель предполагает, что скорость потребления жертв пропорциональна произведению $xy$. В действительности хищник насыщается - при большом числе жертв прирост потребления замедляется (кривая Холлинга типа II и III).

3. Однородная среда: все хищники и все жертвы встречаются с равной вероятностью, пространственная структура отсутствует.

4. Два вида: реальные пищевые сети включают десятки уровней и видов, конкурентов, паразитов.

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

Исторический контекст: данные по рыси и зайцу

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

Наиболее известные полевые данные - архивы Компании Гудзонова залива о заготовке шкур рыси и зайца-беляка за 1845-1935 годы. Численность обоих видов колебалась с периодом около 10 лет, причём пики популяции рыси следовали за пиками зайца с запаздыванием 1-2 года. Это в точности воспроизводит предсказание модели.

Рассчитать для своей задачи

Модификации: учёт ёмкости среды

Наиболее распространённое усовершенствование - замена экспоненциального роста жертвы логистическим (уравнение Розенцвейга-Макартура):

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

Здесь $K$ - ёмкость среды (максимальная численность жертвы без хищников). При добавлении этого члена поведение системы меняется: орбиты больше не замкнуты, а спиралью стягиваются к устойчивому равновесию. Для малых $K$ и достаточно «жадного» хищника возможен парадокс обогащения (Rosenzweig, 1971): улучшение ресурсной базы жертвы дестабилизирует систему вплоть до вымирания хищника.

Связанные концепции устойчивости популяций - в статье о генетическом грузе популяции.

Анализ устойчивости линеаризацией

Для определения характера точки равновесия $(x^*, y^*)$ применяют матрицу Якоби:

$J = \begin{pmatrix} 0 & -\beta x^* \\ \delta y^* & 0 \end{pmatrix} = \begin{pmatrix} 0 & -\gamma/\alpha \cdot \alpha \\ \delta \cdot \alpha/\beta & 0 \end{pmatrix}$

Собственные значения матрицы чисто мнимые: $\lambda_{1,2} = \pm i\sqrt{\alpha\gamma}$. Это и означает нейтральную устойчивость: малые возмущения не затухают и не нарастают, а приводят к осцилляциям с частотой $\omega = \sqrt{\alpha\gamma}$.

Реальный период колебаний: $T = 2\pi/\sqrt{\alpha\gamma}$. Для рыси и зайца: $\alpha \approx 1$ год$^{-1}$, $\gamma \approx 0{,}4$ год$^{-1}$, что даёт $T \approx 10$ лет - близко к наблюдаемым 9-11 годам.

Применение в экологических расчётах

Уравнения Лотки-Вольтерры применяются не только для описания «хищник-жертва», но и в более широком контексте:

• Конкуренция между видами - при замене знаков в уравнениях оба вида подавляют друг друга;
• Мутуализм - при смене знаков оба стимулируют рост партнёра;
• Эпидемиология - модели SIR описывают динамику болезни аналогичными уравнениями;
• Экономика - циклы «бум-спад» в рыночных моделях формально аналогичны уравнениям Вольтерры.

Математически близкую структуру имеет экологическая ниша Хатчинсона: оба подхода описывают пространство допустимых состояний популяций.

Рассчитать для своей задачи

Частые ошибки

• Путать $\delta$ и $\beta$. $\beta$ - коэффициент убывания жертвы при встрече с хищником; $\delta$ - эффективность конверсии этой встречи в прирост хищника. Обычно $\delta < \beta$ (не каждая встреча продуктивна).
• Считать равновесие устойчивым. Точка $(x^*, y^*)$ нейтрально устойчива: при классических уравнениях Лотки-Вольтерры траектории вечно ходят по кругу, а не сходятся к точке. Устойчивость появляется при добавлении логистического члена.
• Игнорировать начальные условия. Орбита (амплитуда колебаний) полностью определяется начальными значениями $x_0$ и $y_0$. Разные начальные условия - разные орбиты, хотя равновесие одно.
• Ожидать точного совпадения с природными данными. Модель качественная: она воспроизводит цикличность и запаздывание, но не точные числа. Реальные системы включают случайность, пространственную структуру, возрастной состав.
• Забывать о единицах $\beta$ и $\delta$. Если $x$ и $y$ - число особей, то $\beta$ имеет размерность [особь$^{-1}$·время$^{-1}$], а $\delta$ - безразмерный (или единицы/особь).

FAQ

Почему хищник достигает пика позже жертвы? Потому что рост хищника требует накопления жертв: сначала должна увеличиться кормовая база, затем с запаздыванием (вынашивание потомства, несколько поколений) вырастет популяция хищника. Это запаздывание - около четверти периода колебаний - прямо следует из фазового сдвига в уравнениях.

Возможно ли вымирание в модели Лотки-Вольтерры? В классической модели - нет. Орбиты замкнуты и никогда не достигают нулевых значений. Однако при добавлении стохастики или малых популяций случайные флуктуации могут привести к вымиранию. На практике популяция с численностью менее 10-20 особей уязвима к вымиранию по случайным причинам даже при «благоприятной» динамике.

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

Коротко

Уравнения Лотки-Вольтерры описывают взаимодействие хищника и жертвы через систему двух дифференциальных уравнений с четырьмя параметрами ($\alpha$, $\beta$, $\delta$, $\gamma$). Главный результат - устойчивые периодические колебания обеих популяций с фазовым сдвигом между ними. Равновесие $(x^* = \gamma/\delta,\ y^* = \alpha/\beta)$ нейтрально устойчиво: траектории замкнуты, а не сходятся к точке. Реальные экосистемы более сложны, но качественное предсказание цикличности подтверждается полевыми данными - прежде всего знаменитыми архивами рыси и зайца-беляка.