
No.1ベストアンサー
- 回答日時:
ピリオドグラムというのは,つまり,確率的な信号についての,FFTの平均値みたいなものです.
普通,信号解析をする場合,多かれ少なかれ,外部ノイズなどの確率的な現象を含む信号(時系列)を対象にしています.
なんでスペクトル解析なんかをするときにも,確率的な取り扱いが入ってきます.
例えば,なんかの観測を無限回,繰り返したとします.
1回目は,時系列 f1(t) を観測,
2回目は,時系列 f2(t) を観測,
…
f1(t),f2(t),…は,だいたい似ているわけですが,細かく見るとノイズなんかの影響でちょっとずつ異なります.
なんで,f1(t),f2(t),…をそれぞれ,フーリエ変換(あるいはFFT)したF1(ω),F2(ω),…もちょっとずつ違います.
で,ピリオドグラムというのは,こいつらの各ωごとの平均のことです.つまり
Periodogram(ω) = E[ Fn(ω) ]_{n=1,2,…}
です.(つまり.確率空間全体にわたる平均のこと)
あと,パワースペクトルというのもあります.これは Fn(ω)の絶対値の2乗の平均のことです.
信号解析では,各周波数成分ごとのパワーのほうに興味があることが多いので,こっちのほうがよく使われます.
つまり,
PowerSpectrum(ω)= E[ |Fn(ω)|^2 ]_{n=1,2,…}
です,(これも,ほんとは,確率空間全体にわたる平均のこと)
注意しないといけないのは,
PowerSpectrum(ω) ≠ |Periodogram(ω)|^2
ということです.これは,ある確率変数Xについて
E[X^2] ≠ (E[X])^2
てことを考えればわかりますね.
なんで,パワースペクトルを計算するには,普通は観測した時系列を直接FFTするのではなくて,
時系列の自己相関をとってそれをFFTします.(ウィーナー・ヒンチンの定理)
定義上は,各観測時系列直接FFTして,それぞれ絶対値の2乗を計算してから,平均すればいいわけですが,
これだと,収束性がものすごく悪いです.
この回答へのお礼
お礼日時:2006/10/03 20:46
かみくだいてご説明いただき,助かりました。
FFT とピリオドグラムとで結果を比較してみました。
たしかにピリオドグラムは,ノイズ成分が非常に少ない!
でも拾いたい(高次の)周波数まで,低減されてました(笑
今後,ピリオドグラムを使うことがあれば,活用したいです。
とてもよい勉強になり今後,ツールの一つとして覚えておきたいです。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
±4σに入る確率について教えてく...
-
4人がじゃんけんしてあいこにな...
-
8頭身あって10人に1人くらいの...
-
数Aの質問です。
-
確率が重複する場合の計算方法
-
【確率の問題】
-
Cp値
-
条件付き確率で、Pa(B)とP(A∩B)...
-
1から13までの数字が1つずつ書...
-
ピリオドグラムって…
-
相関係数についてくるP値とは何...
-
発生確率0と見なせるのは?
-
確率の問題です。回答お願い致...
-
[数学] 無限大÷無限大の答えは?
-
同じ名前、生年月日の人同士が...
-
標準正規分布の確率を求める時...
-
可能性は「高い」?「大きい」?
-
論理的にこの確率の問題がわか...
-
確率はゼロだが可能性がないと...
-
丁半バクチの確率
おすすめ情報