Что можно приготовить из кальмаров: быстро и вкусно

Допущения:

1. Среда однородная.

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

3. Пренебрегаем случайными флуктуациями.

4. Взаимодействие мгновенное.

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

Однако их результаты зачастую противоречат друг другу:

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

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

Построенная в 20-х годах ХХ века Вито Вольтера модель сообщества «хищник-жертва» объясняет многие из этих особенностей.

Это первый успех математической экологии.

При рассмотрении этой системы рассмотрим вопросы устойчивости: условия устойчивости и механизмы устойчивости.

Классическая модель Вольтерра

Численность жертвы,

Численность хищников.

Дополнительные допущения.

1. Единственным лимитирующим фактором, ограничивающим размножение жертв, является давление на них со стороны хищников. Ограниченность ресурсов среды для жертвы не учитывается (как в модели Мальтуса).

2. Размножение хищников ограничивается количеством добытой им пищи (количеством жертв).

− коэффициент естественного прироста жертвы;

− коэффициент естественной смертности хищника;

− количество (биомасса) жертв, потребляемых одним хищником за единицу времени (трофическая функция);

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

Уравнения системы «хищник-жертва»

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

Этот тип характерен для беспозвоночных и некоторых видов хищных рыб.

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

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

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

Классическая модель Вольтерра:

Начальные условия

Система (2) является автономной, т.к. не имеет в правой части. Изменение состояния системы изображается на фазовой плоскости и является решением уравнения

Найдем точки покоя системы (2).

Нетривиальная точка покоя системы (4) имеет вид

Определим характер точки покоя (5).

Сделаем замену

Раскроем скобки и получим систему

Отбросив нелинейные члены, получим систему

Характеристическое уравнение имеет вид

Корни - чисто мнимые числа. Точка покоя – центр. В исходных переменных фазовые траектории имеют вид

Стрелки указывают направление изменения состояния системы со временем.

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

Фазовый портрет решения имеет вид спирали:

В системе «хищник-жертва» возникают затухающие колебания. Численности жертв и хищников стремятся к своим равновесным значениям (8).

Графики зависимости численностей видов.

Модель ситуации типа «хищник-жертва»

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

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

• 1. Жертва может найти достаточно пищи для пропитания;
• 2. При каждой встрече жертвы с хищником последний убивает жертву.

Будем для определённости называть их карасями и щуками. Пусть

состояние системы определяется величинами x(t) и y(t) - количеством карасей и щук в момент г. Чтобы получить математические уравнения, которые приближенно описывают динамику (изменение во времени) популяции, поступим следующим образом.

Как и в предыдущей модели роста популяции (см. п. 1.1) для жертв имеем уравнение

где а > 0 (рождаемость превышает смертность)

Коэффициент а прироста жертв зависит от количества хищников (уменьшается с их увеличением). В простейшем случае а- а - fjy (а>0,р>0). Тогда для численности популяции жертв имеем дифференциальное уравнение

Для популяции хищников имеем уравнение

где b >0 (смертность превышает рождаемость).

Коэффициент b убывания хищников уменьшается, если имеются жертвы, которыми можно питаться. В простейшем случае можно принять b - у -Sx {у > 0, S > 0). Тогда для численности популяции хищников получим дифференциальное уравнение

Таким образом, уравнения (1.5) и (1.6) представляют собой математическую модель рассматриваемой задачи взаимодействия популяций. В этой модели переменные х,у - состояние системы, а коэффициенты характеризуют ее структуру. Нелинейная система (1.5), (1.6) и есть модель Вольтера-Лотки.

Уравнения (1.5) и (1.6) следует дополнить начальными условиями - заданными значениями начальных популяций.

Проведем теперь анализ построенной математической модели.

Посгроим фазовый портрет системы (1.5), (1.6) (по смыслу задачи х > 0, v >0). Разделив уравнение (1.5) на уравнение (1.6), получим уравнение с разделяющимися переменными

Игшлрируя это уравнение, будем иметь

Соотношение (1.7) даёт уравнение фазовых траекторий в неявном виде. Система (1.5), (1.6) имеет стационарное состояние определяемое из

Из уравнений (1.8) получим (т.к. л* Ф 0, у* Ф 0)

Равенства (1.9) определяют на фазовой плоскости положение равновесия (точку О) (Рис. 1.6).

