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

12.1 Belousov–Zhabotinskii 反应中的运动学波（Kinematic Waves in the Belousov–Zhabotinskii Reaction）

本章讨论两类与第 11 章所介绍的扩散/趋化/对流机制"看起来像波"但本质上完全不同的波现象。共同的特征是空间图样来自"生物振荡器"的相位或周期在空间上的变化，而不是任何物质的实际输运。这类波被称作运动学波（kinematic waves），也可称为"伪波"（pseudowaves）。本章两节分别讨论 Belousov–Zhabotinskii（BZ）反应的水平条带（12.1）和七鳃鳗/狗鱼游动的神经控制（12.2–12.4），所发展的数学工具适用于一般的空间分布振荡器。

BZ 反应中当催化剂采用铁（Fe²⁺↔Fe³⁺）时，反应液在红色（橙红）与蓝色之间发生振荡，加入合适染料可清晰显示。在垂直放置、未搅拌的反应管中会形成水平方向的红蓝相间条带，这些条带通常从底部开始向上移动，且后续条带的移动速度逐渐变慢；最终管内充满条带，且越靠近底部条带越密集。整个图样的生成和传播过程中扩散不起作用，这与第 11 章及以后将讨论的反应-扩散波有本质区别。Beck 和 Váradi（1972）给出了这种图样的运动学解释，本节沿用 Kopell 与 Howard（1973）的分析方法。

物理直觉：一排挂在同一水平杆上的单摆，摆长略有梯度（即周期略有梯度），如果同时开始摆动，很快就会看到沿摆列传播的"波"，且波长越来越短。把管内每个位置看作独立的振荡器，周期 T 可以是位置的函数；如果各振荡器相位不同或频率不同，仅由相位/频率的空间变化就会产生空间图样。实验上，BZ 管内产生条带的原因是某种化学浓度梯度（Kopell-Howard 用硫酸浓度梯度）或温度梯度（Thoenes 1973 用温度梯度）使得振荡周期 T(z) 沿轴单调变化。在管中插入任何化学物质都无法通过的屏障既不影响"伪波"的传播也不影响条带的密度：空间输运过程在该图样的生成中完全不起作用。

把管底设为 z=0、管顶归一化为 z=1。设 T(z) 在 0 ≤ z ≤ 1 上有定义；用 2π-周期函数 φ(z,t) 表示振荡器的相位（对 BZ 反应而言，相位前沿即蓝前沿）。设初始相位分布为 φ₀(z)，令 ψ(z,t) 满足 φ(z,t) = ψ(z,t) + φ₀(z)，且 ψ(z,t+T(z)) 比 ψ(z,t) 多 2π。定义参考相位点 φ=0，定义 t(z) 为位置 z 处相位重新为零的时刻，即 0 = ψ(z,t(z)) + φ₀(z)。则在 t = t(z) + nT(z) 时 φ = 2nπ，故相位 2nπ 在 (z,t) 平面上沿曲线 t = t(z) + nT(z) 移动。

