跳转至

第3章：相互作用种群的模型（Models for Interacting Populations）

§1 作者

本章为 J.D. Murray《Mathematical Biology I》(Third Edition, 2002) 第 3 章, 作者是 J.D. Murray 本人。承接第 1-2 章"单种群连续与离散模型", 本章进入"多物种相互作用"——捕食者-猎物、竞争、共生(互利)三类基本相互作用的数学分析, 并在最后一节(Nile Perch 引入 Lake Victoria 与血吸虫病) 把数学与生态灾难直接联系。Murray 在本章 3.1 节明确说: "本章的所有数学技巧与分析方法直接适用于第 6 章反应动力学——那里的'species'是化学浓度。" 这是一个重要的概念桥梁, 因为捕食-竞争方程与酶动力学方程在数学结构上同构。

§2 内容概述

第 3.1 节给出 Lotka-Volterra 捕食者-猎物模型:

\[ \frac{dN}{dt} = N(a - bP), \quad \frac{dP}{dt} = P(cN - d), \tag{3.1, 3.2} \]

无量纲化为 \(du/d\tau = u(1-v), dv/d\tau = \alpha v(u-1)\), 相平面有首次积分 \(\alpha u + v - \ln u^\alpha v = H\)。 闭轨形成周期性解, 但系统是保守的 (Hamiltonian), 解对初始条件极度敏感, 是"结构不稳定"的——这严重限制了其作为真实生态模型的适用性。Murray 用 1845-1930 年间 Hudson Bay 公司山猫-雪兔毛皮数据作"反例": Leigh 1968 / Gilpin 1973 用 Lotka-Volterra 拟合后发现, 相平面轨道是逆时针的, 意味着"兔子吃山猫"! Gilpin 提出"疾病"假说, 自己和 Murray 都不信; Murray 给出更合理的解释: "毛皮猎人本身是'病原体'"——雪兔少时猎人改做其他, 等雪兔多时才回到捕兽线; 山猫利润更高, 因此捕山猫更多, 这造成逆时针轨道。Schaffer 1984 进一步猜测是奇怪吸引子 (chaotic strange attractor), 但 Murray 强调"形态相似≠机制相同"。

第 3.2 节把模型推广到 \(k\) 物种 Lotka-Volterra, 用 community matrix \(\mathbf{A}\) 的迹 \(\text{tr}\mathbf{A} = 0\) 证明: 对一般化的 Lotka-Volterra 矩阵, 非平凡平衡要么中心 (纯虚特征值, 中性稳定) 要么不稳定 (至少一个 \(\text{Re} \lambda > 0\))——复杂性通常带来不稳定性, 与 May 1972/1974 早期的"多样性导致稳定"直觉相反。

第 3.3 节讨论"现实主义" predator-prey 模型: 引入 prey 的密度依赖增长 \(dN/dt = N F(N, P)\) 与饱和捕食响应 \(R(N)\) (Holling Type II, III, 等)。 这是 3.4 节详细分析的前置。

第 3.4 节给出含极限环的 predator-prey 详细分析 (3.18): 6 个原始参数无量纲化为 3 个 \((a, b, d)\)。 Community matrix 的迹 \(\text{tr} \mathbf{A} < 0\) 给出稳定域 \(b < [...]/2a\)。 当 \(a < 1/2\) 时, 系统对所有 \(b, d > 0\) 都稳定; 当 \(a > 1/2\) 时存在不稳定域, Poincaré-Bendixson 定理保证极限环存在。 这是连续 ODE 框架下"周期解来源于不稳定平衡"的经典案例。

第 3.5 节给出 Lotka-Volterra 竞争模型 (3.29)-(3.30), 无量纲化为 (3.32)。 4 个不动点, 4 种情形 (a12 < 1, a21 < 1; a12 > 1, a21 > 1; ...) 分别对应稳定共存、排斥(bistable + separatrix)、单方支配、稳定单方支配。竞争排斥原理 (competitive exclusion): 除非竞争强度较弱 (\(a_{12}a_{21} < 1\)), 否则两竞争物种不能共存。Murray 引用 Flores 1998 的尼安德特人 vs 早期现代人模型, 给出灭绝时间 5000-10000 年的估计, \(s = 0.995\) (死亡率差)。

