hettestrun.asv

计算动力学系统的分岔图

global cds hetds
%C=[0.53604733513342,  0.5155,  .5354,  .5021,  .5573,  .4643,  .6167,  .3560,  .7689,  .0380,  1.0657,  -.8020,  .8098,  -1.0941,  .0027,  -.7955,   -.3299,   -.6318 ,  -.4505,   -.5641,  -.4949,  -.5377, -.5114, -.5280,  -.5178,   -0.50877460786069 ;0.53604733513342,.5354,.5021,.5572,.4641,.6167,.3560 ,.7689,.0380, 1.0657, -.8019 ,.8098,-1.0941,.0027,-0.7955,-.3299,-.6318, -.4505 ,-.5641 ,-.4950, -.5379, -.5113, -.5275, -.5178,-.5238,-0.50877460786069];
 %p=[0.3;-1.03;-0.1];ap=[1];%p=[a,b,r]
 %C=[0.46661702380495, 0.52096866,0.48858445,0.483630993,0.6123,0.8412,1.229,0.642,0.1348,-0.1434,-0.28495812,-0.35634007,-0.3929237, -0.42861702380495; 0.46661702380495, 0.375056655,0.43497755,0.4598092143,0.2067,-0.2762,-1.333,-1.093,-0.7998,-0.6233,-0.52893911,-0.480307989,-0.44762207,-0.42861702380495];
 %C=[0.46661702380495, 0.52096866,0.48858445,0.483630993,0.6123,0.841195,1.22990435,0.6419816,0.1347307,-0.14345697,-0.28495812,-0.35634007,-0.3929237, -0.42861702380495; 0.46661702380495, 0.375056655,0.43497755,0.4598092143,0.2067,-0.276064,-1.3327006,-1.0930203,-0.7998431,-0.6233856,-0.52893911,-0.480307989,-0.44762207,-0.42861702380495];
 %C=[0.46661702380495,0.6123,0.841195,1.22990435,0.6419816,0.1347307,-0.14345697, -0.28495812,-0.35634007,-0.3929237,-0.42861702380495; 0.46661702380495, 0.2067,-0.276064,-1.3327006,-1.0930203,-0.7998431,-0.6233856,-0.52893911,-0.480307989,-0.44762207,-0.42861702380495];
 C=[0.46661702380495,0.51950188769804 ,0.48317182408537 ,0.47311690109635 ,0.6123,0.841195,1.22990435,0.6419816,0.1347307,-0.14345697, -0.33379938281268,-0.38062126993587,-0.40424194997464,-0.41621339121235,-0.42230324914752,-0.42861702380495;0.46661702380495,0.37639405084666,0.43915841488454, 0.45696220552142, 0.2067,-0.276064,-1.3327006,-1.0930203,-0.7998431,-0.6233856,-0.49516248970159,-0.46255049612568,-0.44591635361118,-0.43743735402657,-0.43311132511221,-0.42861702380495];
 %C=[0.46661702380495,0.51950188769804 ,0.48317182408537 ,0.47311690109635 ,0.6123,0.841195,1.22990435,0.6419816,0.1347307,-0.14345697, -0.28481547435852,-0.35598743914033,-0.39177671013039 ,-0.41000328504934 ,-0.41912366327007, -0.42861702380495;0.46661702380495,0.37639405084666,0.43915841488454, 0.45696220552142, 0.2067,-0.276064,-1.3327006,-1.0930203,-0.7998431,-0.6233856,-0.52879374519471,-0.47976672178178,-0.45471025132465,-0.44183990874129,-0.43537107690723,-0.42861702380495];

%C=[0.440593, 0.448949, 0.462354, 0.497145 ,0.585534, 0.794673, 1.143674,  0.794261,0.278157,-0.041412,  -0.213356,-0.303417,-0.350506,  -0.375192, -0.388171,-0.402593 ; 0.440593, 0.427069, 0.404981, 0.346314 ,0.187883 ,-0.241496,-1.159085,-1.129627,-0.848210 ,-0.649767,-0.535571 , 0.473342, -0.440062, -0.422400,-0.413066, -0.402593];
%   load Cit1, 
%   xx=xxhet1(:,350);
%   p=[xx(35);-1.057;-0.5;0];%p=[alpha,beta,R,S]
%    C=[xx(1),xx(3),xx(5),xx(7),xx(9),xx(11),xx(13),xx(15),xx(17),xx(19),xx(21),xx(23),xx(25),xx(27),xx(29),xx(31); xx(2),xx(4),xx(6),xx(8),xx(10),xx(12),xx(14),xx(16),xx(18),xx(20),xx(22),xx(24),xx(26),xx(28),xx(30),xx(32)];
p=[0.3;-1.057;-0.5;0];
 ap=[];%size(C),pause
 %
 opt = contset;
 %opt=contset(opt,'MinStepsize',0.0001);
 %opt=contset(opt,'MaxStepsize',0.05);
 opt=contset(opt,'FunTolerance',1e-6);%0.1
 opt=contset(opt,'VarTolerance',1e-6);%0.01
 opt=contset(opt,'TestTolerance',1e-4);
 

disp('>> [x0,v0]=init_HE_HE(@Ghmap,C, p, ap,2));');%(mapfile, C, p,ap,J, eps0,eps1)
[x0,v0]=init_HE_HE(@Ghmap,C, p, ap,2);

opt=contset(opt,'MaxNumpoints',200);
opt=contset(opt,'Singularities',1);
opt=contset(opt,'Backward',1);
opt = contset(opt,'AutDerivative',0);

opt=contset(opt,'MinStepsize',1);
[xhetr,vhet,shet,hhet,fhet]=cont(@heteroclinic,x0,[],opt);
% for i=1:16
     %cpl(xhet,vhet,shet,[35 2*(i-1)+1])
     cpl(xhetr,vhet,shet,[35 1])
 %end