第12章：运动学波（Kinematic Waves）

章节概述

本章探讨了两类由振子产生的波动现象：一类是Belousov-Zhabotinskii（BZ）反应中的运动学波，另一类是鱼类游泳中枢模式发生器（Central Pattern Generator, CPG）中的相位耦合波。核心观点是：空间模式的形成可以仅由振子的相位或频率的空间变化引起，而不需要扩散、对流或趋化性的作用。

12.1 Belousov-Zhabotinskii反应中的运动学波

12.1.1 实验现象

当在垂直圆柱管中进行BZ反应时，若使用铁催化剂并配合适当染料，反应液会呈现红色（含Fe²⁺）和蓝色（含Fe³⁺）之间的振荡色彩变化。在未搅拌的垂直圆柱中，从底部开始会形成水平的蓝色和红色条带。这些条带从底部向上缓慢移动，后续条带移动得越来越慢。最终圆柱被这些条带填充，但密度不均匀——越靠近底部，波的堆积越密集。

关键发现：扩散在条带的形成和传播中可忽略不计，这与后续章节讨论的扩散驱动波完全不同。

12.1.2 理论分析

将垂直圆柱的底部设为 $z = 0$，顶部标准化为 $z = 1$。由于初始浓度梯度，振子周期 $T(z)$ 是位置的函数。

相位表示：定义振荡器的相为 $\varphi(z, t)$，初始相位分布为 $\varphi_0(z)$，则： $$\varphi(z, t) = \psi(z, t) + \varphi_0(z), \quad \psi(z, 0) = 0 \qquad(12.1)$$

其中 $\psi(z, t)$ 是一个函数，当时间增加一个周期 $T(z)$ 时增加 $2\pi$： $$\varphi(z, t + T(z)) = 2\pi + \psi(z, t) + \varphi_0(z)$$

波前到达时间：定义 $t^*(z)$ 为位置 $z$ 处相为零的时间，满足： $$0 = \varphi(z, t^*(z)) = \psi(z, t^*(z)) + \varphi_0(z) \qquad(12.2)$$

对于任意整数 $n$，在时间 $t = t^*(z) + nT(z)$ 有： $$\varphi(z, t^* + nT) = 2n\pi \qquad(12.3)$$

因此，在 $(z, t)$ 平面中，相位 $2n\pi$ 对应的点满足： $$t = t^*(z) + nT(z) \qquad(12.4)$$

12.1.3 特定示例：$T(z) = 1 + z$，$\varphi_0(z) = 0$

令 $t_n(z) = nT(z) = n(1 + z)$，则： $$\varphi(z, t_n(z)) = 2n\pi \qquad(12.6)$$

