请教关于风速功率谱计算的问题

huaking

关于风速信号的功率谱的定义是什么？怎样计算？
我用傅里叶变换求得的功率谱整体偏高，程序如下：
clear
s1=load('*.dat');
n=length(s1);
%直接用FFT变换与它的转置求能量
y1=fft(s1,n);
Y1=y1.*conj(y1)/n;
Sn1=zeros(n/2,2);fc=50;
Sn1(:,2)=Y1(2:n/2+1);
%换算频率
for i=1:n/2
     Sn1(i,1)=fc*i/n;
end
dlmwrite('*.dat',Sn1);
这是什么问题呢？

CFDerSUN

你的纵坐标是幅值，不是功率谱。
你将纵坐标修改为Y1=abs(y1).^2/n/fc试一下