Система Лотки-Вольтерры. Принцип Вольтерры: различия между версиями

Материал из sawiki
Перейти к навигации Перейти к поиску
(Прошу прощения, пока черновой вариант, я не нашел где и как тут более цивильно его сохранить кроме как засабмитить сюда)
 
 
(не показано 17 промежуточных версий этого же участника)
Строка 1: Строка 1:
== Популяционная система Коломогорова ==
+
== Система Лотки-Вольтерры "хищник — жертва" ==
В данном разделе нами рассматриваются межвидовые отношения типа "хищник - жертва". Опишем формально динамическую систему с непрерывным временем, моделирующую динамику численностей взаимодействующих популяций. Пусть динамика отдельной изолированной популяции задается уравнением:
+
Одной из первых рассмотренных математических моделей взаимодействующих популяций является система ОДУ, именуемая '''моделью [https://ru.wikipedia.org/wiki/Лотка,_Альфред_Джеймс Лотки]-[https://ru.wikipedia.org/wiki/Вольтерра,_Вито Вольтерры]''', которая описывает взаимодействие двух видов, один из которых является хищником, а другой — жертвой (например, экологическая система волки-зайцы).
<math>\dot u =  uf(u).</math>
 
Здесь $$u(t)$$ - численность популяции в момент времени $$t$$, $$f(u)$$ - коэффициент прироста, который зависит от текущей численности. Если имеется несколько взаимодействующих популяций, то коэффициент прироста каждой их них в самой общей ситуации зависит от численностей всех рассматриваемых популяций, и мы получаем систему дифференциальных уравнений:
 
\begin{equation}\label{syst1}
 
\begin{cases}
 
  \dot u_1 = u_1f_1(u_1, u_2, \dots, u_n),\\
 
  \dot u_2 = u_2f_2(u_1, u_2, \dots, u_n),\\
 
  \dots\\
 
  \dot u_n = u_nf_n(u_1, u_2, \dots, u_n),\\
 
\end{cases}
 
\end{equation}
 
где $$u = (u_1, u_2, \dots, u_n)$$ --- вектор численностей (или плотностей) взаимодействующих популяций, а $$f(u)$$ --- вектор-функция, описывающая систему, каждая компонента которой представляет собой коэффициент прироста соответствующего вида. Динамическую систему \eqref{syst1} будем называть '''системой популяционной динамики''', либо, как на неё часто ссылаются, '''популяционной системой Колмогорова'''.
 
  
== Система Лотки-Вольтерры "хищник-жертва" ==
+
[[Файл:Заяц-волк.png|мини| Модель Заяц &mdash; Волк попадалась в художественных произведениях уже полвека тому назад]]
Одной из первых математических моделей взаимодействующих популяций является система ОДУ, предложенная Вольтеррой в связи с попыткой объяснить колебания улова рыбы в Адриатическом море. Та же сама система была предложена Лоткой немного ранее. Модель Лотки-Вольтерры описывает вазимодействие двух видов, один из которых является хищником, а другой --- жертвой (например, экологическая система волки-зайцы)
 
  
[[Файл:Заяц-волк.png|обрамить|центр|Модель Заяц-Волк попадалась в художественных произведениях уже полвека тому назад]]
+
Если $$N(t)$$ &mdash; численность жертв, $$P(t)$$ &mdash; численность хищников в момент времени $$t$$, то модель Лотки-Вольтерры имеет вид:  
 
 
Если $$N(t)$$ - численность жертв, $$P(t)$$ - численность хищников в момент времени $$t$$, то модель Лотки-Вольтерры имеет вид:  
 
  
 
\begin{equation}\label{syst2}
 
\begin{equation}\label{syst2}
Строка 27: Строка 13:
 
\end{equation}
 
\end{equation}
  
где a, b, c, d --- положительные постоянные.
+
где $$a, ~b, ~c, ~d$$ &mdash; положительные постоянные.
\newline
+
 
