123456789101112131415161718192021 |
- import matplotlib.pyplot as plt
- import numpy
- from brian2.units import *
- from scripts.spatial_network.head_direction_index_over_noise_scale import create_head_direction_input
- tunings = numpy.linspace(-numpy.pi, numpy.pi, 100, endpoint=False)
- baseline = 0.05 * nA
- amplitude = 0.6 * nA
- sharpness = 2.5
- direction = 0
- input_to_excitatory_population = create_head_direction_input(baseline, tunings,
- sharpness,
- amplitude, direction)
- plt.plot(tunings / numpy.pi * 180, input_to_excitatory_population / nA)
- plt.xlabel("$\\theta-\\theta_p$")
- plt.ylabel("I (nA)")
- plt.ylim(0, numpy.max(input_to_excitatory_population / nA) * 1.05)
- plt.savefig('../../figures/evolve_spatial_layout/broad_input.png', dpi=300)
|