Направление движения по фазовой траектории можно определить из таких соображений. Пусть карасей мало. г.е. х ~ 0, тогда из уравнения (1.6) у

Все фазовые траектории (за исключением точки 0) замкнутые кривые, охватывающие положение равновесия. Состоянию равновесия соответствует неизменное количество х« и у« карасей и щук. Караси размножаются, щуки их едят, вымирают, но число тех и дрч их не меняется. "Замкнутым фазовым траекториям соответствует периодическое изменение численности карасей и щук. Причём то, по какой траектории движется фазовая точка, зависит от начальных условий. Рассмотрим, как меняется состояние вдоль фазовой траектории. Пусть точка находится в положении А (рис. 1.6). Здесь карасей мало, щук много; щукам есть нечего, и они постепенно вымирают и почти

совсем исчезают. Но и количество карасей тоже уменьшается почти до нуля и

только потом, когда щук стало меньше, чем у , начинается прирост количества карасей; скорость их прироста увеличивается и их число увеличивается - так происходит примерно до точки В. Но увеличение числа карасей приводит к торможению процесса вымирания шук и их число начинает расти (пищи стало больше) - участок ВС. Далее щук много, они едят карасей и почти всех съедают (участок CD). После этого щуки снова начинают вымирать и процесс повторяется с периодом примерно в 5-7 лет. На рис. 1.7 качественно построены кривые изменения численности карасей и щук в зависимости от времени. Максимумы кривых чередуются, причём максимумы численности щук отстают от максимумов карасей.

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

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

Отметим еще один интересный вывод, который можно сделать из Рис. 1.6. Если в точке Р происходит быстрый отлов щук (в другой терминологии - отстрел волков), то система «перепрыгивает » в точку Q, и дальнейшее движение происходит по замкнутой траектории меньшего размера, что интуитивно ожидаемо. Если же уменьшить число щук в точке R, то система перейдет в точку S, и дальнейшее движение будет происходить по траектории большего размера. Амплитуды колебаний увеличатся. Это противоречит интуиции, но как раз объясняет такое явление: в результате отстрела волков их численность увеличивается со временем. Таким образом, важным в этом случае является выбор момента отстрела.

Предположим, что две популяции насекомых (например, тля и божья коровка, которая есть тлю) находились в естественном равновесии х-х*,у = у* (точка О на Рис. 1.6). Рассмотрим влияние разового применения инсектицида, который убивает х> 0 из жертв и у > 0 из хищников, не уничтожая их полностью. Уменьшение численности обеих популяций приводит к тому, что изображающая точка из положения О «перескочит» ближе к началу координат, где х > 0, у 0 (Рис. 1.6) Отсюда следует, что в результате действия инсектицида, призванного уничтожать жертв (тлю), число жертв (тли) увеличивается, а число хищников (божьих коровок) уменьшается. Получается, что численность хищников может стать настолько малой, что им будет фозить полное исчезновение но другим причинам (засуха, болезни и т.д.). Таким образом, применение инсектицидов (если только они не уничтожают вредных насекомых практически полностью) в конечном счёте приводит к увеличению популяции тех насекомых, численность которых находилась под контролем других насекомых-хищников. Такие случаи описаны в книгах по биологии.

В общем случае коэффициент прироста количества жертв а зависит и от Л" и от у: а = а(х, у) (из-за наличия хищников и из-за ограничений на пищу).

При малом изменении модели (1.5), (1.6) к правым частям уравнений добавляются малые члены (учитывающие, например, конкуренцию карасей за пищу и щук за карасей)

здесь 0 f.i « 1.

В таком случае вывод о периодичности процесса (возвращении системы к исходному состоянию), справедливый для модели (1.5), (1.6), теряет силу. В зависимости от вида малых поправок/ и g возможны ситуации, показанные на Рис. 1.8.

В случае (1) равновесное состояние О устойчиво. При любых других начальных условиях через достаточно большое время устанавливается именно оно.

В случае (2) система «идёт в разнос». Стационарное состояние неустойчиво. Такая система в конце концов попадает в такую область значений х и у, что модель перестаёт быть применимой.

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

В.И. Арнольд называет модель Вольтерра-Лотки жесткой, т.к. её малое изменение может привести к выводам, отличным от приведенных выше. Для суждения о том, какая из ситуаций, указанных на Рис. 1.8, реализуется в данной системе, совершенно необходима дополнительная информация о системе (о виде малых поправок/ и g ).

