印ライター：mpcsp079さん（最終更新日時:2016/5/6）投稿日：2012/12/1
.
■
デジタルフィルタで炉はする方法でーーーーーーーす。
■
ローパスフィルタ
Ｈ（ｓ）＝ａ／（ｓ＋ａ）
で炉はすることを考えます。ａはカットオフ周波数をｆ０とすれば２πｆ０です。
ｓ＝（２／T)（１－ｚ）／（１＋ｚ）
と、代入することでデジタルフィルタとなります。Ｔはサンプル周期です
Ｈ（ｚ）＝ａ／（（２／T)（１－ｚ）／（１＋ｚ）＋ａ）
＝（ａ＊（１＋ｚ））／（（２／Ｔ）＊（１－ｚ）＋ａ＊（１＋ｚ））
＝（ａ＋ａｚ）／（（２／Ｔ＋ａ）＋（ａー２／Ｔ）＊ｚ）
＝（１＋ｚ）／（（２／（Ｔａ）＋１）＋（１ー２／（Ｔａ））＊ｚ）ーー（１）
■
これをエクセルデータにあてはめます。
（１）式でＴａがエクセルデータ上でどうなるかがわかればいいのです。
エクセルの１コマをＴとし、正弦波の一周期のデータの数をＮとします。
すると、
周波数Ｆ＝１／（Ｔ＊Ｎ）
カットオフ周波数ｆ０をＦの５倍にしたければ、
ａ＝２π＊５／（Ｔ＊Ｎ）
Ｔａ＝２π＊５／Ｎ
（１）式にこのＴａ値をいれればＯＫです。
■
（１）式のエクセルでの計算法
Ｙ＊（（２／（Ｔａ）＋１）＋（１ー２／（Ｔａ））＊ｚ）
＝Ｘ＊（１＋ｚ）
より、
Ｙ（ｎ－１）＊（１ー２／（Ｔａ））＋Ｙ（ｎ）＊（２／（Ｔａ）＋１）
＝Ｘ（ｎ）＋Ｘ（ｎ－１）
Ｙ（ｎ）＝（Ｘ（ｎ）＋Ｘ（ｎ－１）－Ｙ（ｎ－１）＊（１ー２／（Ｔａ））／（２／（Ｔａ）＋１）
この式が答です。
Ｘ（ｎ）という数列は、測定データ列です。Ｙ（ｎ）は雑音の炉はされた結果
です
■
エクセルで検証しました。
左上元の正弦波、周期１００コマ
左下、ノイズ入り、周期３コマ
右下、上のデジタルフィルタ入れたもの。
ａ＝元正弦波周波数の５倍の周波数＊２π
です
