本课导读
项目课:前六课的兵器全体出列。我们沿着人口建模的历史脉络打四场战役——马尔萨斯模型(指数增长的噩梦)、带收割的马尔萨斯(第一道生死线)、Logistic 模型(S 形曲线正式加冕)、带收割的 Logistic(渔场为什么会在"看起来还行"的捕捞量下突然崩溃)。每一场都是同一套组合拳:建模 → 求平衡点定稳定性 → 必要时求精确解 → 用人话说出结论。预计 55 分钟。
说明
前置工具箱:分离变量(第 2 课)、伯努利换元(第 4 课,y′=y−y2 那题就是本课主角的排练)、相线图与导数判别法(第 5-6 课)、一元二次方程求根公式(今天的隐藏主角)。
🔧 30 秒回顾 · 本课要用的旧手艺(熟练可跳过)
① 求根公式与判别式:aP2+bP+c=0⇒P=2a−b±b2−4ac;判别式缩到 0 ⇔ 两根相向而行、合并成一个。微例:P2−1000P+187500=0⇒P=21000±500——例题 3 原式。⚠️ 翻车点:b 本身带负号时 −b 的符号。
② 用 ln 解指数方程:ekt=A⇔kt=lnA;ln(ab)=blna。微例:e3k=2⇒k=3ln2——任务 1 原式。⚠️ 翻车点:ln(a+b) 没有展开式,别硬拆。
③ 本单元兵器速览:分离变量(第 2 课)、伯努利换元 v=y1−n(第 4 课)、相线图与 f′ 判别法(第 5-6 课)——四幕战役全靠它们,生疏哪个回哪课。
第一幕 · 马尔萨斯模型:增长率与人口成正比
最朴素的假设:人越多,生得越多——增长率正比于当前人口 P:
dtdP=kP(k>0)
第 2 课例题 1 的翻版,分离变量立得 P=P0ekt——指数增长。马尔萨斯 1798 年据此推出"末日论"(人口指数增长、粮食至多线性增长,饥荒是数学必然),这个模型因此也叫末日模型。
它短期内相当准(资源充足的细菌、疫情初期的感染数都走这条曲线),长期必然破产:指数增长意味着"每过固定时间翻一番",任何有限的地球都装不下。模型的失败之处恰恰标出了改进方向——资源上限,第三幕处理它。相线图一秒画完:唯一平衡点 P∗=0,f′(0)=k>0,不稳定——只要有一点人口,就一路奔向无穷。
第二幕 · 末日模型 + 收割:第一道生死线
给种群加上恒定速率的收割(捕捞、狩猎、移民流出),每单位时间拿走 h:
dtdP=kP−h
例题 1 · 定性 + 定量双杀
第 1 步 · 平衡点与稳定性
f(P)=kP−h=0 给出 P∗=kh;f′(P∗)=k>0——不稳定,山顶的球。相线:P>h/k 时 f>0 ↑(增长跑赢收割,奔向无穷);P<h/k 时 f<0 ↓(收割吃穿老本,坠向灭绝)。
P∗=h/k 是一条生死线:种群命运由初值站在线的哪一侧决定,而且这个平衡自身留不住——第 6 课"不稳定平衡是纸面状态"的活例。
第 2 步 · 精确解验证图景
方程是一阶线性(也可分离)。用积分因子 μ=e−kt:(e−ktP)′=−he−kt,积分得 e−ktP=khe−kt+C:
P(t)=kh+(P0−kh)ekt读解:偏差 (P0−kh) 被 ekt 指数放大(对照咖啡的 e−kt 指数遗忘——符号一换,天堂地狱)。P0 略低于生死线时,括号为负、放大到 P=0:有限时间灭绝。定性图景逐条兑现。
第三幕 · Logistic 模型:给增长装上天花板
修正马尔萨斯:环境有承载力 K(能养活的最大数量),越接近 K,增长越乏力。最简洁的实现——给 kP 乘一个"剩余空间比例"因子:
dtdP=rP(1−KP)
这就是 Logistic 方程。两个极端自检:P≪K 时括号 ≈1,退化为马尔萨斯(空旷时代随便长);P→K 时括号 →0,增长熄火(满员了)。
例题 2 · 三板斧齐下
第 1 步 · 定性分析(30 秒出结论)
f(P)=rP(1−P/K)=0:平衡点 P∗=0 与 P∗=K。求导 f′(P)=r−K2rP:f′(0)=r>0 不稳定;f′(K)=−r<0 稳定。
结论已经能写报告了:任何正初值的种群终将稳定在承载力 K。灭绝态 0 是送客的,碰不得一点生机。
第 2 步 · 精确解(第 4 课的排练正式登台)
这是伯努利方程(n=2):展开 P′−rP=−KrP2,令 v=1/P 走一遍第 4 课流程(那里解的 y′=y−y2 就是 r=K=1 的特例),得
P(t)=1+Ce−rtK,C=P0K−P0P0<K 时 C>0:标准 S 形曲线——起步慢(基数小)、中段快(火力全开)、封顶慢(空间挤)。且 t→∞ 时 P→K,与第 1 步的定性结论互相盖章。
第四幕 · Logistic + 收割:渔场崩溃之谜
压轴戏。给 Logistic 种群加恒定捕捞量 h:
dtdP=rP(1−KP)−h
例题 3 · 一个渔场的体检报告
某渔场 r=0.4/年,K=1000 吨,年捕捞 h=75 吨。这个方案可持续吗?安全边际有多大?
第 1 步 · 求平衡点(二次方程登场)
0.4P(1−1000P)−75=0,两边乘 0.41000=2500 并整理(🔧 工具箱 ①:求根公式):
P2−1000P+187500=0⇒P=21000±106−750000=21000±500两个平衡点:P∗=250 和 P∗=750。
第 2 步 · 判稳定性,读体检报告
f′(P)=0.4−0.0008P:f′(750)=−0.2<0 稳定;f′(250)=+0.2>0 不稳定。
报告:鱼群会稳定运行在 750 吨(可持续 ✓);但 250 吨是暗礁——一场赤潮、一次偷捕把存量砸破 250,系统就滑进 f<0 的深渊,再也回不来。收割把 Logistic 原本"任何正初值都得救"的宽厚,改成了"跌破门槛就灭绝"的冷酷。
第 3 步 · 临界捕捞量:崩溃为什么是“突然”的
调大 h,判别式 K2−r4hK 缩小——两个平衡点相向而行(稳定点降、暗礁升),在判别式归零处相撞湮灭(第 6 课练习 3 的剧情原封重演)。临界值:
hc=4rK=40.4×1000=100 吨/年h 越过 100 的一瞬间,平衡点消失、f(P)<0 处处成立:无论存量多雄厚,灭绝是唯一结局。最阴险的是崩溃前毫无渐变的预警——h=99 时系统看起来还稳稳地运行着(平衡点仍在),h=101 时大势已去。这正是历史上多个渔业(如纽芬兰鳕鱼场)"前一年还有鱼、后一年就没了"的数学机制。
✍️ 练习(顾问任务)
任务 1 · 细菌的马尔萨斯时刻
培养皿中细菌每 3 小时翻一番。写出模型并求:初始数量增至 10 倍需要多久?
✅ 参考答案
P=P0ekt,翻倍条件 e3k=2 给 k=3ln2。求 ekt=10:t=kln10=ln23ln10≈9.97 小时。指数增长的直觉冲击:翻一番要 3 小时,翻十倍只要不到 10 小时(因为 10≈23.32——十倍不过是三番多一点)。
任务 2 · Logistic 实战
池塘 r=1/年、K=100 条,投放 P0=10 条鱼。写出 P(t),并求鱼量达到 50 条(半程)的时间。
✅ 参考答案
C=P0K−P0=1090=9,故 P(t)=1+9e−t100。令 P=50:1+9e−t=2,e−t=91,t=ln9≈2.2 年。附赠一个 S 形的小秘密:P=K/2 恰是增长率 f(P) 的最大值点(对 f 求导可证)——半满时长势最猛。
任务 3 · 新渔场评估
渔场 r=0.6/年、K=2000 吨、拟定年捕捞 h=288 吨。求平衡点、稳定性、临界捕捞量,并给一句话结论。
✅ 参考答案
平衡方程两边乘 0.62000:P2−2000P+960000=0,判别式 4×106−3.84×106=160000,⋅=400:P∗=800(不稳定暗礁)与 P∗=1200(稳定运行位,f′=0.6−0.0006P,f′(1200)=−0.12<0 ✓)。临界值 hc=40.6×2000=300 吨/年。
结论:方案可行但逼近红线(288 对 300,安全边际 4%),暗礁高达 800 吨——建议下调捕捞量为系统留出扰动缓冲。
📝 随堂测验
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′=kP:指数增长,短期神器、长期破产;加收割后出现不稳定生死线 h/k
- Logistic P′=rP(1−P/K):0 不稳、K 稳,S 形曲线;P=K/2 处长势最猛
- Logistic + 收割:稳定运行位 + 不稳定暗礁并存;hc=rK/4 一过,两点湮灭、必然灭绝、没有渐变预警
- 定性分析主导全程,精确解只在需要"多快"时出手——工具各司其职
下一课换一个舞台但不换武功:混合问题——水箱、盐水、进进出出。你会发现药物代谢、房间通风、湖泊治污全是同一道题换了戏服,而主角又是积分因子。