Еще в 20-х гг. А. Лотка, а несколько позднее независимо от него В.Вольтера предложили математические модели, описывающие сопряженные колебания численности хищника и жертвы.

Модель состоит из двух компонентов:

С – численность хищников; N – численность жертв;

Предположим, что в отсутствии хищников популяция жертвы будет расти экспоненциально: dN/dt = rN. Но жертвы уничтожаются хищниками со скоростью, которая определяется частотой встреч хищника с жертвой, а частота встреч возрастает по мере увеличения численности хищника (С) и жертвы (N). Точное число встреченных и успешно съеденных жертв будет зависеть от эффективности, с которой хищник находит и ловит жертву, т.е. от а’ – «эффективности поиска» или «частоты нападений». Т.о., частота «успешных» встреч хищника с жертвой и, следовательно, скорость поедания жертв будет равна а’СN и в целом: dN/dt = rN – a’CN (1*).

В отсутствии пищи отдельные особи хищника теряют в весе, голодают и гибнут. Предположим, что в рассматриваемой модели численность популяции хищника в отсутствии пищи вследствие голодания будет уменьшаться экспоненциально: dC/dt = - qC, где q – смертность. Гибель компенсируется рождением новых особей со скоростью, которая, как полагают в этой модели, зависит от двух обстоятельств:

1) скорости потребления пищи, a’CN;

2)эффективности (f), с которой эта пища переходит в потомство хищника.

Т.о., рождаемость хищника равна fa’CN и в целом: dC/dt = fa’CN – qC (2*). Уравнения 1* и 2* составляют модель Лотки-Вольтера. Свойства этой модели можно исследовать, построить изоклины линии, соответствующие постоянной численности популяции, с помощью таких изоклин определяют поведение взаимодействующих популяций хищник-жертва.

В случае с популяцией жертвы: dN/dt = 0, rN = a’CN, или C = r/a’. Т.к. r и a’ = const, изоклиной для жертвы будет линия, для которой величина С является постоянной:

При низкой плотности хищника (С), численность жертвы (N) возрастает, наоборот – снижается.

Точно также для хищников (уравнение 2*) при dC/dt = 0, fa’CN = qC, или N = q/fa’, т.е. изоклиной для хищника будет линия, вдоль которой N постоянна: При высокой плотности жертвы, численность популяции хищника растет, а при низкой – снижается.

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

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

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

5) ПОПУЛЯЦИОННЫЕ СТРАТЕГИИ ХИЩНИКА И ЖЕРТВЫ

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

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

б) истинные хищники . Хищник преследует и убивает жертву;

в) пастбищники . Эти хищники используют жертву многократно, например, оводы или слепни.

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

Истинные хищники делятся на

«жнецов», которые питаются обильными ресурсами (н, планктонные рыбы и даже усатый кит), и «охотников», которые добывают менее обильный корм. В свою очередь

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

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

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

Приобретать защитную окраску («притворяться» листьями или сучками) или, напротив, яркий цвет, Н.: красный, предупреждающий хищника о горьком вкусе. Общеизвестно изменение окраски зайца в разные времена года, что позволяет ему маскироваться летом в траве, а зимой на фоне белого снега. Адаптивное изменение окраски может происходить в разных стадиях онтогенеза: детеныши тюленей белые (цвет снега), а взрослые особи черные (цвет скалистого побережья);

Распространяться группами, что делает их поиск и промысел для хищника более энергоемким;

Прятаться в укрытия;

Переходить к мерам активной обороны (травоядные, имеющие рога, колючие рыбы), иногда совместной (овцебыки могут занимать «круговую оборону» от волков и т.д.).

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

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

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

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

Между биомассами хищников и жертв устанавливается достаточно устойчивое соотношение. Так, Р. Риклефс приводит данные о том, что соотношение биомасс хищника и жертвы колеблется в пределах 1:150 - 1:300. В разных экосистемах умеренной зоны США на одного волка приходится 300 мелких белохвостых оленей (вес 60 кг), 100 крупных оленей вапити (вес 300 кг) или 30 лосей (вес 350). Такая же закономерность выявлена в саваннах.

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

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

1.Постоянный и благоприятный уровень температуры и влажности.

2.Обилие пищи.

3.Защита от неблагоприятных факторов.

4.Агрессивный химический состав среды обитания (пищеварительные соки).

