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

Материал из sawiki
Перейти к навигации Перейти к поиску
 
(не показано 8 промежуточных версий этого же участника)
Строка 1: Строка 1:
== Система Лотки-Вольтерры "хищник-жертва" ==
+
== Система Лотки-Вольтерры "хищник — жертва" ==
Одной из первых рассмотренных математических моделей взаимодействующих популяций является система ОДУ, именуемая '''моделью [https://ru.wikipedia.org/wiki/Лотка,_Альфред_Джеймс Лотки]-[https://ru.wikipedia.org/wiki/Вольтерра,_Вито Вольтерры]''', которая описывает взаимодействие двух видов, один из которых является хищником, а другой — жертвой (например, экологическая система волки-зайцы)
+
Одной из первых рассмотренных математических моделей взаимодействующих популяций является система ОДУ, именуемая '''моделью [https://ru.wikipedia.org/wiki/Лотка,_Альфред_Джеймс Лотки]-[https://ru.wikipedia.org/wiki/Вольтерра,_Вито Вольтерры]''', которая описывает взаимодействие двух видов, один из которых является хищником, а другой — жертвой (например, экологическая система волки-зайцы).
  
[[Файл:Заяц-волк.png|мини| Модель Заяц-Волк попадалась в художественных произведениях уже полвека тому назад]]
+
[[Файл:Заяц-волк.png|мини| Модель Заяц — Волк попадалась в художественных произведениях уже полвека тому назад]]
  
 
Если $$N(t)$$ — численность жертв, $$P(t)$$ — численность хищников в момент времени $$t$$, то модель Лотки-Вольтерры имеет вид:  
 
Если $$N(t)$$ — численность жертв, $$P(t)$$ — численность хищников в момент времени $$t$$, то модель Лотки-Вольтерры имеет вид:  
Строка 13: Строка 13:
 
\end{equation}
 
\end{equation}
  
где a, b, c, d — положительные постоянные.
+
где $$a, ~b, ~c, ~d$$ — положительные постоянные.
  
 
В основу системы \eqref{syst2} положены предположения, характеризующиеся следующими гипотезами:
 
В основу системы \eqref{syst2} положены предположения, характеризующиеся следующими гипотезами:
#в отсутствиии хищников жертвы размножаются неограниченно ($$\dot N = aN$$);
+
#В отсутствиии хищников жертвы размножаются неограниченно ($$\dot N = aN$$).
#хищники в отсутствии жертв вымирают ($$\dot P = -dP$$);
+
#Хищники в отсутствии жертв вымирают ($$\dot P = -dP$$).
#слагаемые, пропорциональные NP, рассматриваются как превращение энергии одного источника в энергию другого (эффект влияния популяции хищников на популяцию жертв заключается в уменьшении относительной скорости прироста численности жертв на величину, пропорциональную численности хищников).
+
#Слагаемые, пропорциональные $$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 фазовыми кривыми].  
 
Так как мы смотрим на систему \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 фазовыми кривыми].  
Строка 35: Строка 35:
  
 
== Неподвижные точки системы  ==
 
