%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