diff --git a/Plots.m b/Plots.m index 5ff4726..2a17d78 100644 --- a/Plots.m +++ b/Plots.m @@ -4,9 +4,11 @@ function Plots(T) if nargin == 0 - Con = [10; % N_tr - 10; % N_rt - 40]; % N_rr + Con = [ 2; % sigma + 5; % N_tr + 5; % N_rt + 50]; % N_rr + var_stim = [ 0; % strength of the stimulus in Hz (spikes per second) 0; % time between stimuli in s diff --git a/Thalamic_Column.h b/Thalamic_Column.h index bd8e0d5..7bd4b40 100644 --- a/Thalamic_Column.h +++ b/Thalamic_Column.h @@ -54,7 +54,8 @@ class Thalamic_Column { /* Constructor for simulation */ Thalamic_Column(double* Con) - : N_tr (Con[0]), N_rt (Con[1]), N_rr (Con[2]) + : sigma_t (Con[0]), sigma_r (Con[0]), + N_tr (Con[1]), N_rt (Con[2]), N_rr (Con[3]) {set_RNG();} /* Set strength of input */ @@ -131,8 +132,8 @@ class Thalamic_Column { /* Declaration and Initialization of parameters */ /* Membrane time in ms */ - const double tau_t = 30; - const double tau_r = 30; + const double tau_t = 20; + const double tau_r = 20; /* Maximum firing rate in ms^-1 */ const double Qt_max = 400.E-3; @@ -143,8 +144,8 @@ class Thalamic_Column { const double theta_r = -58.6; /* Sigmoid gain in mV */ - const double sigma_t = 9; - const double sigma_r = 9; + const double sigma_t = 3; + const double sigma_r = 3; /* Scaling parameter for sigmoidal mapping (dimensionless) */ const double C1 = (3.14159265/sqrt(3)); @@ -160,14 +161,14 @@ class Thalamic_Column { /* Potassium leak current */ const double g_LK_t = 0.02; - const double g_LK_r = 0.02; + const double g_LK_r = 0.012; /* T current */ const double g_T_t = 3; const double g_T_r = 2; /* h current */ - const double g_h = 0.04; + const double g_h = 0.07; /* Reversal potentials in mV */ /* Synaptic */ @@ -206,9 +207,9 @@ class Thalamic_Column { double input = 0.0; /* Connectivities (dimensionless) */ - const double N_tr = 2; - const double N_rt = 5.5; - const double N_rr = 5; + const double N_tr = 10; + const double N_rt = 10; + const double N_rr = 40; friend class Stim; }; diff --git a/Thalamus.cpp b/Thalamus.cpp index f1ed9d0..0d9d672 100644 --- a/Thalamus.cpp +++ b/Thalamus.cpp @@ -38,7 +38,7 @@ /****************************************************************************************************/ /* Fixed simulation settings */ /****************************************************************************************************/ -extern const int onset = 5; /* time until data is stored in s */ +extern const int onset = 15; /* time until data is stored in s */ extern const int res = 1E4; /* number of iteration steps per s */ extern const int red = res/100; /* number of iterations that is saved */ extern const double dt = 1E3/res; /* duration of a timestep in ms */