跳转至

第十二章：心脏读书笔记

作者

本章作者为 James P. Keener 与 James Sneyd（犹他大学数学系 / 奥克兰大学数学系），合著 Mathematical Physiology, Vol. II: Systems Physiology（Springer IAM 8/II, 2009, 第二版）。本章是第二卷最长的一章（104 页），系统建立了从离子通道 → 单细胞动作电位 → 一维纤维传播 → 二维/三维心肌组织 → 起搏器同步的层次化心脏电生理数学模型。其特点是把 Vol I Ch 5 (Hodgkin-Huxley 兴奋性) 和 Ch 6 (行波传播) 的细胞级工具扩展到临床相关的器官级——ECG、心律失常、除颤、起搏器同步——每一步都给出显式方程、参数表（Table 12.1-12.5）和数值预言。

内容概述

本章按 6 大主题组织：(1) 12.1 心电图（ECG）——从单导联标量 ECG 到 12 导联向量 ECG（VCG），Einthoven 三角 + Burger 三角 + Wilson 中心端；(2) 12.2 心脏细胞——Noble 1962（Purkinje 纤维, 3 离子电流）、MNT 1975（9 离子电流）、YNI 1980（SA 节点, Ca²⁺ 振荡）、Beeler-Reuter 1977（心室肌, 4 电流）、Luo-Rudy 1994（15+ 电流的现代标准模型）；(3) 12.3 细胞耦合——一维纤维的离散/连续电缆方程、gap junction 跳跃条件、传播失败的临界条件；(4) 12.4-12.5 Bidomain 模型与各向异性——心肌组织作为双连续介质、平面波速的各向异性（纵 0.5 m/s vs 横 0.17 m/s）、虚拟电极（virtual electrodes）和 dogbone 现象；(5) 12.6 起搏器同步——SA 节点作为异质振荡器网络，同步与电紧张负荷（electrotonic loading）；(6) 12.7-12.9 心律失常——折返（reentry）、螺旋波、室颤、除颤的虚拟电极机制。

最深的数学内容：(a) 12.2.5 的 Common-Pool vs Local-Control Ca²⁺ 释放模型——Stern 1992 证明 common-pool 模型不能同时有高增益和分级释放，必须引入 diadic cleft（仅 15 nm 宽、含 0.2 个 Ca²⁺ 离子的亚细胞区室）作为局部控制单元；(b) 12.3.3 的 Bidomain 模型 + 12.4.2 的 各向异性平面波速 \(c = \Upsilon(n) c_0 / \sqrt{C_m R_m}\)，其中 \(\Upsilon^2 = R_m(n \cdot \sigma_i n)(n \cdot \sigma_e n) / [\chi n \cdot (\sigma_i + \sigma_e) n]\)；(c) 12.5 的 除颤虚拟电极机制——强电场在 unequal-anisotropy 心肌中产生 dogbone 极化图样，可同时去极化和超极化，从而终止螺旋波。

核心方程与概念

12.1 心电图 — 体积导体理论

电流源 Poisson 方程（12.2）：\(\nabla \cdot (\sigma \nabla \phi) = -S\)，其中 \(S\) 为电流源（心脏动作电位波前处的偶极子层）。前向问题（12.3）\(\phi_B(t) = T \cdot S(t)\) 难解（需要全身电导张量），逆问题（12.4）数值不稳定。单偶极子近似（12.5）\(H(t) = \int_V J \, dV\)（心向量），无限均匀导体下表面电位（12.6）\(\phi(x, t) = H(t) \cdot x / (4\pi |x|^3)\)。

Einthoven 三角（12.8-12.10）：\(V_I = \phi_{LA} - \phi_{RA}\)，\(V_{II} = \phi_{LL} - \phi_{RA}\)，\(V_{III} = \phi_{LL} - \phi_{LA}\)，且 \(V_I + V_{III} = V_{II}\)（KVL）。三个导联向量假设为额面内等边三角形。Burger 三角（实测）更准确：\(l_I = (0.923, -0.298, 0.241)\)，\(l_{II} = (0.202, 0.972, -0.121)\)，并非真正额面。

