如何在MATLAB上用二维数组绘制Lyapunov指数图

% 定义一个简单的Logistic映射
f = @(x, r) r * x * (1 - x);

% 参数设置
r_values = 2.8:0.01:4; % 控制参数r的范围
N = 1000; % 迭代次数
x0 = 0.5; % 初始条件

% 存储Lyapunov指数的数组
lyapunov_exponents = zeros(size(r_values));

% 计算Lyapunov指数
for i = 1:length(r_values)
    r = r_values(i);
    x = x0;
    lyapunov_sum = 0;
    
    for j = 1:N
        x_next = f(x, r);
        % 计算局部Lyapunov指数
        lyapunov_sum = lyapunov_sum + log(abs(r - 2*r*x));
        x = x_next;
    end
    
    % 平均Lyapunov指数
    lyapunov_exponents(i) = lyapunov_sum / N;
end

% 绘制Lyapunov指数图
figure;
plot(r_values, lyapunov_exponents, 'b-');
xlabel('r');
ylabel('Lyapunov Exponent');
title('Lyapunov Exponents for Logistic Map');
grid on;