1.Наличие двух сред обитания: среда первого порядка - организм хозяина, среда второго порядка - внешняя среда.

Федеральное агентство по образованию

Государственное образовательное учреждение

высшего профессионального образования

«Ижевский государственный технический университет»

Факультет «Прикладная математика»

Кафедра «Математическое моделирование процессов и технологий»

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

по дисциплине «Дифференциальные уравнения»

Тема: «Качественное исследование модели хищник-жертва»

Ижевск 2010

ВВЕДЕНИЕ

1. ПАРАМЕТРЫ И ОСНОВНОЕ УРАВНЕНИЕ МОДЕЛИ «ХИЩНИК-ЖЕРТВА»

2.2 Обобщенные модели Вольтера типа «хищник-жертва».

3. ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ МОДЕЛИ «ХИЩНИК-ЖЕРТВА»

ЗАКЛЮЧЕНИЕ

СПИСОК ЛИТЕРАТУРЫ

ВВЕДЕНИЕ

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

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

Одной из первых работ в области математической экологии была работа А.Д. Лотки (1880 - 1949), который первый описал взаимодействие различных популяций, связанных отношениями хищник - жертва. Большой вклад в исследование модели хищник -жертва внесли В. Вольтерра (1860 - 1940), В.А. Костицин (1883-1963) В настоящее время уравнения описывающие взаимодействие популяций, называются уравнениями Лотки - Вольтерра.

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

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

1. ПАРАМЕТРЫ И ОСНОВНОЕ УРАВНЕНИЕ МОДЕЛИ ХИЩНИК-ЖЕРТВА

Попытки математического моделирования динамики как отдельных биологических популяций, так и сообществ, включающих взаимодействующие популяции различных видов, предпринимались давно. Одна из первых моделей роста изолированной популяции (2.1) была предложена еще в 1798 г. Томасом Мальтусом:

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

N - численность популяции;

Разность между коэффициентами рождаемости и смертности.

Интегрируя это уравнение получаем:

, (1.2)

где N(0) – численность популяции в момент t = 0. Очевидно, что модель Мальтуса при > 0 дает бесконечный рост численности, что никогда не наблюдается в природных популяциях, где ресурсы, обеспечивающие этот рост, всегда ограничены. Изменения численности популяций растительного и животного мира нельзя описывать простым законом Мальтуса, на динамику роста влияют многие взаимосвязанные причины – в частности, размножение каждого вида саморегулируется и видоизменяется так, чтобы этот вид сохранялся в процессе эволюции.

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

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

,

где - численность популяции;

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

2. КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ ЭЛЕМЕНТАРНОЙ МОДЕЛИ «ХИЩНИК-ЖЕРТВА»

2.1 Модель трофического взаимодействия по типу «хищник-жертва»

Рассмотрим модель трофического взаимодействия по типу «хищник-жертва», построенную В. Вольтерром. Пусть имеется система, состоящая из двух видов, из которых один поедает другой.

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

Коэффициент прироста жертвы;

Коэффициент прироста хищника;

Численность популяции жертвы;

Численность популяции хищника;

Коэффициент естественного прироста жертвы;

Скорость потребления жертвы хищником;

Коэффициент смертности хищника в отсутствие жертвы;

Коэффициент «переработки» хищником биомассы жертвы в собственную биомассу.

Тогда динамика численности популяций в системе хищник - жертва будет описываться системой дифференциальных уравнений (2.1):

(2.1)

где все коэффициенты положительные и постоянные.

Модель имеет равновесное решение (2.2):

По модели (2.1) доля хищников в общей массе животных выражается формулой (2.3):

(2.3)

Анализ устойчивости состояния равновесия по отношению к малым возмущениям показал, что особая точка (2.2) является «нейтрально» устойчивой (типа «центр»), т. е. любые отклонения от равновесия не затухают, но переводят систему в колебательный режим с амплитудой, зависящей от величины возмущения. Траектории системы на фазовой плоскости имеют вид замкнутых кривых, расположенных на различных расстояниях от точки равновесия (рис. 1).

Рис. 1 – Фазовый «портрет» классической вольтерровой системы «хищник-жертва»

Разделив первое уравнение системы (2.1) на второе, получим дифференциальное уравнение (2.4) для кривой на фазовой плоскости .

(2.4)

Интегрируя данное уравнение получим:

(2.5)

где - постоянная интегрирования, где

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

(2.6)

