Lyapunov指数法判断是否具备混沌特性，写出matlab代码

以下是一个简单的matlab代码，用于计算一维动力系统的Lyapunov指数：

% 定义一维动力系统
function x_dot = dyn_sys(x)
    a = 2.6;
    b = 0.1;
    x_dot = -a*x + b*sin(x);

% 初始化参数
x0 = 0.1; % 初始状态
n = 10000; % 迭代次数
dt = 0.01; % 时间步长

% 初始化Lyapunov指数
L = 0;

% 算法主体
for i = 1:n
    % 计算系统状态
    x = dyn_sys(x0);
    
    % 计算Lyapunov指数
    L = L + log(abs(dyn_sys(x0) * dt));
    
    % 更新状态
    x0 = x;
end

% 计算平均Lyapunov指数
L = L / n;

% 输出结果
fprintf('Lyapunov指数为: %f\n', L);

这个代码中的动力系统是一个简单的一维系统，但是它可以很容易地被修改以适应其他系统。