FitzHugh-Nagumo模型 - 版本历史

Matthew：/* 空间推广与应用 */

2026-10-10T09:01:31Z

空间推广与应用

Matthew：/* 神经科学意义与典型现象 */

2026-10-10T09:01:03Z

神经科学意义与典型现象

Matthew：/* 动力学机制 */

2026-10-10T09:00:07Z

动力学机制

Matthew：/* Hopf 分岔与周期发放 */

2026-10-10T08:58:15Z

Hopf 分岔与周期发放

Matthew：/* Hopf 分岔与周期发放 */

2026-10-10T08:55:55Z

Hopf 分岔与周期发放

Matthew：/* 优势、局限与其他模型的关系 */

2026-10-10T08:53:21Z

优势、局限与其他模型的关系

Matthew：/* 动力学机制 */

2026-10-10T08:49:37Z

动力学机制

2026年10月10日 (六) 08:46 Matthew

2026-10-10T08:46:58Z

2026年10月10日 (六) 08:04 Matthew

2026-10-10T08:04:42Z

Matthew：/* Hopf 分岔与周期发放 */

2026-10-10T07:09:22Z

Hopf 分岔与周期发放

@@ 第123行： / 第123行： @@
 == 空间推广与应用 ==
-单个 FHN 单元只描述随时间变化的局部状态。若要研究兴奋信号在介质中的传播，可以加入扩散项，得到反应–扩散形式：<ref name="CebrianLacasa_2024" />
+单个 FHN 单元只描述随时间变化的局部状态。若要研究兴奋信号在介质中的传播，可以加入扩散项，得到反应–扩散形式<ref name="CebrianLacasa_2024" />：
 {{NumBlk|:|<math>\begin{aligned}
@@ 第130行： / 第130行： @@
 \end{aligned}</math>|{{EquationRef|7}}}}
-式（7）中的 <math>v=v(\mathbf{r},t)</math>、<math>w=w(\mathbf{r},t)</math> 为标量场，<math>\mathbf{r}</math> 为位置向量，<math>\nabla^2</math> 为[[拉普拉斯算子]]。<math>D_v>0</math>、<math>D_w\geq0</math> 为[[扩散系数]]。神经电生理的常见简化取 <math>D_w=0</math>，即只有兴奋变量发生空间扩散。具体解还取决于参数、空间维度、初始条件和边界条件。<ref name="CebrianLacasa_2024" />
+式（7）中的 <math>v=v(\mathbf{r},t)</math>、<math>w=w(\mathbf{r},t)</math> 为标量场，<math>\mathbf{r}</math> 为位置向量，<math>\nabla^2</math> 为[[拉普拉斯算子]]。<math>D_v>0</math>、<math>D_w\geq0</math> 为[[扩散系数]]。神经电生理的常见简化取 <math>D_w=0</math>，即只有兴奋变量发生空间扩散。具体解还取决于参数、空间维度、初始条件和边界条件<ref name="CebrianLacasa_2024" />。
-在适当条件下，空间模型能够产生传播脉冲、波列和二维[[螺旋波]]。离散耦合的 FHN 单元则可用于研究神经网络中的[[同步]]及其他集体动力学。模型及其变体也用于心肌兴奋、化学可兴奋介质等问题，主要提供动力学层面的定性描述。空间脉冲的传播不要求每个未耦合单元本身处于持续振荡状态。<ref name="Nagumo_1962" /><ref name="CebrianLacasa_2024" />
+在适当条件下，空间模型能够产生传播脉冲、波列和二维[[螺旋波]]。离散耦合的 FHN 单元则可用于研究神经网络中的[[同步]]及其他集体动力学。模型及其变体也用于心肌兴奋、化学可兴奋介质等问题，主要提供动力学层面的定性描述。空间脉冲的传播不要求每个未耦合单元本身处于持续振荡状态<ref name="Nagumo_1962" /><ref name="CebrianLacasa_2024" />。
 == 优势、局限与其他模型的关系 ==

@@ 第102行： / 第102行： @@
 == 神经科学意义与典型现象 ==
 === 可兴奋响应与恢复 ===
-在可兴奋参数区域，系统具有稳定静息态。小扰动只产生较小响应，适当的较大扰动可以触发一次明显的相平面远足，随后回到静息态。恢复变量的慢变化会影响短时间内再次激发的难易程度。有限时间尺度分离下，响应振幅可随刺激连续变化，出现部分振幅的脉冲，因此不宜把 FHN 模型概括为严格的“全或无”模型。<ref name="Izhikevich_FitzHugh_2006" />
+在可兴奋参数区域，系统具有稳定静息态。小扰动只产生较小响应，适当的较大扰动可以触发一次明显的相平面远足，随后回到静息态。恢复变量的慢变化会影响短时间内再次激发的难易程度。有限时间尺度分离下，响应振幅可随刺激连续变化，出现部分振幅的脉冲，因此不宜把 FHN 模型概括为严格的“全或无”模型<ref name="Izhikevich_FitzHugh_2006" />。
 === 持续发放与[[兴奋阻滞]] ===
