跳转至

第5章：兴奋性

作者

本章作者与 Ch1-Ch4 相同——James P. Keener 与 James Sneyd。本章是 I 卷"细胞生理学"中承前启后的核心章——上承 Ch2/Ch3 的电生理基础与单通道 Markov，下启 Ch6 行波传播、Ch9 神经内分泌、Ch10 基因调控、Ch12 心脏电生理、Ch19 视网膜。本章的标志性论断（§5.1 开头）："Hodgkin-Huxley 方程是整个生理学文献中最重要的模型"。本章的角色是：（1）完整推导 HH 模型（§5.1）——从电压钳数据到 4 维 ODE；（2）相平面分析（§5.1.3）——快相平面 + 快-慢相平面 + 阈值与振荡的几何机制；（3）FitzHugh-Nagumo 简化（§5.2）——把 4 维 ODE 化简为 2 维，揭示兴奋性的"通用"结构；（4）Morris-Lecar、Pushchino、Belousov-Zhabotinsky 等变体——展示兴奋性是跨学科的普适现象。

内容概述

本章围绕"什么是兴奋性、如何分析兴奋系统、兴奋性如何与振荡耦合"这三个核心问题展开。

第一部分（§5.1）HH 模型的完整推导。§5.1.1 简述 Hodgkin-Huxley 1952 的历史——Cole-Hodgkin 1939 电压钳 + 1952 数据拟合 + Huxley 的"3 周手摇计算器"完成动作电位的数值预测。§5.1.2 详细描述门控变量 \(n, m, h\) 的识别过程： - K⁺ 门控（5.11）\(g_K = \bar g_K n^4\)（4 个独立门控粒子，每粒子有 open/closed 两态，故开放概率 \(\propto n^4\)） - Na⁺ 门控（5.18）\(g_{Na} = \bar g_{Na} m^3 h\)（3 个激活粒子 + 1 个失活粒子） - 门控变量动力学（5.13）\(dw/dt = \alpha_w(1-w) - \beta_w w\)（first-order Markov 状态图） - 经验拟合函数 \(\alpha_w(v), \beta_w(v)\)（5.24）–(5.29)——形式选择是经验拟合，无生理学基础

HH 完整方程（5.20）–(5.23)：

\[ C_m \frac{dv}{dt} = -\bar g_K n^4 (v - v_K) - \bar g_{Na} m^3 h (v - v_{Na}) - \bar g_L (v - v_L) + I_{app} \]
\[ \frac{dm}{dt} = \alpha_m(1-m) - \beta_m m, \quad \frac{dn}{dt} = \alpha_n(1-n) - \beta_n n, \quad \frac{dh}{dt} = \alpha_h(1-h) - \beta_h h \]

动作电位的四阶段（§5.1.2 末段）：upstroke (m 快速上升) → excited (\(g_{Na}\) 饱和) → refractory (\(h\) 失活, \(n\) 激活) → recovery (\(n, h\) 返回静息)。关键是时间尺度分离：\(\tau_m \ll \tau_n, \tau_h\)——\(m\) 是"快变量"，\(n, h\) 是"慢变量"。

振荡分析（§5.1.2 末段）：\(I_{app}\) 作为分岔参数。在 \(I_{app} = 9.78\) 发生 subcritical Hopf 分岔（静息失稳），中间范围有稳定极限环振荡，再次 Hopf 回到稳态。生物医学意义：高胞外 K⁺ 通过提升 \(V_K\) 让静息超阈——这是某些心律失常的机制。

第二部分（§5.1.3）相平面分析。快相平面（固定 \(n = n_0, h = h_0\)）：(5.31)-(5.32) 给出 2 维 \((v, m)\) 系统。三个伪稳态 \(v_r, v_s, v_e\)（\(v_r, v_e\) 稳定，\(v_s\) 鞍点）——\(v_s\) 的稳定流形（1 维）构成阈值的几何表现。任何越过稳定流形的扰动会触发动作电位；未越过的扰动会衰减。快-慢相平面（用 \(m = m_\infty(v)\) 准稳态 + \(h + n \approx 0.8\) 经验关系）：(5.35)-(5.36) 化为 2 维 \((v, n)\) 系统。核心几何：\(v\) 的零线 \(f(v, n) = 0\) 是立方曲线（cubic），三个分支——左/右稳定、中不稳定。阈值 = 中分支作为不可穿越的"虚墙"。动作电位的几何图像：\(v\) 跳到上分支 → 沿上分支缓慢上升（\(n\) 增加）→ 达到折点 → 跳到下分支 → 沿下分支缓慢下降 → 回到静息。

