Chaos Behavior and Estimation of the Unknown Parameters of Stochastic Lattice Gas for Prey-Predator Model with Pair-Approximation

Show more

1. Introduction

A lattice model usually represents the motion of a network of particles, where the motion is produced by forces acting between the neighboring particles. Lattice models also are used to simulate the structure of polymers and can exhibit its dynamic behaviors.

The interaction between particles of the systems, which is the subject of this study, is continuous-time Masrkov process on certain spaces of configurations of particles. These systems began as a branch of probability theory in the 1960’s. Most of the inputs came from the work of Spitzer in United States [1] and of Dobrushin in the Soviet Union [2] . From a mathematical point of view, interacting particles systems represent a natural departure from the established theory of Markov processes. The lattice can be in one, two, three or complex dimensions. The one-dimensional lattices consider a linear arrangement between the sequences of particles, and each particle connected to the next one by a spring depends upon a condition. There are many classifications for the two-dimension lattice, such as a centered rectangular lattice, a hexagonal lattice, square lattice and so on (see Figure 1).

Many of models of interacting particles system have shown a chaos behavior. The problems of estimating and controlling stochastic systems are far from solved, and a considerable amount of research is under way. Estimation of the internal states of a stochastic dynamical system is a topic with important applications in different fields such as physics, biology and medicine [3] - [5] .

A spatial stochastic model to discuss strategies to control the epidemic was introduced by Schinazi [6] . El-Gohary has proposed a stochastic model to study the problem of optimal controlling the epidemic [3] [4] . El-Gohary and Al-Ruziza have suggested a non-linear stochastic model to investigate the optimal control of a non-homogenous prey-predator model. They have derived the feedback control law as non-linear functions of the population densities [7] [8] . El-Gohary has studied the problems of chaos and optimal control cancer model with complete unknown parameters [9] . Al-Mahdi and Khirallah have studied stability and bifurcation analysis of a model of cancer [10] [11] . Alwan and El-Gohary have studied the chaos, estimation and optimal control of habitat destruction model and genital herpes epidemic models with uncertain parameters [12] [13] .

A stochastic lattice gas model is proposed to describe the dynamics of two animals populations, one being a prey and the other a predator [14] . El-Gohary and Alwan have discussed the problem of chaos and control of a stochastic lattice gas model for prey- predator when one-site approximation is used. They have studied stability of the system and derived the optimal control inputs. They have also derived the estimators of the unknown parameters and probabilities [15] . This paper is considered as an extension of the paper [15] where it will discuss the stability and estimation of the unknown parameters of the stochastic lattice gas model for prey-predator when two-site approximation is used.

Figure 1. Linear, square and simple 3D lattices.

This paper has the following structure. In section 2, the stochastic probability model and its proposed rules will be discussed, and also the pair approximation mathematical model and its analytical solution will be presented. In section 3, stability analysis of the system will be studied and presented graphically. Estimation of the unknown parameters and the updating rules are derived in section 4. In section 5, numerical solutions are derived and presented graphically. Finally, conclusions are provided in Section 6.

2. Stochastic Probability Model

In this section, we will describe the stochastic rules of the proposed stochastic lattice gas model for prey-predator.

The lattice gas models describing special chemical reaction. Let us consider a lattice of N sites, every site can be either empty (O) or occupied by a prey (1) or occupied by a predator (2). At any time step a site is randomly chosen. For that site, we suppose that and are the numbers of the nearest neighbors of that site occupied by prey or predator, respectively and S is the total number of nearest neighbors of this site. The transition matrix of the model is given in Table 1 or graphically in Figure 2.

This Markov process contains three probability parameters, and, which are associated to the process: birth of prey process, death of prey and simultaneous birth of predator, and spontaneous death of the predator [14] - [17] . The parameters, and satisfy the condition:

Table 1. The transition matrix of the model.

Figure 2. The transitions of the model.

. (1)

