?? homtestrun.m
字號:
% C=[-1.69449184777905 -1.6837643 -1.65949674 -1.581150109 -1.337009429 -0.659661307 0.547365455 0.062065530 -1.009338894 -1.480864536 -1.632164287;-1.69449184777905 -1.65948946 -1.581161931 -1.337017143 -0.659669162 0.547301706 0.062074908 -1.009383240 -1.480879089 -1.632160728 -1.676628792];
% p=[-0.4;1.1;-0.1];ap=[2];%p=[a,b,r] *perfect
% C,pause
% % C=[-0.3144706727040 -1.04920163606 -0.640540809357 0.392687493818 0.324301394255 -0.718447548185 -1.02365494748 -0.588597191062 0.437429350429 0.234771666962 -0.711565274823 -0.929086704063 -0.403516106938 0.547353044603 -0.0659223493914 -0.759037581779 -0.713901008241 0.010375914355 0.528812106379 -0.490774367408 -0.954295085512 -0.656923612169 0.279142457533 0.410479099965 -0.664676612665 -1.03320009562 -0.658206686637 0.352622134037 0.370238275605 -0.709806378525 -1.05518838715 -0.6643179710 ;0.755102259927 -0.640540809357 0.392687493818 0.324301394255 -0.718447548185 -1.02365494748 -0.588597191062 0.437429350429 0.234771666962 -0.711565274823 -0.929086704063 -0.403516106938 0.547353044603 -0.0659223493914 -0.759037581779 -0.713901008241 0.010375914355 0.528812106379 -0.490774367408 -0.954295085512 -0.656923612169 0.279142457533 0.410479099965 -0.664676612665 -1.03320009562 -0.658206686637 0.352622134037 0.370238275605 -0.709806378525 -1.05518838715 -0.66431797106 0.364875727378 ];
% % C=[-0.3144706727040 -1.04920163606 -1.02365494748 -0.929086704063 -0.713901008241 0.279142457533 0.352622134037 ;0.755102259927 -0.640540809357 -0.588597191062 -0.403516106938 0.010375914355 0.410479099965 0.370238275605 ];
%
%
% disp('>> opt = contset; ');
% opt = contset;
% disp('>> [x0,v0]=init_HO_HO(@Henon,C, p, ap,J);');
% [x0,v0]=init_HO_HO(@Henon,C, p, ap,1);
% opt=contset; opt=contset(opt,'MaxNumpoints',41);
% opt=contset(opt,'Singularities',1);
% %opt=contset(opt,'MinStepsize',0.001);
% opt=contset(opt,'Backward',1);
% opt = contset(opt,'AutDerivative',1);
% [xhom,vhom,shom,hhom,fhom]=cont(@homoclinic,x0,[],opt);
%
% for i=1:1
% %cpl(xhom,vhom,shom,[25 2*(i-1)+1])
% cpl(xhom,vhom,shom,[25 11])
% end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%46 NOV
C=[-1.62114638486435 -1.562 -1.443 -1.0987856 -0.27959455684736 0.62285460449110 -0.48255078907068 -1.24433960857123 -1.51139945471518 -1.59072793547249;-1.62114638486435 -1.443 -1.0987856 -0.27959455684736 0.62285460449110 -0.48255078907068 -1.24433960857123 -1.51139945471518 -1.59072793547249 -1.61409645976328];
p=[-0.4;1.03;-0.1];ap=[2];%p=[a,b,r]
disp('>> opt = contset; ');
opt = contset;
disp('>> [x0,v0]=init_HO_HO(@Henon,C, p, ap,J);');
[x0,v0]=init_HO_HO(@Henon,C, p, ap,1);
opt=contset; opt=contset(opt,'MaxNumpoints',35);
opt=contset(opt,'Singularities',1);
%opt=contset(opt,'MinStepsize',0.001);
opt=contset(opt,'Backward',1);
opt = contset(opt,'AutDerivative',1);
[xhom,vhom,shom,hhom,fhom]=cont(@homoclinic,x0,[],opt);
for i=1:1
%cpl(xhom,vhom,shom,[25 2*(i-1)+1])
cpl(xhom,vhom,shom,[23 11])
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -