冷杉フィルターによって生じる遅延の補正

信号をフィルター処理すると遅延が生じます。つまり入力に対し,出力信号では時間のシフトが行われます。この例では,この影響を抑える方法を説明します。

有限インパルス応答フィルターでは,しばしばすべての周波数成分が同じ量だけ遅延します。このため,信号を時間シフトすれば,遅延を簡単に修正できます。

1秒間に500 Hzでサンプルされた心電図の読み取りを例に取ります。ランダムノイズを追加します。再現性が必要な場合は,乱数発生器をリセットします。

Fs = 500;N = 500;rng默认的xn = ecg(N)+0.25*randn([1 N]);tn = (0: n - 1) / Fs;

75赫兹より高い周波数を阻止するフィルターを使用して,ノイズの一部を除去します。designfiltを使用して次数70のフィルターを設計します。

nfilt = 70;置= 75;d = designfilt (“lowpassfir”，“FilterOrder”nfilt,.．.“CutoffFrequency”置,“SampleRate”Fs);

信号をフィルター処理し,プロットします。結果は元の信号より滑らかですが,遅れます。

xf =过滤器(d, xn);情节(tn, xn)在情节(tn xf,“- r”，“线宽”, 1.5)从标题“心电图”包含“时间(s)”传奇(原始信号的，“过滤信号”）

grpdelayを使用して,フィルターによって発生する遅延がフィルター次数の半分に等しいことを確認します。

grpdelay (d, N, Fs)

延迟=意味着(grpdelay (d))

延迟= 35

フィルター処理された信号をシフトして,データを整列させます。最初の延迟サンプルを削除します。元の信号および時間ベクトルの最後の延迟サンプルを削除します。

tt = tn (1: end-delay);sn = xn (1: end-delay);科幻小说= xf;科幻小说(1:延迟)= [];

信号をプロットし,それらが整列していることを確認します。

情节(tt, sn)在情节(tt,科幻,“- r”，“线宽”, 1.5)从标题“心电图”包含(“时间(s)”),传说(原始信号的，“过滤信号转移”）

参考