def question_external_poisson_population():
rate_monitor_excit, spike_monitor_excit, voltage_monitor_excit, idx_monitored_neurons_excit, rate_monitor_inhib, \
spike_monitor_inhib, voltage_monitor_inhib, idx_monitored_neurons_inhib, w_profile = \
wm_model.simulate_wm(
sim_time=800. * b2.ms, poisson_firing_rate=2.2 * b2.Hz, sigma_weight_profile=20., Jpos_excit2excit=1.6)
plot_tools.plot_network_activity(rate_monitor_excit,
spike_monitor_excit,
voltage_monitor_excit,
t_min=0. * b2.ms)
plt.show()
def question_role_of_inhib_population():
rate_monitor_excit, spike_monitor_excit, voltage_monitor_excit, idx_monitored_neurons_excit, rate_monitor_inhib, spike_monitor_inhib, voltage_monitor_inhib, idx_monitored_neurons_inhib, w_profile = \
wm_model.simulate_wm(N_excitatory=1024, N_inhibitory=1, sigma_weight_profile=20,
stimulus_center_deg=120, stimulus_width_deg=30, stimulus_strength=0.5 * b2.namp,
t_stimulus_start=100 * b2.ms,
t_stimulus_duration=100 * b2.ms, sim_time=500. * b2.ms, Jpos_excit2excit=1.6,
poisson_firing_rate=1.5 * b2.Hz)
fig, ax_raster, ax_rate, ax_voltage = plot_tools.plot_network_activity(
rate_monitor_excit,
spike_monitor_excit,
voltage_monitor_excit,
t_min=0. * b2.ms)
plt.show()
def question_weight_profile():
rate_monitor_excit, spike_monitor_excit, voltage_monitor_excit, idx_monitored_neurons_excit, rate_monitor_inhib, \
spike_monitor_inhib, voltage_monitor_inhib, idx_monitored_neurons_inhib, weight_profile_45 = \
wm_model.simulate_wm(
sim_time=800. * b2.ms, poisson_firing_rate=2.3 * b2.Hz, sigma_weight_profile=5., Jpos_excit2excit=6)
plot_tools.plot_network_activity(rate_monitor_excit,
spike_monitor_excit,
voltage_monitor_excit,
t_min=0. * b2.ms)
plt.show()
plt.figure()
plt.plot(weight_profile_45)
plt.show()
def distractor_at_same_time():
rate_monitor_excit, spike_monitor_excit, voltage_monitor_excit, idx_monitored_neurons_excit, rate_monitor_inhib, spike_monitor_inhib, voltage_monitor_inhib, idx_monitored_neurons_inhib, w_profile = \
wm_model.simulate_wm(N_excitatory=1024, N_inhibitory=256, sigma_weight_profile=60,
stimulus_center_deg=40, stimulus_width_deg=20, stimulus_strength=0.5 * b2.namp,
t_stimulus_start=100 * b2.ms,G_excit2inhib=0.355*b2.nS,
t_stimulus_duration=100 * b2.ms, sim_time=500. * b2.ms, Jpos_excit2excit=1.6,
poisson_firing_rate=1.5 * b2.Hz, distractor_center_deg=300, distractor_width_deg=20,
distractor_strength=0.5 * b2.namp,
t_distractor_start=110 * b2.ms, t_distractor_duration=100 * b2.ms)
fig, ax_raster, ax_rate, ax_voltage = plot_tools.plot_network_activity(
rate_monitor_excit,
spike_monitor_excit,
voltage_monitor_excit,
t_min=0. * b2.ms)
plt.show()
def question_integration_of_input():
rate_monitor_excit, spike_monitor_excit, voltage_monitor_excit, idx_monitored_neurons_excit, rate_monitor_inhib, spike_monitor_inhib, voltage_monitor_inhib, idx_monitored_neurons_inhib, w_profile = \
wm_model.simulate_wm(
stimulus_center_deg=120, stimulus_width_deg=60, stimulus_strength=0.5 * b2.namp,
t_stimulus_start=100 * b2.ms,
t_stimulus_duration=200 * b2.ms, sim_time=500. * b2.ms)
fig, ax_raster, ax_rate, ax_voltage = plot_tools.plot_network_activity(
rate_monitor_excit,
spike_monitor_excit,
voltage_monitor_excit,
t_min=0. * b2.ms)
plt.show()
rate_monitor_excit, spike_monitor_excit, voltage_monitor_excit, idx_monitored_neurons_excit, rate_monitor_inhib, spike_monitor_inhib, voltage_monitor_inhib, idx_monitored_neurons_inhib, w_profile = \
wm_model.simulate_wm(
stimulus_center_deg=120, stimulus_width_deg=30, stimulus_strength=0.5 * b2.namp,
t_stimulus_start=100 * b2.ms,
t_stimulus_duration=200 * b2.ms, sim_time=500. * b2.ms)
fig, ax_raster, ax_rate, ax_voltage = plot_tools.plot_network_activity(
rate_monitor_excit,
spike_monitor_excit,
voltage_monitor_excit,
t_min=0. * b2.ms)
plt.show()
def test_woking_memory_sim():
"""Test if the working memory circuit is initialized and simulated for 1ms"""
b2.defaultclock.dt = 0.2 * b2.ms
wm_model.simulate_wm(N_excitatory=40, N_inhibitory=10, sim_time=1. * b2.ms)