if_const.m

用matlab实现模拟神经元模型

% integrate-and-fire neuron
clear; clf;

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