lesliegower28.m

计算动力学系统的分岔图

disp('>> global  x221 v221 s221  opt cds fpmds  lpmds')
global  x24 v24 s24  opt cds fpmds  lpmds

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;

[x0,vO]=init_FPm_FPm(@LeslieGower,[11.17052012;11.94625068;20.99523834;15.21083594], p, ap,1);

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

[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,'Backward',0);
x1=x221(1:4,s221(2).index);p1=p;p1(fpmds.ActiveParams)=x221(5,s221(2).index);
opt=contset(opt,'MaxNumPoints',1300);
opt=contset(opt,'Backward',1);
opt=contset(opt,'Singularities',1);
[x2,v2]=init_PDm_FP2m(@LeslieGower,x1,p1,s221(2),0.01,1);
[x24,v24,s24,h24,f24]=cont(@fixedPointmap,x2,[],opt);
clc
disp('>>>>>>>>>>>>>>fold curve <<<<<<<<<<<');
disp('>> p1(fpmds.ActiveParams)=x24(end,s24(3).index);');
p1(fpmds.ActiveParams)=x24(end,s24(3).index);
disp('>>x1=x24(1:4,s24(3).index);');
x1=x24(1:4,s24(3).index);
disp('>> opt=contset(opt,''Maxnumpoints'',90);');
opt=contset(opt,'Maxnumpoints',90);
disp('>> [x0,v0]=init_LPm_LPm(@LeslieGower,x1,p1,[11 13],2); ');
[x0,v0]=init_LPm_LPm(@LeslieGower,x1,p1,[11 13],2);
disp('>>opt=contset(opt,''Backward'',0); ');
opt=contset(opt,'Backward',0);
disp('>> [x1,v1,s1,h1,f1]=cont(@limitpointmap,x0,v0,opt); ');
[x1,v1,s1,h1,f1]=cont(@limitpointmap,x0,v0,opt);
disp('>> cpl(x1,v1,s1,[6 5]); ');
cpl(x1,v1,s1,[6 5]);
disp('>> opt=contset(opt,''Backward'',1); ');
opt=contset(opt,'Backward',1);
%opt=contset(opt,'IgnoreSingularity',[2]);
disp('>>[x2,v2,s2,h2,f2]=cont(@limitpointmap,x0,v0,opt); ');
[x2,v2,s2,h2,f2]=cont(@limitpointmap,x0,v0,opt);
disp('>> cpl(x2,v2,s2,[6 5]); ');
cpl(x2,v2,s2,[6 5]);