From 4f9ea29db6fcc15e86457a8f31546bb7a408e814 Mon Sep 17 00:00:00 2001 From: Schellenberger Date: Wed, 13 Aug 2014 11:29:59 +0200 Subject: [PATCH] Changed Parameters. Model now shows nice spindle oscillations with a refactory time of slightly above 4s. Can be entrained by stimuli between 5s and 6s. Looses entrainment around 7s. Signed-off-by: Schellenberger --- Plots.m | 18 +++++++++--------- Thalamic_Column.h | 16 ++++++++-------- Thalamus.cpp | 2 +- 3 files changed, 18 insertions(+), 18 deletions(-) diff --git a/Plots.m b/Plots.m index 3540e35..af7c3b9 100644 --- a/Plots.m +++ b/Plots.m @@ -4,10 +4,10 @@ function Plots(T) if nargin == 0 - Con = [ 2; % sigma - 5; % N_tr - 5; % N_rt - 50]; % N_rr + Con = [ 4; % sigma + 6; % N_tr + 5; % N_rt + 100]; % N_rr % stimulation parameters @@ -17,10 +17,10 @@ function Plots(T) % 2 == phase dependend up state % 3 == phase dependend down state - var_stim = [ 1 % mode of stimulation - 30; % strength of the stimulus in Hz (spikes per second) - 50; % duration of the stimulus in ms - 5; % time between stimuli in s + var_stim = [ 0; % mode of stimulation + 200; % strength of the stimulus in Hz (spikes per second) + 100; % duration of the stimulus in ms + 6; % time between stimuli in s 1]; % time until stimuli after min in ms T = 30; % duration of the simulation end @@ -40,4 +40,4 @@ function Plots(T) title('Thalamic reticular membrane voltage'), xlabel('time in s'), ylabel('Vr in mV') subplot(313), plot(f(1:n),Pxx(1:n)) title('Powerspectrum of Steyn-Ross model with pwelch'), xlabel('frequency in Hz'), ylabel('Power') -%save('Thalamus.mat','Vt','Vr') +save('Thalamus.mat','Vt','Vr') diff --git a/Thalamic_Column.h b/Thalamic_Column.h index 4a71e31..f42bfc5 100644 --- a/Thalamic_Column.h +++ b/Thalamic_Column.h @@ -152,7 +152,7 @@ class Thalamic_Column { /* PSP rise time in ms^-1 */ const double gamma_e = 70E-3; - const double gamma_i = 58.6E-3; + const double gamma_i = 100E-3; /* Conductivities in mS/cm^-2 */ /* Leak current */ @@ -161,14 +161,14 @@ class Thalamic_Column { /* Potassium leak current */ const double g_LK_t = 0.02; - const double g_LK_r = 0.012; + const double g_LK_r = 0.02; /* T current */ const double g_T_t = 3; const double g_T_r = 2; /* h current */ - const double g_h = 0.07; + const double g_h = 0.08; /* Reversal potentials in mV */ /* Synaptic */ @@ -186,7 +186,7 @@ class Thalamic_Column { const double E_Ca = 120; /* I_h current */ - const double E_h = -40; + const double E_h = -43; /* Calcium parameters */ const double alpha_Ca = -50E-6; /* influx per spike in nmol */ @@ -203,13 +203,13 @@ class Thalamic_Column { /* Noise parameters in ms^-1 */ const double mphi = 0E-3; - const double dphi = 10E-3;; + const double dphi = 10E-3; double input = 0.0; /* Connectivities (dimensionless) */ - const double N_tr = 10; - const double N_rt = 10; - const double N_rr = 40; + const double N_tr = 5; + const double N_rt = 5; + const double N_rr = 50; friend class Stim; }; diff --git a/Thalamus.cpp b/Thalamus.cpp index a4bc480..0d9d672 100644 --- a/Thalamus.cpp +++ b/Thalamus.cpp @@ -38,7 +38,7 @@ /****************************************************************************************************/ /* Fixed simulation settings */ /****************************************************************************************************/ -extern const int onset = 10; /* 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 */