== Неподвижные точки системы  ==
Система \eqref{syst3} имеет две [https://sawiki.cs.msu.ru/index.php/Неподвижные_точки_системы неподвижные точки]: $$A(0,0)$$ и $$B(1,1)$$.  Стандартный линейный анализ показывает, что точка $$A(0,0)$$ - [https://ru.wikipedia.org/wiki/Седловая_точка седло], а для $$B(1,1)$$ [https://sawiki.cs.msu.ru/index.php?title=Матрица_Якоби._Лемма_о_выпрямлении_векторного_поля&action=edit&redlink=1 матрица Якоби] имеет вид:  
+
Найдём [https://sawiki.cs.msu.ru/index.php?title=Матрица_Якоби._Лемма_о_выпрямлении_векторного_поля&action=edit&redlink=1 матрицу Якоби] для нашей системы ОДУ:
  
 
\begin{equation*}
 
\begin{equation*}
J\vert_{(1,1)}= \left(
+
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)$$  — [https://ru.wikipedia.org/wiki/Особая_точка_(дифференциальные_уравнения) седло].
 +
 
 +
$$\quad\quad$$ 2. Для $$B(1,1)$$ матрица Якоби принимает вид:
 +
\begin{equation*}
 +
J(1,1)= \left(
 
\begin{array}{cc}
 
\begin{array}{cc}
 
     0 & -1 \\
 
     0 & -1 \\
 
     \gamma & 0
 
     \gamma & 0
 
\end{array}
 
\end{array}
\right)
+
\right),
 
\end{equation*}
 
\end{equation*}
  
[[Файл:Contour h.png|мини| Фазовые кривые системы \eqref{syst3} (линии уровня $$H$$) замкнуты.]]
+
[[Файл:ContH.png|мини| Фазовые кривые системы \eqref{syst3} (линии уровня $$H$$) замкнуты.]]
  
а её собственные значения $$\lambda_{1,2} = \pm i \sqrt{\gamma}.$$ Другими словами, положение равновесия (1,1) негиперболическое, и линейный анализ не позволяет сделать вывод об его устойчивости.   
+
а её собственные значения $$\lambda_{1,2} = \pm i \sqrt{\gamma}.$$ Другими словами, положение равновесия $$(1,1)$$ негиперболическое, и линейный анализ не позволяет сделать вывод об его устойчивости.   
  
 
Фазовые кривые системы \eqref{syst3} являются интегральными кривыми уравнения:
 
Фазовые кривые системы \eqref{syst3} являются интегральными кривыми уравнения:
Строка 55: Строка 78:
  
 
\begin{equation}\label{eq4}
 
\begin{equation}\label{eq4}
     \gamma u + v - \ln{u^\gamma v} = H = H(u,v),
+
     \gamma u + v - \ln{(u^\gamma v)} = H = H(u,v),
 
\end{equation}
 
\end{equation}
  
 
где $$H > H_{min} = 1 + \gamma$$, где $$H_{min}$$ достигается в точке $$u=1, ~v=1$$.  
 
где $$H > H_{min} = 1 + \gamma$$, где $$H_{min}$$ достигается в точке $$u=1, ~v=1$$.  
  
Действительно,  
+
Посчитаем частные производные данной функции:
\[H_u(1,1) = 0, ~H_v(1,1) = 0, ~\text{и} ~H_{uu}(1,1) < 0, ~H_{uv}(1,1) = 0, ~H_{vv}(1,1) < 0.\]
+
 
Легко проверить, что $$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}$$ они являются замкнутыми кривыми (см. '''рис.2'''), в общем случае анализ линий уровня функции $$H(u, v)$$ достаточно сложен.  
+
\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)$$ достаточно сложен.  
  
 
== Неустойчивость системы  ==
 
== Неустойчивость системы  ==
Строка 68: Строка 111:
 
<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>
 
<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>
  
[[Файл:Solution 1.png|мини| Решение системы \eqref{syst3}, начальные условия взяты близко к положению равновесия.]]
+
[[Файл:Solv1.png|мини| Решение системы \eqref{syst3}, начальные условия взяты близко к положению равновесия.]]
[[Файл:Solution 2.png|мини| Решение системы \eqref{syst3}, начальные условия взяты далеко от положения равновесия.]]
+
[[Файл:Solv2.png|мини| Решение системы \eqref{syst3}, начальные условия взяты далеко от положения равновесия.]]
 
Последнее равенство означает, что в окрестности точки $$B(1,1)$$ существует обратимая дифференцируемая замена переменных, которая переводит $$H(u,v)$$ в функцию:
 
Последнее равенство означает, что в окрестности точки $$B(1,1)$$ существует обратимая дифференцируемая замена переменных, которая переводит $$H(u,v)$$ в функцию:
 
<math> G(u,v) = 1 + \gamma + \frac{(x-1)^2}{2} + \frac{(y-1)^2}{2},</math>
 
<math> G(u,v) = 1 + \gamma + \frac{(x-1)^2}{2} + \frac{(y-1)^2}{2},</math>
линии уровня которой - окружности. Следовательно, по крайней мере в окрестности точки $$(1,1)$$ линии уровня иисходной функции $$H(u,v)$$ представляют собой замкнутые кривые. Замкнутые фазовые траектории на плоскости $$u,~v$$ отвечают периодическим решениям системы \eqref{syst3}. Начальные условия $$u(0) = u_0, ~v(0) = v_0$$ определяют константу H в выражении \eqref{eq4}. Типичные решения системы \eqref{syst3} приведены на '''рис.3''' и '''рис.4'''.
+
линии уровня которой &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://ru.wikipedia.org/wiki/Устойчивость_(динамические_системы) '''структурную неустойчивость''']: малое изменение правых частей \eqref{syst3} в метрике соответствующего пространства функций может приводить к качественному изменению поведения решений.  
  