第三部分（§5.2）FitzHugh-Nagumo 方程。FitzHugh 1960, 1961, 1969 把 HH 快-慢相平面的本质提炼为 2 维通用模型（5.38）–(5.40)：

\[ \epsilon \frac{dv}{dt} = f(v) - w + I_{app}, \quad \frac{dw}{dt} = v - \gamma w \]
\[ f(v) = v(1-v)(v-\alpha), \quad 0 < \alpha < 1, \quad \epsilon \ll 1 \]

关键解释（§5.2 开头）：\(v\) 是激发变量（fast），\(w\) 是恢复变量（slow）。这是从 HH 4 维 ODE 到 2 维 ODE 的"模型降维"——保留了兴奋性的核心几何结构。§5.2 末段给出其他变体： - McKean 1970 开关 + 线性（5.41）\(f(v) = H(v-\alpha) - v\)（Heaviside 开关） - Pushchino 模型（5.42）——Krinsky, Panfilov, Pertsov, Zykov 开发，心室动作电位的简化 - van der Pol 振子（5.48）\(v'' + a(v^2-1)v' + v = 0\)—— van der Pol & van der Mark 1928 提出的心脏起搏器模型

广义 FitzHugh-Nagumo 理论（§5.2.1–§5.2.2）：单稳态在左/右分支时系统兴奋（Fig. 5.15）；单稳态在中分支时系统振荡（Fig. 5.16）——两种行为是同一 ODE 的不同参数区域。奇异性摄动分析（§5.2.2 末段）：\(\epsilon \ll 1\) → \(v\) 是快变量，\(w\) 是慢变量；\(v\) 沿 \(f(v,w) = 0\) 的稳定分支滑动，直到 \(w\) 达到折点 \(W^*\) → 跳到另一分支 → 周期 \(T\) 由 (5.59) 给出 \(T = \int (1/G^+(w) - 1/G^-(w)) dw\)。

分岔图（§5.2.2 末段 + Fig. 5.17）：\(I_{app} = 0.1\) 第一次 Hopf 出现稳定振荡，\(I_{app} = 1.24\) 第二次 Hopf 消失——两个 Hopf 之间是振荡区域。

第四部分（§5.2.3–§5.3）模型变体与电路实现。§5.2.3 讨论 van der Pol、McKean、Pushchino 等变体的几何与生理意义。Morris-Lecar 模型（Exercise 12）——藤壶肌纤维 2 维简化（\(V\) + 单一门控 \(W\)）。Belousov-Zhabotinsky 反应（Exercise 14）——非生理兴奋系统（柠檬酸溴酸氧化）的化学振荡。§5.3 Exercises 15 的电路实现——Keener 1983 用 LM 741 op-amp 实现了真正的 FHN 模拟电路——这是 §5.2 方程不是"纯数学抽象"而是"可以物理实现"的明证。

主要数学结果： 1. HH 完整方程 (5.20)–(5.23) + 拟合函数 (5.24)–(5.29) 2. 快相平面三伪稳态 + \(v_s\) 稳定流形 = 阈值 3. 快-慢相平面 cubic 零线 + 中分支作为不可穿越墙 4. 动作电位四阶段（upstroke, excited, refractory, recovery）对应 \((m, n, h)\) 的时间尺度分离 5. FitzHugh-Nagumo 通用模型 (5.38)–(5.40)：\(f(v) = v(1-v)(v-\alpha)\) + \(\epsilon \ll 1\) 6. Hopf 分岔：\(I_{app} = 0.1\)（振荡出现）→ \(I_{app} = 1.24\)（消失） 7. 奇异性摄动周期 (5.59)：\(T = \int (1/G^+ - 1/G^-) dw\)

前置知识：ODE 稳定性（线性化、Hopf 分岔）、相平面分析（nullcline、稳定流形）、奇异性摄动（Ch1.6 已奠定基础）、\(g_K = \bar g_K n^4\) 的"独立门控粒子"假设（Ch3 Markov 状态图的具体应用）。

核心方程与概念

1. HH 完整方程

4 维 ODE（5.20）–(5.23)：

\[ C_m \frac{dv}{dt} = -\bar g_K n^4 (v - v_K) - \bar g_{Na} m^3 h (v - v_{Na}) - \bar g_L (v - v_L) + I_{app} \]

经验拟合（5.24）–(5.29)：

\[ \alpha_m = \frac{0.1(25 - v)}{\exp((25-v)/10) - 1}, \quad \beta_m = 4 e^{-v/18} \]

（\(\alpha_n, \beta_n, \alpha_h, \beta_h\) 形式类似，参数不同。）

参数（5.30）：\(\bar g_{Na} = 120\) mS/cm², \(\bar g_K = 36\) mS/cm², \(\bar g_L = 0.3\) mS/cm², \(C_m = 1\) μF/cm², \(v_{Na} = 115\) mV（shifted）, \(v_K = -12\) mV, \(v_L = 10.6\) mV。温度 6.3°C。

注意（§5.1.2 末段）：\(g_{Na} = \bar g_{Na} m^3 h\) 中 \(m^3\) 和 \(h^1\) 的指数是经验拟合，不是从通道结构推导。Ch3.4.2 给出 Markov 状态图的微观解释：3 个独立激活粒子 + 1 个失活粒子，每粒子有 2 态 → 开放概率 \(m^3 h\)。

2. 时间尺度分离与动作电位的四阶段

关键观察（§5.1.2 末段 + Fig. 5.5）：\(\tau_m(v) \ll \tau_n(v), \tau_h(v)\)——\(m\) 是快变量，\(n, h\) 是慢变量。

动作电位的四阶段： - Upstroke（~1 ms）：\(v\) 上升 → \(m\) 快速上升（\(\tau_m\) 短）→ \(g_{Na}\) 上升 → \(v\) 进一步上升（正反馈） - Excited：\(m \to 1\), \(g_{Na}\) 饱和 → \(v \to v_{Na}\)，但 \(h\) 仍在下降（\(\tau_h\) 长）→ \(g_{Na}\) 减少 - Refractory：\(h\) 下降到接近 0 → \(g_{Na} \to 0\)；同时 \(n\) 上升（\(\tau_n\) 长）→ \(g_K\) 上升 → \(v\) 下降 - Recovery：\(h, n\) 缓慢返回静息值 → \(v\) 回到静息

为什么有阈值（§5.1.3）：\(v\) 需要越过 \(m\) 零线上的临界点 \(v_{th}\)，使 \(m_\infty(v_{th}) > m\) 静息值——一旦越过，正反馈让 \(v\) 持续上升。

3. 相平面分析：快相平面

2 维简化（5.31）–(5.32)：固定 \(n = n_0, h = h_0\)，得 \(v\) 和 \(m\) 的 ODE。

三个伪稳态（§5.1.3 中段）：\(v_r\)（静息，\(v\) 接近 \(v_K\)）, \(v_s\)（鞍点）, \(v_e\)（激发，\(v\) 接近 \(v_{Na}\)）。\(v_s\) 的稳定流形是 1 维曲线——把 \((m, v)\) 平面分成"会回到 \(v_r\)"和"会跑到 \(v_e\)"两个区域。

阈值 = 稳定流形。这是几何阈值——不是某个抽象电压值。

4. 相平面分析：快-慢相平面

2 维简化（5.35）–(5.36)：\(m = m_\infty(v)\)（准稳态）+ \(h + n \approx 0.8\)（经验关系）。

核心几何（§5.1.3 + Fig. 5.11）：\(v\) 零线 \(f(v, n) = 0\) 是 cubic： - 左分支（\(v < V^-\)）：稳定 - 中分支（\(V^- < v < V^+\)）：不稳定 - 右分支（\(v > V^+\)）：稳定

动作电位几何路径：\(v\) 从静息（\(v = v_r\)，左分支）→ \(v\) 越过中分支（"虚墙"）→ 跳到右分支（\(v\) 上升到 \(V^+\)）→ 沿右分支缓慢上升（\(n\) 增加）→ 到达折点 → 跳到左分支 → 沿左分支缓慢下降（\(n\) 减少）→ 回到静息。

振荡（§5.1.3 末段 + Fig. 5.12）：当 \(I_{app}\) 足够大时，cubic 零线整体上移，使得 \(n_\infty(v) = n\) 零线与 cubic 中分支相交——稳态变成不稳定，轨迹永远在两稳定分支间跳转——弛豫振荡（relaxation oscillation）。

5. FitzHugh-Nagumo 通用模型

核心方程（5.38）–(5.40）：

\[ \epsilon \frac{dv}{dt} = v(1-v)(v-\alpha) - w + I_{app}, \quad \frac{dw}{dt} = v - \gamma w \]

注意（§5.2 开头）：FHN 不是 HH 的精确简化，而是 HH 快-慢相平面的"提炼"——保留了兴奋性的核心几何（cubic 零线 + 单调 \(w\) 零线），但放弃了 HH 的所有具体参数。

两类行为（Fig. 5.15 vs 5.16）： - Excitable（\(I_{app} = 0\)）：单稳态在左/右分支，稳定。小扰动衰减；超阈扰动触发一次大振荡 - Oscillatory（\(I_{app} = 0.5\)）：单稳态在中分支，不稳定。弛豫振荡——\(v\) 周期性在两稳定分支间跳转

奇异性摄动分析（§5.2.2 末段）：\(\epsilon \ll 1\) → \(v\) 沿 cubic 零线的稳定分支滑动。周期（5.59）：

\[ T = \int_{W_*}^{W^*} \left(\frac{1}{G^+(w)} - \frac{1}{G^-(w)}\right) dw \]

其中 \(G^\pm(w) = g(V^\pm(w), w)\)——这是沿"兴奋分支"和"恢复分支"的时间积分之差。

6. 模型变体：Morris-Lecar / Pushchino / BZ

Morris-Lecar 1981（Exercise 12）：藤壶肌纤维 2 维简化，\(V\) + 单一 K⁺ 门控 \(W\)（无 Na⁺ 失活变量 \(h\)）。优点：只有 1 个时间尺度（K⁺ 激活）→ 几何更清晰。

Pushchino 模型（Exercise 13, Panfilov 等）：心室动作电位的简化。分段线性 \(F(v)\)（5.80）+ 不连续 \(\tau(v)\)（5.81）——这是 §5.2 的 FHN + 心脏电生理（Ch12）的桥梁。

Belousov-Zhabotinsky 反应（Exercise 14, Tyson-Fife 1980）：化学振荡的经典模型。变量 \(u\) = 溴酸浓度，\(v\) = 氧化催化剂金属浓度。这是"兴奋性"超越生理学的明证——BZ 反应是纯化学反应，但具有与神经动作电位相同的几何结构。

7. 物理实现：op-amp 电路

Keener 1983（Exercise 15）用 LM 741 op-amp 实现了 FHN 模拟电路。核心：op-amp 的分段线性输出（rail voltage ±）正好对应 FHN 的分段线性 \(f(v)\)。\(f(v)\) 形式（5.89）： - \(v > \alpha V_{r+}\)：\(F(v) = v - V_{r+}\)（饱和上分支） - \(\alpha V_{r-} \le v \le \alpha V_{r+}\)：\(F(v) = -(R_1/R_2) v\)（中间线性） - \(v < \alpha V_{r-}\)：\(F(v) = v - V_{r-}\)（饱和下分支）

意义：FHN 不仅是数学抽象，它是可物理实现的。这为神经形态计算（neuromorphic computing）提供了理论基础。

关键结论

  1. HH 方程是"最成功的生物物理模型"（§5.1 开头）：从 1952 的 4 篇论文预测动作电位的形状、传导速度、阈值、不应期——之后 70 年的所有电生理学教科书都建立在此基础上。

  2. 兴奋性的几何本质是 cubic 零线（§5.1.3 + §5.2）：动作电位的阈值 = cubic 中分支作为"虚墙"；振荡 = 单稳态在中分支 → 弛豫振荡。这是兴奋性的普适几何——Bénard 对流、BZ 反应、动作电位都共享此结构。

  3. 时间尺度分离是兴奋性的核心（§5.1.2 末段 + §5.2.2）：\(m\) 的快动力学 + \(n, h\) 的慢动力学让 cubic 几何自然涌现。HH 的 4 维 ODE 可被约化为 2 维 FHN——这是奇异性摄动在生物物理的经典胜利。

  4. FHN 是"通用模型"而非"具体模型"（§5.2 开头）：FHN 不描述任何具体生物现象，但抓住了兴奋性的所有核心数学结构。这是"理论模型 vs 现象模型"分工的典范。

  5. 振荡 = 兴奋性的"超阈自持续"（§5.1.2 末段 + §5.2.2）：当 \(I_{app}\) 足够大时，单次动作电位后系统不能回到静息（因为 \(h\) 失活未恢复 + \(n\) 激活未衰退）→ 自持续振荡。这是心律失常、自持续神经元、簇发放的共同机制。

  6. 温度对兴奋性的影响是非平凡的（Exercise 10）：HH 方程从 6.3°C 到 30°C，兴奋性可能变高或变低——取决于哪些参数是温度敏感的。\(Q_{10}\) 效应让单一模型在不同温度下表现不同——这给实验拟合带来额外复杂度。

  7. m, n, h 的指数（3, 4, 1）没有生理学基础（§5.1.2 中段）：HH 选择这些指数是经验拟合而非通道结构推导。Ch3.4.2 给出 Markov 状态图的微观对应：3 个独立激活粒子 + 1 个失活粒子 = \(m^3 h\)。K⁺ 通道的 \(n^4\) 对应 4 个独立门控粒子——真实 K⁺ 通道 4 亚基结构（MacKinnon 1998）的预测。

  8. FHN 可物理实现（Exercise 15）：op-amp 电路的 rail voltage 饱和 + 中间线性 = 分段线性 \(F(v)\)。这是神经形态计算（neuromorphic computing）的理论基础——Mead 1989 开创的"硅神经元"就是 FHN 电路的工程实现。

挑战和开放性问题

  1. m, n, h 的"经验性"（§5.1.2 中段）：HH 选择 \(m^3 h\) 和 \(n^4\) 是"最小可接受指数"——但真实钠通道 4 个同源结构域（DI-DIV）独立控制激活和失活，\(m^3 h\) 是否太简化？现代钠通道模型（Payandeh 2011 晶体结构）显示激活/失活机制远比 HH 复杂。

  2. 温度依赖（Exercise 10）：HH 原始方程 6.3°C，温度敏感性如何标度？Nernst 方程的 \(RT/F\) 因子、速率常数 \(\alpha, \beta\) 的 Q10 效应、电容的温度依赖——所有这些让 HH 在不同温度下是不同模型。

  3. 细胞类型差异（§5.1.2 末段 + §5.1.1 末段）：HH 适用于鱿鱼巨大轴突；脊椎动物神经元有不同 I-V 曲线（Frankenhaeuser 1960）、不同离子通道（特别是 Ca²⁺ 通道）；心肌细胞有完全不同的门控（Ch12）。HH 的"普适性"有限。

  4. 簇发放 (bursting)（Ch9 主题）：HH + 慢变量（如 Ca²⁺ 激活的 K⁺ 通道）可产生多个动作电位的簇——但 Ch5 没有详细讨论。这是 Ch5 → Ch9 的自然过渡。

  5. 2 维 vs 4 维（§5.1.3 + §5.2）：FHN 2 维保留兴奋性的核心，但丢失了 HH 的具体离子电流结构。当需要描述具体药物作用（如 TTX 阻断 Na⁺、TEA 阻断 K⁺）时，必须回到 HH 4 维。这是"模型保真度 vs 可分析性"的根本张力。

  6. \(h + n \approx 0.8\) 经验关系（§5.1.3 中段）：FHN 简化依赖此关系。真实动作电位中 \(h + n\) 不是常数——Ch9 簇发放中 \(h, n\) 完全解耦。FHN 简化对单个动作电位有效，对长期动力学可能失效。

  7. 多离子种类的兴奋性（Ch5 隐含 + Ch12）：HH 只考虑 Na⁺ + K⁺；心肌细胞还有 Ca²⁺ 通道（特别是 L 型 Ca²⁺ 通道的"钙诱导钙释放"）——这让心肌动作电位的"平台相"持续 200-300 ms，远长于 HH 的 1-2 ms。Ch5 模型对心肌严重不足。

  8. 随机性（§5.1 + Ch3 衔接）：HH 是确定性 ODE。真实单通道（Ch3）是离散的 Markov 跳变。对小型神经元（如 C. elegans 的某些神经元），单通道随机性显著影响动作电位的发放时机——这需要 Langevin 扩展或 Monte Carlo 单通道模拟（Ch2.9 Gillespie）。

  9. FHN 的"通用性"边界（§5.2.1 末段）：FHN 是"cubic 零线 + 单调恢复零线"的两类系统。真实生物系统有无 cubic 形状但多变量（如胰岛 β 细胞的簇发放，Ch9）——这些"高维兴奋系统"需要新的分析工具。

  10. 神经形态计算的工程化（Exercise 15 隐含）：Keener 1983 给出 FHN 电路，Mead 1989 实现硅神经元。现代神经形态芯片（如 Intel Loihi、IBM TrueNorth）规模达百万神经元——但单个神经元的简化模型与 FHN 的关系如何？是否需要"更逼真"的模型？

个人反思与批判性分析

建模哲学层面：

Ch5 是 I 卷"建模哲学的巅峰"——HH 1952 是生物数学建模的金标准：（1）从数据出发（Hodgkin 的电压钳数据），（2）建立最小可拟合模型（4 维 ODE），（3）预测新现象（动作电位传导速度），（4）70 年持续验证。没有其他生物模型有这样的成功史。

值得专门提出的是：FitzHugh-Nagumo 简化（§5.2）是"降维建模"的典范——HH 4 维 ODE 化简为 FHN 2 维 ODE，保留了兴奋性的核心结构（cubic 零线 + 时间尺度分离）。这种"几何建模"是动力学系统理论的胜利——不需要求解具体方程，只看相平面结构就能预测行为。

与后续章节的衔接：

  • Ch6（行波传播）：HH + Ch4 电缆方程 = 神经冲动沿轴突传播。FHN + 电缆方程给出行波的解析条件——奇异性摄动 + 速度公式。
  • Ch7（钙动力学）：Ca²⁺ 通道是 Ch3 Markov 模型的多离子例子；Ch5 的兴奋性 + 慢 Ca²⁺ 动力学 = 簇发放 (bursting)。
  • Ch8（细胞间通讯）：突触释放是门控变量的快速开关——Ch5 的 \(m, n, h\) 时间尺度在 Ch8 有直接应用。
  • Ch9（神经内分泌）：簇发放（多个动作电位的串）= HH 兴奋性 + 慢 Ca²⁺ 反馈。这是 Ch5 + Ch7 的综合。
  • Ch10（细胞功能调控）：基因调控网络的双稳态 = 慢变量的"立方"几何——FHN 的"通用性"扩展到生化网络。
  • Ch12（心脏电生理）：心肌动作电位的"平台相"需要 Ca²⁺ 通道——HH 不够，需要扩展模型。

值得重新推导的方程：

  • (5.20)–(5.23) HH 完整方程：从 4 个门控粒子假设（\(m^3 h\) 激活 + \(h\) 失活；\(n^4\) 激活）出发，写下开放概率 → 写出 \(g_{Na}, g_K\) → 写出 4 维 ODE。
  • (5.31)–(5.32) 快相平面：固定 \(n = n_0, h = h_0\)，\(v\) 和 \(m\) 各自 ODE。关键：(5.34) 给出 \(v\) 零线是关于 \(m\) 的三次曲线——三个伪稳态 + 一个鞍点 + 阈值几何。
  • (5.35)–(5.36) 快-慢相平面：\(m = m_\infty(v)\) 准稳态 + \(h + n \approx 0.8\) → 2 维 \((v, n)\) 系统。关键：(5.37) 立方零线的几何意义。
  • (5.59) 奇异性摄动周期：从 \(G^\pm(w)\) 沿稳定分支的慢动力学 + \(W_*\) 折点处快跳变 → 积分。关键：周期 \(T\) 不依赖 \(\epsilon\)（只要 \(\epsilon\) 足够小）。
  • (5.84) op-amp 动力学：从输入电压差 \(v_+ - v_-\) 到输出 \(v_0\)。关键：\(\epsilon_s \to 0\) 时 \(v_0\) 跳到 rail voltage——分段线性 FHN 的物理实现。

对作者建模选择的评价：

  • 优点：（1）HH 完整推导的实验+理论结合（§5.1.1 历史 + §5.1.2 数据拟合）极为清晰；（2）相平面分析（§5.1.3）把 4 维 ODE 化简为 2 维 + 几何直觉——这是动力系统理论的标准应用；（3）FHN（§5.2）的"通用模型"提炼抓住了兴奋性的本质；（4）§5.2.2 奇异性摄动 + 周期公式 (5.59) 是经典的工程数学；（5）§5.3 Exercise 15 的物理电路实现是"理论联系实际"的范例。
  • 不足：（1）HH 4 维 ODE 的完整 Hopf 分岔分析缺失——§5.1.2 末段只给出 \(I_{app} = 9.78\) 的数值，没有给出解析判据；（2）温度依赖只在 Exercise 10 提到——HH 在不同温度下的参数调整是应用的关键；（3）细胞类型差异（心肌、骨骼肌、分泌细胞）只在 §5.1 末段一句话提及；（4）随机性（Ch3 Markov + Ch5 ODE 的衔接）没有明确讨论——单通道随机对小型神经元动作电位的影响。

若能与作者对话，我会问：

  1. HH 的现代扩展：Hodgkin-Huxley 模型的核心是"m³h 激活 + n⁴ 激活"——但现代钠通道（Payandeh 2011 晶体结构）显示 4 个同源结构域，激活机制远比 \(m^3 h\) 复杂。HH 是否应该被"参数化"为 4 维独立门控？
  2. FHN 的边界：FHN 2 维模型抓住了兴奋性核心，但对簇发放（Ch9）失效——需要 3 维或更高维。是否有一个"统一的 3 维 FHN"？
  3. 神经形态计算：现代神经形态芯片（Intel Loihi 2、IBM TrueNorth）规模达百万神经元，但每个神经元用简化的 LIF 模型（LIF = leaky integrate-and-fire，无 cubic 零线）。LIF vs FHN：为何 LIF 在工程上更流行？是否因为 LIF 更易数字实现？
  4. 多离子扩展：心肌动作电位的"平台相"需要 Ca²⁺ 通道（特别是 L 型）。HH + Ca²⁺ 通道的标准扩展是什么？Hodgkin-Huxley + Luo-Rudy 1991？还是其他？
  5. 随机性的工程化：真实神经元的离子通道数有限（小型神经元可能只有 100-1000 个 Na⁺ 通道）→ 单通道随机性显著。HH 的随机版本（stochastic HH）如何实现？是否能让单通道 Markov（Ch3）与集总 ODE（Ch5）形成统一的随机-确定性谱？

重要参考文献

[X1] Hodgkin AL, Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 1952;117:500–544. （HH 模型原始论文——本章 §5.1 的全部基础。）

[X2] Hodgkin AL, Huxley AF. The components of membrane conductance in the giant axon of Loligo. J Physiol 1952a;116:473–496.

[X3] Hodgkin AL, Huxley AF. Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo. J Physiol 1952b;116:449–472.

[X4] Hodgkin AL, Huxley AF. The dual effect of membrane potential on sodium conductance in the giant axon of Loligo. J Physiol 1952c;116:497–506.

[X5] Hodgkin AL, Huxley AF. A quantitative description of the membrane current and its application to conduction and excitation in nerve. J Physiol 1952d;117:500–544. （"5 篇文章"中的最后一篇——拟合函数和动作电位预测。）

[X6] Hodgkin AL, Katz B. The effect of sodium ions on the electrical activity of the giant axon of the squid. J Physiol 1949;108:37–77. （Ch2 已引——HH 早期识别 Na⁺ 重要性。）

[X7] Cole KS, Curtis HJ. Electric impedance of the squid giant axon during activity. J Gen Physiol 1939;22:649–670. （膜电导增加的首次实验观察。）

[X8] FitzHugh R. Impulses and physiological states in theoretical models of nerve membrane. Biophys J 1961;1:445–466. （FitzHugh 简化——§5.2 的源头。）

[X9] Nagumo J, Arimoto S, Yoshizawa S. An active pulse transmission line simulating nerve axon. Proc IRE 1964;50:2061–2070. （Nagumo 电路实现——Exercise 15 的前身。）

[X10] van der Pol B, van der Mark J. The heartbeat considered as a relaxation oscillation, and an electrical model of the heart. Philos Mag 1928;6:763–775. （van der Pol 振子的心脏起搏器模型。）

[X11] Rinzel J. Discussion: One-dimensional and two-dimensional models of nerve activity. In: Mathematical Approaches to Cardiac Arrhythmias, ed. by Zipes DP, Jalife J, Saunders, 1990. （§5.1.1 历史的来源。）

[X12] Morris C, Lecar H. Voltage oscillations in the barnacle giant muscle fiber. Biophys J 1981;35:193–213. （Morris-Lecar 模型。）

[X13] Rinzel J. On the repetitive activity of excitable membranes. Bull Math Biol 1978;40:503–514. （HH 振荡分析。）

[X14] Troy WC. The bifurcation of periodic solutions in the Hodgkin-Huxley equations. Q Appl Math 1978;36:73–83.

[X15] McKean HP. Nagumo's equation. Adv Math 1970;4:209–223. （分段线性 FHN——(5.41)。）

[X16] Tyson JJ, Fife PC. Target patterns in a realistic model of the Belousov-Zhabotinskii reaction. J Chem Phys 1980;73:2224–2237. （BZ 反应——Exercise 14。）

[X17] Krinsky VI, Panfilov AV, Pertsov AM. A mathematical model of the mechanism of the heart muscle contraction. Biofizika 1972;17:839–845. （Pushchino 模型——心室动作电位简化。）

[X18] Panfilov AV, Holden AV (eds.). Computational Biology of the Heart. Wiley, 1997.

[X19] Strogatz SH. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering. Addison-Wesley, 1994. （Hopf 分岔理论的入门标准教材。）

[X20] Keener JP. Analog circuitry for the van der Pol and FitzHugh-Nagumo equations. IEEE Trans Syst Man Cybern 1983;13:1010–1014. （§5.3 Exercise 15 的电路实现。）

[X21] Mead C. Analog VLSI and Neural Systems. Addison-Wesley, 1989. （神经形态计算的开山之作。）

[X22] Armstrong CM, Bezanilla F. Inactivation of the sodium channel. II. Gating current experiments. J Gen Physiol 1977;70:567–590. （Na⁺ 失活门控电流的实验。）

[X23] Payandeh J, Scheuer T, Zheng N, Catterall WA. The crystal structure of a voltage-gated sodium channel. Nature 2011;475:353–358. （钠通道晶体结构——HH 模型现代化的基础。）

[X24] MacKinnon R. Determination of the subunit stoichiometry of a voltage-activated potassium channel. Nature 1991;350:232–235. （K⁺ 通道 4 亚基结构——HH \(n^4\) 的物理基础。）

[X25] Cole KS. Membranes, Ions and Impulses. University of California Press, 1968. （Cole 早期工作的回顾——§5.1.1 历史的另一来源。）