2017-06-14 54 views
0

任何人都可以請建議理想的窗口大小和重疊的樣本在Matlab中的pwelch函數。我有幾個200 ms的EEG信號,採樣率爲1000(信號長度或採樣數量= 200),以評估頻譜功率。默認情況下,pwelch使用漢明窗口並將數據分成8段,重疊度爲50%。對於只有200個樣本的信號,這些默認值是否正確?默認值工作正常,它給我一個PSD圖。我想確保我所做的是在概念上是正確的,如果有人可以提出任何更好的方法來做到這一點?這是對200 ms持續時間的ERP刺激反應的研究。我想比較不同頻段的頻譜功率。窗口大小和重疊PSD的評價matlab pwelch函數

感謝您的幫助!

回答

0

考慮到您的時間窗口只有200毫秒,我建議您使用pwelch,並使用一個200毫秒的窗口並且不要重疊。基礎傅里葉變換的頻率精度與所用數據段的長度有關。具體而言,傅立葉變換塊之間的間距與數據段的長度相關爲(Fs/N),其中Fs是數據的採樣率,並且N是段的長度。因此,例如採樣率爲1000 Hz和200 ms的數據段,頻率變換的有效分辨率將爲4 Hz。 (1000/200)。

還看到:https://electronics.stackexchange.com/questions/12407/what-is-the-relation-between-fft-length-and-frequency-resolution

如果你分析數據的一長段(比如,幾秒鐘),然後使用pwelch具有重疊窗口會更合適。

對基於頻率的分析的一個很好的介紹是M.X的文章'分析神經時間序列數據'。科恩。

+0

是否有任何拇指規則選擇窗口大小和重疊樣本?我們可以選擇頻率分辨率。因此,我們可以提供DFT點的數量爲500或1000或者我們想要的數量,以獲得所需數量的頻率點。所以我不明白Fs/N箱的意義。 – user8020776

+0

增加FFT點將填零數據,這將增加返回的FFT單元的數量,但不會像增加數據記錄的長度一樣增加解析相鄰頻率的能力。請參閱以下答案:https://dsp.stackexchange.com/questions/11560/why-do-we-say-that-zero-padding-doesnt-really-increase-frequency-resolution。根據一般規則,它取決於你想要分析的頻率以及你想要做出的結論。在時間和頻率分辨率之間進行權衡。 – DMR

+0

我試圖比較0-4Hz,4-8Hz,8-12Hz和12-30Hz的不同頻帶的譜功率,用於我的200點數據。 – user8020776