Одой из причин того, что система \eqref{syst3} является структурно неустойчивой, является существование первого интеграла. Системы, обладающие первым интегралом, называются [https://ru.wikipedia.org/wiki/Консервативная_система#:~:text=Консервативная%20система%20(от%20лат.,системы%20верен%20закон%20сохранения%20энергии. '''консервативными'''], и в силу своей структурной неустойчивости редко используются в качестве математических моделей такого рода. Считается, что первые интегралы таких систем описывают некоторые «законы сохранения».
+
Одой из причин того, что система \eqref{syst3} является структурно неустойчивой, является существование первого интеграла. Системы, обладающие первым интегралом, называются [https://sawiki.cs.msu.ru/index.php/Гамильтоновы_системы '''консервативными'''], и в силу своей структурной неустойчивости редко используются в качестве математических моделей такого рода. Считается, что первые интегралы таких систем описывают некоторые «законы сохранения».
  
Изучение консервативности динамической системы неразделимо со свойством сохранения фазового объема.  [https://ru.wikipedia.org/wiki/Теорема_Лиувилля_о_сохранении_фазового_объёма Любая консервативная система сохраняет фазовый объем] (быть может после умножения правых частей на общий скалярный множитель).
+
Изучение консервативности динамической системы неразделимо со свойством сохранения [https://sawiki.cs.msu.ru/index.php/Фазовый_объём._Теорема_Лиувилля фазового объема].  Любая консервативная система сохраняет фазовый объем (быть может после умножения правых частей на общий скалярный множитель).
 
Сделав замену переменных $$u = \ln{N}, ~v = \ln{P}$$ в системе \eqref{syst2}, получим:
 
Сделав замену переменных $$u = \ln{N}, ~v = \ln{P}$$ в системе \eqref{syst2}, получим:
 
\begin{equation}\label{syst5}
 
\begin{equation}\label{syst5}
Строка 92: Строка 135:
 
== Принцип Вольтерры ==
 
== Принцип Вольтерры ==
  
'''Утверждение 1 '''(Принцип Вольтерры) Если в системе хищник–жертва, описываемой моделью \eqref{syst2}, оба вида истребляются равномерно и пропорционально числу их индивидуумов, то среднее число жертв возрастает, а среднее число хищников убывает.
+
'''Утверждение 1 '''(Принцип Вольтерры)  
 +
 
 +
Если в системе хищник &mdash; жертва, описываемой моделью \eqref{syst2}, оба вида истребляются равномерно и пропорционально числу их индивидуумов, то среднее число жертв возрастает, а среднее число хищников убывает.
 +
 
 +
''Доказательство.''
  
''Доказательство.'' Пусть $$N(t), ~P(t)$$ — периодические решения системы \eqref{syst2} с периодом $$T$$. Из \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>
 
<math> \frac{d}{dt}\ln{N} = a - bP, \quad \frac{d}{dt}\ln{P} = -d + cN.</math>
  
Строка 106: Строка 153:
 
Описанный выше эффект наблюдается в природе. Например, во время первой мировой войны лов рыбы в Адриатическом море был сильно сокращен, что, к удивлению биологов, привело к увеличению числа хищников и уменьшению числа жертв (указанный факт стал одним из побудительных мотивов записать и проанализировать модель \eqref{syst2}).  
 
Описанный выше эффект наблюдается в природе. Например, во время первой мировой войны лов рыбы в Адриатическом море был сильно сокращен, что, к удивлению биологов, привело к увеличению числа хищников и уменьшению числа жертв (указанный факт стал одним из побудительных мотивов записать и проанализировать модель \eqref{syst2}).  
  
Кроме всего прочего, принцип Вольтерры показывает двойственный характер применения средств от насекомых (инсектицидов) для сохранения урожая на полях. Почти все такие химические вещества действуют не только на вредителей, но и на их естественных врагов, что зачастую приводит к увеличению числа вредителей и уменьшению, например, числа птиц, питающихся этими вредителями. Отметим также, что принцип Вольтерры впервые теоретически показал, что в экосистеме "хищник–жертва" популяция жертв более чувствительна к процессу пропорционального уменьшения особей в популяции.
+
Кроме всего прочего, принцип Вольтерры показывает двойственный характер применения средств от насекомых (инсектицидов) для сохранения урожая на полях. Почти все такие химические вещества действуют не только на вредителей, но и на их естественных врагов, что зачастую приводит к увеличению числа вредителей и уменьшению, например, числа птиц, питающихся этими вредителями. Отметим также, что принцип Вольтерры впервые теоретически показал, что в экосистеме "хищник &mdash; жертва" популяция жертв более чувствительна к процессу пропорционального уменьшения особей в популяции.
  
 
==Список литературы==
 
==Список литературы==

Текущая версия на 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.