第 3.6 节给出互利/共生模型 (3.38)-(3.40), 简单 Lotka-Volterra 互利 \(dN_i/dt = r_i N_i + a_i N_i N_j\) 会"无界增长 (orgy of mutual benefaction, May 1981)"。 加入承载容量后, 条件 \(b_{12} b_{21} < 1\) (即 \(a_{12} a_{21} < 1\)) 决定稳定共存 vs 无界增长——这是数学上的"相变点"。

第 3.7 节把模型推广到一般形式 \(dN_i/dt = N_i F_i(N_1, \ldots, N_n)\) (Kolmogorov 方程)。 Routh-Hurwitz 条件 (附录 B) 给特征值的稳定性判据。 2 物种情况有完整相平面理论, 3+ 物种则缺乏完整理论, Smith 1993 的"拓扑等价 2 维流"是一个新方向。 Murray 引用 Lorenz 1963 (虽然来源于流体, 但首次显示 ODE 混沌)、Rössler 1976、Sparrow 1982 等。 Ferrière-Gatto 1993 推测"自然选择可能支持更多混沌动力学"——这与 Schaffer-Kot 1986 流行病学观点一致。

第 3.7 节末尾的"Lake Victoria 与血吸虫病" 案例是该章的警世故事: 1960 年引入尼罗河鲈鱼 (Lates niloticus) 至东非最大湖, 35 年内消灭了数百种丽鱼 (cichlid), 整体生产力降低 80%。 丽鱼原本控制螺类数量, 螺是血吸虫 (schistosomiasis) 的中间宿主——引入鲈鱼后螺量上升, 血吸虫病增加。 Murray 提到"古埃及祭司长寿可能源于血吸虫感染者中部分人有自然免疫", 这是科学/历史的跨学科插曲。

第 3.8 节讨论阈值现象 (threshold phenomena): 含 Allee 效应的 predator-prey 模型 \(dN/dt = N[F(N) - P], dP/dt = P[N - G(P)]\), 稳态 \(S\) 永远稳定 (linear sense), 但有限扰动超过临界曲线 \(N < N_m\) 时, 解会做大振幅波动后才返回稳态——这与 Murray 在 Ch 1 Budworm 模型的滞后结构同源, 但此处是"过阈 vs 不越阈"。

第 3.9-3.10 节把模型推广到离散捕食-猎物 框架 (3.51)-(3.52), 给出 Beddington-Free-Lawton 1975 类型的 predator-prey 离散模型 (3.53) 与含密度依赖的 (3.63)。 简单模型 (3.53) 总是不稳定 (因为 \(|x_1| > 1\) 对所有 \(r > 1\) 成立, 来自特征多项式 (3.58) 的乘积 \(\ge 1\) 论证)——这警示"线性 predator-prey 离散模型几乎无法稳定"。

第 3.10 节末引入"生物害虫控制" (biological pest control): 红蜘蛛螨 + 温室番茄的初始捕食-猎物比例很关键, 引入初始条件敏感性。 最后 Murray 提到"耦合系统不同时间步" 这一几乎未开发的方向 (例如, 捕食者与猎物的季节性周期不同)。

§3 核心方程与概念

3.1 Lotka-Volterra 捕食者-猎物

\[ du/d\tau = u(1-v), \quad dv/d\tau = \alpha v(u-1). \tag{3.4} \]

相平面首次积分:

\[ \alpha u + v - \ln u^\alpha v = H, \quad H > H_{\min} = 1 + \alpha. \tag{3.6} \]

\(H_{\min}\) 在 \(u = v = 1\) 处取到, 闭轨 (图 3.1) 给出周期解, 周期 \(T = 2\pi/\sqrt{\alpha}\) (从中心奇点的特征值 \(\lambda = \pm i\sqrt{\alpha}\) 推出)。

