Модель конкурентного взаимодействия
  
     Описание
Качественный анализ
Задачи и упражнения
Литература

Исследуем модель. Построим уравнения, в которых не будет F.
Умножим первое уравнение на γ2, второе - на γ1 и вычтем из первого второе.
(3)
Подчиним начальным условиям:
Будем предполагать, что
(4)
Сформулируем свойства функции 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)
    В силу монотонности решение (5) существует всегда. Значит, вторая точка покоя - P1(0, N2*), где N2* - корень уравнения (5).
    3. N2* = 0, N1*:
    (6)
    Тоже всегда есть решение. Точка покоя - P2(N1*, 0), где N1* - корень уравнения (6).
    4.
    Вывод: система (1) имеет три точки покоя.
    Характер устойчивости
    Исследование на асимптотическую устойчивость проведем с помощью линеаризации. Рассмотрим построение линеаризованной системы для (1) в окрестности произвольного положения равновесия P*.
    Правые части (1) необходимо заменить линеаризованными функциями:
    Аналогично для правой второго уравнения системы (1).
    Получаем:
    Отклонение от положения равновесия по соответствующей координате:
    Линеаризованная система имеет вид:
    (7)
    Итак, имеем три точки покоя и линеаризованную систему (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, то