lyapunov指数的绘制

最新推荐文章于 2025-06-20 16:07:24 发布

原创

最新推荐文章于 2025-06-20 16:07:24 发布 · 552 阅读

2 ·

CC 4.0 BY-SA版权

文章标签：

#图像处理 #图像加密 #混沌系统

有如下方程：

%% 方程式
%  x(n+1)=1+y(n)-a*x(n)^2
%  y(n+1)=b*x(n)

绘制其对应的lyapunov指数。

在这里插入图片描述

MATLAB实现方式：


clc;
clearvars;
close all;

%% 方程式
%  x(n+1)=1+y(n)-a*x(n)^2
%  y(n+1)=b*x(n)



%% 代码
N = 1000; 
a = (0:0.001:1.4)'; 
b = 0.3; 
na = length(a); 
LE1 = zeros(na,1); 
LE2 = zeros(na,1); 
x = 0.2; y = 0.3; 
for i = 1:na 
    LCEvector = zeros(2,1); 
    Q = eye(2); 
    for j=1:N 
        xprev = x; 
        yprev = y; 
        x = 1-a(i)*xprev*xprev+yprev; 
        y = b*xprev; 
        Ji = [-2*a(i)*x,1;b 0];
        B = Ji*Q;
        [Q,R] = qr