请教关于风速功率谱计算的问题
关于风速信号的功率谱的定义是什么?怎样计算?我用傅里叶变换求得的功率谱整体偏高,程序如下:
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);
这是什么问题呢? 你的纵坐标是幅值,不是功率谱。
你将纵坐标修改为Y1=abs(y1).^2/n/fc试一下
页:
[1]