lorenz的分岔图，大家快看一下吧

wangjianbao11 · 发表于 2011-4-29 22:45

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

我做的lorenz分岔图，大家帮忙分析一下吧

meiyongyuandeze · 发表于 2011-4-29 23:10

回复 1 # wangjianbao11 的帖子

用的是什么方法呢，洛伦茨系统不会出现这样的间断吧！

wangjianbao11 · 发表于 2011-4-30 08:01

我用的是极大值法，得到了上边的图像，a=10，b=8/3,c=[0,30],我做的理论推导，也发现是对的，具体是什么图形，您能帮忙分析一下吗？谢谢了大哥哥

meiyongyuandeze · 发表于 2011-4-30 19:13

本帖最后由 meiyongyuandeze 于 2011-4-30 19:42 编辑

回复 3 # wangjianbao11 的帖子

感觉不应该出现那样的间断的情况啊！论坛中有很多关于分岔图的做法你可以参考下,洛伦兹系统是比较典型的非线性系统，应该程序比较的多，你试下庞加莱界面方法，截取像平面X=Y试试！
【总结帖】分岔图绘制不同方法的总结、比较
http://forum.vibunion.com/forum-viewthread-tid-60104-fromuid-77817.html

wangjianbao11 · 发表于 2011-4-30 19:26

大哥，我看了您推荐的方法，但我的程序有什么缺陷，或者说它有什么不对，麻烦您一下，能帮我看一下吗？我的程序运行结果就是如此。

meiyongyuandeze · 发表于 2011-4-30 19:42

你把程序贴出来，让大家一起帮你看看吧！

wangjianbao11 · 发表于 2011-4-30 21:21

function dx=lorze1(t,X)
x=X(1);

y=X(2);

z=X(3);

a=10;

b=8/3;

global c;

dx=zeros(3,1);

dx(1)=-a*(x-y);

dx(2)=-x*z+c*x-y;

dx(3)=-b*z+x*y;

clear;
global c;
range=[0:0.1:30];
k=0;
YY=[];
for c=range
c
   y0=[0.05;0.05;0.05]
   %rand(1,3);
   k=k+1;
tspan=[0:0.1:50];
   [t,Y]=ode45('lorze1',tspan,y0);
count=find(t>0);
Y=Y(count,:);
% 画x的分岔图。
j=1;
   n=length(Y(:,1));
for i=2:n-1
      if Y(i-1,1)+eps<Y(i,1) && Y(i,1)>Y(i+1,1)+eps
% 简单的取出局部最大值。
         YY(k,j)=Y(i,1);
         j=j+1;
      end
   end
if j>1
      plot(c,YY(k,[1:j-1]),'k.','markersize',1);
   end
hold on;
index(k)=j-1;
end
xlabel('r');
ylabel('x max');

wangjianbao11 · 发表于 2011-5-3 09:17

给位帮忙看一下吧

meiyongyuandeze · 发表于 2011-5-3 09:48

用最大方法可能会导致画出的图不连续。可以采用固定截面方法来取点！

wangjianbao11 · 发表于 2011-5-3 09:54

谢谢了师哥，但是用庞加莱截面法我连续做了好几天，还是做不成，师哥您能帮我一下吗，教我一下程序该怎么写吗

账号		自动登录	找回密码
密码			我要加入

[动力学和稳定性] lorenz的分岔图，大家快看一下吧

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。