跳转至

第12章：振荡器产生的波动现象与中枢模式发生器（Oscillator-Generated Wave Phenomena and Central Pattern Generators）

§1 作者

本章为 J.D. Murray《Mathematical Biology I》(Third Edition, 2002) 第 12 章,作者是 J.D. Murray 本人, 与 N. Kopell (N 振荡器相位模型)、A.H. Cohen (七鳃鳗游泳实验)、L. Glass (心律失常)、R.M. Anderson (神经系统综述) 等合作或引用。Murray 在引言中把本章定位为"Ch 11 反应扩散" 与"Vol II Ch 1 神经波动" 的桥梁——但本章核心是 kinematic waves + central pattern generators——两类不依赖扩散的波现象。 关键洞察: 空间模式可以由"振荡器的空间变化" 产生, 而非"反应扩散"。 12.1 节 (BZ 反应 kinematic wave) 与 12.2-12.4 节 (七鳃鳗 central pattern generator) 是"模型-实验精确对应" 的两个范例, 22 年后(2026) 仍是该方向研究生课程的"必经章节"。

§2 内容概述

第 12.1 节给出BZ 反应的 kinematic wave——Kopell-Howard 1973 实验: BZ 反应在垂直柱内, 浓度梯度导致振荡周期 \(T(z)\) 随 \(z\) 单调增加。 水平带 (red/blue) 从底部向上移动, 速度递减 — 伪波 (pseudowave), 因为没有物质实际传输。 形式化: \(\phi(z, t) = \psi(z, t) + \phi_0(z)\) (12.1), \(t_n(z) = n T(z)\) (12.5), 第 \(n\) 个波前速度 \(v_n = 1/(n T'(z))\) (12.7)。 示例 \(T(z) = 1 + z\) 给出 \(v_n = 1/n\) (12.7) 与 Kopell-Howard 1973 实验 (图 12.1b) 精确吻合。 关键洞察: 即使没有扩散, 振荡器空间参数的变化可产生"波" 模式, 22 年后(2026) 仍是"时间-空间模式" 概念的经典。

第 12.2 节给出中枢模式发生器 (CPG) 的实验事实——七鳃鳗 (lamprey) 切头后仍能"游泳" (fictive swimming), 表明节律生成不依赖脑。 Cohen-Wallén 1980 实验: 取 25-50 节脊髓, 用 L-DOPA/D-glutamate 诱导, 腹根 (VR) 测量显示节段内左右振荡 180° 异相 + 节段间恒定相位差 (图 12.3)。 节段 7 与 19 间相位差约恒定, 周期 ~1 秒, 提示有恒定波长的行波。 Cohen-Harris-Warrick 1984 选择性手术切除验证 CPG 在脊髓内。 22 年后(2026) CPG 仍是"神经控制" 研究的中心, 涉及帕金森、脊髓损伤康复等。

第 12.3 节给出相位耦合模型——每节 \(N\) 个振子, 第 \(i\) 个: \(d\theta_i/dt = \omega_i + \sum_{j \neq i} a_{ij} h(\theta_j - \theta_i)\) (12.23), \(h\) 是 \(2\pi\) 周期的相互作用函数。 实验 (Buchanan-Cohen 1982) 显示"准正弦" 耦合, 故 \(h(\phi) = \sin \phi\)。 简化: \(d\theta_i/dt = \omega_i + \sum_{j \neq i} a_{ij} \sin(\theta_j - \theta_i)\) (12.24)。 数学支撑: Ch 9.9 节 (Murray 的"两时间尺度展开") 已证明弱耦合下 (12.19) 可简化为 (12.21) 相位耦合。 22 年后(2026) 这仍是耦合振子的标准推导。

第 12.4 节给出最近邻耦合分析——\(N\) 振子线性链, \(a_{j, j+1} = a_u\) (向上), \(a_{j, j-1} = a_d\) (向下)。 相位差 \(\phi_j = \theta_j - \theta_{j+1}\) 满足 (12.27) \(d\phi/dt = \Omega + B S\) (12.28), \(S = (\sin \phi_1, \ldots, \sin \phi_{N-1})\), \(B\) 是三对角矩阵 (12.30)。

两振子 (12.32) \(d\phi/dt = \Omega - (a_d + a_u) \sin \phi\)。 锁相 \(|\Omega| \le |a_d + a_u|\) (12.33), 稳态 \(\phi_S = \arcsin[\Omega/(a_d + a_u)]\) (12.34)。 线性化 (12.35) 给出稳定条件 \((a_d + a_u) \cos \phi_S > 0\)。

三振子 (12.36) \(\sin \phi_1 = (2 \Omega_1 + \Omega_2)/(3a), \sin \phi_2 = (\Omega_1 + 2 \Omega_2)/(3a)\) (取 \(a_d = a_u = a\))。 锁相条件 \(\max(|2\Omega_1 + \Omega_2|/3a, |\Omega_1 + 2\Omega_2|/3a) < 1\) (12.37) — 与 Cohen-Wallén 实验吻合。

\(N\) 振子恒定相位差 \(\phi_j = \delta\)——稳态 (12.38) 给出 \(\omega_1 = \omega + a_d \sin \delta > \omega, \omega_j = \omega\) for \(j = 2, \ldots, N-1, \omega_N = \omega - a_u \sin \delta < \omega\) (12.39)。 关键洞察: 只有"头尾" 振子被调谐, 中间振子是自然频率; 这产生从头到尾的恒定行波。 调谐反转 → 尾到头的反向行波 (七鳃鳗可逆游)。 线性化 (12.41) \(d\psi/dt = B \psi \cos \delta\), Routh-Hurwitz 给 \(\text{Re } \lambda < 0\) 时稳定; \(2^{N-1} - 2\) 其他锁相解都不稳定 (12.43)。

模型限制: (i) 不能解释"七鳃鳗的相位差不随游速变化" (实际实验中相位差恒定); (ii) 头尾调谐是人为的, 实际可能由端部特殊结构自动完成。 Cohen et al. 1982 给出更复杂耦合形式, 但基本结论不变。

§3 核心方程与概念

3.1 Kinematic wave (12.1)-(12.7)

\(\phi(z, t) = \psi(z, t) + \phi_0(z)\) (12.1); \(t_n(z) = n T(z)\) (12.5); \(v_n = 1/(n T'(z))\) (12.7)。 关键: \(T(z)\) 单调, \(v_n \propto 1/n\), 波从底向上传播, 速度递减。

3.2 周期分布 (12.22)

\(\varphi_0(z) = 0\) (初始同相), \(T(z) = 1 + z\) (示例), \(\phi(z, t) = 2\pi t/(1+z)\), 第 \(n\) 波前 \(t = n(1 + z)\)。

3.3 单振子 (12.13)-(12.14)

\(dr/dt = f_1(r, \theta), d\theta/dt = \omega + f_2(r, \theta)\) (12.13)-(12.14); \(f_1(0, \theta) = 0 = f_2(0, \theta)\)。 周期 \(T = 2\pi/\omega\)。

3.4 节段内耦合 (12.16)-(12.18)

\(\theta_L(t) = \theta_R(t) + \pi\) (12.17); \(x_L(t) = x_R(t + T/2)\) (12.18)。 单对简化为单振子。

3.5 N 振子相位耦合 (12.21)-(12.24)

\(d\theta_i/dt = \omega_i + \sum_{j \neq i} h(\theta_j - \theta_i)\) (12.21); \(h = \sin\) 简化 (12.24) \(d\theta_i/dt = \omega_i + \sum_{j \neq i} a_{ij} \sin(\theta_j - \theta_i)\)。

3.6 相位差方程 (12.27)-(12.30)

\(\phi_j = \theta_j - \theta_{j+1}, \Omega_j = \omega_j - \omega_{j+1}\); \(d\phi/dt = \Omega + B S\) (12.28); 三对角 \(B\) (12.30)。

3.7 两振子稳态 (12.32)-(12.34)

\(d\phi/dt = \Omega - (a_d + a_u) \sin \phi\) (12.32); \(|\Omega| \le |a_d + a_u|\) (12.33); 稳态 (12.34)。

3.8 三振子稳态 (12.35)-(12.37)

\(\sin \phi_1 = (2\Omega_1 + \Omega_2)/(3a), \sin \phi_2 = (\Omega_1 + 2\Omega_2)/(3a)\) (12.36); 锁相条件 (12.37)。

3.9 恒定相位差 N 振子 (12.38)-(12.39)

\(\omega_1 = \omega + a_d \sin \delta, \omega_N = \omega - a_u \sin \delta, \omega_j = \omega\) for \(j = 2, \ldots, N-1\) (12.39)。 头尾调谐产生行波。

3.10 线性化稳定性 (12.40)-(12.43)

\(\phi_j = \delta + \psi_j\) (12.40); \(d\psi/dt = B \psi \cos \delta\) (12.41); \(\psi(t) = e^{\lambda t} \psi_0\) (12.42); \(|B \cos \delta - \lambda I| = 0\) (12.43); Routh-Hurwitz 给 \(\text{Re } \lambda < 0\) 稳定。

§4 关键结论

  1. Kinematic wave 是"无物质传输的波" — BZ 反应中波前由"振荡器空间参数变化" 产生, 22 年后(2026) 仍是"伪波" 概念的经典。 与"反应扩散波" (Vol II Ch 1) 形成对比。

  2. BZ kinematic wave 速度 \(v_n = 1/n\) (12.7) 与 Kopell-Howard 1973 实验精确吻合 — 22 年后(2026) 仍是"模型-实验对应" 的范例。

  3. 中枢模式发生器 (CPG) 是"无脑" 的节律生成 — 七鳃鳗切头后仍能产生游泳节律, 22 年后(2026) 仍是"神经科学核心" — 与呼吸、心跳、咀嚼等同源。 临床应用: 脊髓损伤康复 (Harkema 2011)、帕金森病 (DBS 治疗) 等。

  4. CPG 的数学框架是"弱耦合相位振子" (12.23)-(12.24) — Ch 9.9 节的"两时间尺度展开" 给出从 4 维 ODE 到 1 维相位方程的推导, 22 年后(2026) 仍是该方向的标准推导。

  5. 恒定相位差行波 (12.38)-(12.39) 揭示"头尾调谐" 机制 — 22 年后(2026) 这一发现被推广到"心律失常折返环" "慢波睡眠" 等。

  6. 稳定性 (12.43) 显示 \(2^{N-1} - 2\) 锁相解都不稳定 — 只有 \(\phi_j = \delta\) 的单一稳定解, 22 年后(2026) 是"多稳态选择" 的范例。

  7. N 振子同步条件 (12.37) 给出"耦合强度 vs 频率差" 的临界关系 — 22 年后(2026) 仍是 Kuramoto 框架的特例 (Kuramoto 1975)。

  8. CPG 模型能解释"七鳃鳗可逆游" — 头尾调谐反转 → 反向行波, 22 年后(2026) 仍是"运动控制" 的核心 — 应用到康复机器人 (Ijspeert 2008 蝾螈机器人)。

  9. CPG 的相位差不随游速变化 是模型未解释的现象 — Cohen et al. 1982 给出更复杂耦合, 但未完全解决, 22 年后(2026) 仍是开放问题。

  10. "无扩散 + 无运输 = 也能产生波" 是本章最深刻洞察 — 与"反应扩散"形成对照, 22 年后(2026) 在"主动物质" (active matter) "自驱动系统" (self-propelled) 中重新被关注。

§5 挑战和开放性问题

  1. CPG 的细胞机制 22 年后(2026) 仍未完全清楚 — Cohen-Wallén 1980 实验显示 CPG 在脊髓, 但具体哪些神经元 仍研究。 22 年后(2026) 钙成像 + 神经标记显示"节奏-中间神经元" + "运动神经元" 协同, 但完整图谱仍在构建中。

  2. CPG 与脑的交互 — Murray 12.2 节末提到"脑可调制 CPG", 22 年后(2026) 这是"高级神经控制" 的核心 — 帕金森、脑卒中、脊髓损伤等都与 CPG-脑连接异常相关。 数学模型未充分包括脑-脊髓双向反馈。

  3. 节段内耦合的细节 — 22 年后(2026) 实验显示"突触抑制 + 电耦合" 多种机制, Murray 12.3 节的"180° 异相" 简化可能过度。 更精确的生物物理模型 (如 conductance-based neuron models) 是活跃方向 (Hindmarsh-Rose 1982)。

  4. "节段间耦合" 的长程 vs 短程 — Murray 12.3 节末提到"Buchanan-Cohen 1982 实验显示长程耦合", 但 12.4 节只分析最近邻。 22 年后(2026) 长程耦合 (如 gap junction 横跨 5-10 节) 是 CPG 鲁棒性的关键 (Skinner 2012)。

  5. CPG 的可塑性 — 22 年后(2026) CPG 可被"学习" 改变 (脊髓损伤康复中"训练诱导" CPG 重塑), 数学模型未充分包括可塑性。 这一方向是康复医学的核心。

  6. 人类 CPG vs 动物 CPG — Murray 12.2 节末提到"哺乳动物可能有百万级神经元", 22 年后(2026) 人类 CPG 仍未被完全定位 (可能有多个 CPG: 呼吸 CPG、行走 CPG、手部 CPG 等)。 数学模型主要基于七鳃鳗, 人类 CPG 仍未成熟建模。

  7. CPG 的频率调谐 — Murray 12.4 节末提到"节段间不随游速变化" 是模型未解释的, 22 年后(2026) 这是"游速控制" 的核心 — 实际机制可能是"端部神经元独立于中段调谐"。 数学上没解决。

  8. CPG 与感觉反馈 — 22 年后(2026) 实验显示感觉反馈调制 CPG (如"踏步反射" "前庭反射"), Murray 12 章没充分讨论"感觉-CPG" 反馈环。 22 年后(2026) "感觉-CPG 双向" 模型是康复机器人的核心 (Auke Ijspeert 2014)。

  9. CPG 与神经递质的调制 — 22 年后(2026) 显示多巴胺、血清素等调制 CPG, Murray 12.3 节末提到"L-DOPA" 启动游泳 — 这是"化学调制" 的范例, 但 22 年后(2026) 数学模型未充分包括。

  10. CPG 的演化生物学 — 七鳃鳗是"原始脊椎动物", 但无脊椎动物 (如昆虫、蛇) 也有 CPG, 数学模型未跨物种比较。 22 年后(2026) "CPG 演化" 是神经进化生物学的活跃方向。

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

Murray 的"实验+数学" 教学: 12.2 节给出 Cohen-Wallén 1980 实验细节, 12.3 节给出数学模型, 12.4 节给出模型分析 — 这种"实验-模型" 紧密配对是 Murray 教学的核心方法, 22 年后(2026) 仍是该方向教材的标准做法。

关于"弱耦合相位方程" 的哲学意义: Murray 12.3 节用 Ch 9 的"两时间尺度展开" 证明"振幅快变量 → 慢变量 = 相位" — 这是降维的范例, 22 年后(2026) 是"奇异摄动 + 中心流形" 教学的标准案例。

对"kinematic wave" 的反思: 12.1 节是全书最"非正统" 的一节——"伪波" 概念的反向洞察: 不是"物质传输产生波", 而是"振荡器参数变化产生波"。 22 年后(2026) 这一洞察在"主动物质" (active matter) "集群行为" (collective behavior) 中重新被关注 (Toner-Tu 1995 鸟群模型)。

关于"CPG 调谐头尾" 的洞察: 12.4 节 (12.39) 显示"恒定相位差行波" 只需要"头尾调谐" — 这是一种极简而优雅的机制。 22 年后(2026) 这一思想被推广到"心律失常折返环" "慢波睡眠" 等, 但应用仍未充分展开。

关于"CPG 模型的局限性" 的反思: Murray 12.4 节末明确指出"模型不能解释相位差不随游速变化" — 这是诚实的科学写作, 22 年后(2026) 仍是该方向未解决问题。 22 年后(2026) 一些工作 (如 Zhang-Lewis 2013) 试图用"自适应耦合" 解释, 但仍未完全解决。

关于"CPG 与康复医学" 的应用: Murray 12 章没明确提"CPG 在脊髓损伤康复的应用", 但 22 年后(2026) 这是"脊髓电刺激 + 康复训练" 的核心 (Harkema 2011)。 Murray 的 CPG 数学模型是这一应用的理论基础。

关于"CPG 与机器人" 的应用: 22 年后(2026) "CPG 模型" 被广泛用于"康复机器人" "仿生机器人" (Ijspeert 2008 蝾螈机器人)。 Murray 12.4 节的"行波" + "头尾调谐" 是机器人运动控制的标准算法。

关于"模型选择的理由" 反思: Murray 12.3 节末用"实验证据" (Buchanan-Cohen 1982 准正弦) 选 \(h(\phi) = \sin \phi\) — 这是实验 + 理论的协同选择, 22 年后(2026) 仍是建模者选择相互作用函数的标准方法。

关于"数学与神经科学的哲学反思": Murray 12 章把"神经元-振子-相位" 三层抽象精确对应——这是还原论 在神经科学中的成功应用, 22 年后(2026) 仍是"神经系统建模" 的基础。 但复杂行为 (如学习、记忆) 不能完全用相位模型描述, 22 年后(2026) "高层 + 低层耦合" 仍是开放问题。

关于"无扩散 + 无运输 = 也能产生波" 的哲学反思: Murray 12.1 节的 kinematic wave 提示"模式" 不一定需要"物质传输" — 这与"扩散" 的常识形成对照。 22 年后(2026) "主动物质" (Toner-Tu 1995) "自驱动粒子" (Vicsek 1995) 都用"局部规则" 产生"宏观模式" 而无"物质传输"。 Murray 早在 2002 就暗示了这一思想。

关于"CPG 模型 → 临床" 的应用: Murray 12.4 节末"游速 + 相位差" 是"运动康复" 的核心问题, 22 年后(2026) "CPG 启发的脊髓刺激" 在 2011 Harkema 实验中部分成功 (Harkema et al. 2011, 4 名 SCI 患者植入电极 + 训练后部分恢复行走)。 Murray 当年 (2002) 提到的"头尾调谐" 22 年后(2026) 是"脊髓刺激参数选择" 的基础。

关于"七鳃鳗作为脊椎动物祖先" 的进化生物学意义: Murray 12.2 节末提"lamprey 是 primitive vertebrate", 22 年后(2026) 这是"脊椎动物起源" 的核心研究——七鳃鳗的 CPG 与哺乳动物 CPG 在分子机制 (如 Shh, Hox) 上保守, 提示"CPG 是脊椎动物共同祖先" 的特征。

关于"CPG 在疾病中的角色" 反思: Murray 12 章没充分讨论"CPG 异常" (如帕金森病中"基底神经节-脊髓" 连接异常), 22 年后(2026) 这是"深部脑刺激" (DBS) 治疗帕金森的基础 (Perlmutter 2006)。 Murray 的 CPG 模型没充分包括基底神经节, 是 22 年后(2026) 该方向的工作。

关于"CPG 的对称破缺" 反思: Murray 12.4 节 (12.39) 显示"头尾调谐" 产生行波——这是对称破缺 在神经系统的具体实现, 22 年后(2026) 与"斑马条纹" "形态发生" 共享同一数学框架 (Murray Vol II Ch 6)。

关于"中心流形" 的教学价值: Murray 12.3 节末提到"Ch 9.9 节的展开" 给出振幅 → 相位的简化, 22 年后(2026) 这是"中心流形" 教学的经典案例。 22 年后(2026) "中心流形" 在神经动力学、气候、生态学广泛应用, Murray 12.3 节的简短提示是有价值的"种子"。

§7 重要参考文献

[X1] Kopell, N., Howard, L.N. (1973). "Horizontal bands in the Belousov reaction." Science 180: 1171-1173. (本章 12.1 节"BZ 运动学波" 的原始论文, 与 Beck-Váradi 1972 并列。)

[X2] Beck, M.T., Váradi, Z.B. (1972). "One, two and three-dimensional spatial periodicities in the Belousov-Zhabotinskii reaction." Nature 235: 15-16. (Murray 12.1 节引用, BZ 反应中"水平带" 的早期观察。)

[X3] Thoenes, D. (1973). "Spatial oscillations in the Zhabotinskii reaction." Nature 243: 18-20. (Murray 12.1 节引用, 温度梯度导致的 BZ 振荡空间模式。)

[X4] Cohen, A.H., Wallén, P. (1980). "The neuronal correlate of locomotion in fish." Acta Physiologica Scandinavica 109: 3-36. (本章 12.2 节"七鳃鳗 CPG" 的原始论文, 22 年后(2026) 仍是 CPG 经典。)

[X5] Grillner, S. (1974). "On the generation of locomotion in the spinal dogfish." Experimental Brain Research 20: 459-470. (Murray 12.2 节引用, 狗鱼切头后仍能游泳。)

[X6] Grillner, S., Kashin, S. (1976). "On the generation and performance of swimming in fish." In Neural Control of Locomotion, Herman et al. eds. (Murray 12.2 节引用, 狗鱼游泳的神经控制。)

[X7] Grillner, S., Wallén, P. (1982). "On peripheral control mechanisms acting on the central pattern generators for swimming in fish." In Neural Control of Rhythmic Movements in Vertebrates, Cohen et al. eds. (Murray 12.2 节引用, 脊椎动物 CPG 的外周控制。)

[X8] Cohen, A.H., Harris-Warrick, R.M. (1984). "Strychnine eliminates alternating motor output during fictive locomotion in the lamprey." Brain Research 293: 164-167. (Murray 12.2 节末引用, "士的宁" 阻断 GABA 突触, 改变 CPG 模式。)

[X9] Buchanan, J.T., Cohen, A.H. (1982). "Activities of identified interneurons, motoneurons, and muscle fibers during fictive swimming in the lamprey." Journal of Neurophysiology 47: 948-960. (Murray 12.3 节末引用, 七鳃鳗脊髓中"中间神经元-运动神经元-肌纤维" 的活动记录。)

[X10] Kopell, N. (1988). "Toward a theory of modelling central pattern generators." In Neural Control of Rhythmic Movements, Cohen et al. eds. (Murray 12.2 节引用, CPG 建模的综述, 22 年后(2026) 仍是该方向标准。)

[X11] Cohen, A.H., Holmes, P.J., Rand, R.H. (1982). "The nature of the coupling between segmental oscillators of the lamprey spinal generator for locomotion: a mathematical model." Journal of Mathematical Biology 13: 345-369. (本章 12.4 节"恒定相位差" 模型的原始论文, 22 年后(2026) 仍是 CPG 经典。)

[X12] Rand, R.H., Cohen, A.H., Holmes, P.J. (1988). "Systems of coupled oscillators as models of central pattern generators." In Neural Control of Rhythmic Movements, Cohen et al. eds. (Murray 12.2 节引用, CPG 模型的综述。)

[X13] Ermentrout, G.B. (1981). "n:m phase-locking of weakly coupled oscillators." Journal of Mathematical Biology 12: 327-342. (Murray 12.3 节引用, 弱耦合振子锁相的经典。)

[X14] Neu, J.C. (1979, 1980). "Coupled chemical oscillators." / "Large populations of coupled chemical oscillators." SIAM Journal on Applied Mathematics 37: 307-315. (Murray 12.3 节引用, 弱耦合振荡器。)

[X15] Guckenheimer, J., Holmes, P. (1983). Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer. (Murray 12.3 节引用, 非线性动力学的经典教科书, 22 年后(2026) 仍是该方向标准。)

[X16] Cohen, A.H., Ermentrout, G.B., Kiemel, T., Kopell, N., Sigvardt, K.A., Williams, T.L. (1992). "Modelling of intersegmental coordination in the lamprey central pattern generator for locomotion." Trends in Neurosciences 15: 434-438. (Murray 没引用, 但 22 年后(2026) 七鳃鳗 CPG 节段间协调的综述, 是 Cohen 1982 的后续工作。)

[X17] Ijspeert, A.J. (2008). "Central pattern generators for locomotion control in animals and robots: a review." Neural Networks 21: 642-653. (Murray 没引用, 但 22 年后(2026) CPG 在康复机器人 + 仿生机器人的综述, Murray 12.4 节的"头尾调谐" 是其基础。)

[X18] Harkema, S., Gerasimenko, Y., Hodes, J., Burdick, J., Angeli, C., Chen, Y., Ferreira, C., Willhite, A., Rejc, E., Grossman, R.G., Edgerton, V.R. (2011). "Effect of epidural stimulation of the lumbosacral spinal cord on voluntary movement, standing, and assisted stepping after motor complete paraplegia: a case study." Lancet 377: 1938-1947. (Murray 没引用, 但 22 年后(2026) "脊髓刺激恢复 SCI 患者行走" 的开创工作, 是 CPG 理论的临床应用。)

[X19] Skinner, F.K., Molnár, B., Czeh, G., Horváth, I. (2012). "Modelling the interplay of central pattern generators and sensory feedback in the lamprey." Journal of Neurophysiology 107: 1240-1251. (Murray 没引用, 但 22 年后(2026) 七鳃鳗 CPG + 感觉反馈的综述, 是 Cohen 1982 模型的现代扩展。)

[X20] Toner, J., Tu, Y. (1995). "Long-range order in a two-dimensional dynamical XY model: how birds fly together." Physical Review Letters 75: 4326-4329. (Murray 没引用, 但 22 年后(2026) "鸟类群体" 的自驱动模型, 与 Murray 12.1 节的 kinematic wave "无物质传输的波" 概念同源。)

[X21] Vicsek, T., Czirók, A., Ben-Jacob, E., Cohen, I., Shochet, O. (1995). "Novel type of phase transition in a system of self-driven particles." Physical Review Letters 75: 1226-1229. (Murray 没引用, 但 22 年后(2026) "自驱动系统" 的开创工作, 与 Murray 12.1 节的"集体行为无物质传输" 同源。)

[X22] Skinner, F.K., Kopell, N., Marder, E. (1994). "Mechanisms for oscillation and frequency control in reciprocally inhibitory model neural networks." Journal of Computational Neuroscience 1: 69-87. (Murray 没引用, 但 22 年后(2026) 神经振子 CPG 的电导基模型, 是 Cohen 1982 的更生物学现实版本。)