Let us consider the system state as where according whether the state i empty, occupied by prey or occupied by predator respectively. Let be the probability of the state μ at time t. The time evolution of the probability of the state μand the mathematical system is given in details in [14] . The final mathematical expression of the model is a hierarchical system of equations, where the one-site correlations involve the two-site correlations and the two-site correlations involve three-site correlation and so on, where and take one values of 0, 1, or 2 (see Equation (6) in [14] ). The mean field approximation theory is used to derive a closed set of equations. The one-site approximation of this hierarchical system of Equation (6) in [14] is accomplished by writing the two-site correlation as the product in the hierarchical system. The two-site approximation model which is the subject of this paper will be discussed in the next section.

Pair Approximation Mathematical Model

In this section, the two-site approximation will be presented and the final mathematical model will be written.

In the hierarchical system Equation (6) in [14] , the mean field theory is suggest an approximation to derive a closed set of equations as: The thee-site correlations can be expressed in terms of two-site and one-site correlations as follows

(2)

(3)

where sites i and k are the nearest neighbors of site j. We also seek for spatially homogeneous and isotropic solutions of Equation (6) in [14] . In this case we may drop the indexes in and. We then have three one-site correlations, and nine two-site correlations, which are presented in the transition matrix in Table 1. However, only five of them are independent. We choose them to be (the prey density), (the predator density), , , and. Let us consider also (the vacuum site density). The equations for these variables and the final mathematical system according this pair approximation are given by:

(4a)

, (4b)

, (4c)

, (4d)

, (4e)

where

(5)

and S is the total number of nearest neighbors of the site. This is a nonlinear system of differential equations. The analytical solution for this system is given by solving the system of equation:

. (6)

By Maple program and [14] , when, the system in Equation (6) has the following general solution

, (7a)

, (7b)

, (7c)

, (7d)

, (7e)

where,

and

(for more details, see [14] ), and the following two trivial fixed stationary states as special solutions: corresponds to the vacuum-absorbing state, and corresponds to the prey-absorbing state.

3. Stability Analysis

Study of the stability and the chaos of the system will be discussed in this section, also some of the equilibrium point will be presented.

For the case of, the system has only special stationary solution which is prey- absorbing state [14] . This solution is unstable when, and the general stationary solutions:, , is unstable also, when as in [15] . For the case, the system (3) has two trivial fixed stationary states, that are given by:

, that correspond to the vacuum-absorbing state, and, that correspond to the prey-absorbing state. The Jacobian matrix of the system (4) is given by

(8)

where

The Jacobian matrix in Equation (8) of the system in Equation (4) evaluated at the vacuum-absorbing state is converges to

(9)

and its eigenvalues just are the elements of the main diameter, which are

. (10)

Using the linear stability analysis, since is a positive eigenvalue at least, hence this stationary solution is absolutely unstable.

Similarly, we get the Jacobian matrix in Equation (8) of the system in Equation (4) evaluated at the prey-absorbing state converges to. Accordingly has the same eigenvalues. Therefore, the stationary state is also absolutely unstable. The linear stability analysis for these two stationary states indicates sufficiently that, the system of stochastic lattice gas of prey-predator system according to the pair-appro- ximation is absolutely unstable at least in two dimensions. But the stability conditions are requiring more study.

Such behavior for the system in prey, predator and vacuum densities in Figures 3(a)-(c) respectively, represents an oscillatory behavior. For limit-cycle that appear in Figures 3(d)-(f) where all the neighboring trajectories tend to a limit-cycle at time tends to infinity, causing the so-called a stable limit-cycle, which indicates that the system stochastic lattice gas of prey-predator system according to the pair-approximation has an oscillatory behavior.

4. Estimations of the Unknown Probabilities

In this section we will derive the dynamic estimators of the unknown probabilities and from the conditions of the asymptotic stability of the system in Equation (4) about its stationary states assistance of some feedback variables.

At the beginning, let us assume the modified model with unknown probabilities in Equation (4) to become as follows

(11a)

(11b)

(11c)

(11d)

(11e)