Основные предположения, положенные в основу системы \eqref{syst2} характеризуются следующими гипотезами: в отсутствиии хищников жертвы размножаются неограниченно ($$\dot N = aN$$); хищники в отсутствии жертв вымирают ($$\dot P = -dP$$); слагаемые, пропорциональные NP, рассматриваются как превращение энергии одного источника в энергию другого (эффект влияния популяции хищников на популяцию жертв заключается в уменьшении относительной скорости прироста численности жертв на величину, пропорциональную численности хищников).
+
В основу системы \eqref{syst2} положены предположения, характеризующиеся следующими гипотезами:
\newline
+
отсутствиии хищников жертвы размножаются неограниченно ($$\dot N = aN$$).
Рассматривая систему \eqref{syst2} в качестве математической модели взаимодействущих популяций естественно считать фазовым пространством множество $$\mathbb{R}^2_+$$, которое является инвариантным, так как любая траектория, начинающаяся в $$\mathbb{R}^2_+$$ не может пересечь линии $$N = 0$$ и $$P = 0$$, являющиеся фазовыми кривыми.  
+
#Хищники в отсутствии жертв вымирают ($$\dot P = -dP$$).
\newline
+
#Слагаемые, пропорциональные $$NP$$, рассматриваются как превращение энергии одного источника в энергию другого (эффект влияния популяции хищников на популяцию жертв заключается в уменьшении относительной скорости прироста численности жертв на величину, пропорциональную численности хищников).
 +
 
 +
