psd_multitaper与psd_welch的频率步进

Question

psd_welch返回值中的频率步进为1.953125，而在psd_multitaper中为0.00164441 ~3个数量级差：

import mne
psds,freqs = mne.time_frequency.psd_multitaper(edf)
np.diff(freqs)
==> array([0.00164441, 0.00164441, 0.00164441, 0.00164441, 0.00164441, ...])
psds,freqs = mne.time_frequency.psd_welch(edf)
np.diff(freqs)
==> array([1.953125, 1.953125, 1.953125, 1.953125, 1.953125, 1.953125, ...])

为什么？我该怎么控制它？

PS。另见讨论。

Answer 1

psd_multitaper -不去

频率由edf.info["sfreq"]确定，只有edf.resample()才能改变频率，并修改了频率上限edf.info["lowpass"]。

psd_welch -简单

使用n_fft参数。默认为256。因此，如果要将步长从1.953125减半到0.9765625，请传递n_fft=512。