On a limit cycle curve the following bifurcations can occur
The test function for the Period Doubling bifurcation is defined by
the following system
|
here j and y are so-called bordering vector-functions [10], see [7] for details on the implementation. The system is discretized using orthogonal collocation and solved using the standard MATLAB sparse system solver. The solution component G Î R of this system is the test function and equals zero when there is a Period Doubling bifurcation.
The Fold bifurcation is detected in the same way as the Fold bifurcation of equilibria, the last component of the tangent vector (the a component) is used as the test function.
The Neimark-Sacker bifurcation is detected by monitoring the eigenvalues of the monodromy matrix for the cycle. The monodromy matrix is computed like in auto by a block elimination in the discretized form of the Jacobian of (50).
BPC cycles are not generic in families of limit cycles, but they are common in the case of symmetries, if the
branching parameter is also the continuation parameter. CL_MatCont uses a strategy that requires only the solution of linear systems;
it is based on the fact that in a symmetry-breaking BPC cycle MD has rank defect two, where MD is the square matrix MD,
obtained from the discretized form of the Jacobian of (50).
To be precise, if h Î C1([0,1],IRn), then
|
and
|
where ()dm and ()dc denote discretization in mesh points
and in collocation points, respectively.
Therefore we border MD with two additional rows and columns to obtain
|
so that MDbb is nonsingular in the BPC cycle. Then we solve the systems
|
where y11,y12 have (Nm+1)n components, and gBPC11, gBPC12, gBPC21, and gBPC22, are scalar test functions for the BPC. In the BPC cycle they all vanish.
The singularity matrix is
|
The first row corresponds to the BPC. It contains 4 zeros which indicates that gBPC11, gBPC12, gBPC21, and gBPC22 should vanish. The last row corresponds to the NS. Because we have to exclude that all four testfunctions of the BPC are zeros, we introduce an extra testfunction which corresponds to the norm of these four testfunctions. A NS is detected if this norm is nonzero,the testfunction for the fold is nonzero and the testfunction for the NS is equal to zero.