Модель Мальтуса: экспоненциальный рост популяции

Модель Мальтуса экспоненциального роста - простейшее описание того, как меняется численность популяции, если ресурсы безграничны, а каждая особь размножается с постоянной скоростью. Сформулированная ещё в 1798 году Томасом Мальтусом, она лежит в основе всей математической экологии: с неё начинают разбор динамики популяций, прежде чем переходить к более реалистичным моделям. Ниже разберём вывод уравнения, смысл коэффициента $r$, расчёт времени удвоения и границы применимости. Чтобы сразу подставить свои числа, воспользуйтесь калькулятором ниже.

Суть модели Мальтуса

Основное допущение Мальтуса: скорость прироста численности пропорциональна самой численности. Чем больше особей, тем больше потомства появляется за единицу времени. Если обозначить численность популяции через $N$, а время через $t$, это допущение записывается дифференциальным уравнением:

$\frac{dN}{dt} = rN$

Здесь $r$ - удельная (мальтузианская) скорость роста, равная разности коэффициентов рождаемости $b$ и смертности $d$: $r = b - d$. Левая часть - мгновенная скорость изменения численности, правая говорит, что эта скорость линейно зависит от текущего размера популяции. Никаких внешних ограничений в уравнении нет: ни нехватки пищи, ни хищников, ни конкуренции. Именно поэтому модель называют моделью неограниченного, или экспоненциального, роста.

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

Решение уравнения и формула роста

Уравнение $\frac{dN}{dt} = rN$ решается разделением переменных. Перенесём $N$ влево, $dt$ вправо и проинтегрируем:

$\int \frac{dN}{N} = \int r\, dt \quad\Rightarrow\quad \ln N = rt + C$

Потенцируя и подставляя начальное условие $N(0) = N_0$, получаем итоговую формулу модели Мальтуса экспоненциального роста:

$N(t) = N_0 \, e^{rt}$

где $N_0$ - начальная численность, $e$ - основание натурального логарифма, $r$ - удельная скорость роста, $t$ - время. При $r > 0$ популяция растёт по экспоненте, при $r < 0$ - убывает (вымирает), при $r = 0$ численность неизменна. График при $r > 0$ имеет характерную J-образную форму: сначала рост едва заметен, затем кривая уходит почти вертикально вверх.

Биологический смысл коэффициента r

Параметр $r$ - ключевая величина всей модели. Он показывает, на сколько долей увеличивается популяция за единицу времени в расчёте на одну особь. Размерность $r$ - обратное время (например, $\text{год}^{-1}$ или $\text{сут}^{-1}$). У бактерий в идеальной среде $r$ велик, у крупных млекопитающих с долгим периодом созревания - мал. Иногда вместо непрерывного $r$ используют коэффициент конечной скорости роста $\lambda = e^{r}$: тогда за каждый шаг численность умножается на $\lambda$. Связь экологической ниши и реализуемого $r$ подробнее раскрывается в материале про экологическую нишу Хатчинсона: фактическая скорость роста зависит от того, насколько условия среды попадают в диапазон толерантности вида.

Время удвоения популяции

Удобная характеристика экспоненциального роста - время удвоения $T_2$, за которое численность вырастает вдвое. Из условия $N(T_2) = 2N_0$ следует $2 = e^{rT_2}$, откуда:

$T_2 = \frac{\ln 2}{r} \approx \frac{0{,}693}{r}$

Чем больше $r$, тем короче время удвоения. Эта формула отражает важное свойство экспоненты: время удвоения постоянно и не зависит от текущей численности - популяция удваивается через равные промежутки времени независимо от того, насколько она уже выросла. Аналогичным образом период полураспада постоянен в радиоактивном распаде, который описывается тем же уравнением с отрицательным $r$.

Разбор типовой задачи

Покажем, как модель работает на конкретных числах. Пусть в озеро запустили популяцию из $N_0 = 200$ особей рыбы, а удельная скорость роста оценена как $r = 0{,}15\ \text{год}^{-1}$. Требуется найти численность через 10 лет и время удвоения.

Подставляем в формулу $N(t) = N_0 e^{rt}$ значение $t = 10$:

$N(10) = 200 \cdot e^{0{,}15 \cdot 10} = 200 \cdot e^{1{,}5} \approx 200 \cdot 4{,}48 \approx 896$

За десять лет популяция выросла почти в 4,5 раза. Время удвоения находим отдельно:

$T_2 = \frac{\ln 2}{r} = \frac{0{,}693}{0{,}15} \approx 4{,}6\ \text{года}$

