代碼
%%
a_32_384_sum=zeros(32,384);
%創(chuàng)建0矩陣用于相加
for i=1:1:40
a_32_384_sum=a_32_384_sum+EEG.data(:,:,i);
%訪問a的40頁數(shù)據(jù)加起來
end
a_ave=a_32_384_sum/40;
%除以40頁求平均套媚；這個時候就把數(shù)據(jù)轉(zhuǎn)成 32*384 的二維平均後的數(shù)據(jù)了
%%
%如果3,4,5列代表前額窃肠，把第3蛔屹、4、5行數(shù)據(jù)加起來
frontal=(a_ave(3,:)+a_ave(4,:)+a_ave(5,:))/3;
%由于時間是[-1,2]青团，采樣率是128 Hz，所以
time=-1:1/128:2;
%注意這個時候 time 的 size 是 1*385咖楣，所以我們需要在frontal數(shù)據(jù)后面添一個0
frontal=[frontal 0];
%這個時候 time-frontal 的數(shù)據(jù)維度就一致了
%%
%如果21,22,23列代表頂區(qū)督笆，把第21,22,23行數(shù)據(jù)加起來
parietal=(a_ave(21,:)+a_ave(22,:)+a_ave(23,:))/3;
parietal=[parietal 0];
%%
%開始畫圖
plot(time,frontal,'b-',time,parietal,'g-');
title('前額和頂區(qū)在condition1下的比較');
%添上Xlabel和Ylabel
xlabel('t = -1s to 2s');
ylabel('values of ERP');
%設定x軸步長和步長的label
set(gca,'xticklabel',-1000:500:2000);
%設定坐標軸字體和字體大小
set(gca,'fontname','楷體');
set(gca,'FontSize',25);
%但是我們還需要標注0點，再畫一根垂直于x軸的直線
%畫一條垂直x軸的線诱贿，前面是2個點的橫坐標娃肿，后面是縱坐標
line([0,0],[-20,40],'color','m');
%添上圖例,只要對應寫上先后畫出的函數(shù)就好了；加上title
legend('前額ERP','頂區(qū)ERP','時間0點');
輸出
untitled.jpg