Wilson 中心端 + 6 个胸导联（V1-V6）提供水平面信息。临床要点：(a) QRS 振幅反映心室肌质量，左偏提示左室肥大；(b) QRS 与 T 波同向（lead I、II、aVR），但二者对应动作电位上升支（去极化）和下降支（复极化），方向相反——"最近激活的组织最先复极"这一 retrograde 现象的机制未完全理解；(c) 12 导联可确定心向量方向，从而定位传导阻滞（LBBB、RBBB）和心肌梗死部位。

12.2 心脏细胞模型

12.2.1 Noble 模型（1962，Purkinje 纤维）

通用形式（12.12）\(I = g(V - V_{eq})\)。Noble 模型（12.13）\(C_m dV/dt + g_{Na}(V - V_{Na}) + (g_{K_1} + g_{K_2})(V - V_K) + g_{an}(V - V_{an}) = I_{app}\)，其中 \(V_{Na} = 40\)、\(V_K = -100\)、\(V_{an} = -60\) mV。时间依赖 K⁺ 电导（12.14）\(g_{K_2} = 1.2 n^4\)（100× 慢于 HH 模型，产生动作电位平台）。瞬时 K⁺ 电导（12.15）经验公式。Na⁺ 电导（12.16）\(g_{Na} = 400 m^3 h + g_i\)，\(g_i = 0.14\) 提供平台期恒定内向电流。门控变量动力学（12.17）\(dw/dt = \alpha_w(1-w) - \beta_w w\)，速率常数形式（12.18）。

问题：Noble 时代还没有 Ca²⁺ 数据，模型把内向 Na⁺ 电流（带 \(g_i\) 偏置）赋予"上升支 + 平台"双重角色。电容 \(C_m = 12\) μF/cm² 异常大（典型 1 μF/cm²），Noble 解释为：实验用的是一小束细胞，有效面积小于总膜面积。

12.2.2 MNT 模型（McAllister-Noble-Tsien 1975）

