/math-e1e6230760a8fc9ebc5517bbe56756ea.png "\left\{f_i| i = 1,2,...,N\right\}")
は入力データポイント、/math-cbf2a8ae7d00ceaa3a438a13f080810a.png "\left\{g_i|i = 1,2,...,N\right\}")
は出力データポイントを表します。各/math-91d8947f9447751708c1f4107beb7ab3.png "g_i")
は/math-2d0012d8e86e7519e3d6bcbcc68b10db.png "\left\{f_m|i - floor(npts/2) < m < i + floor(npts/2)\right\}")
から計算されます。
スムージング法が移動平均、Savitzky-Golay、パーセンタイルフィルタの場合、各スムージングポイントは、移動ウィンドウ内データポイントで計算されます。 は入力データポイント、 は出力データポイントを表します。各はから計算されます。
ここでnptsは、変数のウィンドウのポイント数です。
しかし、スムージング法がFFTフィルタの場合、移動ウィンドウは使われません。代わりに、信号全体が処理されます。
隣接平均法は、簡単な平均計算を行います。各は移動ウィンドウ内のデータポイントの平均です。重み付け平均オプションにチェックが付いている場合、平均は加重平均を使って計算されます。この場合、パラボリック加重が使われ、重みの領域は1に正規化されます。中心が i にあるようなウィンドウ(つまり、i 番目の平均点を計算します)に対して、 j 番目のポイント(j=0, 1, …, npts-1) に対応する重みは次のようになります。
/math-fd200c5c25b7b7843b83b864f39daef2.png "w_j=1-(\frac{(j-i)}{(N+1)/2})^2")
ここでNは、ウィンドウのポイント数です。
Savitzky-Golay法によるフィルタリングは、移動ウィンドウ内のデータポイントに多項式回帰を実行します。そしては、i の位置で多項式の値として計算します。
パーセンタイルフィルタに対しては、移動ウィンドウ内のデータポイントのp番目の四分位がとして割り当てられます。 p は、 パーセンタイルパラメータで指定されます。 p番目の四分位(または100 p番目のパーセンタイル)は、平均を持つ経験分布関数から計算されます。
/math-7ae8e6c1c298971d77efadc3dad1752c.png "npts\cdot p/100=j+g")
のようにします。
ここで j は、 /math-a0c7f13b2f9f06d3210a4130a92e4caa.png "npts\cdot p/100")
の整数部であり、 g は の小数部分です。
そして、p番目の四分位を計算でき、これは次式のようにyで表されます。
/math-4aee7d1b3d1dbd0413eea75e0db28d4f.png "y=
\begin{cases}
x_j, & \mbox{if }g=0 \\
x_{j+1}, & \mbox{if }g>0
\end{cases}")
ここで、/math-1f89889020cdc84d9e1c35237cb62f65.png "x_j")
は、移動ウィンドウの j(j=0, 1, ...npts-1)番目のポイント
FFTフィルタ法が選択されていると、Originは次のように実行します。
/math-8a6ec61a846152b8540507019900a780.png "f_{cutoff} = \frac {1}{2n\Delta t}")
/math-5a72f1304af0783657605aed0e38201a.png "\Delta t")
は、任意の2点の時間(または、横軸)です。nを大きな値に設定すると、カットオフ周波数は低くなり、スムージングの値は大きくなります。変換されたデータは片側のウィンドウと乗算されるので、上の式はさらに2で除算し、両側ウィンドウを考慮します。|
Note:
|
Lowess と Loess は"locally weighted scatterplot smoothing" （局所加重散布図スムージング）と"locally weighted least squares"（局所加重最小二乗）の略語です。ここで「局所」と使うのは、ある範囲内でを隣り合ったポイントを元に各スムージング値を計算するためです。この方法は一般的に次の手順で実行されます。
/math-1ba8aaab47179b3d3e24b0ccea9f4e30.png "x_i")
の重みを計算します。次に立方体積比重関数で、その範囲内にある隣り合ったポイント全ての重みを計算します。
/math-d2de986db2dcdfead9130a94deacc948.png "w_i(x)=(1-(\frac{|x-x_i|}{d_i})^3)^3")
/math-9dd4e461268c8034f5c8564e155c67a6.png "x")
は現在の中心点に関連した範囲内の隣り合ったポイントを、/math-d247f594c78d0d2be10fc6d82512cc4e.png "d_i")
はその範囲内のから最も距離がある隣り合ったポイントの横座標距離を表わします。
/math-aeb3b53d344586a9c3e0ec0829777ae7.png "(x_i,\hat{y}_i)")
から、を導きます。
/math-84ceb2513893acb7c9ea001cc9053995.png "x_{i+1}")
まで移動し、ステップ1から3を実行して推定値/math-d845615daf118601c7a2d6719f7d165e.png "(x_{i+1},\hat{y}_{i+1})")
を入手します。全てのポイントの計算が終わると計算は終了します。
2項フィルタは重み付け隣接平均フィルタで、 /math-e447c37303aa003c6e1b9992045ddfda.png "{x_n}")
が入力ソースデータ、/math-83c66063e0d9b8864bd422ee8386d59f.png "{y_n}")
がスムージング出力データです。
/math-ceebc9dfaa7e3ac79481ea2adf795b92.png "y_n=\sum_{k=-N_p}^{Np}b_kx_{n-k}")
スムージング係数の順番/math-e23e5c6178fcede4de9a9673bf51cfd7.png "b_k")
は下記のように与えられます。
/math-da59a8657b40db814c19a032fb1f0f9b.png "b_k=\begin{pmatrix} 2N_p\\ N_p+k \end{pmatrix}/4^{N_p}\; \; (k=0,1,...N_p)")
および、
/math-d32b55686a350cbec96b05c8c8be3bbb.png "b_{-k}=b_{k}")
/math-766b775dcac93fa1cf6f6b365212665f.png "N_p")
は、次数です。
カットオフ周波数/math-e05fe30750d3ea262a610d17ebc07019.png "fc")
は次式で計算されます。
/math-42d1e7ee57aef355b30eddda869213a7.png "fc=\frac{2}{\pi}arccos(Ac^{1/2N_p})\frac{fs}{2}")
/math-54c355936eabc80c15d1dcbc2e697d0d.png "fs=\frac{1}{dt}")
/math-3017d911efceb27d1de6a92b70979795.png "dt")
はサンプリング間隔です。/math-8819a837186918b90b59cc316f36b1e1.png "Ac")
は、-6dBのカットオフ振幅で/math-57caf3aa6efe351c51c785be280b5c0d.png "Ac=0.5")
です。次数が増加すると、カットオフ周波数は減少します。