逻辑斯蒂模型与微分方程（仅BC要求）

AP 微积分 BC· AP 微积分 BC CED — 微分方程· 14 分钟阅读

1. 逻辑斯蒂微分方程与平衡解★★☆☆☆⏱ 4 min

逻辑斯蒂模型修正了指数增长不切实际的无限增长假设，加入了有限资源限制的考量。它用于建模人口、疾病传播、技术普及以及其他趋近于最大可持续规模（称为环境容纳量 $K$ ）的量。该主题占AP微积分BC考试总分的2-3%，会出现在选择题和自由问答题两个部分中。

📘 定义

标准逻辑斯蒂微分方程

\frac{dP}{dt} = kP\left(1 - \frac{P}{K}\right)

其中 $P (t)$ 是 $t$ 时刻的研究对象数量， $k > 0$ 是内禀增长率（种群数量极小时的人均增长率）， $K > 0$ 是环境容纳量（ $P$ 的最大可持续规模）。

例:

当 $P$ 非常小时， $\frac{d P}{d t} \approx k P$ ，符合指数增长。当 $P = K$ 时，增长完全停止；当 $P > K$ 时， $P$ 会向 $K$ 减小。

平衡解是增长率 $\frac{d P}{d t} = 0$ 的常数解。对于逻辑斯蒂方程，两个平衡解是 $P = 0$ （不稳定：任何小的正 $P$ 都会增长偏离0）和 $P = K$ （稳定：任何非零 $P$ 在 $t \to \infty$ 时都会趋近于 $K$ ）。增长率 $\frac{d P}{d t}$ 是关于 $P$ 的开口向下的二次函数，因此其最大值出现在 $P = \frac{K}{2}$ ，即当数量达到环境容纳量的一半时增长最快。

📐 例题

给定微分方程 $\frac{d P}{d t} = 0.3 P (1 - \frac{P}{1800})$ ，(a) 识别环境容纳量和内禀增长率，(b) 求出所有平衡解并分类其稳定性，(c) 求出使增长率最大的 $P$ 值。

1
直接将给定方程与标准形式 $\frac{d P}{d t} = k P (1 - \frac{P}{K})$ 对比得到参数：
2
$k=0.3 \quad (\text{intrinsic growth rate}), \quad K=1800 \quad (\text{carrying capacity})$
3
令 $\frac{d P}{d t} = 0$ 求平衡解：
4
$0.3P\left(1-\frac{P}{1800}\right)=0 \implies P=0 \text{ or } P=1800$
5
当 $0 < P < 1800$ 时， $\frac{d P}{d t} > 0$ ，因此 $P$ 会偏离0向1800移动。当 $P > 1800$ 时， $\frac{d P}{d t} < 0$ ，因此 $P$ 会向1800移动。因此， $P = 0$ 是不稳定的， $P = 1800$ 是稳定的。
6
根为 $P = 0$ 和 $P = 1800$ 的二次函数的最大值出现在根的中点：
7
$P = \frac{0 + 1800}{2} = 900$

Exam tip:

识别 $K$ 之前，一定要先将逻辑斯蒂方程改写为括号内首项为1的标准形式。不要从非标准因式分解形式中猜测 $K$ ——这是选择题中最常见的错误。

2. 求解逻辑斯蒂微分方程★★★☆☆⏱ 4 min

逻辑斯蒂微分方程是可分离的，因此我们可以使用变量分离法和部分分式分解求出 $P (t)$ 的显式解。

🔬 推导

目标:

推导逻辑斯蒂微分方程的通解和特解

起点:

标准逻辑斯蒂微分方程： $\frac{d P}{d t} = k P (1 - \frac{P}{K})$

