hettesthenon.asv

计算动力学系统的分岔图

x1=[-0.31622776601684;-0.31622776601684];x2=[-0.31560902084518;-0.31539331317194];
x3=[-0.31539355715830;-0.31508178332301];x4=[-0.31508467959730; -0.31466796144265];    
x5=[-0.31466813859739; -0.31410047552221]; x6=[-0.31410046307950; -0.31332725645617];  
x7=[-0.31332723597227;-0.31227444777060];x8=[-0.31227448507563;-0.31084253925511];
x9=[-0.31084264307434;-0.30889749276922];x10=[-0.30889771353950;-0.30626014122335];
x11=[-0.30626046517850;-0.30269274794841]; x12=[-0.30269317475526; -0.29788304781852];
x13=[-0.29788367908864;-0.29142701352582];x14=[-0.29142818352847;-0.28281250314794];
x15=[-0.28281516993790;-0.27140906914917]; x16=[-0.27141565244602;-0.25647298255308];
x17=[-0.25648860532045;-0.23718191017136]; x18=[-0.23720688901728;-0.21272406111122];
x19=[-0.21275866615506;-0.18243024189220]; x20=[-0.18248039606343;-0.14599724084430];
x21=[-0.14607278783116;-0.10372937367153]; x22=[-0.10384307676712;-0.05672813648452];
x23=[-0.05689064603087;-0.00690403479402]; x24= [-0.00711578678428;0.04327756910813];   
x25=[0.04303188174671; 0.09127496361279];  x26=[0.09102308870982; 0.13498411331315];   
x27=[0.13475485581872; 0.17307101848360];  x28=[0.17288318591903; 0.20502897288831];   
x29=[0.20488782105962;0.23102215898824];   x30=[0.23092288757316; 0.25164295697730];
x31=[0.25157637798561;0.26768489795363];  x32=[0.26764275584108; 0.27997672537471];
x33=[0.27995211355481;0.28928851172323];  x34=[0.28927483876242; 0.29628256787354];
x35=[0.29627498723939;0.30150167316507];  x36=[0.30149744431039; 0.30537743637787];
x37=[0.30537504839377;0.30824529794061]; x38=[0.30824393010609; 0.31036174159792];
x39=[0.31036097611244; 0.31192056006836]; x40=[0.31192016442221;0.31306707871867];
x41=[0.31306687092959; 0.31390945162499]; x42=[0.31390934025118; 0.31452787887399];
x43=[0.31452781777021; 0.31498163696176]; x44=[0.31498160261154;0.31531443299027];
x45=[0.31531441398504; 0.31555843688580]; x46=[0.31555842716090; 0.31573730016370];
x47=[0.31573729509531; 0.31586839136228]; x48=[0.31586838866002; 0.31596445812121];
x49=[0.31596445664381; 0.31603485188048]; x50=[0.31603485105268; 0.31608643017051];
x51=[0.316086429714;  0.316124220337];  x52=[0.31622776601684;0.31622776601684];

C(:,1)=x1;C(:,2)=x2;C(:,3)=x3;C(:,4)=x4;,C(:,5)=x5;C(:,6)=x6;C(:,7)=x7;C(:,8)=x8;
C(:,9)=x9;C(:,10)=x10;C(:,11)=x11;C(:,12)=x12;,C(:,13)=x13;C(:,14)=x14;C(:,15)=x15;C(:,16)=x16;
C(:,17)=x17;C(:,18)=x18;C(:,19)=x19;C(:,20)=x20;,C(:,21)=x21;C(:,22)=x22;C(:,23)=x23;C(:,24)=x24;
C(:,25)=x25;C(:,26)=x26;C(:,27)=x27;C(:,28)=x28;,C(:,29)=x29;C(:,30)=x30;C(:,31)=x31;C(:,32)=x32;
C(:,33)=x33;C(:,34)=x34;C(:,35)=x35;C(:,36)=x36;,C(:,37)=x37;C(:,38)=x38;C(:,39)=x39;C(:,40)=x40;
C(:,41)=x41;C(:,42)=x42;C(:,43)=x43;C(:,44)=x44;,C(:,45)=x45;C(:,46)=x46;C(:,47)=x47;C(:,48)=x48;
C(:,49)=x49;C(:,50)=x50;C(:,51)=x51;C(:,52)=x52;
ap=[2];%p=[a,b,r]
 
 disp('>> opt = contset; ');
 opt = contset;p=[0.1;-1];
disp('>> [x0,v0]=init_HE_HE(@Henon,C, p, ap,J,eps0,eps1);');%(mapfile, C, p,ap,J, eps0,eps1
opt=contset(opt,'MaxNumpoints',);
opt=contset(opt,'Singularities',1);
opt=contset(opt,'Backward',0);
opt = contset(opt,'AutDerivative',1);
opt=contset(opt,'FunTolerance',1e-3);%-1
opt=contset(opt,'VarTolerance',2e-2);
opt=contset(opt,'MinStepsize',1e-6);
[x0,v0]=init_HE_HE(@SHmap,C, p, ap,1);
[xhet,vhet,shet,hhet,fhet]=cont(@heteroclinic,x0,[],opt);
grid,
for i=1:52
    cpl(xhet,vhet,shet,[107 2*(i-1)+1])
    
end