where are the estimators of the unknown probabilities and are the control inputs that will be derived to make the trajectory of the system specified by the steady-states, and the general solution in Equation (7) to any of these states. If, then the system in Equation (11) has an unstable special solution:

Figure 3. The densities of the prey, predator and vacuum sites and their limit cycle, respectively for and, at the initial densities and

(12)

where are the steady-states of the uncontrolled system in Equation (4) that should be stabilized by finding the controllers

To solve the problem of this stabilization, we will use the Liapunov stability technique. Constructive upon this, let us consider the following positive definite form of Liapunov function

(13)

then

. (14)

By substituting from Equations (11) and choosing the following controllers inputs

(15a)

(15b)

(15c)

(15d)

(15e)

and the following update rules of the unknown probabilities

(16a)

(16b)

(16c)

the total time derivative of the Liapunov function in takes the form:

(17)

where and are non-negative control constants. In this case is a negative definite if and and a negative semi-definite if and. This implies that, under the action of the controllers in Equation (15) and updating rules in Equation (16) of the unknown system probabilities, the solution in Equation (7) of the systems in Equation (4) and Equation (17) is asymptotic stable in the Liapunov sense if and and only stable but not necessarily asymptotic stable if and. Since V is radially unbounded, and widen over time, therefore the solution in Equation (12) is globally asymptotically stable which completes the proof.

By substituting Equation (15) in the modified controlled system in Equation (11), in addition the update rules in Equation (16) we get the final system as follows

(18a)

(18b)

(18c)

(18d)

(18e)

(18f)

(18g)

(18h)

where and as defined in Equation (5). This system of nonlinear differential equations will be solved numerically and displayed some graphical solutions as examples.

5. Numerical Solution

This section presents some numerical solutions of the controlled nonlinear system of the stochastic lattice gas of prey-predator model in Equation (18) and the estimators of the system unknown probabilities to show how the control for this system is possible. Also, numerical examples for controlled stochastic lattice gas of prey-predator model were carried out for various probabilities values and different initial densities. For illustration purpose, we display the numerical solutions of the system graphically. Furthermore, the percentage error in the estimate for real values of the parameters will be calculated. The percentage error of the estimator of the parameter is given by the following rule:

. (19)

The following figures display two examples of numerical solutions of the non-linear system in Equation (18). The first solution is shown below.

Clearly, the densities to in Figures 4(b)-(e) respectively, tend to the steady states. While, the density in Figure 4(a) converges to the steady state. Also, the estimator in Figure 4(f) tends to the real value. While the estimators and in Figure 4(g) and Figure 4(h) converge to the real values and respectively. The steady state and the real values are represented by the dotted lines.

The second numerical example of solution is given below.

The densities and in Figure 5(a), Figure 5(d) and Figure 5(e) respectively, tend to the steady states. While, the densities and in Figure 5(b) and Figure 5(c) respectively, converge to the steady state. Also, the estimators

Figure 4. The controlled densities and the estimators of the unknown parameters of the stochastic lattice gas of prey-predator model with pair approximation for some values of parameters and initial densities as follows:

Figure 5. The controlled densities and the estimators of the unknown parameters of the stochastic lattice gas of prey-predator model with pair approximation for some values of parameters and initial densities as follows:

Table 2. The estimated values of the unknown parameters against the real values, and the percentage errors for some real values.

and in Figure 5(f) and Figure 5(g) tends to the real values and respectively. While the estimator in Figure 5(h) converges to the real values. The steady states and the real values are represented by the dotted lines. The previous figures show the approaching of the trajectories of the system to its steady states and approaching of the estimated values to the real values of the system unknown parameters over time.

In the following, numerical calculation for the percentage error in each estimate for different real values of the system unknown parameters.

Table 2 represents the comparing between the real values of the unknown parameters and its estimated values, in addition to the infinitesimal values of the percentage errors indicate to a strong convergence between the assumed real values and its estimated values so they are almost equally. These so good results are shown a high efficiency of the proposed Liapunov technique in the estimation process.

6. Conclusion