9 个离子电流 + 9 门控变量。新电流：\(I_{si} = (0.8 df + 0.04 d')(V - V_{si})\)（Ca²⁺ 慢内向，\(V_{si} = 70\) mV）（12.20）。3 个 K⁺ 电流：\(I_{K_2}\)（起搏，pacemaker）、\(I_{x_1}\)、\(I_{x_2}\)（平台期）（12.22-12.25）。Cl⁻ 电流 \(I_{Cl} = 2.5 qr(V - V_{Cl})\)。模型困境：MNT 比 Noble 更复杂但行为相似——这一对比明确显示"模型越复杂 ≠ 越有用"。

12.2.3 YNI 模型（Yanagihara-Noma-Irisawa 1980，SA 节点）

SA 节点的关键差异：没有快速 Na⁺ 电流，phase 0 去极化由慢 Ca²⁺ 电流承担。YNI 模型（12.30-12.40）含 5 个电流：\(I_{Na}\)（快但非 HH 式）、\(I_K\)、\(I_l\)（漏）、\(I_s\)（慢内向）、\(I_h\)（超极化激活的内向电流）。振荡机制：复极化后 \(I_s\) 缓慢去极化至阈值（约 -40 mV）触发新动作电位。简化：因 SA 节点动作电位无初始尖峰，可用 2 变量 FitzHugh-Nagumo 方程定性重现（图 12.13）。

12.2.4 Beeler-Reuter 模型（1977，心室肌）

4 个电流：\(I_{Na}\)（12.41，带再激活变量 \(j\)，因 \(j\) 比 \(h\) 慢得多，HH 单变量 \(h\) 不够）、\(I_K\)（时间无关）（12.42）、\(I_x\)（时间激活外向）（12.43）、\(I_s\)（Ca²⁺ 慢内向）（12.44）。Ca²⁺ 平衡（12.45）\(dc/dt = 0.07(1-c) - I_s\)（\(c = 10^7 [\text{Ca}^{2+}]_i\)）。

12.2.5 Ca²⁺ 释放：Common-Pool vs Local-Control

Stern 1992 关键定理：common-pool 模型（myoplasm 或 diadic cleft 作为 well-mixed 隔室）不能同时具有高增益和分级释放。Local-Control 模型（Stern 1992）假设 L 型通道与 RyR 在 diadic cleft 中"紧邻"，每个 cleft 独立控制，统计招募产生分级响应。

diadic cleft 几何：宽 ~15 nm、半径 ~200 nm、体积 ~\(2 \times 10^{-18}\) L。静息状态 cleft 内平均仅 0.2 个 Ca²⁺ 离子——经典确定性连续模型在 cleft 内失效，需要随机模拟（Greenstein & Winslow 2002、Sobie et al. 2002）。

12.2.6 L 型 Ca²⁺ 通道（Jafri 1998 Markov 模型）

Monod-Wyman-Changeux 类型：两个构象（normal、Ca²⁺），仅 O 状态开放。全模型（12 状态）+ 简化版（6 状态）+ 最简版（3 状态）。激活率 \(\alpha = 2 e^{0.012(V-35)}\)，Ca²⁺ 失活率 \(\gamma = 0.44 / [\text{Ca}^{2+}]_{cleft}\)。

12.2.7 RyR（Stern 1999 模型）

4 状态（Fig. 12.20）：R（关）→ O（开，Ca²⁺ 激活）→ I（失活，Ca²⁺ 诱导）。参数：\(k_1 = 35\) μM⁻² ms⁻¹、\(k_{-1} = 0.06\) ms⁻¹、\(k_2 = 0.5\) μM⁻¹ ms⁻¹、\(k_{-2} = 0.005\) ms⁻¹。争议：RyR 关闭的机制（cleft Ca²⁺ 失活、SR 耗竭、随机关闭、适应）至今无共识。

12.2.8 Na⁺-Ca²⁺ 交换（NCX）

3 Na⁺ 内流 / 1 Ca²⁺ 外排，可逆、电生。模型（12.46）来自 Vol I §2.4.3，Ca²⁺ 外流通量 \(J\) 由 17 个项的复杂分子动力学决定。通量方向（12.47）：当 \([\text{Ca}^{2+}]_i / [\text{Ca}^{2+}]_e > e^{F V_m / RT} ([\text{Na}^+]_i / [\text{Na}^+]_e)^3\) 时 \(J\) 为正（Ca²⁺ 外排）。功能：动作电位起始时 NCX 反转（Ca²⁺ 内流），后期正向工作（Ca²⁺ 外排）。

12.3 细胞耦合 — 1D 电缆方程

带 gap junction 跳跃条件的电缆方程（12.48-12.49）： $\(p\left(C_m \frac{\partial V}{\partial t} + I_{ion}\right) = \frac{\partial}{\partial x}\left(\frac{1}{r_c} \frac{\partial V_i}{\partial x}\right) = -\frac{\partial}{\partial x}\left(\frac{1}{r_e} \frac{\partial V_e}{\partial x}\right)\)$ $\([V_i]/r_g = (1/r_c) \partial V_i/\partial x\)$

特征方程（12.58）：解 \(\mu < 1\) 满足 \(R_g = 2(\mu - 1/E)(\mu - E)/[\mu(E - 1/E)]\)，其中 \(E = e^{\lambda L}\)、\(\lambda^2 = p/(R_m(r_c + r_e))\)。大空间常数极限（12.60）：\(L^2/\lambda^2_g \approx R_g + L^2 p/[R_m(r_c + r_e)]\)，即"电阻沿电缆可加"。

生理学悖论：教科书常说"gap junction 是低电阻"，但 Chapman & Fry 1978 在蛙心室肌的实测：\(L = 131\) μm，\(\lambda_g = 0.328\) cm，得到 48% 电阻来自细胞质、52% 来自 gap junction——gap junction 实际承担总电阻的一半。

传播类型（12.79-12.80）：Fast & Kleber 1993 实测——细胞内 30 μm 需 38 μs，跨细胞 30 μm 需 118 μs（gap junction 延迟 80 μs）。传播是跳跃式（saltatory）的，不是连续的。推导 \(r_g/L = 1.78(r_e + r_c)\)，即 gap junction 占总电阻约 65%。

12.3.2 传播失败（Propagation Failure）

条件（12.86）：当 \(A_1 p_1 / (A_2 p_2) \geq 1 - F(\beta)/F(\alpha)\) 时（\(F\) 是 \(f\) 的原函数），传播方向上 \(pA\) 突增足够大将导致传播阻滞。bundle branch block 机制：左束支分叉处总膜周长 \(p\) 与截面积 \(A\) 同时变化，若满足 (12.86) 条件，AP 不能穿越分叉。

gap junction 介导的阻滞（12.93）：存在临界 \(R_g^* = R_g/(\lambda L) \geq (1-2\alpha)[\text{某函数}]\)，gap junction 电阻大于此值时阻滞发生。关键结论（12.97-12.98）：\(R_g^*\) 随 \(\alpha\) 减小（兴奋性 ↑）而增大，更兴奋的组织更不易阻滞；长细胞比短细胞更不易阻滞。

12.3.3 Bidomain 模型

双连续介质（12.99-12.103）：心肌作为 intracellular + extracellular 两个互穿的连续介质，各有电导张量 \(\sigma_i\)、\(\sigma_e\)。总电流守恒 \(\nabla \cdot i_t = 0\)，跨膜电流 \(i_T = \chi(C_m \partial V/\partial t + I_{ion}) = \nabla \cdot (\sigma_i \nabla V_i)\)。

单域简化（12.106）：当 \(\sigma_i = \alpha \sigma_e\)（各向异性比相等）时，bidomain 化为 monodomain。实际心肌不满足此条件——12.4.2 实测各向异性比不等。

12.3.4 平面波速的各向异性

波速公式（12.114）：\(c = \Upsilon(n) c_0 / \sqrt{C_m R_m}\)，其中 \(\Upsilon(n)^2 = R_m (n \cdot \sigma_i n)(n \cdot \sigma_e n) / [\chi n \cdot (\sigma_i + \sigma_e) n]\) 是方向依赖空间常数。实测（Roberts & Scher 1982，狗心室肌）：\(\Delta V_{eL} = 74 \pm 7\) mV，\(\Delta V_{eT} = 43 \pm 6\) mV，\(\Delta V = 100\) mV → \(\sigma_{eL}/\sigma_{iL} = 0.35\)、\(\sigma_{eT}/\sigma_{iT} = 1.33\)。结合传播速度比 3:1 → \(\sigma_{iT}/\sigma_{iL} = 0.05\)、\(\sigma_{eT}/\sigma_{iL} = 0.07\)。

12.4 起搏器同步

SA 节点是异质自振细胞网络：空间频率梯度 + 随机偏离。关键问题：(1) 集体频率由谁决定？(2) firing 顺序由谁主导？(3) sinus node dysfunction 的发生条件？

反直觉结论：highest intrinsic frequency 的细胞不一定是 leader——gap junction 耦合产生"频率折中"（mutual entrainment），整个网络同步到中间频率。电紧张负荷（electrotonic loading）：fast 细胞被慢细胞拖累，slow 细胞被 fast 推动。

关键结论

  1. ECG 标量/向量理论：单导联 → Einthoven 三角 → 12 导联 → 心向量 \(H(t)\)。心向量 85% 信息集中在 3 导联，但临床仍用 12 导联以提高几何分辨率和冗余度。
  2. 细胞模型演化：Noble 1962（3 电流、\(C_m=12\) μF/cm²）→ MNT 1975（9 电流、含 Ca²⁺）→ Beeler-Reuter 1977（4 电流、j 再激活变量）→ Luo-Rudy 1994（15+ 电流的现代标准）。"模型越复杂 ≠ 越有用"——MNT 比 Noble 复杂但行为相似，是建模者的两难。
  3. SA 节点无 Na⁺ 电流：phase 0 由慢 Ca²⁺ 电流驱动，故可被 2 变量 FHN 方程定性描述。
  4. diadic cleft 的 0.2 离子问题：cleft 体积 \(2 \times 10^{-18}\) L，静息仅 0.2 个 Ca²⁺ 离子——经典确定性连续模型在 cleft 内失效，必须用随机模型。
  5. Common-Pool vs Local-Control 定理（Stern 1992）：common-pool 模型不能同时具有高增益和分级释放；高增益+分级释放需要 Local-Control（diadic cleft 单元独立 + 统计招募）。
  6. 传播是跳跃式（saltatory）的——Fast & Kleber 1993 实测：gap junction 占总电阻约 65%（远高于教科书的"低电阻"说法）。
  7. Bundle branch block 条件（12.86）：分叉处 \(A_1 p_1 / (A_2 p_2) \geq 1 - F(\beta)/F(\alpha)\)。兴奋性低的组织在分叉处更易阻滞（解释了缺血/低兴奋状态下束支阻滞高发）。
  8. Bidomain 模型的各向异性：心室肌的 \(\sigma_{eL}/\sigma_{iL} = 0.35\)、\(\sigma_{eT}/\sigma_{iT} = 1.33\)（不相等），导致非椭圆极化图样（dogbone）和虚拟电极。
  9. 平面波速方向依赖：\(c(n) \propto \sqrt{(n \cdot \sigma_i n)(n \cdot \sigma_e n) / [n \cdot (\sigma_i + \sigma_e) n]}\)，纵向 0.5 m/s vs 横向 0.17 m/s，比值 ~3。
  10. 虚拟电极（virtual electrodes）：强电场 + unequal anisotropy 在电场正交方向产生反极性区，是除颤的物理基础。

挑战和开放性问题

  1. RyR 关闭机制仍未达成共识——Stern 1999 等 4 状态模型只是候选之一，cleft Ca²⁺ 失活、SR 耗竭、随机关闭、适应 4 种机制都有支持证据，"there is not enough space here to discuss ... all the hypothetical models and the arguments for and against them"（12.2.7）。
  2. diadic cleft 的 0.2 离子问题——"traditional deterministic and continuous models may not even be applicable"。这意味着心室肌兴奋-收缩耦合的最核心物理过程（cleft 内 Ca²⁺ 动力学）尚无可信的第一性原理模型。
  3. ECG 逆问题（12.4）——从体表电位重建心内电流源，numerically unstable。临床上只能用单偶极子 + Burger 三角近似，无法定位小尺度梗死或精确传导路径。
  4. 电导张量实测的不一致性（12.4.2 末尾）——"one can find a variety of estimates of the anisotropy ratios in the literature"，原因是：(a) 测量方法（阻抗 vs 扩散张量成像）；(b) 心脏几何（心外膜 vs 心内膜 vs 中层）；(c) 疾病状态（缺血、心衰、肥厚）。
  5. Bidomain 模型的计算成本（12.3.3 末尾）——"is computationally expensive and there are not many situations where it matters"。心脏整体（10⁸ 细胞）的精细 bidomain 模拟在 2026 年仍需要 supercomputer，限制了临床转化。
  6. 激活序列的"非椭圆"传播（Fig. 12.27）——纤维旋转 + 曲面几何使最快传播方向不沿纤维方向（仅 60° 偏转），3D 整体心脏模拟的纤维方向数据（Hunter 1991, Nielsen 1991）仍不完整。
  7. Mayer 波在心脏中的对应——Ch 11 的 baroreflex 延迟微分方程解释了动脉压振荡，但心率本身的振荡（respiratory sinus arrhythmia、Cheyne-Stokes）的完整数学模型仍未统一。
  8. 胎儿-新生儿循环的分子细节——12 章未涉及，但新生儿窦房结功能不全（sinus node dysfunction）的离子机制与成人不同，可能涉及 \(I_f\)（"funny" current）和 HCN 通道的发育差异。

个人反思与批判性分析

本章的最大优点在于"模型演化的编年史"视角——把 Noble 1962 → MNT 1975 → BR 1977 → LR 1994 的连续 4 代模型并列展示，每个都给出参数表、方程、Fig 和局限。这种"边写边批判"的叙述方式使读者能直观感受离子通道模型的发展节奏——大约每 15-20 年一次大迭代，每次都因新电压钳数据而推翻前代。

本章的最大缺点在于"图多公式少"——12.1 ECG 节几乎全部是临床图（Fig. 12.3 各种心律失常、Fig. 12.5 Einthoven 三角、Fig. 12.6 6 导联），没有可解的方程。这反映了一个根本问题：ECG 的临床应用主要是模式识别而非方程求解——心脏科医生靠眼睛和经验，不靠 ODE。数学家写这一节时只能转录临床事实，缺乏分析性贡献。

12.2.5 节的 Common-Pool vs Local-Control 是本章最深的洞察。Stern 1992 的定理（"不能同时高增益+分级释放"）实际上是信息论层面的约束——单隔室系统的 SNR 不足以同时支持放大和保真。多隔室系统（每个 cleft 独立）通过统计多样性（diversity gain）突破这一约束。这与控制论中的"模块化系统优于集成系统"思想完全一致——Stern 的工作本质上是心电生理的"模块化定理"。

12.3.3 的 Bidomain 模型 + 12.4.2 的各向异性平面波速 一起构成了现代心脏电生理仿真的数学基础。然而本章对这一基础的处理偏教科书化——没有讨论：(a) bidomain 解的唯一性和正则性（PDE 理论）；(b) Stokes-Darcy 耦合的边界条件（心外膜与心包之间）；(c) 纤维旋转的微结构建模（Rule-based 方法）。这些缺失使得本章的 bidomain 部分对做实际仿真的人不够用，对做理论分析的人又太浅。

对我自己研究血管生物力学的启发：(1) 12.3.2 的传播失败条件 (12.86) 本质是几何不连续性导致的失稳阈值——这与血管分叉处的血流分离、动脉瘤入口的壁剪切应力集中、动脉硬化斑块肩部的应力集中在数学上同构——都是"局部不均匀性破坏全局稳态"的标准问题；(2) 12.4.2 的各向异性电导张量与血管壁的各向异性力学（周向 vs 轴向弹性模量差异 ~10 倍）在 PDE 结构上完全相同——一个是对 \(\sigma \nabla V\) 求解，一个是对 \(\sigma \nabla u\) 求解；(3) 12.2.5 的 diadic cleft 0.2 离子问题是任何亚细胞钙动力学的根本性挑战——血管平滑肌细胞中也有类似的"微域"（microdomain），cytosol 与膜附近的 Ca²⁺ 浓度可差 10-100 倍，血管钙化的研究必须面对同样的统计性问题。

值得指出，本章的 12.5 节（除颤虚拟电极机制）虽然占篇幅不大，但实际上是临床除颤器设计的理论基础——双向波（biphasic waveform）比单向波更有效，正是因为反向电流可以"擦除"狗骨头图样中的虚拟电极。这一节是 Vol I+II 二十章中临床转化最成功的数学例子。

重要参考文献

[X1] Noble, D. (1962). A modification of the Hodgkin–Huxley equations applicable to Purkinje fibre action and pacemaker potentials. Journal of Physiology, 160, 317-352. — 12.2.1 节 Purkinje 纤维的 HH 式模型原型；3 离子电流、\(C_m=12\) μF/cm²。 [X2] McAllister, R. E., Noble, D., & Tsien, R. W. (1975). Reconstruction of the electrical activity of cardiac Purkinje fibres. Journal of Physiology, 251, 1-59. — 12.2.2 节 MNT 9 电流模型；首次包含 Ca²⁺ 慢内向电流 \(I_{si}\)。 [X3] Yanagihara, K., Noma, A., & Irisawa, H. (1980). Reconstruction of sino-atrial node pacemaker potential based on the voltage clamp experiments. Japanese Journal of Physiology, 30, 841-857. — 12.2.3 节 SA 节点 YNI 模型；首次明确 SA 节点无快 Na⁺ 电流。 [X4] Beeler, G. W., & Reuter, H. (1977). Reconstruction of the action potential of ventricular myocardial fibres. Journal of Physiology, 268, 177-210. — 12.2.4 节心室肌 BR 模型；4 电流 + j 再激活变量。 [X5] Luo, C. H., & Rudy, Y. (1994a,b). A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes; II. Afterdepolarizations, triggered activity, and potentiation. Circulation Research, 74, 1071-1096, 1097-1113. — 12.2.4 末尾提及的 LR-I 和 LR-II 模型；现代心室电生理标准模型。 [X6] Stern, M. D. (1992). Theory of excitation-contraction coupling in cardiac muscle. Biophysical Journal, 63, 497-517. — 12.2.5 节 Local-Control 理论的奠基；common-pool 模型不能同时高增益+分级释放的定理。 [X7] Jafri, M. S., Rice, J. J., & Winslow, R. L. (1998). Cardiac Ca²⁺ dynamics: The roles of ryanodine receptor adaptation and sarcoplasmic reticulum load. Biophysical Journal, 74, 1149-1168. — 12.2.6 节 L 型 Ca²⁺ 通道的 Monod-Wyman-Changeux Markov 模型。 [X8] Bers, D. M. (2001). Excitation-Contraction Coupling and Cardiac Contractile Force (2nd ed.). Kluwer. — 12.2.4 节 EC 耦合实验文献的最全面综述；Fig. 12.16 心肌 Ca²⁺ 通量图出处。 [X9] Tung, L. (1978). A bidomain model for describing ischemic myocardial D-C potentials. PhD thesis, MIT. — 12.3.3 节 Bidomain 模型的首次提出。 [X10] Henriquez, C. S. (1993). Simulating the electrical behavior of cardiac tissue using the bidomain model. Critical Reviews in Biomedical Engineering, 21, 1-77. — Bidomain 模型的综述；现已被广泛接受为心肌电活动的标准模型。 [X11] Roth, B. J. (1992). How the anisotropy of the intracellular and extracellular conductivities influences stimulation of cardiac muscle. Journal of Mathematical Biology, 30, 633-646. — 12.4 节虚拟电极（virtual electrodes）的首次提出。 [X12] Wikswo, J. P., Lin, S. F., & Abbas, R. A. (1995). Virtual electrodes in cardiac tissue: A common mechanism for anodal and cathodal stimulation. Biophysical Journal, 69, 2195-2210. — 虚拟电极的统一理论；阴极和阳极刺激都通过虚拟电极起作用。 [X13] Stern, M. D., Song, L. S., Cheng, H., et al. (1999). Local control models of cardiac excitation-contraction coupling: A possible role for allosteric interactions between ryanodine receptors. Journal of General Physiology, 113, 469-489. — 12.2.7 节 RyR 的 4 状态 Markov 模型。 [X14] Fast, V. G., & Kleber, A. G. (1993). Microscopic conduction in cultured strands of neonatal rat heart cells measured with voltage-sensitive dyes. Circulation Research, 73, 914-925. — 12.3.1 节 saltatory 传播的实测证据；gap junction 延迟 80 μs。 [X15] Roberts, D. E., & Scher, A. M. (1982). Effect of tissue anisotropy on extracellular potential fields in canine myocardium in situ. Circulation Research, 50, 342-351. — 12.4.2 节心室肌电导张量的各向异性实测；\(\sigma_{eL}/\sigma_{iL} = 0.35\)、\(\sigma_{eT}/\sigma_{iT} = 1.33\)。 [X16] Rushmer, R. F. (1976). Cardiovascular Dynamics (4th ed.). W. B. Saunders. — 12.1 ECG 节各临床图（Fig. 12.2-12.10）的主要出处。 [X17] DiFrancesco, D., & Noble, D. (1985). A model of cardiac electrical activity incorporating ionic pumps and concentration changes. Philosophical Transactions of the Royal Society B, 307, 353-398. — 12.2.4 末尾提及的 MNT 后继 Purkinje 模型；首次将离子泵动力学整合入细胞模型。