Nengo Example: Integrator

This demo implements a one-dimensional neural integrator.

This is the first example of a recurrent network in the demos. It shows how neurons can be used to implement stable dynamics. Such dynamics are important for memory, noise cleanup, statistical inference, and many other dynamic transformations.

When you run this demo, it will automatically put in some step functions on the input, so you can see that the output is integrating (i.e. summing over time) the input. You can also input your own values. Note that since the integrator constantly sums its input, it will saturate quickly if you leave the input non-zero. This makes it clear that neurons have a finite range of representation. Such saturation effects can be exploited to perform useful computations (e.g. soft normalization).

Step 1: Create the neural populations

Our model consists of one recurrently connected ensemble, and an input population.


In [ ]:
import nengo

model = nengo.Network(label='Integrator')
with model:
    # Our ensemble consists of 100 leaky integrate-and-fire neurons,
    # representing a one-dimensional signal
    A = nengo.Ensemble(nengo.LIF(100), dimensions=1, label='Integrator')

Step 2: Create input for the model

We will use a piecewise step function as input, so we can see the effects of recurrence.


In [ ]:
# Create a piecewise step function for input
from nengo.utils.functions import piecewise
with model:
    input = nengo.Node(piecewise({0: 0, 0.2: 1, 1: 0, 2: -2, 3: 0, 4: 1, 5: 0}), label='Piecewise input')

Step 3: Connect the network elements


In [ ]:
with model:
    # Connect the population to itself
    tau = 0.1
    nengo.Connection(A, A, transform=[[1]], synapse=tau) # Using a long time constant for stability
    
    # Connect the input
    nengo.Connection(input, A, transform=[[tau]], synapse=tau) # The same time constant as recurrent to make it more 'ideal'

Step 4: Probe outputs

Anything that is probed will collect the data it produces over time, allowing us to analyze and visualize it later.


In [ ]:
with model:
    # Add probes
    input_probe = nengo.Probe(input, 'output')
    A_probe = nengo.Probe(A, 'decoded_output', synapse=0.01)

Step 5: Run the model


In [ ]:
# Create our simulator
sim = nengo.Simulator(model)
# Run it for 6 seconds
sim.run(6)

Step 6: Plot the results


In [ ]:
import matplotlib.pyplot as plt

# Plot the decoded output of the ensemble
plt.plot(sim.trange(), sim.data[input_probe], label="Input")
plt.plot(sim.trange(), sim.data[A_probe], 'k', label="Integrator output")
plt.legend()

The graph shows the response to the input by the integrator. Because it is implemented in neurons, it will not be perfect (i.e. there will be drift). Running several times will give a sense of the kinds of drift you might expect. Drift can be reduced by increasing the number of neurons.