SNN脉冲神经网络中的IF神经元脉冲仿真-源码

%V_reset =0 value to reset voltage to after a spike [mV] %t_ref %refractory time constant [ms] function AveRate_vect = IF_UFun(V_res,t_ref) %DEFINE PARAMETERS dt = 0.1; %time step [ms] t_time =1000; %total time of run [ms] V_th =10; %spike threshold [mV] R_m =10; %membrane resistance [MOhm] tau =10; %membrane time constant [ms] V0 = 1; %resting membrane potential [mV] %t_ref=5; %refractory time constant [ms] %DEFINE INITIAL VALUES AND VECTORS TO HOLD RESULTS t_vect = 0:dt:t_time; %will hold vector of times V_vect = zeros(1,length(t_vect)); %initialize the voltage vector %INTEGRATE THE EQUATION tau*dV/dt = -V + E_L + I_e*R_m I_Num=0; I_vect_long = 1.00000000001:0.0001:3; %vector of injected current for producing theory plot AveRate_vect= zeros(0,length(I_vect_long)); for I_In =I_vect_long; %loop over different I_In values I_Num = I_Num + 1; NumSpikes = 0; %holds number of spikes that have occurred i = 1; % index denoting which element of V is being assigned spike_flag=0; Times=t_ref/dt;%max times V_vect(i)= V0; %first element of V, i.e. value of V at t=0 I_e_vect = I_In*ones(1,1+t_time/dt); for t=dt:dt:t_time %loop through values of t in steps of dt ms %INTEGRATE THE EQUATION tau*dV/dt = -V + IR RI =I_e_vect(i)*R_m; if (V_vect(i) > V_th) %cell spiked spike_flag=1; Times=0; NumSpikes = NumSpikes + 1; %add 1 to the total spike count end if (spike_flag==1) V_vect(i+1) = V_res; %set voltage back to V_res Times=Times+1; else V_vect(i+1) = RI*dt/tau + (1-dt/tau)*V_vect(i); end if (Times==(t_ref/dt+1)) spike_flag=0; end i = i + 1; %add 1 to index, corresponding to moving forward 1 time step end AveRate_vect(I_Num) =NumSpikes/t_time; %gives average firing rate in [#/sec = Hz] end end