Main Content

このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。

filtfilt

ゼロ位相デジタル フィルター処理

ページ内をすべて折りたたむ

構文

y = filtfilt(b,a,x)
y = filtfilt(sos,g,x)
y = filtfilt(d,x)

説明

y = filtfilt(b,a,x) は、入力データ x を順方向と逆方向の両方で処理することにより、ゼロ位相デジタル フィルター処理を実行します。データを順方向でフィルター処理した後、関数は初期条件を一致させて起動時と終了時の過渡現象を最小限に抑え、フィルター処理したシーケンスを反転し、その反転したシーケンスを逆方向にフィルター処理します。結果には、以下の特性があります。

  • ゼロ位相の歪み

  • 元のフィルター伝達関数の振幅の 2 乗と等しい、フィルター伝達関数

  • b と a により指定されるフィルターの次数の倍のフィルター次数

filtfilt は、Gustafsson による推奨アルゴリズム[1]を実装します。

微分器フィルターやヒルベルト FIR フィルターの演算は位相応答に大きく依存しているため、関数 filtfilt をこれらのフィルターと共に使うことは避けてください。

y = filtfilt(sos,g,x) ゼロ位相は、行列 sos およびスケール値 g で表される 2 次セクション (双二次) フィルターを使用して入力データ x をフィルター処理します。

例

y = filtfilt(d,x) ゼロ位相は、デジタル フィルター d を使用して入力データ x をフィルター処理します。d を周波数応答仕様に基づいて生成するには、関数 designfilt を使用します。

例

すべて折りたたむ

ライブ スクリプトを開く

ゼロ位相フィルター処理は、フィルター処理された時間波形の特徴を、フィルター処理されていない信号で現れる場所にそのまま保持する上で役立ちます。

合成心電図 (ECG) 波形をゼロ位相フィルター処理します。波形を生成する関数は、例の最後にあります。QRS 群は心電図波形の重要な特徴です。この例では時点 160 あたりから始まっています。

wform = ecg(500);

plot(wform)
axis([0 500 -1.25 1.25])
text(155,-0.4,"Q")
text(180,1.1,"R")
text(205,-1,"S")

Figure contains an axes object. The axes object contains 4 objects of type line, text.

心電図を加法性ノイズで乱します。再現可能な結果が必要な場合は、乱数発生器をリセットします。ローパス FIR 等リップル フィルターを作成し、ゼロ位相のフィルターと通常のフィルターの両方を使用してノイズの多い波形をフィルター処理します。

rng default

x = wform' + 0.25*randn(500,1);
d = designfilt("lowpassfir", ...
    PassbandFrequency=0.15,StopbandFrequency=0.2, ...
    PassbandRipple=1,StopbandAttenuation=60, ...
    DesignMethod="equiripple");
y = filtfilt(d,x);
y1 = filter(d,x);

subplot(2,1,1)
plot([y y1])
title("Filtered Waveforms")
legend("Zero-phase Filtering","Conventional Filtering")

subplot(2,1,2)
plot(wform)
title("Original Waveform")

Figure contains 2 axes objects. Axes object 1 with title Filtered Waveforms contains 2 objects of type line. These objects represent Zero-phase Filtering, Conventional Filtering. Axes object 2 with title Original Waveform contains an object of type line.

ゼロ位相フィルター処理では信号のノイズが低減し、QRS 群は元々の発生と同じ時点に保持されます。通常のフィルター処理では信号のノイズは低減しますが、QRS 群に遅延が生じます。

バタワース 2 次セクション フィルターを使用してフィルター処理を繰り返します。

d1 = designfilt("lowpassiir",FilterOrder=12, ...
    HalfPowerFrequency=0.15,DesignMethod="butter");
y = filtfilt(d1,x);

subplot(1,1,1)
plot(x)
hold on
plot(y,LineWidth=3)
legend("Noisy ECG","Zero-Phase Filtering")

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Noisy ECG, Zero-Phase Filtering.

この関数は ECG 波形を生成します。

