Колебания системы хищник-жертва (модель Лотки-Вольтера). Равновесие хищник—жертва Модель системы хищник жертва

позволяющее сформулировать условие сосуществования видов:

Произведение коэффициентов межпопуляционного взаимодействия меньше произведения коэффициентов внутри популяционного взаимодействия.

Действительно, пусть естественные скорости роста двух рассматриваемых видов a 1 , a 2 одинаковы. Тогда необходимым для устойчивости условием будет

c 2 > b 12 ,c 1 > b 21 .

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

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

x 1 (a 1 –c 1 x 1 – b 12 x 2) = 0 (dx 1 /dt = 0),

x 2 (a 2 –b 21 x 1 – c 2 x 2) = 0 (dx 2 /dt = 0),

При этом получим уравнения для главных изоклин системы

x 2 = – b 21 x 1 / c 2 +a 2 /c 2 , x 2 = 0

– уравнения изоклин вертикальных касательных.

x 2 = – c 1 x 1 / b 12 +a 1 /b 12 , x 1 = 0

– уравнения изоклин вертикальных касательных. Точки попарного пересечения изоклин вертикальных и горизонтальных касательных систем представляют собой стационарные решения системы уравнений (9.2.), а их координаты суть стационарные численности конкурирующих видов.

Возможное расположение главных изоклин в системе (9.2) изображено на рис.9.1. Рис. 9.1 а соответствует выживанию вида x 1 , рис. 9.1 б – выживанию вида x 2 , рис. 9.1 в – сосуществованию видов при выполнении условия (9.6). Рисунок 9.1 г демонстрирует триггерную систему. Здесь исход конкуренции зависит от начальных условий. Ненулевое для обоих видов стационарное состояние (9.5) – неустойчивое. Это – седло, через которое проходит сепаратриса, отделяющая области выживания каждого из видов.

Рис. 9.1. Расположение главных изоклин на фазовом портрете вольтерровской системы конкуренции двух видов (9.2) при разном соотношении параметров. Пояснения в тексте.

Для изучения конкуренции видов ставились эксперименты на самых различных организмах. Обычно выбирают два близкородственных вида и выращивают их вместе и по отдельности в строго контролируемых условиях. Через определенные промежутки времени проводят полный или выборочный учет численности популяции. Регистрируют данные по нескольким повторным экспериментам и анализируют. Исследования проводили на простейших (в частности, инфузориях), многих видах жуков рода Tribolium, дрозофиллах, пресноводных ракообразных (дафниях). Много экспериментов проводилось на микробных популяциях (см. лекция 11). В природе также проводили эксперименты, в том числена планариях (Рейнольдс), двух видах муравьев (Понтин) и др. На рис. 9.2. изображены кривые роста диатомовых водорослей, использующих один и тот же ресурс (занимающих одну и ту же экологическую нишу). При выращивании в монокультуре Asterionella Formosa выходит на постоянный уровень плотности и поддержвает концентрацию ресурса (силиката) на постоянно низком уровне. Б. При выращивании в монокультуре Synedrauina ведет себя сходным образом и поддерживает концентрацию силиката на еще более низком уровне. В. При совместном культивировании (в двух повторностях) Synedrauina вытесняет Asterionella Formosa . По-видимому, Synedra

Рис. 9.2. Конкуренция у диатомовых водорослей. а ‑ при выращивании в монокультуре Asterionella Formosa выходит на постоянный уровень плотности и поддерживает концентрацию ресурса (силиката) на постоянно низком уровне.б ‑ при выращивании в монокультуре Synedrauina ведет себя сходным образом и поддерживает концентрацию силиката на еще более низком уровне. в ‑ при совместном культивировании (в двух повторностях) Synedruina вытесняет Asterionella Formosa . По-видимому, Synedra выигрывает конкуренцию благодаря своей способности к более полному использованию субстрата (см. также Лекцию 11).