第 $n$ 个波前的速度： $$v_n(z) = \frac{dz}{dt}\bigg|_{\varphi=2n\pi} = \frac{1}{d t_n(z)/dz} = \frac{1}{nT'(z)} = \frac{1}{n} \cdot \frac{1}{T(z)} \qquad(12.7)$$

取 $T(z) = 1 + z$（单调递增函数），则 $T'(z) = 1$，$v_n(z) = 1/n$。

物理意义： - 第一个波前（$n=1$）以速度 $v_1 = 1$ 传播 - 第二个波前（$n=2$）以速度 $v_2 = 1/2$ 传播 - 第 $n$ 个波前以 $1/n$ 的速度传播

由于速度随 $n$ 减小，越来越多的波被限制在 $0 \leq z \leq 1$ 区间内，导致波在底部堆积。

12.1.4 生物学意义

这一分析表明，生物振子和生物钟可能导致随时间变化的空间模式，而不必假设存在反应扩散机制。空间模式可以简单地由振子参数的空间变化引起。

12.2 中枢模式发生器：鱼类游泳的实验事实

12.2.1 鱼的游泳机制

鱼通过从头部向尾部传播的行波推动自身在水中前进，其速度是波频率的函数。控制肌肉运动产生实际波浪并协调它们的是沿背部排列的神经元网络。

12.2.2 脊髓的概念

在哺乳动物中，通常认为节律性波动的控制中心位于大脑。然而，在许多动物中，脊髓横断后游泳仍然能够发生。在狗鱼身上，这一现象自19世纪末就已为人所知。脊髓横断后观察到的游泳运动表现出正确的节段间肌肉协调。

中枢模式发生器（CPG）：脊髓中能够产生时间信号序列的神经网络，不需要外部输入控制即可产生所需的节律性输出。

12.2.3 七鳃鳗（Lamprey）实验

实验动物选择：七鳃鳗是一种较原始的脊椎动物，具有相对较少的细胞（约100个节段），但仍具有基本的脊椎动物组织结构，并能表现出拟游泳行为。

实验设置： - 使用25-50个节段的脊髓片段 - 在生理盐溶液中培养 - 通过L-DOPA或D-谷氨酸化学诱导拟游泳 - 通过放置在单个节段两侧腹根上的电极监测运动神经元活动

关键实验发现： 1. 单个节段的左、右腹根（VR）放电呈180°反相，即相位锁定 2. 不同节段之间存在近似恒定的相位滞后 3. 仅约10个节段的脊髓片段就能产生稳定的神经拟游泳输出

12.2.4 重要的生物学特征

节段内耦合（Intrasegmental coupling）：同一节段的左右VR的振荡器180°反相锁定，这种耦合非常稳定
节段间耦合（Intersegmental coupling）：相邻节段之间存在近似恒定的相位滞后

12.3 中枢模式发生器的数学模型

12.3.1 模型基本假设

根据实验观察，模型的核心假设为： 1. 每个节段关联一对神经元振子，每个振子都表现出稳定的极限环周期振荡 2. 振子的振幅仅依赖于内部参数，通常不受药物或电刺激等外部因素影响 3. 每个振子与其最近邻耦合，也可能存在远程耦合

12.3.2 单个振子及其表示

设极限环变量的向量为： $$\mathbf{x}(t) = (x_1(t), x_2(t), \ldots, x_n(t)) \qquad(12.8)$$

极限环由以下向量微分方程控制： $$\frac{d\mathbf{x}}{dt} = \mathbf{f}(\mathbf{x}) \qquad(12.9)$$

相表示：将周期极限环用相 $\theta$ 参数化，$\theta$ 从 $0$ 到 $2\pi$ 完成一个完整回路。选择坐标系使得沿极限环的解的速度 $d\theta/dt$ 为常数。

简单示例——二维系统： $$\frac{dx_1}{dt} = x_1(1 - \rho) - \omega x_2$$ $$\frac{dx_2}{dt} = x_2(1 - \rho) + \omega x_1 \qquad(12.10)$$ $$\rho = (x_1^2 + x_2^2)^{1/2}$$

转换为极坐标： $$\frac{d\rho}{dt} = \rho(1 - \rho)$$ $$\frac{d\theta}{dt} = \omega \qquad(12.12)$$

极限环为 $\rho = 1$，是渐近稳定的。

12.3.3 一般振子系统的相变量表示

对于一般系统 (12.9)，引入局部坐标 $(r, \theta)$，其中 $r$ 垂直于极限环的扰动，$\theta$ 是相。系统可写为： $$\frac{dr}{dt} = f_1(r, \theta) \qquad(12.13)$$ $$\frac{d\theta}{dt} = \omega + f_2(r, \theta) \qquad(12.14)$$

其中 $f_1(0, \theta) = 0 = f_2(0, \theta)$，周期 $T = 2\pi/\omega$。

无扰动时： $$r = 0, \quad \frac{d\theta}{dt} = \omega \quad \Rightarrow \quad \theta(t) = \theta(0) + \omega t \qquad(12.15)$$

12.3.4 振子对模型

每节段有一对耦合振子，左、右振子输出 $\mathbf{x}_R(t)$ 和 $\mathbf{x}_L(t)$ 满足： $$\frac{d\mathbf{x}_R}{dt} = \mathbf{f}_R(\mathbf{x}_R, \mathbf{x}_L)$$ $$\frac{d\mathbf{x}_L}{dt} = \mathbf{f}_L(\mathbf{x}_L, \mathbf{x}_R) \qquad(12.16)$$

实验观察表明左右振子180°反相： $$\theta_L(t) = \theta_R(t) + \pi \qquad(12.17)$$

这意味着： $$\mathbf{x}_L(t) = \mathbf{x}_R(t + T/2) \qquad(12.18)$$

12.3.5 耦合振子序列

考虑由 $N$ 个振子组成的系统，每个振子满足： $$\frac{dx_j}{dt} = f_j(x_j) + g_j(x_1, \ldots, x_N, \mathbf{c}), \quad j = 1, \ldots, N \qquad(12.19)$$

其中 $\mathbf{g}$ 表示其他振子的耦合效应，$\mathbf{c}$ 是耦合参数向量。

弱耦合假设：当 $|\mathbf{c}| \ll 1$ 且 $|\mathbf{g}| \ll |\mathbf{f}|$ 时，使用微扰方法，可将系统简化为仅涉及相位的方程。

12.3.6 相位耦合模型

通过微扰分析，最终得到相位耦合系统： $$\frac{d\theta_i}{dt} = \omega_i + \sum_{j=1}^{N} a_{ij} h(\theta_j - \theta_i), \quad i = 1, \ldots, N \qquad(12.23)$$

对于线性耦合，具体形式为： $$\frac{dx_i}{dt} = f_i(x_i) + \sum_{j=1, j \neq i}^{N} A_{ij} x_j \qquad(12.22)$$

取 $h(\varphi) = \sin\varphi$，得到最终模型： $$\frac{d\theta_i}{dt} = \omega_i + \sum_{j=1, j \neq i}^{N} a_{ij} \sin(\theta_j - \theta_i), \quad i = 1, \ldots, N \qquad(12.24)$$

耦合系数的意义： - $a_{ij} > 0$：兴奋性耦合，$\theta_j$ 倾向于将 $\theta_i$ 拉向其值 - $a_{ij} < 0$：抑制性耦合，$\theta_j$ 倾向于增加与 $\theta_i$ 的差异

12.4 相位耦合模型系统的分析

12.4.1 最近邻耦合模型

假设每个振子仅与最近邻耦合，方程 (12.24) 变为： $$\frac{d\theta_1}{dt} = \omega_1 + a_{12} \sin(\theta_2 - \theta_1)$$ $$\frac{d\theta_2}{dt} = \omega_2 + a_{21} \sin(\theta_1 - \theta_2) + a_{23} \sin(\theta_3 - \theta_2)$$ $$\vdots$$ $$\frac{d\theta_j}{dt} = \omega_j + a_{j,j-1} \sin(\theta_{j-1} - \theta_j) + a_{j,j+1} \sin(\theta_{j+1} - \theta_j)$$ $$\vdots$$ $$\frac{d\theta_N}{dt} = \omega_N + a_{N,N-1} \sin(\theta_{N-1} - \theta_N) \qquad(12.25)$$

12.4.2 相位差变量

引入相位差 $\varphi_j = \theta_j - \theta_{j+1}$ 和频率差 $\Delta_j = \omega_j - \omega_{j+1}$： $$\varphi_j = \theta_j - \theta_{j+1}, \quad \Delta_j = \omega_j - \omega_{j+1} \qquad(12.26)$$

通过两两相减，得到 $N-1$ 个方程： $$\frac{d\varphi_1}{dt} = \Delta_1 - (a_{12} + a_{21}) \sin\varphi_1 + a_{23} \sin\varphi_2$$ $$\frac{d\varphi_2}{dt} = \Delta_2 + a_{21} \sin\varphi_1 - (a_{23} + a_{32}) \sin\varphi_2 + a_{34} \sin\varphi_3$$ $$\vdots$$ $$\frac{d\varphi_j}{dt} = \Delta_j + a_{j,j-1} \sin\varphi_{j-1} - (a_{j,j+1} + a_{j+1,j}) \sin\varphi_j + a_{j+1,j+2} \sin\varphi_{j+1}$$ $$\vdots$$ $$\frac{d\varphi_{N-1}}{dt} = \Delta_{N-1} + a_{N-1,N-2} \sin\varphi_{N-2} - (a_{N-1,N} + a_{N,N-1}) \sin\varphi_{N-1} \qquad(12.27)$$

向量形式为： $$\frac{d\varphi}{dt} = \Delta + B\mathbf{S} \qquad(12.28)$$

其中 $\mathbf{S} = (\sin\varphi_1, \ldots, \sin\varphi_{N-1})^T$，$B$ 为 $(N-1) \times (N-1)$ 矩阵。

12.4.3 恒定耦合系数

令所有向上（从头到尾）耦合系数 $a_{j,j+1} = a_u$，所有向下耦合系数 $a_{j,j-1} = a_d$，则 $B$ 矩阵为： $$B = \begin{pmatrix} -(a_d + a_u) & a_u & \cdot & \cdot \\ a_d & -(a_d + a_u) & a_u & \cdot \\ \cdot & \cdot & \cdot & \cdot \\ \cdot & \cdot & a_d & -(a_d + a_u) \end{pmatrix} \qquad(12.30)$$

12.4.4 相位锁定条件

对于拟游泳，我们需要所有振子具有相同周期，即所有振子间的相位差 $\varphi_j$ 恒定。这意味着 $\frac{d\varphi_j}{dt} = 0$。

平衡解满足： $$0 = \Delta + B\mathbf{S} \quad \Rightarrow \quad \mathbf{S} = -B^{-1}\Delta \qquad(12.31)$$

相位锁定解存在当且仅当 $B^{-1}\Delta$ 的所有元素介于 $\pm 1$ 之间。

12.4.5 两振子系统

对于 $N=2$，单相方程为： $$\frac{d\varphi}{dt} = \Delta - (a_d + a_u) \sin\varphi \qquad(12.32)$$

相位锁定解（$\frac{d\varphi}{dt} = 0$）存在当且仅当： $$|\Delta| \leq |a_d + a_u| \qquad(12.33)$$

解的稳定性通过线性化分析确定。

12.4.6 三振子系统

平衡态满足（设 $a_d = a_u = a$）： $$\sin\varphi_1 = (2\Delta_1 + \Delta_2)/3a, \quad \sin\varphi_2 = (\Delta_1 + 2\Delta_2)/3a \qquad(12.36)$$

相位锁定解存在条件： $$\max\{|2\Delta_1 + \Delta_2|/3a, |\Delta_1 + 2\Delta_2|/3a\} < 1 \qquad(12.37)$$

12.4.7 常相位滞后系统（$N$ 振子）

假设相邻振子间相位差恒定：$\varphi_j = \theta_j - \theta_{j+1} = \delta$（正常数）。

稳态方程给出： $$\omega_j = \omega, \quad j = 2, \ldots, N-1$$ $$\omega_1 = \omega + a_d \sin\delta > \omega$$ $$\omega_N = \omega - a_u \sin\delta < \omega \qquad(12.39)$$

物理解释：除头尾振子外，所有振子保持其孤立时的自然频率。头部振子被调高频率，尾部振子被调低频率。

波的传播方向： - 当 $a_d, a_u > 0$（兴奋性耦合）时，波从头向尾传播 - 若头部振子调低、尾部振子调高，则波从尾向头传播（七鳃鳗可以倒游）

12.4.8 稳定性分析

设 $\varphi_j = \delta + \psi_j$，线性化后： $$\frac{d\psi}{dt} = B\psi \cos\delta \qquad(12.41)$$

特征值方程为： $$|B\cos\delta - \lambda I| = 0 \qquad(12.43)$$

可以证明（使用Routh-Hurwitz条件），当 $0 < \cos\delta < 1$ 时，$\text{Re}\lambda < 0$，故常相位滞后解是线性稳定的。

12.4.9 模型局限性

孤立脊髓片段可以向前和向后"游泳"，需要假设某种机制自动调节末端节段
模型似乎无法解释实验观察到的相位滞后不随游泳速度变化的现象
存在其他可能的耦合方式值得考虑

公式汇总表

编号	公式	描述
(12.1)	$\varphi(z, t) = \psi(z, t) + \varphi_0(z)$	相位的空间-时间表示
(12.4)	$t = t^*(z) + nT(z)$	相位 $2n\pi$ 的波前轨迹
(12.5)	$t_n(z) = nT(z)$	$n$ 个波前通过位置 $z$ 的时间
(12.7)	$v_n(z) = \frac{dz}{dt}\\|_{\varphi=2n\pi} = \frac{1}{nT(z)}$	第 $n$ 个波前的速度
(12.8)	$\mathbf{x}(t) = (x_1(t), x_2(t), \ldots, x_n(t))$	极限环变量向量
(12.9)	$\frac{d\mathbf{x}}{dt} = \mathbf{f}(\mathbf{x})$	极限环微分方程
(12.12)	$\frac{d\rho}{dt} = \rho(1-\rho), \quad \frac{d\theta}{dt} = \omega$	极坐标下的极限环系统
(12.13)-(12.14)	$\frac{dr}{dt} = f_1(r,\theta), \quad \frac{d\theta}{dt} = \omega + f_2(r,\theta)$	一般振子系统的相变量表示
(12.17)	$\theta_L(t) = \theta_R(t) + \pi$	左右振子180°反相
(12.19)	$\frac{dx_j}{dt} = f_j(x_j) + g_j(x_1,\ldots,x_N,\mathbf{c})$	耦合振子序列一般形式
(12.23)	$\frac{d\theta_i}{dt} = \omega_i + \sum_{j=1}^N a_{ij}h(\theta_j - \theta_i)$	一般相位耦合方程
(12.24)	$\frac{d\theta_i}{dt} = \omega_i + \sum_{j=1, j\neq i}^N a_{ij}\sin(\theta_j - \theta_i)$	正弦相位耦合模型
(12.25)	最近邻耦合振子系统方程	最近邻耦合形式
(12.26)	$\varphi_j = \theta_j - \theta_{j+1}, \quad \Delta_j = \omega_j - \omega_{j+1}$	相位差和频率差定义
(12.28)	$\frac{d\varphi}{dt} = \Delta + B\mathbf{S}$	相位差的向量形式
(12.30)	$B$ 矩阵（带状结构）	耦合系数矩阵
(12.31)	$\mathbf{S} = -B^{-1}\Delta$	相位锁定平衡条件
(12.32)	$\frac{d\varphi}{dt} = \Delta - (a_d + a_u)\sin\varphi$	两振子系统的相位方程
(12.33)	$	\Delta
(12.36)	$\sin\varphi_1 = (2\Delta_1 + \Delta_2)/3a, \quad \sin\varphi_2 = (\Delta_1 + 2\Delta_2)/3a$	三振子系统（$a_d = a_u = a$）
(12.37)	$\max{	2\Delta_1 + \Delta_2
(12.39)	$\omega_1 = \omega + a_d\sin\delta, \quad \omega_N = \omega - a_u\sin\delta$	常相位滞后系统的频率调节
(12.41)	$\frac{d\psi}{dt} = B\psi\cos\delta$	线性化稳定性方程
(12.43)	$	B\cos\delta - \lambda I

关键概念总结

运动学波 vs 反应扩散波

特征	运动学波	反应扩散波
产生机制	振子相位/频率的空间变化	化学物质扩散与反应
传播速度	多波长，速度随序数变化	由扩散系数决定
传输介质	无实际物质传输（伪波）	实际物质浓度波
扩散作用	可忽略	关键作用

中枢模式发生器模型要点

极限环振子：每个节段包含一对180°反相的极限环振子
弱耦合假设：节段间耦合足够弱，可用微扰方法处理
相位模型：振幅动力学被忽略，仅考虑相位演化
正弦耦合：$h(\varphi) = \sin\varphi$ 作为相互作用函数
相位锁定：所有振子必须以相同频率振荡
行波：相邻振子间的恒定相位差产生行波

练习提示（第12章）

四振子系统：验证 $N$ 振子系统中当频率差 $\varepsilon \leq a/2$ 时相位锁定解存在
一般相互作用函数：对于一般 $h(\delta)$ 推导等效于 (12.38) 的稳态频率差方程

参考来源：Murray, J.D. Mathematical Biology I: An Introduction, 第12章 Kinematic Waves