学径XUEJING · 个人学习平台通知设置

课程库 Calculus 2 S5 · 一阶常微分方程 › 第 7

🎯 人口模型专题 —— 从指数噩梦到渔场崩溃

项目课55 分钟

本课导读

项目课:前六课的兵器全体出列。我们沿着人口建模的历史脉络打四场战役——马尔萨斯模型(指数增长的噩梦)、带收割的马尔萨斯(第一道生死线)、Logistic 模型(S 形曲线正式加冕)、带收割的 Logistic(渔场为什么会在"看起来还行"的捕捞量下突然崩溃)。每一场都是同一套组合拳:建模 → 求平衡点定稳定性 → 必要时求精确解 → 用人话说出结论。预计 55 分钟。

说明

前置工具箱:分离变量(第 2 课)、伯努利换元(第 4 课,y=yy2y' = y - y^2 那题就是本课主角的排练)、相线图与导数判别法(第 5-6 课)、一元二次方程求根公式(今天的隐藏主角)。

🔧 30 秒回顾 · 本课要用的旧手艺(熟练可跳过)

① 求根公式与判别式aP2+bP+c=0P=b±b24ac2aaP^2 + bP + c = 0 \Rightarrow P = \dfrac{-b \pm \sqrt{b^2 - 4ac}}{2a};判别式缩到 0 ⇔ 两根相向而行、合并成一个。微例:P21000P+187500=0P=1000±5002P^2 - 1000P + 187500 = 0 \Rightarrow P = \dfrac{1000 \pm 500}{2}——例题 3 原式。⚠️ 翻车点:bb 本身带负号时 b-b 的符号。

② 用 ln\ln 解指数方程ekt=Akt=lnAe^{kt} = A \Leftrightarrow kt = \ln Aln(ab)=blna\ln(a^b) = b\ln a。微例:e3k=2k=ln23e^{3k} = 2 \Rightarrow k = \dfrac{\ln 2}{3}——任务 1 原式。⚠️ 翻车点:ln(a+b)\ln(a+b) 没有展开式,别硬拆。

③ 本单元兵器速览:分离变量(第 2 课)、伯努利换元 v=y1nv = y^{1-n}(第 4 课)、相线图与 ff' 判别法(第 5-6 课)——四幕战役全靠它们,生疏哪个回哪课。

第一幕 · 马尔萨斯模型:增长率与人口成正比

最朴素的假设:人越多,生得越多——增长率正比于当前人口 PP

dPdt=kP(k>0)\frac{dP}{dt} = kP \qquad (k > 0)

第 2 课例题 1 的翻版,分离变量立得 P=P0ektP = P_0 e^{kt}——指数增长。马尔萨斯 1798 年据此推出"末日论"(人口指数增长、粮食至多线性增长,饥荒是数学必然),这个模型因此也叫末日模型

它短期内相当准(资源充足的细菌、疫情初期的感染数都走这条曲线),长期必然破产:指数增长意味着"每过固定时间翻一番",任何有限的地球都装不下。模型的失败之处恰恰标出了改进方向——资源上限,第三幕处理它。相线图一秒画完:唯一平衡点 P=0P^* = 0f(0)=k>0f'(0) = k > 0,不稳定——只要有一点人口,就一路奔向无穷。

第二幕 · 末日模型 + 收割:第一道生死线

给种群加上恒定速率的收割(捕捞、狩猎、移民流出),每单位时间拿走 hh

dPdt=kPh\frac{dP}{dt} = kP - h

例题 1 · 定性 + 定量双杀

第 1 步 · 平衡点与稳定性

f(P)=kPh=0f(P) = kP - h = 0 给出 P=hkP^* = \dfrac{h}{k}f(P)=k>0f'(P^*) = k > 0——不稳定,山顶的球。相线:P>h/kP > h/kf>0f > 0 ↑(增长跑赢收割,奔向无穷);P<h/kP < h/kf<0f < 0 ↓(收割吃穿老本,坠向灭绝)。

P=h/kP^* = h/k 是一条生死线:种群命运由初值站在线的哪一侧决定,而且这个平衡自身留不住——第 6 课"不稳定平衡是纸面状态"的活例。

第 2 步 · 精确解验证图景

方程是一阶线性(也可分离)。用积分因子 μ=ekt\mu = e^{-kt}(ektP)=hekt(e^{-kt}P)' = -he^{-kt},积分得 ektP=hkekt+Ce^{-kt}P = \dfrac{h}{k}e^{-kt} + C

P(t)=hk+(P0hk)ektP(t) = \frac{h}{k} + \left(P_0 - \frac{h}{k}\right)e^{kt}

读解:偏差 (P0hk)\left(P_0 - \frac{h}{k}\right)ekte^{kt} 指数放大(对照咖啡的 ekte^{-kt} 指数遗忘——符号一换,天堂地狱)。P0P_0 略低于生死线时,括号为负、放大到 P=0P = 0有限时间灭绝。定性图景逐条兑现。

第三幕 · Logistic 模型:给增长装上天花板

修正马尔萨斯:环境有承载力 KK(能养活的最大数量),越接近 KK,增长越乏力。最简洁的实现——给 kPkP 乘一个"剩余空间比例"因子:

dPdt=rP(1PK)\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right)

