第9章：受扰与耦合振子及黑洞（Perturbed and Coupled Oscillators and Black Holes）

§1 作者

本章为 J.D. Murray《Mathematical Biology I》(Third Edition, 2002) 第 9 章,作者是 J.D. Murray 本人, 受 A.T. Winfree (Arizona State) 关于果蝇 eclosion 节律、心脏起搏、神经节律等系列工作启发。Murray 在引言中把本章作为 Ch 7-8 (振荡器) 的扩展——"外部扰动 + 耦合"对极限环振荡的影响。Winfree 1980 提出的"奇点 = 黑洞" (singular points = black holes) 是该章的核心概念, Murray 称之为"对内生振荡器理论的革命性贡献"——能在 (I, θ) 平面上找到使振荡完全消失 的精确刺激+相位。

§2 内容概述

第 9.1 节给出相位重置 (phase resetting) 的几何基础——把极限环振荡看作"环上的点匀速运动", 周期 \(T\), 相位 \(\theta\)。刺激 \(I\) 改变点位置/速度, 重置相位。 Murray 用单摆作为类比——"在正确相位施加正确冲量, 可以完全停止摆动"。

第 9.2 节给出Type 1 vs Type 0 相位重置曲线——

Type 1 (9.5): 刺激只改变角速度, \(d\theta/dt = 1 + I \cos 2\theta\)。新相位 \(\varphi\) 是旧相位 \(\theta\) 的单调函数 (\(d\varphi/d\theta > 0\)); 任何新相位都可达, 平均斜率 = 1。显式解 (9.8) 用 \(\tan \varphi = A \tan[TB + \tan^{-1}(A^{-1} \tan \theta)]\) 给出。

Type 0 (9.10): 刺激将点推出极限环, 位置变为 \(\rho \neq 1\)。 \(\rho \cos \varphi = \cos \theta, \rho \sin \varphi + I = \sin \theta\), 消去 \(\rho\) 得 (9.11) \(\tan \varphi = \tan \theta - I/\cos \theta\)。当 \(|I| > 1\), \(d\varphi/d\theta < 0\) 在某范围, 新相位 \(\varphi\) 不能覆盖完整 \([0, 2\pi)\)——平均斜率 = 0。 Type 0 与 Type 1 在 \(I = 1\) 处有分岔 (图 9.6)。

第 9.3 节给出黑洞 (black holes)——Type 0 出现的 \(I = 1\) 是奇点。在 \((I, \theta)\) 平面上, 等 \(\varphi\) 曲线从 \(I = 1, \theta = \pi/2\) 与 \(\theta = 3\pi/2\) 两个奇点汇出 (图 9.7)。在奇点处, 任何 \(0 \le \varphi \le 2\pi\) 都可达——刺激结果是"不可预测"。 Murray 类比"在钟摆正过最低点时施以正确冲量, 可以让它停死"——理论上精确可计算, 实际上极难实现。实际实验中只能找到"黑洞附近的小区域", 新相位在该区域内不确定。

第 9.4 节给出真实生物振荡器中的黑洞验证——

Best 1979 用 Hodgkin-Huxley 模型 (Ch 7.5) 数值实验: 弱刺激 (2 mV) → Type 1 曲线 (图 9.8a), 强刺激 (60 mV) → Type 0 曲线 (图 9.8b); 黑洞出现在中间电压区域 (图 9.9); 关键刺激 (5 mV 在 0.45 周期) 完美消除发放 (图 9.10)。
Jalife-Antzelevitch 1979 心脏起搏细胞: 短刺激 (10 ms) → Type 1, 中刺激 (50 ms) → Type 0; 130 ms 处的刺激完全抑制起搏 (图 9.12a)——是临床"心律失常终止" 的理论基础。
临床意义: Winfree 1983b 提出"突发心源性死亡" 与刺激-相位拓扑学有关; 纤颤可能是"螺旋波 + 起搏异常" 的组合。

