跳转至

7.7.3 Solving Logistic Equations(求解逻辑斯蒂方程)

运用分离变量法求解逻辑斯蒂微分方程,得到解析解 P(t) = L/(1+Ae^(-kt))

定义

逻辑斯蒂方程(Logistic Equation)是描述种群增长的常微分方程,考虑了环境容纳量的限制。标准形式为 \(\frac{dP}{dt} = kP(L-P)\),其中 \(P(t)\) 表示时刻 \(t\) 的种群数量,\(k > 0\) 是增长率常数,\(L\) 是环境容纳量(承载量)。通过分离变量法求解该微分方程,可得到解析解 \(P(t) = \frac{L}{1+Ae^{-kLt}}\),其中 \(A\) 是由初始条件 \(P(0) = P_0\) 确定的常数,\(A = \frac{L-P_0}{P_0}\)。该解描述了一条S形曲线,初期近似指数增长,后期趋向于承载量 \(L\)

核心公式

  • \(\frac{dP}{dt} = kP(L-P)\)
  • \(P(t) = \frac{L}{1+Ae^{-kLt}}\)
  • \(A = \frac{L-P_0}{P_0}\),其中 \(P_0 = P(0)\) 是初始种群数量
  • \(\lim_{t \to \infty} P(t) = L\)
  • \(\frac{dP}{dt}\bigg|_{P=L/2} = \frac{kL^2}{4}\)(最大增长率出现在 \(P = L/2\) 处)

易错点

  • ⚠️ 混淆逻辑斯蒂方程的两种形式:\(\frac{dP}{dt} = kP(L-P)\)\(\frac{dP}{dt} = rP(1-\frac{P}{L})\)(两者等价但参数含义不同,\(r=kL\)),导致在求解或代入初值时出错
  • ⚠️ 在分离变量后进行部分分式分解时出错,特别是忘记处理 \(\frac{1}{P(L-P)} = \frac{1}{L}(\frac{1}{P} + \frac{1}{L-P})\) 中的系数 \(\frac{1}{L}\),或在积分后忘记指数化
  • ⚠️ 错误地确定常数 \(A\) 的值,例如直接用 \(P_0\) 代替 \(A\),或在计算 \(A = \frac{L-P_0}{P_0}\) 时分子分母颠倒
  • ⚠️ 误解解的长期行为,认为 \(P(t)\) 会超过 \(L\) 或无限振荡,而实际上 \(P(t)\) 单调趋向于 \(L\)(当 \(0 < P_0 < L\) 时)