Continuous dynamical systems that involve differential equations mostly contain parameters. It can happen that a slight variation in a parameter can have significant impact on the solution. The main questions of interest in this chapter are: How to continue equilibria and periodic orbits of dynamical systems with respect to a parameter? How to compute stability boundaries of equilibria and limit cycles in the parameter space? How to predict qualitative changes in system’s behavior (bifurcations) occurring at these equilibrium points? This chapter will also cover the classification of bifurcations in terms of equilibria and periodic orbits. Especially it will present the specific bifurcation called ” Hopf bifurcation” which refers to the development of periodic orbits from stable equilibrium point, as a bifurcation parameter crosses a critical value. Since the theory of bifurcation from equilibria based on center manifold reduction and Poincaré-Normal forms, the direction of bifurcations for the mathematical models will also be explained using this theory. Finally, by introducing several software packages and numerical methods this chapter will also cover the techniques to determine and continue in some control parameters all local bifurcations of periodic orbits of dynamical systems and relevant normal form computations combined with the center manifold theorem, including periodic normal forms for periodic orbits.
In general, in a dynamical system, a parameter is allowed to vary, then the differential system may change. An equilibrium can become unstable and a periodic solution may appear or a new stable equilibrium may appear making the previous equilibrium unstable. The value of parameter at which these changes occur is known as ”bifurcation value” and the parameter that is varied is known as the ”bifurcation parameter”.
In this chapter, we also discuss several types of bifurcations, saddle node, transcritical, pitchfork and Hopf bifurcation. Among these types, we especially focus on Hopf bifurcation. The first three types of bifurcation occur in scalar and in systems of differential equations. The fourth type called Hopf bifurcation does not occur in scalar differential equations because this type of bifurcation involves a change to a periodic solution. Scalar autonomous differential equations can not have periodic solutions. Hopf bifurcation occurs in systems of differential equations consisting of two or more equations. This type is also referred to as a ”Poincare-Andronov-Hopf bifurcation”.
For a given system of differential equations first we shall consider the stability and the local Hopf bifurcation. By using the Hopf bifurcation theorem we prove the occurrence of the Hopf bifurcation. And then, based on the normal form method and the center manifold reduction introduced by Hassard et al.,, we derive the formulae determining the direction, stability and the period of the bifurcating periodic solution at the critical value of the bifurcation parameter. To verify the theoretical analysis, numerical simulations for bifurcation analysis are given in this chapter. For references see -.
We also introduce the Hopf bifurcation for continuous dynamical systems and state the Hopf bifurcation theorem for these models. As it is well known, Hopf bifurcations occur when a conjugated complex pair of eigenvalues crosses the boundary of stability. In the time-continuous case, a limit cycle bifurcates. It has an angular frequency which is given by the imaginary part of the crossing pair. In the discrete case, the bifurcating orbit is generally quasi-periodic, except that the argument of the crossing pair times an integer gives just . If we consider an ordinary differential equation (ODE) that depends on one or more parameters
where, for simplicity, we assume to be the only parameter. There is the possibility that under variation of nothing interesting happens to Equation (1). There is only a quantitatively different behavior. Let us define Equation (1) to be structurally stable in the case there are no qualitative changes occurring. However, the ODE (Ordinary Differential Equation) might change qualitatively. At that point, bifurcations will have occurred.
Many of the basic principles for one dimensional systems apply also for two-dimensional systems. Let us define a two-dimensional system
where biologically we mostly interpret x as prey or resource and y as predator or consumer. Equilibria can be found by taking the equations equal to zero, i.e.,
We have three possibilities for the stability of an equilibrium. Next to the stable and unstable equilibrium, there is the saddle equilibrium. A two-dimensional stable equilibrium is attracting in two directions, while a two-dimensional unstable equilibrium is repelling in two directions. A saddle point is attracting in one direction and repelling in the other direction. In the less formal literature saddles are often considered just unstable equilibria. A second remark is that also the dynamics of the system around the equilibria can differ.
The attracting or repelling can occur via straight orbits (a node) or via spiralling orbits (a spiral or focus). Note that, it is not possible to have a saddle focus in two dimensions. It is possible though in three of higher dimensional systems.
To prove the existence of Hopf bifurcation, we first obtain the Hopf bifurcation theorem hypothesis, i.e., the existence of purely imaginary eigenvalues of the corresponding characteristic equation with respect to the parameter and also we prove the transversality condition
at the critical value where Hopf bifurcation occurs. Then based on the normal form approach and the center manifold theory introduced by Hassard et. al,, we derive the formula for determining the properties of Hopf bifurcation of the model.
Finally in this chapter, to support these theoretical results, we illustrate them by numerical simulations. In numerical analysis, generally MATLAB solver packages are used to analyze the dynamics of nonlinear models. In these solvers, differential equation systems are simulated by difference equations. In this chapter, we also give some examples from biology (such as well known predator-prey models with time delay) with numerical simulations and by graphing the solutions in two or three dimensions, we illustrate the occurrence of periodic solutions. For some examples see references by Çelik, ,  and .
2. Basic concepts of bifurcation analysis
As it is stated above, in dynamical systems, a bifurcation occurs when a small smooth change made to the parameter values (the bifurcation parameters) of a system causes a sudden ” qualitative" or topological change in its behavior. Generally, at a bifurcation, the local stability properties of equilibria, periodic orbits or other invariant sets changes. It has two types;
2.1. Equilibrium points
In dynamical systems, only the solutions of linear systems may be found explicitly. The problem is that in general real life problems may only be modeled by nonlinear systems. The main idea is to approximate a nonlinear system by a linear one (around the equilibrium point). Of course, we do hope that the behavior of the solutions of the linear system will be the same as the nonlinear one. But this is not always true. Before the linear stability analysis, we give some basic definitions below.
where is a function mapping A point is called an equilibrium point if there is a specific such that
Suppose is an equilibrium point (with the input ). Consider the initial condition , and applying the input for all , then resulting solution satisfies
for all . That is why it is called an equilibrium point or solution.
where denotes the population density at time , and are positive constants, is the carrying capacity. Then by setting right hand side function equal to zero,
we obtain two equilibrium points and
2.2. Linear stability analysis
Linear stability of dynamical equations can be analyzed in two parts: one for scalar equations and the other for two dimensional systems;
2.2.1. Linear stability analysis for scalar equations
To analyze the ODE
locally about the equilibrium point , we expand the function in a Taylor series about the equilibrium point . To emphasize that we are doing a local analysis, it is customary to make a change of variables from the dependent variable to a local variable. Now let
where it is assumed that , so that we can justify dropping all terms of order two and higher in the expansion. Substituting into the RHS of the ODE yields;
and dropping higher order terms, we obtain
Note that dropping these higher order terms is valid since Now substituting into the LHS of the ODE,
The goal is to determine if we have growing or decaying solutions. If the solutions grows, then the equilibrium point is unstable. If the solution decays, then the fixed point is stable.
To determine whether or not the solution is stable or unstable we simply solve the ODE and get the solution as
where is a constant. Hence, the solution is growing if and decaying if As a result, the equilibrium point is stable if unstable if as it is stated in the following theorem.
the derivative function is continuous on an open interval where the equilibrium point . Then the equilibrium point is locally stable if and it is unstable if
If the equilibrium point is stable and in addition
then it is called asymptotically stable equilibrium point.
We first find the equilibrium points by setting
Next we compute and evaluate it at the equilibrium points.
2.2.2. Linear stability analysis for systems
Consider the two dimensional nonlinear system
and suppose that is a steady state (equilibrium point), i.e.,
Now let’s consider a small perturbation from the steady state
where and are understood to be small as and Is is natural to ask whether and are growing or decaying so that and will move away form the steady state or move towards the steady states. If it moves away, it is called unstable equilibrium point, if it moves towards the equilibrium point, then it is called stable equilibrium point.As in scalar equations, by expanding the Taylor’s series for and
Since and are assumed to be small, the higher order terms are extremely small, we can neglect the higher order terms and obtain the following linear system of equations governingthe evolution of the perturbations and
where the matrix is called Jacobian matrix of the nonlinear system.The above linear system for and has the trivial steady state , and the stability of this trivial steady state is determined by the eigenvalues of the Jacobian matrix, as follows:
As a summary,
Hyperbolic equilibria are robust(i.e., the system is structurally stable): Small perturbations of order do not change qualitatively the phase portrait near the equilibria. Moreover, local phase portrait of a hyperbolic equilibrium of a nonlinear system is equivalent to that of its linearization. This statement has a mathematically precise form known as the Hartman-Grobman. This theorem guarantees that the stability of the steady state of the nonlinear system is the same as the stability of the trivial steady state of the linearized system.
Non-hyperbolic equilibria are not robust (i.e., the system is not structurally stable). Small perturbations can result in a local bifurcation of a non-hyperbolic equilibrium, i.e., it can change stability, disappear, or split into many equilibria. Some refer to such an equilibrium by the name of the bifurcation (See the section below).
The equilibria are the points and and the Jacobian matrix is
We compute the Jacobian at the equilibrium point (0,0) where which implies that the eigenvalues are purely imaginary
by solving the characteristic equation
Since the point is non-hyperbolic, the linearized system can not tell about the stability. Later on we will show that this is a center.
For the equilibrium point the Jacobian thus the point is locally unstable (as
where one of the eigenvalues is strictly positive). Since it is a hyperbolic equilibrium point, the stability of fixed point is the same as in the linearized system. So it is also unstable.
3. Bifurcation analysis
Consider a family of ODE’s that depend on one parameter
where is analytic for , . Let be a family of equilibrium points of equation (5), i.e., Now let’s set
whereLet be the eigenvalues of If, for some , changes sign at we say that is a bifurcation point of equation (5).
3.1. Bifurcation in one dimension
We may assume that so that and is a real valued analytic function of provided
Therefore, the equilibrium point is asymptotically stable if and unstable if This implies that is a bifurcation point if Hence, bifurcation points are solutions of
The most common bifurcation types are illustrated by the following examples.
3.1.1. Saddle-node bifurcation
A saddle-node bifurcation or tangent bifurcation is a collision and disappearance of two equilibria in dynamical systems. In autonomous systems, this occurs when the critical equilibrium has one zero eigenvalue. This phenomenon is also called fold or limit point bifurcation.
Now consider the dynamical system defined by
An equilibrium solution (where is simply Therefore,if then we have no real solutions,if then we have two real solutions.
We now consider each of the two solutions for , and examine their linear stability in the usual way. First, we add a small perturbation:
Substituting this into the equation yields
and since the term in brackets on the RHS is trivially zero, therefore
which has the solution
From this, we see that for , as (linear stability); for ,as (linear stability).
As sketched in the “ bifurcation diagram” below, therefore, the saddle node bifurcation at corresponds to the creation of two new solution branches. One of these is linearly stable, the other is linearly unstable.
3.1.2. Transcritical bifurcation
In a transcritical bifurcation, two families of fixed points collide and exchange their stability properties. The family that was stable before the bifurcation is unstable after it. The other fixed point goes from being unstable to being stable.
Now consider the dynamical system
Again, and are control parameters. We can find two steady states to this system
We now examine the linear stability of each of these states in turn, following the usual procedure.
For the state , we add a small perturbation
with the linearized form
has the solution
Therefore, perturbations grow for a > 0 and decay for a < 0. So
Now for the state we add a small perturbation
which yields the linearized form
has the solution
Therefore, perturbations grow for a > 0 and decay for a < 0. So
It can be easily seen that the bifurcation point corresponds to an exchange of stabilities between the two solution branches.
3.1.3. The pitchfork bifurcation
In pitchfork bifurcation one family of fixed points transfers its stability properties to two families after or before the bifurcation point. If this occurs after the bifurcation point, then pitchfork bifurcation is called supercritical. Similarly, a pitchfork bifurcation is called subcritical if the nontrivial fixed points occur for values of the parameter lower than the bifurcation value. In other words, the cases in which the emerging nontrivial equilibria are stable are called supercritical whereas the cases in which these equilibria are called subcritical.
Consider the dynamical system
As usual, and are external control parameters. Steady states, for which are as follows:
Note that the equilibrium points and only exist when if and for if .
As usual, we now examine the linear stability of each of these steady states in turn. (This can be done for a general ). First we write the perturbation for
that yields the linearized equation
with the solution
So we see that
For the states and , setting
with the solution
Thus it is obvious that
3.2. Bifurcation in two dimension
3.2.1. Hopf bifurcation
3.3. Hopf bifurcation theorem
Consider the two dimensional system
where is the parameter and suppose that is the equilibrium point and are the eigenvalues of the Jacobian matrix which is evaluated at the equilibrium point.
In addition let’s assume that the change in the stability of the equilibrium point occurs at where
First the system is transformed so that the equilibrium is at the origin and the parameter at gives purely imaginary eigenvalues. System (7) is rewritten as follows;
The linearization of the system (7) about the origin is given by , where and
is the Jacobian matrix evaluated at origin.
Let and , in system (8) have continuous third order partial derivatives in and Suppose that the origin is an equilibrium point of (8) and that the Jacobian matrix as above, is valid for all sufficiently small Moreover, assume that the eigenvalues of matrix are where such that the eigenvalues cross the imaginary axis with nonzero speed, i.e.,
Then in any open set containing the origin in and for any there exists a value such that the system of differential equations (8) has a periodic solution for in . (Allen, L.J.S).
In a supercritical Hopf bifurcation, the limit cycle grows out of the equilibrium point. In other words, right at the parameters of the Hopf bifurcation, the limit cycle has zero amplitude, and this amplitude grows as the parameters move further into the limit-cycle. (See the figure below)
However in a subcritical Hopf bifurcation, there is an unstable limit cycle surrounding the equilibrium point, and a stable limit cycle surrounding that. The unstable limit cycle shrinks down to the equilibrium point, which becomes unstable in the process. For systems started near the equilibrium point, the result is a sudden change in behavior from approach to a stable focus, to large-amplitude oscillations.(See the figure below).
where is a parameter. When we compute the equilibrium points depending on parameterIf then there is no x-nullclines, hence the system has no equilibrium points.If then the system has exactly one equilibrium point at If then the system has two equilibrium points and
Then the Jacobian matrix is
where at the equilibrium points andFor , there will be a line equilibrium (since one of the eigenvalues is zero) and for , the point is a sink and is a saddle point so that is the bifurcation point for this differential equation system.
where is the bifurcation parameter. We can easily show that the conditions of the Hopf Bifurcation theorem hold. In this system and are zero. Then the Jacobian matrix is
for which the eigenvalues are where Reand the imaginary part ImIt follows that Reand Imand also
Hence, we conclude that there exists a periodic solution for in every neighborhood of origin.
4. Center Manifold teoremLet , for the dynamical system
and let the eigenvalues of the Jacobian matrix be Suppose that, the real parts of the eigenvalues are zero and if not, suppose there are numbers of eigenvalues with number of eigenvalues with and number of eigenvalues with Let be the eigenspace on imaginary axis corresponding to eigenvalues. The eigenvalues on the imaginary axis are called the critical eigenvalues as on the eigenspace And suppose the function denote the flow corresponding to the equation .
With these assumption, we state the Center Manifold theorem as follows;
There exists a locally invariant center manifold such that
such that the dynamics of the system
(where and are are the blocks in the canonical form whose diagonals contain the eigenvalues with and ; respectively) restricted to the center manifold are given by
And the manifold is called center manifold.
Remark: Center manifolds are not unique.
4.1. Center Manifold reduction for two dimensional systems
The only equilibrium point is , we linearize around that and obtain
Now we look for Then,
and on the other hand,
Comparing the two expressions, we deduce that and the center manifold reduction takes the form
Hence for the last equation is asymptotically stable and therefore is asymptotically stable for the original system.
Again is an equilibrium point and the jacobian matrix for the linearized system is
Consider the transformation which leadsand
from which we deduce that and and
For this reduced equation, is unstable and hence, is also unstable for the original system.
4.2. Center Manifold reduction for Hopf bifurcation
The aim of this section is to give a formal framework for the analytical bifurcation analysis of Hopf bifurcations in delay differential equations
with a single fixed time delay to be chosen as a bifurcation parameter. Characteristic equations of the delay differential equation form (14) are often studied in order to understand changes in the local stability of equilibria of certain delay differential equations. It is therefore important to determine the values of the delay at which there are roots with zero part. We give a general formalization of these calculations and determine closed form algebraic equations where the stability and amplitude of periodic solutions close to bifurcation can be calculated.
We shall determine the direction of Hopf bifurcation and the stability of the bifurcating periodic solutions by applying the normal form theory and the center manifold theorem by Hassard et al.,, and throughout this section, we assume that the three dimensional system of delay differential equations (14) undergoes Hopf bifurcations at the positive equilibrium at and is the corresponding purely imaginary root of the characteristic equation at the positive equilibrium . For the sake of simplicity, we use the notation for .
We first consider the system (14) by the transformation
which is equivalent to the following Functional Differential Equation(FDE) system in ,
where and are given respectively, byand
By Riesz representation theorem, there exists a function of bounded variation for such that
Indeed we may take
where is the Dirac delta function. For defineand
Then the system (15) is equivalent to
where for For define
and a bilinear inner product
where Then and are adjoint operators. Suppose that and are eigenvectors of and corresponding to and respectively. Then suppose that is the eigenvector of corresponding to then Then in the following, we use the theory by Hassard et al.,, to compute the coordinates describing center manifold at Define
On the center manifold , we have
where and are local coordinates for center manifold in the direction of and Note that is real if is real. We consider only real solutions. For the solution since and (15), we have
By using (17), we have , and
From the definition of , we have
and we evaluate
To determine , we need to compute and . By (15) and (18), we have
Note that on the center manifold near to the origin,
Thus we obtain,
By using (19), for
Comparing the coefficients with (20), we obtain the following
From (22) and (24) and the definition of we get
Noticing we evaluate by is a constant vector. From the definition of and (22), we obtain
where Next we compute and from (25) and (26) and determine the following values to investigate the qualities of bifurcating periodic solution in the center manifold at the critical value For this purpose, we express the direction of Hopf bifurcation in terms of and eigenvalues . And then we can evaluate the following values;
and we state this as in the following theorem.
In the following section, we shall give a numerical example to verify the theoretical results.
4.2.1. Numerical example of Center Manifold reduction
Consider the following system with discrete time delay ;
Now we present some numerical simulations by using MATLAB(7.6.0) programming (Çelik-3). We simulate the predator-prey system (28) by choosing the parameters , , , , and , i.e., we consider the following system
which has only one positive equilibrium . By algorithms in the previous sections, we obtain , ,and which leads to . So by the theorem above, the equilibrium point is asymptotically stable when and unstable when and also Hopf bifurcation occurs at as it is illustrated by computer simulations.
By the theory of Hassard et al.,, as it is discussed in previous section, we also determine the direction of Hopf bifurcation and the other properties of bifurcating periodic solutions. From the formulae in Section 5.2 we evaluate the values of , and as
from which we conclude that Hopf bifurcation of system (29) occurring at is supercritical and the bifurcating periodic solution exists when crosses to the right, and also the bifurcating periodic solution is stable.
In computer simulations, the initial conditions are taken as and MATLAB DDE (Delay Differential Equations) solver is used to simulate the system (29). We first take and plot the density functions , and in Figs.9,10,11 respectively which shows the positive equilibrium is asymptotically stable for . Moreover in Fig.12, we illustrate the asymptotic stability in three dimension.
However in Figs.13,14,15 and 16 below, we take sufficiently close to which illustrates the existence of bifurcating periodic solutions from the equilibrium point .
1. Consider the Van der Pol equation
and convert this into a system and check the stability of fixed points.
2. Let be an equilibrium point of the equationLet be a positive function. Deduce the stability of as an equilibrium solution of the equation
from its stability as an equilibrium solution of the equation Repeat the same question if is a negative function.
3. For the following nonlinear system
determine the stability of fixed points.
4. For the following nonlinear system
classify the fixed points.
5. Analyze the bifurcation properties of the following problems choosing as bifurcation parameter,
6. Find the equilibrium points and identify the bifurcation in the following system, and sketch the appropriate bifurcation diagram and phase portraits:
Then compute the extended center manifold near the bifurcation point by choosing as bifurcation parameter.
7. Show that the following system is structurally unstable,
and the following system is structural stable. Explain our reason.
for the above system, classify the fixed points.
9. For , consider the differential equation
on the real line.
Show that is a fixed point for any value of the parameter r, and determine its stability. Hence identify a bifurcation point
Show that for certain values of the parameter there are additional fixed points. For which values of do these fixed points exist? Determine their stability and identify a further bifurcation points
Using a Taylor expansion of the differential equation above, determine the normal form of the bifurcation at . What type of bifurcation takes place.
Similarly, determine the normal from of the bifurcation at What type bifurcation takes place?
Sketch the bifurcation diagram for all values of and . (Use a full line to denote a curve of stable fixed points, and a dashed line for a curve of unstable fixed points).
10. Consider the system of differential equations
in the plane.
Determine all fixed points of the system.
Sketch the phase portrait in the plane, including trajectories through and Whichfixed point does the trajectory through approach?
11. Consider the Lorenz system (the model of heat convection by Rayleigh-Benard occurring in the earth’s atmosphere) in three dimension as follows;
where are constants. Pom the stability analysis of this nonlinear system.
12. For the Holling-Tanner type Predator-Prey model
where and are positive constants. Find the equilibria and classify them.
13. Consider the Delayed Predator-Prey model
with time delay and positive constants and By choosing as bifurcation parameter, check whether bifurcating periodic solutions occur around the equilibrium points or not.