6/13 - tkyswiki @ ウィキ

load eog_left2.mat
load eog_right2.mat
load move_data.mat

加算平均したときの誤差

plot(-1:.001:1,move_data(eog_right2(1)-1000:eog_right2(1)+1000,21))
1回目の波を表示
(ｙ：マイクロV　　ｘ：時間)

---

加算平均

加算平均10個

1~10個目の加算平均

signal = zeros(2001,1);
for i = 1:10
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
figure
plot(-1:.001:1,signal/10)
ylim([-50 50]);
signal1=signal;

signal = zeros(2001,1);
for i = 11:20
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/10)
ylim([-50 50]);
signal2=signal;

signal = zeros(2001,1);
for i = 21:30
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/10)
ylim([-50 50]);
signal3=signal;

signal = zeros(2001,1);
for i = 31:40
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/10)
ylim([-50 50]);
signal4=signal;

signal = zeros(2001,1);
for i = 41:50
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/10)
ylim([-50 50]);
signal5=signal;

plot

plot(-1:.001:1,signal1/10,'b')
hold on
plot(-1:.001:1,signal2/10,'r')
plot(-1:.001:1,signal3/10,'g')
plot(-1:.001:1,signal4/10,'m')
plot(-1:.001:1,signal5/10,'c')

---

20個

1~20個目の加算平均

signal = zeros(2001,1);
for i = 1:20
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/20)
ylim([-50 50]);
signal21=signal;

11~30個目の加算平均

signal = zeros(2001,1);
for i = 11:30
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/20)
ylim([-50 50]);
signal22=signal;

21~40個目の加算平均

signal = zeros(2001,1);
for i = 21:40
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/20)
ylim([-50 50]);
signal23=signal;

31~50個目の加算平均

signal = zeros(2001,1);
for i = 31:50
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/20)
ylim([-50 50]);
signal24=signal;

---

30個

1~30個目の加算平均

signal = zeros(2001,1);
for i = 1:30
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/30)
ylim([-50 50]);
signal31=signal;

21~50個目の加算平均

signal = zeros(2001,1);
for i = 11:40
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/30)
ylim([-50 50]);
signal32=signal;

21~50個目の加算平均

signal = zeros(2001,1);
for i = 21:50
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/30)
ylim([-50 50]);
signal33=signal;

---

フーリエ変換

fft_signal = fft(signal);
log_signal = log10(abs(fft_signal));
ifft_signal = ifft(log_signal);
figure
plot(real(ifft_signal))
ylim([-0.5 3]);
figure
plot(f/t:f/t:f/2,abs(fft_signal(2:1001)));
ylim([-1500 100000]);

---

octave:42> zero_y1 = fft_y1; octave:43> length（fft_y1 ） ans = 378874 octave:44> length（fft_y1）; octave:45> length（fft_y1）; octave:46> length（fft_y1） ans = 378874 octave:47> length（fft_y1）/388747 ans = 0.97460 octave:48> length（fft_y1）/378874 ans = 1 octave:49> length（fft_y1）*10/378874 ans = 10 octave:50> length（fft_y1）*10/51200 ans = 73.999 octave:51> zero_y1（1:74）=zeros; octave:52> zero_y1（length（zero_y1）-74:length（zero_y1））=zeros; octave:53> plot（real（zero_y1））; octave:54> plot（real（zero_y1））; octave:55> ifft_y1 = ifft（zero_y1）; octave:56> plot（real（ifft_y1））; octave:57> zero_y1 = fft_y1; octave:58> zero_y1（1:74*5）=zeros; octave:59> zero_y1（length（zero_y1）-（74*5）:length（zero_y1））=zeros; octave:60> ifft_y1 = ifft（zero_y1）; octave:61> plot（real（ifft_y1））;

signal = zeros(2001,1);
for i = 1:50
signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21);
end
plot(-1:.001:1,signal/50)
hold off
plot(-1:.001:1,signal/50)
signal_S = signal;
signal_N101 = signal-signal1;
signal_SN = signal_S/signal_N101;
signal_S = signal/50;
signal_N101 = signal_S - signal1/10;
signal_SN = signal_S/signal_N101;
plot(signal_SN )

signal2 = zeros(2001,1);
for i = 1:10
signal2 = signal2 + move_data(eog_right2(i)-1000:eog_right2(i)+1000,8);
end
plot(-1:.001:1,signal2/10)
ylim([-30 30])

plot(-1:.001:1,signal3/10)
signal4 = zeros(2001,1);
for i = 1:10
signal4 = signal4 + move_data(eog_left2(i)-1000:eog_left2(i)+1000,8);
end
figure
plot(-1:.001:1,signal4/10)
ylim([-40 40])