跳转至

第1章：单物种连续种群模型（Continuous Population Models for Single Species）

§1 作者

本章是 J.D. Murray《Mathematical Biology I. An Introduction》(Third Edition, Springer 2002, IAM Vol. 17) 的开篇章节,作者为 J.D. Murray 本人,并非特邀作者综述。Murray 是牛津大学(及后续宾大、UCL、Oxford 数学研究所)长期从事应用数学与数学生物学研究的学者,本书是其在生物数学领域多年教学与科研的结晶。该章是全书"单种群动力学"四章(Ch. 1-4)的引入章节,扮演"建模哲学 + 工具入门"的双重角色——既给出 Malthus / Verhulst 等经典模型,又通过 Budworm 害虫爆发、Cheyne-Stokes 呼吸、血液生成等实例,展示微分方程、延迟微分方程(DDE)、偏微分方程(PDE)与非线性动力学(含分岔与混沌)如何服务于生物学建模。

§2 内容概述

本章围绕"单一物种连续时间种群动力学"展开,从最朴素的指数增长出发,逐步引入密度依赖、捕食、年龄结构等修正因素,展示了为什么"真实种群"几乎从不服从简单的 Malthus 律。

第 1.1 节建立 \(dN/dt = f(N)\) 的一般框架,讨论 Malthus 模型 \(dN/dt = (b-d)N\) 与 Verhulst 物流模型 \(dN/dt = rN(1 - N/K)\)。Murray 明确指出,物流模型不应被视为"宇宙普适定律",而是一个数学上便利的密度依赖调节隐喻,这一点借 Pearl 1925 年拟合美、法人口数据失败的史实得到警示。

第 1.2 节给出"云杉芽虫(Spruce Budworm)"作为双稳态的经典案例:三个正平衡点中 \(u_1\) 为"避难(Refuge)",\(u_3\) 为"爆发(Outbreak)",中间 \(u_2\) 不稳定,参数 \((r,q)\) 平面上出现尖点(cusp)区域。这段是本章核心——它把"灾难性跳变"从纯数学引申到视觉感知: Fisher 1967 那组从男人脸到女人的图片序列,57 名学生实验数据展示出单次感知滞后——与生物系统的可重复迟滞形成微妙对比。

第 1.3-1.4 节引入延迟微分方程 \(dN/dt = rN(t)[1 - N(t-T)/K]\),用 Nicholson 澳大利亚羊蝇实验数据说明延迟如何引入极限环周期解;线性化导出超越特征方程 \(\lambda = -e^{-\lambda T}\),得到 Hopf 分岔点 \(rT_c = \pi/2\),并对 \(rT\) 在 \(\pi/2\) 附近作扰动展开预测周期约为 \(4T\)。

