if_vary1.m

用matlab实现模拟神经元模型

clear; clf;
%parameter
tau=10;         %ms
Rm=10;          %K Ohs*cm2
Vth=-50;        %mV
El=-65;         %mV
Vreset=-65;     %mV
% construct noisy external input
Ie=-0.5+4*rand(100,1);      %uA/cm2
%initial membrane voltage
V=-65;          %mV
%integration
for t=1:1:100
    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);