То есть каждые примерно 4,6 года численность удваивается: 200, затем около 400, около 800 и так далее. Обратите внимание, что результат экспоненциального роста сильно чувствителен к точности оценки $r$: ошибка в коэффициенте на 0,02 за десять лет меняет прогноз на сотни особей. Поэтому в полевых исследованиях $r$ стараются оценивать по нескольким замерам, а не по двум точкам.

Допущения и идеализации модели

Модель Мальтуса экспоненциального роста опирается на ряд сильных упрощений, которые важно осознавать при её применении:

• Неограниченные ресурсы. Предполагается, что пищи, пространства и других ресурсов хватает на любую численность - отсюда отсутствие верхнего предела.
• Постоянство $r$. Коэффициент роста считается неизменным во времени, хотя в реальности он колеблется из-за сезонности, болезней и изменений среды.
• Однородность популяции. Все особи считаются одинаковыми по плодовитости и выживаемости, без учёта возрастной и половой структуры.
• Отсутствие взаимодействий. Нет хищников, паразитов и межвидовой конкуренции, влияющих на динамику.

Эти допущения делают модель прозрачной для анализа, но именно из-за них она применима лишь на коротких интервалах или при низкой плотности. Снятие первого допущения - конечность ресурсов - и приводит к логистической модели, рассмотренной ниже.

Сравнение с логистической моделью

Главный недостаток модели Мальтуса - бесконечный рост: формула $N_0 e^{rt}$ при $t \to \infty$ даёт неограниченно большую численность, что физически невозможно. Реальные ресурсы среды конечны, и при росте плотности усиливается конкуренция, падает рождаемость, растёт смертность. Эти эффекты учитывает логистическая модель Ферхюльста, в которой скорость роста замедляется по мере приближения к ёмкости среды $K$:

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

При малых $N$ множитель $(1 - N/K)$ близок к единице, и логистическая кривая практически совпадает с мальтузианской экспонентой - то есть модель Мальтуса работает как хорошее приближение на ранней стадии роста. По мере приближения к $K$ рост тормозится, и кривая выходит на горизонтальную асимптоту.

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

Где применяется модель

Несмотря на упрощённость, модель Мальтуса экспоненциального роста остаётся рабочим инструментом в нескольких ситуациях:

• Начальная фаза роста любой популяции, когда плотность ещё низкая и ресурсы не лимитируют. Колония бактерий в свежей питательной среде, вид-вселенец на новой территории.
• Демография - модель применима к человечеству на коротких интервалах и исторически послужила Мальтусу основой прогноза о нехватке продовольствия.
• Радиоактивный распад и химическая кинетика первого порядка - то же уравнение с $r < 0$.
• Финансовая математика - непрерывное начисление сложных процентов описывается формулой $A_0 e^{rt}$, идентичной мальтузианской.

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

• Путают $r$ и $\lambda$. Удельная скорость $r$ входит в показатель экспоненты, а $\lambda = e^r$ - множитель за шаг. Подстановка $\lambda$ вместо $r$ в формулу $N_0 e^{rt}$ даёт завышенный результат.
• Считают рост линейным. За равные интервалы прирост не одинаков по абсолютной величине: он сам растёт пропорционально численности.
• Забывают про размерность времени. Если $r$ задан в год$^{-1}$, то и $t$ должно быть в годах, иначе результат бессмысленен.
• Применяют модель к большим $t$. При длительном прогнозе экспонента даёт нереально огромные числа - нужно переходить к логистической модели.

FAQ

Чем модель Мальтуса отличается от геометрической прогрессии? Это две формы одной идеи. Геометрическая прогрессия описывает рост в дискретном времени ($N_{t+1} = \lambda N_t$), а модель Мальтуса - в непрерывном ($N(t) = N_0 e^{rt}$). При малом шаге обе дают близкие значения.

Что означает отрицательный коэффициент $r$? Отрицательный $r$ возникает, когда смертность превышает рождаемость ($d > b$). Тогда $N(t) = N_0 e^{rt}$ убывает к нулю - популяция экспоненциально вымирает.

Можно ли по двум замерам численности найти $r$? Да. Из $N_1$ и $N_2$, измеренных через интервал $\Delta t$, скорость роста равна $r = \frac{1}{\Delta t}\ln\frac{N_2}{N_1}$. Эта оценка лежит в основе обработки данных полевых учётов.

Коротко

Модель Мальтуса экспоненциального роста описывает популяцию уравнением $\frac{dN}{dt} = rN$ с решением $N(t) = N_0 e^{rt}$, где $r$ - удельная скорость роста. Численность удваивается за постоянное время $T_2 = \ln 2 / r$ и формально стремится к бесконечности. Это идеализация для условий неограниченных ресурсов: она точна на ранней стадии и служит отправной точкой, от которой переходят к логистической модели с конечной ёмкостью среды.