lesliegower34.m

计算动力学系统的分岔图

disp('>> global  x15 v15 s15  opt cds fpmds')
global  x15 v15 s15  opt cds fpmds
disp('>> opt = contset; ');
opt = contset;
disp('>> opt = contset(opt,''Multipliers'',1); ');
opt = contset(opt,'Multipliers',1);
disp('>> global opt cds fpmds');
global opt cds fpds 
disp('>> opt = contset; ');
opt = contset;
disp('>> opt = contset(opt,''Multipliers'',1); ');
opt = contset(opt,'Multipliers',1);
disp('>> ap=11; p=[20;18;0.72;0.23;0.29;0.98;0.36;0.55;0.18;0.26;0;0.23;0;.08];');
p = [20;18;0.72;0.23;0.29;0.98;0.36;0.55;0.18;0.26;0;0.23;0;.08];ap=11;
disp('>> [x0,vO]=init_FPm_FPm(@@LeslieGower,[0;0;32.68698060;23.68138391], p, ap,1);');
[x0,vO]=init_FPm_FPm(@LeslieGower,[0;0;32.68698060;23.68138391], p, ap,1);
disp('>> opt=contset;opt=contset(opt,''MaxNumPoints'',200);');
opt=contset(opt,'MaxNumPoints',200);
disp('>> opt=contset(opt,''Singularities'',1);');
opt=contset(opt,'Singularities',1);
disp('>> [x15,v15,s15,h15,f15]=cont(@fixedpointmap,x0,[],opt);');
[x15,v15,s15,h15,f15]=cont(@fixedPointmap,x0,[],opt,1);

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

clc
disp('>>>>>>> flip curve <<<<<<<');
disp('>> global x15 v15 s15 opt cds fpmds pdmds');
global x15 v15 s15 opt cds fpmds pdmds
disp('>> p1(fpmds.ActiveParams)=x15(end,s15(3).index);');
p1=p;
p1(fpmds.ActiveParams)=x15(end,s15(3).index);
disp('>> x1=x15(1:4,s15(3).index);');
x1=x15(1:4,s15(3).index);
disp('>> opt=contset(opt,''Maxnumpoints'',40);');
opt=contset(opt,'Maxnumpoints',40);
opt = contset(opt,'Multipliers',1);
disp('>> [x0,v0]=init_PDm_PDm(@LeslieGower,x1,p1,[11 13],1);');
[x0,v0]=init_PDm_PDm(@LeslieGower,x1,p1,[11 13],1);
disp('>> opt=contset(opt,''Backward'',0);');

opt=contset(opt,'Backward',0);
disp('>> [x341,v341,s341,h341,f341]=cont(@perioddoublingmap,x0,v0,opt);');

[x341,v341,s341,h341,f341]=cont(@perioddoublingmap,x0,v0,opt);
disp('>> cpl(x341,v341,s341,[6 5]);');

cpl(x341,v341,s341,[6 5]);
disp('>> opt=contset(opt,''Backward'',1);');

opt=contset(opt,'Backward',1);
disp('>> opt=contset(opt,''Maxnumpoints'',20); ');

opt=contset(opt,'Maxnumpoints',20); 
%opt=contset(opt,'IgnoreSingularity',[2]);
disp('>> [x342,v342,s342,h342,f342]=cont(@perioddoublingmap,x0,v0,opt);');

[x342,v342,s342,h342,f342]=cont(@perioddoublingmap,x0,v0,opt);
disp('cpl(x342,v342,s342,[6 5]);')
cpl(x342,v342,s342,[6 5]);