第 1.5 节把延迟建模推向两个生理学应用: (a) Cheyne-Stokes 呼吸——用 Hill 函数描述 \(\text{CO}_2\) 对通气量的非线性响应,延迟来自脑干化学感受器,临界条件是通气梯度 \(V_0'\) 超过 \(\pi/(2\alpha x_0 T)\); (b) 血液生成调节—— Mackey-Glass 方程,调节参数 \(m\) 经历级联周期倍化→混沌的路径,与慢性白血病白细胞计数的实验观察惊人相似。

第 1.6 节把单种群模型用于可再生资源捕捞:Beddington-May 物流+恒定努力捕捞给出最大持续产量 \(Y_M = rK/4\),相应的稳定平衡 \(K/2\);并指出努力策略比恒定产量策略更稳健——后者在临界点附近恢复时间发散。Murray 在末尾引用了 Ludwig 关于鳕鱼、北大西洋渔业崩溃的告诫,严肃指出即使有科学知识,政治与短期经济仍可能阻碍合理管理。

第 1.7 节给出年龄结构建模: Von Foerster 方程 \(\partial n/\partial t + \partial n/\partial a = -\mu(a) n\) 加边界 \(n(t,0) = \int_0^\infty b(a) n(t,a) da\),通过相似解 \(n(t,a) = e^{\gamma t} r(a)\) 得到阈值参数 \(S = \int_0^\infty b(a) \exp[-\int_0^a \mu(s) ds] da\),\(S>1\) 增长,\(S<1\) 衰减。

§3 核心方程与概念

3.1 连续增长与线性稳定性

最简守恒律:

\[ \frac{dN}{dt} = \text{births} - \text{deaths} + \text{migration} = f(N). \tag{1.4} \]

Malthus 指数律: \(dN/dt = bN - dN\),\(N(t) = N_0 e^{(b-d) t}\)。\(b>d\) 增长,\(b<d\) 死亡。

Verhulst 物流律:

\[ \frac{dN}{dt} = rN\left(1 - \frac{N}{K}\right), \quad r, K > 0. \tag{1.2} \]

显式解:

\[ N(t) = \frac{N_0 K e^{rt}}{K + N_0(e^{rt} - 1)} \to K, \quad t \to \infty. \tag{1.3} \]
  • \(r\): 线性出生率与线性死亡率之差,代表"恢复时间尺度的倒数",\(1/r\) 是种群对扰动的特征响应时间;
  • \(K\): 承载容量,反映环境资源约束;
  • \(N = 0\) 不稳定 (\(f'(0) = r > 0\)); \(N = K\) 稳定 (\(f'(K) = -r < 0\))。

Murray 强调"物流只是隐喻"——Pearl 1925 年用其拟合美法人口都只拟合了部分曲线,后续预测大幅偏离。

3.2 Spruce Budworm 模型: 尖点突变

加入 sigmoid 捕食 \(p(N) = BN^2/(A^2 + N^2)\):

\[ \frac{dN}{dt} = r_B N\left(1 - \frac{N}{K_B}\right) - \frac{BN^2}{A^2 + N^2}. \tag{1.6} \]

无量纲化 \(u = N/A\), \(\tau = Bt/A\), \(r = Ar_B/B\), \(q = K_B/A\):

\[ \frac{du}{d\tau} = r u\left(1 - \frac{u}{q}\right) - \frac{u^2}{1 + u^2} = f(u; r, q). \tag{1.8} \]

平衡点 \(f(u; r, q)=0\) 等价于 \(r(1 - u/q) = u/(1 + u^2)\)。当 \(r,q\) 落入"折叠区"时,系统有 3 个正平衡 \((u_1, u_2, u_3)\): \(u_1\), \(u_3\) 稳定, \(u_2\) 不稳定。

折叠边界参数化(习题 1):

\[ r(a) = \frac{2a^3}{(1 + a^2)^2}, \quad q(a) = \frac{2a^3}{a^2 - 1}, \quad a \ge \sqrt{3}. \tag{1.9} \]
  • \(r\) 增加, 系统从 A→B→C→D 走, 在 C 点 \(u_1\) 跳到 \(u_3\) (向上跳跃);
  • \(r\) 减小, 在 B 点 \(u_3\) 跳回 \(u_1\) (向下跳跃);
  • 形成尖点突变, 这是该模型作为害虫管理理论工具的核心。

3.3 感知突变: 单次滞后

Fisher 1967 图片序列 (男人脸→女人) 实验: 57 名学生, 三轮播放, 记录每次的"突变点"。

  • 第一轮 (1→8): 突变多发生于 7 (延迟锁在原感知);
  • 第二轮 (8→1): 突变多发生于 4-5 (中点附近);
  • 第三轮 (1→8): 突变多发生于 5 附近。

关键差异: 视觉感知的"滞后"只发生一次, 而 Budworm 的迟滞在每次越过临界点都会重现。

3.4 延迟模型与 Hopf 分岔

延迟 Logistic 方程:

\[ \frac{dN}{dt} = rN(t)\left[1 - \frac{N(t-T)}{K}\right]. \tag{1.12} \]

无量纲化 \(N^* = N/K\), \(t^* = rt\), \(T^* = rT\):

\[ \frac{dN}{dt} = N(t)\left[1 - N(t - T)\right]. \tag{1.15} \]

在 \(N = 1\) 附近线性化 \(N = 1 + n\), 得:

\[ \frac{dn(t)}{dt} \approx -n(t - T). \tag{1.16} \]

设 \(n(t) = ce^{\lambda t}\), 特征方程:

\[ \lambda = -e^{-\lambda T}. \tag{1.17} \]

实虚部展开: \(\mu = -e^{-\mu T}\cos\omega T\), \(\omega = e^{-\mu T}\sin\omega T\)。 令 \(\mu = 0\), 得临界:

\[ rT_c = \pi/2, \quad \omega_c = 1. \tag{1.19} \]

周期估计: 在 \(T = \pi/2 + \varepsilon\) 附近扰动, 解 \(\lambda = \mu + i\omega\) 给出:

\[ N(t) \approx 1 + \text{Re}\left\{ c \exp\left[\frac{\varepsilon t}{1 + \pi^2/4}\right] \exp\left[it\left(1 - \frac{\varepsilon\pi/2}{1 + \pi^2/4}\right)\right] \right\}. \tag{1.23} \]
  • 周期 \(t_p \approx 2\pi/[1 - \varepsilon\pi/2(1 + \pi^2/4)] \approx 2\pi\) (无量纲), 量纲下 \(4T\);
  • 当 \(\varepsilon = 0.029\) (\(rT = 1.6\)): 公式预测 \(4.05T\), 数值 \(4.03T\) (一致);
  • 当 \(\varepsilon = 0.53\) (\(rT = 2.1\)): 公式预测 \(5.26T\), 数值 \(4.54T\) (需高阶展开)。

Lyapunov 函数 (van den Driessche-Zou 1998) 对一般线性 DDE \(dy/dt = a y(t) + b y(t-\tau)\):

\[ L[y(t)] = y^2(t) + |b| \int_{t-\tau}^{t} y^2(s) ds. \tag{1.25} \]

\(dL/dt \le 0\) 当 \(a < -|b|\), 此时零解全局稳定 (不仅是线性稳定)。

3.5 生理学延迟模型

Cheyne-Stokes 呼吸 — \(\text{CO}_2\) 控制:

  • 通气量: \(V = V_{\max} c^m(t-T) / [a^m + c^m(t-T)]\) (Hill 函数);
  • 动力学: \(dc/dt = p - bV c\);
  • 无量纲化得 \(x' = 1 - \alpha x V(x(t-T))\), \(V(x) = x^m/(1 + x^m)\);
  • 线性化超越方程: \(\lambda = -A - B e^{-\lambda T}\) 其中 \(A = \alpha V_0\), \(B = \alpha x_0 V_0'\);
  • 临界条件: \(\alpha x_0 V_0' T < [(\alpha V_0 T)^2 + s_1^2]^{1/2}\), \(s_1 \cot s_1 = -\alpha V_0 T\);
  • 小 \(AT\) 近似: \(V_0' < \pi/(2\alpha x_0 T)\)。

正常人类参数 (\(c_0 = 40\) mmHg, \(T = 0.25\) min): \(\alpha V_0 T = 0.0375\), 临界 \(V_0' = 7.44\) litre/min·mmHg, 实测周期 \(\approx 1\) min \(= 4T\)。

血液生成 (Mackey-Glass 1977):

\[ \frac{dc}{dt} = \frac{\lambda_0 a^m c^m(t-T)}{a^m + c^m(t-T)} - gc. \tag{1.43} \]

调节 Hill 系数 \(m\):

  • \(m=7\): 简单极限环;
  • \(m=10, T=6\) 天: 低幅周期解, 周期 \(\approx 20\) 天 (与图 1.15a 数值一致);
  • \(m=10, T=20\) 天: 混沌;
  • \(m=20\): 复杂混沌相图 (图 1.16e)。

血液生成模型展示了级联周期倍化通往混沌的完整 Feigenbaum 路径, 与下章离散模型的理论分析互相印证。

3.6 捕捞模型

Beddington-May 恒定努力模型:

\[ \frac{dN}{dt} = rN\left(1 - \frac{N}{K}\right) - EN. \tag{1.44} \]

新平衡 \(N_h(E) = K(1 - E/r)\) (\(E < r\))。最大持续产量:

\[ Y_M = \left. EN_h\right|_{E = r/2} = \frac{rK}{4}, \quad N_h\big|_{Y_M} = \frac{K}{2}. \tag{1.47} \]

恢复时间比:

\[ \frac{T_R(E)}{T_R(0)} = \frac{1}{1 - E/r}. \tag{1.50} \]
  • \(E = r/2\) (最优点): \(T_R = 2 T_R(0)\);
  • \(E \to r\): \(T_R \to \infty\) (系统变脆)。

Brauer-Sanchez 恒定产量模型:

\[ \frac{dN}{dt} = rN(1 - N/K) - Y_0. \tag{1.52} \]

两个正平衡, 其中稳定的 \(N_2(Y_0) = (K/2)[1 + (1 - 4Y_0/(Kr))^{1/2}]\)。 当 \(Y_0 \to Y_M\) 时 \(T_R \to \infty\), 因此恒定产量策略对扰动更敏感。

3.7 年龄结构与 Von Foerster 方程

\[ \frac{\partial n}{\partial t} + \frac{\partial n}{\partial a} = -\mu(a) n, \quad n(t, 0) = \int_0^\infty b(a) n(t, a) da. \tag{1.54, 1.56} \]

相似解 \(n(t, a) = e^{\gamma t} r(a)\) 代入:

\[ r(a) = r(0) \exp\left[-\gamma a - \int_0^a \mu(s) ds\right]. \tag{1.63} \]

边界条件化为:

\[ \phi(\gamma) = \int_0^\infty b(a) \exp\left[-\gamma a - \int_0^a \mu(s) ds\right] da = 1. \tag{1.64} \]

阈值 \(S = \phi(0)\): \(S > 1\) 增长, \(S < 1\) 衰减。这是 Lotka 的"净生殖率"推广。

§4 关键结论

  1. Malthus 律仅在无约束环境下短期有效——17 至 21 世纪世界人口数据大致符合指数增长, 但承载容量 \(K\) 与密度依赖机制必然介入。物流模型不是普适定律, 而是密度依赖的"原型"。

  2. Spruce Budworm 是双稳态 + 突变的生物物理原型—— 临界尖点 \((r_c, q_c)\) 在 \((\sqrt{3}/8, 2\sqrt{3})\) 附近 (从参数化方程), 三个正平衡 \(u_1 < u_2 < u_3\) 解释"避难"与"爆发"两种状态。害虫管理的核心目标是让 \((r, q)\) 离开折叠区——可减少 \(K_B\) (喷洒)、降 \(r_B\) (不育技术)、增 \(A\) (招引天敌)。

  3. 延迟引入极限环—— 单标量 ODE \(dN/dt = f(N)\) 不可能有限环 (能量积分), 引入 \(T > 0\) 即可。Hopf 分岔点 \(rT_c = \pi/2\) 在 \(N = K\) 稳定态产生振荡, 周期 \(\approx 4T\)。Nicholson 羊蝇实验数据 (35-40 天) 与该预测高度一致。

  4. Lyapunov 方法给出全局稳定域——对线性 DDE \(dy/dt = a y(t) + b y(t-\tau)\), \(L[y] = y^2 + |b|\int_{t-\tau}^t y^2(s) ds\) 是 Lyapunov 函数, 全局稳定条件 \(a < -|b|\) 比精确谱分析得到的"线性稳定"略小, 但同型。

  5. Cheyne-Stokes 呼吸源于反馈时延—— 通气响应曲线越陡 (Hill 系数 \(m\) 越大) 越易不稳定。Murray 给出的临界 \(V_0' = 7.44\) L/min·mmHg 是与生理实验一致的"硬数字", 而 Fowler-Kalamangalam (2000) 的替代解释 (心输出量降低) 提示该模型并不唯一。

  6. Mackey-Glass 方程再现周期倍化→混沌—— 在 \(m\) 增加时, 极限环经过分岔变成双环、四环、… 进入混沌相, 与慢性白血病的白细胞计数 (准周期、72 天) 形态相似 (图 1.15b vs c)。Murray 强调"形态相似≠机制相同"——建模只提供假设, 必须有医学协作。

  7. 恒定努力捕捞比恒定产量捕捞更稳健—— \(T_R\) 在 \(Y_0 \to Y_M\) 时发散, 而努力策略下 \(T_R(E)\) 始终有限。这对渔业政策有直接含义: 配额制比总产量上限制更安全 (尽管两者都依赖准确的参数估计)。

  8. 年龄结构模型有显式阈值—— Lotka 阈值 \(S = \int_0^\infty b(a) \exp[-\int_0^a \mu(s) ds] da\), \(S > 1\) 增长。这是单纯出生率/死亡率数据的长期预测, 不需要解整个 PDE。

§5 挑战和开放性问题

  1. Spruce Budworm 的空间维度被本章完全忽略。Murray 明确说"第 II 卷第 2 章再讨论扩散"。问题: 当芽虫扩散尺度与森林斑块尺度可比时, 避难状态是否能被邻域爆发态"传染"? 实际林业管理需要空间显式模型。

  2. Lyapunov 函数的"过保守"。van den Driessche-Zou 1998 的 \(L[y]\) 给出 \(a < -|b|\), 但谱分析给出的精确稳定域更大。能否构造"恰好"刻画边界的 Lyapunov 函数, 仍是 DDE 稳定性理论的开放问题。

  3. Cheyne-Stokes 呼吸的延迟参数 \(T\) 难以测量。Murray 自己承认实测周期 (2-3 倍延迟) 与模型预测 (4 倍) 差距过大。Fowler-Kalamangalam 2000 提出的"心输出量降低"机制与 Murray 模型互斥而非互补, 哪个对? 这是临床医学仍有争议的开放问题。

  4. Mackey-Glass 混沌的临床意义。Murray 承认图 1.15b 与 c 的相似不应推断"延迟是白血病的主因"。问题是: 在已知有反馈时延的生理系统中, 哪些是真的混沌? 真正的"白细胞混沌"实验需要高时间分辨率的长期数据, 这在临床上很难获得。

  5. 捕捞模型的"开环"假设。Beddington-May 不含经济反馈, 但 Murray 引用 Clark 1976/1985/1990 指出真实渔业管理需要经济反馈 (价格、补贴、监管) 才能稳定。问题: 在"自上而下"的渔业配额政策中, 反馈时延本身就是另一个延迟——是 Hopf 分岔源。

  6. Von Foerster 相似解的初始条件不匹配。Murray 在 §1.7 末尾承认 \(n(t,a) = e^{\gamma t} r(a)\) 只能刻画"长时间后"行为, 初始瞬态需要数值或迭代 (1.61) 处理。短时间人口预测 (例如 5-10 年) 的可靠性受此限制。

  7. 延迟高阶展开何时失效。Murray 给出 \(rT = 2.1\) (\(\varepsilon = 0.53\)) 时一阶展开预测 \(5.26T\), 数值 \(4.54T\)。二阶展开能修复, 但需要重新作系统的双时间尺度匹配, 是一类还未在数学生物学界完全自动化的"形式摄动"工作。

  8. 季节性 (周期性参数)。所有本章模型都假设参数恒定, 但 Budworm 真实出生率是季节驱动, 通气量受昼夜节律调制。引入时变参数 \(r(t) = r_0(1 + \varepsilon \cos\omega t)\) 会让分岔分析进入 Floquet 理论领域, 这是该章未触及的方向。

§6 个人反思与批判性分析

对 Murray 建模哲学的评价: Murray 的最大贡献是把"建模的诚实"放在"解的精巧"之前——他在 1.1 节就引用 Pearl 1925 的失败例子, 1.5 节末尾又强调"形态相似≠机制相同"。这种反"过度数学化"的姿态在 1990-2000 年代的应用数学界并不常见, 是该书能成为经典的核心原因。

本章的"叙事跳跃": 1.1 物流 → 1.2 害虫尖点 → 1.3-1.4 延迟 Hopf → 1.5 生理学混沌 → 1.6 渔业 → 1.7 年龄结构, 七节之间没有递进关系, 几乎是"建模工具箱巡礼"。这既是优点 (展示多样性), 也是缺点 (读者很难形成"主线索")。如果重写, 建议在第 1.1 节末尾就预告"非线性反馈→分岔→混沌"这条贯穿全章的逻辑链。

对 Budworm 模型的具体批评: 该模型用 sigmoid 捕食 \(BN^2/(A^2 + N^2)\) 是 ad-hoc 的; 更现代的捕食函数 (Beddington-DeAngelis, Crowley-Martin, ratio-dependent) 给出不同的分岔结构。Ludwig 1978 的原始论文只给出"定性讨论", 后续工作 (如 Dwyer 1994) 证明加入空间维度后, 害虫循环的周期由扩散而非延迟决定, 这是 Murray 没有引用的。

与自身研究的相关性: 作为计算生物力学研究者, 我最感兴趣的是 1.5 节的 Hill 函数——它在血管平滑肌钙动力学、动脉壁重塑、应力-生长律中反复出现。Murray 在 1.5 节末提到"Keener & Sneyd 1998", 而该书正是我做血管生物力学时的参考。Murray 强调的"延迟如何引起 Hopf 分岔"在心血管系统中是实在的——比如动脉压力感受器反射弧 (\(\approx 1-2\) 秒延迟) 在某些病理条件下引发 Mayer 波 (周期约 10 秒, \(\approx 5-10 T\)), 与本章 1.4 节的预测框架完全同构。

未在章末习题中得到引导的关键问题: 习题 1-8 中, 习题 2 (Malthus 周期解的反例) 是哲学性提示, 习题 5 (须鲸延迟模型) 直接对应国际捕鲸委员会现实问题, 习题 6/7 是 Cheyne-Stokes 的进一步推导。但没有习题引导读者思考空间效应——这是第 II 卷 (MurrayII) 弥补的方向。

关于 Cheyne-Stokes 呼吸, 我有一个具体疑问: Murray 把 \(\text{CO}_2\) 当作唯一调控变量, 但现代呼吸生理学公认 \(\text{O}_2\) 也参与 (尤其在高原/心衰患者中)。如果把 \(\text{CO}_2\) 和 \(\text{O}_2\) 同时建模, 会得到 2 维 DDE 系统, 其分岔结构会从 Hopf 变成 Bautin (即存在弱振幅与强振幅极限环的共存) 或双 Hopf。这是 Murray 1993 年第一版未深入的方向, 2002 年第三版也没有大幅更新。

关于血液生成模型的统计推断: Mackey-Glass (1.43) 的参数 (\(\lambda_0, a, m, g, T\)) 都从临床数据拟合, 但患者个体差异大。Murray 没讨论参数不确定性如何影响混沌预测——这是一个"鲁棒混沌"概念, 在 Hassell 1999 的复杂芽虫模型中才出现。我倾向于把这一类"参数识别 + 不确定性传播"问题作为后续与计算生物学家合作的接口。

§7 重要参考文献

[X1] Malthus, T.R. (1798). An Essay on the Principle of Population. J. Johnson, London. (首次系统阐述指数增长, 1798 年匿名出版, 后续引发 Darwin 等进化论先驱的思考。)

[X2] Verhulst, P.F. (1838, 1845). "Correspondance mathématique et physique" (十年间两篇通讯, 提出物流方程 \(dN/dt = rN(1 - N/K)\) 的原始形式。)

[X3] Pearl, R. (1925). The Biology of Population Growth. Williams and Wilkins, Baltimore. (用物流方程拟合美、法人口数据, 后续预测失败, 警示"参数少≠模型好"。)

[X4] Ludwig, D., Jones, D.D., Holling, C.S. (1978). "Qualitative analysis of insect outbreak systems: the spruce budworm and forest." Journal of Animal Ecology 47: 315-332. (本章 Budworm 模型的原始来源, 把害虫循环放到生态系统层面。)

[X5] Fisher, R.A. (1967). The Theory of Confidence Intervals and the Use of Maximum Likelihood Estimates. (Murray 引用其图片序列作为感知突变的视觉证据, 后续 Zeeman 1972 用尖点突变理论建模。)

[X6] Zeeman, E.C. (1982). "Catastrophe theory in psychology." Behavioral Science 27: 39-49. (把 Budworm 尖点结构与感知突变连接, 区分"可重现迟滞"与"单次滞后"。)

[X7] Nicholson, A.J. (1957). "The self-adjustment of populations to change." Cold Spring Harbor Symposia on Quantitative Biology 22: 153-173. (澳大利亚羊蝇 2 年实验, 35-40 天周期, 验证延迟 Logistic 1.12 的 \(4T\) 周期预测。)

[X8] Mackey, M.C., Glass, L. (1977). "Oscillation and chaos in physiological control systems." Science 197: 287-289. (本章 1.5 血液生成 + Cheyne-Stokes 的原始模型来源, 引入"动态疾病"概念。)

[X9] Glass, L., Mackey, M.C. (1988). From Clocks to Chaos: The Rhythms of Life. Princeton University Press. (生理节律与混沌的科普专著, Murray 在 1.5 节中作为后续阅读推荐。)

[X10] Beddington, J.R., May, R.M. (1977). "Harvesting natural populations in a randomly fluctuating environment." Science 197: 463-465. (本章 1.6 节恒定努力捕捞模型的原始论文, 给出"最大持续产量 \(Y_M = rK/4\)"的结论。)

[X11] van den Driessche, P., Zou, X. (1998). "Global attractivity in delayed Hopfield neural network models." SIAM Journal on Applied Mathematics 58: 1878-1890. (给出 1.4 节 Lyapunov 函数 1.25 的构造法, 适用一切线性 DDE。)

[X12] Fowler, A.C., Kalamangalam, G.P. (2000). "A mathematical analysis of the Grodins model of respiratory control." Mathematical Medicine and Biology 17: 67-82. (对 Cheyne-Stokes 提出"心输出量降低"替代机制, 与 Murray 模型形成对比。)

[X13] Ludwig, D. (1993, 1994, 1995, 1996a, 1996b); Ludwig, D., Hilborn, R., Walters, C. (1993). "Uncertainty, resource exploitation, and conservation: lessons from history." Science 260: 17-36. (Murray 在 1.6 节末尾引用的渔业崩溃系列论文, 强调"科学知识≠政策实施"。)

[X14] Leslie, P.H. (1945). "On the use of matrices in certain population mathematics." Biometrika 33: 183-212. (本章 1.7 提及的年龄结构矩阵模型开创者, 与 Von Foerster 偏微分方程互为补充。)

[X15] von Foerster, H. (1959). "Some remarks on changing populations." In The Kinetics of Cellular Proliferation, Grune and Stratton. (本章 1.7 节 PDE 年龄模型的原始文献, 后来广泛用于细胞动力学。)

[X16] Keener, J.P., Sneyd, J. (1998). Mathematical Physiology. Springer, Interdisciplinary Applied Mathematics Vol. 8. (Murray 在 1.5 节提到的"生理学数学"参考书, 与本书构成"动力学"与"生理"双线。)