カットオフ周波数を決定するための式

for Fc = 1:50

Wn = ((Fc)/(Fs/2));

[b,a] = butter(order,Wn,'low');

raw = filtfilt(b,a,data);

rev_raw = flipud(raw);

filrev_raw = filtfilt(b,a,rev_raw);

alter_filrev_raw(Fc,:) = flipud(filrev_raw);

end

ノイズを除去するためのバターワースフィルタのカットオフ周波数を決定するためにこのコードを使用しているのですが、「左辺のインデックスが右辺とサイズが適合しないため、代入は実行できません」というエラーが出てしまいます。

この場合の原因はなんでしょうか？詳しい方に教えていただけると助かります。

[takemoto]

エラーメッセージを見ると、

alter_filrev_raw(Fc,:) = flipud(filrev_raw);

の代入がうまく処理できていない様です。filrev_rawが行ベクトルとなっていない可能性がありますので、まずはfilrev_rawのサイズを確認されてはいかがでしょうか。

[光貴 川島]

ありがとうございます！filver_rawは170行×116列の行列でした。行ベクトルに変換して試してみます！