In this paper, we have introduced the mathematical model of stochastic lattice gas of prey-predator model with pair-approximation. The stability for this model has been discussed and it is found that this system has a chaos behavior. The estimators of the unknown parameters and the updating rules are derived according the conditions of the asymptotic stability. Some numerical solutions to show the stable system are presented graphically.

References

[1] Dobrushin, R.L. (1971) Markov Processes with a Large Number of Locally Interacting Components: Existence of a Limit Process and Its Ergodicity. Problems of Information Transmission, 7, 149-164.

[2] Dobrushin, R.L. (1971) Markov Processes with Many Locally Interacting Components-Reversible Cases and Some Generalization. Problems of Information Transmission, 7, 235-241.

[3] El-Gohary, A. (2001) Optimal Control of the Genital Herpes Epidemic. Chaos, Solitons Fractals, 12, 1817-1822.

http://dx.doi.org/10.1016/S0960-0779(00)00012-6

[4] El-Gohary, A. (2009) Chaos and Optimal Control of Equilibrium States of Tumor System with Drug. Chaos, Solitons and Fractals, 41, 425-435.

http://dx.doi.org/10.1016/j.chaos.2008.02.003

[5] El-Gohary, A. and Alwasel, I.A. (2009) The Chaos and Optimal Control of Cancer Model with Complete Unknown Parameters. Chaos, Solitons and Fractals, 42, 2865-2874.

http://dx.doi.org/10.1016/j.chaos.2009.04.028

[6] Schinazi, R. (1997) Predatorprey and Host Parasite Spatial Stochastic Models. Annals of Applied Probability, 7, 1-9.

http://dx.doi.org/10.1214/aoap/1034625250

[7] El-Gohary, A. and Al-Ruziza, A. (2003) Optimal Control of Non-Homogenous Prey Predator Models During Finite and Infinite Time Intervals. Applied Mathematics and Computation, 146, 495-508.

http://dx.doi.org/10.1016/S0096-3003(02)00601-X

[8] El-Gohary, A. and Al-Ruzaiza, A. (2007) Chaos and Adaptive Control in Two Prey, One Predator System with Nonlinear Feedback. Chaos, Solitons and Fractals, 34, 443-453.

http://dx.doi.org/10.1016/j.chaos.2006.03.101

[9] El-Gohary, A. (2008) Chaos and Optimal Control of Cancer Self-Remission and Tumor System Steady States. Chaos, Solitons and Fractals, 37, 1305-1316.

http://dx.doi.org/10.1016/j.chaos.2006.10.060

[10] Al-Mahdi, A. and Khirallah, M. (2016) Bifurcation Analysis of a Model of Cancer. European Scientific Journal, 12, 68-83.

[11] Khirallah, M. and Al-Mahdi, A. (2014) Stability Analysis of a Model of Cancer Treatment by Immunotherapy. European Journal of Scientific Research, 121, 161-173.

[12] Alwan, S. and El-Gohary, A. (2011) Chaos, Estimation and Optimal Control of Habitat Destruction Model with Uncertain Parameters. Computers and Mathematics with Applications, 62, 4089-4099.

[13] El-Gohary, A. and Alwan, S. (2011) Estimation of Parameters and Optimal Control of the Genital Herpes Epidemic Model. Canadian Journal on Science and Engineering Mathematics, 2, 31-41.

[14] Satulovsky, J. E and Tome, T. (1994) Stochastic Lattice Gas Model for a Predator-Prey System. Physical Renew E, 49, 5073-5079.

[15] El-Gohary, A. and Alwan, S. (2011) Chaos and Optimal Control of Stochastic Lattice of Prey-Predator with Unknown Probabilities. International Journal of Mathematical Sciences and Applications, 1, 301-315.

[16] El-Gohary, A. (2005) Optimal Control of Stochastic Lattice of Prey-Predator Models. Applied Mathematics and Computation, 160, 15-28.

[17] van Kampen, N.G. (1981) Stochastic Processes in Physics and Chemistry. North-Holland Publishing Company Amsterdam, New York, Oxford.