关于Izhikevich神经元模型的MATLAB代码

1. 在神经元模型中，Izhikevich模型相较于HH模型和IF模型，具有准确度高（相比IF），计算简单（相比HH）等优点，其数学表达式如下：

v‘ =0.04v^2 + 5v + 140 - u + I 
u' =a(bv - u) 

if v > 30 mV;

then
v=c
u=u + d:

现在网络上和各大论坛所有的大部分都是多神经元模型，有关于单个神经元模型的代码很少，这不利于我们对于其特性的探索。

因此，我将附上Izhikevich的单神经元matlab代码

具体代码如下：

a=0.02;
b=0.2;
c=-65;
d=8;
v=-65;
u=b*v;
dt=0.001;
t=0:dt:0.15;
n=length(t);
V=[];
U=[];
for i=1:n
    I=50;
    v=0.5*(0.04*v^2+5*v+140-u+I);
    u=u+a*(b*v-u);
    V(i)=v;
    U(i)=u;
    if v>=30
        v=c;
        u=u+d;
    else
        v=v;
        u=u;
    end
end
    plot(t,V);

所得最终结果如下所示：

。。。弄不上去，总之能跑出来就对了