Интернет-журнал "Домашняя лаборатория", 2008 №4
Шрифт:
Преобразуем наше выражение в следующий вид:
(N(t + ?t) — N(t))/?t = (b — d)N(t)
Далее видим, что в левой части стоит производная N при ?t, устремленной к нулю:
Значит, имеем формулу:
dN/dt = (b — d)N(t)
Или
dN/dt = rN(t)
Таким образом, мы получили основную формулу для модели экспоненциального
На следующих графиках представлена зависимость плотности популяции индеек (англ. turkeys на графике) в штате Мичиган от времени, по оси 0Y отложены N в левом и In N во правом графиках. К концу IX века всех индеек истребили, но скоро они опять появились, сначала это было небольшое количество особей, но затем популяция стала расти экспоненциально, что можно заметить на графиках.
Обычно, так ведет себя сравнительно небольшая популяция, расселившаяся на обширном пространстве. Но на самом деле в природе ни одна популяция по крайней мере долгое время экспоненциально не растет, то есть эта модель слишком упрощенная. Поэтому мы попробуем добавлять в наше уравнение некоторые ограничения, которые будут приближать нашу модель к реальной.
Заметим, что r не может быть константой, это функция от числа особей в популяции, ведь, как мы помним, r = b — d, a b и d зависит от числа особей, конкуренции внутри популяции, хищников и др. факторов, ведь чем больше особей в популяции, тем больше животные конкурируют за еду, пространство и т. д., и тем труднее становится им выживать. Чем больше N, тем меньше r, поэтому имеет место формула:
Здесь мы сделали самое простое предположение, что r зависит линейно от N. Величины r0 и К — начальные параметры, некоторые константы, причем К — емкость среды, то есть некоторое стабилизирующее значение числа особей, при котором устанавливается равновесие, то есть популяция перестает расти. Если же N больше К, то популяция наоборот начинает убывать, r0 — это удельная скорость популяционного роста при малых N. Подставляя это значение r в модель экспоненциального роста, мы получаем логистическую модель, которая выражает второй закон популяционной динамики, говорящий о том, что всегда существует предел экспоненциального роста, накладываемый средой.
Третья модель — это модель Лотки-Вольтерра. Но сначала поговорим о такой важной концепции, которая называется функциональной зависимостью. Это зависимость между скоростью, с которой хищник убивает жертв, и плотностью количества жертв.
Для наглядности проведем параллель с химической кинетикой. Представим, что у нас есть пробирка с находящимся внутри раствором какого-то реагента А. Бросим внутрь одну молекулу вещества В. Частицы будут участвовать в броуновском движении и, рано или поздно, молекула В встретится с молекулой А. Вероятность этого события за единицу времени, очевидно, пропорциональна концентрации [А] вещества А. Это и есть пример функциональной зависимости. Экологи провели массу экспериментов, в которых вещество А «представляли» жертвы, а вещество В — хищники (к примеру, запускала в аквариум головастиков и одного тритона и смотрели, с какой скоростью тритон будет их поедать, в зависимости от количества головастиков). Итак, функциональная зависимость —
В экологии мы говорим о ресурсах и потребителях, поскольку не только хищники и жертвы, но и патогенны и их хозяева, паразиты и их хозяева, травоядные и их пища — это все примеры взаимодействия типа «хищник-жертва». В описывающей эту функциональную зависимость формуле f(R) = aR буквой R обозначены ресурсы. Эта зависимость — линейная. Это самый простой вид функциональной зависимости первого типа. Проблема заключается в том, что даже в химических реакциях это очень неточное приближение. В популяционной динамике это очень плохое приближение, потому что, когда мы начинаем увеличивать количество жертв, достаточно быстро наступает момент, когда хищник может убить жертву, но не может ее съесть. Есть некоторый потолок, выше которого хищник не будет убивать жертв, потому что ему это не будет нужно. Функциональная зависимость этот факт учитываклцая, называется функциональной зависимостью второго типа. Она описывается такой формулой:
f(R) = cR/(d + R)
Плотность жертв
В формуле с — это как раз тот самый потолок, к которому асимптотически приближается наша гипербола; d показывает, с какой скоростью она к нему приближается (количество жертв, которых убивает хищник со скоростью, равной с/2).
Рассматривая функциональные зависимости, мы постепенно приближались к модели Лотки-Вольтерра.
Эту модель мы будем выводить, используя так называемый «коробочный» подход (он проиллюстрирован на рисунке).
В этой системе слева записаны плотности ресурсов и потребителей. Справа вторая «коробка» у жертв отвечает за популяционный рост, у хищников — за уровень смерти, а первая — за взаимодействие ресурсов и потребителей. В эту общую модель мы можем подставлять разные функциональные зависимости. Хищники жертв поедают, поэтому перед первой коробкой в первом уравнении стоит знак минус, а во втором — плюс, т. к. при поедании жертв улучшается возможность увеличения количества хищников. Начнем с простого, перейдя потом к сложному. Подставив в эту модель самую простую функциональную зависимость (экспоненциальную), получим следующую систему уравнений:
dR/dt = rR — aRP
dP/dt = caRP — dP
где
• R — плотность жертвы (в килограммах мяса на квадратный километр),
• Р — плотность хищника
• r — удельная скорость роста популяции при малом числе хищников, т. е. удельная скорость, с которой хищники убивают жертв,
• а — эффективность хищника (удельная скорость, с которой хищник убивает жертв),
• с — эффективность превращения хищниками биомассы жертв в себя и потомство,
• d — удельный коэффициент смертности, с которой хищники вымирают, когда нет жертв.
Эта простейшая модель взаимодействия хищника и жертвы, или ресурсов и потребителей и она называется моделью Лотки-Вольтерра (в русском языке почему-то делают ударение на последнем слоге, ВольтеррА, наверное, потому что его первая статья вышла на французском языке, а он итальянец, и произносить надо ВольтЕрра). Лотка был американец польского происхождения. Они изобрели свои модели почти одновременно, один в 1925 году, а другой — в 1926.