计算湍流谱

2024-06-26 16:51:23 发布

您现在位置:Python中文网/ 问答频道 /正文

我在评估热线信号时遇到了麻烦。我正在风洞中记录一个热线信号(电压),它测量速度,更重要的是速度波动。所以根据这个电压信号,我可以计算出速度信号-到目前为止没问题。可以用能量密度谱对湍流进行重要的描述。数学上它被定义为

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

其中X(f)是FFT,X*(f)是共轭复数FFT,N是我样本的长度。所以我所做的是:

^{pr2}$

因为这不会导致正确的频谱(与给定的数据集相比),我尝试了不同的方法,如:

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

每个版本都给出了不同的结果,但没有一个是我需要的频谱。有没有人以前做过这件事,能帮上忙吗?差异如下图所示。在

提前非常感谢!在

“右”谱的信号和数据可以在这里找到:Link to Files with Signal and Spectrum

Difference between my try and the "right" solution.


Tags: 数据fft信号np记录速度频谱能量