lesliegower30.m

计算动力学系统的分岔图

disp('global x222 v222 s222opt cds fpmds pdmds')
global x222 v222 s222  opt cds fpmds pdmds
disp('>> opt = contset; ');
opt = contset;
disp('>> opt = contset(opt,''Multipliers'',1); ');
opt = contset(opt,'Multipliers',1);
p = [20;18;0.72;0.23;0.29;0.98;0.36;0.55;0.18;0.26;0.3;0.23;0.3;.08];ap=11;
x=[11.17052012;11.94625068;20.99523834;15.21083594];
[x0,vO]=init_FPm_FPm(@LeslieGower,x, p, ap,1);
opt = contset;
opt=contset;opt=contset(opt,'MaxNumPoints',900);

opt=contset(opt,'Singularities',1);

disp('>>>>>>> FIXEDPOINTS <<<<<<<');

[x221,v221,s221,h221,f221]=cont(@fixedPointmap,x0,[],opt,1);

opt=contset(opt,'Backward',1);

[x222,v222,s222,h222,f222]=cont(@fixedPointmap,x0,[],opt,1);

opt=contset(opt,'Maxnumpoints',250);

ap1=11; p1=[20;18;0.72;0.23;0.29;0.98;0.36;0.55;0.18;0.26;0.3;0.23;0.3;.08];


p1(ap1)=x221(end,s221(2).index);x1=x221(1:4,s221(2).index);

clc
opt=contset(opt,'Maxnumpoints',290);
disp('>>>>>>> PERIODDOUBLING <<<<<<<');

[x0,v0]=init_PDm_PDm(@LeslieGower,x1,p1,[11 13],1);

opt=contset(opt,'Backward',0);

[x1,v1,s1,h1,f1]=cont(@perioddoublingmap,x0,v0,opt);

cpl(x1,v1,s1,[6 5]);


opt=contset(opt,'Backward',1);
opt=contset(opt,'Maxnumpoints',290); %260
%opt=contset(opt,'IgnoreSingularity',[2]);

[x2,v2,s2,h2,f2]=cont(@perioddoublingmap,x0,v0,opt);

cpl(x2,v2,s2,[6 5]);

  
clc
disp('global x12 v12 s12 opt cds fpmds pdmds')
global opt cds fpmds pdmds
disp('>>>>>>> SWITCHING <<<<<<<'); 
disp('opt=contset(opt,''Maxnumpoints'',250); '); 


opt=contset(opt,'Maxnumpoints',200); %300
disp('ap = [11 13];y1=x1(1:4,s1(2).index);p1(ap)=x1([5 6],s1(2).index); '); 


ap = [11 13];y1=x1(1:4,s1(2).index);p1(ap)=x1([5 6],s1(2).index); 
disp('[x0,v0]=init_GPD_LP2m(@LeslieGower,1e-3,y1,p1,[11 13],1);'); 


[x0,v0]=init_GPD_LP2m(@LeslieGower,1e-3,y1,p1,[11 13],1);
disp('opt=contset(opt,''Backward'',1);'); 


opt=contset(opt,'Backward',1);
disp('[x7,v7,s7,h7,f7]=cont(@limitpointmap,x0,[],opt);'); 

[x7,v7,s7,h7,f7]=cont(@limitpointmap,x0,[],opt);
disp('cpl(x7,v7,s7,[6 5]);'); 


cpl(x7,v7,s7,[6 5]);
disp('y1=x2(1:4,s2(2).index);p1(ap)=x2([5 6],s2(2).index);'); 


y1=x2(1:4,s2(2).index);p1(ap)=x2([5 6],s2(2).index);
disp('opt=contset(opt,''Maxnumpoints'',260);'); 

opt=contset(opt,'Maxnumpoints',250);
disp('[x0,v0]=init_GPD_LP2m(@LeslieGower,1e-1,y1,p1,[11 13],1);'); 

[x0,v0]=init_GPD_LP2m(@LeslieGower,1e-1,y1,p1,[11 13],1);
disp('[x8,v8,s8,h8,f8]=cont(@limitpointmap,x0,[],opt);'); 

opt=contset(opt,'Backward',0);
[x8,v8,s8,h8,f8]=cont(@limitpointmap,x0,[],opt);
disp('cpl(x8,v8,s8,[6 5]);'); 


cpl(x8,v8,s8,[6 5]);