|
@@ -1,48 +0,0 @@
|
|
|
-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',[]);
|
|
|
-
|