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

4 次查看(过去 30 天)
光貴 川島
光貴 川島 2021-9-1
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
ノイズを除去するためのバターワースフィルタのカットオフ周波数を決定するためにこのコードを使用しているのですが、「左辺のインデックスが右辺とサイズが適合しないため、代入は実行できません」というエラーが出てしまいます。
この場合の原因はなんでしょうか?詳しい方に教えていただけると助かります。
  2 个评论
takemoto
takemoto 2021-9-2
エラーメッセージを見ると、
alter_filrev_raw(Fc,:) = flipud(filrev_raw);
の代入がうまく処理できていない様です。filrev_rawが行ベクトルとなっていない可能性がありますので、まずはfilrev_rawのサイズを確認されてはいかがでしょうか。
光貴 川島
光貴 川島 2021-9-3
ありがとうございます!filver_rawは170行×116列の行列でした。行ベクトルに変換して試してみます!

请先登录,再进行评论。

回答(0 个)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!