ave_t.m

小数据量法求Lyapunov指数_龙格库塔

function P=ave_T(x,N)
Y=fft(x);
N=length(Y);
L=round(N/2);
Y(1)=[];
power=abs(Y(1:L)).^2;%%幅度的平方为能量
nyquist=1/2;
freq=(1:L)/(L)*nyquist;
%subplot(2,1,1);plot(freq,power);
%title('freq-power gram');
period=1./freq;
%subplot(2,1,2);plot(period,power);
%title('period-power gram');
[mp,index]=max(power);
P=period(index);
%end