Please cite "Y. Burak and I. R. Fiete. Accurate path integration in continuous attractor network models of grid cells. PLoS Comp. Biol. 5(2) (2009)" If this code is used.
Notion 笔记：https://www.notion.so/Burak-Fiete-2009-280ca7f3face81f2ba8cd04dee126334?source=copy_link
- gc_dynamics.m
- gc_periodic.m
- gc_non_periodic.m
- gc_dynamics_1d.m
The associated code models periodic and aperiodic continuous attractor networks of Grid Cells with random or recorded trajectory data. Neuron spiking is simulated by a Poisson Process. For the 2d model a single neuron response for the simulation is recorded and returned. During the simulation, the network population response and single neuron spiking with respect to the trajectory are displayed in figures.
1d Simulation -
- To run the simulation call gc_dynamics_1d(periodic,n) from Matlab or Octave where periodic is 0 or 1 and n is the population size. ex. gc_dynamics_1d(0,128) simulates an aperiodic network with 128 neurons.
2d Simulation -
-
To run the simulation call gc_dynamics('filename', dt, useSpiking, simulate_aperiodic) from Matlab or Octave where : filename is your trajectory filename. If you're not using a trajectory, then use ''. dt is your recorded timestep in ms. If you're not using a trajectory, then use 0. useSpiking is {0,1}. simulate_aperiodic is {0,1}
ex. gc_dynamics('test.mat', .5, 1,0) Simulates a periodic network with spiking that reads data from test.mat at a timestep of .5 ms.
-
If you are using trajectory data, move it to the same directory as gc_dynamics.m. It must be a .mat file with two columns: position_x in centemeters position_y in centemeters you must change dt to your recording timestep, measured in milliseconds.
-
If spiking is used, you will have a variable named 'spikes' in your environment after the simulation which is the spiking of a random single neuron.
- The simulation places the figures at specific locations on the screen. If the figures are not being displayed properly, try adjusting 'set(fig,'Position',[50,1000,450,900])' in gc_periodic.m and gc_non_periodic.m