123456789101112131415161718192021222324252627282930313233343536373839404142434445464748 |
- 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',[]);
|