Широко известны эксперименты по изучению конкуренции Г. Гаузе, продемонстрировавшие выживание одного из конкурирующих видов и позволившие ему сформулировать «закон конкурентного исключения». Закон гласит, что в одной экологической нише может существовать только один вид. На рис. 9.3. приведены результаты экспериментов Гаузе для двух видов Parametium , занимающих одну экологическую нишу (рис. 9.3 а, б) и видами, занимающими разные экологические ниши (рис. 9.3. в).

Рис. 9.3. а - Кривые роста популяций двух видов Parametium в одновидовыхкультурах. Черные кружки – P Aurelia , белые кружки – P . Caudatum

б - Кривые роста P Aurelia и P . Caudatum в смешанной культуре.

По Gause , 1934

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

Система ХИЩНИК+ЖЕРТВА

(9.7)

Здесь, в отличие от (9.2) знаки b 12 и b 21 - разные. Как и в случае конкуренции, начало координат

(9.8)

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

,(9.9)

(9.10)

(9.11)

Таким образом, возможно выживание только жертвы (9.10), только хищника (9.9) (если у него имеются и другие источники питания) и сосуществование обоих видов (9.11). Последний вариант уже был рассмотрен нами в лекции 5. Возможные типы фазовых портретов для системы хищник-жертва представлены на рис. 9.4.

Изоклины горизонтальных касательных представляют собой прямые

x 2 = – b 21 х 1 /c 2 + a 1 /c 2 , х 2 = 0,

а изоклины вертикальных касательных – прямые

x 2 = ­– c 1 х 1 /b 12 + a 2 /b 12 , х 1 = 0.

Стационарные точки лежат на пересечении изоклин вертикальных и горизонтальных касательных.

Из рис. 9.4 видно следующее. Система хищник – жертва (9.7) может иметь устойчивое положение равновесия, в кот o ром популяция жертв полностью вымерла ( ) и остались только хищники (точка 2 на рис. 9.4 а). Очевидно, такая ситуация может реализоваться лишь в случае, если кроме рассматриваемого вида жертв х 1 хищник х 2 – имеет дополнительные источники питания. Этот факт в модели отражается положительным членом в правой части уравнения для х 2 . Особые точки (1) и (3) (рис. 9.4 а ) являются неустойчивыми. Вторая возможность – устойчивое стационарное состояние, в котором популяция хищников полностью вымерла и остались одни жертвы – устойчивая точка (3) (рис. 9.4 6 ). Здесь особая точка (1) – также неустойчивый узел.

Наконец, третья возможность – устойчивое сосуществование популяций хищника и жертвы (рис. 9.4 в ), стационарные численности которых выражаются формулами (9.11).

Как и в случае одной популяции (см. Лекция 3), для модели (9.7) можно разработать стохастическую модель, но для нее нельзя получить решение в явном виде. Поэтому мы ограничимся общими рассуждениями. Допустим, например, что точка равновесия находится на некотором расстоянии от каждой из осей. Тогда для фазовых траекторий, на которых значения x 1 , x 2 остаются достаточно большими, вполне удовлетворительной будет детерминистическая модель. Но если в некоторой точке фазовой траектории какая–либо переменная не очень велика, то существенное значение могут приобрести случайные флюктуации. Они приводят к тому, что изображающая точка переместится на одну из осей, что означает вымирание соответствующего вида.

Таким образом, стохастическая модель оказывается неустойчивой, так как стохастический “дрейф” рано или поздно приводит к вымиранию одного из видов. В такого рода модели хищник в конечном счете вымирает, это может произойти либо случайно, либо вследствие того, что сначала элиминируется популяция его жертвы. Стохастическая модель системы хищник – жертва хорошо объясняет эксперименты Гаузе (Гаузе, 1934), в которых инфузория Paramettum candatum служила жертвой для другой инфузории Didinium nasatum – хищника. Ожидавшиеся согласно детерминистическим уравнениям (9.7) равновесные численности в этих экспериментах составляли примерно всего по пять особей каждого вида, так что нет ничего удивительного в том, что в каждом повторном эксперименте довольно быстро вымирали либо хищники, либо жертвы (а за ними и хищники) Результаты экспериментов представлены на рис. 9.5.

