最大lyapunov指数图

bohua1208

最大lyapunov指数图
不管改哪个参数
最大lyapunov指数图总是趋于0呢，一直没有大于0的情况怎么回事啊，我想要的是混沌运动的最大lyapunov指数图

function dX = DF(t,X)
a0=0.05;f1=5;f2=10;mu=0.1;w0=1;beta=0.1;ks=2;alafa=0.67;
b0=(f1-(f1^2-4*mu^2*w0^2*a0^2)^0.5)/(2*mu*w0);
omega=1;
x=X(1); y=X(2);
Y=[X(3),X(5);
   X(4),X(6)];   
dX = zeros(6,1);
dX(1)=y;
dX(2)=-(w0^2*x-ks*beta*x^2+ks*alafa*x^3+ks*mu*y-ks*f1*cos(omega*t)-ks*f2*cos(2*omega*t));      
J=[    0,    1;         
   -w0^2+2*ks*beta*x-3*ks*alafa*x^2,  -mu*ks];
dX(3:6) = J*Y;
% 计算Rossler吸引子的Lyapunov指数
clc;clear all;
a0=0.05;f1=5;f2=10;mu=0.1;w0=1;beta=0.1;ks=2;alafa=0.67;
b0=(f1-(f1^2-4*mu^2*w0^2*a0^2)^0.5)/(2*mu*w0);
yinit = [a0,b0];
orthmatrix = [1 0 ;
              0 1 ];   
y = zeros(6,1);
% 初始化输入
y(1:2) = yinit;
y(3:6) = orthmatrix;
tstart = 0; % 时间初始值
tstep = 1e-3; % 时间步长
wholetimes = 20000; % 总的循环次数
steps = 100; % 每次演化的步数
iteratetimes = wholetimes/steps; % 演化的次数
mod = zeros(2,1);
lp = zeros(2,1);
Lyapunov1 = zeros(iteratetimes,1);
Lyapunov2 = zeros(iteratetimes,1);
for i=1:iteratetimes
    tspan = tstart:tstep:(tstart + tstep*steps);   
    [T,Y] = ode45('DF', tspan, y);
    % 取积分得到的最后一个时刻的值
    y = Y(size(Y,1),:);
    % 重新定义起始时刻
    tstart = tstart + tstep*steps;
    y0 = [y(3) y(5);
          y(4) y(6)];   
    %正交化
    y0 = TwoGS(y0);
    mod(1) = sqrt(y0(:,1)'*y0(:,1));
    mod(2) = sqrt(y0(:,2)'*y0(:,2));
    y0(:,1) = y0(:,1)/mod(1);
    y0(:,2) = y0(:,2)/mod(2);
   
    lp = lp+log(abs(mod));
  
    Lyapunov1(i) = lp(1)/(tstart);
    Lyapunov2(i) = lp(2)/(tstart);
  
        y(3:6) = y0';
end
% 作Lyapunov指数谱图
i = 1:iteratetimes;
plot(i,Lyapunov1,i,Lyapunov2)

shengjunjie

为何，我用这个程序计算出错呢？??? Undefined function or method 'TwoGS' for input arguments of type 'double'.

Error in ==> Lyapunov at 31
     y0 = TwoGS(y0);

lihaitao123

shengjunjie 发表于 2012-6-21 11:01

                               
登录/注册后可看大图

为何，我用这个程序计算出错呢？??? Undefined function or method 'TwoGS' for input arguments of type ' ...

http://www.chinavib.com/forum.ph ... mp;page=2#pid323869