lesliegower31.m

来自「计算动力学系统的分岔图」· M 代码 · 共 38 行

disp('>> global opt cds fpmds pdmds ')
global opt cds fpmds pdmds
LeslieGower12;
disp('>> opt = contset; ');
opt = contset;
disp('>> opt = contset(opt,''Multipliers'',1); ');
opt = contset(opt,'Multipliers',1);
opt=contset(opt,'Singularities',1);
clc
disp('global x12 v12 s12 opt cds fpmds pdmds')
global x12 v12 s12 opt cds fpmds pdmds
disp('>>>>>>> flip curve <<<<<<<');
disp('p1(fpmds.ActiveParams)=x12(end,s12(4).index);');
p1(fpmds.ActiveParams)=x12(end,s12(4).index);
disp('x1=x12(1:4,s12(4).index);');
x1=x12(1:4,s12(4).index);
disp('opt=contset(opt,''Maxnumpoints'',15);');
opt=contset(opt,'Maxnumpoints',15);
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('[x1,v1,s1,h1,f1]=cont(@perioddoublingmap,x0,v0,opt);');
[x1,v1,s1,h1,f1]=cont(@perioddoublingmap,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);
disp('opt=contset(opt,''Maxnumpoints'',20);');
opt=contset(opt,'Maxnumpoints',20); %260
disp('[x2,v2,s2,h2,f2]=cont(@perioddoublingmap,x0,v0,opt);');
[x2,v2,s2,h2,f2]=cont(@perioddoublingmap,x0,v0,opt);
disp('cpl(x2,v2,s2,[6 5]);');
cpl(x2,v2,s2,[6 5]);