Chapter 1: Multi-Species Waves and Practical Applications
《数学生物学 II》J.D. Murray 读书笔记
1.1 直观预期 (Intuitive Expectations)
概述
本章开篇介绍反应扩散系统中行波（travelling waves）的广泛应用背景。Murray指出，多物种反应扩散系统的空间波现象远比单物种模型丰富得多。书中以三个经典的生物学实例引入：
-
Belousov-Zhabotinskii (BZ) 反应：可观察到靶向图案（target patterns）和螺旋波（spiral waves），分别如图1.1(a)和(b)所示。靶向波由起搏点（pacemaker）发出同心圆波；螺旋波则可通过轻轻搅拌试剂人为制造，旋转周期约2分钟。
-
盘基网柄菌（Dictyostelium discoideum）：这种黏菌在发育过程中会发出周期性的cAMP化学信号（化学趋向剂），某些起搏细胞可引发靶向波或螺旋波。细胞移动时明亮，静止时暗淡，形成明暗相间的条纹。
-
神经冲动传播：FitzHugh-Nagumo模型描述的动作电位传播也是行波的一个例子。
反应扩散系统的基本框架
多物种反应扩散系统的基本形式为：
其中 \(\mathbf{u} = (u_1, u_2, \ldots, u_n)^T\) 是物种浓度向量，\(\mathbf{f}(\mathbf{u})\) 是反应动力学项，\(D\) 是扩散系数矩阵。
核心概念
- 起搏子（Pacemaker）：持续产生振荡输出的区域
- 靶向波：从中心点向外传播的同心圆波
- 螺旋波：具有旋转核心的行波，在二维空间中表现为螺旋图案
- 行波脉冲（Pulse wave）：代表从稳态出发再回到稳态的单次 excursion
- 行波列（Wave train）：连续的周期性波形
三维螺旋波
在三维空间中，螺旋波具有惊人的拓扑结构——二维螺旋的每一部分本身又是一个螺旋。Winfree (1974) 和 Welsh et al. (1983) 提供了三维波的照片。Mimura及其同事也对球面波进行了大量解析和数值研究。
医学重要性
螺旋波在医学上具有重要意义，特别是在心脏病学（心律失常、房扑、室颤）和神经生物学中。Winfree的工作表明，心脏猝死的前兆可能是旋转的电生理波。
1.2 捕食者-猎物系统中的追逃波 (Waves of Pursuit and Evasion in Predator-Prey Systems)
引言
当捕食者和猎物呈空间分布时，由于捕食者追逐猎物、猎物逃避捕食者，必然产生时空种群波动。海洋浮游生物、动物迁徙、真菌和植被等都能观察到行波带。
模型建立
考虑一个修正的Lotka-Volterra系统，包含猎物的逻辑斯谛增长和捕食者的扩散：
其中 \(U\) 是猎物，\(V\) 是捕食者，\(A, B, C, D, K\) 为正常数，\(D_1, D_2\) 为扩散系数。
无量纲化
令：
得到一维形式：
定常态分析
系统有三个定常态： - \((0, 0)\)：不稳定结点 - \((1, 0)\)：不稳定结点（无捕食者，猎物达环境容纳量） - \((b, 1-b)\)：稳定定常态（共存），条件为 \(b < 1\)
共存定常态 \((b, 1-b)\) 的稳定性取决于 \(4a\) 与 \(b/(1-b)\) 的关系： - \(4a \leq b/(1-b)\) 时为稳定结点 - \(4a > b/(1-b)\) 时为稳定焦点（螺线）
行波解的存在性
设行波解为：
其中 \(c > 0\) 为波速。代入后得常微分方程组：
D = 0 的简化分析
假设猎物扩散远小于捕食者（\(D \to 0\)），将系统化为一阶方程组：
定常态 \((1,0,0)\) 和 \((b,1-b,0)\) 分别为不稳定和稳定的。
线性化分析给出特征方程的根：
行波存在的必要条件：
两类波行为
存在临界值 \(a^*\)，使得： - 当 \(a > a^*\) 时：\((U,V)\) 以振荡阻尼方式趋近 \((b, 1-b)\)（图1.3a） - 当 \(a \leq a^*\) 时：\((U,V)\) 单调趋近 \((b, 1-b)\)（图1.3b）
对流捕食-猎物模型
另一种完全不同的"追逃"模型：猎物的逃跑速度和捕食者的追逐速度与对方的梯度成正比：
展开后得到交叉扩散项：
其中 \(h_1 uv_{xx}\) 和 \(-h_2 vu_{xx}\) 为交叉扩散项。
1.3 灰松鼠在英国空间扩散的竞争模型 (Competition Model for the Spatial Spread of the Grey Squirrel in Britain)
背景
20世纪初，北美灰松鼠（Sciurus carolinensis）被引入英国东南部，此后成功扩散到苏格兰低地，同时本地红松鼠（Sciurus vulgaris）从这些地区消失。Lloyd (1983) 注意到灰松鼠入侵的区域，红松鼠在短短几年重叠期后就会消失。
竞争机制
灰松鼠具有竞争优势的原因： - 在北美原产地，灰松鼠在强烈的种间竞争下进化，偏好过度繁殖 - 红松鼠在无种间竞争的环境中进化，繁殖水平较低 - 灰松鼠体型约为红松鼠的两倍 - 在英国，红松鼠已适应在硬木森林中生活，与灰松鼠产生资源竞争
竞争模型
设 \(S_1(X,T)\) 和 \(S_2(X,T)\) 分别为灰松鼠和红松鼠的种群密度：
竞争条件：\(b_2 > c_1\)，\(c_2 > b_1\)（灰松鼠占优）
无量纲化
令：
得到：
行波解
设行波形式：
边界条件：
特殊情况的解析结果
当 \(\kappa = \alpha = 1\)，\(\gamma_1 + \gamma_2 = 2\) 时，可得：
这是经典的Fisher-Kolmogoroff方程，其最小波速：
对于红松鼠：
量纲最小波速：
参数估计
| 参数 | 含义 | 估计值 |
|---|---|---|
| \(a_1\) | 灰松鼠内禀净增长率 | 0.82/年 |
| \(a_2\) | 红松鼠内禀净增长率 | 0.61/年 |
| \(1/b_1\) | 灰松鼠环境容纳量 | 10只/公顷 |
| \(1/b_2\) | 红松鼠环境容纳量 | 0.75只/公顷 |
扩散系数估计
基于 woodlands 之间扩散的启发式估计：
| \(l\) (km) | \(D_1\) (km²/年) | \(C_{min}\) (km/年) |
|---|---|---|
| 1 | 0.179 | 0.77 |
| 5 | 4.46 | 3.82 |
| 10 | 17.9 | 7.66 |
| 20 | 71.4 | 15.3 |
观测到的扩散速率约为 7.7 km/年，与理论值吻合良好。
1.4 转基因生物的扩散 (Spread of Genetically Engineered Organisms)
引言
重组DNA技术用于改良动植物，但存在生态系统破坏的潜在风险。核心问题是：转基因生物可能以多快的速度和范围扩散？能否通过地理屏障实现控制？
模型框架
使用与灰松鼠相同的Lotka-Volterra竞争模型，但推广到空间异质环境：
\(E\) 为转基因微生物，\(N\) 为自然微生物。
周期性斑块环境
设环境由长度为 \(l_1\) 的适宜斑块（Patch 1）和长度为 \(l_2\) 的不适宜斑块（Patch 2）交替组成，\(l = l_1 + l_2\)：
Patch 1 为适宜环境：\(D_1 \geq D_2\)，\(d_1 \geq d_2\)，\(1 \geq G_2\)，\(1 \geq g_2\)
入侵条件
线性稳定性分析给出特征方程（Hill方程）：
稳定性条件：
包含条件：
三种情况
- \(\gamma_e > 1 > G_2\)：原生定常态在两类斑块中均稳定
- \(1 > G_2 > \gamma_e\)：原生定常态在两类斑块中均不稳定，必然入侵
- \(1 > \gamma_e > G_2\)：原生定常态在适宜斑块不稳定但在不适斑块稳定，结果取决于斑块相对大小
临界斑块长度
临界长度 \(l_1^*\) 决定稳定性边界：
控制策略
当 \(\gamma_e < 1\) 时，确保控制的必要条件：
最简单的控制策略是使 \(\gamma_e > 1\)。
1.5 BZ反应中的行波前 (Travelling Fronts in the Belousov-Zhabotinskii Reaction)
引言
BZ反应是研究化学波的经典范例，可作为生物模式形成的范式，并用于理解心脏等器官的复杂波行为。
反应机理
简化后的反应序列：
\(X\) 为亚溴酸 (\(HBrO_2\))，\(Y\) 为溴离子 (\(Br^-\))，\(A\) 为常量。
反应扩散方程
应用质量守恒定律并包含扩散：
无量纲化
得到：
其中 \(L \approx M \approx O(10^{-4})\)，\(b = O(1)\)，\(r\) 可在5到50之间变化。
前缘波模型
忽略 \(L\) 和 \(M\) 项（\(L, M \ll 1\)），得到前缘波模型：
行波解
设行波形式 \(u(s,t) = f(z)\)，\(v(s,t) = g(z)\)，\(z = s + ct\)，边界条件：
Fisher-Kolmogoroff极限
当 \(v = \frac{1-b}{r}(1-u)\)，\(b \neq 1\) 时，系统约化为Fisher-Kolmogoroff方程：
波速 \(c \geq 2\sqrt{b}\)。从初始条件 \(u(s,0) \sim O(e^{-\beta s})\) 出发：
Murray得到的波速bounds
量纲分析
波速的量纲形式：
典型值：\(c_D \approx 4.5 \times 10^{-3}\) cm/s \(\approx\) 0.27 cm/min，与实验值一致。
关键发现：反应扩散波比纯扩散快几个数量级——传播1 cm约需4分钟，而纯扩散需约850分钟。
1.6 可激发介质中的波 (Waves in Excitable Media)
可激发系统概述
可激发介质最著名的例子是神经细胞通过电信号进行的通讯。FitzHugh-Nagumo (FHN) 模型是 Hodgkin-Huxley 模型的数学简化。
FHN模型
其中 \(0 < a < 1\)，\(b > 0\)，\(\gamma > 0\)。
行波脉冲解
设行波坐标 \(z = x - ct\)，边界条件：
阈值行为
- 若扰动的最大值 \(u_{max} < u_A\)（阈值），解将衰减回原点
- 若 \(u_{max} > u_A\)，则触发大 excursion（ABCDA轨道）
小参数分析
设 \(b = \varepsilon L\)，\(\gamma = \varepsilon M\)，\(\varepsilon \ll 1\)：
前缘波速
忽略 \(O(\varepsilon)\) 项，得到：
这是带三次非线性源的扩散方程，前缘波速：
要求 \(a < 1/2\) 以保证 \(c > 0\)。
尾缘波速
尾缘波速（从 \(C\) 到 \(D\)）：
通过令前缘与尾缘波速相等来确定 \(v_C\)。
一般可激发系统
更一般的快-慢系统：
其中 \(0 < \varepsilon \ll 1\)，\(f\) 和 \(g\) 的零息线呈立方形状。
奇异扰动分析
在 \(O(\varepsilon)\) 尺度下，前缘的 \(O(1)\) 系统为：
这与第13章研究的方程本质相同，存在连接 \(u_S\) 和 \(u_A\) 的单调行波解。
1.7 振荡动力学R-D系统中的行波列 (Travelling Wave Trains in R-D Systems with Oscillatory Kinetics)
背景
对于Hopf分叉产生的极限环动力学，扩散效应会产生行波列。
一般框架
假设空间均匀系统 \(\frac{d\mathbf{u}}{dt} = \mathbf{f}(\mathbf{u}; \gamma)\) 在 \(\gamma < \gamma_c\) 时有稳定定常态，在 \(\gamma > \gamma_c\) 时通过Hopf分叉产生稳定极限环。
行波列形式
其中 \(\sigma > 0\) 是频率，\(k\) 是波数，波长 \(w = 2\pi/k\)，波速 \(c = \sigma/k\)。
代入得常微分方程：
λ-ω 系统
两反应物的λ-ω模型：
其中 \(r^2 = u^2 + v^2\)。
极坐标变换
设 \(u = r\cos\theta\)，\(v = r\sin\theta\)：
若 \(\lambda'(r_0) < 0\) 且 \(\omega(r_0) \neq 0\)，存在极限环 \(r = r_0\)，\(\theta = \theta_0 + \omega(r_0)t\)。
行波列解
设 \(r = \alpha\)，\(\theta = \sigma t - kx\)，得到：
因此：
波速 \(c = \frac{\sigma}{k} = \frac{\omega(\alpha)}{\sqrt{\lambda(\alpha)}}\)。
具体例子
设 \(\omega(r) \equiv 1\)，\(\lambda(r) = \gamma - r^2\)：
小振幅行波列：
Hopf分叉的一般理论
一般两物种系统通过适当的坐标变换可化为λ-ω形式。设：
在Hopf分叉点 \(\gamma = \gamma_c\)，\(Tr M = 0\)，\(det M > 0\)，可选择矩阵 \(N\) 使 \(N^{-1}MN\) 具有标准形式。
这意味着λ-ω系统的分析方法可推广到一般的反应扩散系统。
1.8 螺旋波 (Spiral Waves)
引言
旋转螺旋波广泛存在于生物、生理和化学情境中。BZ反应中的螺旋波是研究最充分的例子之一。
实验观察
- BZ反应：Müller等(1985, 1986)用光吸收技术定量突出显示实际浓度水平
- 脑皮层传播性抑郁波：Shibata和Bures(1974)在鼠脑皮层实验观察到螺旋波
- 心脏组织：Allessie等(1977)在兔心房肌诱导出旋转波
- 黏菌Dictyostelium：Newell (1983)拍摄到令人印象深刻的螺旋信号模式
心脏猝死
Winfree (1994a, 1995) 提出假说：心脏猝死可能涉及电活动的三维旋子（螺旋波），当心脏厚度超过某临界值时突然失稳。心脏纤维化时，小区域独立收缩，导致心室颤动。
数学描述
在极坐标 \((r,\theta)\) 中，简单旋转螺旋由相位的周期函数描述：
其中 \(\Omega\) 是频率，\(m\) 是螺旋臂数，\(\psi(r)\) 描述螺旋类型。
Archimedes螺旋和对数螺旋
带中心核的形式：
螺旋波的数学表达
设 \(w = u + iv\)，螺旋波解可写为：
在极坐标下：
振幅和相位方程：
边界条件
在原点需正则且有界：
远处需满足色散关系：
螺旋波尖端的运动
Yagisita等(1998) 研究了球面上 excitable 系统中的螺旋波，表明螺旋尖端是旋转的。
1.9 λ-ω反应扩散系统的螺旋波解 (Spiral Wave Solutions of λ-ω Reaction Diffusion Systems)
λ-ω系统的螺旋波
对于一般的λ-ω系统：
其中 \(A^2 = u^2 + v^2\)。
化为复数形式：
Koga的研究
Koga (1982) 研究了：
计算了1臂和2臂螺旋波解。
归一化形式
通过变量变换可化为标准形式：
Kuramoto和Koga (1981) 发现当 \(|\beta|\) 较小时，产生稳定旋转螺旋波；增大 \(|\beta|\) 时螺旋波失稳并出现混沌。
一维类比
螺旋波的一维类比是从核心周期性发射的脉冲。方程(1.148)给出：
一个特解为：
公式汇总表
第1章 核心公式一览
| 编号 | 公式名称 | 方程 | 适用场景 |
|---|---|---|---|
| (1.1) | 反应扩散基本方程 | \(\mathbf{u}_t = \mathbf{f}(\mathbf{u}) + D\nabla^2\mathbf{u}\) | 多物种反应扩散系统 |
| (1.2) | 捕食者-猎物模型 | \(\partial_T U = AU(1-U/K) - BUV + D_1\nabla^2 U\) | Lotka-Volterra扩散系统 |
| (1.3) | 无量纲捕食者-猎物 | \(u_t = u(1-u-v) + Du_{xx}\)，\(v_t = av(u-b) + v_{xx}\) | D=0简化分析 |
| (1.4) | 行波变换 | \(z = x + ct\) | 行波求解 |
| (1.9) | 特征根 | \(\lambda_{2,3} = \frac{1}{2}\{c \pm [c^2-4a(1-b)]^{1/2}\}\) | 线性稳定性分析 |
| (1.10) | 捕食者-猎物波速下界 | \(c \geq [4a(1-b)]^{1/2}\) | 波存在的必要条件 |
| (1.12) | 对流模型 | \(u_t - [(c_1+h_1v_x)u]_x = f(u,v)\) | 追逃模型 |
| (1.15) | 竞争模型 | \(\partial_T S_i = D_i\nabla^2 S_i + a_i S_i(1-b_i S_i - c_i S_j)\) | 灰松鼠入侵 |
| (1.18) | 无量纲竞争模型 | \(\theta_{1t} = \nabla^2\theta_1 + \theta_1(1-\theta_1-\gamma_1\theta_2)\) | 竞争行波分析 |
| (1.27) | 灰松鼠最小波速 | \(c_{min} = 2(1-\gamma_1)^{1/2}\) | 入侵波速 |
| (1.30) | 量纲波速 | \(C_{min} = 2\sqrt{a_1 D_1(1-c_1/b_2)}\) | 实际扩散速率估计 |
| (1.31) | 扩散系数估计 | \(D_1 = l^2/(4\times 1.4) = l^2/5.6\) | 斑块间扩散 |
| (1.32) | 转基因入侵模型 | \(E_t = [D(x)E_x]_x + r_E E[G(x)-a_E E-b_E N]\) | 空间异质环境 |
| (1.55) | Hill方程 | \([D(x)f']' + [G(x)-\gamma_e+\lambda]f = 0\) | 入侵稳定性 |
| (1.65) | 入侵包含条件 | \((1-\gamma_e)l_1 \geq (\gamma_e-G_2)l_2\) | 转基因生物控制 |
| (1.71) | 临界斑块长度 | \(l_1^* = \frac{2}{\sqrt{1-\gamma_e}}\arctan\sqrt{\cdots}\) | 稳定性边界 |
| (1.80) | BZ反应机理 | \(A+Y\to X+P\)，\(X+Y\to 2P\)，\(A+X\to 2X\)，\(2X\to P+A\) | 化学振荡 |
| (1.85) | 前缘波模型 | \(u_t = u(1-u-rv) + u_{xx}\)，\(v_t = -buv + v_{xx}\) | BZ反应前缘 |
| (1.88) | Fisher-K波速公式 | \(c = \beta + b/\beta\)（小\(\beta\)）或 \(2\sqrt{b}\)（大\(\beta\)） | 行波速度 |
| (1.93) | BZ波速bounds | \(\frac{r}{2}+\frac{b}{3r}-[2(b+2r)]^{-1/2} \leq c \leq 2\) | Murray的界限 |
| (1.94) | FHN方程 | \(u_t = f(u)-v+Du_{xx}\)，\(v_t = bu-\gamma v\) | 可激发介质 |
| (1.99) | FHN前缘波速 | \(c = \sqrt{D/2}(1-2a)\) | 神经脉冲 |
| (1.101) | 尾缘波速 | \(c = \sqrt{D/2}(u_C-2u_P+u_D)\) | 脉冲尾缘 |
| (1.106) | 振荡R-D系统 | \(u_t = f(u) + u_{xx}\) | 行波列 |
| (1.108) | 行波列 | \(z = \sigma t - kx\) | 周期波 |
| (1.110) | λ-ω系统 | \(w_t = (\lambda+i\omega)w + \nabla^2 w\) | 极限环系统 |
| (1.116) | λ-ω波速关系 | \(\sigma = \omega(\alpha)\)，\(k^2 = \lambda(\alpha)\) | 行波列存在条件 |
| (1.119) | Hopf分叉示例 | \(\omega(r)=1\)，\(\lambda(r)=\gamma-r^2\) | 简单振荡 |
| (1.131) | 螺旋波相位 | \(\phi = \Omega t \pm m\theta + \psi(r)\) | 螺旋波描述 |
| (1.136) | 复数λ-ω方程 | \(w_t = (\lambda+i\omega)w + D\nabla^2 w\) | 螺旋波分析 |
| (1.139) | 螺旋波形式 | \(A=A(r)\)，\(\phi=\Omega t+m\theta+\psi(r)\) | 旋转波 |
| (1.143) | 色散关系 | \(\psi'(\infty)=\sqrt{\lambda(A_\infty)/D}\)，\(\Omega=\omega(A_\infty)\) | 螺旋波远场 |
| (1.147) | 归一化Swift-Hohenberg | $w_t = w-(1+i\beta) | w |
关键概念总结
行波类型
- 行波前（Wavefront）：连接两个不同空间均匀态
- 行波脉冲（Wave pulse）：从定常态出发再回到同一定常态
- 行波列（Wave train）：无限延伸的周期波形
稳定性条件
- Fisher-Kolmogoroff：最小波速 \(c_{min} = 2\sqrt{Df'(0)}\)
- Hopf分叉：\(Tr M = 0\)，\(det M > 0\) 产生极限环
- 螺旋波：需要相位奇点的存在
实际应用
- 生态入侵：灰松鼠取代红松鼠的扩散速率预测
- 基因工程控制：通过斑块设计控制转基因生物扩散
- 心脏病理：螺旋波与心室颤动的关系
- 化学工程：BZ反应作为生物模式形成的范式
分析方法
- 行波变换：将偏微分方程化为常微分方程
- 相平面分析：研究行波解的定性行为
- 奇异扰动理论：处理快-慢系统
- Lyapunov函数：证明全局稳定性
- 特征值分析：确定波速下界
读书备注：本章系统介绍了多物种反应扩散系统中各种行波现象的数学建模与分析方法。从捕食者-猎物追逃波、竞争物种的入侵波、化学振荡反应中的行波列，到螺旋波，Murray展示了反应扩散方程在生物学中的广泛应用。核心数学工具包括行波变换、相平面分析、Hopf分叉理论和奇异扰动方法。这些内容为了解心脏电活动、胚胎发育、种群入侵等复杂时空动力学提供了坚实的理论基础。