关于小数据量法求Lyapunov指数？？

yufeiyfyf

下面是一个小数据量法求Lyapunov指数程序：
clear
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h=0.01;
k=8000;
x(1)=2.2;
y(1)=2.4;
z(1)=28;
for i=1:k
    kx1=fx(x(i),y(i));
    ky1=fy(x(i),y(i),z(i));
    kz1=fz(x(i),y(i),z(i));
   
    kx2=fx(x(i)+(1/2)*h*kx1,y(i)+(1/2)*h*ky1);
    ky2=fy(x(i)+(1/2)*h*kx1,y(i)+(1/2)*h*ky1,z(i)+(1/2)*h*kz1);
    kz2=fz(x(i)+(1/2)*h*kx1,y(i)+(1/2)*h*ky1,z(i)+(1/2)*h*kz1);
   
    kx3=fx(x(i)+(1/2)*h*kx2,y(i)+(1/2)*h*ky2);
    ky3=fy(x(i)+(1/2)*h*kx2,y(i)+(1/2)*h*ky2,z(i)+(1/2)*h*kz2);
    kz3=fz(x(i)+(1/2)*h*kx2,y(i)+(1/2)*h*ky2,z(i)+(1/2)*h*kz2);
   
    kx4=fx(x(i)+h*kx3,y(i)+h*ky3);
    ky4=fy(x(i)+h*kx3,y(i)+h*ky3,z(i)+h*kz3);
    kz4=fz(x(i)+h*kx3,y(i)+h*ky3,z(i)+h*kz3);

    x(i+1)=x(i)+(h/6)*(kx1+2*kx2+2*kx3+kx4);
    y(i+1)=y(i)+(h/6)*(ky1+2*ky2+2*ky3+ky4);
    z(i+1)=z(i)+(h/6)*(kz1+2*kz2+2*kz3+kz4);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Data = x(1000:1:6000);
m=8;
tau=10;
N=5001;
P=ave_T(Data,N);
[Y_num,Y]=small_data_sets_Lyapunov1(Data,N,m,tau,P,h);
XX=1:Y_num;
plot(XX,Y(XX));hold on;
linear=80;%线性长度的变化引起最大Lapunov指数的变化，线性长度和最大Lapunov指数也受步长h的影响
XX1=1:linear;
pp=polyfit(XX1,Y(XX1),1);
lambda_1=pp(1)
yp=polyval(pp,XX1);
plot(XX1,yp,'r');

想问下这个线性长度是怎么确定的？上面程序里为什么选取linear=80？？说它的变化引起最大Lapunov指数的变化，线性长度和最大Lapunov指数也受步长h的影响？？是什么意思？？

cqupenghao

回复 1 # yufeiyfyf 的帖子

小数据量法求Lyapunov指数的过程中，在最后是用的最小二乘法对所得到的有些数据进行拟合，从而得到一个直线，直线的斜率就是Lyapunov指数的值，所以那个线性区域的选择对Lyapunov指数的值有很大影响……

yufeiyfyf

回复 2 # cqupenghao 的帖子

那这个线性长度是怎么确定的？

octopussheng

根据你最后计算的结果来定线性区域。