matlab画虫口模型分岔图,[转载]Matlab——混沌分岔图

经过近期的研究发现，目前对于系统单参数分岔图的计算共有以下的几种方法：

1)最大值法

即对系统微分方程(组)进行求解，对求解的结果用getmax函数进行取点，并绘图。

2)Poincare截面法

对系统参数的每一次取值，绘制其Poincare截面，进而得到其分岔图。

这种方法需要注意的是，自治系统的Poincare截面是选取一超平面，平面上点的分布即构成一Poincare截面，非自治系统的Poincare截面则是根据系统激励的频率进行取点并绘图。

本帖将以Lorenz系统为例，对这两种方法进行比较

首先对第二种方法进行阐述。

编程如下(matlab)

Lorenz系统：

function dy = Lorenz(t,y)

% Lorenz系统

% 系统微分方程：

%

dx/dt = -a(x-y)

%

dy/dt = x(r-z)-y

%

dz/dt = xy-bz

%

a=y(4)

%    r=y(5)

%

b=y(6)

dy=zeros(6,1);

dy(1)=-y(4)*(y(1)-y(2));

dy(2)=y(1)*(y(5)-y(3))-y(2);

dy(3)=y(1)*y(2)-y(6)*y(3);

dy(4)=0;

dy(5)=0;

dy(6)=0;​

随r的