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

飞翼

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

ChaChing

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

飞翼

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

                               
登录/注册后可看大图

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

function fga
y0=[-0.05;0;0];
options=odeset;options.RelTol=1e-7;options.AbsTol=1e-8;
ts=[0:0.04*pi:3000];
[t,x]=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=[x(2);-2*u*x(2)-(1+108/35*g^2)*x(1)-108/35*g*x(1)^2-36/35*x(1)^3-g-36/35*g^3+B*cos(x(3));w];
else
   dx=[x(2);-2*u*x(2)-(1+108/35*g^2)*x(1)-108/35*g*x(1)^2-36/35*x(1)^3-g-36/35*g^3+B*cos(x(3))-a*x(1);w];
end

function fga
y0=[-0.05;0;0];
options=odeset;options.RelTol=1e-7;options.AbsTol=1e-8;
ts=[0:0.04*pi:3000];
[t,x]=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=[x(2);-2*u*x(2)-(1+108/35*g^2)*x(1)-108/35*g*x(1)^2-36/35*x(1)^3-g-36/35*g^3+B*cos(x(3));w];
else
   dx=[x(2);-2*u*x(2)-(1+108/35*g^2)*x(1)-108/35*g*x(1)^2-36/35*x(1)^3-g-36/35*g^3+B*cos(x(3))-a*x(1);w];
end

飞翼

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

                               
登录/注册后可看大图

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

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

ChaChing

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

                               
登录/注册后可看大图

就像这样，g分别取这两个数时，都这个相图，是不是最后两位精度不够显示不出来？不知道怎么设置，正常来说 ...

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

飞翼

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

                               
登录/注册后可看大图

matlab变数是双精确度, 不致於有"最后两位精度不够"的问题!
试下aaa=0.202531-0.2025不就知道!?
不过不 ...

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

飞翼

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

                               
登录/注册后可看大图

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

⊙﹏⊙b汗，各种贴吧论坛上的少……