http://w.atwiki.jp/tkyswiki/ tkyswiki @ ウィキ ja 2013-11-28T15:27:33+09:00 1385620053 https://w.atwiki.jp/tkyswiki/pages/34.html URL: -[[iwahori>http://www.cvl.cs.chubu.ac.jp/~iwahori/Ogl/]] 2013-11-28T15:27:33+09:00 1385620053 https://w.atwiki.jp/tkyswiki/pages/33.html octave:122> max(abs(signal_n1)) ans = 7.1833 octave:123> max(abs(signal_n2)) ans = 5.3990 octave:124> max(abs(signal_n3)) ans = 6.0399 octave:125> max(abs(signal_n4)) ans = 8.3160 octave:126> max(abs(signal_n5)) ans = 6.9596 octave:127> max(abs(signal_n21)) ans = 5.1727 octave:128> max(abs(signal_n22)) ans = 3.6067 octave:129> max(abs(signal_n23)) ans = 3.6014 octave:130> max(abs(signal_n24)) ans = 5.7461 2013-07-18T11:52:00+09:00 1374115920 https://w.atwiki.jp/tkyswiki/pages/32.html set(gca,'FontName','Helvetica','Fontsize',12) 2013-07-11T17:53:39+09:00 1373532819 https://w.atwiki.jp/tkyswiki/pages/31.html signal = zeros(2001,1); for i = 1:10 signal = signal + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21); end signal1 = zeros(2001,1); for i = 1:50 signal1 = signal1 + move_data(eog_right2(i)-1000:eog_right2(i)+1000,21); end hold on plot(-1:.001:1,move_data(eog_right2(1)-1000:eog_right2(1)+1000,21),'g') plot(-1:.001:1,signal/10) plot(-1:.001:1,signal1/50,'r') #ref(zemi_kasan1.jpg) 縦軸：μV　横軸：眼球を動かした時間を基準に秒 2013-06-27T19:05:24+09:00 1372327524 https://w.atwiki.jp/tkyswiki/pages/30.html -加算平均が少ない->多い --波形がどうなるか --SN比が何回やると上がるか(SN比の検討) --スペクトルがどうなるか 2013-06-27T19:20:58+09:00 1372328458 https://w.atwiki.jp/tkyswiki/pages/29.html あほだ 2014-05-20T20:20:44+09:00 1400584844 https://w.atwiki.jp/tkyswiki/pages/28.html 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　　ｘ：時間) ---- #region(加算平均) #region(加算平均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; //xlim([1 2001]); //11~20個目の加算平均 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; //21~30個目の加算平均 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; //31~40個目の加算平均 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) yli 2013-10-24T19:45:59+09:00 1382611559 https://w.atwiki.jp/tkyswiki/pages/27.html [[sig1]] = wavread("oc1.wav"); [[sig2]] = wavread("oc2.wav"); sig3 = wavread("oc3.wav"); plot(sig3); xlim([5000 11000]); sig3a=sig3(7001:7512); plot(sig3); xlim([23000 28000]); sig3i=sig3(25001:25512); plot(sig3); xlim([42000 48000]); sig3u=sig3(44001:44512); plot(sig3); xlim([70000 77000]); sig3e=sig3(72501:73012); plot(sig3); xlim([90000 100000]); sig3o=sig3(94001:94512); fft_sig3a = fft(sig3a); log_sig3a = log(abs(fft_sig3a)); ifft_sig3a = ifft(log_sig3a); zero_sig3a = ifft_sig3a; zero_sig3a(12:512-12) = 0; fft2_sig3a = fft(zero_sig3a); plot(8000/512:8000/512:8000,fft2_sig3a,"r"); fft_sig3i = fft(sig3i); log_sig3i = log(abs(fft_sig3i)); ifft_sig3i = ifft(log_sig3i); zero_sig3i = ifft_sig3i; zero_sig3i(12:512-12) = 0; fft2_sig3i = fft(zero_sig3i); plot(8000/512:8000/512:8000,fft2_sig3i,"r"); fft_sig3u = fft(sig3u); log_sig3u = log(abs(fft_sig3u)); ifft_sig3u = ifft(log_sig3u); zero_sig3u = ifft_sig3u; ze 2013-06-06T19:27:15+09:00 1370514435 https://w.atwiki.jp/tkyswiki/pages/26.html [[sig1]] = wavread("oc1.wav"); sig2 = wavread("oc2.wav"); [[sig3]] = wavread("oc3.wav"); plot(sig2) xlim([9000 15000]); sig2a=sig2(10001:10512); plot(sig2) xlim([22000 31000]); sig2i=sig2(25501:26012); plot(sig2) xlim([43000 48000]); sig2u=sig2(45001:45512); plot(sig2) xlim([62000 67000]); sig2e=sig2(64001:64512); plot(sig2) xlim([90000 95000]); sig2o=sig2(93001:93512); fft_sig2a = fft(sig2a); log_sig2a = log(abs(fft_sig2a)); ifft_sig2a = ifft(log_sig2a); zero_sig2a = ifft_sig2a; zero_sig2a(12:512-12) = 0; fft2_sig2a = fft(zero_sig2a); plot(8000/512:8000/512:8000,fft2_sig2a,"r"); fft_sig2i = fft(sig2i); log_sig2i = log(abs(fft_sig2i)); ifft_sig2i = ifft(log_sig2i); zero_sig2i = ifft_sig2i; zero_sig2i(12:512-12) = 0; fft2_sig2i = fft(zero_sig2i); plot(8000/512:8000/512:8000,fft2_sig2i,"r"); fft_sig2u = fft(sig2u); log_sig2u = log10(abs(fft_sig2u)); ifft_sig2u = ifft(log_sig2u); zero_sig2u = ifft_sig2u; zero_sig2 2013-06-06T19:21:22+09:00 1370514082 https://w.atwiki.jp/tkyswiki/pages/25.html sig1 = wavread("oc1.wav"); [[sig2]] = wavread("oc2.wav"); [[sig3]] = wavread("oc3.wav"); plot(sig1) xlim([5000 10000]); sig1a=sig1(7001:7512); plot(sig1) xlim([20000 25000]); sig1i=sig1(21501:22012); plot(sig1) xlim([35000 40000]); sig1u=sig1(36001:36512); plot(sig1) xlim([47000 53000]); sig1e=sig1(48501:49012); plot(sig1) xlim([60000 65000]); sig1o=sig1(62001:62512); fft_sig1a = fft(sig1a); log_sig1a = log10(abs(fft_sig1a)); ifft_sig1a = ifft(log_sig1a); zero_sig1a = ifft_sig1a; zero_sig1a(12:512-12) = 0; fft2_sig1a = fft(zero_sig1a); plot(8000/512:8000/512:8000,fft2_sig1a,"r"); fft_sig1i = fft(sig1i); log_sig1i = log10(abs(fft_sig1i)); ifft_sig1i = ifft(log_sig1i); zero_sig1i = ifft_sig1i; zero_sig1i(12:512-12) = 0; fft2_sig1i = fft(zero_sig1i); plot(8000/512:8000/512:8000,fft2_sig1i,"r"); fft_sig1u = fft(sig1u); log_sig1u = log10(abs(fft_sig1u)); ifft_sig1u = ifft(log_sig1u); zero_sig1u = ifft_sig1u; zero_sig1u(12:512-1 2013-06-06T19:10:53+09:00 1370513453