关于fft计算平均周期
function T_mean=period_mean_fft(data)delt_t = 0.01;
t1 = 0:delt_t:500;
=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;
=max(xPower);
T_mean=period(index)
这个程序的问题存在哪里？蓝色部分对吗？计算结果应该约为4，按照上述程序运行为250。
页:
[1]