В таком случае, подставив значения системы (2.6) в систему (2.1), будем иметь:

(2.7)

Умножив первое уравнение на , а второе - на и сложив их, получим:

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

Величина , как видно из (4.9), всегда больше нуля. Таким образом, не меняет знака, и вращение все время идет в одну сторону.

Интегрируя (2.9) найдем период:

Когда мало, то уравнения (2.8) и (2.9) переходят в уравнения эллипса. Период обращения в этом случае равен:

(2.11)

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

(2.12)

и проинтегрируем по периоду:

(2.13)

Так как подстановки от и в силу периодичности обращаются в нуль, средние по периоду оказываются равными стационарным состояниям (2.14):

(2.14)

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

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

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

(2.15)

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

(2.16)

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

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

Начиная с некоторого момента времени t, при котором , решение второго уравнения (2.15) также начинает убывать и при стремится к нулю. Таким образом, в случае оба вида исчезают.

2.1 Обобщенные модели Вольтера типа «хищник-жертва»

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

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

В 1936 г. А.Н. Колмогоров предложил использовать для описания динамики системы хищник - жертва следующую систему уравнении:

, (2.17)

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

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

Позднее в своей работе, Колмогоров исследовал подробно менее общую модель:

(2.18)

Различные частные случаи системы дифференциальных уравнений (2.18) исследовались многими авторами. В таблице приведены различные частные случаи функций , , .

Таблица 1 - Различные модели сообщества «хищник-жертва»

			Авторы
			Вольтерра-Лотка
			Гаузе
			Пислоу
			Холинг
			Ивлев
			Рояма
			Шимазу
			Мэй

математическое моделирование хищник жертва

3. ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ МОДЕЛИ ХИЩНИК-ЖЕРТВА

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

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

Заданы следующие начальные показатели:

Со временем число карасей и щук меняется, но так как рыбы в пруду много, то не будем различать 1020 карасей или 1021 и поэтому будем считать и непрерывными функциями времени t. Будем называть пару чисел (,) состоянием модели.

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

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

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

Проинтегрируем оба уравнения систему по t, которое будет изменяться от - начального момента времени, до , где T – период, за который происходят изменения в экосистеме. Пусть в нашем случае период равен 1 году. Тогда система принимает следующий вид:

;

;

Принимая = и = приведем подобные слагаемые, получим систему, состоящую из двух уравнений:

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

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

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

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

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

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

Интересно узнать, как именно в данной модели зависит динамика популяций от начальных параметров и насколько это соответствует действительности и здравому смыслу, причём увидеть это графически, не прибегая к сложным расчётам. Для этой цели на основе модели Вольтерра была создана программа в среде Mathcad14.

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

Полученные результаты отражают теоретические: хищники постепенно вымирают(Рис.1), а численность жертвы неограниченно возрастает(Рис.2).

Рис.1 Зависимость числа хищников от времени при отсутствии жертвы

Рис.2 Зависимость числа жертв от времени при отсутствии хищников

Как видно, в данных случаях система соответствует математической модели.

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

Таблица 1. Коэффициенты колебательного режима системы

Рис.3 Система при значении параметров из Таблицы 1

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

Положим теперь начальную численность хищника равную 200 особей при сохранении остальных параметров(Рис.4).

Таблица 2. Коэффициенты колебательного режима системы

Рис.4 Система при значении параметров из Таблицы 2

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

Зададим начальную численность хищника равную 100 особей, численность жертв 200 при сохранении остальных параметров(Рис.5).

Таблица 3. Коэффициенты колебательного режима системы

Рис.5 Система при значении параметров из Таблицы 3

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

Рассмотрим изменение других параметров системы. Пусть начальные численности соответствуют второму случаю. Увеличим коэффициент размножения жертв (Рис.6).

Таблица 4. Коэффициенты колебательного режима системы

Рис.6 Система при значении параметров из Таблицы 4

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

Теперь увеличим коэффициент прироста хищников (Рис.7).

Таблица 5. Коэффициенты колебательного режима системы

Рис.7 Система при значении параметров из Таблицы 5

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

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

Таблица 6. Коэффициенты колебательного режима системы

Рис.8 Система при значении параметров из Таблицы 6

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

И теперь уменьшим частоту поедания хищниками жертв:

Таблица 7. Коэффициенты колебательного режима системы

Рис.9 Система при значении параметров из Таблицы 7

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

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

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

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