FitzHugh-Nagumo 模型连续-MATLAB仿真

1. FitzHugh-Nagumo 模型微分方程、

function dYdt = fhn_model(Y, a, b, c, I)
    V = Y(1);
    W = Y(2);
    dVdt = V - V^3 / 3 - W + I; 
    %表示膜电位的变化率。这个方程包含了两个部分，第一部分 V - V^3 / 3 描述了膜电位的变化，第二部分 - W + I 描述了外部刺激对膜电位的影响。
    dWdt = c * (V + a - b*W);
 %表示激活变量的变化率。这个方程也包含了两个部分，第一部分 V + a 描述了膜电位对激活变量的影响，第二部分 - b*W 描述了激活变量自身的耗散，其中 b 是耗散系数，c 是一个比例系数。
    dYdt = [dVdt; dWdt];
end

2. FitzHugh-Nagumo 模型脉冲响应微分方程

    function dYdt = fhn_model_impulse(Y, a, b, c, I)
        V = Y(1);
        W = Y(2);    
        dVdt = V - V^3 / 3 - W + I;
        dWdt = c * (V + a - b*W);
        dYdt = [dVdt; dWdt];
        if rand < 0.01 % 添加脉冲激励
            dYdt = dYdt + [1; 0];
        end
    end
    %该模型模拟神经元的放电行为。

3. MATLAB仿真测试

clc;clear
% FitzHugh-Nagumo 模型参数