Limitcycle initialization
For limit cycles the same problems occur as with equilibria, a limit cycle continuation can't be done by just calling the continuer as
[x,v,s,h,f]=cont(@limitcycle, x0, v0, opt)
The limit cycle curve file has to know
- which ode file to use,
- which parameter is active,
- the values of all parameters,
- the number of mesh and collocation points to use for the discretization.
Also an initial cycle x0 has to be known. All this information can be supplied using any of the following three starting functions. All three
return an initial cycle x0 as well as its tangent vector v0.
- [x0,v0]=init_H_LC(@odefile, x, p, ap, h, ntst, ncol)
Calculates an initial cycle from a Hopf point detected on an equilibrium curve. Here odefile is the ode-file to be used. x is a vector containing the values of the state variables returned by a previous equilibrium curve continuation. p is the vector containing the current values of the parameters. ap is the active parameter, h contains the value of the initial amplitude and ntst and ncol are the number of mesh and collocation points to be used for the discretization.
- [x0,v0]=init_BPC_LC(@odefile, x, v, s, ap, ntst, ncol,h)
Calculates an initial cycle for starting a secondary cycle from a BPC detected on a previous calculated limit cycle. odefile, ap, ntst and ncol are the same as for init_H_LC. x, v and s are here the x, v and s as returned by a previous limit cycle continuation and h contains the value of the initial amplitude.
- [x0,v0]=init_LC_LC(@odefile, x, v, s, par, ap, ntst, ncol)
This starter can be used to start a limit cycle continuation from a previous limit cycle continuation. odefile, ap, ntst and ncol are the same as for init_H_LC. x and v are here the x and v as returned by a previous limit cycle continuation. s is the special point structure for the
point from where to start the new continuation. par is the parameter vector that is to be used for the new continuation. If this is the same as in the cycle returned by the continuer (arguably the natural situation), then this is the same as
s.data.parametervalues
- [x0,v0]=init_PD_LC(@odefile, x, s, ntst, ncol,h)
This starter calculates an initial double period cycle from a period doubling bifurcation point. Here h is a distance parameter which can be set to zero. We note also that there is no active parameter ap
input argument.