trials_1 = xlsread('Block1_3');
saveloc = '/Users/tinahocher/Desktop/MATLAB_EEGLAB_DATA /Phase_stop/';
t = [1/fsample:1/fsample:protlength];
beta{1} = sin(IBF*2*pi*[1/fsample:1/fsample:protlength]);
stimuli = zeros(5,protlength.*fsample);
stimuli(3,in) = beta{it}(in).*linspace(0,1, length(in));
save([saveloc 'try_stim_' subject], 'stimuli');
start_trial(i) = ((start_trial(i))*fsample) +(ramptime*fsample);
for j = 1:length(start_trial)
stimuli(5,round(start_trial(j)-5):round(start_trial(j)+5)) = 16;
stop_time(i) = ((stop_time(i))*fsample);
stop_trigger(i) = stop_time(i) + start_trial(i);
else stop_trigger(i) = NaN;
stop_trigger_final = stop_trigger(~isnan(stop_trigger));
stop_phase = x(~isnan(stop_trigger),3);
for j = 1:length(stop_trigger_final)
beta_stop(j) = stimuli(3,round(stop_trigger_final(j)));
loc_phase = stimuli(3,[round(stop_trigger_final(j)-250):round(stop_trigger_final(j)+250)])==1;
loc_phase = find(loc_phase==1);
sample_loc_phase(j) = (round(stop_trigger_final(j)-[250-loc_phase]));
elseif stop_phase(j) == 2
data = stimuli(3,[round(stop_trigger_final(j)-250):round(stop_trigger_final(j)+250)]);
zerocross = find(diff(sign(data)));
for z = 1:length(zerocross)
if data(zerocross(z)+1) < 0
sample_loc_phase(j) = (round(stop_trigger_final(j)-[250-zerocross(z)]));
else sample_loc_phase(j) = (round(stop_trigger_final(j)-[250-zerocross]));
elseif stop_phase(j) == 3
loc_phase = stimuli(3,[round(stop_trigger_final(j)-250):round(stop_trigger_final(j)+250)])==-1;
loc_phase = find(loc_phase==1);
sample_loc_phase(j) = (round(stop_trigger_final(j)-[250-loc_phase]));
elseif stop_phase(j) == 4
data = stimuli(3,[round(stop_trigger_final(j)-250):round(stop_trigger_final(j)+250)]);
zerocross = find(diff(sign(data)));
for z = 1:length(zerocross)
if data(zerocross(z)+1) > 0
sample_loc_phase(j) = (round(stop_trigger_final(j)-[250-zerocross(z)]));
else sample_loc_phase(j) = (round(stop_trigger_final(j)-[250-zerocross]));
else sample_loc_phase(j) = 0;
for j = 1:length(sample_loc_phase)
stimuli(5,round(sample_loc_phase(j)-5):round(sample_loc_phase(j)+5)) = 32;
save([saveloc 'stop_phase_2_4_' subject], 'stimuli');
for j = 1:length(sample_loc_phase)
plot(t(sample_loc_phase(j)-8000:sample_loc_phase(j)+1000), stimuli(3,sample_loc_phase(j)-8000:sample_loc_phase(j)+1000))
plot(t(sample_loc_phase(j)-8000:sample_loc_phase(j)+1000), stimuli(5,sample_loc_phase(j)-8000:sample_loc_phase(j)+1000)/10)