From 37f038e81c7620aac091f4d1059ab2dd7f97777e Mon Sep 17 00:00:00 2001 From: Schellenberger Date: Mon, 10 Mar 2014 16:00:05 +0100 Subject: [PATCH] Updated the matlab script --- Plots.m | 18 +++--------------- 1 file changed, 3 insertions(+), 15 deletions(-) diff --git a/Plots.m b/Plots.m index 1134f8f..49df712 100644 --- a/Plots.m +++ b/Plots.m @@ -15,26 +15,14 @@ function Plots(T, onset) T = 30; % duration of the simulation end -[Vt, Vr] = Thalamus(T, Con, var_stim); +[Vt] = Thalamus(T, Con, var_stim); L = max(size(Vt)); timeaxis = linspace(0,T,L); figure(1) subplot(211), plot(timeaxis,Vt) -title('Exitatory membrane voltage'), xlabel('time in s'), ylabel('Vt in mV') +title('Thalamic relay membrane voltage'), xlabel('time in s'), ylabel('Vt in mV') subplot(212), plot(timeaxis,Vr) -title('Inhibitory membrane voltage'), xlabel('time in s'), ylabel('Vr in mV') - -% -fs = L/T; -NFFT = 2^nextpow2(L); - -[Pxx,f] = pwelch(Vt-mean(Vt), [], [], [], fs,'onesided'); -n = find(f<=60, 1, 'last' ); - -figure(3) -plot(f(1:n),Pxx(1:n)) -title('Powerspectrum of Steyn-Ross model with pwelch'), xlabel('frequency in Hz'), ylabel('Power') - +title('Thalamic reticular membrane voltage'), xlabel('time in s'), ylabel('Vr in mV') %save('Thalamus.mat','Vt','Vr')