Parameter Estimation of the Lotka–Volterra Model with Fractional Order Based on the Modulation Function and Its Application

Hao, Ying; Guo, Mingshun

doi:https://doi.org/10.1155/2021/6645059

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 6645059 | https://doi.org/10.1155/2021/6645059

Parameter Estimation of the Lotka–Volterra Model with Fractional Order Based on the Modulation Function and Its Application

Ying Hao^1,2and Mingshun Guo¹

Academic Editor: Yakov Strelniker

Received25 Nov 2020

Revised25 May 2021

Accepted30 Jun 2021

Published14 Jul 2021

Abstract

The Lotka–Volterra model is widely applied in various fields, and parameter estimation is important in its application. In this study, the Lotka–Volterra model with universal applicability is established by introducing the fractional order. Modulation function is multiplied by both sides of the Lotka–Volterra model, and the model is converted into linear equations with parameters to be estimated by the fractional integration method. The parameters are obtained by solving the equations. The state of the system is estimated by shifted Chebyshev polynomial. Last, the implementation program of the model is compiled. The concrete implementation method of the improved model is proposed by an example in this study.

1. Introduction

In the 1940s, Lotka (1925) and Volterra (1926) jointly laid the theoretical foundation for the Lotka–Volterra model. As an extension of the logistic model, this model was used to study the competition and cooperation between different populations in the process of biological evolution. The Lotka–Volterra model is widely applied in the financial industry, resource development, Internet of Things, currency transactions, national economy, food, and other fields [1–8].

To make the Lotka–Volterra model better adapted to applications, a stochastic Lotka–Volterra model disturbed by G-Brownian motion in the framework of nonlinear expectation was proposed by He et al. [9], which considered the uncertainty of the noise variance. Vadillo [10] studied this problem using three different stochastic models built from a single Lotka–Volterra deterministic model. More concretely, they study their mean-extinction time which satisfies the backward Kolmogorov differential equation, a linear second-order partial differential equation with variable coefficients. An extended Lotka–Volterra model was proposed by Stefan et al. [11], and they studied the microbial dynamics of two interactions in the chemostat. Tahara et al. [12] demonstrated that the Lotka–Volterra system can be stabilized by inducing small immigration of the prey or predator population, and they also checked the effect of the inclusion of the nonlinear interaction term to the stability of the Lotka–Volterra system. Stomov et al. [13] extended the existing N-species impulsive competitive models to an uncertainty case, and the main results are obtained by using Lyapunov-type functions and a comparison principle. Gavina et al. [14] built a Lotka–Volterra competition equation with a nonlinear crowding effect.

Estimating the unknown parameters of the Lotka–Volterra model is a critical issue. A Lotka–Volterra competitive model parameter estimation method based on the grey direct modeling method was proposed by Whu and Wanghai-Jun [15]. Using the least squares method, parameter estimation of a stochastic Lotka–Volterra model driven by small Levy noises from discrete observations was proposed by Wei [16]. A polynomial approximation method for the solutions of continuous-time delayed population models was presented by Yuzbasi and Karasayir [17]. They transformed nonlinear algebraic equations for the solution and proposed a technique of residual correction. A regularization method was applied to parameter inversion of the Lotka–Volterra cooperative system with the periodic boundary conditions by Xu et al. [18]. The unknown parameters of the Lotka–Volterra system can be estimated using the approximate Bayesian computation method by Skvortsov et al. [19]. Perturbed Lyapunov function methods were used to solve multispecies Lotka–Volterra models with regime switching modulated by a continuous-time Markov chain involving a small parameter in Wang et al.’s study [20]. The Lindstedt–Poincare method was applied to the Lotka–Volterra model by Amore and Fernandez [21]. They utilized an efficient systematic algorithm, there were an unprecedented number of perturbation corrections for the two dynamical variables, and the frequency was obtained.

Empirical formulas for complex dynamical processes in some fields are often expressed in the form of power-law functions, and the standard “gradient” law is not satisfied by the corresponding mechanical constitutive relation. Mechanical and physical processes with historical memory and spatial global correlation can be described by fractional differential operators concisely and accurately. In this study, a universal applicability model is established by implanting a fractional-order model into the Lotka–Volterra model. The antinoise property of the algorithm is improved by introducing integrals into the algebraic formula of the variable of the modulation function. The algebraic expression of the algorithm makes it nonasymptotically convergent. Unknown parameter estimation in the Lotka–Volterra model is accurately realized.

