Exemplo n.º 1
0
 def test_pair(self, pre_duration=2.0, stimulus_duration=10.0, post_duration=2.0, delay=0.01):
   total_duration = pre_duration + stimulus_duration + post_duration
   stimulus_begin = pre_duration
   stimulus_end = pre_duration + stimulus_duration
   
   startTime = time()
   timeNow = 0.0
   
   n1 = Neuron((0.0, 0.0, 0.0), timeNow)
   n2 = Neuron((0.0, 0.0, 1.0), timeNow)
   n1.synapseWith(n2)
   
   # Set up plotting
   figure(figsize = (12, 9))
   hold(True)  # [graph]
   subplot(211)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.02)  # [graph]
   setup_neuron_plot("Presynaptic neuron", None, "Membrane potential (V)")
   subplot(212)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.02)  # [graph]
   setup_neuron_plot("Postsynaptic neuron", "Time (s)", "Membrane potential (V)")
   subplots_adjust(hspace = 0.33)
   
   while timeNow <= total_duration:
     timeNow = time() - startTime
     if stimulus_begin <= timeNow <= stimulus_end:
       n1.accumulate(np.random.normal(0.0035, 0.0005))  # TODO accumulate value based on deltaTime (since last accumulate)
     n1.update(timeNow)
     #print n1.id, n1.timeCurrent, n1.potential  # [log: potential]
     n2.update(timeNow)
     #print n2.id, n2.timeCurrent, n2.potential  # [log: potential]
     
     subplot(211)  # [graph]
     n1.plot()  # [graph]
     subplot(212)  # [graph]
     n2.plot()  # [graph]
     pause(delay)  # [graph]
     #sleep(delay)
   
   if holdPlots:
     show()  # [graph]
Exemplo n.º 2
0
 def test_single(self, duration=10.0, delay=0.01):
   startTime = time()
   timeNow = 0.0
   n = Neuron((0.0, 0.0, 0.0), timeNow)
   
   # Set up plotting
   figure(figsize = (12, 9))
   hold(True)  # [graph]
   xlim(0.0, duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.02)  # [graph]
   setup_neuron_plot("Stimulated neuron", None, "Membrane potential (V)")
   
   while timeNow <= duration:
     timeNow = time() - startTime
     n.accumulate(np.random.normal(0.001, 0.00025))
     n.update(timeNow)
     n.plot()  # [graph]
     pause(delay)  # [graph]
     #sleep(delay)
   
   if holdPlots:
     # Show plot (and wait till it's closed)
     show()  # [graph]
Exemplo n.º 3
0
 def test_gatekeeper(self, total_duration=14.0, stimulus_period=(2.0, 12.0), gate_period=(5.0, 8.0), delay=0.01):
   startTime = time()
   timeNow = 0.0
   
   n1 = Neuron((0.0, 0.0, 0.0), timeNow)
   n2 = Neuron((0.0, 0.0, 1.0), timeNow)
   g1 = Neuron((-1.0, -1.0, 1.0), timeNow)  # gatekeeper
   n1.synapseWith(n2, None, g1)  # auto-initialize synaptic strength
   
   # Set up plotting
   figure(figsize = (12, 9))
   hold(True)  # [graph]
   ax = subplot(311)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.01)  # [graph]
   setup_neuron_plot("Presynaptic neuron", None, None)
   ax.get_xaxis().set_ticklabels([])
   ax = subplot(312)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.01)  # [graph]
   setup_neuron_plot("Postsynaptic neuron", None, "Membrane potential (V)")
   ax.get_xaxis().set_ticklabels([])
   subplot(313)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.01)  # [graph]
   setup_neuron_plot("Gatekeeper neuron", "Time (s)", None)
   subplots_adjust(hspace = 0.33)
   
   while timeNow <= total_duration:
     timeNow = time() - startTime
     
     if stimulus_period[0] <= timeNow <= stimulus_period[1]:
       n1.accumulate(np.random.normal(0.004, 0.0005))  # TODO accumulate value based on deltaTime (since last accumulate)
     
     if gate_period[0] <= timeNow <= gate_period[1]:
       g1.accumulate(np.random.normal(0.0035, 0.0005))
     
     n1.update(timeNow)
     #print n1.id, n1.timeCurrent, n1.potential  # [log: potential]
     n2.update(timeNow)
     #print n2.id, n2.timeCurrent, n2.potential  # [log: potential]
     g1.update(timeNow)
     
     subplot(311)  # [graph]
     n1.plot()  # [graph]
     subplot(312)  # [graph]
     n2.plot()  # [graph]
     subplot(313)  # [graph]
     g1.plot()  # [graph]
     pause(delay)  # [graph]
     #sleep(delay)
   
   if holdPlots:
     show()  # [graph]
Exemplo n.º 4
0
 def test_inhibitor(self, total_duration=14.0, stimulus_period=(2.0, 12.0), inhibition_period=(5.0, 8.0), delay=0.01):
   startTime = time()
   timeNow = 0.0
   
   n1 = Neuron((0.0, 0.0, 0.0), timeNow)
   n2 = Neuron((0.0, 0.0, 1.0), timeNow)
   i1 = Neuron((-1.0, -1.0, 1.0), timeNow)  # inhibitor
   n1.synapseWith(n2)  # no synaptic gating
   i1.synapseWith(n2, -np.random.normal(synaptic_strength.mu, synaptic_strength.sigma))  # inhibitory synapse
   
   # Set up plotting
   figure(figsize = (12, 9))
   hold(True)  # [graph]
   subplot(311)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.01)  # [graph]
   setup_neuron_plot("Neuron " + str(n1.id))
   subplot(312)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.01)  # [graph]
   setup_neuron_plot("Neuron " + str(n2.id))
   subplot(313)  # [graph]
   xlim(0.0, total_duration)  # [graph]
   ylim(action_potential_trough.mu - 0.01, action_potential_peak + 0.01)  # [graph]
   setup_neuron_plot("Neuron " + str(i1.id) + " (inhibitor)")
   subplots_adjust(hspace = 0.33)
   
   while timeNow <= total_duration:
     timeNow = time() - startTime
     
     if stimulus_period[0] <= timeNow <= stimulus_period[1]:
       n1.accumulate(np.random.normal(0.004, 0.0005))  # TODO accumulate value based on deltaTime (since last accumulate)
     
     if inhibition_period[0] <= timeNow <= inhibition_period[1]:
       i1.accumulate(np.random.normal(0.0035, 0.0005))
     
     n1.update(timeNow)
     #print n1.id, n1.timeCurrent, n1.potential  # [log: potential]
     n2.update(timeNow)
     #print n2.id, n2.timeCurrent, n2.potential  # [log: potential]
     i1.update(timeNow)
     
     subplot(311)  # [graph]
     n1.plot()  # [graph]
     subplot(312)  # [graph]
     n2.plot()  # [graph]
     subplot(313)  # [graph]
     i1.plot()  # [graph]
     pause(delay)  # [graph]
     #sleep(delay)
   
   if holdPlots:
     show()  # [graph]