Так как мы смотрим на систему \eqref{syst2} как на математическую модель взаимодействущих популяций, естественно будет считать [https://sawiki.cs.msu.ru/index.php?title=Фазовые_и_интегральные_кривые._Фазовое_пространство&action=edit&redlink=1 фазовым пространством] множество $$\mathbb{R}^2_+$$, которое является [https://ru.wikipedia.org/wiki/Инвариантное_подпространство инвариантным], так как любая траектория, начинающаяся в $$\mathbb{R}^2_+$$ не может пересечь линии $$N = 0$$ и $$P = 0$$, являющиеся [https://sawiki.cs.msu.ru/index.php?title=Фазовые_и_интегральные_кривые._Фазовое_пространство&action=edit&redlink=1 фазовыми кривыми].  
 +
 
 
В безразмерных переменных система принимает вид:  
 
В безразмерных переменных система принимает вид:  
  
Строка 44: Строка 33:
 
где:
 
где:
 
<math>u(\tau) = \frac{d}{c}N(t), ~v(\tau) = \frac{bP(t)}{a}, ~\tau = at, ~\gamma = \frac{c}{a}.</math>  
 
<math>u(\tau) = \frac{d}{c}N(t), ~v(\tau) = \frac{bP(t)}{a}, ~\tau = at, ~\gamma = \frac{c}{a}.</math>  
Система \eqref{syst3} имеет две неподвижные точки: (0,0), (1,1). Стандартный линейный анализ показывает, что точка (0,0) - седло, а для (1,1) матрица Якоби имеет вид:  
+
 
<math>$$J = \begin{bmatrix}
+
== Неподвижные точки системы  ==
0 & -1 \\
+
Найдём [https://sawiki.cs.msu.ru/index.php?title=Матрица_Якоби._Лемма_о_выпрямлении_векторного_поля&action=edit&redlink=1 матрицу Якоби] для нашей системы ОДУ:
\gamma & 0
+
 
\end{bmatrix}. $$</math>
+
\begin{equation*}
 +
J(u, v)= \left(
 +
\begin{array}{cc}
 +
    1-v & -u \\
 +
    \gamma v & \gamma(u-1)
 +
\end{array}
 +
\right)
 +
\end{equation*}
 +
 
 +
Система \eqref{syst3} имеет две [https://sawiki.cs.msu.ru/index.php/Неподвижные_точки_системы неподвижные точки]: $$A(0,0)$$ и $$B(1,1)$$, исследуем их по отдельности:
 +
 
 +
$$\quad\quad$$ 1. В точке $$A(0,0)$$ имеем
 +
\begin{equation*}
 +
J(0,0)= \left(
 +
\begin{array}{cc}
 +
    1 & 0 \\
 +
    0 & -\gamma
 +
\end{array}
 +
\right),
 +
\end{equation*}
 +
Следовательно, $$\lambda_1 = 1$$, $$\lambda_2 = -\gamma$$. Знаем, что $$\gamma = \frac{c}{a} > 0$$, следовательно собственные значения $$\lambda_{1,2}$$ представляют из себя действительные числа противоположных знаков, что означает, что точка $$A(0,0)$$  &mdash; [https://ru.wikipedia.org/wiki/Особая_точка_(дифференциальные_уравнения) седло].
 +
 
 +
$$\quad\quad$$ 2. Для $$B(1,1)$$ матрица Якоби принимает вид:
 +
\begin{equation*}
 +
J(1,1)= \left(
 +
\begin{array}{cc}
 +
    0 & -1 \\
 +
    \gamma & 0
 +
\end{array}
 +
\right),
 +
\end{equation*}
 +
 
 +
[[Файл:ContH.png|мини| Фазовые кривые системы \eqref{syst3} (линии уровня $$H$$) замкнуты.]]
 +
 
 +
а её собственные значения $$\lambda_{1,2} = \pm i \sqrt{\gamma}.$$ Другими словами, положение равновесия $$(1,1)$$ негиперболическое, и линейный анализ не позволяет сделать вывод об его устойчивости. 
 +
 
 +
Фазовые кривые системы \eqref{syst3} являются интегральными кривыми уравнения:
 +
<math>\frac{dv}{dv} = \gamma\frac{v(u-1)}{u(1-v)}, </math>
 +
решение которого:
 +
 
 +
\begin{equation}\label{eq4}
 +
    \gamma u + v - \ln{(u^\gamma v)} = H = H(u,v),
 +
\end{equation}
 +
 
 +
где $$H > H_{min} = 1 + \gamma$$, где $$H_{min}$$ достигается в точке $$u=1, ~v=1$$.
 +
 
 +
Посчитаем частные производные данной функции:
 +
 
 +
\begin{equation}
 +
H_u = \frac{\gamma(u-1)}{u}, ~H_v = \frac{v-1}{v}, ~H_{uu} = \frac{\gamma}{u^2}, ~H_{uv} = 0, ~H_{vv} = \frac{1}{v^2}.
 +
\end{equation}
 +
 
 +
Действительно, обе частные производные первого порядка $$H_u(1,1) = 0, ~H_v(1,1) = 0$$ в ней обнуляются, в то время как производные второго порядка $$H_{uu}(1,1) = A > 0, ~H_{uv} = B = 0, ~H_{vv}(1,1) = C > 0,$$ удовлетворяют неравенству $$AC - B^2 > 0$$, что указывает на достижение минимума в данной точке.
 +
 
 +
'''Определение 1 '''(Производная Ли)
 +
 
 +
Пусть задана динамическая система
 +
\begin{equation}\label{eq5}
 +
    \dot{u} = f(u), \quad u \in \mathbb{R}^n, \quad f:\mathbb{R}^n \rightarrow \mathbb{R}^n
 +
\end{equation}
 +
и дифференцируемая функция $$V:\mathbb{R}^n \rightarrow \mathbb{R}^n$$.
 +
 
 +
Производной Ли или производной вдоль траекторий системы \eqref{eq5} называется выражение
 +
\begin{equation}
 +
    L_{t}V = \langle grad~V, \dot{u} \rangle = \sum \frac{\partial V}{\partial u_{i}}\dot{u_{i}} = \sum \frac{\partial V}{\partial u_{i}}f_{i}(u) = \langle grad~V, f \rangle.
 +
\end{equation}
 +
 
 +
 +
Легко проверить, что $$L_tH(u,v) = 0$$, и, следовательно, функция $$H(u,v)$$ задает [https://ru.wikipedia.org/wiki/Первый_интеграл первый интеграл] системы \eqref{syst3}. Анализируя [https://ru.wikipedia.org/wiki/Изолиния линии уровня] фукции $$H(u,v)$$ можно показать, что для любых $$H > H_{min}$$ они являются замкнутыми кривыми, в общем случае анализ линий уровня функции $$H(u, v)$$ достаточно сложен.
 +
 
 +
== Неустойчивость системы  ==
 +
Если разложить в окрестности точки $$B(1,1)$$ функцию $$H(u, v)$$ в [https://ru.wikipedia.org/wiki/Ряд_Тейлора ряд Тейлора], то получим:
 +
<math> H(u,v) = 1 + \gamma + \gamma\frac{(u-1)^2}{2} + \gamma\frac{(v-1)^2}{2} + \o\left((v-1)^2\right) + \o\left((u-1)^2\right).</math>
 +
 
 +
[[Файл:Solv1.png|мини| Решение системы \eqref{syst3}, начальные условия взяты близко к положению равновесия.]]
 +
[[Файл:Solv2.png|мини| Решение системы \eqref{syst3}, начальные условия взяты далеко от положения равновесия.]]
 +
Последнее равенство означает, что в окрестности точки $$B(1,1)$$ существует обратимая дифференцируемая замена переменных, которая переводит $$H(u,v)$$ в функцию:
 +
<math> G(u,v) = 1 + \gamma + \frac{(x-1)^2}{2} + \frac{(y-1)^2}{2},</math>
 +
линии уровня которой &mdash; окружности. Следовательно, по крайней мере в окрестности точки $$(1,1)$$ линии уровня иисходной функции $$H(u,v)$$ представляют собой замкнутые кривые. Замкнутые фазовые траектории на плоскости $$u,~v$$ отвечают периодическим решениям системы \eqref{syst3}. Начальные условия $$u(0) = u_0, ~v(0) = v_0$$ определяют константу  $$H$$ в выражении \eqref{eq4}. Типичные решения системы \eqref{syst3} приведены на приложенных иллюстрациях.
 +
 
 +
 
 +
Главным же недостатком системы \eqref{syst3} как математической модели экологической системы принято считать её [https://ru.wikipedia.org/wiki/Устойчивость_(динамические_системы) '''структурную неустойчивость''']: малое изменение правых частей \eqref{syst3} в метрике соответствующего пространства функций может приводить к качественному изменению поведения решений.
 +
 
 +
Одой из причин того, что система \eqref{syst3} является структурно неустойчивой, является существование первого интеграла. Системы, обладающие первым интегралом, называются [https://sawiki.cs.msu.ru/index.php/Гамильтоновы_системы '''консервативными'''], и в силу своей структурной неустойчивости редко используются в качестве математических моделей такого рода. Считается, что первые интегралы таких систем описывают некоторые «законы сохранения».
 +
 
 +
Изучение консервативности динамической системы неразделимо со свойством сохранения [https://sawiki.cs.msu.ru/index.php/Фазовый_объём._Теорема_Лиувилля фазового объема].  Любая консервативная система сохраняет фазовый объем (быть может после умножения правых частей на общий скалярный множитель).
 +
Сделав замену переменных $$u = \ln{N}, ~v = \ln{P}$$ в системе \eqref{syst2}, получим:
 +
\begin{equation}\label{syst5}
 +
\begin{cases}
 +
  \dot u = a - be^v,\\
 +
  \dot v = -d + ce^u.\\
 +
\end{cases}
 +
\end{equation}
 +
 
 +
Можно проверить, что для правых частей системы выполняется условие $$\div{f} = 0$$, т.е. система Лотки-Вольтерры сохраняет фазовый объем в $$\mathbb{R}^2_+$$.
 +
 
 +
== Принцип Вольтерры ==
 +
 
 +
'''Утверждение 1 '''(Принцип Вольтерры)
 +
 
 +
Если в системе хищник &mdash; жертва, описываемой моделью \eqref{syst2}, оба вида истребляются равномерно и пропорционально числу их индивидуумов, то среднее число жертв возрастает, а среднее число хищников убывает.
 +
 
 +
''Доказательство.''
 +
 
 +
Пусть $$N(t), ~P(t)$$ — периодические решения системы \eqref{syst2} с периодом $$T$$. Из \eqref{syst2} следует, что:
 +
<math> \frac{d}{dt}\ln{N} = a - bP, \quad \frac{d}{dt}\ln{P} = -d + cN.</math>
 +
 
 +
Интегрируя последние равенства по $$t$$ в промежутке от $$0$$ до $$T$$, получим:
 +
<math> \frac{1}{T}\int_{0}^{T} P(t)dt= \frac{a}{b}, \quad \frac{1}{T}\int_{0}^{T} N(t)dt= \frac{d}{c},</math>
 +
 
 +
так как $$N(T) = N(0), ~P(T) = P(0)$$. То есть среднее число жертв и хищников остается постоянным и равным координатам нетривиального положения равновесия \eqref{syst2}.
 +
 
 +
Если оба вида истребляются равномерно и пропорционально числу их индивидуумов, то у жертв уменьшается коэффициент рождаемости $$a$$, который становится равным $$a - \delta_2$$, а у хищников увеличивается коэффциент смертности $$d$$, который становится равен $$d+\delta_1$$, другими словами, среднее число жертв равно $$\frac{d+\delta_1}{c}$$, а среднее число хищников &mdash; $$\frac{a-\delta_2}{b}$$.<math>\blacksquare</math>  
 +
 
 +
Описанный выше эффект наблюдается в природе. Например, во время первой мировой войны лов рыбы в Адриатическом море был сильно сокращен, что, к удивлению биологов, привело к увеличению числа хищников и уменьшению числа жертв (указанный факт стал одним из побудительных мотивов записать и проанализировать модель \eqref{syst2}).
 +
 
 +
Кроме всего прочего, принцип Вольтерры показывает двойственный характер применения средств от насекомых (инсектицидов) для сохранения урожая на полях. Почти все такие химические вещества действуют не только на вредителей, но и на их естественных врагов, что зачастую приводит к увеличению числа вредителей и уменьшению, например, числа птиц, питающихся этими вредителями. Отметим также, что принцип Вольтерры впервые теоретически показал, что в экосистеме "хищник &mdash; жертва" популяция жертв более чувствительна к процессу пропорционального уменьшения особей в популяции.
 +
 
 +
==Список литературы==
 +
 
 +
1. Братусь А.С., Новожилов А.С., Платонов А.П. Динамические системы и модели биологии 2011.
 +
 
 +
2. Абрамова В.В. Лекции по курсу "Динамические системы и биоматематика", 2023.

Текущая версия на 17:43, 19 декабря 2023

Система Лотки-Вольтерры "хищник — жертва"

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

Модель Заяц — Волк попадалась в художественных произведениях уже полвека тому назад

Если $$N(t)$$ — численность жертв, $$P(t)$$ — численность хищников в момент времени $$t$$, то модель Лотки-Вольтерры имеет вид:

\begin{equation}\label{syst2} \begin{cases} \dot N = aN - bNP,\\ \dot P = -dP + cNP,\\ \end{cases} \end{equation}

где $$a, ~b, ~c, ~d$$ — положительные постоянные.

В основу системы \eqref{syst2} положены предположения, характеризующиеся следующими гипотезами:

  1. В отсутствиии хищников жертвы размножаются неограниченно ($$\dot N = aN$$).
  2. Хищники в отсутствии жертв вымирают ($$\dot P = -dP$$).
  3. Слагаемые, пропорциональные $$NP$$, рассматриваются как превращение энергии одного источника в энергию другого (эффект влияния популяции хищников на популяцию жертв заключается в уменьшении относительной скорости прироста численности жертв на величину, пропорциональную численности хищников).

Так как мы смотрим на систему \eqref{syst2} как на математическую модель взаимодействущих популяций, естественно будет считать фазовым пространством множество $$\mathbb{R}^2_+$$, которое является инвариантным, так как любая траектория, начинающаяся в $$\mathbb{R}^2_+$$ не может пересечь линии $$N = 0$$ и $$P = 0$$, являющиеся фазовыми кривыми.

В безразмерных переменных система принимает вид:

\begin{equation}\label{syst3} \begin{cases} \dot u = u(1-v),\\ \dot v = \gamma v(u-1),\\ \end{cases} \end{equation}

где\[u(\tau) = \frac{d}{c}N(t), ~v(\tau) = \frac{bP(t)}{a}, ~\tau = at, ~\gamma = \frac{c}{a}.\]

Неподвижные точки системы

Найдём матрицу Якоби для нашей системы ОДУ:

\begin{equation*} J(u, v)= \left( \begin{array}{cc} 1-v & -u \\ \gamma v & \gamma(u-1) \end{array} \right) \end{equation*}

Система \eqref{syst3} имеет две неподвижные точки: $$A(0,0)$$ и $$B(1,1)$$, исследуем их по отдельности:

$$\quad\quad$$ 1. В точке $$A(0,0)$$ имеем \begin{equation*} J(0,0)= \left( \begin{array}{cc} 1 & 0 \\ 0 & -\gamma \end{array} \right), \end{equation*} Следовательно, $$\lambda_1 = 1$$, $$\lambda_2 = -\gamma$$. Знаем, что $$\gamma = \frac{c}{a} > 0$$, следовательно собственные значения $$\lambda_{1,2}$$ представляют из себя действительные числа противоположных знаков, что означает, что точка $$A(0,0)$$  — седло.

$$\quad\quad$$ 2. Для $$B(1,1)$$ матрица Якоби принимает вид: \begin{equation*} J(1,1)= \left( \begin{array}{cc} 0 & -1 \\ \gamma & 0 \end{array} \right), \end{equation*}

Фазовые кривые системы \eqref{syst3} (линии уровня $$H$$) замкнуты.

а её собственные значения $$\lambda_{1,2} = \pm i \sqrt{\gamma}.$$ Другими словами, положение равновесия $$(1,1)$$ негиперболическое, и линейный анализ не позволяет сделать вывод об его устойчивости.

Фазовые кривые системы \eqref{syst3} являются интегральными кривыми уравнения\[\frac{dv}{dv} = \gamma\frac{v(u-1)}{u(1-v)}, \] решение которого:

\begin{equation}\label{eq4} \gamma u + v - \ln{(u^\gamma v)} = H = H(u,v), \end{equation}

где $$H > H_{min} = 1 + \gamma$$, где $$H_{min}$$ достигается в точке $$u=1, ~v=1$$.

Посчитаем частные производные данной функции:

\begin{equation} H_u = \frac{\gamma(u-1)}{u}, ~H_v = \frac{v-1}{v}, ~H_{uu} = \frac{\gamma}{u^2}, ~H_{uv} = 0, ~H_{vv} = \frac{1}{v^2}. \end{equation}

Действительно, обе частные производные первого порядка $$H_u(1,1) = 0, ~H_v(1,1) = 0$$ в ней обнуляются, в то время как производные второго порядка $$H_{uu}(1,1) = A > 0, ~H_{uv} = B = 0, ~H_{vv}(1,1) = C > 0,$$ удовлетворяют неравенству $$AC - B^2 > 0$$, что указывает на достижение минимума в данной точке.

Определение 1 (Производная Ли)

Пусть задана динамическая система \begin{equation}\label{eq5} \dot{u} = f(u), \quad u \in \mathbb{R}^n, \quad f:\mathbb{R}^n \rightarrow \mathbb{R}^n \end{equation} и дифференцируемая функция $$V:\mathbb{R}^n \rightarrow \mathbb{R}^n$$.

Производной Ли или производной вдоль траекторий системы \eqref{eq5} называется выражение \begin{equation} L_{t}V = \langle grad~V, \dot{u} \rangle = \sum \frac{\partial V}{\partial u_{i}}\dot{u_{i}} = \sum \frac{\partial V}{\partial u_{i}}f_{i}(u) = \langle grad~V, f \rangle. \end{equation}


Легко проверить, что $$L_tH(u,v) = 0$$, и, следовательно, функция $$H(u,v)$$ задает первый интеграл системы \eqref{syst3}. Анализируя линии уровня фукции $$H(u,v)$$ можно показать, что для любых $$H > H_{min}$$ они являются замкнутыми кривыми, в общем случае анализ линий уровня функции $$H(u, v)$$ достаточно сложен.

Неустойчивость системы

Если разложить в окрестности точки $$B(1,1)$$ функцию $$H(u, v)$$ в ряд Тейлора, то получим\[ H(u,v) = 1 + \gamma + \gamma\frac{(u-1)^2}{2} + \gamma\frac{(v-1)^2}{2} + \o\left((v-1)^2\right) + \o\left((u-1)^2\right).\]

Решение системы \eqref{syst3}, начальные условия взяты близко к положению равновесия.
Решение системы \eqref{syst3}, начальные условия взяты далеко от положения равновесия.

Последнее равенство означает, что в окрестности точки $$B(1,1)$$ существует обратимая дифференцируемая замена переменных, которая переводит $$H(u,v)$$ в функцию\[ G(u,v) = 1 + \gamma + \frac{(x-1)^2}{2} + \frac{(y-1)^2}{2},\] линии уровня которой — окружности. Следовательно, по крайней мере в окрестности точки $$(1,1)$$ линии уровня иисходной функции $$H(u,v)$$ представляют собой замкнутые кривые. Замкнутые фазовые траектории на плоскости $$u,~v$$ отвечают периодическим решениям системы \eqref{syst3}. Начальные условия $$u(0) = u_0, ~v(0) = v_0$$ определяют константу  $$H$$ в выражении \eqref{eq4}. Типичные решения системы \eqref{syst3} приведены на приложенных иллюстрациях.


Главным же недостатком системы \eqref{syst3} как математической модели экологической системы принято считать её структурную неустойчивость: малое изменение правых частей \eqref{syst3} в метрике соответствующего пространства функций может приводить к качественному изменению поведения решений.

Одой из причин того, что система \eqref{syst3} является структурно неустойчивой, является существование первого интеграла. Системы, обладающие первым интегралом, называются консервативными, и в силу своей структурной неустойчивости редко используются в качестве математических моделей такого рода. Считается, что первые интегралы таких систем описывают некоторые «законы сохранения».

Изучение консервативности динамической системы неразделимо со свойством сохранения фазового объема. Любая консервативная система сохраняет фазовый объем (быть может после умножения правых частей на общий скалярный множитель). Сделав замену переменных $$u = \ln{N}, ~v = \ln{P}$$ в системе \eqref{syst2}, получим: \begin{equation}\label{syst5} \begin{cases} \dot u = a - be^v,\\ \dot v = -d + ce^u.\\ \end{cases} \end{equation}

Можно проверить, что для правых частей системы выполняется условие $$\div{f} = 0$$, т.е. система Лотки-Вольтерры сохраняет фазовый объем в $$\mathbb{R}^2_+$$.

Принцип Вольтерры

Утверждение 1 (Принцип Вольтерры)

Если в системе хищник — жертва, описываемой моделью \eqref{syst2}, оба вида истребляются равномерно и пропорционально числу их индивидуумов, то среднее число жертв возрастает, а среднее число хищников убывает.

Доказательство.

Пусть $$N(t), ~P(t)$$ — периодические решения системы \eqref{syst2} с периодом $$T$$. Из \eqref{syst2} следует, что\[ \frac{d}{dt}\ln{N} = a - bP, \quad \frac{d}{dt}\ln{P} = -d + cN.\]

Интегрируя последние равенства по $$t$$ в промежутке от $$0$$ до $$T$$, получим\[ \frac{1}{T}\int_{0}^{T} P(t)dt= \frac{a}{b}, \quad \frac{1}{T}\int_{0}^{T} N(t)dt= \frac{d}{c},\]

так как $$N(T) = N(0), ~P(T) = P(0)$$. То есть среднее число жертв и хищников остается постоянным и равным координатам нетривиального положения равновесия \eqref{syst2}.

Если оба вида истребляются равномерно и пропорционально числу их индивидуумов, то у жертв уменьшается коэффициент рождаемости $$a$$, который становится равным $$a - \delta_2$$, а у хищников увеличивается коэффциент смертности $$d$$, который становится равен $$d+\delta_1$$, другими словами, среднее число жертв равно $$\frac{d+\delta_1}{c}$$, а среднее число хищников — $$\frac{a-\delta_2}{b}$$.\(\blacksquare\)

Описанный выше эффект наблюдается в природе. Например, во время первой мировой войны лов рыбы в Адриатическом море был сильно сокращен, что, к удивлению биологов, привело к увеличению числа хищников и уменьшению числа жертв (указанный факт стал одним из побудительных мотивов записать и проанализировать модель \eqref{syst2}).

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

Список литературы

1. Братусь А.С., Новожилов А.С., Платонов А.П. Динамические системы и модели биологии 2011.

2. Абрамова В.В. Лекции по курсу "Динамические системы и биоматематика", 2023.