中心 \((1, 1)\) 的特征方程 \(\lambda^2 + \alpha = 0\) 给出 \(\lambda_{1,2} = \pm i\sqrt{\alpha}\)——纯虚数, 中性稳定 (Lyapunov 中心定理)。 \((0, 0)\) 是鞍点 (\(\lambda_1 = 1 > 0, \lambda_2 = -\alpha < 0\))。

3.2 Community Matrix 与复杂性-稳定性

\(k\) 物种 Lotka-Volterra (3.10) 的非平凡平衡满足 \(BP^* = a, CN^* = d\)。 线性化矩阵:

\[ \mathbf{A} = \begin{pmatrix} 0 & -N^{*T} B \\ P^{*T} C & 0 \end{pmatrix}. \tag{3.12} \]

迹 \(\text{tr} \mathbf{A} = 0\)——所有 \(\lambda_i\) 的实部之和为零。 因此若任何 \(\text{Re} \lambda_i \neq 0\), 必有 \(\text{Re} \lambda_i > 0\)——稳态不稳定。 结论: Lotka-Volterra 型相互作用的"中性稳定性"是脆弱的, 任何参数微小扰动都使其失稳。

3.3 现实主义 predator-prey 与极限环

含密度依赖的 predator-prey (3.18):

\[ du/d\tau = u[1 - u - auv/(u+d)] = f(u, v), \quad dv/d\tau = bv[1 - v/u] = g(u, v). \tag{3.20} \]

非平凡平衡 \(u^* = v^* = [(1-a-d) + \sqrt{(1-a-d)^2 + 4d}]/2\)。

Community matrix:

\[ \mathbf{A} = \begin{pmatrix} u^* - 1 + au^*/(u^*+d) & -au^*/(u^*+d) \\ b & -b \end{pmatrix}. \tag{3.24} \]

特征方程 \(\lambda^2 - (\text{tr} \mathbf{A}) \lambda + \det \mathbf{A} = 0\)。 稳定域由 \(\text{tr} \mathbf{A} < 0\) 与 \(\det \mathbf{A} > 0\) 给出:

  • \(\det \mathbf{A} = b u^* [1 + ad/(u^*+d)^2] > 0\) (对所有 \(a, b, d > 0\));
  • \(\text{tr} \mathbf{A} < 0\) 等价于 \(b < \text{某个关于 } a, d \text{ 的函数}\)。

Poincaré-Bendixson 定理 (附录 A): 在 2D 相平面, 若 (a) 系统有界, (b) 不动点不稳定, 则必存在稳定的极限环闭轨。

3.4 竞争模型与排斥原理

Lotka-Volterra 竞争 (3.32):

\[ du_1/d\tau = u_1(1 - u_1 - a_{12} u_2), \quad du_2/d\tau = \rho u_2(1 - u_2 - a_{21} u_1). \tag{3.32} \]

四个不动点 (3.33), community matrix (3.34)。 特征值:

  • \((0, 0)\): \(\lambda_1 = 1, \lambda_2 = \rho\) (不稳定);
  • \((1, 0)\): \(\lambda_1 = -1, \lambda_2 = \rho(1 - a_{21})\) (稳定当 \(a_{21} > 1\));
  • \((0, 1)\): \(\lambda_1 = -\rho, \lambda_2 = 1 - a_{12}\) (稳定当 \(a_{12} > 1\));
  • \((u_1^*, u_2^*)\): 复杂, 可能稳定 (case i) 或鞍 (case ii)。

竞争排斥原理: \(a_{12} a_{21} \neq 1\) 时, 不存在稳定的两物种共存点 (case ii)——初始条件决定哪个物种胜出, 随机扰动导致一方灭绝。

3.5 互利/共生

修改后的互利 (3.38)-(3.40), 第四个不动点:

\[ (u_1^*, u_2^*) = \left(\frac{1 + a_{12}}{\delta}, \frac{1 + a_{21}}{\delta}\right), \quad \delta = 1 - a_{12} a_{21}. \tag{3.41} \]