-短时刺激触发的单次兴奋响应，应与恒定输入支持的持续周期发放区分。某些参数下，增大恒定输入可使系统从静息进入振荡。进一步增大输入又可使系统停留在较高兴奋水平的稳定平衡态，持续发放终止，这称为兴奋阻滞。因此，输入越强并不意味着始终产生更多脉冲。<ref name="Izhikevich_FitzHugh_2006" />
+短时刺激触发的单次兴奋响应，应与恒定输入支持的持续周期发放区分。某些参数下，增大恒定输入可使系统从静息进入振荡。进一步增大输入又可使系统停留在较高兴奋水平的稳定平衡态，持续发放终止，这称为兴奋阻滞。因此，输入越强并不意味着始终产生更多脉冲<ref name="Izhikevich_FitzHugh_2006" />。
 === 参数示例 ===
-为说明式（3）和式（5）的用法，取 <math>a=0.7</math>、<math>b=0.8</math>、<math>\varepsilon=0.08</math>。这是一组常见的示例参数。<ref name="Izhikevich_FitzHugh_2006" />按本文方程计算得到以下结果，数值保留四位小数：
+为说明式（3）和式（5）的用法，取 <math>a=0.7</math>、<math>b=0.8</math>、<math>\varepsilon=0.08</math>。这是一组常见的示例参数<ref name="Izhikevich_FitzHugh_2006" />。按本文方程计算得到以下结果，数值保留四位小数：
 {| class="wikitable"

@@ 第94行： / 第94行： @@
 \end{aligned}</math>|{{EquationRef|6}}}}
-式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的[[非退化 Hopf 分岔]]，还需检查[[特征值]]实部随控制参数穿过零的[[横截条件]]及相应非线性条件。新生[[周期轨道]]是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述<ref name="CebrianLacasa_2024" />。下图是该动力学模型的[[分岔图]]。
+式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的[[非退化 Hopf 分岔]]，还需检查[[特征值]]实部随控制参数穿过零的[[横截条件]]及相应非线性条件。新生[[周期轨道]]是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述<ref name="CebrianLacasa_2024" />。下图是该动力学模型的[[分岔图]]，引用自教材<ref>Keener, J., & Sneyd, J. (2009). Mathematical physiology: II. Systems physiology (2nd ed.). Springer. https://doi.org/10.1007/978-0-387-79388-7</ref>。
 [[文件:分岔图.png|无框|507x507像素]]
 横坐标为[[外部刺激电流]] <math>I_{\mathrm{app}}</math>，纵坐标为[[兴奋变量]] <math>v</math>。实线和虚线分别表示稳定与不稳定平衡解，其中 ss 表示[[稳态解]]（steady state）；osc max 和 osc min 分别表示[[稳定周期振荡]]过程中兴奋变量的最大值和最小值；HB 表示 [[Hopf 分岔点]]。随着外部刺激电流增大，系统在第一个 Hopf 分岔点附近由[[稳定平衡态]]转入[[周期振荡态]]，并在第二个 Hopf 分岔点附近重新获得稳定平衡态。图示参数为 <math>\alpha=0.1、\gamma=0.5、\varepsilon=0.01</math>，两个 Hopf 分岔点分别位于 <math>I_{\mathrm{app}}\approx0.105</math> 和 <math>I_{\mathrm{app}}\approx1.238</math>。
 == 神经科学意义与典型现象 ==

@@ 第98行： / 第98行： @@
 [[文件:分岔图.png|无框|507x507像素]]
