# 逻辑斯蒂模型与微分方程（仅BC要求） > AP 微积分 BC · AP 微积分 BC CED 第7单元 > 来源: https://www.owlsprep.com/zh/study/ap-calculus-bc-u7-logistic-models-with-differential-equations/ 这是AP微积分BC专属的子主题，讲解使用逻辑斯蒂微分方程进行有界增长建模。你将学习如何识别模型参数、求解可分离的逻辑斯蒂方程、分析增长率、凹凸性、拐点，并在AP考试题目中解读结果。 **先修:** [一阶微分方程的变量分离法](https://www.owlsprep.com/zh/study/ap-calculus-bc-u7-separation-of-variables/); 指数增长建模; 对数法则和部分分式分解 ## 学习目标 - 从逻辑斯蒂微分方程中识别环境容纳量和内禀增长率 - 求解并分类逻辑斯蒂微分方程的平衡解 - 求解逻辑斯蒂微分方程的通解和特解 - 分析逻辑斯蒂解的凹凸性并确定拐点位置 - 在实际情境问题中解释逻辑斯蒂模型 ## 逻辑斯蒂微分方程与平衡解 逻辑斯蒂模型修正了指数增长不切实际的无限增长假设，加入了有限资源限制的考量。它用于建模人口、疾病传播、技术普及以及其他趋近于最大可持续规模（称为环境容纳量$K$）的量。该主题占AP微积分BC考试总分的2-3%，会出现在选择题和自由问答题两个部分中。 **标准逻辑斯蒂微分方程** — 其中$P(t)$是$t$时刻的研究对象数量，$k>0$是内禀增长率（种群数量极小时的人均增长率），$K>0$是环境容纳量（$P$的最大可持续规模）。 *记法:* \frac{dP}{dt} = kP\left(1 - \frac{P}{K}\right) *例:* 当$P$非常小时，$\frac{dP}{dt} \approx kP$，符合指数增长。当$P=K$时，增长完全停止；当$P>K$时，$P$会向$K$减小。 平衡解是增长率$\frac{dP}{dt}=0$的常数解。对于逻辑斯蒂方程，两个平衡解是$P=0$（不稳定：任何小的正$P$都会增长偏离0）和$P=K$（稳定：任何非零$P$在$t \to \infty$时都会趋近于$K$）。增长率$\frac{dP}{dt}$是关于$P$的开口向下的二次函数，因此其最大值出现在$P = \frac{K}{2}$，即当数量达到环境容纳量的一半时增长最快。 **例题:** 给定微分方程$\frac{dP}{dt} = 0.3P\left(1 - \frac{P}{1800}\right)$，(a) 识别环境容纳量和内禀增长率，(b) 求出所有平衡解并分类其稳定性，(c) 求出使增长率最大的$P$值。 1. 直接将给定方程与标准形式$\frac{dP}{dt}=kP\left(1-\frac{P}{K}\right)$对比得到参数： 2. $$k=0.3 \quad (\text{intrinsic growth rate}), \quad K=1800 \quad (\text{carrying capacity})$$ 3. 令$\frac{dP}{dt}=0$求平衡解： 4. $$0.3P\left(1-\frac{P}{1800}\right)=0 \implies P=0 \text{ or } P=1800$$ 5. 当$00$，因此$P$会偏离0向1800移动。当$P>1800$时，$\frac{dP}{dt}<0$，因此$P$会向1800移动。因此，$P=0$是不稳定的，$P=1800$是稳定的。 6. 根为$P=0$和$P=1800$的二次函数的最大值出现在根的中点： 7. $$P = \frac{0 + 1800}{2} = 900$$ > **考试提示:** 识别$K$之前，一定要先将逻辑斯蒂方程改写为括号内首项为1的标准形式。不要从非标准因式分解形式中猜测$K$——这是选择题中最常见的错误。 ## 求解逻辑斯蒂微分方程 逻辑斯蒂微分方程是可分离的，因此我们可以使用变量分离法和部分分式分解求出$P(t)$的显式解。 **推导:** 推导逻辑斯蒂微分方程的通解和特解 *起点:* 标准逻辑斯蒂微分方程：$\frac{dP}{dt} = kP\left(1 - \frac{P}{K}\right)$ 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{KP_0}{P_0 + (K-P_0)e^{-kt}}$。当$t \to \infty$时，$P(t) \to K$，验证了$P=K$是稳定平衡。 **例题:** 一个种群按照$\frac{dP}{dt} = 0.2P\left(1 - \frac{P}{600}\right)$增长，初始种群数量为$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. 这与给定的初始值一致，确认解正确。 > **考试提示:** 一定要代入$t=0$验证你的特解，确认得到$P_0$。这能发现90%整理项时常见的代数错误。 ## 逻辑斯蒂解的凹凸性与拐点 为了分析逻辑斯蒂解曲线$P(t)$的形状，我们使用链式法则计算二阶导数$\frac{d^2P}{dt^2}$，以识别凹凸性和拐点。 **推导:** 求逻辑斯蒂微分方程中$P(t)$的二阶导数 *起点:* $\frac{dP}{dt} = kP\left(1 - \frac{P}{K}\right)$ 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^2P}{dt^2} = k\left(1-\frac{2P}{K}\right)\frac{dP}{dt} = k^2P\left(1-\frac{P}{K}\right)\left(1-\frac{2P}{K}\right)$ 当$00$，因此$\frac{d^2P}{dt^2}$的符号与$\left(1-\frac{2P}{K}\right)$的符号一致。当$P<\frac{K}{2}$时，$\frac{d^2P}{dt^2}>0$，因此$P(t)$是凹向上的。当$P>\frac{K}{2}$时，$\frac{d^2P}{dt^2}<0$，因此$P(t)$是凹向下的。这说明逻辑斯蒂曲线的拐点（凹凸性改变的点）恰好出现在$P=\frac{K}{2}$处，这里也是增长率最大的位置。如果初始种群$P_0>K$，$P(t)$会向$K$减小，永远不会穿过$P=\frac{K}{2}$，因此$t>0$时没有拐点。 **例题:** 对于逻辑斯蒂模型$\frac{dP}{dt} = 0.3P\left(1-\frac{P}{900}\right)$，初始条件$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. 根据二阶导数的符号说明凹凸区间：当$00$，$P(t)$凹向上。当$t>10\ln 2$时，$P>450$，因此$\frac{d^2P}{dt^2}<0$，$P(t)$凹向下。 > **考试提示:** 记住我们求的是关于$t$的函数$P(t)$的凹凸性，不是关于$P$的函数$\frac{dP}{dt}$的凹凸性。一定要使用链式法则求对$t$的二阶导数，不要停留在求$\frac{dP}{dt}$对$P$的导数。 ## AP风格例题 **例题:** 选择题：一个种群按照$\frac{dP}{dt} = 0.04P(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{dC}{dt} = 0.5C\left(1-\frac{C}{800}\right)$，初始条件$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人尝试过这款新的季节性饮品。 ## 常见错误 - **错误做法:** 对于$\frac{dP}{dt} = 0.2P(1000 - 2P)$，将环境容纳量算为$K=1000$，而非正确的$K=500$。 - 原因: 学生没有将方程改写为标准形式，直接认为因式分解二项式中的常数项就是$K$。 - 正确做法: 一定要从二项式中提取$P$的系数，得到首项为常数的形式，此时$P$下方的分母就是$K$。本例中$0.2P(1000-2P) = 200P\left(1-\frac{P}{500}\right)$，因此$K=500$。 - **错误做法:** 积分后，错误地将$\ln|P| - \ln|K-P|$合并为$\ln|P(K-P)|$。 - 原因: 学生在快速完成部分分式步骤时记错了对数法则。 - 正确做法: 取指数前先写出对数法则$\ln A - \ln B = \ln\left(\frac{A}{B}\right)$确认。 - **错误做法:** 求解常数$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{kK}{4}$，也就是$P=K/2$处$\frac{dP}{dt}$的值，不是$K$。 ## 速查表 | 类别 | 公式 | 说明 | | --- | --- | --- | | 标准逻辑斯蒂微分方程 | \frac{dP}{dt} = kP\left(1-\frac{P}{K}\right) | $P(t)$ = 时刻$t$的数量，$k$ = 内禀增长率，$K$ = 环境容纳量 | | 平衡解 | $P=0$ (unstable), $P=K$ (stable) | 出现在$\frac{dP}{dt}=0$处；所有非零解在$t\to\infty$时趋近于$P=K$ | | 最大增长率 | Maximum $\frac{dP}{dt} = \frac{kK}{4}$ at $P=\frac{K}{2}$ | 仅适用于$0K$则无拐点 | | 长期行为 | $\lim_{t\to\infty} P(t) = K$ | 对于所有$0