# 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).

In [1]:

import matplotlib.pyplot as plt
%matplotlib inline

import nengo


## Step 1: Create the neural populations¶

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

In [2]:

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(100, dimensions=1)


## 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 [3]:

# Create a piecewise step function for input
from nengo.processes import Piecewise
with model:
input = nengo.Node(
Piecewise({
0: 0,
0.2: 1,
1: 0,
2: -2,
3: 0,
4: 1,
5: 0
}))


## Step 3: Connect the network elements¶

In [4]:

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 [5]:

with model:
input_probe = nengo.Probe(input)
A_probe = nengo.Probe(A, synapse=0.01)


## Step 5: Run the model¶

In [6]:

# Create our simulator
with nengo.Simulator(model) as sim:
# Run it for 6 seconds
sim.run(6)


## Step 6: Plot the results¶

In [7]:

# Plot the decoded output of the ensemble
plt.figure()
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.