thenon8.asv

计算动力学系统的分岔图

global cds homTds homds xhom vhom shom

p=[-0.4;1.03;-0.1];
ap=[1 2];
x=xhom(1:end-1,shom(2).index);
p2=p;p2(homds.ActiveParams)=xhom(end,shom(2).index);
nu=homds.nu;
ns=homds.ns;
nphase=homds.nphase;
[x0,v0]=init_homT_homT(@Henon,x,nphase,nu,ns, p2, ap,1);
opt=contset(opt,'MaxNumpoints',53);
opt=contset(opt,'Singularities',0);
opt=contset(opt,'Backward',0);
%opt = contset(opt,'AutDerivative',1);
[xhomT,vhomT,shomT,homT,fhomT]=cont(@homoclinicT,x0,[],opt);
 cpl(xhomT,vhomT,shomT,[24 23])