Уравнение экспоненциального роста и бактерии

Я микробиолог, но я преподаю экологию на вводном курсе, поэтому, когда мы дошли до роста популяции, я подумал, что воспользуюсь примером микробной популяции. Но я столкнулся со странной проблемой, которую, как я думал, может помочь мне понять эколог:

Представьте себе популяцию бактерий, которые могут делиться каждые 30 минут. За час каждая клетка производит четыре клетки. Поскольку на самом деле никто не умирает, внутренняя скорость роста (r) равна 4. Уравнение экспоненциального роста dN/dt = rN прекрасно показывает рост популяции: начиная с одной клетки, через час это 4, затем через два часов rN = 4*4 = 16, через три часа rN = 16*4 = 64 и так далее. В 16 часов мы получаем около 4 миллиардов бактерий, чего и ожидает микробиолог.

Но затем мы говорим учащимся, что для облегчения предсказания будущих чисел мы можем немного посчитать и получить N(t) = N(0) e^rt. Если N(0) = 1, t = 16 часов и r = 4 бактерии/час/клетка, это дает e^64, что составляет примерно 6x10^27. Ой! Вроде далеко не 4 млрд.

Итак... почему это не работает? Мне кажется, что либо я неправильно понимаю r (думаю, чтобы это сработало, оно должно быть ниже примерно на 0,6), либо, может быть, у уравнения есть предел, и оно просто не работает для абсурдно высокого r бактерии?

Спасибо за помощь...

Ответы (3)

Учитывая ваше предположение:

Я просто смотрю на экспоненциальную часть, где работает простое экспоненциальное уравнение. Если мы предположим, что питательных веществ достаточно для бесконтрольного роста бактерий в течение нескольких часов (более или менее верно в реальной культуре),

В вашей исходной модели вы используете дискретные состояния и фиксированные временные шаги. Итак, если 30 минут - это один временной шаг, то после n шагов у вас есть 2 н или же 4 ( н / 2 ) клетки. Вы в основном считаете двойные шаги, и это нормально. Для непрерывного приближения вы должны оценить скорость следующим образом.

У вас есть уравнение роста первого порядка:

г Н г т знак равно р Н

Когда вы интегрируете это уравнение (определенный интеграл), вы получите:

л н ( Н ф Н я ) знак равно р т

Где Н ф конечное количество ячеек, Н я - начальное количество ячеек и т это интервал времени.


Когда номер ячейки удваивается, то Н ф Н я знак равно 2 а также т знак равно 30   мин (время удвоения).


Тогда ваша константа скорости будет: л н ( 2 ) 30 0,023  мин 1

Через 16 часов (16×30 минут) ваше количество клеток будет: е 0,023 × 16 × 60 4,27 × 10 9 (около 4 млрд).

Так что же в принципе не так? Вы не изменили масштаб своей ставки. Поскольку показатель степени е и ты это знаешь а Икс знак равно е Икс . л н ( а ) , все, что вам нужно сделать, это масштабировать константу скорости с помощью л н ( а ) (а=2). Более того, в непрерывной модели у вас нет фиксированных временных шагов. Таким образом, вы также масштабируете константу со временем удвоения.

@Jon Если этот ответ касается вашей проблемы , подумайте о том, чтобы принять его , нажав на галочку слева от ответа, сделав его зеленым. Это помечает вопрос как решенный к вашему удовлетворению и присуждает репутацию как вам, так и человеку, который ответил.

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

Логистическая модель

Здесь я быстро представляю одну стандартную модель логистического роста для непрерывного (а не дискретного) времени.

г Н г т знак равно р Н ( т ) ( 1 Н ( т ) К )

, куда К является грузоподъемностью. Из приведенной выше модели видно, что когда Н ( т ) является низким по сравнению с К , то дифференциальное уравнение хорошо аппроксимируется описанной вами экспоненциальной моделью г Н г т р Н ( т ) . Параметр г Н г т знак равно 0 показывает, что равновесие достигается, когда Н ( т ) знак равно 0 (неустойчивое равновесие) или Н ( т ) знак равно К (устойчивое равновесие).

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

Более продвинутые модели

Есть и другие более реалистичные модели. Например, реальные популяции бактерий могут также потреблять свои ресурсы со скоростью, превышающей скорость, с которой производятся ресурсы, и тогда размер популяции будет уменьшаться и в конечном итоге достигнет исчезновения.

Здесь вы найдете описание модели роста популяции, включающей популяции хищников.

Создайте собственную модель роста населения

Эти модели не так сложно собрать в соответствии с вашими потребностями. Первые главы « Руководства биолога по математическому моделированию в экологии и эволюции» (от Otto and Day) помогут вам построить и проанализировать собственные модели популяционной биологии, если вам это интересно.

Это полезно с точки зрения понимания общей картины, и я ценю ссылку на то, что выглядит действительно полезной книгой. Но для этого примера я просто смотрю на экспоненциальную часть, где работает простое экспоненциальное уравнение. Если мы предположим, что питательных веществ достаточно для бесконтрольного роста бактерий в течение нескольких часов (более или менее верно для реальной культуры), то я думаю, что уравнение dN/dt и уравнение er^rt должны давать ОДИНАКОВОЕ приближение популяции. размера в данный момент времени. Почему они дают такие разные результаты?
Либо ваш р слишком велик (но я думаю, что это не слишком плохое приближение) или нарушены предположения об отсутствии конкуренции за доступ к ресурсу.

Количество р должна быть частота, измеряемая в часах 1 . Поэтому его нельзя измерить в бактериях/час/клетку. Если количество бактерий умножается на 4 каждый час, это означает, что опыт р знак равно 4 , или же р знак равно п 4 , или 1, 39 в час.