leaky integrate and fire MATLAB model

LIF is considered as one of the simpler spiking neuron models, and I’ve been doing a bit of work with it in MATLAB:

Leaky intergrate-and-fire graph

Leaky intergrate-and-fire graph

The graph shows the membrane potential of a neuron over time. There’s an inhibitory current injection at 0.04 seconds, followed by an excitatory injection at 0.05 and 0.09 seconds. The second current injection causes the voltage to exceed a set threshold (shown by the red line), and so an activation potential is generated (the neuron spikes). A final excitatory injection is given at 0.11 seconds, however this is within the time of the absolute refractory period: although the voltage crosses the threshold again, the neuron doesn’t spike.

The code is over on my GitHub.

Feng, J. (2001). Is the integrate-and-fire model good enough?–a review. Neural networks : the official journal of the International Neural Network Society, 14(6-7), 955–75.