请教脉动风的验证

chanbill

:@) 我已经得出了0.25秒到256秒的脉动风速v,时间间隔为0.25秒，v有1024个点。现在要验证模拟的脉动风是否正确。即由脉动风得出功率谱，与目标功率谱比较。
第一种方法：
d=xcorr(v,'unbiased');%求自相关函数
s=fft(d,1024);
Pyy=abs(s)/512;
f=4*(0:511)/1024;
plot(f,Pyy(1:512))

第二种方法：
[power,freq]=psd(v,1024,4,boxcar(1024),0,'mean');
power=power*2*0.25;%规一化修正
plot(freq,power)
两种方法得出的图形都与目标功率谱有很大的差别，请问高手我的程序哪里有问题？谢谢！