Mean-Field of 2D Spiking Neurons

We consider the following McKean-Vlasov equation: dVtμ=b(Vtμ)dt+JE[Xtμf(Vtμ)]dtVtμR+1{zf(Vtμ)}N(dt,dz),dXtμ=1XtμτdtUXtμR+1{zf(Vtμ)}N(dt,dz), where b(v)=mλv, and f(v)=v2 or f(v)=v1/β or f(v)=1β1v1.

Simulation & Initial Conditions

0501001502002503003504000200400
Neural Raster Plot (U = 0.1, J = 3.1177)Time (s)Neuron Index
05010015020025030035040001234
Interaction over time: J E[X f(V)]Time (s)ValueInteractionMean Rate
00.511.522.500.0010.0020.003
Power Spectrum of Mean Rate E[f(V_t)] (Peak: f = 0.01 Hz, ω = 0.08 rad/s, T = 81.92 s)Angular Frequency ω (rad/s)Power Spectral DensityFrequency f (Hz)Angular Freq ω

Calculate Invariant Measures

Interactive: Click anywhere on the computed curve to automatically select an α value and jump to Stability Analysis!

0.050.10.150.20.250.30.350.40.450.500.511.522.533.5
1D model2D modelInvariant Measures (Click curve to select α)α (alpha)J (Coupling)

Stability Integral

Computing Complex ODEs...