ロジスティック方程式
Contents
ロジスティック方程式
ロジスティック方程式 (英:logistic equation) とは、マルサスモデルのマルサス係数に、個体数に応じて増加率が減少するモデル $r\left(1-\frac{P}{K}\right)$ を適用したもの。
\[ \frac{dP}{dt} = rP\left(1-\frac{P}{K}\right) \]
パラメタ | 値域 | 説明 |
---|---|---|
$P$ | $P\in\R, ~0\le P\le K$ | 個体数 |
$K$ | $K\in\R, ~P\le K$ | 環境収納力 |
$r$ | $r\in\R$ | 内的自然増加率 |
$t$ | $t\in\R$ | 時間 |
ロジスティック方程式の解
一般に $t=0$ 時の $P$ の値である $P_0$ の値をパラメタとした次の解がロジスティック方程式の解として知られている。
\[ P = \frac{P_0Ke^{rt}}{K-P_0+P_0e^{rt}} \]
ロジスティック方程式の解の導出: \( \begin{aligned} \frac{dP}{dt} &= rP\left(1-\frac{P}{K}\right) \cr rdt &= \frac{K}{P(K-P)}dP \end{aligned} \)
部分分数分解により、
\[ \begin{aligned} rdt &= \frac{K}{P(K-P)}dP \cr &= \left(\frac{1}{P}+\frac{1}{K-P}\right)dP \end{aligned} \]
両辺を積分すると、
\[ \begin{aligned} \int{rdt} &= \int{\left(\frac{1}{P}+\frac{1}{K-P}\right)dP} \cr &= \int{\frac{1}{P}dP} + \int{\frac{1}{K-P}dP} \cr \cr rt &= \ln P - \ln(K-P) + C \quad \because P,K>0, ~P\le K \cr rt - C &= \ln\frac{P}{K-P} \cr \exp(rt-C) &= \frac{P}{K-P} \cr \cr P &= \exp(rt-C)(K-P) \cr &= K\exp(rt-C)-P\exp(rt-C) \cr \cr K\exp(rt-C) &= P(1+\exp(rt-C)) \cr P &= \frac{K\exp(rt-C)}{1+\exp(rt-C)} \cr &= \frac{K}{1+\exp(-rt+C)} \end{aligned} \]
$P_0$ を計算すると、
\[ \begin{aligned} P_0 &= \left.\frac{K}{1+\exp(-rt+C)}\right|_{t=0} \cr &= \frac{K}{1+\exp(C)} \cr 1+\exp(C) &= \frac{K}{P_0} \cr \exp(C) &= \frac{K-P_0}{P_0} \end{aligned} \tag{1} \]
よって、
\[ \begin{aligned} P &= \frac{K}{1+\exp(-rt+C)} \cr &= \frac{K}{1+\exp(C)\exp(-rt)} \cr &= \frac{K}{1+\frac{K-P_0}{P_0}\exp(-rt)} \quad \because (1) \cr &= \frac{KP_0\exp(rt)}{\left[1+\frac{K-P_0}{P_0}\exp(-rt)\right]P_0\exp(rt)} \cr &= \frac{P_0K\exp(rt)}{K-P_0+P_0\exp(rt)} \cr \cr \therefore P &= \frac{P_0Ke^{rt}}{K-P_0+P_0e^{rt}} \end{aligned} \]