用于 PSD 评估的 MATLAB 的 pwelch 函数中的窗口大小和重叠



任何人都可以建议理想的窗口大小和 Matlab 中 pwelch 函数的重叠样本。我有几个采样率为1000(信号长度或样本数= 200)的200毫秒EEG信号来评估频谱功率。默认情况下,pwelch 使用汉明窗口并将数据分成 8 个重叠 50% 的段。对于只有 200 个样本的信号,这些默认值可以吗?默认值工作正常,它给了我一个 PSD 图。我想确保我正在做的事情在概念上是正确的,如果有人可以提出更好的方法来做到这一点?这是一项关于ERP在200毫秒持续时间内对刺激的反应的研究。我想比较不同频段的频谱功率。

感谢您的帮助!

考虑到您的时间窗口只有 200 毫秒,我建议使用具有单个 200 毫秒窗口且没有重叠的pwelch。基础傅里叶变换的频率精度与所用数据段的长度有关。具体来说,傅里叶变换箱之间的间距与数据段的长度相关,如(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. Cohen的"分析神经时间序列数据"。

最新更新