2. Preliminary Knowledge

2.1. Fractional Calculus

The following notations and definitions are performed before the parameter and system state estimation proof.

Mark: represents the n-dimensional Euclidean space, and represents the set of real matrices. Given the column vector , the upper corner represents the transpose. Given that , . Besides, represents an n-order diagonal matrix, and represents the determinant of matrix A. Given and , represents the Kronecker product of matrix and matrix . Interval , . Moreover, , and represent the largest integer which is less than or equal to and the smallest integer which is more than or equal to , respectively_. The square integrable space of interval is .

Definition 1. The left-sided Riemann–Liouville fractional integral of a function is defined bywhere is the gamma function.
It should be noted that is also used to refer to the -order Riemann–Liouville fractional derivative of .

Definition 2. The -order left-sided Riemann–Liouville fractional derivative of a function is defined by

Definition 3. The -order left-sided Caputo fractional derivative of a function is defined by

Definition 4. The right-sided Riemann–Liouville fractional integral of a function is defined by

Definition 5. The -order right-sided Riemann–Liouville fractional derivative of a function is defined by

Definition 6. The -order right-sided Caputo fractional derivative of a function is defined by

Lemma 1. For any interval , the following formula holds [22]:

2.2. Modulation Function

Definition 7. Given , , if , then function meets the criteria , , , and the function is called an -order modulation function on the interval .

2.3. Shifted Chebyshev Polynomial

Shifted Chebyshev polynomial is defined by the following recursive formula:

Equation (8) is equivalent to the following form:

Here, and .

Any could be represented as a shifted Chebyshev polynomial form as follows:where , . Especially, the preceding paragraph is truncated, and the following can be obtained:

3. Problem Description and Main Conclusions

3.1. Problem Description

Assume that multiple species in the same natural environment must affect the survival and development of each other. The fractional-order model is introduced to establish the general applicability Lotka–Volterra model in this study, which is as follows:where , is the fractional order, , is the system state, is the known function, refers to the initialization function of the state , and and are the system parameters to be estimated.

Furthermore, the numerical estimates of in need to be given, .

3.2. Parameter and System State Estimation

As a robust and fast estimation method, the modulating functions method has been widely used in the fields of signal processing and automatic control for parameter identification [23], derivative estimation [24, 25], and integral estimation [26]. Recently, it has been extended to state estimation problems for fractional-order systems [27–30]. However, the models considered in previous works are linear. In this study, the modulating functions method is applied for both parameter and state estimation in the nonlinear fractional-order Lotka–Volterra model. First, according to the known initial value function, the unknown parameters and of system (12) are estimated, where . On this basis, the approximate estimation in is given.

3.2.1. Parameter Estimation

Before the conclusion is given, _,. Since system (12) is known at , estimation of parameters and can be transformed into a linear system of equations about these parameters by the modulation function method. Then, can be obtained by selecting the appropriate modulation function.

Proposition 1. Given the known state of system (12), , _, and if there is a modulation function , thenThe system parameters and can be given by the following formula:where , , , , , , and .

Proof. Multiply in both sides of equation (12) and the integrals from to . Obtain and .According to Lemma 1, the left side of equation (16) becomesAccording to Definition 7, from formula (17), the following equation can be obtained:where and . Therefore, according to formula (18),Thus, Proposition 1 can be proved.

3.2.2. State Estimation

Suppose

According to the Chebyshev polynomial, can be used to estimate the corresponding system state in approximately, that is, estimates of the system state , , can be obtained by solving the algebraic equation ( can be obtained based on the known system state in ).

Proposition 2. If there are time points, then becomeswhere .
Then, the state of system (12) , can be approximately given by , and can be given by the following formula:where .

Proof. According to formula (11), the numerical solution of the system state , can be obtained only by solving the corresponding coefficient .
According to the known conditions and formula (13), the following can be obtained:The first term is used to approximate , and then, system (23) isthat is,Thus, according to formula (21), the proposition is proved.

4. Example Analysis

The survival and development of different populations in the same natural environment must influence each other. The degree of interaction among multiple populations can be analyzed quantitatively by the Lotka–Volterra model under the premise of a certain population size (n). The Lotka–Volterra model with fractional-order multicompetition relations in this study has more general application. The corresponding examples are discussed for n = 3 and n = 2, respectively.

