关于fft计算平均周期

wyk_1216

function T_mean=period_mean_fft(data)
delt_t = 0.01;
t1 = 0:delt_t:500;
[tt1,y1]=ode45(@Lorenz,t1,[-1,0,1]);
xx1 = y1(:,1)';
x1 = spline(tt1, xx1, t1);
data = x1(30001:1:40000);
x=data;
xPower=abs(fft(x));
NN=length(xPower);
xPower(1)=[];
NN=floor(NN/2);
xPower=xPower(1:NN);
freq=(1:NN)/(NN*0.5);
period=1./freq;
[mP,index]=max(xPower);
T_mean=period(index)

这个程序的问题存在哪里？蓝色部分对吗？  计算结果应该约为4，按照上述程序运行为250。