load('data_fig_1.db','-mat');


subplot(2,1,1)

%% theta frequency
frwin=3;
for a=1:3
    rx=size(data.unwrapped_phase_right{a},1);
    freq_right{a}=([data.unwrapped_phase_right{a}(:,1+frwin:end) nan(rx,frwin)]-[nan(rx,frwin) data.unwrapped_phase_right{a}(:,1:end-frwin)])/2/pi/0.02/frwin/2;
    rx=size(data.unwrapped_phase_left{a},1);
    freq_left{a}=([data.unwrapped_phase_left{a}(:,1+frwin:end) nan(rx,frwin)]-[nan(rx,frwin) data.unwrapped_phase_left{a}(:,1:end-frwin)])/2/pi/0.02/frwin/2;    
end

%% acceleration
win=10;    %gfr=fr;


    rx=size(x,1);
    acc_right=-([nan(rx,win) x(:,1:end-win)]+[x(:,1+win:end) nan(rx,win)]-2*x)/0.02/0.02/win/win;


axlim=0; fct=0; th=0; 
th=ths(r); axlim=[-120  220]; %fct=0.0127;
new_figfig
% set(AX(2),'XTickLabel',[]);
% set(AX(1),'XTickLabel',[]);

set(AX(2),'XTick',[]);
% L=legend('acceleration','CA1','MEC','Location', 'NorthEast');
% set(L,'Color','none')
% set(L,'box','off')
% set(L,'FontSize',10,'FontName','Arial')
subplot(2,1,2)
% win=5;
if r==3
    gfr=data.unwrapped_phase_left; gacc=lacc;
    rx=size(lx,1);
    gacc=([nan(rx,win) lx(:,1:end-win)]+[lx(:,1+win:end) nan(rx,win)]-2*lx)/0.02/0.02/win/win;    
else
    gfr=fr; gacc=acc; 
    rx=size(x,1);
    gacc=-([nan(rx,win) x(:,1:end-win)]+[x(:,1+win:end) nan(rx,win)]-2*x)/0.02/0.02/win/win;    
end
% axlim=[-160 260];
new_figfig
set(AX(2),'XTick',[]);