win10のMATLAB R2017aで開くと下記エラーの為開けず、エラーの出る状況を再現出来ませんでした。
エラー: load
MAT ファイル csp_dataset_test.mat を読み込むことができません。ファイルが破損している可能性があります。
FLT (96735 × 25 double) の「位置1のインデックスが配列範囲を超えています (96735 を超えてはならない)」ですから、下記の計算結果が96735 を超える時のmrkがどのような値になっているか調べれば原因を特定できると思います。データの欠損などによりNaNが紛れ込む可能性も考えましたがfind関数やnnz関数ではその様な事にならないようです。
% repmat(find(mrk == k),length(wnd),1) + repmat(wnd',1,nnz(mrk == k))
- 上記の値が96735 を超えた事実がある事は間違いない
- mrk==kの合致数は最大24個だがおそらく常に1個
- なのでwnd'(125~875の1×751 double)にmrk==kがヒットしたインデックス(1~24)が加算される(最大値は899)
- kの範囲を1:2より増やしても上記は同様、どう考えても96735どころか899を超えない
下記の様に計算が通る適当な値を入れると、エラーはでません。
flt = @(f)(f>7&f<30).*(1-cos((f-(7+30)/2)/(7-30)*pi*4));
EEG = single(ones(96735,25));
mrk = sparse(1,1:24,((-1:3/24:1.9)+3)/2);
FLT = real(ifft(fft(EEG).*repmat(flt(250*(0:96735-1)/96735)',1,25)));
wnd = round(250*0.5):round(250*3.5)
for k=1:2
EPO{k} = FLT(repmat(find(mrk == k),length(wnd),1) + repmat(wnd',1,nnz(mrk == k)),:)
end