Example 1. The Lotka–Volterra competition model is a dynamic model, and the cooperative relationship between enterprises can be reflected by this model [31]. Three companies with more invention patents in the new energy vehicle industry are chosen by this study. The number of patent holdings in the future can be predicted by the Lotka–Volterra model with the fractional order. In order to meet the requirements of the algorithm for sample data, the patent dates of enterprises are fitted to make the data continuous. Equation (12) can be written aswhere is the patent number of enterprise , is the patent number of enterprise , is the change rate of each enterprise’s own patents, and is the influence intensity of enterprise on enterprise .
The distribution of the patent holdings of the three enterprises in the past five years is given in Table 1.
According to the parameter identification (as shown in Figure 1), the number of patent granted to the Zhejiang Geely Holding Group Co., Ltd. is influenced mostly by BYD Company Limited, and the number of patent granted to Beijing New Energy Automobile Co., Ltd. is influenced least by BYD Company Limited. The number change of patent granted of the three enterprises is less affected by the enterprises themselves and more affected by other enterprises in the same system. In addition, , , , and ; it shows that in terms of patent granted number, BYD Company Limited and Zhejiang Geely Holding Group Co., Ltd. are in the predator-prey system and BYD Company Limited and Beijing New Energy Automobile Co. Ltd. are in the mutualism system.
Figure 2 shows that the growth trend of the patent holdings of BYD Company Limited has a significant slowing trend while patent holdings of Zhejiang Geely Holding Group Co., Ltd. and Beijing New Energy Automobile Co., Ltd. are in a period of rapid growth after 10 years. Patent holdings of BYD Company Limited are influenced by another two enterprises increasingly seriously. Moreover, the patent holdings of Beijing New Energy Automobile Co. Ltd. will surpass Zhejiang Geely Holding Group Co., Ltd. in five years. The gap between patent holdings of the two enterprises will widen as time goes on.

Example 2. Research on parameter estimation of patent application number of two competitive enterprises is in order to reflect the relationship between patent applied number and patent granted number more objectively. The insufficiency of decision condition in Example 1 is supplemented. Equation (12) can be written aswhere is the patent applied number of enterprise , is the patent applied number of enterprise , is the patent application change rate of each enterprise, and is the influence intensity of enterprise on enterprise .
Patents applied number of the two enterprises from 2014 to 2018 is given in Table 2.
According to the parameter identification (as shown in Figure 3), in terms of the patent application number, BYD Company Limited and Zhejiang Geely Holding Group Co., Ltd. are in the predator-prey system. The change of the patent applied number of Zhejiang Geely Holding Group Co., Ltd. is influenced mostly by BYD Company Limited, while the influence of it owns much less. It shows that the innovation ability of Zhejiang Geely Holding Group Co., Ltd. is slightly lower than BYD Company Limited, and BYD Company Limited has a leading edge in technology.

5. Conclusion

In this study, the modulating functions method was applied to estimate unknown parameters, and the state of the Lotka–Volterra model was defined by a fractional-order nonlinear system. First, by multiplying the modulating functions to the considered model and applying fractional-order integration by parts formulas, the differential equations can be transformed into a set of integral equations where the unknown parameters are linear concerning the integrals. Then, the unknowns are estimated by solving the linear system of the obtained integral equations. Second, by taking the truncated expansion of the state based on the shifted Chebyshev polynomials, the unknown coefficients are obtained by solving a linear system using the initialization function of the state. Then, these estimated coefficients are used to predict the future behavior of the state.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Humanities and Social Sciences Fund of Department of Education of Liaoning Province, China (LGZY2019004), the Special Project of Serving Local Economic and Social Development Decision-Making of Shenyang University of Technology (FWDFGD2020010), and the Humanities and Social Sciences Fund of Liaoning Engineering Vocational College, China (ZYS201901).

References