Рис . 9.5. Рост Parametium caudatum и хищной инфузории Dadinium nasutum . Из : Gause G.F. The struggle for existence. Baltimore , 1934

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

Обобщенные модели взаимодействия двух видов

Было предложено большое число моделей, описывающих взаимодействие видов, правые части уравнений которых представляли собой функции численностей взаимодействующих популяций. Решался вопрос о выработке общих критериев, позволяющих установить, какого вида функции могут описать особенности поведения временного численности популяции, в том числе устойчивые колебания. Наиболее известные из этих моделей принадлежат Колмогорову (1935, переработанная статья ‑ 1972) и Розенцвейгу (1963).

(9.12)

В модель заложены следующие предположения:

1) Хищники не взаимодействуют друг с другом, т.е. коэффициент размножения хищников k 2 и число жертв L , истребляемых в единицу времени одним хищником, не зависит от y.

2) Прирост числа жертв при наличии хищников равен приросту в отсутствие хищников минус число жертв, истребляемых хищниками. Функции k 1 (x ), k 2 (x ), L (x ), ‑ непрерывны и определены на положительной полуоси x , y ³ 0.

3) dk 1 /dx< 0. Это означает, что коэффициент размножения жертв в отсутствие хищника монотонно убывает с возрастанием численности жертв, что отражает ограниченность пищевых и иных ресурсов.

4) dk 2 / dx> 0, k 2 (0) < 0 < k 2 (¥ ). С ростом численности жертв коэффициент размножения хищников монотонно убывает с возрастанием численности жертв, переходя от отрицательных значений, (когда нечего есть) к положительным.

5) Число жертв, истребляемых одним хищником в единицу времени L (x )> 0 при N> 0; L (0)=0.

Возможные типы фазовых портретов системы (9.12) представлены на рис. 9.6:

Рис. 9.6. Фазовые портреты системы Колмогорова (9.12), описывающей взаимодействие двух видов при разных соотношениях параметров. Пояснения в тексте.

Стационарные решения (их два или три) имеют следующие координаты:

(1). ` x=0; ` y=0 .

Начало координат при любых значениях параметров представляет собой седло(рис. 9.6 а-г).

(2). ` x=A, ` y=0. (9.13)

A определяется из уравнения:

k 1 (A )=0.

Стационарное решение (9.13) ‑ седло, если B< A (рис. 9.6 а , б , г ), B определяется из уравнения

k 2 (B )=0

Точка (9.13) помещается в положительном квадранте, еслиB>A. Это ‑ устойчивый узел.

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

(3). ` x=B, ` y=C. (9.14)

Величина С определяется из уравнений:

Точка (9.14) – фокус (рис.9.6 а ) или узел (рис.9.6 г ), устойчивость которых зависит от знакавеличины s

s 2 = – k 1 (B) – k 1 (B )B+L (B )C.

Если s >0, точка устойчива, если s <0 ‑ точка неустойчива, и вокруг нее могут существовать предельные циклы (рис. 9.6 б )

В зарубежной литературе чаще рассматривается сходная модель, предложенная Розенцвейгом и Макартуром (1963):

(9.15)

где f (x ) - скорость изменения численности жертв x в отсутствие хищников, Ф(x,y ) ‑ интенсивность хищничества, k ‑ коэффициент, характеризующий эффективность переработки биомассы жертвы в биомассу хищника, e - смертность хищника.

Модель (9.15) сводится к частному случаю модели Колмогорова (9.12) при следующих предположениях:

1) численность хищника ограничивается только численностью жертвы,

2) скорость, с которой данная особь хищника поедает жертву, зависит только от плотности популяции жертв и не зависит от плотности популяции хищников.

Тогда уравнения (9.15) принимают вид.

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

Модельвзаимодействиядвух видов насекомых (M acArthur, 1971)

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

(9.16)

Здесь x,y - биомассы двух видов насекомых. Трофические взаимодействия описываемых в этой модели видов – весьма сложные. Этим определяется вид многочленов в правых частях уравнений.

