thenon4.m

来自「计算动力学系统的分岔图」· M 代码 · 共 28 行

global cds hetTds hetds xhet vhet shet

p=[0.3;-1.057;-0.5;0];
ap=[1 2];
x=xhet(1:end-1,shet(3).index);
p2=p;p2(hetds.ActiveParams)=xhet(end,shet(3).index);
nu=hetds.nu;
ns=hetds.ns;
nphase=hetds.nphase;

[x0,v0]=init_HetT_HetT(@Ghmap,x,nphase,nu,ns, p2, ap,2);

opt=contset(opt,'MaxNumpoints',200);
opt=contset(opt,'Singularities',0);
opt=contset(opt,'Backward',0);
% opt=contset(opt,'FunTolerance',1e-8);
% opt=contset(opt,'VarTolerance',1e-8);
 opt=contset(opt,'MaxStepsize',1e1);
 opt=contset(opt,'MinStepsize',1e-5);
%MaxStepsize
opt = contset(opt,'AutDerivative',1);

[xhetT,vhetT,shetT,hetT,fhetT]=cont(@heteroclinicT,x0,[],opt);
 cpl(xhetT,vhetT,shetT,[36 35]), pause
 opt=contset(opt,'Backward',1);
 [xhetT,vhetT,shetT,hetT,fhetT]=cont(@heteroclinicT,x0,[],opt);
 cpl(xhetT,vhetT,shetT,[36 35])