-横坐标为[[外部刺激电流]] \(I_{\mathrm{app}}\)，纵坐标为[[兴奋变量]] \(v\)。实线和虚线分别表示稳定与不稳定平衡解，其中 ss 表示[[稳态解]]（steady state）；osc max 和 osc min 分别表示[[稳定周期振荡]]过程中兴奋变量的最大值和最小值；HB 表示 [[Hopf 分岔点]]。随着外部刺激电流增大，系统在第一个 Hopf 分岔点附近由[[稳定平衡态]]转入[[周期振荡态]]，并在第二个 Hopf 分岔点附近重新获得稳定平衡态。图示参数为 \(\alpha=0.1\)、\(\gamma=0.5\)、\(\varepsilon=0.01\)，两个 Hopf 分岔点分别位于 \(I_{\mathrm{app}}\approx0.105\) 和 \(I_{\mathrm{app}}\approx1.238\)。
+横坐标为[[外部刺激电流]] <math>I_{\mathrm{app}}</math>，纵坐标为[[兴奋变量]] <math>v</math>。实线和虚线分别表示稳定与不稳定平衡解，其中 ss 表示[[稳态解]]（steady state）；osc max 和 osc min 分别表示[[稳定周期振荡]]过程中兴奋变量的最大值和最小值；HB 表示 [[Hopf 分岔点]]。随着外部刺激电流增大，系统在第一个 Hopf 分岔点附近由[[稳定平衡态]]转入[[周期振荡态]]，并在第二个 Hopf 分岔点附近重新获得稳定平衡态。图示参数为 <math>\alpha=0.1、\gamma=0.5、\varepsilon=0.01</math>，两个 Hopf 分岔点分别位于 <math>I_{\mathrm{app}}\approx0.105</math> 和 <math>I_{\mathrm{app}}\approx1.238</math>。
 === 快慢运动与配图 ===

@@ 第85行： / 第85行： @@
 固定 <math>w</math> 考察快子系统时，三次零倾线的外支 <math>|v|>1</math> 为吸引支，中支 <math>|v|<1</math> 为排斥支。这一结论依据快方程的导数 <math>1-v^2</math>，不能代替完整二维系统的稳定性判据。尤其是，交点位于中支并不必然意味着平衡点不稳定，也不单独保证存在稳定[[极限环]]。
-=== [[Hopf 分岔]]与周期发放 ===
+=== [[Hopf 分岔]]与[[周期发放]] ===
-当控制参数变化时，平衡点可能发生 Hopf 分岔。由式（5）的零迹和正行列式条件，得到候选位置
+当控制参数变化时，平衡点可能发生 [[Hopf 分岔]]。由式（5）的[[零迹]]和[[正行列式]]条件，得到[[分岔]]点：
 {{NumBlk|:|<math>\begin{aligned}
@@ 第94行： / 第94行： @@
 \end{aligned}</math>|{{EquationRef|6}}}}
-式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的非退化 Hopf 分岔，还需检查[[特征值]]实部随控制参数穿过零的横截条件及相应非线性条件。新生周期轨道是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述。<ref name="CebrianLacasa_2024" />
+式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的[[非退化 Hopf 分岔]]，还需检查[[特征值]]实部随控制参数穿过零的[[横截条件]]及相应非线性条件。新生[[周期轨道]]是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述<ref name="CebrianLacasa_2024" />。下图是该动力学模型的[[分岔图]]。
 [[文件:分岔图.png|无框|507x507像素]]
 === 快慢运动与配图 ===

@@ 第135行： / 第135行： @@
 == 优势、局限与其他模型的关系 ==
-模型中的恢复变量不等同于某个单独的离子通道门控概率。标准无量纲形式也不直接给出真实[[神经元]]的离子电流、[[动作电位]]持续时间或不应期长度。它可以表现某些恢复与不应性的定性特征，但定量预测仍需标定，并检验所选变体是否适合目标细胞或组织。与 HH 模型相比，FHN 保留的生物物理细节较少。<ref name="Izhikevich_FitzHugh_2006" />
+FHN 模型变量少，便于用相平面分析兴奋和恢复过程，也便于构建空间系统或耦合网络。它适合研究静息、兴奋、振荡之间的转换及其几何机制<ref name="FitzHugh_1961" />。
-从数学历史看，FHN 与范德波尔振荡器的联系体现在三次非线性、恢复动力学和弛豫振荡。从神经建模看，它与 HH 模型的联系是对兴奋行为的简化描述。仅凭其他模型同样具有兴奋或抑制变量，不能将其认定为 FHN 的直接离散化或严格等价形式。
+它可以表现某些恢复与不应性的定性特征，但定量预测仍需标定，并检验所选变体是否适合目标细胞或组织。与 HH 模型相比，FHN 保留的生物物理细节较少<ref name="Izhikevich_FitzHugh_2006" />。
 == 参考文献 ==

@@ 第61行： / 第61行： @@
 当 <math>0<b<1</math> 时，式（3）左端对 <math>v_*</math> 的导数为 <math>bv_*^2+1-b>0</math>，因此每个恒定输入对应唯一平衡点。当 <math>b>1</math> 时，某些 <math>a,I</math> 可对应三个不同平衡点。平衡点的个数本身不决定稳定性。