Рассмотрим правую часть первого уравнения. Насекомые вида х поедают личинок вида у (член + k 3 y), но взрослые особи вида у поедают личинок вида х при условии высокой численности видов х или у или обоих видов (члены – k 4 xy, – у 2 ). При малых х смертность вида х выше, чем его естественный прирост (1–k 1 +k 2 x–x 2 < 0 при малых х). Во втором уравнении член k 5 отражает естественный прирост вида y; –k 6 y – самоограничение этого вида, –k 7 x – поедание личинок вида у насекомыми вида х, k 8 xy – прирост биомассы вида у за счет поедания взрослыми насекомыми вида у личинок вида х.

На рис. 9.7 представлен предельный цикл, являющийся траекторией устойчивого периодического решения системы (9.16).

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

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

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

Модель А.Д.Базыкина

Теоретический анализ моделей взаимодействий видов наиболее исчерпывающе проведен в книге А.Д.Базыкина “Биофизика взаимодействующих популяций” (М., Наука, 1985).

Рассмотрим одну из изученных в этой книге моделей типа хищник-жертва.

(9.17)

Система (9.17) является обобщением простейшей модели хищник-жертва Вольтерра (5.17) с учетом эффекта насыщения хищников. В модели (5.17) предполагается, что интенсивность выедания жертв линейно растет с ростом плотности жертв, что при больших плотностях жертв не соответствует реальности. Для описания зависимости рациона хищника от плотности жертв могут быть выбраны разные функции. Наиболее существенно, чтобы выбранная функция с ростом x асимптотически стремилась к постоянному значению. В модели (9.6) использовалась логистическая зависимость. В модели Базыкина в роли такой функции выбрана гипербола x /(1+px ). Вспомним, что такой вид имеет формула Моно, описывающая зависимость скорость роста микроорганизмов от концентрации субстрата. Здесь в роли субстрата выступает жертва, а в роли микроорганизмов – хищник.

Система (9.17) зависит от семи параметров. Число параметров можно уменьшить с помощью замены переменных:

x ® (A/D )x; y ® (A/D )/y;

t ® (1/A )t; g (9.18)

и зависит от четырех параметров.

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

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

Построение полного параметрического портрета производится в виде набора “срезов” (проекций) параметрического портрета малой размерности при фиксированных значениях некоторых из параметров.

Параметрический портрет системы (9.18) при фиксированных g и малых e представлен на рис.9.8. Портрет содержит 10 областей с различным типом поведения фазовых траекторий.

Рис. 9.8. Параметрический портрет системы (9.18) при фиксированных g

и малых e

Поведение системы при различных соотношениях параметров может быть существенно различным (рис.9.9). В системе возможны:

1) одно устойчивое равновесие (области 1 и 5);

2) один устойчивый предельный цикл (области 3 и 8);

3) два устойчивых равновесия (область 2)

4) устойчивый предельный цикл и неустойчивое равновесие внутри него (области 6, 7, 9, 10)

5) устойчивый предельный цикл и устойчивое равновесие вне его (область 4).

В параметрических областях 7, 9, 10 область притяжения равновесия ограничивается неустойчивым предельным циклом, лежащим внутри устойчивого. Наиболее интересно устроен фазовый портрет, соответствующий области 6 на параметрическом портрете. Детально он изображен на рис. 9.10.

Область притяжения равновесия В 2 (заштрихована) представляет собой “улитку”, скручивающуюся с неустойчивого фокуса В 1 . Если известно, что в начальный момент времени система находилась в окрестности В 1 , то судить о том, придет ли соответствующая траектория в равновесие В 2 или на устойчивый предельный цикл, окружающий три точки равновесия С (седло), В 1 и В 2 можно лишь на основе вероятностных соображений.

Рис.9.10. Фазовый портрет системы 9.18 для параметрической области 6. Область притяжения В 2 заштрихована

На параметрическом портрете (9.7) имеются 22 различные бифуркационные границы, которые образуют 7 различных типов бифуркаций. Их изучение позволяет выявить возможные типы поведения системы при изменении ее параметров. Например, при переходе из области 1 в область 3 происходит рождение малого предельного цикла, или мягкое рождение автоколебаний вокруг единственного равновесия В. Аналогичное мягкое рождение автоколебаний, но вокруг одного из равновесий, а именно B 1 , происходит при пересечении границы областей 2 и 4. При переходе из области 4 в область 5 устойчивый предельный цикл вокруг точки B 1 “лопается” на петле сепаратрис и единственной притягивающей точкой остается равновесие B 2 и т.д.

