二维离散动力学系统的混沌研究【基于matlab的动力学模型学习笔记_9】

歪卜巴比

已于 2022-07-29 11:26:56 修改

阅读量5.1k

点赞数 9

分类专栏：基于matlab的动力学模型研究文章标签： matlab 线性代数矩阵动力学

于 2021-11-10 19:43:54 首次发布

本文链接：https://blog.csdn.net/L1234X/article/details/121254801

版权

本文通过MATLAB研究了二维动力学系统，包括厄农映射和Lotka-Volterra捕食模型的混沌行为。详细探讨了最大Lyapunov指数及正平衡点的求解，展示了如何从一维系统过渡到二维混沌系统，并通过实例分析系统参数变化对混沌状态的影响。

摘要由CSDN通过智能技术生成

摘要：混沌（Chaos）是指发生在确定系统中的貌似随机的不规则运动，本文将基于经典的二维系统，然后根据动力学方程研究其混沌产生过程以及相对应的MATLAB仿真，再讨论Lyapunov指数以及正平衡点。

上一篇中介绍了一维系统，这次我们将维数提升到二。

/*仅当作学习笔记，若有纰漏欢迎友好交流指正，此外若能提供一点帮助将会十分荣幸*/

2 Lotka-Volterra捕食模型的混沌

0 引言

对于一维离散混沌系统来说，其具有形式简单、产生混沌序列时间短等优点。但其缺点是空间太小，难以处理较为复杂的情况。因此本文将从经典的二维混沌系统入手以此为基础，研究动力学系统的发展、混沌的产生以及相应的MATLAB仿真，借此建立一个对二维混沌系统有更为完整的认识与理解。

1 厄农（Henon）映射

1.1系统的混沌

厄农映射是由法国数学家米歇尔·厄农提出的，以此作为洛伦茨模型的庞加莱截面的简化模型。对经典厄农映射而言，任意初始点或趋向厄农奇异吸引子，或发散至无穷大。厄农吸引子具有分形结构，其在一个方向上连续，另一个方向上则为一个康托尔集，其方程动力学微分方程为：

其中：a、b为系统的参数分别取为0.4、0.3，其对应MATLAB程序如下：

%厄农映射
%x(j+1)=1+yj-a*xj*xj
%y(j+1)=b*xj
clc
clear all
close all
a=0.4;b=0.3;%系统参数设置
xjL=[];yjL=[];
for k=1:140%迭代次数
    a=0.5+(k-1)*0.01;
xj=[];yj=[];
    x=0.1;y=0.1;
for j=1:200
        x=1+y-a*x*x;y=b*x;%系统函数
xj=[xj;j x];  
yj=[yj;j y];
if j>=120
xjL=[xjL;a b x];    
yjL=[yjL;a b y]; 
end
end
end
plot(xjL(:,1),xjL(:,3),'k.')
%坐标轴设置
xlabel('\itt','FontSize',20,'FontName','Times New Roman');
ylabel('\itx','FontSize',20,'FontName','Times New Roman');
grid on;

仿真结果：