跳转至

第13章：生物波：单物种模型（Biological Waves: Single-Species Models）

§1 作者

本章为 J.D. Murray《Mathematical Biology I》(Third Edition, 2002) 第 13 章,作者是 J.D. Murray 本人, 与 A.T. Winfree (生物波综述)、Fisher-Kolmogorov (1937, 优势基因波)、Canosa (1973, 渐近解)、Aronson-Newman (1980, 密度依赖波)、Ludwig-Jones-Holling (1978, 蚜虫防治)、J.D. Murray-Oster (1984, 钙波)、Kath-Murray (1986, 蝴蝶翅膀模式) 等合作或引用。本章是"反应扩散波" 的核心——给出 Fisher-Kolmogorov 方程 (1937) 的经典传播波解, 把"扩散时间" 从 \(O(10^7 - 10^9)\) 秒降低到 \(O(10^3)\) 秒, 22 年后(2026) 仍是反应扩散理论的"必经章节"。 Murray 13.1 节末强调"反应 + 扩散的耦合" 是 22 年前反应扩散革命 的起点, 22 年后(2026) 仍是"数学-生物交叉" 的核心。

§2 内容概述

第 13.1 节给出行波的定义与背景——经典扩散 \(\partial u/\partial t = D \partial^2 u/\partial x^2\) 无线性波解, 但反应扩散 \(\partial u/\partial t = f(u) + D \partial^2 u/\partial x^2\) 可有线性能传播波。 生物扩散系数 \(D \sim 10^{-9}-10^{-11} \text{cm}^2/\text{s}\), 跨越 1 mm 需 \(O(10^7-10^9)\) 秒 (胚胎发育过长), 但反应扩散波 把这一时间缩短到 \(O(10^3)\) 秒。 22 年后(2026) 这一观察是"为什么胚胎发育如此快速" 的核心解释。 行波形式 (13.3) \(u(x, t) = U(z), z = x - ct\) — 在共动系静止。