-零倾线并不是轨迹必须沿着运动的曲线。在零倾线上的某一瞬间，只有对应的速度分量为零。对于小 <math>\varepsilon</math>，部分轨迹会在三次零倾线的吸引支附近缓慢运动，这是一种[[快慢系统]]近似。
+零倾线并不是轨迹必须沿着运动的曲线。在零倾线上的某一瞬间，只有对应的速度分量为零。对于小 <math>\varepsilon</math>，参数支持稳定周期轨道且时间尺度分离较强，轨迹通常在三次零倾线的两条外支附近缓慢运动，并在接近折点时快速跳转，形成弛豫振荡。对式（1），快子系统的折点位于 <math>v=\pm1</math>。下图是时间序列上的图像。
 === 平衡点的局部稳定性 ===
@@ 第96行： / 第99行： @@
 === 快慢运动与配图 ===
 == 神经科学意义与典型现象 ==

@@ 第1行： / 第1行： @@
-'''菲茨休–南云模型'''（FitzHugh–Nagumo model，简称 '''FHN 模型'''）是描述可兴奋系统的二维非线性模型。它用一个快速的兴奋变量和一个较慢的恢复变量，刻画神经膜的兴奋、恢复及周期发放。作为[[霍奇金-赫胥黎模型]]的定性简化，FHN 模型保留了部分关键动力学机制，但不逐一描述钠、钾等离子通道的电流。<ref name="FitzHugh_1961">{{cite journal|author1=Richard FitzHugh|title=Impulses and Physiological States in Theoretical Models of Nerve Membrane|journal=Biophysical Journal|year=1961|volume=1|issue=6|pages=445–466|doi=10.1016/S0006-3495(61)86902-6|url=https://doi.org/10.1016/S0006-3495(61)86902-6}}</ref>
+'''菲茨休–南云模型'''（FitzHugh–Nagumo model，简称 '''FHN 模型'''）是描述可兴奋系统的二维非线性模型。它用一个快速的兴奋变量和一个较慢的恢复变量，刻画神经膜的兴奋、恢复及周期发放。作为[[霍奇金-赫胥黎模型]]的定性简化，FHN 模型保留了部分关键动力学机制，但不逐一描述钠、钾等离子通道的电流<ref name="FitzHugh_1961">{{cite journal|author1=Richard FitzHugh|title=Impulses and Physiological States in Theoretical Models of Nerve Membrane|journal=Biophysical Journal|year=1961|volume=1|issue=6|pages=445–466|doi=10.1016/S0006-3495(61)86902-6|url=https://doi.org/10.1016/S0006-3495(61)86902-6}}</ref>。
 == 基本信息 ==
@@ 第21行： / 第21行： @@
 == 历史溯源 ==
-有一种可以描述具有自激振荡特征的系统的动力学模型，名为[[范德波尔振荡器]]。1961 年，理查德·菲茨休（Richard FitzHugh）在[[范德波尔振荡器]]的基础上加入恢复项和刺激项，构造了能够表现稳定静息态、兴奋响应与重复发放的二维模型。他将其称为 Bonhoeffer–van der Pol 模型，简称 [[BVP 模型]]。这一工作的目的，是用较简单的数学结构分析神经膜的兴奋与恢复，而非证明具体离子通道对真实神经活动不重要。<ref name="FitzHugh_1961" />
+有一种可以描述具有自激振荡特征的系统的动力学模型，名为[[范德波尔振荡器]]。1961 年，理查德·菲茨休（Richard FitzHugh）在[[范德波尔振荡器]]的基础上加入恢复项和刺激项，构造了能够表现稳定静息态、兴奋响应与重复发放的二维模型。他将其称为 Bonhoeffer–van der Pol 模型，简称 [[BVP 模型]]。这一工作的目的，是用较简单的数学结构分析神经膜的兴奋与恢复，而非证明具体离子通道对真实神经活动不重要<ref name="FitzHugh_1961" />。
-年，南云仁一（J. Nagumo）、S. Arimoto 和 S. Yoshizawa 用[[隧道二极管]]等元件构建了模拟神经轴突的[[有源脉冲传输线]]。该电路提供了相关兴奋动力学及脉冲传播的物理实现，后来这一模型通常被称为 FitzHugh–Nagumo（FHN）模型。<ref name="Nagumo_1962">{{cite journal|author1=J. Nagumo|author2=S. Arimoto|author3=S. Yoshizawa|title=An Active Pulse Transmission Line Simulating Nerve Axon|journal=Proceedings of the IRE|year=1962|volume=50|issue=10|pages=2061–2070|doi=10.1109/JRPROC.1962.288235|url=https://doi.org/10.1109/JRPROC.1962.288235}}</ref>
+年，南云仁一（J. Nagumo）、S. Arimoto 和 S. Yoshizawa 用[[隧道二极管]]等元件构建了模拟神经轴突的[[有源脉冲传输线]]。该电路提供了相关兴奋动力学及脉冲传播的物理实现，后来这一模型通常被称为 FitzHugh–Nagumo（FHN）模型<ref name="Nagumo_1962">{{cite journal|author1=J. Nagumo|author2=S. Arimoto|author3=S. Yoshizawa|title=An Active Pulse Transmission Line Simulating Nerve Axon|journal=Proceedings of the IRE|year=1962|volume=50|issue=10|pages=2061–2070|doi=10.1109/JRPROC.1962.288235|url=https://doi.org/10.1109/JRPROC.1962.288235}}</ref>。
 在早期研究中，FHN模型主要用于解释神经元的兴奋阈值、动作电位产生、不应期以及神经冲动传播。后来，它从单个神经元动力学模型，逐渐成为研究一般兴奋介质的理论模型。
@@ 第45行： / 第45行： @@
 以下平衡点和[[分岔]]分析假定输入恒定，即 <math>I(t)=I</math>。我们可以定义两条[[零倾线]]，即某一变量的时间导数为零。
-[[文件:相平面.png|thumb|left|520px|图1：原词条相平面示意图，图中标注恒定输入 <math>I=0.8</math>。横轴为兴奋变量 <math>v</math>，纵轴为恢复变量 <math>w</math>。红色三次曲线和蓝色直线分别为式（2）的两条零倾线，灰色箭头表示向量场，绿色闭合轨迹表示图示参数下的周期运动。图中的恢复零倾线使用 <math>w=(v+a)/b</math>，与本文约定一致。“中支不稳定”指固定恢复变量时快子系统的性质。图源：站内原文件，由 Matthew 上传。]]
+[[文件:相平面.png|thumb|left|520px|图1：相平面示意图，图中标注恒定输入 <math>I=0.8</math>。横轴为兴奋变量 <math>v</math>，纵轴为恢复变量 <math>w</math>。红色三次曲线和蓝色直线分别为式（2）的两条零倾线，灰色箭头表示向量场，绿色闭合轨迹表示图示参数下的周期运动。图中的恢复零倾线使用 <math>w=(v+a)/b</math>，与本文约定一致。“中支不稳定”指固定恢复变量时快子系统的性质。]]
 <div style="clear:both;"></div>

@@ 第42行： / 第42行： @@
 == 动力学机制 ==
 === [[零倾线]]与[[平衡点]] ===
-以下平衡点和[[分岔]]分析假定输入恒定，即 <math>I(t)=I</math>。我们可以定义两条[[零倾线]]，即某一变量的时间导数为零。由式（1）得到
 {{NumBlk|:|<math>\begin{aligned}
@@ 第92行： / 第98行： @@
 当参数支持稳定周期轨道且时间尺度分离较强时，轨迹通常在三次零倾线的两条外支附近缓慢运动，并在接近折点时快速跳转，形成弛豫振荡。对式（1），快子系统的折点位于 <math>v=\pm1</math>。有限 <math>\varepsilon</math> 下的实际轨迹及跳转位置不必与零倾线或折点完全重合。
 [[文件:时间响应.png|thumb|left|760px|图2：原词条时间响应示意图。横轴为无量纲时间 <math>t</math>，蓝线为 <math>v(t)</math>，红线为 <math>w(t)</math>。左图标注恒定输入 <math>I=0.1</math>，显示轨迹趋于静息态。右图标注恒定输入 <math>I=0.8</math>，显示持续振荡。这是两个恒定输入情形的比较，不能直接解释为短时阈下、阈上脉冲的比较。图源：站内原文件，由 Matthew 上传。]]

←上一版本		2026年10月10日 (六) 07:09的版本
第86行：		第86行：

	式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的非退化 Hopf 分岔，还需检查[[特征值]]实部随控制参数穿过零的横截条件及相应非线性条件。新生周期轨道是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述。<ref name="CebrianLacasa_2024" />		式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的非退化 Hopf 分岔，还需检查[[特征值]]实部随控制参数穿过零的横截条件及相应非线性条件。新生周期轨道是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述。<ref name="CebrianLacasa_2024" />
		+
		+	[[文件:分岔图.png\|无框\|507x507像素]]

	=== 快慢运动与配图 ===		=== 快慢运动与配图 ===