I have a system with 11 states. I use Cl_matcont5p4 to define my system and perform continuation from an equilibrium point. I get a series of limit points, with the active parameter almost unvarying and then it stops after some steps saying "closed curve detected". Basically, at every value of the active parameter, I get a fixed point.
This is the code I used to continue:
X0 = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]; % state corresponding to an equilibrium
p = [1.2977]; % parameter corresponding to initial equilibrium state
ap = [1]; % I just use one active parameter [x0,v0]=init_EP_EP(@my_system, X0', p, ap);
opt=contset;
opt=contset(opt,'MaxNumPoints',600);
opt=contset(opt,'Singularities',1); [x,v,s,h,f]=cont(@equilibrium,x0,[],opt);
It obviously depends on the system. But could there be a rookie matcont user error that is leading to this?
Last edit: Anonymous 2019-07-13
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
I get a series of limit points, with the active parameter almost unvarying and then it stops after some steps saying "closed curve detected"
After 50 steps (default option) it starts to check whether you are back at the original point. If the points are really similar, then this happens. Set the value to 0 to have MatCont ignore this check. All other things indeed depend on your system very much.
Strangely, the fixed points plotted, correspond to the time histories of the states.
That's a strange statement. How would you get "history"?
Basically, at every value of the active parameter, I get a fixed point.
Of course, by definition, as you continue equilibria.
From: Rohini V rohini-v@users.sourceforge.net
Sent: Thursday, July 11, 2019 6:48 AM
To: [matcont:discussion]
Subject: [matcont:discussion] Matcont finds a limit point at every point n the curve.
I have a system with 11 states. I use Cl_matcont5p4 to define my system and perform continuation from an equilibrium point. I get a series of limit points, with the active parameter almost unvarying and then it stops after some steps saying "closed curve detected". Strangely, the fixed points plotted, correspond to the time histories of the states. Basically, at every value of the active parameter, I get a fixed point.
This is the code I used to continue:
X0 = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]; % state corresponding to an equilibrium
p = [1.2977]; % parameter corresponding to initial equilibrium state
ap = [1]; % I just use one active parameter [x0,v0]=init_EP_EP(@my_system, X0', p, ap);
opt=contset;
opt=contset(opt,'MaxNumPoints',600);
opt=contset(opt,'Singularities',1); [x,v,s,h,f]=cont(@equilibrium,x0,[],opt);
It obviously depends on the system. But could there be a rookie matcont user error that is leading to this?
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Anonymous
-
2019-07-11
That's a strange statement. How would you get "history"?
By time history I mean time integration of the states. It results in a particular periodic orbit and the limit points seem to follow the same orbit.
After 50 steps (default option) it starts to check whether you are back at the original point. If the >>points are really similar, then this happens. Set the value to 0 to have MatCont ignore this check. >>All other things indeed depend on your system very much.
Thanks! I will try that out.
Last edit: Anonymous 2019-07-11
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
That's a strange statement. How would you get "history"?
By time history I mean time integration of the states. It results in a particular periodic orbit and the limit points seem to follow the same orbit.
That suggests your system has some additional features... MatCont focusses on generic dissipative systems. If your system has more structure you should take that into account yourself.
I have a system with 11 states. I use Cl_matcont5p4 to define my system and perform continuation from an equilibrium point. I get a series of limit points, with the active parameter almost unvarying and then it stops after some steps saying "closed curve detected". Basically, at every value of the active parameter, I get a fixed point.
This is the code I used to continue:
X0 = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]; % state corresponding to an equilibrium
p = [1.2977]; % parameter corresponding to initial equilibrium state
ap = [1]; % I just use one active parameter
[x0,v0]=init_EP_EP(@my_system, X0', p, ap);
opt=contset;
opt=contset(opt,'MaxNumPoints',600);
opt=contset(opt,'Singularities',1);
[x,v,s,h,f]=cont(@equilibrium,x0,[],opt);
It obviously depends on the system. But could there be a rookie matcont user error that is leading to this?
Last edit: Anonymous 2019-07-13
After 50 steps (default option) it starts to check whether you are back at the original point. If the points are really similar, then this happens. Set the value to 0 to have MatCont ignore this check. All other things indeed depend on your system very much.
That's a strange statement. How would you get "history"?
Of course, by definition, as you continue equilibria.
From: Rohini V rohini-v@users.sourceforge.net
Sent: Thursday, July 11, 2019 6:48 AM
To: [matcont:discussion]
Subject: [matcont:discussion] Matcont finds a limit point at every point n the curve.
I have a system with 11 states. I use Cl_matcont5p4 to define my system and perform continuation from an equilibrium point. I get a series of limit points, with the active parameter almost unvarying and then it stops after some steps saying "closed curve detected". Strangely, the fixed points plotted, correspond to the time histories of the states. Basically, at every value of the active parameter, I get a fixed point.
This is the code I used to continue:
X0 = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]; % state corresponding to an equilibrium
p = [1.2977]; % parameter corresponding to initial equilibrium state
ap = [1]; % I just use one active parameter
[x0,v0]=init_EP_EP(@my_system, X0', p, ap);
opt=contset;
opt=contset(opt,'MaxNumPoints',600);
opt=contset(opt,'Singularities',1);
[x,v,s,h,f]=cont(@equilibrium,x0,[],opt);
It obviously depends on the system. But could there be a rookie matcont user error that is leading to this?
Matcont finds a limit point at every point n the curve.
Sent from sourceforge.net because you indicated interest in https://sourceforge.net/p/matcont/discussion/762214/
To unsubscribe from further messages, please visit https://sourceforge.net/auth/subscriptions/
By time history I mean time integration of the states. It results in a particular periodic orbit and the limit points seem to follow the same orbit.
Thanks! I will try that out.
Last edit: Anonymous 2019-07-11
That suggests your system has some additional features... MatCont focusses on generic dissipative systems. If your system has more structure you should take that into account yourself.
Thanks for the tip! But could you be a little more clear about what exactly you mean by 'more structure'?