%lorenz系统的映射分岔图 速度还可以
function Lorenz_bifur_c
Z=[];
for r=linspace(1,500,500);
% 舍弃前面迭带的结果，用后面的结果画图
[T,Y]=ode45('Lorenz',1,[1;1;1;10;r;8/3]);
[T,Y]=ode45('Lorenz',50,Y(length(Y),:));
Y(:,1)=Y(:,2)-Y(:,1);
% 对计算结果进行判断，如果点满足x=y，则取点
for k=2:length(Y)
f=k-1;
if Y(k,1)<0
if Y(f,1)>0
y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));
Z=[Z ,r+abs(y)*i];
end
else
if Y(f,1)<0
y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));
Z=[Z ,r+abs(y)*i];
end
end
end
end
plot(Z,'.','markersize',1)
title('Lorenz映射分岔图')
xlabel('c'),ylabel('|y| where x=y')
Bifurcations.zip
需积分: 9 49 浏览量
2021-11-15
16:45:31
上传
评论
收藏 3KB ZIP 举报
Bifurcations.zip （6个子文件） 
分岔图程序 getmax法 画图速度较慢 
Lorenz.M 305B getmax.m 773B lorenz映射的分岔图 lorenz7.m 819B Lorenz.M 305B 最大值法 画图速度较快 lorenz6.m 696B Lorenz.M 305B
评论0