Дифференциальные уравнения · Уравнения Лотки — Вольтерры

Уравнения Лотки — Вольтерры

Уравнения Лотки — Вольтерры описывают популяции жертв и хищников: жертвы размножаются сами и служат пищей; хищники размножаются, питаясь жертвами, а без них вымирают. Их численности циклически растут и падают вокруг точки равновесия.

Уравнения Лотки — Вольтерры в этой визуализации

Здесь xx — число жертв, а yy — число хищников, оба в тысячах. Темпы изменения задаются на одну особь: жертвы растут с темпом α\alpha и гибнут от хищников с интенсивностью, пропорциональной числу хищников. Хищники умирают с интенсивностью γ\gamma и растут пропорционально числу жертв. Фазовая плоскость отслеживает точку (x,y)(x, y) от выбранного начала, а график времени снизу показывает обе численности год за годом.

Модель без вылова и предела
x′=x(α−βy),y′=y(δx−γ)x' = x(\alpha - \beta y), \qquad y' = y(\delta x - \gamma)

Темп изменения жертв α−βy\alpha - \beta y равен нулю на прямой y=α/βy = \alpha/\beta, а хищников δx−γ\delta x - \gamma — на прямой x=γ/δx = \gamma/\delta. Эти пунктирные линии нулевого роста пересекаются в точке равновесия (x∗,y∗)=(γ/δ,α/β)(x^*, y^*) = (\gamma/\delta, \alpha/\beta) и делят плоскость на четыре области. Ниже y∗y^* жертвы прибавляются, выше — убывают. Справа от x∗x^* хищники прибавляются, слева — убывают. Поэтому траектория идёт против часовой стрелки. Пик хищников наступает, когда они перестают расти: к этому моменту жертвы уже прошли свой пик и снова опустились до x∗x^*. Именно этот интервал времени показан как запаздывание.

Сам по себе обход не возвращает траекторию к началу: спираль тоже делает обороты. Замыкание обеспечивает величина, постоянная вдоль каждого решения: V=δx−γln⁡x+βy−αln⁡y,V′=(δx−γ)(α−βy)+(βy−α)(δx−γ)=0.\begin{aligned} V &= \delta x - \gamma \ln x + \beta y - \alpha \ln y, \\ V' &= (\delta x - \gamma)(\alpha - \beta y) + (\beta y - \alpha)(\delta x - \gamma) = 0. \end{aligned} Каждый уровень VV выше минимума образует замкнутую кривую вокруг точки равновесия. Траектория следует по ней к началу и повторяется. При пределе KK высокая плотность замедляет рост жертв, и VV уже не постоянно. В примере с пределом траектория закручивается к равновесию вместо замыкания.

Как бы велики ни были колебания, среднее за каждый полный цикл точно совпадает с точкой равновесия. Разделив уравнение хищников на yy, получим (ln⁡y)′=δx−γ(\ln y)' = \delta x - \gamma. За цикл yy возвращается к начальному значению, поэтому интеграл равен 0, а среднее xx равно γ/δ=x∗\gamma/\delta = x^*. Уравнение жертв таким же образом даёт среднее yy: α/β=y∗\alpha/\beta = y^*. Большие колебания меняют крайние значения и период, но не средние.

Непрерывный вылов обеих популяций с одинаковой интенсивностью на особь ε\varepsilon заменяет α\alpha на α−ε\alpha - \varepsilon, а γ\gamma — на γ+ε\gamma + \varepsilon. Без предела точка равновесия сдвигается к x∗=(γ+ε)/δx^* = (\gamma + \varepsilon)/\delta, y∗=(α−ε)/βy^* = (\alpha - \varepsilon)/\beta: в среднем больше жертв, меньше хищников. Средняя численность жертв — это число, при котором прирост и гибель хищников точно уравновешиваются. Поэтому оно растёт, когда вылов ускоряет гибель хищников.

На что обратить внимание

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

Управление

  • Начальная точка не может лежать на оси: популяция, начавшая с 0, останется нулевой.

Смотрите также