Izhikevich

代码说明：

U(1)=round(u0);dV(1)=0;dU(1)=0;I=18*ones(1,1500);I(1:300)=0;I(1201:1500)=0;               %Initialization of current II=I+(40+b)*ones(1,1500);  for k=2:1500    V(k)=V(k-1)+round((22/(512*512))*V(k-1)*V(k-1)-(20/32)*V(k-1)-(1/(512))*U(k-1)+20*I(k-1));    U(k)=U(k-1)+round(a*b*(1/512)*V(k-1)-a*(1/512)*U(k-1));    dV(k)=round(2*((1/512)*(1/512)*V(k-1)*V(k-1))-15*V(k-1)-(26/128)*U(k-1)+200*I(k-1));    dU(k)=round(a*b*(26/128)*V(k-1)-a*(1/512)*U(k-1));