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 = []