| 1 | ''' |
|---|
| 2 | A very simple example Brian script to show how to implement |
|---|
| 3 | a leaky integrate and fire model. In this example, we also |
|---|
| 4 | drive the single leaky integrate and fire neuron with |
|---|
| 5 | regularly spaced spikes from the :class:`SpikeGeneratorGroup`. |
|---|
| 6 | ''' |
|---|
| 7 | from brian import * |
|---|
| 8 | |
|---|
| 9 | tau = 10*ms |
|---|
| 10 | Vr = -70*mV |
|---|
| 11 | Vt = -55*mV |
|---|
| 12 | |
|---|
| 13 | model = Model(equations=''' |
|---|
| 14 | dV/dt = -(V-Vr)/tau : volt |
|---|
| 15 | ''', threshold=Vt, reset=Vr) |
|---|
| 16 | G = NeuronGroup(1,model) |
|---|
| 17 | |
|---|
| 18 | spikes = linspace(10*ms,100*ms,25) |
|---|
| 19 | input = MultipleSpikeGeneratorGroup([spikes]) |
|---|
| 20 | |
|---|
| 21 | C = Connection(input, G) |
|---|
| 22 | C[0,0] = 5*mV |
|---|
| 23 | |
|---|
| 24 | M = StateMonitor(G, 'V', record=True) |
|---|
| 25 | |
|---|
| 26 | G.V = Vr |
|---|
| 27 | run(100*ms) |
|---|
| 28 | plot(M.times/ms,M[0]/mV) |
|---|
| 29 | show() |
|---|
