计算湍流谱

Question

我对热线信号的评估遇到了麻烦。我正在风洞中记录一个热线信号（电压）来测量速度，更重要的是速度波动。所以从这个电压信号我可以计算速度信号 - 到目前为止没问题。关于流动中的湍流的重要陈述可以用能量密度谱完成。在数学上它被定义为

PSD = 2.0/N X*(f) X(f). 

X(f)是FFT，X*(f)是共轭复数FFT，N是我的样本长度。所以我做的是：

def powerspec(x):
    s = np.fft.fft(x)
    return 2.0/N*np.real(s*np.conjugate(s))

因为这不是导致正确的光谱（与给定的数据集相比），我尝试了不同的方法，如：

Pxx_den =np.abs(np.fft.fft(x))**2
f, Pxx_den = signal.periodogram(x, fs)

每个版本都会给出不同的结果，但没有一个是我需要的频谱。有没有人之前做过这件事可以提供帮助？差异如下图所示。

非常感谢您提前！

信号和＆＃34;右＆＃34;的数据可以在此处找到频谱：Link to Files with Signal and Spectrum