Логистическая модель роста популяции Ферхюльста

Логистическая модель роста популяции Ферхюльста описывает, как численность живых организмов сначала растёт почти экспоненциально, а затем замедляется и выходит на плато, упираясь в ограниченные ресурсы среды. Это базовая модель математической экологии: она исправляет нереалистичный неограниченный рост модели Мальтуса, добавляя «тормоз» - внутривидовую конкуренцию за пищу, пространство и свет. Ниже разберём дифференциальное уравнение Ферхюльста, его точное решение, смысл ёмкости среды и точки перегиба, а интерактивный калькулятор покажет, как параметры лепят S-образную кривую.

Дифференциальное уравнение Ферхюльста

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

$\frac{dN}{dt} = r\,N\left(1 - \frac{N}{K}\right),$

где $N(t)$ - численность популяции в момент времени $t$, $r$ - удельная (мальтузианская) скорость роста при малой плотности, а $K$ - ёмкость среды, то есть максимальная численность, которую территория способна прокормить. Множитель $\left(1 - \tfrac{N}{K}\right)$ - это и есть ключевая идея Ферхюльста: пока $N \ll K$, он близок к единице и рост почти экспоненциальный; когда $N \to K$, он стремится к нулю и рост останавливается.

Рассчитать для своей задачи

Связь с моделью Мальтуса

Если убрать ограничивающий множитель, уравнение превращается в $\tfrac{dN}{dt} = rN$ - это модель Мальтуса экспоненциального роста, дающая неограниченное удвоение численности. Модель Ферхюльста - это её естественное продолжение: при $N \ll K$ обе кривые практически совпадают, и логистический рост неотличим от экспоненциального. Расхождение начинается, когда популяция заполняет заметную долю ёмкости: мальтузианская кривая уходит в бесконечность, а логистическая сворачивает к горизонтали $N = K$. Поэтому модель Мальтуса корректна лишь на ранней стадии - как начальный участок S-образной кривой Ферхюльста.

Решение уравнения и S-образная кривая

Логистическое уравнение разделяется, и его точное решение - логистическая функция:

$N(t) = \frac{K}{1 + \left(\dfrac{K - N_0}{N_0}\right)e^{-rt}},$

где $N_0 = N(0)$ - начальная численность. График $N(t)$ имеет характерную S-образную (сигмоидную) форму: медленный старт при малой численности, бурный рост в середине и плавный выход на плато $K$ при $t \to \infty$. Эта кривая - визитная карточка модели: по ней узнают логистический рост в данных по бактериям, дрожжам, оленям на острове и даже в распространении технологий.

Рассчитать для своей задачи

Ёмкость среды и устойчивое равновесие

Ёмкость среды $K$ - это устойчивая стационарная точка уравнения. Производная $\tfrac{dN}{dt}$ обращается в ноль при $N = 0$ и при $N = K$. Точка $N = 0$ неустойчива: малейшее отклонение вверх запускает рост. Точка $N = K$ устойчива: если численность случайно превысит $K$, множитель $\left(1 - \tfrac{N}{K}\right)$ станет отрицательным и популяция начнёт убывать обратно к $K$; если упадёт ниже - снова дорастёт до $K$. Поэтому ёмкость среды работает как аттрактор: куда бы ни стартовала популяция (выше или ниже $K$), она сходится к этому уровню.

Точка перегиба и максимальная скорость роста

Скорость прироста $\tfrac{dN}{dt}$ как функция от $N$ - это парабола с максимумом ровно посередине, при $N = K/2$. Именно в этот момент популяция растёт быстрее всего: ресурсов ещё хватает, а особей уже много. На S-образной кривой это точка перегиба - здесь рост из ускоряющегося превращается в замедляющийся. Максимальная скорость равна

$\left.\frac{dN}{dt}\right|_{\max} = \frac{rK}{4}.$

Этот результат важен в практике: в рыболовстве и лесном хозяйстве уровень $N = K/2$ соответствует максимальному устойчивому вылову (MSY) - именно при половине ёмкости среды популяция восстанавливается быстрее всего.

Биологический смысл параметров

Параметр $r$ задаёт темп: чем он больше, тем круче S и тем быстрее популяция доходит до плато. Биологически $r$ зависит от плодовитости и смертности вида при изобилии ресурсов. Параметр $K$ задаёт уровень плато и определяется ёмкостью местообитания: площадью, кормовой базой, доступностью укрытий. Начальная численность $N_0$ сдвигает кривую по времени, но не меняет её форму и финальный уровень. Эта триада ($r$, $K$, $N_0$) полностью определяет траекторию - тесно связанная с динамикой видов тема разобрана в материале про уравнение Лотки-Вольтерры хищник-жертва, где взаимодействие двух популяций надстраивается над логистической базой.

Рассчитать для своей задачи

Где модель работает, а где нет

Логистическая модель отлично описывает рост в замкнутой среде с постоянными ресурсами: культуры микроорганизмов в пробирке, дрожжи в питательном растворе, изолированные популяции на островах. Однако реальные популяции часто отклоняются от идеальной S-кривой: при перенаселении возможен «проскок» через ёмкость с последующим обвалом, сезонные колебания ресурсов меняют $K$ во времени, а запаздывание реакции на плотность порождает колебания вокруг $K$. Эти эффекты моделируют усложнениями - логистическим уравнением с запаздыванием или дискретным логистическим отображением, где при больших $r$ возникает даже хаос.

Частые ошибки

• Путают $r$ и фактическую скорость роста. $r$ - это удельная скорость при малой плотности; реальный прирост $\tfrac{dN}{dt}$ зависит ещё и от того, насколько $N$ близко к $K$.
• Считают, что популяция растёт быстрее всего в начале. Максимум скорости приходится на $N = K/2$, а не на старт: при малом $N$ особей слишком мало, при $N$ близком к $K$ мешает конкуренция.
• Забывают про устойчивость $K$. Если в задаче $N_0 > K$, популяция не растёт, а убывает к ёмкости среды - кривая идёт сверху вниз.
• Подставляют $N_0 = 0$. При нулевой начальной численности популяция остаётся нулём: точка $N = 0$ - стационарная.
• Смешивают непрерывную и дискретную модель. Сигмоида Ферхюльста - это непрерывное уравнение; дискретное логистическое отображение ведёт себя иначе и может давать колебания и хаос.

FAQ

Чем модель Ферхюльста отличается от модели Мальтуса? Мальтус даёт неограниченный экспоненциальный рост $N = N_0 e^{rt}$. Ферхюльст добавляет множитель $\left(1 - \tfrac{N}{K}\right)$, который тормозит рост по мере заполнения среды, и кривая выходит на плато $K$.

Что такое ёмкость среды $K$? Это максимальная численность популяции, которую ресурсы территории способны поддерживать неограниченно долго. В модели $K$ - устойчивое равновесие, к которому популяция сходится сверху или снизу.

Где находится точка перегиба логистической кривой? При $N = K/2$. В этой точке скорость роста максимальна и равна $rK/4$, а S-образная кривая переходит от ускоряющегося роста к замедляющемуся.

Коротко

Логистическая модель Ферхюльста описывает рост популяции уравнением $\tfrac{dN}{dt} = rN\left(1 - \tfrac{N}{K}\right)$, решение которого - S-образная логистическая функция с плато на уровне ёмкости среды $K$. Параметр $r$ задаёт темп, $K$ - потолок, а точка перегиба при $N = K/2$ отмечает момент максимального прироста $rK/4$. Модель обобщает экспоненциальный рост Мальтуса, добавляя внутривидовую конкуренцию, и служит фундаментом математической экологии.