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 <schellenberger@inb.uni-luebeck.de>

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')

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;
 };

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			*/