From f20dabf6c712bdaaf4f4f2b58bb482d1b17f8814 Mon Sep 17 00:00:00 2001 From: Schellenberger Date: Wed, 14 May 2014 17:33:06 +0200 Subject: [PATCH] Changed parameters so that the model produces spindle oscillations --- Plots.m | 18 ++++++++++++------ Thalamic_Column.h | 29 +++++++++++++---------------- 2 files changed, 25 insertions(+), 22 deletions(-) diff --git a/Plots.m b/Plots.m index 8c2457c..5ff4726 100644 --- a/Plots.m +++ b/Plots.m @@ -4,24 +4,30 @@ function Plots(T) if nargin == 0 - Con = [2; % N_tr - 5.5; % N_rt - 5]; % N_rr + Con = [10; % N_tr + 10; % N_rt + 40]; % N_rr var_stim = [ 0; % strength of the stimulus in Hz (spikes per second) 0; % time between stimuli in s 0; % time until first stimuli in s 0]; % duration of the stimulus in ms - T = 60; % duration of the simulation + T = 30; % duration of the simulation end [Vt, Vr] = Thalamus(T, Con, var_stim); L = max(size(Vt)); timeaxis = linspace(0,T,L); +fs = L/T; +[Pxx,f] = pwelch(Vt-mean(Vt), [], [], [], fs,'onesided'); +n = find(f<=60, 1, 'last' ); + figure(1) -subplot(211), plot(timeaxis,Vt) +subplot(311), plot(timeaxis,Vt) title('Thalamic relay membrane voltage'), xlabel('time in s'), ylabel('Vt in mV') -subplot(212), plot(timeaxis,Vr) +subplot(312), plot(timeaxis,Vr) 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') diff --git a/Thalamic_Column.h b/Thalamic_Column.h index 1b8575d..bd8e0d5 100644 --- a/Thalamic_Column.h +++ b/Thalamic_Column.h @@ -139,12 +139,12 @@ class Thalamic_Column { const double Qr_max = 400.E-3; /* Sigmoid threshold in mV */ - const double theta_t = -45; - const double theta_r = -45; + const double theta_t = -58.6; + const double theta_r = -58.6; /* Sigmoid gain in mV */ - const double sigma_t = 3; - const double sigma_r = 3; + const double sigma_t = 9; + const double sigma_r = 9; /* Scaling parameter for sigmoidal mapping (dimensionless) */ const double C1 = (3.14159265/sqrt(3)); @@ -153,24 +153,21 @@ class Thalamic_Column { const double gamma_e = 70E-3; const double gamma_i = 58.6E-3; - /* axonal flux time constant */ - const double nu = 120E-3; - /* Conductivities in mS/cm^-2 */ /* Leak current */ - const double g_L_t = 0.6; - const double g_L_r = 1.5; + const double g_L_t = 1; + const double g_L_r = 1; /* Potassium leak current */ const double g_LK_t = 0.02; const double g_LK_r = 0.02; /* T current */ - const double g_T_t = 2.2; + 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.04; /* Reversal potentials in mV */ /* Synaptic */ @@ -178,8 +175,8 @@ class Thalamic_Column { const double E_GABA = -70; /* Leak */ - const double E_L_t = -64; - const double E_L_r = -64; + const double E_L_t = -70; + const double E_L_r = -70; /* Potassium */ const double E_K = -100; @@ -188,11 +185,11 @@ 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 */ - const double tau_Ca = 5; /* calcium time constant in ms */ + const double tau_Ca = 10; /* calcium time constant in ms */ const double Ca_0 = 2E-4; /* resting concentration */ /* I_h activation parameters */ @@ -205,7 +202,7 @@ class Thalamic_Column { /* Noise parameters in ms^-1 */ const double mphi = 0E-3; - const double dphi = 2E-3;; + const double dphi = 10E-3;; double input = 0.0; /* Connectivities (dimensionless) */