1
分离变量：
2
$\frac{dP}{P\left(1-\frac{P}{K}\right)} = k \, dt$
3
使用部分分式分解改写左侧：
4
$\left(\frac{1}{P} + \frac{1}{K-P}\right)dP = k \, dt$
5
对两边积分并合并对数：
6
$\ln|P| - \ln|K-P| = kt + C \implies \ln\left|\frac{P}{K-P}\right| = kt + C$
7
两边取指数，令 $A = e^{C}$ ：
8
$\frac{P}{K-P} = Ae^{kt}$
9
整理得到 $P (t)$ ，得到通解：
10
$P(t) = \frac{KAe^{kt}}{1 + Ae^{kt}}$
11
代入初始条件 $P (0) = P_{0}$ 求解得 $A = \frac{P _{0}}{K - P _{0}}$

结果:

初始条件 $P (0) = P_{0}$ 对应的特解为： $P (t) = \frac{K P _{0}}{P _{0} + ( K - P _{0} ) e ^{- k t}}$ 。当 $t \to \infty$ 时， $P (t) \to K$ ，验证了 $P = K$ 是稳定平衡。

📐 例题

一个种群按照 $\frac{d P}{d t} = 0.2 P (1 - \frac{P}{600})$ 增长，初始种群数量为 $P (0) = 150$ ，求特解 $P (t)$ 。

1
从微分方程和初始条件中识别参数：
2
$k=0.2, \quad K=600, \quad P_0=150$
3
代入特解公式：
4
$P(t) = \frac{KP_0}{P_0 + (K-P_0)e^{-kt}} = \frac{600 \cdot 150}{150 + (600-150)e^{-0.2t}}$
5
分子分母同除以150化简表达式：
6
$P(t) = \frac{90000}{150 + 450e^{-0.2t}} = \frac{600}{1 + 3e^{-0.2t}}$
7
验证初始条件检查是否存在错误：
8
$P(0) = \frac{600}{1+3} = 150$
9
这与给定的初始值一致，确认解正确。

Exam tip:

一定要代入 $t = 0$ 验证你的特解，确认得到 $P_{0}$ 。这能发现90%整理项时常见的代数错误。

3. 逻辑斯蒂解的凹凸性与拐点★★★☆☆⏱ 3 min

为了分析逻辑斯蒂解曲线 $P (t)$ 的形状，我们使用链式法则计算二阶导数 $\frac{d ^{2} P}{d t ^{2}}$ ，以识别凹凸性和拐点。

🔬 推导

目标:

求逻辑斯蒂微分方程中 $P (t)$ 的二阶导数

起点:

$\frac{d P}{d t} = k P (1 - \frac{P}{K})$

1
使用链式法则对 $t$ 两边求导：
2
$\frac{d^2P}{dt^2} = \frac{d}{dP}\left(kP\left(1-\frac{P}{K}\right)\right) \cdot \frac{dP}{dt}$
3
计算对 $P$ 的导数：
4
$\frac{d}{dP}\left(kP - \frac{kP^2}{K}\right) = k\left(1 - \frac{2P}{K}\right)$

结果:

$\frac{d ^{2} P}{d t ^{2}} = k (1 - \frac{2 P}{K}) \frac{d P}{d t} = k^{2} P (1 - \frac{P}{K}) (1 - \frac{2 P}{K})$

当 $0 < P < K$ 时， $\frac{d P}{d t} > 0$ ，因此 $\frac{d ^{2} P}{d t ^{2}}$ 的符号与 $(1 - \frac{2 P}{K})$ 的符号一致。当 $P < \frac{K}{2}$ 时， $\frac{d ^{2} P}{d t ^{2}} > 0$ ，因此 $P (t)$ 是凹向上的。当 $P > \frac{K}{2}$ 时， $\frac{d ^{2} P}{d t ^{2}} < 0$ ，因此 $P (t)$ 是凹向下的。这说明逻辑斯蒂曲线的拐点（凹凸性改变的点）恰好出现在 $P = \frac{K}{2}$ 处，这里也是增长率最大的位置。如果初始种群 $P_{0} > K$ ， $P (t)$ 会向 $K$ 减小，永远不会穿过 $P = \frac{K}{2}$ ，因此 $t > 0$ 时没有拐点。

📐 例题

