Complexity

Volume 2019, Article ID 3426974, 14 pages

https://doi.org/10.1155/2019/3426974

## Numerical Simulation of a Class of Three-Dimensional Kolmogorov Model with Chaotic Dynamic Behavior by Using Barycentric Interpolation Collocation Method

^{1}Institute of Computer Information Management, Inner Mongolia University of Finance and Economics, Hohhot 010070, China^{2}Department of Mathematics, Inner Mongolia University of Technology, Hohhot 010051, China^{3}Institute of Economics and Management, Jining Normal University, Jining 012000, Inner Mongolia, China

Correspondence should be addressed to Wei Zhang; moc.361@wzyxfsnj

Received 5 January 2019; Revised 13 February 2019; Accepted 3 March 2019; Published 9 April 2019

Academic Editor: Giacomo Innocenti

Copyright © 2019 Mingjing Du et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper numerically simulates three-dimensional Kolmogorov model with chaotic dynamic behavior by barycentric Lagrange interpolation collocation method. Some numerical examples are studied for finding some new chaotic behaviors and demonstrating some existing chaotic dynamic behaviors of the Kolmogorov model. Results obtained by the present method indicate that the method has merits of small operations and good numerical stability.

#### 1. Introduction

J.P. Berrut [1, 2] introduced barycentric Lagrange interpolation and studied its numerical stability and convergence. Barycentric Lagrange interpolation is unconditionally stable at the Chebyshev points. S.P. Li, Z.Q. Wang [3, 4] gave some algorithms of barycentric Lagrange interpolation collocation method (BLICM). Some authors [3, 5–8] solved all sorts of equations and showed the BLICM has merits of small operations and high precision (see [3, 4, 9]). This paper numerically simulates some three-dimensional Kolmogorov models with chaotic dynamic behavior. The purpose of this paper is to find some new chaotic behaviors and verify the existing chaotic dynamic behaviors by the BLICM.

Three-dimensional Kolmogorov models comprise a significant class of ecological models that are used widely in ecology to represent the dynamic behavior of prey and predators, which are expressed in the following form:where represents the population density of the species and represents the per capita growth rate of the species.

In ecology, the most frequently used model is the Lotka-Volterra system; that is, each per capita growth function is affine and chosen as the logistic growth. In this circumstance, model (1) reads as

Model (2) is a totally competitive system if all parameters are positive.

There are several famous functional responses in the Kolmogorov model which are referred to as Holling type I, type II, type III, type IV, Monod-Haldane type, Hassel-Verley type, Beddington-DeAngelis type functional response, etc.

In this paper, we consider the following three-dimensional Kolmogorov model with functional response:with the following initial condition:where represent the population density of the species and is the known functional response. are known constants. are unknown functions of time .

Many researchers [10–20] studied the dynamics of three dimensional Kolmogorov model with different type of functional response in theory. They found some chaotic dynamics [13, 18–23] of the three-dimensional Kolmogorov model. Chaos and hyperchaos exist in many natural processes and are one of the main contents of nonlinear science research. Although many kinds of numerical methods of the Kolmogorov model have been announced, simple and efficient methods have always been the direction that scholars strive to pursue. This paper suggests the BLICM to solve the three-dimensional Kolmogorov model. Model (3) is adopted as an example to elucidate the solution process.

#### 2. Barycentric Lagrange Interpolation Collocation Method

First of all, we discrete computational interval by Chebyshev points into and construct following linear iterative format of model (3):

The format (5) is convergent; then , , .

Next, we transform format (5) into the following linear algebraic equations.where is order matrix. is order unit matrix and is a symbol of diagonal matrix composed of vectors. is, respectively, barycentric interpolation primary function and is center of gravity Lagrange interpolation weight. The vector

The vector

The last and the first line of equations (6) are replaced separately by the equation of the initial condition (9) in turn.

So, we can get a numerical solution of (3) and (4).

#### 3. Numerical Experiments

In this section, some numerical examples are studied to demonstrate the accuracy of the present method. The examples are computed using MatlabR2017a. In numerical experiments, the number of nodes . The accuracy of iteration control and the initial iteration value .

*Experiment 1. *We consider the following three-species food chain model [18]:where is real parameters, which satisfy the initial condition Results of Experiment 1 are given in Figures 1-2.

Figure 1 is obtained by using the current method with . Among them, is the time series plot; is the phase diagram of ; is the phase diagram of ; is the graph projected on -plane; is the three-dimensional space graph. Figure 2 is obtained by using the current method with . is the time series plot; is the phase diagram of ; is the three-dimensional space graph; is the graph projected on -plane; is the the graph projected on -plane; is the graph projected on -plane.

Figure 1 gets some new chaotic behaviors. Figure 2 verifies the existing chaotic dynamic behaviors [18]. Our study suggests that model (10) will go chaotic when the rate of the self-reproduction of the prey is large. Our numerical results are in good agreement with the theory [18].