function x = ecg(L)
%ECG Electrocardiogram (ECG) signal generator.
%   ECG(L) generates a piecewise linear ECG signal of length L.
%
%   EXAMPLE:
%   x = ecg(500).';
%   y = sgolayfilt(x,0,3); % Typical values are: d=0 and F=3,5,9, etc. 
%   y5 = sgolayfilt(x,0,5); 
%   y15 = sgolayfilt(x,0,15); 
%   plot(1:length(x),[x y y5 y15]);

%   Copyright 1988-2002 The MathWorks, Inc.

a0 = [0,1,40,1,0,-34,118,-99,0,2,21,2,0,0,0]; % Template
d0 = [0,27,59,91,131,141,163,185,195,275,307,339,357,390,440];
a = a0/max(a0);
d = round(d0*L/d0(15)); % Scale them to fit in length L
d(15)=L;

for i=1:14
       m = d(i):d(i+1)-1;
       slope = (a(i+1)-a(i))/(d(i+1)-d(i));
       x(m+1) = a(i)+slope*(m-d(i));
end

end

入力引数

すべて折りたたむ

伝達関数の係数。ベクトルで指定します。全極フィルターを使用する場合は、b に対し 1 を入力します。全零 (FIR) フィルターを使用する場合は、a に対し 1 を入力します。

例: b = [1 3 3 1]/6 および a = [3 0 1 0]/3 は、正規化された 3 dB の周波数 0.5π ラジアン/サンプルをもつ 3 次のバタワース フィルターを指定します。

データ型: single | double

入力信号。実数値または複素数値のベクトル、行列、または N 次元配列として指定します。x は有限の値でなければなりません。x の長さは、max(length(B)-1, length(A)-1) として定義されるフィルター次数の 3 倍より大きくなければなりません。関数は、x が行ベクトルでなければ、x の最初の配列次元に沿って動作します。x が行ベクトルである場合は、2 番目の次元に沿って動作します。

例: cos(pi/4*(0:159))+randn(1,160) は単一チャネルの行ベクトル信号です。

例: cos(pi./[4;2]*(0:159))'+randn(160,2) は 2 チャネル信号です。

データ型: single | double
複素数のサポート: あり

2 次セクションの係数。行列として指定します。sos は K 行 6 列の行列で、セクション数 K が 2 以上でなければなりません。セクション数が 2 未満の場合、関数は入力を分子ベクトルとして扱います。sos の各行は 2 次 (双二次) フィルターの係数に対応しています。sos の i 行目は [bi(1) bi(2) bi(3) ai(1) ai(2) ai(3)] に対応しています。

例: s = [2 4 2 6 0 2;3 3 0 6 0 0] は、正規化された 3 dB の周波数 0.5π ラジアン/サンプルをもつ 3 次のバタワース フィルターを指定します。

データ型: single | double

スケール係数。ベクトルとして指定します。

データ型: single | double

デジタル フィルター。digitalFilter オブジェクトで指定します。デジタル フィルターを周波数応答仕様に基づいて生成するには、関数 designfilt を使用します。

例: d = designfilt("lowpassiir",FilterOrder=3,HalfPowerFrequency=0.5) は、正規化された 3 dB の周波数 0.5π ラジアン/サンプルをもつ 3 次のバタワース フィルターを指定します。

出力引数

すべて折りたたむ

フィルター処理された信号。ベクトル、行列、または N 次元配列として返されます。

filtfilt への入力が単精度の場合、関数は出力 y を単精度で返します。

参照

[1] Gustafsson, F. “Determining the initial states in forward-backward filtering.” IEEE® Transactions on Signal Processing. Vol. 44, April 1996, pp. 988–992. https://doi.org/10.1109/78.492552.

[2] Mitra, Sanjit K. Digital Signal Processing. 2nd Ed. New York: McGraw-Hill, 2001.

[3] Oppenheim, Alan V., and Ronald W. Schafer, with John R. Buck. Discrete-Time Signal Processing. 2nd Ed. Upper Saddle River, NJ: Prentice Hall, 1999.

拡張機能

バージョン履歴

R2006a より前に導入

参考

| | | |

トピック