Исследуем модель. Построим уравнения, в которых не будет F.
Умножим первое уравнение на γ2, второе - на γ1 и вычтем из первого второе.
Подчиним начальным условиям:
Будем предполагать, что
Сформулируем свойства функции F, которыми будем пользоваться при дальнейшем анализе:
F = 0, если N1 = N2 = 0.
F > 0 для любых N1, N2 > 0
При неограниченном возрастании N1 и N2 F тоже неограниченно возрастает.
F монотонно возрастает по каждому своему аргументу:
Пусть N1* - корень уравнения
В момент, когда достигается численность первого вида, равная N1*, численность второго вида равна N2. Может ли N1 увеличиваться дальше?
Значит, не может. Значит, N1 - ограниченная функция. Аналогично, N2 - ограниченная. Следовательно, невозможен неограниченный рост численности.
Правая часть (3) стремится к бесконечности в силу (4). Левая часть - отношение двух ограниченных функций, значит, знаменатель стремится к нулю, т.е. первый вид вытесняет второй.
Положения равновесия
Решение системы (1) вида
называется стационарным решением. (N1*, N2*) - положение равновесия. Найдем их из решений системы
Система распадается на 4 системы:
1. N1* = N2* = 0. Обозначим эту точку покоя - P0(0 ,0).
2. N1* = 0, N2*:
В силу монотонности решение (5) существует всегда. Значит, вторая точка покоя - P1(0, N2*), где N2* - корень уравнения (5).
3. N2* = 0, N1*:
Тоже всегда есть решение. Точка покоя - P2(N1*, 0), где N1* - корень уравнения (6).
Вывод: система (1) имеет три точки покоя.
Характер устойчивости
Исследование на асимптотическую устойчивость проведем с помощью линеаризации.
Рассмотрим построение линеаризованной системы для (1) в окрестности произвольного положения равновесия P*.
Правые части (1) необходимо заменить линеаризованными функциями:
Аналогично для правой второго уравнения системы (1).
Получаем:
Отклонение от положения равновесия по соответствующей координате:
Линеаризованная система имеет вид:
Итак, имеем три точки покоя и линеаризованную систему (7). Исследуем на устойчивость каждую точку отдельно.
Первая точка покоя. P0(0, 0). Подставляя в (7), получаем:
Решение:
x0, y0 - начальное отклонение. Точка P0 неустойчива, малые отклонения от положения равновесия со временем нарастают. Т.к. собственные значения системы
 , |
то P0 - неустойчивый узел.
Вторая точка покоя. P1(0, N2*). Соответствующая линеаризованная система:
 , |
Матрица системы имеет треугольный вид, значит собственные значения совпадают с коэффициентами главной диагонали:
 , |
Определим знак первого собственного значения λ1:
 , |
Получили, что собственные значения имеют разный знак. P1 - "седло". "Седло" всегда неустойчиво.
У "седла" две сепаратрисы - прямые, найдем их. Будем искать в виде: y = kx.
Выразим k1:
Числитель дроби отрицателен, знаменатель положителен. Значит k1 < 0. Нашли только одно конечное значение k, значит
Получаем две сепаратрисы:
Сепаратрисы разбивают фазовую плоскость на 4 области. Вся плоскость заполняется траекториями. Как определить направление движения по траектории? Пусть x = 0, как изменится y? Подставим x = 0 во второе уравнение. Получим решением экспоненциальную функцию, стремящуюся к нулю. Значит y убывает.
Третья точка покоя. P2(N1*, 0). Соответствующая линеаризованная система имеет вид:
Матрица системы имеет треугольный вид. Ниже главной диагонали - нули, значит, собственные значения системы имеют вид:
Получаем, что собственные значения вещественны и отрицательны. Значит, P2 - устойчивый узел.
Первый вид (если он есть) вытесняет второй. Количество особей первого вида стремится к некоторой константе. Чем больше величина ε/γ, тем меньше нехватка пищи сказывается на поведении особей. В нашем случае
Условие (4) может быть нарушено. Например, ε1γ2 - γ1ε2 < 0, но это тоже самое, что изменить нумерацию видов. Ничего не изменится, только второй вид, наоборот, будет вытеснять первый.
Рассмотрим случай, когда ε1γ2 = γ1ε2, т.е. оба вида реагируют на нехватку пищи одинаково. Изменится соотношение (3):
- уравнение фазовой траектории.
Вспомним, как мы искали положения равновесия. Четвертый случай был
- это целая линия точек.
Рассмотрим случай, когда F(N1, N2) = N1 + N2. Эта функция удовлетворяет условиям.
Нулевое положение так и останется неустойчивым. Весь отрезок состоит из устойчивых точек, но нет асимптотической устойчивости.
Пересечение вертикальных и горизонтальных изоклин дает точку покоя. Вертикальные изоклины найдем из условия N1* = 0, горизонтальные - N2* = 0:
Если F = N1 + N2, то
|