Особый интерес для практики представляет, конечно, выработка критериев близости системы к бифуркационным границам. Действительно, биологам хорошо известно свойство “буферности”, или “гибкости”, природных экологических систем. Этими терминами обычно обозначают способность системы как бы поглощать внешние воздействия. Пока интенсивность внешнего воздействия не превышает некоторой критической величины, поведение системы не претерпевает качественных изменений. На фазовой плоскости это соответствует возвращению системы в устойчивое состояние равновесия или на устойчивый предельный цикл, параметры которого не сильно отличаются от первоначального. Когда же интенсивность воздействия превышает допустимую, система “ломается”, переходит в качественно иной режим динамического поведения, например просто вымирает. Это явление соответствует бифуркационному переходу.

Каждый тип бифуркационных переходов имеет свои отличительные особенности, позволяющие судить об опасности такого перехода для экосистемы. Приведем некоторые общие критерии, свидетельствующие о близости опасной границы. Как и в случае одного вида, если при уменьшении численности одного из видов происходит “застревание” системы вблизи неустойчивой седловой точки, что выражается в очень медленном восстановлении численности к начальному значению, значит, система находится вблизи критической границы. Индикатором опасности служит также изменение формы колебаний численностей хищника и жертвы. Если из близких к гармоническим колебания становятся релаксационными, причем амплитуда колебаний увеличивается, это может привести к потере устойчивости системы и вымиранию одного из видов.

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

Литература.

Колмогоров А.Н. Качественное изучение математических моделей динамики популяций. // Пороблемы кибернетики. М ., 1972, Вып .5.

MacArtur R. Graphycal analysis of ecological systems// Division of biology report Perinceton University. 1971

А.Д.Базыкин “Биофизика взаимодействующих популяций”. М., Наука, 1985.

В.Вольтерра: «Математическая теория борьбы за существование». М .. Наука , 1976

Gause G.F. The struggle for existence. Baltimore, 1934.

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

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

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

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

Пусть в отсутствие хищников популяция жертвы растет экспоненциаль­но:

dN/dt =r 1 N 1

где N -численность, а r, - удельная мгновенная скорость роста популя­ции жертвы. Если же хищники присутствуют, то они уничтожают особей жертвы со скоростью, которая определяется, во-первых, частотой встреч хищников и жертв, возрастающей по мере увеличения их численностей, и, во-вторых, эффективностью, с которой хищник обнаруживает и ловит свою жертву при встрече. Число жертв, встреченных и съеденных одним хищником N с, пропорционально эффективности охоты, которую мы выразим через коэффициент С 1; численности (плотности) жертвы N и времени, затраченному на поиски Т:

N C =C 1 NT (1)

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

В рассматриваемой модели С 1 является константой. Это означает, что число жертв, изъятых хищниками из популяции, линейно возрастает с увеличением ее плотности (так называемый функциональный ответ типа 1). Ясно, что общая скорость поедания жертв всеми особями хищника составит:

(3)

где Р - численность популяции хищника. Теперь мы можем записать уравнение роста популяции жертвы следующим образом:

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

(5)

где r 2 - удельная мгновенная смертность в популяции хищника.

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

Поскольку С 1 и s - константы, их произведение - это также константа, которую мы обозначим как С 2 . Тогда скорость роста популяции хищника будет определяться балансом рождаемости и смертности в соответствии с уравнением:

(6)

Уравнения 4 и 6 вместе образуют модель Лотки-Вольтерра.

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

Для популяции жертвы: при откуда

Таким образом, поскольку r, и С 1 , - константы, изоклиной для жертвы будет линия, на которой численность хищника (Р) является постоянной, т.е. параллельная оси абсцисс и пересекающая ось ординат в точке Р =r 1 / С 1 . Выше этой линии численность жертвы будет уменьшаться, а ниже- возрастать.