S. J. Lee, D. J. Lee, and H. S. Oh, “Technological forecasting at the Korean stock market: a dynamic competition analysis using Lotka-Volterra model,” Technological Forecasting and Social Change, vol. 72, pp. 1004–1057, 2015.
View at: Google Scholar
L. Lambertini and G. Leitmann, “On the attainment of the maximum sustainable yield in the Verhulst-Lotka-Volterra model,” Automatica, vol. 110, Article ID 108555, 2019.
View at: Publisher Site | Google Scholar
R. Casado-Vara, A. Canal-Alonso, M. D. Rey et al., “Smart buildings ioT networks accuracy evolution prediction to improve their reliability using a Lotka-Volterra Ecosystem Model,” Sensors, vol. 19, no. 21, Article ID 4642, 2019.
View at: Publisher Site | Google Scholar
P. Gatabazi, J. C. Mba, and E. Pindza, “Modeling cryptocurrencies transaction counts using variable-order fractional grey Lotka-Volterra dynamical system,” Chaos, Solitons and Fractals, vol. 127, pp. 283–290, 2019.
View at: Publisher Site | Google Scholar
I. Ivanova, O. Strand, and L. Leydesdorff, “An eco-systems approach to constructing economic complexity measures: endogenization of the technological dimension using Lotka-Volterra equations,” Advances in Complex Systems, vol. 22, no. 1, Article ID 1850023, 2019.
View at: Publisher Site | Google Scholar
L. Shi, Y. Xu, L. Zhan et al., “Enhancing vitamin B12 content in co-fermented soy-milk via a Lotka Volterra model,” Turkish Journal of Biochemistry, vol. 43, no. 6, pp. 671–678, 2018.
View at: Publisher Site | Google Scholar
Z. Cao, Q. Zhang, and Y. Zhao, “The principle of competitive exclusion about a stochastic lotka-volterra model with two predators competing for one prey,” Discrete Dynamics in Nature and Society, vol. 2018, Article ID 7312581, 13 pages, 2018.
View at: Publisher Site | Google Scholar
H.-C. Hung, Y.-C. Chiu, H.-C. Huang et al., “An enhanced application of Lotka-Volterra model to forecast the sales of two competing retail formats,” Computers and Industrial Engineering, vol. 109, pp. 325–334, 2017.
View at: Publisher Site | Google Scholar
P. He, Y. Ren, and D. Zhang, “Asymptotic moment estimation for stochastic Lotka–Volterra model driven by G -Brownian motion,” Stochastics An International Journal of Probability and Stochastic Processes, vol. 2020, Article ID 3705325, 18 pages, 2020.
View at: Publisher Site | Google Scholar
F. Vadillo, “Comparing stochastic Lotka-Volterra predator-prey models,” Applied Mathematics and Computation, vol. 360, pp. 181–189, 2019.
View at: Publisher Site | Google Scholar
V. Stefan, D. B. Sophie, F. Karoline et al., “Bistability in a system of two species interacting through mutualism as well as competition: chemostat vs. Lotka-Volterra equations,” PLoS One, vol. 13, no. 6, Article ID e0197462, 2018.
View at: Publisher Site | Google Scholar
T. Tahara, M. K. A. Gavina, T. Kawano et al., “Asymptotic stability of a modified Lotka-Volterra model with small immigrations,” Entific Reports, vol. 8, Article ID 7029, 2018.
View at: Publisher Site | Google Scholar
G. T. Stamov, S. Simeonov, and I. M. Stamova, “Uncertain impulsive Lotka-Volterra competitive systems: robust stability of almost periodic solutions,” Chaos, Solitons and Fractals: Applications in Science and Engineering: An Interdisciplinary Journal of Nonlinear Science, vol. 110, pp. 178–184, 2018.
View at: Publisher Site | Google Scholar
M. K. A. Gavina, T. Tahara, K. I. Tainaka et al., “Multi-species coexistence in Lotka-Volterra competitive systems with crowding effects,” Scientific Reports, vol. 8, Article ID 1198, 2018.
View at: Publisher Site | Google Scholar
L. F. Wu and Y. N. WangHai-jun, “Estimation the parameters of Lotka-Volterra model based on the grey direct modeling and its application,” Expert Systems with Applications, vol. 38, no. 6, pp. 6412–6416, 2011.
View at: Publisher Site | Google Scholar
C. Wei, “Parameter estimation for stochastic Lotka-Volterra model driven by small Levy noises from discrete observations,” Communication in Statistics-Theory and Methods, vol. 12, pp. 1–10, 2020.
View at: Publisher Site | Google Scholar
S. Yuzbasi and M. Karacayir, “A numerical method for solutions of Lotka-Volterra predator-prey model with time-delay,” International Journal of Biomathematics, vol. 11, no. 2, Article ID 1850028, 2018.
View at: Publisher Site | Google Scholar
L. Xu, J. Y. Liu, and G. Zhang, “Pattern formation and parameter inversion for a discrete Lotka-Volterra cooperative system,” Chaos Solitons and Fractals, vol. 110, pp. 226–231, 2018.
View at: Publisher Site | Google Scholar
A. Skvortsov, B. Ristic, and A. Kamenev, “Predicting population extinction from early observations of the Lotka-Volterra system,” Applied Mathematics and Computation, vol. 320, pp. 371–379, 2018.
View at: Publisher Site | Google Scholar
R. Wang, X. Y. Li, and G. Yin, “Asymptotic properties of multi-species lotka-volterra models with regime switching involving weak and strong interactions,” Journal of Nonlinear Science, vol. 30, no. 2, pp. 565–601, 2020.
View at: Publisher Site | Google Scholar
P. Amore and F. M. Fernandez, “On the application of the Lindstedt-Poincare method to the Lotka-Volterra system,” Annals of Physics, vol. 396, pp. 293–303, 2018.
View at: Publisher Site | Google Scholar
Y. Q. Podlubny and Y. Chen, “Adjoint fractional differential expressions and operators,” in Proceedings of ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference IDETC/CIE, pp. 4–7, Las Vegas, NV, USA, September 2007.
View at: Publisher Site | Google Scholar
D. Y. Liu, T. M. Laleg-Kirati, O. Gibaru, and W. Perruquetti, “Identification of fractional order systems using modulating functions method,” in Proceedings of the 2013 American Control Conference, pp. 17–19, Washington, DC, USA, June 2013.
View at: Publisher Site | Google Scholar
D. Y. Liu, T. M. Laleg-Kirati, W. Perruquetti, and O. Gibaru, “Non-asymptotic state estimation for a class of linear time-varying systems with unknown inputs,” in Proceedings of the 19th World Congress of the International Federation of Automatic Control, pp. 324–328, Cape Town, South Africa, August 2014.
View at: Publisher Site | Google Scholar
D. Y. Liu and T. M. Laleg-Kirati, “Robust fractional order differentiators using generalized modulating functions method,” Signal Processing, vol. 107, pp. 395–406, 2015.
View at: Publisher Site | Google Scholar
Y. Tian, Y. Q. Wei, D. Y. Liu, and D. Boutat, “Fast and robust estimation for positions and velocities from noisy accelerations using generalized modulating functions method,” Mechanical Systems and Signal Processing, vol. 133, Article ID 106270, 2019.
View at: Publisher Site | Google Scholar
X. Wei, D. Y. Liu, and D. Boutat, “Non-asymptotic state estimation for a class of fractional order linear systems,” IEEE Transactions on Automatic Control, vol. 62, no. 3, pp. 1150–1164, 2016.
View at: Google Scholar
Y. Q. Wei, D. Y. Liu, D. Boutat, and Y. M. Chen, “An improved pseudo-state estimator for a class of commensurate fractional order linear systems based on fractional order modulating functions,” Systems and Control Letters, vol. 118, pp. 29–34, 2018.
View at: Publisher Site | Google Scholar
Y. Q. Wei, D. Y. Liu, and D. Boutat, “Innovative fractional derivative estimation of the pseudo-state for a class of fractional order linear systems,” Automatica, vol. 99, pp. 157–166, 2019.
View at: Publisher Site | Google Scholar
Y. Q. Wei, D. Y. Liu, D. Boutat, H. R Liu, and C. Lv, “Modulating functions based differentiator of the pseudo-state for a class of fractional order linear systems,” Journal of Computational and Applied Mathematics, vol. 384, Article ID 113161, 2021.
View at: Google Scholar
D. Guo, H. Zhang, C. Zheng et al., “Analysis of the future development of Chinese auto energy saving and environmental benefits,” Systems Engineering-- Theory&Practice, vol. 36, no. 6, pp. 1593–1599, 2016.
View at: Google Scholar

Copyright

Copyright © 2021 Ying Hao and Mingshun Guo. 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.

PDF Download Citation

Download other formats

Order printed copies