def ousim(mu_amp, mu_offs, sigma_amp, sigma_offs, freq, V_th): # mu_amp, mu_offs, sigma_amp, sigma_offs, freq, V_th = config if sigma_amp > sigma_offs: sigma_amp = sigma_offs # print("Setting up OU LIF simulation...") ounet = Network() clock.reinit_default_clock() eqs =Equations('dV/dt = mu-(V+V0)/tau + sigma*I/sqrt(dt) : volt') eqs+=Equations('dI/dt = -I/dt + xi/sqrt(dt) : 1') eqs+=Equations('mu = mu_amp*sin(t*freq*2*pi) + mu_offs : volt/second') eqs+=Equations('sigma = sigma_amp*sin(t*freq*2*pi) + sigma_offs :' ' volt/sqrt(second)') eqs.prepare() ounrn = NeuronGroup(1, eqs, threshold=V_th, refractory=t_refr, reset=V_reset) ounet.add(ounrn) ounrn.V = V0 V_mon = StateMonitor(ounrn, 'V', record=True) st_mon = SpikeMonitor(ounrn) ounet.add(V_mon, st_mon) ounet.run(duration) V_mon.insert_spikes(st_mon, value=V_th*2) times = V_mon.times membrane = V_mon[0] return times, st_mon.spiketimes[0], membrane
def lifsim(mu_amp, mu_offs, simga_amp, sigma_offs, freq, V_th): lifnet = Network() clock.reinit_default_clock() eqs = Equations('dV/dt = (-V+V0)/tau : volt') eqs.prepare() lifnrn = NeuronGroup(1, eqs, threshold=V_th, refractory=t_refr, reset=V_reset) lifnet.add(lifnrn) pulse_times = (np.arange(1, duration*freq, 1)+0.25)/freq pulse_spikes = [] Npoiss = 5000 Npulse = 5000 wpoiss = (mu_offs-mu_amp)/(Npoiss*freq) wpulse = mu_amp/(Npulse*freq) sigma = 1/(freq*5) if (wpulse != 0): for pt in pulse_times: pp = PulsePacket(t=pt*second, n=Npulse, sigma=sigma) pulse_spikes.extend(pp.spiketimes) pulse_input = SpikeGeneratorGroup(Npulse, pulse_spikes) pulse_conn = Connection(pulse_input, lifnrn, 'V', weight=wpulse) lifnet.add(pulse_input, pulse_conn) if (wpoiss != 0): poiss_input = PoissonGroup(Npoiss, freq) poiss_conn = Connection(poiss_input, lifnrn, 'V', weight=wpoiss) lifnet.add(poiss_input, poiss_conn) V_mon = StateMonitor(lifnrn, 'V', record=True) st_mon = SpikeMonitor(lifnrn) lifnet.add(V_mon, st_mon) lifnet.run(duration) V_mon.insert_spikes(st_mon, value=V_th*2) times = V_mon.times membrane = V_mon[0] return times, st_mon.spiketimes[0], membrane
nkreuzsamples = 3 Vrest = 0*mV Vth = 20*mV tau = 20*ms Nnrns = 4 Ningroups = 1 Nin_per_group = 50 fin = 20*Hz ingroup_sync = [0.5] sigma = 0*ms weight = 2.0*mV Nallin = Nin_per_group*Ningroups Nin = 25 # total number of connections each cell receives lifeq_exc = Equations("dV/dt = (Vrest-V)/tau : volt") lifeq_exc.prepare() nrngroup = NeuronGroup(Nnrns, lifeq_exc, threshold="V>Vth", reset=Vrest, refractory=2*ms) nrngroup.V = Vrest network.add(nrngroup) print("Setting up inputs and connections ...") ingroups = [] inpconns = [] for ing in range(Ningroups): ingroup = sl.tools.fast_synchronous_input_gen(Nin_per_group, fin, ingroup_sync[ing], sigma, duration, shuffle=False) inpconn = Connection(ingroup, nrngroup, 'V') ingroups.append(ingroup) inpconns.append(inpconn) inputneurons = []