为看清实质，取 φ₀(z) = 0（即初始全同相位），此时 t*(z) = 0；并取 T(z) 单调递增的简单情形。定义 t_n(z) = nT(z)，则 φ(z,t_n(z)) = 2nπ，即 t_n(z) 是第 n 个波前经过 z 的时刻。第 n 个波前的速度为 v_n(z) = |dt_n/dz|⁻¹ = 1/(nT'(z))。由于 T'(z) > 0，第 n 个波前在 t = nT(0) 从 z=0 出发，以第 1 个波速度 1/n 倍的速度向上移动，并在 t = nT(1) 到达 z=1。因为速度按 1/n 递减，相同时间间隔内从 z=0 进入的波数多于从 z=1 离开的波数，管内波数随时间增多。

取具体例子 T(z) = 1+z、φ₀(z) = 0，则 φ(z,t) = 2πt/(1+z)，第 n 个波的速度 v_n = 1/n：在 t=1 时 φ=2π 的波从 z=0 进入，速度 v_1 = 1，t=1 走完全管；t=2 时 φ=4π 的波以 v_2 = 1/2 进入，t=2 走完全管；t=3 时 φ=6π 的波以 v_3 = 1/3 进入，t=3 走完全管……t 越大，管内波数越多，且越靠近 z=0 越密集。实验图（图 12.1(b)）与之定性吻合：约 7 分钟内管内条带密度自下而上递减。若初始相位 φ₀(z) ≠ 0，则 t*(z) ≠ 0，对大 t 仍有 T(z) ∼ t/n，上述分析在 t → ∞ 时仍渐近成立；实际实验装置中 T(z) 不会严格线性，因此实验图样正是这种渐近行为。

生物学含义：生物振荡器与生物钟普遍存在，时变的空间图样可能仅由振荡器参数的空间变化引起，而不是反应-扩散或其他图样生成机制的必然结果。本节图样是连续变化的，下两节讨论一种更"相干"的波图样生成器，即中枢模式发生器（central pattern generator, CPG）。

12.2 中枢模式发生器：鱼类游动的实验事实（Central Pattern Generator: Experimental Facts in the Swimming of Fish）

鱼类游动由沿身体从头到尾的行波驱动，速度取决于波频；沿脊柱分布的神经网络控制肌肉运动并协调行波。一般直觉认为节律控制中心在脑，但许多动物在脊髓与脑断离（脊髓横切，spinal transection）后仍能游动，狗鱼即自 19 世纪末已被记录到这一现象。横切后所观察的游动仍具有正常的节段间肌肉协调。所需节律行为及其节段间协调的神经基础是脊髓中的一组中枢神经网络，能生成时间序列信号以驱动周期性肌肉活动。这类网络被称为中枢模式发生器（CPG），其定义特征是无需任何外部输入即可产生节律输出。

对高等脊椎动物可能有上百万神经元参与，研究必然从尽可能小但仍能呈现横切后活动的脊髓入手；这种神经活动被称为"虚拟游动（fictive swimming / fictive locomotion）"，即使产生实际游动的肌肉被移除，脊髓输出的神经信号与完整游动鱼相同。Grillner（1974）、Grillner 与 Kashin（1976）、Grillner 与 Wallén（1982）提供了狗鱼的实验数据；Kopell（1988）以这些工作作为其理论的案例研究。Cohen 与 Wallén（1980）（另见 Cohen 与 Harris-Warrick 1984）使用七鳃鳗（lamprey）进行了关键实验：七鳃鳗是较为原始的脊椎动物，体长约 13–30 cm，脊髓约 0.3 mm 厚、1.5 mm 宽，神经元数目较少但具有基本的脊椎动物组织，并表现虚拟游动行为。实验取约 25–50 节段长度的离体脊髓（七鳃鳗约有 100 节），动物被去头、暴露脊髓但保留大部分肌肉。

实验中脊髓浸于生理盐水，用 L-DOPA 或 D-谷氨酸诱发虚拟游动，可持续数小时；电极放在同一节段的左、右腹根（ventral root, VR）上记录运动神经元活动。电极显示同一节段左右 VR 之间的脉冲呈交替爆发（bursts），即同一节段的左右两侧是相位相差 180° 的两个振荡器，且这种节段内耦合非常稳定。取自 27 节段离体脊髓第 7 与第 19 节段的记录（图 12.3）显示爆发周期约 1 s/次，并且这两节段之间存在近似恒定的相位滞后（无论左还是右 VR 都近似相同）——这是节段间协调的关键定量特征。仅约 10 节段长度的脊髓片即可产生稳定的虚拟游动输出。

Cohen 与 Wallén（1980）、Cohen 与 Harris-Warrick（1984）的结果与游动鱼的波动运动直接对应（图 12.4）；这种周期性活动由 CPG 控制，本章即对此建模。第 12.3 与 12.4 节详细描述 Cohen 等（1982）提出的模型，它已被成功用于解释脊髓选择性外科损伤相关的实验结果（Cohen 与 Harris-Warrick 1984）。

12.3 中枢模式发生器的数学模型（Mathematical Model for the Central Pattern Generator）

第 12.2 节实验（图 12.3）揭示的关键特征：同一节段的左右 VR 是相位锁定 180° 的两个振荡器；从头到尾存在近似恒定的相位滞后。模型的假设：(i) 每个脊髓节段关联一对神经元振荡器，每个孤立时呈稳定的极限环周期振荡，振幅只取决于内部参数，通常不受药物或电刺激等外部因素影响；(ii) 每个振荡器与其最近邻耦合，但允许长程耦合（Buchanan 与 Cohen 1982 提供实验证据）。

第 7 章曾指出许多生化反应系统可呈现稳定的极限环周期振荡，此处假设节段振荡器的驱动力正是这样一类生物振荡器（或耦合到某种神经元电学性质）。无需知道具体的生物振荡器机制（事实上并不知道），也无需 f 的具体函数形式。用极限环作为局部坐标系，把 γ 上的周期极限环看作相位 θ 从 0 变到 2π 完成一圈闭合轨道；假设电极所观察到的爆发发生在 θ = 0 的某个点 P，并且坐标系选取得使沿 γ 的速率 dθ/dt 为常数。

用一个刻意构造的教学例子（系统 (12.10)–(12.12)）展示该思想：在极坐标 (ρ,θ) 下，系统变为 dρ/dt = ρ(1-ρ)，dθ/dt = ω，从而极限环 ρ = 1 渐近稳定（任何 ρ≠1 的扰动都会被拉回 ρ=1，且因 dθ/dt > 0 逆时针行进）；这里 ρ=1 等价于 γ，θ 即为 0 到 2π 的相位。回到模型系统 (12.9)，取 θ 为 n 个变量之一（值取模 2π），其余 n−1 个变量记为 r，表示对极限环轨道 γ 的垂直扰动；未受扰时 r=0。把系统 (12.9) 重写为 dr/dt = f₁(r,θ)，dθ/dt = ω + f₂(r,θ)，其中 f₁(0,θ)=0=f₂(0,θ)，振荡器周期 T = 2π/ω，且 f₁、f₂ 关于 θ 是 2π 周期的。无外部激发时极限环为 r=0、θ(t) = θ(0) + ωt。

每个节段有两个耦合的振荡器，在无节段间影响时，左右输出 x_R(t)、x_L(t) 简单地以 180° 相位振荡，每个遵循 (12.9) 形式，即 d x_R/dt = f_R(x_R, x_L)，d x_L/dt = f_L(x_L, x_R)。与每个振荡器相关的是相位 θ_R、θ_L 与偏差向量 r_R、r_L（即节段内耦合引起的相对孤立情况的偏差）。实验观察到节段内对是 180° 相位差，故模型假设 θ_L(t) = θ_R(t) + π（这只是一阶近似，因为稍后将加入弱的节段间耦合）；由此 x_L(t) = x_R(t + T/2)，即 (12.16) 可简化为单振荡器方程，回到 (12.13)–(12.14) 形式。

设定单振荡器的相位 θ 增加时，输出变量 x(t) 周期性变化（图 12.7）；设爆发在 x(t) 超过某阈值时开始，并持续到 x(t) 降到该阈值以下；把阈值出现时所对应的相位取为 θ=0。实验观察到爆发仅占周期约 0.4（即 2π × 0.4），设 θ₀ 为爆发停止的相位，则爆发发生在 0 < θ < θ₀ + 2mπ（m=0,1,2,…）区间。

耦合振荡器序列：考虑 N 节段振荡器（每节段一对，但分析只需考虑其中一个），第 j 节段为 dx^j/dt = f^j(x^j) + g^j(x¹,…,x^N,c)，其中 g 表示其他 N−1 振荡器的耦合效应，c 是耦合参数向量；c=0 时退化为各振荡器自身的极限环 γⱼ。系统 (12.19) 一般不可解，需要简化假设或特殊化；耦合振荡器相互作用的研究由来已久，硬性定量结果并非易事，多数研究聚焦于少量振荡器或弱耦合情形，常用奇异扰动方法（第 9 章 9.5、9.10 节）。当耦合弱（即 |c|≪1，|g|≪|f|）时，使用扰动方法处理：意味着孤立振荡器的极限环 γⱼ 仅被耦合轻微扰动，故可继续使用 (12.13)–(12.14) 形式（相位 θⱼ 与对 γⱼ 的偏差 rⱼ），只需加入一个小的耦合项，得到 (12.20)。

实验观察到各孤立振荡器频率 ωⱼ 不同因而周期 Tⱼ = 2π/ωⱼ 不同。关键事实是耦合后振荡器仍执行极限环振荡，所以仍能用相位 θⱼ 刻画；虚拟游动是相位耦合的反映，故只需研究相位耦合模型。把对 (12.19) 在 rⱼ = 0 的扰动约化，得到相位耦合方程 dθⱼ/dt = ωⱼ + hⱼ(θ₁,…,θ_N,c)，j=1,…,N（式 (12.21)），其中 hⱼ 包含其他振荡器的弱耦合效应；该方程不包含振幅。弱耦合下耦合振荡器群体的研究已有相当深度（Neu 1979, 1980；Rand 与 Holmes 1980；Ermentrout 1981；Guckenheimer 与 Holmes 1983）；把 (12.19) 在 rⱼ=0 处扰动最终会导出 (12.21)（见第 9 章）。

考虑线性耦合模型——第 j 个振荡器对第 i 个的作用简单与 xⱼ 成正比——得到 d x_i/dt = f_i(x_i) + Σⱼ≠ᵢ A_{ij} xⱼ（式 (12.22)），A_{ij} 为耦合系数矩阵。(12.22) 同时包含相位与振幅，但论证表明只需考虑相位耦合模型，所以 (12.22) 仍比所需的更复杂。第 9 章 9.9 节（另见 Neu 1979, 1980；Rand 与 Holmes 1980）证明弱耦合下 (12.22) 化为相位耦合方程 dθᵢ/dt = ωᵢ + Σⱼ a_{ij} h(θⱼ−θᵢ)（式 (12.23)），h 是其变量的周期函数；亦可启发式论证：因为只需考虑相位模型，线性耦合应合理地包含相位差的函数；h 的周期性来自任意两个振荡器之间相位差的周期性；Buchanan 与 Cohen（1982）发现节段内运动神经元的慢变势是准正弦的，为该猜想提供实验支持。取简单周期函数 h(φ)=sin φ 作为具体研究对象（更一般的 h 讨论见 Kopell 1988）。

所详细分析的相位耦合模型为 (12.24)：dθᵢ/dt = ωᵢ + Σⱼ≠ᵢ a_{ij} sin(θⱼ−θᵢ)，i=1,…,N。a_{ij} 衡量第 j 振荡器对第 i 振荡器的效应：a_{ij} > 0 时为兴奋性（θⱼ 把 θᵢ 拉向自己），a_{ij} < 0 时为抑制性（θⱼ 加大与 θᵢ 的差异）。取 h(φ)=sin φ 时，最大兴奋效应发生在相位差 π/2，最大抑制在 −π/2，同相时无效应（图 12.8）。建模要求耦合后振荡器仍保持周期特性，且 dθᵢ/dt 始终是 t 的正单调增函数（使极限环 γᵢ 上的运动始终单向，输出呈规则升降），故假设 ωᵢ 与 a_{ij} 的幅值使 θᵢ(t) 满足此行为。

12.4 相位耦合模型系统的分析（Analysis of the Phase Coupled Model System）

需选定 (12.24) 的耦合类型；本节分析最简单情形，即每个振荡器只与最近邻耦合，得到 N 个方程 (12.25)：端点只有一个耦合项 dθ₁/dt = ω₁ + a₁₂ sin(θ₂−θ₁)，dθ_N/dt = ω_N + a_{N,N−1} sin(θ_{N−1}−θ_N)，中间节段 dθⱼ/dt = ωⱼ + a_{j,j−1} sin(θ_{j−1}−θⱼ) + a_{j,j+1} sin(θ_{j+1}−θⱼ)。引入 φⱼ = θⱼ−θ_{j+1}（相位差）与 Ωⱼ = ωⱼ−ω_{j+1}（频率差），两两相减得到 N−1 个相位差方程 (12.27)。寻找与虚拟游动对应的规则周期图样，可对耦合系数 a_{ij} 做假设：设上行耦合（按节段编号方向，即头到尾方向）系数 a_{j,j+1} = a_u，下行系数 a_{j,j−1} = a_d；则 (12.27) 可写成向量形式 dφ/dt = Ω + BS（式 (12.28)），其中 φ = (φ₁,…,φ_{N−1})^T、S = (sin φ₁,…,sin φ_{N−1})^T、Ω = (Ω₁,…,Ω_{N−1})^T，B 是 (N−1)×(N−1) 三对角矩阵 (12.30)。

为把模型应用到七鳃鳗的虚拟游动，需寻找 (12.28) 的相位锁定（phase locked）解：所有振荡器具有同一周期，亦即相位差 φⱼ 恒定（j=1,…,N−1），即 dφⱼ/dt = 0，故求 (12.28) 的平衡解 0 = Ω + BS，即 S = −B⁻¹Ω（式 (12.31)）；因为 S 仅含 sin φⱼ，解存在当且仅当 B⁻¹Ω 的每个分量绝对值 ≤ 1。

两振荡器系统：只剩一个相位差方程 dφ/dt = Ω − (a_d + a_u) sin φ（式 (12.32)）。相位锁定解存在当且仅当 |Ω| ≤ |a_d + a_u|（式 (12.33)）。记解为 φ_S：当 |Ω| = |a_d + a_u| 时 φ_S = π/2 或 3π/2；当 |Ω| < |a_d + a_u| 时 |φ_S| < π/2 或 π/2 < |φ_S| ≤ π（式 (12.34)）。用线性化 dψ/dt ≈ −[(a_d+a_u) cos φ_S] ψ 判定稳定性：(12.34) 的前两个解 cos φ_S = 0 为中性稳定，后两个解中一个稳定、一个不稳定。随着耦合强度 |a_d+a_u| 相对频率差 |Ω| 减小（|Ω|/|a_d+a_u| 增大），稳定与不稳定稳态都趋向中性稳态；当 |a_d+a_u| → |Ω| 时发生分叉，稳定相位锁定解消失。模型中两振荡器对应两对 VR，结合图 12.7 的阈值爆发与每对 VR 的相位差，爆发输出如图 12.9(a)。

三振荡器系统：稳态由 (12.31) 给出矩阵方程 (12.35)；取 a_d = a_u = a 的代数简化情形，得 sin φ₁ = (2Ω₁+Ω₂)/(3a)、sin φ₂ = (Ω₁+2Ω₂)/(3a)（式 (12.36)）。相位锁定解存在当且仅当 max{|2Ω₁+Ω₂|/3a, |Ω₁+2Ω₂|/3a} < 1（式 (12.37)）。一般情形下 φ₁ ≠ φ₂，典型 VR 输出如图 12.9(b)。当振荡器数目增加时代数复杂度迅速失控，但求解过程（得到耦合系数与失谐参数 Ω 之间需满足的条件）的设置方法已经清楚。

N 振荡器恒定相位滞后系统：保留一般 N 振荡器模型的代数处理可能过于繁琐，本节旨在说明耦合振荡器能给出稳定的相位锁定系统，从而可作为 CPG 产生虚拟游动所需 VR 输出的候选。设相邻节段间存在恒定相位差 φⱼ = θⱼ−θ_{j+1} = δ（δ 为正常数），该情形与图 12.3 的实验事实直接对应，这样的振荡器链是均匀行波的典型特征。取 sin δ > 0 代入 (12.31)、(12.29)、(12.30) 后得到稳态频率差方程 (12.38)：Ω₁ = a_d sin δ，Ωⱼ = 0（j=2,…,N−2），Ω_{N−1} = a_u sin δ。回到原频率（因 Ωⱼ = ωⱼ−ω_{j+1}）：中间所有振荡器具有同一自然频率 ω；头部（喙端）振荡器被调高 ω₁ = ω + a_d sin δ > ω；尾部（尾端）振荡器被调低 ω_N = ω − a_u sin δ < ω（式 (12.39)）。该解要求耦合系数 a_d、a_u 为正（即兴奋性耦合），所得行波从头向尾传播。另一解是头部调低、尾部调高，对应行波从尾向头传播（七鳃鳗确实能倒退游动）。

需考察该恒定相位滞后解的稳定性：把 φⱼ = δ + ψⱼ 代入时变方程 (12.28) 并线性化（|ψ|≪1），得到 dψ/dt = B ψ cos δ（式 (12.41)），B 由 (12.30) 给出且只依赖 a_d、a_u，为三对角矩阵。设 ψ(t) = e^{λt} ψ₀，则 λ 是 |B cos δ − λ I| = 0（式 (12.43)）的解。实验观察到的相位滞后 δ 很小，故 0 < cos δ < 1；在此条件下（例如用 Routh-Hurwitz 判据对 λ 的多项式分析）可得 Re λ < 0，由 (12.42) 知 ψ(t) → 0（t → ∞），即该相位锁定解对线性扰动稳定。sin φⱼ = sin δ 还有 2^{N−1} − 2 个其他相位锁定解，但通过特征值矩阵研究可知它们全部线性不稳定（即存在某 λ 使 Re λ > 0）。因此恒定相位差解 (12.39) 是相关解，给出沿脊髓向下（或向上）传播的稳定行波。

需要指出的是：就七鳃鳗而言，离体脊髓片段也能"前游"和"后游"，所以必须假设某种机制自动调节两个末端节段的频率。该简单模型的另一局限是无法解释实验事实——相位滞后在不同游动速度下似乎保持恒定；这是重要的问题，因为相位滞后决定波长，进而决定游动鱼的身形，但其他耦合方式（如 Cohen 等 1982 讨论）可以考虑。第 12.2–12.4 节的目的之一是展示这样一个相对简单的模型如何成为实验观察到的虚拟游动中行波传播的模式生成器；值得强调的是，各种简单的节段间振荡器耦合都能产生稳定行波，即便所分析的简单模型也已足够生成前后游动所需的相位耦合协调——只需重新调谐头尾两个振荡器。虽然仍有问题，结果令人鼓舞；这一切并不证明该模型机制就是 CPG，只说明它是一个可能的候选。

本章个人批注

第 12 章最值得反复咀嚼的核心论点是"运动学波"这个概念。第 11 章建立的直觉是"波 = 某种物质在空间中的输运"，无论是扩散（物质本身）、趋化（细胞主动迁移）还是对流（流体携带）。本章一上来就把这个直觉打碎：在 BZ 反应管里，水平条带的图样看起来完全是波，但实际没有任何东西在管中被输运——每个位置都是一个独立振荡器，位置之间的"波传播"只是因为每个位置达到相位 2nπ 的时刻 t_n(z) = nT(z) 随位置单调变化。这个反直觉的事实让我重新审视了"空间图样"的本质：图样的生成机制未必与图样在时间-空间上的视觉表现对应。

Kopell-Howard 的具体例子 T(z) = 1+z 给我印象极深。公式 v_n = 1/n 简洁得令人发指——第 n 个波的速度严格是第 1 个波的 1/n 倍，这意味着从动力学上看，每个波前都满足同样的几何规律，只是速度按 1/n 衰减。这直接推导出"底部条带更密"这个看似偶然的现象其实是 T(z) 单调递增的必然结果。如果 T(z) 是常数，则所有波速度相同，管内波数恒定；如果 T(z) 单调递减（化学浓度梯度反向），则后到的波反而比先到的快，情形会完全不同。这个简单图像让我意识到"梯度"在振荡器空间分布中扮演的角色远比它在反应-扩散中的角色更直接。

第 12.2–12.4 节则从化学振荡转到神经振荡。CPG 概念的提出让我感受到生物学家与数学家之间方法论的深刻差异。生物学家花了近一个世纪才发现脊髓横切后狗鱼仍能游动，这本身是令人震惊的发现；而数学家把"每个节段一对相位锁定 180° 的振荡器"翻译成 dθⱼ/dt = ωⱼ + Σ a_{ij} sin(θⱼ−θᵢ) 这一简洁的方程组。实验事实中我最看重的是"约 10 节段长度的脊髓片即可产生稳定的虚拟游动输出"——这意味着 CPG 是高度模块化的、可扩展的，每个节段都是相对独立的"组件"，节段间通过弱耦合维持整体协调。这与现代神经科学对模块化神经网络的设计哲学是一致的。

第 12.4 节对相位耦合模型的数学处理有几个让我印象深刻的点：(1) 把 (12.25) 通过 φⱼ = θⱼ−θ_{j+1}、Ωⱼ = ωⱼ−ω_{j+1} 的代换化为关于相位差的方程 (12.27)，这是降维技巧的典范——把绝对相位消去，只保留相对关系；(2) 对两振荡器系统的稳定性分析显示 cos φ_S = 0 给出中性稳定，这对应于分叉点 |Ω| = |a_d+a_u|，是经典的鞍结分叉结构；(3) 对 N 振荡器恒定相位滞后系统的解 (12.39) 显示只有头尾两个振荡器需要被"调谐"，中间所有振荡器都保持自然频率——这正是"边界条件决定全局行为"在振荡器网络中的体现。

我对第 12.4 节末尾的两点坦诚评论印象尤其深刻：(a) 离体脊髓能前后游动，所以"自动调谐头尾"的机制必须存在，但本模型没有给出；(b) 相位滞后与游动速度无关的实验事实，本模型也解释不了。这两点展示了 Murray 作为数学家的科学诚实：他给出了清晰的模型，展示了模型能解释什么，并明确指出了模型不能解释什么，而不是用花言巧语掩盖缺陷。

本章与全书的结构关系上，第 12 章是第 7 章（生物振荡器与开关）和第 9 章（受扰耦合振荡器）的直接应用：第 7 章给出了"什么是生物振荡器"的数学基础，第 9 章给出了"耦合振荡器如何相互作用"的分析工具（弱耦合相位方程 (12.23) 正是第 9 章 9.9 节的核心结果），第 12 章则把这两个工具应用到两个具体的生物现象（BZ 反应运动学波和七鳃鳗虚拟游动）上。这是从抽象到具体的标准数学教材结构。

与上下章的衔接（一段话）

第 11 章"反应-扩散、趋化"主要讨论的是"物质如何在空间中被输运并形成图样"，重点在 Turing 不稳定性、趋化方程、行波解等空间模式形成的机制；本章则把关注点从"物质输运"切换到"振荡器性质在空间上的变化也能产生图样/波"，开辟了一条完全不同的空间图样生成路径。在第 11 章的框架中，"波"意味着某种浓度的传播；本章 12.1 节明确指出 BZ 反应管的条带不是真正的波——没有物质输运，是"伪波"。在第 11 章的图样形成与第 13 章即将讨论的行波之间，本章插入了一个关于"振荡器群体"的中间视角：把第 7 章的生物振荡器与第 9 章的耦合振荡器理论应用到具体生物学现象（BZ 反应、七鳃鳗游动）。本章的"位置"在全书结构上是这样的：第 7–9 章提供生物振荡器及其相互作用的数学工具，第 11 章提供空间图样/输运机制的数学工具，本章则是这两种工具在两类具体生物现象上的应用示范，并为后续章节（特别是第 13 章生物行波、第 14 章分形）提供空间分布振荡器的一般理论框架。