第 9.5 节给出耦合振荡器——Marek-Stuchl 1975 双向 BZ 实验: 两 BZ 反应通过多孔壁耦合, 自主频率差小时, 相位锁定 (phase locking) 到常数差; 频率差大时, 节奏分裂 (rhythm splitting)——长期慢变, 短期快变。弱耦合模型 (9.16) 包含 \(\varepsilon\) 与 \(k\): \(k\) 是质量传递耦合, \(\lambda\) 表征两个振荡器的不相同性。

第 9.6 节给出萤火虫同步——Pteroptyx malaccae 雄萤周期性闪光 (~0.9 秒), 用以吸引雌萤。单个萤火虫能"相位锁定"到外部闪光, 但仅在外部频率 \(\omega_e\) 接近自然频率 \(\omega\) 时。简化模型 (9.19): \(d\theta/dt = \omega + I \sin(\theta_e - \theta)\), \(I\) 是同步强度。相位差 \(\varphi = \theta_e - \theta\) 满足 (9.22) \(\varphi' = \delta - \sin \varphi\), \(\delta = (\omega_e - \omega)/I\)。稳态: \(\delta = \sin \varphi_s\), 锁定条件 \(|\delta| < 1\) 即 \(|\omega_e - \omega| < I\)。 临界 \(\delta_c = 1\): 当 \(|\delta| > 1\) 时无稳态, 出现"节奏漂移" (phase drift), 周期 \(T = 2\pi/[I(\delta^2 - 1)^{1/2}] = 2\pi/[(\omega_e - \omega)^2 - I^2]^{1/2}\) (9.25)。这一模型与 Marek-Stuchl 实验直接对应。

第 9.7-9.8 节给出奇摄动分析——把每个振荡器用局部坐标 \((A_i, \theta_i)\) 表示, \(A_i\) 是垂直于极限环的位移, \(\theta_i\) 是沿极限环的相位。变换 (9.37) 给出新的 (A, θ) 系统, 关键函数 \(V, U\) (9.41, 9.44)。稳定性条件 \(\int_0^T \Phi(\sigma) d\sigma < 0\) (9.47)。

第 9.9 节给出两时间尺度展开——短时间 \(t\) + 长时间 \(\tau = \varepsilon t\), 把 \(A_i, \theta_i\) 展开为 \(A_i \sim {}_0 A_i + \varepsilon {}_1 A_i\) (9.50)。 O(1) 给出 \(0\theta_i = t + \psi_i(\tau)\), \(0 A_i\) 是 \(\Phi\) 周期 ODE 的特解。 O(ε) 给出有界性条件 (9.67), 要求 \(\mu_1 = \mu_2 = 0\), 得相移方程 (9.69):

\[ d\chi/d\tau = P(\chi) - \lambda \beta, \quad P(\chi) = H(-\chi) - H(\chi), \quad \chi = \psi_2 - \psi_1. \tag{9.69} \]

\(P(0) = 0, P'(0) = -2k\) (9.70)。

第 9.10 节分析相移方程——稳态 \(\chi_0\) 满足 \(P(\chi_0) = \lambda \beta\) (9.71), 线性稳定 \(P'(\chi_0) < 0\)。临界条件 (9.73): \(\min P < \lambda \beta < \max P\)。当 \(\lambda \beta\) 落在 \([\min P, \max P]\) 之外时, 出现 rhythm splitting——慢变+快变的"准周期" 行为 (图 9.17), 解释了 Marek-Stuchl 1975 实验观察。频率分裂 (9.78): 振荡器 1 频率 \(1 + \varepsilon[H(\chi) + \lambda \beta]\), 振荡器 2 频率 \(1 + \varepsilon H(-\chi)\)。

§3 核心方程与概念

3.1 极限环的"环" 几何

\((r, \theta)\) 极坐标 (9.4): \(dr/dt = r(1 - r), d\theta/dt = 1\), 极限环 \(r = 1\)。周期 \(T = 2\pi\)。任意点 \((r_0, \theta_0)\): \(r(t) = r_0 e^t/[(1 - r_0) + r_0 e^t]\), \(\theta(t) = \omega t + \theta_0\)。

3.2 Type 1 相位重置 (9.5)-(9.8)

\(d\theta/dt = 1 + I \cos 2\theta\), 刺激时间 \(T\)。显式解 (9.8):

\[ |I| < 1: \tan \varphi = A \tan[TB + \tan^{-1}(A^{-1} \tan \theta)], \quad A = (|1+I|/|1-I|)^{1/2}, \quad B = |1-I^2|^{1/2}. \]

特征: \(d\varphi/d\theta > 0\) 对所有 \(\theta\), 平均斜率 = 1。任何新相位 \(\varphi\) 都可达。

3.3 Type 0 相位重置 (9.10)-(9.11)

刺激 \(I\) 平行于 \(y\) 轴, \((r, \theta) \to (\rho, \varphi)\): \(\rho \cos \varphi = \cos \theta, \rho \sin \varphi + I = \sin \theta\)。消去 \(\rho\):

\[ \tan \varphi = \tan \theta - I/\cos \theta. \tag{9.11} \]

\[ d\varphi/d\theta = 1 - I \sin\theta/\cos^2\theta. \tag{9.12} \]

\(|I| < 1\): \(d\varphi/d\theta > 0\) → Type 1。 \(|I| > 1\): \(d\varphi/d\theta < 0\) 在 \(\sin\theta > 1/|I|\) 范围 → Type 0。在 \(\varphi = \pi/2, 3\pi/2\) 处新相位不能取某些值, 平均斜率 = 0。

3.4 黑洞 (9.13)

投影 \(\varphi = \varphi(\theta, I)\) 到 \((I, \theta)\) 平面: \(I = \sin\theta - \cos\theta \tan\varphi\) (9.13)。在 \(I = 1, \theta = \pi/2\) 与 \(I = 1, \theta = 3\pi/2\) 处, 任何 \(\varphi\) 都汇于一点——奇点 (图 9.7)。 物理意义: 精确的"正确刺激 + 正确相位" 不可预测, 真实系统必有"近奇点小区域" 行为不可预测。

3.5 萤火虫相位锁定 (9.19)-(9.25)

\[ d\theta/dt = \omega + I \sin(\theta_e - \theta), \quad \varphi = \theta - \theta_e, \quad d\varphi/d\tau = \delta - \sin\varphi, \quad \delta = (\omega_e - \omega)/I. \tag{9.19-9.22} \]

稳态: \(\delta = \sin\varphi_s\), 锁定条件 \(-1 \le \delta \le 1\)。锁定时频率 (9.25): \(T = 2\pi/[I(\delta^2 - 1)^{1/2}]\)。临界 \(|\delta| > 1\): 节奏漂移, \(\varphi\) 单调增加但速度变化。

3.6 弱耦合振荡器 (9.16)

\[ dx_i/dt = F(x_i, y_i) + \varepsilon[k(x_j - x_i) + \lambda f(x_i, y_i)], \quad i \neq j. \]

\(0 < \varepsilon \ll 1\) 弱耦合, \(k\) 质量传递系数, \(\lambda\) 表征两振荡器的不相同性。 \(\varepsilon = 0\): 各自极限环; \(\varepsilon \neq 0, \lambda = 0\): 相同耦合; \(\varepsilon \neq 0, \lambda \neq 0\): 不相同耦合。

3.7 局部坐标变换 (9.37)

\[ x_i = X(\theta_i) + \varepsilon A_i Y'(\theta_i), \quad y_i = Y(\theta_i) - \varepsilon A_i X'(\theta_i), \quad i = 1, 2. \tag{9.37} \]

\(A_i\) 是垂直于极限环的位移, \(\theta_i\) 是沿极限环的相位。变换后的方程 (9.45)-(9.46):

\[ dA_i/dt = \Phi(\theta_i) A_i + kU(\theta_i, \theta_j) + \lambda \varphi(\theta_i) + O(\varepsilon), \]

\[ d\theta_i/dt = 1 + \varepsilon[\{\Omega(\theta_i) - \Gamma(\theta_i)\}A_i + \lambda \gamma(\theta_i) + kr(\theta_i) + kV(\theta_i, \theta_j)] + O(\varepsilon). \]

3.8 相移方程 (9.69)

\[ d\chi/d\tau = P(\chi) - \lambda \beta, \quad P(\chi) = H(-\chi) - H(\chi), \quad \chi = \psi_2 - \psi_1. \tag{9.69} \]

\(P(0) = 0, P'(0) = -2k\)。稳态锁定条件: \(\min P(\chi) < \lambda \beta < \max P(\chi)\) (9.73)。落出该区间 → 节奏分裂 (rhythm splitting)。

3.9 频率分裂 (9.78)

\[ d\theta_1/dt = 1 + \varepsilon[H(\chi) + \lambda \beta], \quad d\theta_2/dt = 1 + \varepsilon H(-\chi). \tag{9.78} \]

稳态时两频率相等 (锁定), 节奏分裂时频率突变。

§4 关键结论

Type 1 与 Type 0 相位重置在 \(I = 1\) 处有清晰分岔——弱刺激 (\(I < 1\)) 总能重置相位; 强刺激 (\(I > 1\)) 不能重置到所有相位, 出现"不可达" 区域。这是任何 2 维极限环的普适现象, 22 年后(2026) 仍是神经动力学的基础知识。
黑洞 (Winfree 1970) 在 \((I, \theta)\) 平面存在, 表明某些精确的"刺激+相位"组合完全消除振荡——是"振荡→静止" 的临界点。真实生物系统中只能找到"近黑洞"——但这是 Jalife 1979 实验的"130 ms 抑制起搏" 的理论基础。
Best 1979 在 Hodgkin-Huxley 模型中数值证实了 Type 1 / Type 0 分岔与黑洞存在——这是数学理论先于实验的罕见案例。 5 mV 在 0.45 周期的刺激完全消除神经发放 (图 9.10)。
Jalife-Antzelevitch 1979 心脏起搏实验验证了强刺激 50 ms 导致起搏完全抑制 (图 9.12a, trace 3)——是"刺激终止心律失常" 的临床前理论基础。 22 年后(2026) 这一原理用于植入式心脏除颤器 (ICD) 与抗心动过速起搏 (ATP)。
Marek-Stuchl 1975 的 BZ 耦合实验与 Murray 的耦合振荡器分析 (9.69) 完美对应——相位锁定 vs 节奏分裂取决于 \(\lambda \beta\) 与 \(P(\chi)\) 的相对位置。数学理论精确预测了实验观察的"长期慢变 + 短期快变" 行为。
萤火虫相位锁定模型 (9.19) 是 Ch 7 神经发放 + Ch 9 弱耦合的最小实例。临界 \(|\delta| < 1\) 是同步条件, \(|\delta| > 1\) 节奏漂移。 22 年后(2026) "萤火虫同步" 仍是数学-生物互动的经典教学案例。
相移方程 (9.69) 是该章的最重要结果——Murray 用两时间尺度展开从 4 维 ODE (9.16) 化简为 1 维 ODE, 给出相移的解析演化。这是奇摄动方法应用于弱耦合振荡器的范本, Murray 9.7-9.10 节本身就是"奇摄动应用于振荡器" 的教学示范。
频率分裂 (9.78) 显示, 在相移锁定时两振荡器频率精确同步; 节奏分裂时频率在两值之间突变。这是临床"心律失常+除颤" 的数学基础, 也是神经同步 + 认知节律 研究的现代基础 (Lisman 2010 的 hippocampal theta-gamma 耦合)。

§5 挑战和开放性问题

真实神经元的"刺激响应" 比 Hodgkin-Huxley 复杂得多——Best 1979 只验证了 4 维 HH 简化模型; 真实神经元有 10+ 离子通道, 黑洞位置可能因通道动力学不同而变。 22 年后(2026) 这一方向与深部脑刺激 (DBS) 临床应用紧密相关 (Lozano 2019 综述)。
心律失常的"终止" 阈值——Jalife-Antzelevitch 1979 实验显示 130 ms 处的刺激可终止起搏, 但临床上起搏细胞群有空间分布, 终止需要"临界质量"——这与 Vol II Ch 1 螺旋波终止相关。 22 年后(2026) 这是计算电生理学 的核心问题 (Trayanova 2014 综述)。
萤火虫同步的"非零相位差" —— 简化模型 (9.19) 假设 \(h(\varphi) = I \sin\varphi\), 对称。真实萤火虫 \(h(\varphi)\) 不对称, 这导致同步后仍存在相位差 (Ermentrout 1991)。 22 年后(2026) 这一观察被网络科学重新研究 (Mirollo-Strogatz 1990 的 Kuramoto 模型)。
Kuramoto 模型是 22 年后(2026) 的标准框架——Murray 9.7-9.10 的弱耦合分析可以被 Kuramoto 1975 框架统一。 Murray 当年没明确引用 Kuramoto, 但他的 (9.69) 是 Kuramoto 框架的特例。 22 年后(2026) 这是同步现象 (Strogatz 2000, "Sync") 的统一理论。
Marek-Stuchl 实验的现代扩展——22 年后(2026) "耦合化学振荡器" 在化学计算 (chemical computing) 与微流控 中被广泛研究 (Adamatzky 2001)。实际耦合不再是单一隔膜, 而是网络化——Murray 9.7-9.10 的两振荡器框架被推广到 \(N\) 振荡器网络。
黑洞的"鲁棒性"——理论预测黑洞在精确点, 但实验中只能找到"近黑洞" 小区域。真实生物系统有噪声, 会使精确奇点模糊。 22 年后(2026) 这是"随机动力学" 的标准议题 (Lyapunov 指数 + 噪声)。
心脏除颤 (defibrillation) 的临床应用——Murray 9.4 节末提到"纤颤 = 螺旋波 + 起搏异常"。 22 年后(2026) 这是计算心脏病学的核心议题 (Trayanova-Nash 2006), 但数学还没完全解决"为什么有些刺激能终止纤颤而其他不能"。
强耦合 (即 \(\varepsilon\) 不再小) 的振荡器——Murray 9.5-9.10 严格限制在弱耦合。 22 年后(2026) 强耦合系统的"同步 → 混沌" 转变 (Pisarchik 2019 综述) 仍是活跃方向。

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

Murray 的"几何直觉" 教学: 该章是全书最"几何"的一章——把抽象 ODE 极限环映射为"环上的点", 刺激是"位移/速度冲量", 黑洞是"汇出等 \(\varphi\) 曲线"的点。这种几何直觉 比代数推导更易记, 是 22 年后(2026) 该章仍被广泛引用的原因。

Winfree 的开创性工作: Murray 在该章中多次引用 Winfree 1967-2000, 强调 Winfree 的"环上点" 几何是"对内生振荡器理论的革命性贡献"。 Murray 没明确说"Winfree 1970 提出黑洞"是哪个具体会议论文, 但概念来源是清楚的。 22 年后(2026) Winfree 的工作被视为"生物时间理论" 的核心。

对"模型-实验一致性" 的反思: Best 1979 的 HH 模型数值实验精确证实了 Type 1/0 分岔与黑洞, 这是模型预测先于实验的范例。但 22 年后(2026), 真实神经元的离子通道多样性 (10+ 种 K, Na, Ca 通道) 远超 HH 4 维——黑洞的位置与"消失振荡" 的精确阈值会有定量偏差。 Murray 没明确指出这一"实验-模型差距" 是后续工作方向。

关于"实验→临床" 的链: Jalife-Antzelevitch 1979 → ICD (植入式心脏除颤器) 在 1980s 临床应用, 22 年后(2026) ICD 已成为预防心源性猝死的金标准。 Murray 当年 (2002) 只能引用 Jalife 1979 的基础实验, 但他正确预测了临床应用方向。

对萤火虫模型的反思: (9.19) 是最简单的相位模型, 但 Ermentrout 1991 指出"对称 \(\sin\) 函数" 不足以解释零相位差同步——真实萤火虫 \(h(\varphi)\) 不对称, 同步后仍有小相位差。 Murray 当年 (2002) 没充分讨论这一限制, 22 年后(2026) Mirollo-Strogatz 1990 的Kuramoto 框架 + Ermentrout 的"脉冲耦合"是主流。

与自身研究的相关性: 我做血管生物力学, 不直接做振荡器动力学。但血管脉搏波的传播可视为耦合振荡器——主动脉 + 周围动脉的脉搏反射形成耦合振荡系统。在某些病理条件下 (主动脉瓣反流、动脉瘤) 出现异常同步——Murray 9.6 节的萤火虫模型可以借鉴到"脉搏同步" 的数学建模。这是一个未充分开发的应用方向。

对"rhythm splitting" 的现代修正: Murray 9.10 节给出节奏分裂的确定性 ODE 描述。 22 年后(2026) 真实神经节律有噪声——节奏分裂应被随机微分方程 (SDE) 描述, 且噪声可能"加速"或"延缓" 节奏分裂的临界点 (Lyapunov 指数 + 噪声诱导转变)。 Murray 2002 完全用确定性, 22 年后这是标准扩展。

关于"双时间尺度" 教学: Murray 9.9 节的 (9.50)-(9.68) 推导是"两时间尺度展开应用于弱耦合振荡器" 的范本。但推导非常代数 (几十行展开), 22 年后(2026) 用计算机代数系统 (Mathematica, Maple) 自动化是标准做法, Murray 当年手工推导的精神值得保留但技术已经过时。

关于"窦房结与房颤" 的现代理解: Murray 9.4 节末提到"纤颤 = 螺旋波 + 起搏异常"。 22 年后(2026) 这一理解已被计算电生理学 (Trayanova 2014, Fenton 2002) 详尽研究, 包括多尺度建模 (离子通道 → 细胞 → 组织 → 心脏整体) 与除颤能量阈值的定量预测。 Murray 当年的"心脏起搏 + 刺激" 简化模型是这一多尺度框架的开端。

关于"Topological Singularity" 的哲学意义: Winfree 的"黑洞" 在 \((I, \theta)\) 平面是一个拓扑奇点——在这一点, "新相位" 是"完全不确定的"。这种"完全不可预测" 的概念在确定性 ODE 中是惊人的——它意味着某些数学模型的可预测性是有限的, 即便模型是确定性的。 22 年后(2026) "奇点理论" 在数据科学/机器学习中重新被关注 (神经网络的损失景观中的鞍点/极小值)——Murray 9.3 节的"拓扑奇点" 是这一更广主题的具体实例。

§7 重要参考文献

[X1] Winfree, A.T. (1970, 1980, 1982, 1983a, 1983b, 1987, 2000). The Geometry of Biological Time. Springer. (本章 9.1-9.4 节"相位重置 + 黑洞" 概念的开创者, 是 22 年后(2026) 仍被引用的"生物时间几何" 经典著作。)

[X2] Best, E.N. (1979). "Null space in the Hodgkin-Huxley equations: a small problem for large-scale computation of oscillating neurons." Biophysical Journal 35: 533-546. (本章 9.4 节"HH 模型中黑洞" 的数值验证论文, 是 Type 1/0 分岔与黑洞在真实模型中的首次观察。)

[X3] Jalife, J., Antzelevitch, C. (1979). "Phase resetting and annihilation of pacemaker activity in cardiac tissue." Science 206: 695-697. (本章 9.4 节"心脏起搏细胞" 的实验验证, 是 22 年后(2026) ICD 临床应用的理论基础。)

[X4] Marek, M., Stuchl, I. (1975). "Synchronization of two oscillating chemical reactions." Biophysical Chemistry 3: 263-272. (本章 9.5 节"耦合 BZ 反应" 的原始实验, 是 Murray 9.7-9.10 节数学分析的目标案例。)

[X5] Hanson, F.E. (1978). "Comparative studies of firefly pacemakers." Federation Proceedings 37: 2158-2164. (本章 9.6 节"萤火虫相位锁定" 的实验, 给出"周期 ~0.9 秒" 与"频率调整 ~15%" 的实测数据。)

[X6] Buck, J., Buck, E. (1976). "Synchronous fireflies." Scientific American 234: 74-85. (本章 9.6 节"萤火虫同步" 的科普介绍, 是 Pteroptyx malaccae 现象的经典文献。)

[X7] Mirollo, R.E., Strogatz, S.H. (1990). "Synchronization of pulse-coupled biological oscillators." SIAM Journal on Applied Mathematics 50: 1645-1662. (Murray 没引用, 但 22 年后(2026) 该论文是"脉冲耦合振荡器" 的标准框架, 推广了 Murray 9.6 的连续正弦耦合。)

[X8] Guevara, M.R., Glass, L. (1982). "Phase locking, period-doubling bifurcations and chaos in a mathematical model of a periodically driven oscillator: a theory for the entrainment of biological oscillators and the generation of cardiac dysrhythmias." Journal of Mathematical Biology 14: 1-23. (Murray 9.6 节末"心脏细胞 + 刺激 + 混沌" 的原始工作, 是 22 年后(2026) 心律失常混沌理论的开端。)

[X9] Kuramoto, Y. (1975). "Self-entrainment of a population of coupled non-linear oscillators." In International Symposium on Mathematical Problems in Theoretical Physics, Springer. (Murray 没引用, 但 22 年后(2026) Kuramoto 模型是"耦合振荡器" 的标准框架, Murray 9.7-9.10 的 (9.69) 是其特例。)

[X10] Neu, J.C. (1979). "Coupled chemical oscillators." SIAM Journal on Applied Mathematics 37: 307-315. (本章 9.7-9.9 节"弱耦合 + 奇摄动" 的核心方法学参考, 是 Murray (9.16)-(9.69) 推导的源流。)

[X11] Ermentrout, G.B. (1991). "An adaptive model for synchrony in the firefly Pteroptyx malaccae." Journal of Mathematical Biology 29: 571-585. (Murray 9.6 节末"更一般的影响函数" 的来源, 给出萤火虫"非零相位差同步" 的更复杂模型。)

[X12] Guevara, M.R., Shrier, A., Glass, L. (1981). "Phase resetting of spontaneously beating embryonic ventricular heart cell aggregates." American Journal of Physiology 240: H1-H16. (Murray 9.6 节末"心脏细胞 + 刺激 + 混沌" 的实验, Guevara-Glass 1982 的实验依据。)

[X13] Murray, J.D. (1984). Asymptotic Analysis. Springer. (本章 9.7-9.10 节"两时间尺度展开" 的方法学参考, 是弱耦合振荡器奇摄动分析的基础。)

[X14] Krinsky, V.I. (1978). "Mathematical models of cardiac arrhythmias (spiral waves)." Pharmacology and Therapeutics B 3: 539-555. (Murray 9.1 节末引用——心脏纤颤 + 螺旋波, 是 Vol II Ch 1 的前导。)

[X15] Hoppensteadt, F.C. (1985). An Introduction to the Mathematics of Neurons. Cambridge University Press. (Murray 9.1 节末引用——呼吸 + 步态同步, 是 Ch 12 神经模型的早期工作。)

[X16] Keener, J.P., Glass, L. (1984). "Global bifurcations of a periodically forced oscillator." Journal of Mathematical Biology 21: 175-190. (Murray 9.6 节末"周期强迫振荡器" 的全局分岔, 是 22 年后(2026) 心律失常建模的标准参考。)

[X17] Strogatz, S.H., Stewart, I. (1993). "Coupled oscillators and biological synchronization." Scientific American 269: 102-109. (Murray 9.6 节引用, "耦合振荡器同步" 的科普综述, 是 22 年后(2026) 同步科学的开端。)

[X18] Rinzel, J., Ermentrout, G.B. (1983). "Analysis of neural excitability and oscillations." In Methods in Neuronal Modeling, Koch-Segev eds. (Murray 9.6 节引用, 神经振荡 + 萤火虫的早期模型。)

[X19] Glass, L., Mackey, M.C. (1988). From Clocks to Chaos: The Rhythms of Life. Princeton University Press. (Murray 9.6 节末引用——生物节律与混沌的科普专著, 与本章内容高度互补。)

[X20] Lisman, J. (2010). "Working memory: a binding problem." Science 329: 1059-1060. (Murray 没引用, 但 22 年后(2026) 神经科学中"theta-gamma 频率耦合" 与本章 9.6 节"萤火虫同步" 在数学上同源, 是记忆与认知的频率编码假说。)