求解系统的 Lyapunov 指数谱程序
Lyapunov 指数是描述时序数据所生成的相空间中两个极其相近的初值所产生的轨道,
随时间推移按指数方式分散或收敛的平均变化率。任何一个系统,只要有一个Lyapunov 大
于零,就认为该系统为混沌系统。
李雅普诺夫指数是指在相空间中相互靠近的两条轨线随着时间的推移,按指数分离或聚合的
平均变化速率。
一 chen 系统的 Lyapunov 指数谱
function dX = Chen2(t,X)
% Chen吸引子,用来计算Lyapunov指数
% dx/dt=a*(y-x)
% dy/dt=(c-a)*x+c*y-x*z
% dz/dt=x*y-b*z
global a; % 变量不放入参数表中
global b;
global c;
x=X(1); y=X(2); z=X(3);
% Y的三个列向量为相互正交的单位向量
Y = [X(4), X(7), X(10);
X(5), X(8), X(11);
X(6), X(9), X(12)];
% 输出向量的初始化
dX = zeros(12,1);
% Lorenz吸引子
dX(1) = a*(y-x);
dX(2) = (c-a)*x+c*y-x*z;
dX(3) = x*y-b*z;
% Lorenz吸引子的Jacobi矩阵
Jaco = [-a a 0;
c-a-z c -x;
y x -b];
dX(4:12) = Jaco*Y;
Z1=[];
- 1
- 2
前往页