飞翼 发表于 2013-7-11 11:24

分岔图分岔参数的有效数字问题

在做分岔图时，分岔参数g=0.202531时做的图与g=0.2025的图一样，也就是后两位数字显示不出来，请高手指点

ChaChing 发表于 2013-7-11 20:02

不清楚"g=0.202531时做的图与g=0.2025的图一样"的意思!?

飞翼 发表于 2013-7-11 20:54

ChaChing 发表于 2013-7-11 20:02 static/image/common/back.gif
不清楚"g=0.202531时做的图与g=0.2025的图一样"的意思!?

function fga
y0=[-0.05;0;0];
options=odeset;options.RelTol=1e-7;options.AbsTol=1e-8;
ts=;
=ode45(@dsg,ts,y0,options);
N1=length(x(:,1))-10000;
N2=length(x(:,1));
YYY1=x(N1:N2,1)
YYY2=x(N1:N2,2)
T=t(N1:N2)
figure(1)
plot(YYY1,YYY2)
function dx=dsg(t,x)
u=0.01;
a=20;
B=0.2;
w=0.3;
g=0.2025;
if x(1)<0
dx=;
else
   dx=;
end

function fga
y0=[-0.05;0;0];
options=odeset;options.RelTol=1e-7;options.AbsTol=1e-8;
ts=;
=ode45(@dsg,ts,y0,options);
N1=length(x(:,1))-10000;
N2=length(x(:,1));
YYY1=x(N1:N2,1)
YYY2=x(N1:N2,2)
T=t(N1:N2)
figure(1)
plot(YYY1,YYY2)
function dx=dsg(t,x)
u=0.01;
a=20;
B=0.2;
w=0.3;
g=0.202531;
if x(1)<0
dx=;
else
   dx=;
end

飞翼 发表于 2013-7-11 21:02

ChaChing 发表于 2013-7-11 20:02 static/image/common/back.gif
不清楚"g=0.202531时做的图与g=0.2025的图一样"的意思!?

就像这样，g分别取这两个数时，都这个相图，是不是最后两位精度不够显示不出来？不知道怎么设置，正常来说g=0.202531时，相图应该是稳定的1周期，我是看别人的文章自己仿的，但是g=0.202531时显示不出1周期

ChaChing 发表于 2013-7-11 23:42

飞翼 发表于 2013-7-11 21:02 static/image/common/back.gif
就像这样，g分别取这两个数时，都这个相图，是不是最后两位精度不够显示不出来？不知道怎么设置，正常来说 ...

matlab变数是双精确度, 不致於有"最后两位精度不够"的问题!
试下aaa=0.202531-0.2025不就知道!?
不过不清楚不到0.02%(31/202500)的差异, LZ的结果会有影响吗?

飞翼 发表于 2013-7-12 10:32

ChaChing 发表于 2013-7-11 23:42 static/image/common/back.gif
matlab变数是双精确度, 不致於有"最后两位精度不够"的问题!
试下aaa=0.202531-0.2025不就知道!?
不过不 ...

目前正在初学，还没做Lz图，g取0.2025到0.202531做了次分岔图，结果都是混沌现象，也许是那人的图本身就是错误的

飞翼 发表于 2013-7-12 14:20

飞翼 发表于 2013-7-12 10:32 static/image/common/back.gif
目前正在初学，还没做Lz图，g取0.2025到0.202531做了次分岔图，结果都是混沌现象，也许是那人的图本身就是 ...

⊙﹏⊙b汗，各种贴吧论坛上的少……
页: [1]
查看完整版本: 分岔图分岔参数的有效数字问题