06 > 05

EEG=脳電位エレクトロ～グラム
eog=眼電位(Electro OculoGram)エレクトロオキュログラム

配布されたファイルをロード

load move_eye.mat

サンプリング周波数1000

eog_right

• 中央から右に動いた点

eog_left

• 中央から右に動いた点

---

プロットしてみよう
plot(EEG_O1(eog_left2(1)-1000:eog_left2(1)+1000))

       O(オー)1です注意

• EEG_O1
• 左後頭葉
• EEG_O2
• 右後頭葉

からサッカード開始点±1000点のデータを使って同期加算(加算平均)

同期点：サッカード開始点(eog_left2,eog_right2)が格納されている

課題配分
left担当

• 陳～原

right担当

• 福永～磯谷

加算数：2,5,10,50

• O1,O2の比較(加算数50)
• 加算数の比較(加算数2,5,10,50)

S/N比もできれば？やる
課題説明ここまで
%-----------------------------------------------------
%初期化
x2=EEG_O1(eog_left2(1)-1000:eog_left2(1)+1000);
x2(1:end)=0;
x50=x2;
x10=x2;
x5 =x2;
%
%加算数2回
j=2;
for i=1:j
x2+=EEG_O1(eog_left2(i)-1000:eog_left2(i)+1000);
end
O1_02=x2/j;
%加算数5回
j=5;
for i=1:j
x5+=EEG_O1(eog_left2(i)-1000:eog_left2(i)+1000);
end
O1_05=x5/j;
%加算数10回
j=10;
for i=1:j
x10+=EEG_O1(eog_left2(i)-1000:eog_left2(i)+1000);
end
O1_10=x10/j;
%加算数50回
j=50;
for i=1:j
x50+=EEG_O1(eog_left2(i)-1000:eog_left2(i)+1000);
end
O1_50=x50/j;
%プロット
figure
hold on
plot(-1:.001:1,O1_02,'b')
plot(-1:.001:1,O1_05,'r')
plot(-1:.001:1,O1_10,'m')
plot(-1:.001:1,O1_50,'g')
hold off
title("O1_left")
legend("2times","5times","10times","50times")
xlabel("time(s)")
ylabel("voltage(uV)")
%-------------------------------------------------------
%初期化
x2=EEG_O1(eog_right2(1)-1000:eog_right2(1)+1000);
x2(1:end)=0;
x50=x2;
x10=x2;
x5 =x2;
%
%加算数2回
j=2;
for i=1:j
x2+=EEG_O1(eog_right2(i)-1000:eog_right2(i)+1000);
end
O1_02=x2/j;
%加算数5回
j=5;
for i=1:j
x5+=EEG_O1(eog_right2(i)-1000:eog_right2(i)+1000);
end
O1_05=x5/j;
%加算数10回
j=10;
for i=1:j
x10+=EEG_O1(eog_right2(i)-1000:eog_right2(i)+1000);
end
O1_10=x10/j;
%加算数50回
j=50;
for i=1:j
x50+=EEG_O1(eog_right2(i)-1000:eog_right2(i)+1000);
end
O1_50=x50/j;
%プロット
figure
hold on
plot(-1:.001:1,O1_02,'b')
plot(-1:.001:1,O1_05,'r')
plot(-1:.001:1,O1_10,'m')
plot(-1:.001:1,O1_50,'g')
hold off
title("O1_right")
legend("2times","5times","10times","50times")
xlabel("time(s)")
ylabel("voltage(uV)")
%-------------------------------------------------------
%初期化
x2=EEG_O2(eog_left2(1)-1000:eog_left2(1)+1000);
x2(1:end)=0;
x50=x2;
x10=x2;
x5 =x2;
%
%加算数2回
j=2;
for i=1:j
x2+=EEG_O2(eog_left2(i)-1000:eog_left2(i)+1000);
end
O2_02=x2/j;
%加算数5回
j=5;
for i=1:j
x5+=EEG_O2(eog_left2(i)-1000:eog_left2(i)+1000);
end
O2_05=x5/j;
%加算数10回
j=10;
for i=1:j
x10+=EEG_O2(eog_left2(i)-1000:eog_left2(i)+1000);
end
O2_10=x10/j;
%加算数50回
j=50;
for i=1:j
x50+=EEG_O2(eog_left2(i)-1000:eog_left2(i)+1000);
end
O2_50=x50/j;
%プロット
figure
hold on
plot(-1:.001:1,O2_02,'b')
plot(-1:.001:1,O2_05,'r')
plot(-1:.001:1,O2_10,'m')
plot(-1:.001:1,O2_50,'g')
hold off
title("O2_left")
legend("2times","5times","10times","50times")
xlabel("time(s)")
ylabel("voltage(uV)")
%-------------------------------------------------------
初期化
x2=EEG_O2(eog_right2(1)-1000:eog_right2(1)+1000);
x2(1:end)=0;
x50=x2;
x10=x2;
x5 =x2;
%
%加算数2回
j=2;
for i=1:j
x2+=EEG_O2(eog_right2(i)-1000:eog_right2(i)+1000);
end
O2_02=x2/j;
%加算数5回
j=5;
for i=1:j
x5+=EEG_O2(eog_right2(i)-1000:eog_right2(i)+1000);
end
O2_05=x5/j;
%加算数10回
j=10;
for i=1:j
x10+=EEG_O2(eog_right2(i)-1000:eog_right2(i)+1000);
end
O2_10=x10/j;
%加算数50回
j=50;
for i=1:j
x50+=EEG_O2(eog_right2(i)-1000:eog_right2(i)+1000);
end
O2_50=x50/j;
%プロット
figure
hold on
plot(-1:.001:1,O2_02,'b')
plot(-1:.001:1,O2_05,'r')
plot(-1:.001:1,O2_10,'m')
plot(-1:.001:1,O2_50,'g')
hold off
title("O2_right")
legend("2times","5times","10times","50times")
xlabel("time(s)")
ylabel("voltage(uV)")
%-------------------------------------------------------

%plot(-1:.001:1,x)