对于逻辑斯蒂模型 $\frac{d P}{d t} = 0.3 P (1 - \frac{P}{900})$ ，初始条件 $P (0) = 100$ ，求拐点的 $t$ 坐标并说明凹凸区间。

1
初始值低于 $K$ 的解的拐点出现在 $P = \frac{K}{2}$ ：
2
$P = \frac{900}{2} = 450$
3
使用标准公式写出 $P (t)$ 的特解：
4
$P(t) = \frac{900 \cdot 100}{100 + 800e^{-0.3t}} = \frac{900}{1 + 8e^{-0.3t}}$
5
令 $P (t) = 450$ 求解 $t$ ：
6
$450 = \frac{900}{1+8e^{-0.3t}} \implies 1 + 8e^{-0.3t} = 2 \implies e^{-0.3t} = \frac{1}{8} \implies t = \frac{\ln 8}{0.3} = 10\ln 2 \approx 6.93$
7
根据二阶导数的符号说明凹凸区间：当 $0 < t < 10 ln 2$ 时， $P < 450$ ，因此 $\frac{d ^{2} P}{d t ^{2}} > 0$ ， $P (t)$ 凹向上。当 $t > 10 ln 2$ 时， $P > 450$ ，因此 $\frac{d ^{2} P}{d t ^{2}} < 0$ ， $P (t)$ 凹向下。

Exam tip:

记住我们求的是关于 $t$ 的函数 $P (t)$ 的凹凸性，不是关于 $P$ 的函数 $\frac{d P}{d t}$ 的凹凸性。一定要使用链式法则求对 $t$ 的二阶导数，不要停留在求 $\frac{d P}{d t}$ 对 $P$ 的导数。

4. AP风格例题★★★★☆⏱ 3 min

📐 例题

选择题：一个种群按照 $\frac{d P}{d t} = 0.04 P (800 - P)$ 增长，其环境容纳量是多少？

1
提取二项式中的800，将微分方程改写为标准逻辑斯蒂形式：
2
$\frac{dP}{dt} = 0.04P \cdot 800\left(1-\frac{P}{800}\right) = 32P\left(1-\frac{P}{800}\right)$
3
与标准形式对比， $K = 800$ 。32是内禀增长率 $k$ ，不是 $K$ ，因此正确答案是800。

📐 例题

实际情境自由问答题：一家咖啡店推出了新的季节性饮品，最多共有800名常客会尝试这款饮品。新顾客尝试的速率符合逻辑斯蒂模型 $\frac{d C}{d t} = 0.5 C (1 - \frac{C}{800})$ ，初始条件 $C (0) = 20$ （开业当天有20名顾客尝试了它）。问5天后共有多少顾客尝试过这款饮品，结果四舍五入到整数。

1
识别参数： $K = 800$ ， $k = 0.5$ ， $P_{0} = 20$ ，代入特解公式：
2
$C(t) = \frac{KP_0}{P_0 + (K-P_0)e^{-kt}} = \frac{800 \cdot 20}{20 + 780e^{-0.5t}} = \frac{16000}{20 + 780e^{-0.5t}}$
3
代入 $t = 5$ 计算：
4
$e^{-0.5(5)} = e^{-2.5} \approx 0.082085, \quad \text{denominator} \approx 20 + 780(0.082085) \approx 84.026$
5
求解 $C (5)$ 并解释：
6
$C(5) \approx \frac{16000}{84.026} \approx 190$
7
5天后，800名常客中大约有190人尝试过这款新的季节性饮品。

5. 常见陷阱

错误做法:

对于 $\frac{d P}{d t} = 0.2 P (1000 - 2 P)$ ，将环境容纳量算为 $K = 1000$ ，而非正确的 $K = 500$ 。

原因:

学生没有将方程改写为标准形式，直接认为因式分解二项式中的常数项就是 $K$ 。

正确做法:

一定要从二项式中提取 $P$ 的系数，得到首项为常数的形式，此时 $P$ 下方的分母就是 $K$ 。本例中 $0.2 P (1000 - 2 P) = 200 P (1 - \frac{P}{500})$ ，因此 $K = 500$ 。

