Normal form coefficients

The numerical methods used in MATCONT to compute normal form coefficients of codimension 1 bifurcations of limit cycles are discussed in [28].

The periodic normal form at the Limit Point of Cycles (LPC) bifurcation is

\begin{displaymath}
\left\{\begin{array}{rcl}
\frac{d\tau}{dt}&=&1 - \xi + a...
...
\frac{d\xi}{dt}&=& b \xi^2 + \cdots ,
\end{array}
\right.
\end{displaymath} (51)

where $\tau \in [0,T]$, $\xi$ is a real coordinate on the center manifold that is transverse to the limitcycle, $a,b\in \mathbb{R}$, and dots denote nonautonomous $T$-periodic $O(\xi^3)$-terms. For each detected LPC point the normal form coefficient $b$ is computed. If $b\neq 0$ then the LPC bifurcation is nondegenerate, i.e., the cycle manifold is quadratically tangential to the hyperplane orthogonal to the parameter direction. In practice, by round-off errors $b$ is never zero, so it is better to check if the testfunction changes sign. (for $b=0$ the LPC degenerates to a cusp of cycles (CP) but this will not be detected on a branch of limit cycles since it is a codimension 2 phenomenon.)

The periodic normal form at the Period Doubling (PD) bifurcation is

\begin{displaymath}
\left\{\begin{array}{rcl}
\frac{d\tau}{dt}&=&1 + a \xi^2...
...
\frac{d\xi}{dt}&=& c \xi^3 + \cdots,
\end{array}
\right.
\end{displaymath} (52)

where $\tau \in [0,2T]$, $\xi$ is a real coordinate on the center manifold that is transverse to the limit cycle, $a,c \in \mathbb{R}$, and dots denote nonautonomous $2T$-periodic $O(\xi^4)$-terms. The coefficient $c$ determines the stability of the period doubled cycle in the center manifold and is computed during the processing of each PD point.

If $c<0$ then the PD bifurcation is supercritical, i.e., within the center manifold on one side of the bifurcation only stable cycles exist and on the other side unstable cycles coexist with stable period-doubled cycles. If $c>0$ the PD bifurcation is subcritical, i.e., within the center manifold on one side of the bifurcation stable cycles coexist with unstable period-doubled cycles and on the other side only unstable cycles exist. (for $c=0$ the PD degenerates to Generalized Period Doubling bifurcation (GPD) but this will not be detected on a branch of limit cycles since this is a codimension 2 phenomenon).

The periodic normal form at the Neimark-Sacker (NS) bifurcation is

\begin{displaymath}
\left\{\begin{array}{rcl}
\frac{d\tau}{dt}&=&1 + a \vert...
...} \xi + d \xi \vert\xi\vert^2 + \cdots,
\end{array}
\right.
\end{displaymath} (53)

where $\tau \in [0,T]$, $\xi$ is a complex coordinate on the center manifold that is complementary to $\tau$, $a \in \mathbb{R},d \in \mathbb{C}$, and dots denote nonautonomous $T$-periodic $O(\vert\xi\vert^4)$-terms. The critical coefficient $d$ in the periodic normal form for the NS bifurcation is computed during the processing of a NS point.

If ${\rm Re}~d<0$ then the bifurcation is supercritical, i.e. in the center manifold the limit cycle is stable on one side of the bifurcation point and unstable on the other side and the unstable cycles coexist with stable invariant tori. If ${\rm Re}~d>0$ then the bifurcation is subcritical, i.e. in the center manifold the limit cycle is stable on one side of the bifurcation point and unstable on the other side and the stable cycles coexist with unstable invariant tori. (for ${\rm Re}~d=0$ the NS bifurcation degenerates to a Chenciner bifurcation (CH) but this will not be detected on a branch of limit cycles since it is a codimension 2 phenomenon).