7.7: Limit Cycles
- Page ID
- 91093
So far we have just been concerned with equilibrium solutions and their behavior. However, asymptotically stable fixed points are not the only attractors. There are other types of solutions, known as limit cycles, towards which a solution may tend. In this section we will look at some examples of these periodic solutions.
Such solutions are common in nature. Rayleigh investigated the problem
\[x^{\prime \prime}+c\left(\dfrac{1}{3}\left(x^{\prime}\right)^{2}-1\right) x^{\prime}+x=0 \nonumber \]
in the study of the vibrations of a violin string. Balthasar van der Pol (1889-1959) studied an electrical circuit, modeling this behavior. Others have looked into biological systems, such as neural systems, chemical reactions, such as Michaelis-Menten kinetics, and other chemical systems leading to chemical oscillations. One of the most important models in the historical study of dynamical systems is that of planetary motion and investigating the stability of planetary orbits. As is well known, these orbits are periodic.
Limit cycles are isolated periodic solutions towards which neighboring states might tend when stable. A key example exhibiting a limit cycle is given in the next example.
Find the limit cycle in the system
\[ \begin{aligned} x^{\prime} &=\mu x-y-x\left(x^{2}+y^{2}\right) \\ y^{\prime} &=x+\mu y-y\left(x^{2}+y^{2}\right) \end{aligned} \label{7.48} \]
Solution
It is clear that the origin is a fixed point. The Jacobian matrix is given as
\[D f(0,0)=\left(\begin{array}{cc} \mu & -1 \\ 1 & \mu \end{array}\right) \nonumber \]
The eigenvalues are found to be \(\lambda=\mu \pm i\). For \(\mu=0\) we have a center. For \(\mu<0\) we have a stable spiral and for \(\mu>0\) we have an unstable spiral. However, this spiral does not wander off to infinity. We see in Figure \(\PageIndex{1}\) that the equilibrium point is a spiral. However, in Figure \(\PageIndex{2}\) it is clear that the solution does not spiral out to infinity. It is bounded by a circle.
One can actually find the radius of this circle. This requires rewriting the system in polar form. Recall from Chapter 2 that we can change derivatives of Cartesian coordinates to derivatives of polar coordinates by using the relations
\[r r^{\prime} =x x^{\prime}+y y^{\prime} \nonumber \]
\[r^{2} \theta^{\prime} =x y^{\prime}-y x^{\prime} \nonumber \]
Inserting the system \(\PageIndex{2}\) into these expressions, we have
\[r r^{\prime}=\mu r^{2}-r^{4}, \quad r^{2} \theta^{\prime}=r^{2}. \nonumber \]
This leads to the system
\[ \begin{aligned} r^{\prime} &=\mu r-r^{3}, \\ \theta^{\prime} &=1 . \end{aligned}\label{7.52} \]
Of course, for a circle the radius is constant, \(r=\) const. Therefore, in order to find the limit cycle, we need to look at the equilibrium solutions of Equation \(\PageIndex{6}\). This amounts to finding the constant solutions of \(\mu r-r^{3}=0 .\) The equilibrium solutions are \(r=0, \pm \sqrt{\mu} .\) The limit cycle corresponds to the positive radius solution, \(r=\sqrt{\mu}\).
In Figures \(\PageIndex{1} - \PageIndex{2}\) we take \(\mu=0.4\). In this case we expect a circle with \(r=\sqrt{0.4} \approx 0.63\). From the \(\theta\) equation, we have that \(\theta^{\prime}>0\). This means that we follow the limit cycle in a counterclockwise direction as time increases.
Limit cycles are not always circles. In Figures \(\PageIndex{3} -\PageIndex{4}\) we show the behavior of the Rayleigh system \(\PageIndex{1}\) for \(c=0.4\) and \(c=2.0 .\) In this case we see that solutions tend towards a noncircular limit cycle in a clockwise direction.
A slight change of the Rayleigh system leads to the van der Pol equation:
\[x^{\prime \prime}+c\left(x^{2}-1\right) x^{\prime}+x=0 \nonumber \]
(The van der Pol system). The limit cycle for \(c=2.0\) is shown in Figure \(\PageIndex{5}\).
Can one determine ahead of time if a given nonlinear system will have a limit cycle? In order to answer this question, we will introduce some definitions.
We first describe different trajectories and families of trajectories. A flow on \(R^{2}\) is a function \(\phi\) that satisfies the following
- \(\phi(\mathbf{x}, t)\) is continuous in both arguments.
- \(\phi(\mathbf{x}, 0)=\mathbf{x}\) for all \(\mathbf{x} \in R^{2}\).
- \(\phi\left(\phi\left(\mathbf{x}, t_{1}\right), t_{2}\right)=\phi\left(\mathbf{x}, t_{1}+t_{2}\right)\).
(Orbits and trajectories). The orbit, or trajectory, through \(\mathbf{x}\) is defined as \(\gamma=\{\phi(\mathbf{x}, t) \mid t \in I\} .\) In Figure \(\PageIndex{6}\) we demonstrate these properties. For \(t=0, \phi(\mathbf{x}, 0)=\mathbf{x}\). Increasing \(t\), one follows the trajectory until one reaches the point \(\phi\left(\mathbf{x}, t_{1}\right)\). Continuing \(t_{2}\) further, one is then at \(\phi\left(\phi\left(\mathbf{x}, t_{1}\right), t_{2}\right)\). By the third property, this is the same as going from \(\mathbf{x}\) to \(\phi\left(\mathbf{x}, t_{1}+t_{2}\right)\) for \(t=t_{1}+t_{2}\).
Having defined the orbits, we need to define the asymptotic behavior of the orbit for both positive and negative large times. We define the positive semiorbit through \(\mathbf{x}\) as \(\gamma^{+}=\{\phi(\mathbf{x}, t) \mid t>0\} .\) The negative semiorbit through \(\mathbf{x}\) is defined as \(\gamma^{-}=\{\phi(\mathbf{x}, t) \mid t<0\} .\) Thus, we have \(\gamma=\gamma^{+} \cup \gamma^{-}\).
(Limit sets and limit points). The positive limit set, or w-limit set, of point \(\mathbf{x}\) is defined as
\[\Lambda^{+}=\left\{\mathbf{y} \mid \text { there exists a sequence of } t_{n} \rightarrow \infty \text { such that } \phi\left(\mathbf{x}, t_{n}\right) \rightarrow \mathbf{y}\right\} \nonumber \]
The \(\mathbf{y}^{\prime}\) s are referred to as w-limit points. This is shown in Figure \(\PageIndex{7}\).
Similarly, we define the negative limit set, or the alpha-limit set, of point \(\mathbf{x}\) is defined as
\[\Lambda^{-}=\left\{\mathbf{y} \mid \text { there exists a sequences of } t_{n} \rightarrow-\infty \text { such that } \phi\left(\mathbf{x}, t_{n}\right) \rightarrow \mathbf{y}\right\}\nonumber \]
and the corresponding \(\mathbf{y}^{\prime}\) s are \(\alpha\)-limit points. This is shown in Figure \(\PageIndex{8}\).
(Cycles and periodic orbits). There are several types of orbits that a system might possess. A cycle or periodic orbit is any closed orbit which is not an equilibrium point. A periodic orbit is stable if for every neighborhood of the orbit such that all nearby orbits stay inside the neighborhood. Otherwise, it is unstable. The orbit is asymptotically stable if all nearby orbits converge to the periodic orbit.
A limit cycle is a cycle which is the \(\alpha\) or \(\omega\)-limit set of some trajectory other than the limit cycle. A limit cycle \(\Gamma\) is stable if \(\Lambda^{+}=\Gamma\) for all \(\mathbf{x}\) in some neighborhood of \(\Gamma\). A limit cycle \(\Gamma\) is unstable if \(\Lambda^{-}=\Gamma\) for all \(\mathbf{x}\) in some neighborhood of \(\Gamma\). Finally, a limit cycle is semistable if it is attracting on one side and repelling on the other side. In the previous examples, we saw limit cycles that were stable. Figures \(\PageIndex{7}\) and \(\PageIndex{8}\)depict stable and unstable limit cycles, respectively.
We now state a theorem which describes the type of orbits we might find in our system.
Let \(\gamma^{+}\)be contained in \(a\) bounded region in which there are finitely many critical points. Then \(\Lambda^{+}\)is either
- a single critical point;
- a single closed orbit;
- a set of critical points joined by heteroclinic orbits. [Compare Figures \(\PageIndex{9}\) and \(\PageIndex{10}.\) ]
We are interested in determining when limit cycles may, or may not, exist. A consequence of the Poincaré-Bendixon Theorem is given by the following corollary.
Let \(D\) be a bounded closed set containing no critical points and suppose that \(\gamma^{+} \subset D .\) Then there exists a limit cycle contained in \(D .\)
More specific criteria allow us to determine if there is a limit cycle in a given region. These are given by Dulac’s Criteria and Bendixon’s Criteria.
Consider the autonomous planar system
\[x^{\prime}=f(x, y), \quad y^{\prime}=g(x, y)\nonumber \]
and a continuously differentiable function \(\psi\) def
\(D\) contained in some open set. If
\[\dfrac{\partial}{\partial x}(\psi f)+\dfrac{\partial}{\partial y}(\psi g) \nonumber \]
does not change sign in \(D\), then there is at most one limit cycle contained entirely in \(D .\)
Consider the autonomous planar system
\[x^{\prime}=f(x, y), \quad y^{\prime}=g(x, y)\nonumber \]
defined on a simply connected domain \(D\) such that
\[\dfrac{\partial}{\partial x}(\psi f)+\dfrac{\partial}{\partial y}(\psi g) \neq 0 \nonumber \]
in \(D\). Then, there are no limit cycles entirely in \(D\).
- Proof.
-
These are easily proved using Green’s Theorem in the Plane. (See your calculus text.) We prove Bendixon’s Criteria. Let \(\mathbf{f}=(f, g)\). Assume that \(\Gamma\) is a closed orbit lying in \(D .\) Let \(S\) be the interior of \(\Gamma .\) Then
\[ \begin{aligned} \int_{S} \nabla \cdot \mathbf{f} d x d y &=\oint_{\Gamma}(f d y-g d x) \\ &=\int_{0}^{T}(f \dot{y}-g \dot{x}) d t \\ &=\int_{0}^{T}(f g-g f) d t=0 . \end{aligned} \label{7.54} \]
So, if \(\nabla \cdot \mathbf{f}\) is not identically zero and does not change sign in \(S\), then from the continuity of \(\nabla \cdot \mathbf{f}\) in \(S\) we have that the right side above is either positive or negative. Thus, we have a contradiction and there is no closed orbit lying in \(D\).
Consider the earlier example in \((7 \cdot 48)\) with \(\mu=1\).
\[ \begin{aligned} &x^{\prime}=x-y-x\left(x^{2}+y^{2}\right) \\ &y^{\prime}=x+y-y\left(x^{2}+y^{2}\right) \end{aligned} \label{7.55} \]
We already know that a limit cycle exists at \(x^{2}+y^{2}=1 .\) A simple computation gives that
\[\nabla \cdot \mathbf{f}=2-4 x^{2}-4 y^{2} \nonumber \]
For an arbitrary annulus \(a<x^{2}+y^{2}<b\), we have
\[2-4 b<\nabla \cdot \mathbf{f}<2-4 a \nonumber \]
For \(a=3 / 4\) and \(b=5 / 4,-3<\nabla \cdot \mathbf{f}<-1\). Thus, \(\nabla \cdot \mathbf{f}<0\) in the annulus \(3 / 4<x^{2}+y^{2}<5 / 4\). Therefore, by Dulac’s Criteria there is at most one limit cycle in this annulus.
Consider the system
\[ \begin{aligned} &x^{\prime}=y \\ &y^{\prime}=-a x-b y+c x^{2}+d y^{2} . \end{aligned} \label{7.56} \]
Let \(\psi(x, y)=e^{-2 d x} .\) Then,
\[\dfrac{\partial}{\partial x}(\psi y)+\dfrac{\partial}{\partial y}\left(\psi\left(-a x-b y+c x^{2}+d y^{2}\right)\right)=-b e^{-2 d x} \neq 0 \nonumber \]
We conclude by Bendixon’s Criteria that there are no limit cycles for this system.