Microsoft Excel に、時系列のデータが入っています。
2次微分までのデータが欲しいのですが、ノイズが多く、
きれいな曲線になりませんでした。
このような事情で、Savitzky Golayの方法(名前は知っていた)で
データのスムージングを行おうとしていたところ、あるサイトで、
{0, 0.0323313, 0.00850822, -0.00829007,..., 0.0859867,...,0.0323313}
というような数列のようなものを見つけました。
これを係数に41項からなる和算を行ったところ、全く何も根拠のないまま、
一見したところ目的を達しているかのような曲線が得られました。
これ、全くのデタラメでもないのでしょうか? それに、もう少し項目を
減らしたり増やしたりするには、どんな係数を使えばいいのでしょうか。
そもそも、適切な係数を知るにはどこを調べればいいのでしょうか。
No.1ベストアンサー
- 回答日時:
【原典】
Savitzky-Golay法は次の文献で発表されています。
A.Savitzky,M.J.E.Golay"Smoothing and Differention of Data bye Simplified Least Squares Procedures,"Analytical Chemistry,vol.36,no.8,pp1627-1639,1964
【方法の主眼】
等間隔で得られた観測値から雑音を除去するために最小2乗法
を用いて多項式に当てはめる。
観測点が"等間隔"言うのがミソです。
たとえば多項式を2次式として、平滑化の対象を奇数点の(-2,-1,0,2,1)
という点を選んでスケール変換して計算すると、重み係数(-3,12,17,12,-3)、
正規化定数35が得られます。
【得られる情報源】
科学計測のための波形データ処理 南茂夫 CQ出版社
1986年 初版 ¥1960
この本には、具体的なデータ処理が例示されています。
問われている係数の算出方法についても式が提示されています。
原典よりも新しいですが、まだ出版されているのかな?
【展望】
観測値から雑音を除去すると言うテーマは多くの分野で
あつかわれているようです。
Savitzky-Golay法は、クラシックな部類に属する手法かと思います。
しかし、科学計測の分野では根強く活躍しているようです。
数値フィルタまたはディジタルフィルタという観点から見ると、スペクトル
というものの見方が欠けているのかなと思われます。
最新号のC-MAGAZINEの特集またはその参考文献がきっと
役にたつと思います。
【最後に】
ご質問のテーマは、きっと物理または化学が適切では
なかったかな、と思います。
741さん
『科学計測のための波形データ処理』はオンラインでも入手可能の
ようですので、さっそく購入することにしました。
初めは係数の算出方法を知るだけが目的でしたが、
面白そうな分野を知るいいきっかけになりました。
あまり深く理解すると、他のもっと適切な方法が見えてきて
大変になるのかも知れませんが、少しでも理解して使えるように
がんばってみます。
おっしゃるとおり、質問のカテゴリーをもう少し考えるべきでした。
が、それでも論文と書籍の情報を教えていただけたことは、幸運でした。
貴重なお答え、ありがとうございました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Visual Basic(VBA) 3つのプロシージャをまとめたら実行時エラー発生で対応不能 6 2022/05/17 01:47
- Visual Basic(VBA) vba 等間隔の列に対しての計算 6 2022/05/17 20:15
- Visual Basic(VBA) 【VBA】データを入力後に,同一シート内に履歴として転記するVBAコードを教えていただきたいです。 3 2022/11/16 01:37
- 物理学 微分方程式の物理現象への適用について 3 2023/05/14 12:22
- Excel(エクセル) VBAで重複データを合算したい(時間) 1 2022/12/08 23:06
- その他(IT・Webサービス) 高速処理可能な表計算ソフトについて ExcelやGoogleスプレッドシートのような表計算ソフトで、 2 2023/04/29 16:06
- Excel(エクセル) Countifよりも早く重複数をカウントする方法ありますか? 18 2022/07/04 13:39
- システム CSVファイルのマッピング処理の省力化 1 2022/11/24 00:01
- Excel(エクセル) VLOOKUP が機能しない、その原因は何 ? 8 2022/10/19 12:06
- Visual Basic(VBA) VBAで大量データの処理 3 2022/11/15 21:53
このQ&Aを見た人はこんなQ&Aも見ています
-
新NISA制度は今までと何が変わる?非課税枠の拡大や投資対象の変更などを解説!
少額から投資を行う人のための非課税制度であるNISAが、2024年に改正される。おすすめの銘柄や投資額の目安について教えてもらった。
-
EXCELでSavitzky-Golay法
Excel(エクセル)
-
Excelで、Savitzky Golay- データのノイズ除去 - を行うには?
Excel(エクセル)
-
手軽に曲線の平滑化をしたい
Excel(エクセル)
-
-
4
savitzky-golayのフィルタについて
物理学
-
5
EXCELにてローパスフィルタを作成する
その他(教育・科学・学問)
-
6
エクセルで極大値を拾うには
Excel(エクセル)
-
7
エクセルでノイズ値を除去する方法。
数学
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
e^sinXの展開式について。。。
-
可算個の不連続点をもつ関数の...
-
約数と因数の違い
-
約数と因数の違い(∈N)
-
多項式について質問です。 エク...
-
これがどうしても分かりません❗...
-
データのノイズ除去法 - Savitz...
-
(x-1)(x-2)(x-3)の展開の...
-
エルミート補間について
-
多項式とはなんですか?
-
deg f?
-
3次式と2次式の最大公約数の問題
-
単項式とは
-
原始多項式の求め方
-
問題が理解できません
-
三乗根を含んだ最小多項式
-
(x+y+2z)(2x+3y-z)(4x-y-3z)を...
-
単項式・多項式
-
(1)ニューグローバルβの190につ...
-
テーラー展開を用いて0の近く...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
多項式について質問です。 エク...
-
単項式と分数式の違いについて
-
斉次とは?(漢字と意味)
-
データのノイズ除去法 - Savitz...
-
阪大2014年数学挑戦枠2問からで...
-
余次元って何?
-
(x-1)(x-2)(x-3)の展開の...
-
数学 因数分解 X^3+x^2+x−1 ...
-
(x+y+2z)(2x+3y-z)(4x-y-3z)を...
-
約数と因数の違い(∈N)
-
数を拡張するとはなんですか? ...
-
等差×等比 型の数列の和を求め...
-
arcsinのマクローリン展開について
-
ローラン展開についてです。
-
CRCのアルゴリズムって、どんな...
-
なぜ、2変数以上の多項式を因数...
-
deg f?
-
0は偶関数?
-
原始多項式の求め方
-
テイラー展開の剰余項
おすすめ情報