Stability, Chaos Diagnose and Adaptive Control of Two Dimensional Discrete - Time Dynamical System ()
1. Introduction
In the last two decades, the interest in dynamical systems has increased, because they are an important concept in describing the behavior of many models and in various fields. Some studies and research have focused on discrete dynamical systems in which the systems are described in the form of difference equations [1], as intermittent dynamic systems are more appropriate in describing models, especially when Models have non-overlapping generations, as the computational techniques used in discrete dynamical systems are more effective for numerical simulations, unlike in continuous dynamic systems. Mathematical models have been studied in biology and the environment, which have been described in the form of food chains. Food chains is one of the basic relationships that describe for us the interactions that occur between societies in the field of ecology, which studies the interrelationships between living organisms and the environment in which they are. Malthus [2] presented a study on the interactions that occur in societies that consist of one species, which were described using difference mathematical equation. There are many mathematical models in ecology that have been touched upon by researchers in this field. One of the most important mathematical models that have received wide attention in recent years is the prey and predator model due to its wide scope of application. The interaction between the prey community and the predator was described by Volterra [3], Lutka [4]. This type of model was found in temperate regions due to its appropriate environment seasonality, as Holling presented [5] a more realistic model of prey predator describing for us the interaction of three types of societies, as the relationship in the intermittent dynamic systems between prey and predator societies was and remains one of the important topics in many studies and research because of its importance. This type of the models depends on many assumptions, including that prey societies grow in a limited way when the predator community is absent, and predation societies depend on the presence of prey in order to survive, that is, the rate of the prey community is proportional to the predator community, and finally, we can say that the models in the environmental science are not simple and we cannot rely on the results obtained as stable and perfect results.
Accordingly, the research was arranged as follows: system description, system analysis and finding fixed points [6], stability analysis [7] [8] [9] of the fixed points by (characteristic roots equation, jury criteria [10], Lyapunov function [10] ), numerical behavior study and finding diagrams of bifurcation of the system [11] [12], chaos analysis using maximum Lyapunov exponent and the Lyapunov dimension [13] and the binary test (0 - 1) [14] - [20]. In the final section, adaptive control [21] [22] of the chaotic system is performed, and stability tests and the Lyapunov exponent test are performed.
2. System Description
In this work, a two-dimensional discrete time dynamical system was taken [23] and defined as follows:
(1)
(2)
where
represents the prey society and
represents the predator community in discrete time (t), and that
represent the parameters of the system and are positive, the part:
is represents the rate of increase of the prey community in the absence of the predation community, and that Part
represents the rate of decline of the prey community due to the presence of the predation community, and
represents the predation parameter, while the part
represents the variance in the size of the predation population which depends on the size of the prey community.
3. System Analysis
3.1. Jacobian Matrix
For 2D discrete dynamical system (1) with continuous differentiable transition function
and
, given by
The Jacobian Matrix of system (1) is:
(3)
3.2. Fixed Point
In this section we find the fixed points of system (1), assume that
,
(4)
(5)
from Equations (4) and (5) we get the following fixed points:
,
,
Theorem (1): Let
(6)
A characteristic equation of (3), the following cases are true:
1) If the absolute value of the roots of Equation (6) is less than one, then the fixed point of the system (1) is locally asymptotically stable and is called the sink.
2) If the absolute value of the roots of Equation (6) is greater than one, then the fixed point of the system (1) is unstable and is usually called the source, but if at least one of the values of the roots of Equation (6) is greater than one, then the fixed point is called the Saddle.
3) If the absolute value of the roots of Equation (6) is equal to one, then the fixed point of the system (1) is called the non-hyperbolic point. But if there are no roots of values equal to one, then the fixed point is called the hyperbolic point.
4. Stability Analysis
In this section, the stability of the fixed points of the system (1) using the following criteria:
4.1. Characteristic Equation Roots
Substituting the point
into (3) we get:
And from Equation (2) we get:
(7)
so the roots of quadratic Equation (7) are:
,
since
,
, so by theorem (1) we get system (1) is unstable at
, similarly we test the points
and
, the results shown in Table 1 that the three points
are saddle.
Lemma (1):
Let the characteristic equation of system (1)
Then it’s jury Table 2 is:
Table 1. Stability of fixed points by characteristic equation.
Such that
,
,
,
,
We say that the fixed point of the system (1) is stable if the satisfies following conditions:
,
,
Otherwise, fixed points are unstable.
4.2. Jury Stability Criteria
We test the stability of point
, by using lemma (1) and values from eq. (7) we get the jury Table 3, since the condition
, so the point
is unstable. Similarly, we test the rest of the points
, which shows that they are unstable, so the system (1) is unstable.
4.3. Lyapunov Function Criteria
The fixed points are said to be stable using the criterion of the Lyapunov function if
except for the origin, which is stable, and to study the stability of the fixed points of the system (1) we impose the quadratic equation of the following Lyapunov function:
Where x, y are not equal to zero, and using ∆V we get:
(8)
By substituting (2) in Equation (8) and test the point
we get:
Table 3. Jury table of the point q0.
It is clear that V is positive definite and ∆V is positive, hence the fixed point
is unstable and so system (1) is unstable.
5. The Graphical and Numerical Behavior of System (1)
5.1. Newtons’ - Raphson Method
Newton’s Raphson method is one of the important numerical methods for finding the roots of the difference equations. By using a written programme on MATLAB we get the best result obtained for the system (1) is
with minimum error (0.0001).
5.2. Trajectories of System (1)
In this section the time behavior of the system (1) was studied and the parameters were fixed at the values
,
with the values
and for (1000) iterations, shown in Figure 1. It turns out that the system behaves unstable.
5.3. Phase Space of System (1)
In this section, the phase space of the system (1) was found, the parameters were fixed at the values
,
with the values
and shown in Figure 2 which show that the system behaves chaotic. In Figure 3
Figure 1. The trajectories of state variables
and
, when
.
Figure 2. The phase space of the system (1) at parameter
.
Figure 3. The phase space of the system (1) with some values at parameter
.
Figure 4. Bifurcation of parameter
with variable
.
shows that the system (1) generates trajectories chaotic for
and
at some values parameter
.
5.4. Bifurcation Diagram
In this section the bifurcation diagrams of system (1) are found at the bifurcation parameter
, the parameters
with the values
and the
parameter ranging from 1.8 to 2.7 showing the behavior chaotic of the system (1), also the internal balance of the parameter
on the period [2.6, 2.65], as shown in Figure 4.
6. Lyapunov Exponent and Lyapunov Dimension
The Lyapunov exponent is one of the important tests in detecting the chaotic behavior of dynamic systems, as the system is said to be a chaotic system if one of the values of the Lyapunov exponent is greater than zero, and by using a mathematical program in MATLAB, the Lyapunov exponent of the system (1) was tested and the following values were obtained:
,
Since one of the values of the Lyapunov exponent of system (1) is a positive value, then system is chaotic, show in Figure 5.
To calculate the Lyapunov dimension of the system (1), we use the following law:
7. Binary Test (0 - 1)
In this section, the binary test was used to analysis the chaos of the system (1), a time series
was generated from the regularity (1) at the parameters
,
and the values
For (1000) iterations, and by using a mathematical program in MATLAB, the binary selection of the system (1), we calculate
with
for t = 100 and chose(c) is random value within the period (0, π), show in Figure 6(a), which shows us the behavior aligned to the system and similar to Brownian movement, and the average square displacement
was calculated with time k, shown in Figure 6(b), which shows us that the average square of the displacement
It is a function that grows linearly with time (t), and by finding the mean (k) of the contiguous growth
of the mean of the displacement square
where t = 100 shows that the system behaves in chaotic behavior (k = 0.982) as the mean (k) approaches the one shown in Figure 6(c).
8. Adaptive Control Technique
To achieve the stability of the chaotic system (1) we will design an adaptive control law with the unknown parameter
.
(9)
where
are the controllers for the adaptive and are known as follows:
Figure 5. Lyapunov exponent of the system (1).
Figure 6. The binary test of system (1) with
.
(10)
where
are positive constants and the parameter
is an approximate parameter of the parameter
, and by substituting (10) in (9) we get:
(11)
Let the error for the discretionary parameter be defined as follows:
(12)
Substituting (12) into (11) we get:
(13)
8.1. Numerical Results
In this section, we will test the stability of the fixed points of the system (1) in the controlled system (11) with the values
and
,
and the parameter
is an estimated parameter of the parameter
, and let
.
8.1.1. Characteristic Equation Roots
From (11) the Jacobian matrix is:
(14)
Table 4 shows that the testing of the
of system (11) are stable.
8.1.2. Jury’s Test
We test the stability of the fixed point
, from characteristic equation
at point
we get
,
,
Then we get Jury Table 5, since, all conditions are satisfies lemma (1), therefore the fixed point
is stable. Similarly, we test the rest of the points
, shown that all are stable.
8.1.3. Lyapunov Function Test
From the system (11) we get
Substituting the parameters
we obtain
We test the fixed points
, it is clear that V is positive definite and ∆V is negative definite, consequently the adaptive strategy success to control system (1), and Table 6 illustrates this.
Table 4. Stability of fixed points by characteristic equation for system (11).
Table 5. The results Jury test of q0.
Table 6. Results of the Lyapunov function test of the system (11).
8.2. Lyapunov Exponent
In this section, the Lyapunov exponent test of the system was performed and the values were obtained
,
, accordingly, the system is regular and Figure 7 shows that.
Figure 7. Lyapunov exponent of system (11).
9. Conclusion
In this research, 2D discrete - time dynamical system was taken, the system was analyzed, fixed points were found, and the stability analyzed for fixed points using (roots characteristic equation, Jury test, Lyapunov function test). The roots of the system were found using Newton’s Raphson numerical method, and the dynamic behavior was analyzed and studied. The phase space of the system shows that the system is unstable at parameter
. For the chaos analysis of the system, the bifurcation diagrams of the bifurcation parameter
are found for the system, and the Lyapunov exponent test was used and the value (
) was obtained, which is an indication of the chaos of the system. And Lyapunov dimension was calculated as (
). When using the binary test (0 - 1), it was found that the value of (k = 0.982) with the parameter
which is an indication of the chaos of the system. Finally, the adaptive control of the system and the stability test of the system after the control was performed which showed us that the system is stable, and when we tested the Lyapunov exponent (
,
) shown the regular behavior of the system.
Acknowledgements
We would like sincerely thank and also acknowledge to the university of Mosul and college of computer science and mathematics for the support and encouragement that help us to improve the quality of this work.