错误做法:

积分后，错误地将 $ln ∣ P ∣ - ln ∣ K - P ∣$ 合并为 $ln ∣ P (K - P) ∣$ 。

原因:

学生在快速完成部分分式步骤时记错了对数法则。

正确做法:

取指数前先写出对数法则 $ln A - ln B = ln (\frac{A}{B})$ 确认。

错误做法:

求解常数 $A$ 时得到 $A = \frac{K - P _{0}}{P _{0}}$ ，而非正确的 $A = \frac{P _{0}}{K - P _{0}}$ 。

原因:

整理初始条件方程时简单地把分数颠倒了。

正确做法:

得到 $P (t)$ 后代入 $t = 0$ 确认得到 $P_{0}$ ，就能立刻发现这个错误。

错误做法:

对于初始值 $P_{0} > K$ 的解，断言其在 $P = K /2$ 处存在拐点。

原因:

学生记住了拐点永远在 $K /2$ ，但忘记如果 $P$ 初始值高于 $K$ ，它只会向 $K$ 减小，永远不会穿过 $K /2$ 。

正确做法:

断言存在拐点前，一定要根据初始条件检查 $P (t)$ 的取值范围。

错误做法:

题目问最大增长率，却回答了环境容纳量 $K$ 。

原因:

学生混淆了'最大种群规模'和'种群的最大增长率'。

正确做法:

如果题目问最大增长率，计算 $\frac{k K}{4}$ ，也就是 $P = K /2$ 处 $\frac{d P}{d t}$ 的值，不是 $K$ 。

6. 速查表

类别	公式	说明
标准逻辑斯蒂微分方程	\frac{dP}{dt} = kP\left(1-\frac{P}{K}\right)	$P (t)$ = 时刻 $t$ 的数量， $k$ = 内禀增长率， $K$ = 环境容纳量
平衡解	$P = 0$ (unstable), $P = K$ (stable)	出现在 $\frac{d P}{d t} = 0$ 处；所有非零解在 $t \to \infty$ 时趋近于 $P = K$
最大增长率	Maximum $\frac{d P}{d t} = \frac{k K}{4}$ at $P = \frac{K}{2}$	仅适用于 $0 < P < K$ ；与 $P (t)$ 的拐点位置一致
通解	$P (t) = \frac{K A e ^{k t}}{1 + A e ^{k t}}$	变量分离法和部分分式积分的结果
特解（ $P (0) = P_{0}$ ）	$P (t) = \frac{K P _{0}}{P _{0} + ( K - P _{0} ) e ^{- k t}}$	AP考试题目中最常用的形式
二阶导数（凹凸性）	$\frac{d ^{2} P}{d t ^{2}} = k^{2} P (1 - \frac{P}{K}) (1 - \frac{2 P}{K})$	由链式法则推导得到，用于求凹凸区间
拐点	Inflection at $P = \frac{K}{2}$	仅当解穿过 $P = \frac{K}{2}$ 时存在；若 $P_{0} > K$ 则无拐点
长期行为	$lim_{t \to \infty} P (t) = K$	对于所有 $0 < P_{0} < \infty, P_{0} \neq = 0$ 都成立

真题中的出现

AI 根据考纲规律估算的考点位置,请对照官方真题核实准确性。仅作复习重点参考。

2022 · BC
自由问答题：人口增长建模
2023 · BC MCQ
从微分方程识别环境容纳量

下一步

逻辑斯蒂模型建立在你对变量分离法和微分方程的知识基础上，它是AP微积分BC考试中应用建模里有界增长的核心例子。掌握逻辑斯蒂微分方程对于跨主题的自由问答题至关重要，这类题目会结合微分方程和实际建模，是考试中常见的题型。本主题也为分析平衡解的稳定性打下了基础，这是高等微积分课程的核心技能。掌握本主题后，你就可以学习更一般的自治微分方程、它们的平衡解，以及无法代数求解的方程的其他实用解法了。