LIF神经元模型

clear
clf
% membrane constants
tau_m = 0.030;
Rm = 3e7;
thresh=0.05;
t_spikes=[];
No_spikes=0;
in_abs_ref=0;
abs_ref_count=0;
% time
delta_t = 0.001; % integration time step
duration = 0.5;% duration of the simulation
abs_T=0.003;%absolute refractory period绝对不应期
abs_ref_steps=round（abs_T./delta_t）;%四舍五入
No_steps = round（duration./delta_t）;
times = linspace（0 duration No_steps + 1）; % the time points for update
% membrane potential
V_0= 0.03;
V = zeros（1 No_steps +1）;
V（1）= V_0;

%% injection current
I_0 =2e-9;
I_duration = 0.25;
I_start= 0.05;
I_start_index= round（I_start./ delta_t）;
I_finish_index = round（（I_start+I_duration）./delta_t）;
Iinj=zeros（1No_steps+1）;
Iinj（I_start_index: I_finish_index）= I_0;
%% noise current
std_dev_noise = 1.5e-10;
I_noise