for idx1 = 1:19
for idx2 = 1:19
% 関数でangleを取ってからi番目のchに注目した時の位相差を処理してください
% 19×100が19個出来ると思います。
angle1 = get_angle(record(idx1,:));
angle2 = get_angle(record(idx2,:));
dif1{idx1,idx2} = angle1 - angle2;
end
end
function p = get_angle(data)
srate = 250;
locutoff = 8;
hicutoff = 13;
smoothdata1 = eegfilt(data,srate,locutoff,hicutoff);
y = hilbert(smoothdata1);%ヒルベルト変換
p = angle(y); %位相
end