ウェーブレットコヒー​レンスの一部データを​抽出する方法

2種類の周波数データx,yを下記関数にてウェーブレットコヒーレンス処理を実施した後、一部の周波数帯域（15-30Hz）のデータ抽出を試みています。

しかし、下記手順にて実施すると、下線部にてエラーが出てしまいます。

手順につき、ご教示頂ければ幸いです。

fs=1000; %サンプリング周波数

[wcoh,wcs,f]= wcoherence(x,y,fs)

t=(1:length(wcoh)); %時間

indx_f = (f>=15)&(f<=30); % 周波数15-30Hz

frq = f(indx_f);

wcoh1 = wcoh(frq, :);

surface(t,frq,wcho1); %図

Answer 1

[Naoya]

wcoh1 は 行列ですので、 実際の周波数の値を入れるのではなく

インデックスを表す整数、または論理値 を入れることになります。

以下部分を修正すればよいと思います。

wcoh1 = wcoh(frq, :);

↓

wcoh1 = wcoh(indx_f, :);

[Kenta Otsuka]

ありがとうございました。