这就是 Logistic 方程。两个极端自检:PKP \ll K 时括号 1\approx 1,退化为马尔萨斯(空旷时代随便长);PKP \to K 时括号 0\to 0,增长熄火(满员了)。

例题 2 · 三板斧齐下

第 1 步 · 定性分析(30 秒出结论)

f(P)=rP(1P/K)=0f(P) = rP(1 - P/K) = 0:平衡点 P=0P^* = 0P=KP^* = K。求导 f(P)=r2rPKf'(P) = r - \dfrac{2rP}{K}f(0)=r>0f'(0) = r > 0 不稳定f(K)=r<0f'(K) = -r < 0 稳定

结论已经能写报告了:任何正初值的种群终将稳定在承载力 KK。灭绝态 00 是送客的,碰不得一点生机。

第 2 步 · 精确解(第 4 课的排练正式登台)

这是伯努利方程(n=2n = 2):展开 PrP=rKP2P' - rP = -\dfrac{r}{K}P^2,令 v=1/Pv = 1/P 走一遍第 4 课流程(那里解的 y=yy2y' = y - y^2 就是 r=K=1r = K = 1 的特例),得

P(t)=K1+Cert,C=KP0P0P(t) = \frac{K}{1 + Ce^{-rt}}, \qquad C = \frac{K - P_0}{P_0}

P0<KP_0 < KC>0C > 0:标准 S 形曲线——起步慢(基数小)、中段快(火力全开)、封顶慢(空间挤)。且 tt \to \inftyPKP \to K,与第 1 步的定性结论互相盖章。

第四幕 · Logistic + 收割:渔场崩溃之谜

压轴戏。给 Logistic 种群加恒定捕捞量 hh

dPdt=rP(1PK)h\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right) - h

例题 3 · 一个渔场的体检报告

某渔场 r=0.4r = 0.4/年,K=1000K = 1000 吨,年捕捞 h=75h = 75 吨。这个方案可持续吗?安全边际有多大?

第 1 步 · 求平衡点(二次方程登场)

0.4P(1P1000)75=00.4P\left(1 - \dfrac{P}{1000}\right) - 75 = 0,两边乘 10000.4=2500\dfrac{1000}{0.4} = 2500 并整理(🔧 工具箱 ①:求根公式):

P21000P+187500=0P=1000±1067500002=1000±5002P^2 - 1000P + 187500 = 0 \quad\Rightarrow\quad P = \frac{1000 \pm \sqrt{10^6 - 750000}}{2} = \frac{1000 \pm 500}{2}

两个平衡点:P=250P^* = 250P=750P^* = 750

第 2 步 · 判稳定性,读体检报告

f(P)=0.40.0008Pf'(P) = 0.4 - 0.0008Pf(750)=0.2<0f'(750) = -0.2 < 0 稳定f(250)=+0.2>0f'(250) = +0.2 > 0 不稳定

报告:鱼群会稳定运行在 750 吨(可持续 ✓);但 250 吨是暗礁——一场赤潮、一次偷捕把存量砸破 250,系统就滑进 f<0f < 0 的深渊,再也回不来。收割把 Logistic 原本"任何正初值都得救"的宽厚,改成了"跌破门槛就灭绝"的冷酷。

第 3 步 · 临界捕捞量:崩溃为什么是“突然”的

调大 hh,判别式 K24hKrK^2 - \dfrac{4hK}{r} 缩小——两个平衡点相向而行(稳定点降、暗礁升),在判别式归零处相撞湮灭(第 6 课练习 3 的剧情原封重演)。临界值:

hc=rK4=0.4×10004=100 吨/年h_c = \frac{rK}{4} = \frac{0.4 \times 1000}{4} = 100 \text{ 吨/年}

