Лотка–Вольтерра
Уравнения Лотки–Вольтерры, также известные как модель «хищник-жертва», описывают динамику двух взаимодействующих биологических популяций. Система управляется парой нелинейных дифференциальных уравнений первого порядка: dN/dt = αN - βNP для популяции жертв N(t) и dP/dt = δNP - γP для популяции хищников P(t). Здесь α — собственная скорость роста жертв в отсутствие хищников, β — скорость, с которой хищники ловят жертв, δ — эффективность превращения съеденных жертв в новых хищников, а γ — уровень смертности хищников. Произведения βNP и δNP моделируют ключевое взаимодействие: предполагается, что скорость встреч между видами пропорциональна произведению их численностей. Эта модель упрощает реальные экосистемы, игнорируя пространственные эффекты, возрастную структуру, ёмкость среды для жертв и другие потенциальные источники пищи для хищника. Несмотря на упрощения, модель воспроизводит фундаментальное колебательное поведение, наблюдаемое в некоторых природных системах, где численность хищников отстаёт от численности жертв. Работая с этим симулятором, студенты смогут исследовать фазовую плоскость — график зависимости численности хищников от численности жертв — и наблюдать, как траектории образуют замкнутые орбиты вокруг ненулевой точки равновесия (γ/δ, α/β). Они узнают, как начальные условия влияют на амплитуду и фазу циклов, увидят интегрирование дифференциальных уравнений методом Рунге–Кутты (RK4) и поймут концепцию нейтрально устойчивого центрального равновесия в консервативной системе.
Для кого: Студенты бакалавриата по математике, экологии или курсам динамических систем, изучающие связанные дифференциальные уравнения и биологическое моделирование.
Ключевые понятия
- Уравнения Лотки–Вольтерры
- Модель «хищник-жертва»
- Дифференциальные уравнения
- Фазовая плоскость
- Точка равновесия
- Метод Рунге–Кутты
- Динамика популяций
- Нейтральная устойчивость
Как это работает
Карикатура **хищник–жертва**: **dN/dt = αN − βNP**, **dP/dt = δNP − γP**. Внутреннее равновесие **(N*, P*) = (γ/δ, α/β)**. Траектории в плоскости **(N, P)** для этой модели — **замкнутые** кривые вокруг него (логистический предел для жертвы не вводится). Интегрирование **RK4**; **красная точка** отмечает **(N*, P*)**.
Основные формулы
Часто задаваемые вопросы
- Всегда ли популяции в модели Лотки–Вольтерры колеблются вечно?
- В представленной идеализированной модели — да, колебания являются периодическими и продолжаются бесконечно, образуя замкнутые петли на фазовой плоскости. Это происходит потому, что модель консервативна и не имеет затухания или ограничивающих факторов. В реальных экосистемах такие идеальные незатухающие колебания редки, поскольку модель не учитывает факторы вроде ёмкости среды, которые стабилизировали бы систему в предельный цикл или неподвижную точку.
- Что представляет собой точка равновесия?
- Ненулевая точка равновесия (γ/δ, α/β) представляет собой постоянные размеры популяций, при которых скорости изменения для обоих видов равны нулю. Если система начинает работу точно в этой точке, популяции остаются постоянными. Однако это равновесие является нейтрально устойчивым: любое малое возмущение приводит к устойчивым колебаниям вокруг него, а не к возврату к равновесию или неограниченному расхождению.
- Почему член взаимодействия пропорционален N*P?
- Произведение N*P, известное как член закона действующих масс, предполагает, что скорость встреч хищников и жертв пропорциональна вероятности случайной встречи особей из каждой группы. Это распространённое упрощение в химической кинетике и эпидемиологии (где оно моделирует передачу болезни), основанное на принципе, что большее количество особей обоих типов приводит к более частым взаимодействиям.
- Каковы основные ограничения этой базовой модели?
- Ключевые ограничения включают предположение о неограниченном росте жертв (отсутствие ёмкости среды), зависимость хищника только от одного вида жертв, отсутствие временных задержек (например, в размножении хищника) и игнорирование пространственного перемещения или убежищ для жертв. Более сложные модели, такие как модель Розенцвейга–МакАртура, добавляют ёмкость среды для жертв, чтобы создать более реалистичную динамику.
Ещё из «Визуализация математики»
Другие симуляторы в этой категории — или все 26.
Логистический рост
dN/dt = rN(1−N/K); точная S-образная кривая и ёмкость среды K.
Матрица 2×2 и собственные векторы
Деформация сетки под действием M; стрелки вдоль собственных направлений для вещественных λ.
Сглаживание Савицкого–Голэя
Зашумлённый косинус vs свёртка SG(7,2) — лучше сохраняет пики по сравнению с широким прямоугольным фильтром.
Цепь Маркова (Погода)
Двухсостоянияя цепь Солнце/Дождь: матрица P, стационарное распределение π, эмпирическое и теоретическое.
Градиентный спуск (2D)
Линии уровня функции f(x,y) и траектория (x,y) ← (x,y) − η∇f; чаша или эллиптическая впадина.
Диаграмма Минковского
Световой конус и повёрнутые оси в 1+1 измерениях; γ от v.