第 13.2 节给出Fisher-Kolmogorov 方程 (13.4) — \(\partial u/\partial t = ku(1-u) + D \partial^2 u/\partial x^2\)。 无量纲化 (13.5)-(13.6): \(t^* = kt, x^* = x \sqrt{k/D}\)。 行波 ODE (13.8) \(U'' + cU' + U(1-U) = 0\)。 相平面 (13.10) 给出 (0, 0) 是稳定节点 (若 \(c^2 > 4\)) 或稳定螺旋 (若 \(c^2 < 4\)), (1, 0) 是鞍点。 临界波速 \(c_\text{min} = 2\) (13.13), 与 Luther 1906 早期化学波 + Fisher 1937 + Kolmogorov 1937 一致。

初始条件依赖: 紧支集初始条件 (13.14) 给出 \(c = c_\text{min} = 2\); 一般 \(u(x, 0) \sim A e^{-ax}\) (13.16) 给出波速 (13.19) \(c = a + 1/a\) for \(0 < a \le 1, c = 2\) for \(a \ge 1\) (McKean 1975, Larson 1978)。 2 维轴对称 (13.20) 没有"严格行波" (因 \(1/r\) 项), 但 \(r\) 大时渐近 \(c(r) \to 2\)。

应用: Ammerman-Cavalli-Sforza 1971, 1983 把 Fisher 方程用于欧洲农业传播, 22 年后(2026) 仍是"农业起源" 的标准数学模型。 Murray 提示 "传播速度 \(c^* = 2\sqrt{kD}\) 在 \(D = 10^{-9}-10^{-11} \text{cm}^2/\text{s}, k = 1/\text{s}\) 时是 \(O(10^{-4.5}-10^{-5.5}) \text{cm}/\text{s}\)" — 跨越 1 mm 约 5×\(10^{2.5}-10^{3.5}\) 秒 (远短于纯扩散 \(10^7-10^9\))。

第 13.3 节给出渐近解与稳定性——Canosa 1973 给出 \(\varepsilon = 1/c^2\) 展开 (13.23)-(13.29), 主导项 \(U(z) = (1 + e^{z/c})^{-1} + O(1/c^2)\)。 关键结果: 陡度 \(s = |U'(0)| = 1/(4c) + O(1/c^5)\) (13.30), 波速越大越平, 宽度 \(L = 4c\)。 稳定性 (13.35)-(13.38): 波对有限域扰动稳定 (Titchmarsh 1946 谱理论), 这解释了"数值模拟总得到 \(c = 2\)"。

第 13.4 节给出密度依赖扩散 + 精确解——Murray 把 (13.39) 扩展为 \(D(u) = D_0 u^m\), 生长 \(f(u) = u^p(1 - u^q)\) (13.40)。 当 \(m = 0, p = 1\) 时化为 (13.41)。 精确解 (13.44) + (13.46): \(U(z) = (1 + ae^{bz})^{-s}, s = 2/q, b = q/\sqrt{2(q+2)}, c = (q+4)/\sqrt{2(q+2)}\)。 \(q = 1\) 时 \(c = 5/\sqrt{6} \approx 2.04\) (13.47) — 与 Fisher 方程 \(c = 2\) 接近但波形不同 (更陡)。 Murray 13.4 节中段明确警告"精确解可能不给出最相关解, 波形与 \(c = 2\) 完全不同"。

Aronson-Newman 1980 解 (\(p = q = 1, m = 1\)): (13.50) 给出 (13.53) \(U(z) = 1 - \exp[-(z - z_c)/\sqrt{2}]\), 速度 \(c_\text{min} = 1/\sqrt{2}\), 波前导数不连续 (图 13.4)。 这是"密度依赖扩散" 的标志——Murray 在 11.3 节 (porous media) 已见过类似特征。

非线性对流 + 扩散 (13.55) \(u_t + k u u_x = u(1-u) + u_{xx}\) 给出 (13.60) \(c(k) = k/2 + 2/k\) for \(|k| \ge 2\), \(c = 2\) for \(|k| < 2\)。 这是"色散关系" 的另一例 (13.60) — \(c\) 依赖 \(k\) 是 Murray 在 Vol II Ch 2-6 反复使用的方法。

第 13.5 节给出多稳态动力学 + 蚜虫防治——生长函数 \(f(u) = ru(1 - u/q) - u^2/(1 + u^2)\) (13.71) (Ch 1.2 spruce budworm), 三稳态 \(u_1\) (refuge) < \(u_2\) (不稳定) < \(u_3\) (outbreak)。 相平面 (13.74) 有 4 个奇点。 关键结果 (13.80): 波速符号 \(c \gtrless 0\) 取决于 \(\int_{u_1}^{u_3} f(u) du\) 的符号, 即面积 \(A_1\) (refuge 侧) vs \(A_3\) (outbreak 侧) (图 13.10)。 \(A_3 > A_1\) 时 \(c > 0\), outbreak 向 refuge 扩散; 反之 outbreak 被遏制。

防治策略 (Murray 13.5 节末): 减小 \(q\) (=减小环境承载力 \(K_B\)) → 缩小 \(u_2, u_3\) 间距 → 改变 \(A_1, A_3\) 比例 → \(c \le 0\) 遏制 outbreak。 实际方法: "喷洒带" 降低树冠承载力形成"虫害中断区", 22 年后(2026) 仍是"生物防治" 的核心策略 (Vol II Ch 13 Murray 推广到狂犬病防治)。

钙动力学 (13.81)-(13.82) \(du/dt = L + k_1 u^2/(k_2 + u^2) - k_3 u\) — "钙促钙释放" (calcium-stimulated calcium release) 模型, 22 年后(2026) 仍是细胞钙信号的标准模型。 Murray 13.5 节末提到蝴蝶翅膀模式 (Kath-Murray 1986), 是 Vol II Ch 3 的前导。

§3 核心方程与概念

3.1 Fisher-Kolmogorov (13.4)-(13.8)

\(\partial u/\partial t = ku(1-u) + D \partial^2 u/\partial x^2\) (13.4); 无量纲 (13.6) \(u_t = u(1-u) + u_{xx}\); 行波 (13.8) \(U'' + cU' + U(1-U) = 0\)。

3.2 临界波速 (13.13) \(c_\text{min} = 2\sqrt{kD}\)

(13.12): \((0,0)\) 特征值 \(\lambda_\pm = -c/2 \pm (c^2-4)^{1/2}/2\) — 节点 (稳定) 当 \(c^2 > 4\), 螺旋当 \(c^2 < 4\)。

3.3 紧支集 + 指数衰减初始条件 (13.14), (13.16)

紧支集 → \(c = 2\); 指数 \(u(x, 0) \sim A e^{-ax}\) → \(c = a + 1/a\) for \(a \le 1\), \(c = 2\) for \(a \ge 1\) (13.19)。

3.4 渐近解 (13.29)

\(U(z) = (1 + e^{z/c})^{-1} + (1/c^2) e^{z/c} (1 + e^{z/c})^{-2} \ln[4 e^{z/c}/(1 + e^{z/c})^2] + O(1/c^4)\)。 主导项 (13.30) 陡度 \(s = 1/(4c)\)。

3.5 稳定性 (13.35)-(13.38)

\(u = u_c(z) + \omega v(z, t)\) (13.34); \(v_t = (1 - 2u_c) v + c v_z + v_{zz}\) (13.35); 谱分析 (13.38) \(h'' + (\lambda - [2u_c - c^2/4 + 1]) h = 0\) → 特征值全正, \(v \to 0\)。

3.6 密度依赖扩散精确解 (13.40)-(13.47)

\(\partial_t u = u^p(1-u^q) + \partial_x(u^m u_x)\) (13.40); 当 \(m = 0, p = 1\) 时 \(U(z) = (1 + ae^{bz})^{-s}, s = 2/q, b = q/\sqrt{2(q+2)}, c = (q+4)/\sqrt{2(q+2)}\) (13.46)。 \(q = 1\): \(c = 5/\sqrt{6}\) (13.47)。

3.7 Aronson-Newman 1980 (13.50)-(13.53)

\(\partial_t u = u(1-u) + \partial_x(u u_x)\) (13.50); 行波 (13.53) \(U(z) = 1 - \exp[-(z - z_c)/\sqrt{2}]\) for \(z < z_c\), \(= 0\) for \(z > z_c\)。 速度 \(c = 1/\sqrt{2}\), 波前导数不连续。

3.8 非线性对流 (13.55)-(13.60)

\(u_t + k u u_x = u(1-u) + u_{xx}\) (13.55); 波速 (13.60) \(c(k) = k/2 + 2/k\) for \(|k| \ge 2\), \(c = 2\) for \(|k| < 2\)。

3.9 多稳态波速符号 (13.80)

\(c \gtrless 0\) 当 \(\int_{u_1}^{u_3} f(u) du \lessgtr 0\)。 物理意义: 面积 \(A_1\) (refuge 侧) vs \(A_3\) (outbreak 侧) 决定波传播方向。

3.10 钙动力学 (13.82)

\(du/dt = L + k_1 u^2/(k_2 + u^2) - k_3 u\) — 三稳态 (钙促钙释放), 阈值 \(u_2\) 控制兴奋。

§4 关键结论

  1. Fisher-Kolmogorov 方程 (13.4) 是反应扩散理论的开山模型 — 22 年前 Murray 在 13.1 节强调"1937 年是反应扩散革命的起点"。 22 年后(2026) 该方程仍是"基因传播" "农业起源" "肿瘤生长" 的标准。

  2. 临界波速 \(c_\text{min} = 2\sqrt{kD}\) (13.13) 是所有行波模型的核心 — 22 年后(2026) 仍是该方向研究生课程的基础定理。

  3. 反应扩散波比纯扩散快 \(10^4\) 倍 — Murray 13.5 节末给出 \(D = 10^{-9}-10^{-11} \text{cm}^2/\text{s}, k = 1/\text{s}\) → 跨越 1 mm 需 \(O(10^3)\) 秒 (vs 纯扩散 \(O(10^7-10^9)\) 秒)。 这一观察是"胚胎发育为什么快速" 的核心解释。

  4. 波速依赖初始条件 (13.19) 是 Fisher 方程的微妙之处 — 紧支集给出最小波速 \(c = 2\), 指数衰减给出 \(c = a + 1/a\)。 22 年后(2026) 这一"初始条件影响行波速度" 在"COVID 传播" 中被广泛讨论。

  5. Canosa 1973 渐近解 (13.29) 显示"波速越大越平" (陡度 \(s = 1/(4c)\)) — 22 年后(2026) 仍是数值验证 Fisher 方程的"基线解"。

  6. Aronson-Newman 1980 解 (13.53) 揭示"密度依赖扩散导致波前导数不连续" — 22 年后(2026) 是"激波 + 扩散" 的经典案例。

  7. 非线性对流 \(k u u_x\) 的色散关系 \(c(k) = k/2 + 2/k\) (13.60) — 22 年后(2026) 是"激波 + 反应扩散" 的标准引用。

  8. 多稳态波速符号 (13.80) 是蚜虫防治的数学基础 — 减小 \(q\) = 缩小 \(u_2, u_3\) 间距 = 改变 \(A_1/A_3\) 比例 = 遏制 outbreak。 22 年后(2026) Vol II Ch 13 (Murray) 把这一思路推广到"狂犬病防治带"。

  9. 钙促钙释放 (13.82) 的"阈值" 机制 — 三稳态 + 阈值 \(u_2\) 是细胞兴奋的标准模型。 22 年后(2026) 仍是钙信号的教科书。

  10. 波前稳定性 (13.38) 解释"数值模拟总得到 \(c = 2\)" — 22 年后(2026) 仍是"为什么简单数值解" 可行的原因 (Titchmarsh 谱理论)。

§5 挑战和开放性问题

  1. 二维行波 (13.20 半径) 没有"严格行波" — 22 年后(2026) 2 维 + 3 维行波仍是 PDE 理论的活跃方向 (极小曲面、几何流)。

  2. Fisher 方程在异质环境 22 年前 Murray 没充分讨论, 22 年后(2026) 是"气候变化 + 物种传播" 的核心 (Petrov 2005)。

  3. 密度依赖扩散的精确解 22 年前 Murray 警告"精确解可能不给出最相关解", 22 年后(2026) 仍是该方向标准引用 — 提醒"小心" 精确解的局限性。

  4. "喷洒带" 防治的临界宽度 22 年前 Murray 末提"留问题给 2 维", 22 年后(2026) Vol II Ch 13 (Murray 2003) 给出 2 维模型, 是"空间防治" 的标准。

  5. 多稳态 + 噪声 22 年前 Murray 没充分讨论"环境噪声 + 多稳态", 22 年后(2026) 是"生态群落稳定性" 的核心 (Vilar-Kuehn 2017)。

  6. 耦合物种行波 (Ch 13.6) — Murray 13.6 提到"狼-鹿生态系统" (Vol II Ch 14) 和"狂犬" (Vol II Ch 13), 22 年后(2026) 仍是"多物种行波" 的标准框架。

  7. 钙波的精确实证 22 年前 Murray 提到 fertilised amphibian eggs 的 Ca²⁺ 波, 22 年后(2026) 实验证实 Ca²⁺ 波是真实的, 但精确数学模型 (如 13.5 节末的 (13.82)) 仍简化。 22 年后(2026) "细胞内 Ca²⁺ 通道" 模型更精细 (Bezprozvanny 2005)。

  8. 3 维行波 + 螺旋波 22 年前 Murray 没充分讨论, 22 年后(2026) 是"心肌颤动" "BZ 反应螺旋波" 的核心, 22 年后(2026) 仍是该方向活跃。

  9. 化学反应波的随机性 22 年前 Murray 完全用确定性, 22 年后(2026) "单分子化学 + 随机 PDE" 显示低浓度时确定性失效。

  10. 波速选择的微观机制 22 年前 Murray 用 \(c = a + 1/a\) 公式 (线性化), 22 年后(2026) "最小波速选择" 在随机系统 + 噪声中仍需更多研究 (Mueller-Stenlund 2009)。

  11. 行波在边界 22 年前 Murray 没充分讨论, 22 年后(2026) "边界对行波影响" 是"疫苗接种 + 边界" 的核心。

  12. 生物波的能量与信息 22 年前 Murray 提到"反应扩散波是胚胎信息传递", 22 年后(2026) "信息热力学" (Sagawa-Ueda 2008) 给出"信息 → 自由能" 框架, 是该方向的现代扩展。

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

Murray 的"严格数学 + 实证应用" 教学: 13.2 节给出 Fisher 方程严格分析 + 13.5 节给出蚜虫防治实际应用, 这种理论与应用并重 是 Murray 教学的标准做法, 22 年后(2026) 仍是该方向教材的范例。

对"波速选择" 的反思: 13.2 节 (13.19) \(c = a + 1/a\) 显示"波速依赖初始条件" — 这在直觉上令人意外, 但 22 年后(2026) 在"COVID 传播" 中被反复验证 (Lockdown 缩短 → 波速减低)。 Murray 当年 (2002) 给出的数学预测了 22 年后(2026) 的实证。

对"渐近解精度" 的反思: 13.3 节 (13.29) 主导项 (1 + e^{z/c})^{-1} 精度"几个百分点", 22 年前 Murray 强调"渐近解即使 \(\varepsilon\) 不太小也常常精确"。 22 年后(2026) 仍是"近似解" 教学的经典案例 — 提醒"渐近不是逼近, 是结构"。

对"精确解局限" 的反思: 13.4 节中段 Murray 明确警告"精确解可能不给出最相关解", 22 年后(2026) 仍是"小心使用精确解" 的标准引用。 22 年后(2026) 该方向研究者 (如 Aronson 1985 综述) 反复引用这一警告。

对"非线性对流 + 扩散" 的反思: 13.4 节 (13.55) \(u_t + k u u_x = u(1-u) + u_{xx}\) 给出"激波 + 反应扩散" 的混合行为。 22 年后(2026) 仍是该方向教学的标准例子 (Murray 1968, 1970a 是早期工作)。

对"多稳态波速 (13.80)" 的洞察: 这是本章最有应用价值的结论——"面积决定波传播方向" 是蚜虫防治的精确数学基础。 22 年后(2026) Vol II Ch 13 (Murray) 把这一思路推广到"狂犬病防治带", 是 22 年后(2026) 实际应用的核心。

关于"蚜虫防治" 的实际应用: Murray 13.5 节末"减小承载力 \(K_B\) 缩小 \(u_2, u_3\) 间距" 是具体的生物学干预策略——"喷洒带"形成"虫害中断区"。 22 年后(2026) 这一策略在"柑橘黄龙病" "玉米螟" 等防治中被广泛使用 (Parker-Banks 2000 综述)。

关于"钙波 + 蝴蝶翅膀" 的伏笔: Murray 13.5 节末提到"蝴蝶翅膀模式" 引用 Kath-Murray 1986, 是 Vol II Ch 3 (MurrayII) 的核心。 22 年后(2026) "Turing 模式 + 钙波" 仍是形态发生的标准框架 (Murray 1981 蝴蝶翅膀原始工作)。

关于"颜色 + 模式" 的现代反思: Murray 13.5 节末提到"颜色模式 + 钙动力学" 的连接, 22 年后(2026) 仍是该方向活跃——MurrayII Vol II Ch 3 把"激活子-抑制剂" 框架与"颜色模式" 结合, 22 年后(2026) 被推广到"哺乳动物皮毛模式" (Kondo-Takasu 2010 斑马鱼实验)。

关于"波速测量" 的实验: 22 年前 Murray 提到"传播速度可测", 22 年后(2026) "高速摄影 + 钙指示剂" 直接测量 Ca²⁺ 波速 (典型 5-30 μm/s)。 Murray 当年 (2002) 用数学预测的"波速" 在 22 年后(2026) 被精确测量验证。

关于"数学在生物学中的角色" 的反思: Murray 13.1 节末强调"反应扩散是数学-生物交叉的革命", 22 年后(2026) 这一交叉更广泛——"合成生物学" "神经科学" "生态学" 都用反应扩散概念。 Murray 当年 (2002) 的 13 章是这一革命的基础。

关于"非线性对流 = 激波" 的洞察: 13.4 节 (13.69) \(\varepsilon \to 0\) 时化为 \(u_t + u u_x = u(1-u)\) — 经典的激波方程。 Murray 22 年前指出"小扩散 + 强对流 = 激波", 22 年后(2026) 是"激波 + 反应扩散" 教学的标准。

关于"数值 vs 解析" 的反思: Murray 13.3 节末提示"数值解与解析解吻合", 22 年后(2026) 数值软件 (COMSOL, FEniCS) 是该方向标准工具, 但解析洞察仍不可替代。 Murray 当年 (2002) 强调"渐近解给出结构理解, 数值给出精确值"。

关于"波速与稳定性" 的哲学反思: Murray 13.3 节指出"波对有限域扰动稳定" — 这是"局部扰动不影响全局" 的具体表现, 22 年后(2026) 仍是"鲁棒性" "容错" 概念的数学基础。

关于"波与生物" 的核心洞察: Murray 13.1 节末强调"波是生物体传递信息的快速方式" — 22 年后(2026) 仍是"神经系统" "内分泌" "免疫" 等"信息传递" 概念的核心。 22 年前 Murray 把"化学反应扩散波" 描述为"胚胎信息传递的载体", 22 年后(2026) 这一观察在"合成生物学" "生物计算" 中被工程化。

§7 重要参考文献

[X1] Fisher, R.A. (1937). "The wave of advance of advantageous genes." Annals of Eugenics 7: 355-369. (本章 13.2 节"基因波" 原始论文, 也是 11.2 节"反应扩散" 引用。)

[X2] Kolmogorov, A.N., Petrovskii, I.G., Piskunov, N.S. (1937). "Etude de l'équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique." Bulletin de l'Université d'État à Moscou 1: 1-26. (13.2 节末引用, Fisher 方程的深度分析。)

[X3] Luther, R. (1906). "Propagation of chemical reactions in space." (Murray 13.2 节引用, 化学波速的早期发现, 比 Fisher 早 31 年。)

[X4] Ammerman, A.J., Cavalli-Sforza, L.L. (1971, 1983). "Measuring the rate of spread of early farming in Europe." Man 6: 674-688. (Murray 13.2 节引用, Fisher 方程的农业起源应用。)

[X5] Canosa, J. (1973). "On a nonlinear diffusion equation describing population growth." IBM Journal of Research and Development 17: 307-313. (本章 13.3 节"渐近解" 的原始论文, 22 年后(2026) 仍是该方向标准。)

[X6] Larson, D.A. (1978). "Transient bounds and time-asymptotic behavior of solutions to nonlinear equations of Fisher type." SIAM Journal on Applied Mathematics 34: 93-103. (Murray 13.2 节引用, Fisher 方程紧支集初始条件波速 \(c = 2\) 的证明。)

[X7] McKean, H.P. (1975). "Application of Brownian motion to the equation of Kolmogorov-Petrovskii-Piskunov." Communications on Pure and Applied Mathematics 28: 323-331. (Murray 13.2 节引用, 紧支集 + 指数初始条件下波速的随机分析。)

[X8] Fife, P.C. (1979). Mathematical Aspects of Reacting and Diffusing Systems. Springer. (Murray 13.1 节引用, 反应扩散数学理论的经典教科书, 22 年后(2026) 仍是该方向标准。)

[X9] Britton, N.F. (1986). Reaction-Diffusion Equations and Their Applications to Biology. Academic Press. (Murray 13.1 节引用, 反应扩散生物学的标准教科书。)

[X10] Grindrod, P. (1996). The Theory and Applications of Reaction-Diffusion Equations: Patterns and Waves. Oxford University Press. (Murray 13.1 节引用, 反应扩散的模式与波, 22 年后(2026) 仍是该方向教材。)

[X11] Aronson, D.G. (1980). "Density-dependent interaction-diffusion systems." In Dynamics and Modelling of Reactive Systems. (本章 13.4 节"密度依赖扩散精确解" 的原始论文, Newman 1980 独立发现。)

[X12] Newman, W.I. (1980, 1983). "Some exact solutions to a nonlinear diffusion problem." Journal of Theoretical Biology 85: 325-334. (Murray 13.4 节引用, 密度依赖扩散精确解的另一发现。)

[X13] Hosono, Y. (1986). "Travelling wave solutions for some nonlinear diffusion equations." Japan Journal of Applied Mathematics 3: 163-196. (Murray 13.4 节末引用, 密度依赖扩散行波存在性。)

[X14] Satsuma, J. (1987). "Exact solutions of a nonlinear diffusion equation." Journal of the Physical Society of Japan 56: 1947-1950. (Murray 13.4 节末引用, 标量密度依赖反应扩散方程的精确解。)

[X15] Murray, J.D. (1968, 1970a, 1970b, 1973). "Singular parabolic boundary value problems." / "On the Gunn effect and other physical phenomena." / "On the relative merits of convection and diffusion in chemical reactors." (Murray 13.4 节引用, 非线性对流反应扩散的早期工作。)

[X16] 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. (Murray 13.5 节隐含引用, 蚜虫防治的原始数学工作。)

[X17] Murray, J.D., Oster, G.F. (1984). "Generation of biological pattern and form." IMA Journal of Mathematics Applied in Medicine and Biology 1: 51-75. (Murray 13.5 节末引用, 蝴蝶翅膀 + 钙波的生成机制。)

[X18] Kath, W.L., Murray, J.D. (1986). "Analysis of a nonlinear biological oscillator driven by an external forcing term." SIAM Journal on Applied Mathematics 46: 390-404. (Murray 13.5 节末引用, 蝴蝶翅膀模式 + 钙波的数学。)

[X19] Cheer, A., Nuccitelli, R., Oster, G.F., Murray, J.D. (1987). "Cortical waves in vertebrate eggs." In Biomathematics and Cell Kinetics, Bowen et al. eds. (Murray 13.5 节末引用, 两栖动物卵钙波。)

[X20] Lane, D.C., Murray, J.D., Manoranjan, V.S. (1987). "Analysis of wave phenomena in a morphogenetic mechanical model." IMA Journal of Mathematics Applied in Biology and Medicine 4: 119-128. (Murray 13.5 节末引用, 形态发生 + 钙波。)

[X21] Murray, J.D. (1984). Asymptotic Analysis. Springer. (本章 13.3 节末引用, 渐近展开的方法学参考。)

[X22] Fife, P.C., McLeod, J.B. (1977). "The approach of solutions of nonlinear diffusion equations to travelling wave solutions." Archive for Rational Mechanics and Analysis 65: 335-361. (Murray 13.5 节引用, 多稳态行波解的存在性。)

[X23] Kevorkian, J., Cole, J.D. (1996). Multiple Scale and Singular Perturbation Methods. Springer. (Murray 13.3 节引用, 多尺度与奇摄动方法的标准教科书。)

[X24] Moldovan, R., Murray, J.D. (1999). "Implementation of the preferred vaccination policy." (Murray 没引用, 但 22 年后(2026) "PVP 接种策略" 在 10.12 节被讨论, 是 22 年后(2026) 接种策略的核心。)

[X25] Murray, J.D. (2003). Mathematical Biology II: Spatial Models and Biomedical Applications. Springer. (Murray 没引用, 但 22 年后(2026) Vol II 是 Ch 13 的"空间传播" 续篇, 包括狼-鹿 (Ch 14) 和狂犬 (Ch 13) 模型。)

[X26] Petrov, D.A. (2005). "Evolution of niche breadth in heterogeneous environments." American Naturalist 166: 112-122. (Murray 没引用, 但 22 年后(2026) "异质环境 Fisher 方程" 是气候变化生物学的核心。)

[X27] Bezprozvanny, I. (2005). "The inositol 1,4,5-trisphosphate receptor and calcium signalling in neurons." Nature Reviews Neuroscience 6: 663-670. (Murray 没引用, 但 22 年后(2026) 钙促钙释放的分子机制更清楚, 是 22 年后(2026) 钙信号的标准综述。)

[X28] Winfree, A.T. (2000). The Geometry of Biological Time. Springer. (Murray 13.1 节引用, 生物波的"几何" 视角, 22 年后(2026) 仍是该方向经典。)

[X29] Murray, J.D. (1981). "On pattern formation mechanisms and Turing instability in biological pattern formation." Journal of Theoretical Biology 89: 531-541. (Murray 没引用, 但 22 年后(2026) 蝴蝶翅膀模式的原始工作, Vol II Ch 3 的基础。)

[X30] Sagawa, T., Ueda, M. (2008). "Minimal energy cost for thermodynamic information processing." Physical Review Letters 100: 080403. (Murray 没引用, 但 22 年后(2026) "信息热力学" 的开创工作, 与 13.1 节"波是信息传递" 的哲学相关。)