if_vary2.m

来自「用matlab实现模拟神经元模型」· M 代码 · 共 26 行

%parameter
tau=10;         %ms
Rm=10;          %K Ohs*cm2
Vth=-50;        %mV
El=-65;         %mV
Vreset=-65;     %mV
% constant external input
load vary_Ie        %uA/cm2
%initial membrane voltage
V=-65;          %mV
%integration
for t=1:1:500
    V=El+Rm*Ie(t)+(V-El-Rm*Ie(t))*exp(-1/tau);
    if( V>=Vth)
        plot_V(t)=0;
        V=Vreset;
    else
        plot_V(t)=V;
    end
    plot_t(t)=t;
end
subplot(2,1,1);
plot(plot_t,plot_V);
subplot(2,1,2);
plot(plot_t,Ie);