赋个初值看看:
我们使用雅克比矩阵计算lyapunov exponent图谱看看:
matlab代码如下:
clear;
%% ===========初始化输入==============
yinit = [0.1,0.1,0.1];
orthmatrix = [1 0 0;
0 1 0;
0 0 1];
y = zeros(12,1);
y(1:3) = yinit;
y(4:12) = orthmatrix;
mod = zeros(3,1);
lp = zeros(3,1);
%% ============迭代控制==============
tstart =