Логистическое уравнение и его свойства

Определения

Пусть $$ N(t) $$ - численность изолированной популяции в момент времени $$ t $$. Скорость её изменения может быть представлена в следующем виде

\[ \dot N = \text{рождаемость} - \text{смертность} + \text{миграция} . \]

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

\[ \dot N = bN - cN, \ \ \ \ \ \ \ \Rightarrow \ \ \ \ \ \ \ \dot N = aN, \ a = b - c .\]

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

\[ \dot N = N f(N) , \]

где $$ f(N) $$ - коэффициент скорости роста популяции. Разложим $$ F(N) $$ в ряд Тейлора в окрестности нуля и отбросим все члены, кроме первых двух. Получим

\[ \dot N = N(a + bN), \]

где $$ a,b $$ - некоторые постоянные, причем естественно предположить, что $$ a > 0, \ b < 0 $$. Именно таким образом Альфред Лотка (Alfred Lotka, 1880–1949, один из создателей математической биологии) пришел к уравнению, которое стало известно как логистическое уравнение, которое после некоторых переобозначений запишется в виде:

\[ \dot N = rN \left( 1 - \dfrac{N}{k} \right) . \]

Здесь $$ r $$ и $$K$$ - положительные параметры.

Свойства

Можно заметить, что когда $$ N(t) $$ мало, то $$ \dot N \approx rN $$. Значит, $$N$$ экспоненциально растёт. Параметр $$ K $$ интерпретируется как потенциальная емкость экологической системы, которая определяется доступным наличным количеством ресурсов. Величина $$ K $$ определяет предельное значение численности популяции. Для доказательства этого факта решим логистическое уравнение:

\[ \dfrac{dP}{dt} = rP \left( 1 - \dfrac{P}{K} \right), \]

\[ -\dfrac{1}{r} \int \dfrac{dP}{\frac{P^2}{K} - P} = \int dt , \]

\[ -\dfrac{1}{r} \int \dfrac{dP}{ \left( \frac{P}{\sqrt{K}} - \frac{1}{2} \sqrt{k} \right)^2 - \frac{1}{4}k} = \int dt . \]

Сделаем замену в правой части равенства: $$ z = \frac{P}{\sqrt{K}} - \frac{1}{2} \sqrt{k}, \ dz = \frac{dP}{\sqrt{k}} $$. Получим

\[ - \dfrac{\sqrt{K}}{r} \int \dfrac{dz}{z^2 - \frac{1}{4} k} = \int dt , \]

\[ -\dfrac{1}{r} \ln{ \left| \dfrac{z - \frac{1}{2} \sqrt{k} }{z + \frac{1}{2} \sqrt{k}} \right| } = t + C , \]

\[ \dfrac{z - \frac{1}{2} \sqrt{k} }{z + \frac{1}{2} \sqrt{k}} = Ce^{-rt} . \]

Возвращаясь к P, получим

\[ 1 - \dfrac{K}{P} = C e^{-rt} ,\] \[ P = \dfrac{K}{1 - Ce^{-rt}}. \]

Найдем $$ C $$ из начального условия $$ P(0) = P_0: \ C = 1 - \frac{K}{P_0}. $$ Окончательно получим

\[ P = \dfrac{K}{1 - \left( 1 - \frac{K}{P_0} \right) e^{-rt} } = \dfrac{K e^{-rt} }{e^{rt} - 1 + \frac{K}{P_0}} = \dfrac{P_0 K e^{rt}}{P_0 (e^{rt} - 1) + K} \]