为什么scipy.signal.welch会抑制零频率？

Question

我尝试使用whelch方法，我发现零频率异常

import numpy as np
import scipy.signal as signal
import matplotlib.pyplot as plt

n = 100000
s = np.ones(n)
f, psd = signal.welch(s, return_onesided=False)
plt.plot(f, psd)
plt.show()

因此结果PSD为零。 为什么scipy.signal.welch会抑制零频率？

Answer 1

再看看the docstring for welch特别注意detrend参数。

默认detrend是'constant'，它在计算频谱之前从输入中减去均值。要禁用去除趋势，请使用detrend=False：

In [57]: from scipy.signal import welch

In [58]: from scipy.fftpack import fftshift

In [59]: n = 1000

In [60]: s = np.ones(n)

In [61]: f, psd = signal.welch(s, return_onesided=False, detrend=False)

In [62]: plot(fftshift(f), fftshift(psd))
Out[62]: [<matplotlib.lines.Line2D at 0x10f8b6dd8>]