lorenz系统分岔图的三种算法

混沌Lorenz系统的matalb仿真，采用getmax、映射、最大值三种算法，绘出Lorenz系统的分岔图，分析系统的混沌时刻。

%lorenz系统分岔图 法2 速度很慢
function Lorenz_bifur_r_getmax
clear all
t0=[0 100];
 
%bifurcation
for r=linspace（1500500）; 
    r
    [ty]=ode45（‘Lorenz‘t0[1;1;1;10;r;8/3]）;
    [Xmax]=max（y（:1））;
a=length（y（:1））;
j=1;
k=1;
for l=1:a;
    YY= mod（y2）;
if（YY==1）
    XX（kj）=y（a1）;
k=k+1;
end
end
for i=（a-1）/2:a
    b=（y（i1）-y（i-21））/2;
    c=（y（i1）+y（i-21））/2-y（i-11）;
    if y（i-21）<=y（i-11）&&y（i-11）>=y（i1）&&c==0
        Xmax（j）=y（i-11）;
        j=j+1;
    elseif y（i-21）<=y（i-11）&&y（i-11）>=y（i1）
        Xmax（j）=y（i-11）-b^2/（4*c）;
        j=j+1;
    end
end
if j>1
   plot（rXmax‘b‘‘markersize‘1）
   hold on  
end
end
xlabel（‘随参数c变化‘）;
ylabel（‘x的最大值‘）;
title（‘lorenz bifurcation diagram‘）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        773  2011-04-09 12:44  分岔图程序\getmax法 画图速度较慢\getmax.m

     文件        305  2011-03-23 23:09  分岔图程序\getmax法 画图速度较慢\Lorenz.M

     文件        305  2011-03-23 23:09  分岔图程序\lorenz映射的分岔图\Lorenz.M

     文件        819  2011-04-09 21:31  分岔图程序\lorenz映射的分岔图\lorenz7.m

     文件        305  2011-03-23 23:09  分岔图程序\最大值法 画图速度较快\Lorenz.M

     文件        696  2011-04-09 12:51  分岔图程序\最大值法 画图速度较快\lorenz6.m

     目录          0  2016-10-28 19:46  分岔图程序\getmax法 画图速度较慢

     目录          0  2016-10-28 19:46  分岔图程序\lorenz映射的分岔图

     目录          0  2016-10-28 19:46  分岔图程序\最大值法 画图速度较快

     目录          0  2016-01-21 12:45  分岔图程序

----------- ---------  ---------- -----  ----

                 3203                    10