周波数帯域内のパワースペクトルをどのように計算できますか？

Question

私は離散EEG信号を持っており、各チャンネルから絶対パワーを抽出しようとしています。サンプリングレートは30Hzです。私が混乱していることは、これを正確にプログラムする方法がわからないことです。 は、まず私が行っている：

y = fft(eegdata); 

これは私に複雑な二重を与えます。私は周波数帯域（アルファ）+ハミング窓+ 50％オーバーラップでこれを計算しなければならないが、これをどうやって行うのか分からない。誰もがこれで私を助けることができますか？

Answer 1

Matlabのpwelch関数を試してください。あなたのユニットがU（例えば）の場合、pwelchはU^2 /（Hz s）を出力するので、パワースペクトラル密度をU^2/Hzで求めたい場合、pwelchの結果にサンプル間隔を掛けます、1/30秒）。

デフォルト設定では、50％オーバーラップのハミングウィンドウが適用されます。