Для популяции хищника:

при откуда

Поскольку r 2 и С 2 - константы, изоклиной для хищника будет линия, на которой численность жертвы (N) является постоянной, т.е. перпендикуляр­ная оси ординат и пересекающая ось абсцисс в точке N = r 2 /С 2 . Слева от нее численность хищника будет уменьшаться, а справа - возрастать.

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

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

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

КОМПЬЮТЕРНАЯ МОДЕЛЬ «ХИЩНИК-ЖЕРТВА»

Казачков Игорь Алексеевич 1 , Гусева Елена Николаевна 2
1 Магнитогорский государственный технический университет им. Г.И. Носова, институт строительства, архитектуры и искусства, студент 5 курса
2 Магнитогорский государственный технический университет им. Г.И. Носова, институт энергетики и автоматизированный систем, кандидат педагогических наук, доцент кафедры бизнес-информатики и информационных технологий

COMPUTER MODEL «PREDATOR-VICTIM»

Kazatchkov Igor Alekseevich 1 , Guseva Elena Nikolaevna 2
1 Nosov Magnitogorsk State Technical University, Civil Engineering, Architecture and Arts Institute, student of the 5th course
2 Nosov Magnitogorsk State Technical University, Power Engineering and Automated Systems Institute, PhD in Pedagogical Science, Associate Professor of the Business Computer Science and Information Technologies Department

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

В настоящий момент подобные формы используется для изучения окружающей нас среды, а когда требуется изучить какую-либо из ее областей, то применяют математическое моделирование. Данная модель дает возможность спрогнозировать влияние тех или иных факторов на объект изучения. В свое время был предложен тип «хищник – жертва» такими учеными как: Т. Мальтусом (Malthus 1798, Мальтус 1905), Ферхюльстом (Verhulst 1838), Пирлом (Pearl 1927, 1930), а также А. Лотки (Lotka 1925, 1927) и В. Вольтерры (Volterra 1926).Эти модели воспроизводят периодический колебательный режим, возникающий в результате межвидовых взаимодействий в природе.

Одним из основных методов познания является моделировка. Помимо того, что в нем можно спрогнозировать изменения, происходящие в окружающей среде, к тому же помогает найти оптимальный способ решения проблемы. Уже давно в экологии используют математические модели, для того чтобы установить закономерности, тенденции развития популяций, помогают выделить суть наблюдений. Макет может служить образцом поведения, объекта.

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

Рассмотрим, что же такое объект моделирования. «Объектом является замкнутая среда обитания, где происходит взаимодействие двух биологических популяций: хищников и жертв. Процесс роста, вымирания и размножения происходит непосредственно на поверхности среды обитания. Питание жертв происходит за счет тех ресурсов, которые присутствуют в данной среде, а питание хищников происходит за счет жертв. При этом питательные ресурсы могут быть как возобновляемые, так и не возобновляемые.

В 1931 году Вито Вольтеррой были выведены следующие законы отношения хищник-жертва.

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

Закон сохранения средних величин – средняя численность каждого вида постоянна, независимо от начального уровня, при условии, что специфические скорости увеличения численности популяций, а также эффективность хищничества постоянны.

Закон нарушения средних величин – при сокращении обоих видов пропорционально их числу, средняя численность популяции жертвы растет, а хищников – падает.

Модель хищник-жертва – это особая взаимосвязь хищника с жертвой, в результате которой выигрывают оба. Выживают наиболее здоровые и приспособленные особи к условиям среды обитания, т.е. все это происходит благодаря естественному отбору. В той среде где нет возможности для размножения, хищник рано или поздно уничтожит популяцию жертвы, в последствии чего вымрет и сам» .

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

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

Рассмотрим компьютерную модель размножения млекопитающих на примере модели Лотки-Вольтерры. Пусть на некоторой территории обитают два вида животных: олени и волки. Математическая модель изменения численности популяций в модели Лотки-Вольтерры:

Начальное число жертв - xn, число хищников - yn.

Параметры модели:

P1– вероятность встречи с хищником,

