位相差の求め方

現在１９ｃｈ×100のデータがあります。

これの位相差を求めたいです。

1chに注目した時の１ｃｈと１ｃｈ、１ｃｈと２ｃｈ、１ｃｈと３ｃｈ。。。。。１ｃｈと１９ｃｈ

２ｃｈに注目した時２ｃｈと１ｃｈ、２ｃｈと２ｃｈ。。。。。２ｃｈと１９ｃｈ

これを１９ｃｈに注目した時までの位相を求めたいです

１ｃｈと２ｃｈだけの位相差であれば１ｃｈー２ｃｈをするだけなのですが、

これを全部書くことはかなり時間かかるのでforloopなどで簡単に出来るやり方がありましたら

スクリプトを送っていただきたいです。

よろしくお願いいたします。

下記スクリプトは全体ではないですが５ｃｈー10ｃｈの位相差のスクリプトになります

上の方のでα波取り出しとヒルベルト変換していますが、それぞれ全体の位相差部分のスクリプトだけで大丈夫です。

data=record(5,:);

samplerate=250;

length=231250;

t=(0:length-1)/samplerate;

srate=250;

locutoff=8;

hicutoff=13;

smoothdata1 = eegfilt(data,srate,locutoff,hicutoff);

y1=hilbert(smoothdata1);%ヒルベルト変換

p=angle(y);

hold on

data1=record(10,:);

length=231250;

p2=angle(y1); %位相

dif1=p-p2;%位相差(5-10ch)

Answer 1

[Atsushi Ueno]

for idx1 = 1:19
    for idx2 = 1:19
        % 関数でangleを取ってからi番目のchに注目した時の位相差を処理してください
        % １９×１００が１９個出来ると思います。
        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