Я догоняю, ты убегаешь

Могут ли сложные математические инструменты применяться в биологии? Могут, если биологи изучают сложные динамические системы, например взаимодействие разных видов животных в естественной среде. Американец Альфред Лотка и итальянец Вито Вольтерра разработали модель, позволяющую описывать, как будет меняться поголовье хищников и их травоядных жертв в зависимости от множества привходящих условий. Это наш второй материал о самых интересных дифференциальных уравнениях (с первым можно ознакомиться здесь). Если вы читаете нас с телефона, переключайте страницу на десктопную версию, так вы сможете увидеть интерактивный график целиком.

Изначально Альфред Лотка вообще не планировал создавать никаких математических моделей. Он собирался разработать новую предметную область — «физическую биологию» — и поэтому начиная с 1902 года стал публиковать небольшие статьи, посвященные этой теме.

Параллельно с этим его все более интересовало применение математических методов в биологии. Идеи Лотки, однако, не получили широкого распространения — в то время американский ученый не имел широких связей в научной среде и работал в одиночестве.

Ситуация изменилась в 1920 году, когда статьи Лотки привлекли внимание биолога и статистика Раймонда Пирла, который нашел в них близкие для себя идеи: Пирл интересовался ростом популяции в пределах одного вида.

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

Позже Лотка расширил это наблюдение до общего случая взаимодействия типа «хищник-жертва».

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

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

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

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

Система Лотки-Вольтерры является первоначальной и простейшей системой (усложненные системы будут рассмотрены ниже) для описания модели «хищник-жертва», то есть популяции хищников и популяции жертв, взаимодействующих в какой-то среде: жертвы едят растительность, хищники — жертв:

где

• — численность жертв (травоядных);
• — численность хищников;
• — вероятность того, что травоядные размножатся;
• — вероятность того, что травоядное будет съедено хищником;
• — вероятность того, что хищник умрет от голода;
• — вероятность того, что хищнику хватит еды на дальнейшее размножение.

Из системы сразу следует, что если жертв нет ( = 0), то хищники будут вымирать экспоненциально с неким начальным коэффициентом ( согласно уравнению).

Схожую ситуацию получаем при полном отсутствии хищников ( = 0):

Рост жертв получается экспоненциальным с некой заранее заданной константой (). Стоит отметить, что в данной модели принимаются несколько допущений:

• Количество пищи для травоядных не ограничено;
• Ни жертвы, ни хищники не эмигрируют из среды;
• Никакие другие животные не мигрируют в среду;
• Данная модель не учитывает вымирание животных по причине старения и прочих внешних воздействий.

Ниже можно посмотреть, как будут меняться размеры популяции в зависимости от заданных начальных условий (если вы читаете нас с телефона, переключайте страницу на десктопную версию, так вы сможете увидеть интерактивный график целиком):

Особые точки

Найдем особые точки, которыми обладает система:

Понятно, что при (0) = 0, (0) = 0 особой точкой будет как раз (0, 0), но этот случай не интересен, так как в нулевой момент времени животные обоих видов отсутствуют и, что логично, дальше не появляются.

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

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

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

Миграция животных

Существует усложнение стандартной модели Лотки-Вольтерры, при котором учитывается миграция животных. В такой модели система принимает вид:

где C(), D() — миграция травоядных и хищников соответственно. Причем функции могут задаваться двумя разными способами. В первом случае:

То есть в каждый момент времени особи обеих популяций константно мигрируют.

Второй случай менее примитивен:

То есть функции показывают отношение мигрирующих животных к общей массе. Для обоих случаев верно, что при положительных константах , особи будут прибывать в среду, при отрицательных — покидать ее, а при нулевых миграции не будет.

При данном задании модели возможны различные интересные комбинации миграции двух видов животных. Рассмотрим ниже пару примеров, чтобы было понятно, как это происходит.

Миграция травоядных в среду

Рассмотрим случай, когда мигрируют только жертвы по второму способу задания функций, то есть:

[Найдем особые точки (сразу рассматриваем случай, когда размеры популяций ненулевые):

А теперь исследуем ситуацию на устойчивость: найдем якобиан и собственные значения:

Если  = 0, то получаем особую точку типа центр, иначе — фокус, причем, если  < 0, то точка будет устойчивой, если  > 0, то неустойчивой.

Миграция хищников в среду

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

Опять найдем особые точки (а точнее, одну особую, так как случай, когда жертв нет, не интересен — тогда хищники просто вымрут):

Теперь так же, как и в предыдущем случае, исследуем особую точку:

Приходим к такие же выводам: если  = 0, то получаем особую точку типа центр, иначе — фокус. Если  < 0, то точка будет устойчивой, если  > 0 — неустойчивой.

Ниже можно поиграть с начальными условиями и понаблюдать, как они сказываются на поведение системы:

Многомерный случай

Вито Вольтерра вывел уравнения для -мерного случая, которые записываются в виде:

Поделиться