P2– коэффициент роста хищников за счет жертв,

d – коэффициент смертности хищников,

a – коэффициент прироста численности жертв.

В учебной задаче были заданы такие значения: численность оленей равнялось 500, численности волков равна 10, коэффициент прироста оленей равен 0,02, коэффициент прироста численности волков равен 0,1, вероятность встречи с хищником 0,0026, коэффициент роста хищников за счет жертв 0,000056. Данные рассчитаны на 203 года.

Исследуем влияние коэффициент прироста жертв на развитие двух популяций, остальные параметры оставим без изменений. На схеме 1 наблюдается увеличение численности жертвы и затем, с некоторым опозданием наблюдается прирост хищников. Затем хищники выбивают жертв, число жертв резко падает и вслед за ним уменьшается число хищников (рис. 1).

Рисунок 1. Численность популяций при низкой рождаемости у жертв

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

Рисунок 2.Численность популяций при средней рождаемости у жертв

Рассмотрим как изменится динамика популяций при высоком значении коэффициента рождаемости жертвы а=1,13. На рис. 3 наблюдается резкое увеличение численности обеих популяций с последующим вымиранием, как жертвы, так и хищника. Это происходит за счет того, что численность популяции жертв увеличилось до такого количества, что стали заканчиваться ресурсы, вследствие чего происходит вымирание жертвы. Вымирание хищников происходит из-за того, что сократилось количество жертв и у хищников закончились ресурсы для существования.

Рисунок 3.Численность популяций при высокой рождаемости у жертв

Исходя из анализа данных компьютерного эксперимента, можно сделать выводы о том, что компьютерное моделирование позволяет нам прогнозировать численность популяций, изучать влияние различных факторов на популяционную динамику. В приведенном примере мы исследовали модель «хищник-жертва», влияние коэффициента рождаемости жертв на численность оленей и волков. Небольшой прирост популяции жертв приводит к небольшому увеличению жертв, которую через некоторый период уничтожают хищники. Умеренный прирост популяции жертв приводит к увеличению численности обеих популяций. Высокий прирост популяции жертв приводит сначала к быстрому росту популяции жертв, это влияет на увеличение роста хищников, но затем расплодившиеся хищники быстро уничтожают популяцию оленей. В итоге оба вида вымирают.

к договору от ___.___,20___ г. об оказании платных образовательных услуг

Министерство образования и науки Российской Федерации

Лысьвенский филиал

Пермского Государственного технического университета

Кафедра ЕН

Курсовая работа

по дисциплине «Моделирование систем»

тема: Система хищник-жертва

Выполнил:

Студент гр. БИВТ-06

------------------

Проверил преподаватель:

Шестаков А. П.

Лысьва, 2010 г.

Реферат

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

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

Кроме многоклеточных животных, в роли хищников могут выступать протисты, грибы и высшие растения.

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

В процессе совместной эволюции хищники и жертвы приспосабливаются друг к другу. У хищников появляются и развиваются средства обнаружения и атаки, а у жертв - средства скрытности и защиты. Поэтому наибольший вред жертвам могут нанести новые для них хищники, с которыми те не вступали ещё в «гонку вооружений».

Хищники могут специализироваться на одном–нескольких видах для добычи, это делает их в среднем более успешными в охоте, но повышает зависимость от данных видов.

Система хищник-жертва.

Взаимодействие хищник-жертва - основной тип вертикальных взаимоотношений организмов, при котором по пищевым цепям передаются вещество и энергия.

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

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

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

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

• убегать от хищника (в этом случае в результате соревнования повышается подвижность и жертв, и хищников, что особенно характерно для степных животных, которым негде прятаться от преследователей);
• приобретать защитную окраску (<притворяться> листьями или сучками) или, напротив, яркий (например, красный) цвет, предупреждающий хищника о горьком вкусе;
• прятаться в укрытия;
• переходить к мерам активной обороны (рогатые травоядные, колючие рыбы), часто совместной (птицы-жертвы коллективно отгоняют коршуна, самцы оленей и сайгаков занимают <круговую оборону> от волков и т. д.).