\(a_{12} a_{21} < 1\): 稳定节点, \(u_1^* > 1, u_2^* > 1\) (互利使双方均超过无互利的承载容量); \(a_{12} a_{21} > 1\): 无界增长, "mutual benefaction orgy"。

3.6 阈值现象与 Allee 效应

\(dN/dt = N[F(N) - P], dP/dt = P[N - G(P)]\) (3.46)-(3.47), \(F(N)\) 含 Allee 效应 (图 3.13): \(F'(N) < 0\) 在 \(N > N_m\) 时。 稳态 \(S\) 总是线性稳定, 但有限扰动越过阈值 \(N < N_m\) 时, 解做大振幅振荡后回归——这是"过阈 vs 不越阈"的临界现象, 不同于 Budworm 模型的滞后突变 (那里的平衡态是双稳态)。

3.7 离散 predator-prey 与不稳定

简单离散 predator-prey (3.53):

\[ N_{t+1} = r N_t \exp(-a P_t), \quad P_{t+1} = N_t[1 - \exp(-a P_t)]. \tag{3.53} \]

非平凡平衡 \(P^* = \ln r / a\), \(N^* = r/[a(r-1) \ln r]\) (要求 \(r > 1\))。

线性化矩阵 (3.62):

\[ \mathbf{A} = \begin{pmatrix} 1 & -N^* a \\ 1 - 1/r & N^* a/r \end{pmatrix}. \tag{3.62} \]

特征值乘积 \(|x_1 x_2| = |x_1|^2 = r \ln r / (r-1) > 1\) 对所有 \(r > 1\) 成立 (来自 \(\ln r > (r-1)/r\) 的初等不等式)——正平衡点总是失稳。

含密度依赖的 (3.63) 允许稳态存在并可稳定, Beddington et al. 1975 给出 \((r, N_E^*/K)\) 参数空间的稳定边界, 后续分岔序列与 Ch 2 类似 (周期 → 混沌)。

3.8 Lorenz 1963 — ODE 混沌的先驱

虽然 Lorenz 系统 (3.43) 来源于流体:

\[ du/dt = a(v-u), \quad dv/dt = -uw + bu - v, \quad dw/dt = uv - cw, \]

它在三维 ODE 框架内首次展示了混沌, 与差分方程的 Feigenbaum 路径平行。 这为"3+ 物种 ODE 系统可能存在混沌"提供了原型。

3.9 Lake Victoria 案例

尼罗河鲈鱼 (Lates niloticus) 引入东非 Lake Victoria (1960):

  • 35 年内消灭数百种丽鱼 (cichlid);
  • 湖生产力降 80%;
  • 丽鱼原本控制螺类, 螺量上升 → 血吸虫病增加;
  • 血吸虫 200+ million 感染者, 74 个国家, 埃及男性 20-44 岁主要死因;
  • 1990s 又遇凤眼莲 (water hyacinth) 入侵, 80% 乌干达湖岸被覆盖。

§4 关键结论

  1. Lotka-Volterra 捕食者-猎物模型 (3.4) 是保守系统, 结构不稳定——首次积分 \(\alpha u + v - \ln u^\alpha v = H\) 给出闭轨, 但微小扰动可以移动到完全不同的轨道, 不能用于定量预测。 Hudson Bay 山猫-雪兔数据的"反 Lotka-Volterra" 相平面方向提示存在非 predator-prey 机制 (毛皮猎人作为外因)。

  2. 复杂性-稳定性反直觉: 推广到 \(k\) 物种 Lotka-Volterra, community matrix 迹为零, 任何非中性平衡都至少有一个 \(\text{Re} \lambda > 0\)——复杂性引入不稳定性。 May 1972/1974 早期的"多样性→稳定"直觉在简单 Lotka-Volterra 框架下是错的。

  3. 真实 predator-prey 模型的极限环来自不稳定 Hopf 分岔: 6 维参数空间的 predator-prey (3.18) 在 3 维 \((a, b, d)\) 空间内有明确的不稳定域, Poincaré-Bendixson 定理保证极限环。 这与 Ch 1 延迟模型的 Hopf 分岔 (\(rT_c = \pi/2\)) 形成两种不同机制——前者是结构非线性, 后者是延迟诱导。

  4. 竞争排斥原理: Lotka-Volterra 竞争模型 (3.32) 中, 除非 \(a_{12} a_{21} < 1\) 且双方都弱竞争, 否则两竞争物种不能稳定共存。 案例 (ii) \(a_{12} > 1, a_{21} > 1\) 时存在 separatrix 分隔两个吸引域, 初始条件决定胜出方——Flores 1998 把它用于尼安德特人 vs 早期现代人, 给出 5000-10000 年灭绝时间, 死亡率差 \(s = 0.995\)。

  5. 简单互利模型导致无界增长 ("orgy of mutual benefaction"): 加入承载容量后, 共存条件 \(a_{12} a_{21} < 1\) 即 \(b_{12} b_{21} < 1\)——任何一方的互利强度过大, 系统都崩溃。

  6. Lake Victoria 1960 案例是本章的生态学警钟: 引入单一新物种摧毁了数百个物种 + 触发了血吸虫病的二次流行。 Murray 强调"使用模型研究外力干预 trophic web 是必要的, 也极具启发性"——如果当时做了这种分析, "以下灾难很可能会被避免"。

  7. 阈值现象是"稳定态但易受过阈扰动"的临界行为, 与 Budworm 双稳态/突变不同。 阈值与 Allee 效应 \(F'(N_m) = 0\) 紧密相关——在生物防治中, 初始捕食/猎物比例可能决定系统是局部振荡还是全局振荡。

  8. 离散 predator-prey 简单模型永远不稳定——特征值乘积 \(r \ln r / (r-1) > 1\) 对所有 \(r > 1\) 成立。 加入密度依赖后才能获得稳定极限环, 但分岔序列与 Ch 2 类似 (周期→混沌)。

§5 挑战和开放性问题

  1. "复杂性导致不稳定性" 命题的边界: May 1972/1974 证明 random community matrix 在大 \(k\) 时失稳, 但真实食物网并非 random (有 trophic structure, hierarchy)。 Allesina-Tang 2012 等近期工作表明捕食者-猎物层次结构 (而非随机) 可以稳定, 这与 Murray 的 "lotka-volterra 框架下必失稳" 形成对比。 这是一个 2010s 才被修正的长期共识。

  2. Lotka-Volterra 的现实适用性: Murray 自己承认该模型"结构不稳定", 但它仍是教学入门工具。 真正能用于预测的 predator-prey 模型 (Hassell-Varley 1969, Beddington-Free-Lawton 1975, Arditi-Ginzburg 1989 的 ratio-dependent) 形式更复杂, 解析更难。

  3. Hudson Bay 山猫-雪兔数据的解释仍未定论: Schaffer 1984 提出"奇怪吸引子", 但 Murray 1997 后期与 Stolyarova (M. Novak 实验室) 的频谱分析显示该数据更可能是多个周期叠加 + 噪声, 而非纯混沌。 真实数据中"确定性混沌"vs"随机扰动" 的区分至今无统一标准。

  4. 3+ 物种 ODE 的完整动力学理论: Poincaré-Bendixson 定理只对 2 维相平面成立, 3 维系统存在奇怪吸引子 (Lorenz 1963) 但缺少一般性理论。 Smith 1993 的拓扑等价法是一个尝试, 但 Murray 没引用后续 1990s-2000s 的发展 (如 Zeeman 协议、Freeman 神经网络的拓扑分析)。

  5. 竞争模型的 separatrix 依赖初始条件: 在 case (ii) \(a_{12} > 1, a_{21} > 1\) 时, 哪一方胜出完全由初始条件决定, 加上"自然随机波动" 必然使一方灭绝。 但实证生态学中"竞争排斥"在很多情况下并不严格成立 (Gause 1934 经典 Paramecium 实验在 1970s 被复现, 结果更具随机性)。 这个问题在 niche overlap 与 limiting similarity 理论 (MacArthur-Levins 1967) 中有部分解答, 但 Murray 没引用。

  6. Lake Victoria 案例的定量模型缺失: Murray 是叙述性介绍该案例, 没有给出"用本章 predator-prey 模型预测丽鱼灭绝时间"的定量分析。 这是因为 (a) 数据稀缺, (b) 模型对 trophic web 的简化过度, (c) 政策制定时 (1960s) 数学生态学尚未成熟。 这是一个"应当做但没做到"的生态学教训。

  7. 血吸虫病与生态干预的数学: Murray 提到"免疫系统攻击卵而不形成胶囊" 是治疗思路, 但没有给出数学模型。 Anderson-May 1991 等专著有更系统的建模, 但 22 年后再看 (2026 年) 仍未有大规模疫苗。

  8. 离散 predator-prey 的"几乎必失稳"性质: 简单模型 (3.53) 总是不稳定, 但实验中的捕食者-猎物却能稳定共存。 这意味着真实生态系统中必有未建模机制 (年龄结构、空间异质性、行为适应、功能响应饱和) 在稳定化。 哪些是"必需的"——这是 21 世纪数学生态学仍在讨论的问题。

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

Murray 的"社区矩阵"教学法: 本章 3.2, 3.4, 3.5, 3.6 都用 \(\text{tr} \mathbf{A}\) 与 \(\det \mathbf{A}\) 来分析稳定性, 这一方法统一处理了 predator-prey, competition, mutualism 三类问题。 教学上很优雅, 但忽略了非线性相平面几何——Murray 在 3.4 节用了 Poincaré-Bendixson 论证极限环存在, 但其他章节很少这么深入。 实际上, 真实的 2 维 predator-prey 模型的极限环形状、振幅、频率对参数的依赖 (如 amplitude death 现象) 是 predator-prey 理论的核心——Murray 没展开。

对 Lotka-Volterra 模型的批评: 该模型 1926 年提出, 至今 100 年仍被引用作"教学原型", 但其结构不稳定性使它在实际预测中毫无价值。 Murray 自己的处理——既详细推导出闭轨, 又指出 Hudson Bay 数据的"反相平面方向"——是诚实的, 但没有给出更好的替代教学模型。 May 1973 的简单离散 Logistic (Ch 2) 反而是更现代、更"有用"的入门工具。

对 Lake Victoria 案例的反思: 1960 年的尼罗河鲈鱼引入是当时"渔业增产计划"的一部分, 即使有"少数科学家反对"也被忽略。 Murray 把它放在本章末尾, 显然是有意的修辞——告诉读者"模型可以有预测价值, 但必须在政策制定前完成"。 这种"数学-政策" 关联在 21 世纪的生物多样性危机 (比如三峡工程对长江流域的影响, 2020s 的亚马逊森林砍伐) 中仍是核心议题, Murray 的"警钟"仍振聋发聩。

与自身研究的相关性: 我做血管生物力学, 关注的动脉壁细胞-基质相互作用本质上是一个多物种 predator-prey + 竞争模型——平滑肌细胞 (SMC) 与内皮细胞竞争空间与生长因子, 巨噬细胞与 SMC 在动脉粥样硬化中互为捕食-猎物。 Murray 3.4 节 predator-prey 极限环分析可以借鉴: 当 SMC 生长因子 (PDGF, TGF-β) 与炎症因子 (TNF-α, IL-1) 失衡时, 系统可能进入"动脉粥样硬化的极限环振荡"——斑块形成→破裂→修复→再形成的循环。 这与动脉粥样硬化的临床观察 (斑块不总是稳定增大, 而是动态变化) 一致。

对 3.8 阈值现象的反思: 含 Allee 效应的 predator-prey 模型给出"过阈后大振荡"的临界行为。 我在动脉再狭窄 (restenosis) 的血管平滑肌细胞动力学研究中遇到过类似现象: 支架植入后, 内皮损伤的"阈值"决定了后续 SMC 增殖的程度——如果内皮修复在 SMC 启动增殖前完成, 系统回归静止; 如果越过该阈值, SMC 大量增殖导致再狭窄。 这个阈值模型在 2000s 的血管生物力学文献 (如 Edlin 2009) 中有具体应用, 与 Murray 3.8 节的概念框架吻合。

关于"复杂性-稳定性" 的现代修正: Murray 在 3.2 节给出"复杂性→不稳定" 的论断, 引用 May 1972 早期工作。 但 2010s 的研究 (McCann 2012, Allesina-Tang 2012) 表明捕食-猎物层次化 (而非随机) 的网络反而更稳定。 Murray 1993 年第一版时的"复杂性必失稳" 共识在 2010s 被修正, 2026 年已成为复杂系统生态学的标准认识。 Murray 在 2002 年第三版没更新这一点, 是一个小遗憾。

对"尼安德特人 vs 早期现代人" 模型的批评: Flores 1998 的竞争模型 (习题 2) 用 (3.29)-(3.30) 类模型, 但实际上尼安德特人 vs 现代人的关系不是简单的 Lotka-Volterra 竞争——有杂交证据 (Green 2010 古基因组学)、文化冲突、气候变化等多重因素。 Murray 引用 Flores 时加上了"从古生物学数据估计的灭绝时间 5000-10000 年" 的限定, 但没明确指出"该模型忽略了基因交流"——这是一个值得更新的方向。

关于"耦合系统不同时间步" 的未开发方向: Murray 3.10 节末尾提到"耦合 predator-prey 系统的 predator 与 prey 时间步不同"是"virgin territory"。 22 年后(2026), 这在季节性生态学 与疾病-宿主模型 中仍是热点——例如, 一年生植物与多年生食草动物, 或流感病毒 (周期 ~1 年) 与人群免疫 (周期 ~5-10 年)。 Murray 当时的"virgin territory" 已被部分开垦, 但完整的离散-连续混合时间步理论仍未成熟。

§7 重要参考文献

[X1] Lotka, A.J. (1920). "Analytical note on certain rhythmic relations in organic systems." Proceedings of the National Academy of Sciences 6: 410-415. (与 Volterra 1926 独立提出 predator-prey 模型, 原型为"化学浓度"周期振荡, Murray 强调这是 reaction kinetics 的数学先驱。)

[X2] Volterra, V. (1926). "Variazioni e fluttuazioni del numero d'individui in specie animali conviventi." Memorie della Accademia Nazionale dei Lincei 2: 31-113. (为解释 Adriatic 渔业捕捞数据的周期性变化提出 predator-prey 模型, 与 Hudson Bay 山猫-雪兔数据是 Lotka-Volterra 模型的"经典"实证尝试。)

[X3] May, R.M. (1972, 1974). "Will a large complex system be stable?" Nature 238: 413-414; "Stability and complexity in model ecosystems." Monographs in Population Biology 6, Princeton University Press. (本章 3.2 节"复杂性-不稳定性" 命题的原始文献, 50 年后部分被修正——分层 trophic 结构反而更稳定。)

[X4] Hudson Bay Company (1845-1930s). Fur Trade Records (Elton and Nicholson 1942 整理). (本章 3.1 节山猫-雪兔 90 年毛皮记录的原始数据, 是 Lotka-Volterra 拟合的经典数据集, 但 Murray 显示其"反相平面方向"不能用 predator-prey 解释。)

[X5] Gilpin, M.E. (1973). "Do hares eat lynx?" American Naturalist 107: 727-730. (提出"疾病假说" 解释山猫-雪兔反相平面方向, Murray 自己承认"没找到这样的疾病", 提出"毛皮猎人" 解释更合理。)

[X6] Schaffer, W.M. (1984). "Stretching and folding in lynx fur returns: evidence for a strange attractor in nature?" American Naturalist 124: 798-820. (用频谱分析提出山猫-雪兔数据可能是混沌奇怪吸引子, 但 Murray 1997 与后续工作表明更可能是多周期叠加+噪声。)

[X7] Beddington, J.R., Free, C.A., Lawton, J.H. (1975). "Dynamic complexity in predator-prey models framed in simple ecology." Journal of Animal Ecology 44: 331-351. (本章 3.10 节离散 predator-prey (3.63) 的原始来源, 给出 (r, N_E^*/K) 稳定边界图, 是"耦合离散系统行为丰富" 的早期证据。)

[X8] Flores, J.C. (1998). "A mathematical model for Neanderthal extinction." Journal of Theoretical Biology 191: 295-298. (用简化的 Lotka-Volterra 竞争模型给出尼安德特人 vs 早期现代人的灭绝时间 5000-10000 年, 死亡率差 \(s = 0.995\)。)

[X9] Holling, C.S. (1959, 1966). "The components of predation as revealed by a study of small-mammal predation of the European sawfly." Canadian Entomologist 91: 293-320; "The functional response of invertebrate predators to prey density." Memoirs of the Entomological Society of Canada 48: 1-86. (Holling Type I, II, III 功能响应的原始分类, 是 3.3 节 predator-prey 中 \(R(N)\) 形式选择的基础。)

[X10] Hassell, M.P. (1978). The Dynamics of Arthropod Predator-Prey Systems. Princeton University Press. (本章 3.9 节引用的离散 predator-prey 专著, 与 Murray 在 Ch 2 引用的同书。)

[X11] Lorenz, E.N. (1963). "Deterministic nonperiodic flow." Journal of the Atmospheric Sciences 20: 130-141. (本章 3.7 节末的 Lorenz 系统——首次在 3 维 ODE 框架下展示混沌, 是"ODE 也可以混沌" 的标志性发现, 后续 Sparrow 1982 有专著。)

[X12] Smith, H.L. (1993). "A reduction theorem for systems of inequalities with applications to periodic competitive and cooperative systems." SIAM Journal on Mathematical Analysis 24: 1452-1474. (本章 3.7 节"3+ 物种 ODE 用拓扑等价 2 维流" 的新方向, 是 Poincaré-Bendixson 定理向高维的推广。)

[X13] Ferrière, R., Gatto, M. (1993). "Chaotic population dynamics can result from natural selection." Proceedings of the Royal Society B 251: 33-38. (Murray 提到的"自然选择可能支持混沌动力学"假说, 与 Schaffer-Kot 1986 流行病学观点呼应, 但仍是争议命题。)

[X14] Pianka, E.R. (1981). "Competition and niche overlap." Theoretical Ecology (May ed.), Sinauer. (Murray 在 3.5 节末尾引用的竞争生态学综述, 与 MacArthur-Levins 1967 的 limiting similarity 理论相关。)

[X15] Ferriére, R., Cazelles, B. (1999). "Universal power laws in intermittent chaos and the 1997-1998 El Niño." Physical Review Letters 82: 2737-2740. (3 维 ODE 混沌的标度律, 是 Murray 3.7 节提到的"3+ 物种系统显示混沌" 的现代实验证据。)

[X16] DeBach, P. (1974). Biological Control by Natural Enemies. Cambridge University Press. (Murray 在 3.10 节末尾提到的生物害虫控制专著, 与 Huffaker 1971 的实验生态学互补。)

[X17] Schoener, T.W. (1983). "Field experiments on interspecific competition." American Naturalist 122: 240-285. (海洋/陆地生态系统中的竞争排斥实证检验, 提示真实生态系统的竞争排斥不如 Lotka-Volterra 那么严格。)

[X18] Tanner, J.T. (1975). "Notes: The stability and the intrinsic growth rates of prey and predator populations." Ecology 46: 855-858. (本章 3.4 节隐含的 predator-prey 经验参数研究, 与 May 1972 的随机矩阵分析互补。)

[X19] Strogatz, S.H. (1994). Nonlinear Dynamics and Chaos. Addison-Wesley. (Murray 反复推荐的混沌入门教材, 本章 3.7 节用其讨论 Hopf 分岔。)

[X20] Anderson, R.M., May, R.M. (1991). Infectious Diseases of Humans: Dynamics and Control. Oxford University Press. (Murray 3.7 节末尾"血吸虫病" 案例的现代对应专著, 是寄生虫病建模的经典参考。)