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
真是强人得好好向各位学习啊{:{26}:}
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 = 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 = 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
{:{23}:}最大值法和Poincare截面法哪个方法效果好一些?哪个迭代更快些
mybook86
发表于 2012-6-13 16:05
{:{23}:}
如枫
发表于 2013-1-28 22:28
请教一下,
Z=;
这一行是什么意思?
赵双双
发表于 2013-10-4 09:23
自治系统用哪些方法
echo-echo
发表于 2014-5-21 18:06
个人认为,最大值方法不靠谱,并不是一个周期就对应用一个最值的