hh 越过 100 的一瞬间,平衡点消失、f(P)<0f(P) < 0 处处成立:无论存量多雄厚,灭绝是唯一结局。最阴险的是崩溃前毫无渐变的预警——h=99h = 99 时系统看起来还稳稳地运行着(平衡点仍在),h=101h = 101 时大势已去。这正是历史上多个渔业(如纽芬兰鳕鱼场)"前一年还有鱼、后一年就没了"的数学机制。

✍️ 练习(顾问任务)

任务 1 · 细菌的马尔萨斯时刻

培养皿中细菌每 3 小时翻一番。写出模型并求:初始数量增至 10 倍需要多久?

✅ 参考答案

P=P0ektP = P_0 e^{kt},翻倍条件 e3k=2e^{3k} = 2k=ln23k = \dfrac{\ln 2}{3}。求 ekt=10e^{kt} = 10t=ln10k=3ln10ln29.97t = \dfrac{\ln 10}{k} = \dfrac{3\ln 10}{\ln 2} \approx 9.97 小时。指数增长的直觉冲击:翻一番要 3 小时,翻十倍只要不到 10 小时(因为 1023.3210 \approx 2^{3.32}——十倍不过是三番多一点)。

任务 2 · Logistic 实战

池塘 r=1r = 1/年、K=100K = 100 条,投放 P0=10P_0 = 10 条鱼。写出 P(t)P(t),并求鱼量达到 50 条(半程)的时间。

✅ 参考答案

C=KP0P0=9010=9C = \dfrac{K - P_0}{P_0} = \dfrac{90}{10} = 9,故 P(t)=1001+9etP(t) = \dfrac{100}{1 + 9e^{-t}}。令 P=50P = 501+9et=21 + 9e^{-t} = 2et=19e^{-t} = \dfrac{1}{9}t=ln92.2t = \ln 9 \approx 2.2 年。附赠一个 S 形的小秘密:P=K/2P = K/2 恰是增长率 f(P)f(P) 的最大值点(对 ff 求导可证)——半满时长势最猛。

任务 3 · 新渔场评估

渔场 r=0.6r = 0.6/年、K=2000K = 2000 吨、拟定年捕捞 h=288h = 288 吨。求平衡点、稳定性、临界捕捞量,并给一句话结论。

✅ 参考答案

平衡方程两边乘 20000.6\dfrac{2000}{0.6}P22000P+960000=0P^2 - 2000P + 960000 = 0,判别式 4×1063.84×106=1600004\times10^6 - 3.84\times10^6 = 160000=400\sqrt{\cdot} = 400P=800P^* = 800(不稳定暗礁)与 P=1200P^* = 1200(稳定运行位,f=0.60.0006Pf' = 0.6 - 0.0006Pf(1200)=0.12<0f'(1200) = -0.12 < 0 ✓)。临界值 hc=0.6×20004=300h_c = \dfrac{0.6 \times 2000}{4} = 300 吨/年。

结论:方案可行但逼近红线(288 对 300,安全边际 4%),暗礁高达 800 吨——建议下调捕捞量为系统留出扰动缓冲。

📝 随堂测验

随堂测验0 / 5 题正确

1. 马尔萨斯模型 $P' = kP$ 的解与根本缺陷是?

2. $P' = kP - h$ 的平衡点 $P^* = h/k$ 的性质是?

3. Logistic 方程 $P' = rP(1 - P/K)$ 中因子 $(1 - P/K)$ 的建模含义是?

4. 例题 3 渔场(稳定点 750、暗礁 250)里,“暗礁”的含义是?

5. 捕捞量 $h$ 越过临界值 $h_c = rK/4$ 时发生了什么?

本课小结

  • 组合拳定式:建模 → 平衡点 + 稳定性 → (必要时)精确解 → 人话结论——四幕全按它打
  • 马尔萨斯 P=kPP' = kP:指数增长,短期神器、长期破产;加收割后出现不稳定生死线 h/kh/k
  • Logistic P=rP(1P/K)P' = rP(1-P/K)00 不稳、KK 稳,S 形曲线;P=K/2P = K/2 处长势最猛
  • Logistic + 收割:稳定运行位 + 不稳定暗礁并存;hc=rK/4h_c = rK/4 一过,两点湮灭、必然灭绝、没有渐变预警
  • 定性分析主导全程,精确解只在需要"多快"时出手——工具各司其职

下一课换一个舞台但不换武功:混合问题——水箱、盐水、进进出出。你会发现药物代谢、房间通风、湖泊治污全是同一道题换了戏服,而主角又是积分因子。

划选正文任意文字可高亮、批注或加入复习卡

讨论

载入中…