如何计算频带内的功率谱?

How can I compute the power spectrum within a frequency band?

我有离散的 EEG 信号,我正在尝试从每个通道中提取绝对功率。采样率为 30Hz。我感到困惑的是我不知道如何准确地编程。 首先我做了:

y = fft(eegdata);

这给了我一个复数。我必须在频带 (alpha) + 汉明 window + 50% 重叠范围内计算这个,但我不知道该怎么做。谁能帮我解决这个问题?

试试 Matlab 的 pwelch 函数。如果您的单位是 U(例如),则 pwelch 输出 U^2/(Hz s),因此如果您想要 U^2/Hz 中的功率谱密度,请将 pwelch 的结果乘以您的采样间隔(在你的情况下,1/30 秒)。

默认设置应用具有 50% 重叠的汉明 window。