【总结帖】分岔图绘制不同方法的总结、比较

zhong124 · 发表于 2011-3-27 20:16

回复 7 # liliangbiao 的帖子

我一直都没画出来，你可不可以分享一下啊？

netor · 发表于 2011-3-30 15:54

图挂了……

wangjianbao11 · 发表于 2011-4-29 21:44

这是我做的lorenz分岔图，希望给位大哥哥帮忙看看，看做的对吗？怎么去分析它？

wangjianbao11 · 发表于 2011-4-29 21:50

这是我做的lorenz分岔图，希望给位大哥哥帮忙看看，看做的对吗？怎么去分析它？
我选的是a=10，b=8/3

xiaozh99 · 发表于 2011-11-28 11:46

回复 7 # liliangbiao 的帖子

能否给出相关程序

0810064 · 发表于 2011-12-7 12:58

真是强人得好好向各位学习啊

ldlnihao · 发表于 2011-12-10 21:29

http://forum.vibunion.com/thread-108009-1-1.html

bujin · 发表于 2011-12-19 23:39

好东西，谢谢分享

伤痕累累 · 发表于 2012-3-28 20:49

function [Xmax] = getmax(y)
a=length(y);
j=1;
for i=(a-1)/2:a
b=(y(i,1)-y(i-2,1))/2;
c=(y(i,1)+y(i-2,1))/2-y(i-1,1);
if y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
      Xmax(j)=y(i-1,1);
      j=j+1;
elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
      Xmax(j)=y(i-1,1)-b^2/(4*c);

里面elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
      Xmax(j)=y(i-1,1)-b^2/(4*c);

为什么要减去b^2/(4*c);表示不明白。为什么这里的getmax函数是这么写得。换成其他样的不行吗。

伤痕累累 · 发表于 2012-3-28 20:53

function [Xmax] = getmax(y)
a=length(y);
j=1;
for i=(a-1)/2:a
b=(y(i,1)-y(i-2,1))/2;
c=(y(i,1)+y(i-2,1))/2-y(i-1,1);
if y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
      Xmax(j)=y(i-1,1);
      j=j+1;
elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
      Xmax(j)=y(i-1,1)-b^2/(4*c);

里面elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
      Xmax(j)=y(i-1,1)-b^2/(4*c);

为什么要减去b^2/(4*c);表示不明白。为什么这里的getmax函数是这么写得。换成其他样的不行吗。

yiyewu.mian · 发表于 2012-4-19 15:31

最大值法和Poincare截面法哪个方法效果好一些？哪个迭代更快些

mybook86 · 发表于 2012-6-13 16:05

如枫 · 发表于 2013-1-28 22:28

请教一下，
Z=[Z r+abs(y)*i];
这一行是什么意思？

赵双双 · 发表于 2013-10-4 09:23

自治系统用哪些方法

echo-echo · 发表于 2014-5-21 18:06

个人认为，最大值方法不靠谱，并不是一个周期就对应用一个最值的

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

[分形与混沌] 【总结帖】分岔图绘制不同方法的总结、比较