Psi-Caputo Logistic Population Growth Model

Awadalla, Muath; Yameni Noupoue, Yves Yannick; Asbeh, Kinda Abu

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

Journal of Mathematics

On this page

Abstract Introduction Preliminaries Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

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

Psi-Caputo Logistic Population Growth Model

Muath Awadalla,¹Yves Yannick Yameni Noupoue,²and Kinda Abu Asbeh¹

Academic Editor: Nan-Jing Huang

Received14 Apr 2021

Accepted06 Jul 2021

Published27 Jul 2021

Abstract

This article studies modeling of a population growth by logistic equation when the population carrying capacity K tends to infinity. Results are obtained using fractional calculus theories. A fractional derivative known as psi-Caputo plays a substantial role in the study. We proved existence and uniqueness of the solution to the problem using the psi-Caputo fractional derivative. The Chinese population, whose carrying capacity, K, tends to infinity, is used as evidence to prove that the proposed approach is appropriate and performs better than the usual logistic growth equation for a population with a large carrying capacity. A psi-Caputo logistic model with the kernel function performed the best as it minimized the error rate to 3.20% with a fractional order of derivative = 1.6455.

1. Introduction

Real-life phenomena are often modeled by the mean of differential equations. Differential equations are built using two mains tools which are derivatives and integrals. In the early stage of differential equations, only integer numbers were used as order of a derivative or integral. However, three centuries ago, a new type of derivative called fractional derivative was introduced [1]. This is a new form of derivative and integral tries to generalize the use of any real number or fraction as order of derivative or integral which was the genesis of fractional calculus. In the earlier days following its introduction, fractional calculus was mainly investigated theoretically. Hence, early research works were focused designing fractional models and then proving existence and uniqueness of their solutions. Samples such theoretical studies are found in [2–6] and the references therein. In all what follows, statements as integer order of derivation, classical approach, and classical differential equations are often used to refer to the model in which order of derivatives or integrals are all integers.

In recent decades, fractional calculus has taken a new direction. Indeed, after researchers have well established the theory of fractional differential equations, many are now focusing on their application to solve real-life problems. It has since been proven in many works that the fractional differential equations can be successfully used for modeling some natural phenomena even when the classical differential equations fail to do so. Moreover, in the case where both classical and fractional approaches are applicable to solve a problem, it is common to observe that the fractional approach will minimize the fitting error compared to the classical approach. Fractional differential equations have been successfully used for modeling in finance [7] and epidemiology [8, 9], a few examples just to mention. More generally, samples of works in which fractional calculus has successfully been used for modeling are found in [10–17] and the references therein.

This work contributes in showing the power of fractional differential equations over the classical differential equations in modeling. A fractional model is built, its existence and uniqueness are proven, and the numerical study is carried out with the idea to support theoretical results. More generally, we referred to the model built in this work as the fractional logistic equation. Recent studies on finding solution of fractional logistic equation are found in [18, 19]. The population growth modeling is considered when the carrying capacity is very large. The Chinese population which has such property was considered in simulation. It appeared that the -Caputo model introduced by Almeida [20], with recent development found in [21], is a good tool for solving this problem. Moreover, square root function used as kernel better fits the Chinese population more than the classical logistic population growth model. The work is organized as follows. Preliminaries definitions of fractional calculus are given in Section 2. The logistic growth model is defined in Section 3. The -Caputo is defined in Section 4. The approximation of population with large carrying capacity is given in Section 5, and finally, a simulation with Chinese population is done.

2. Preliminaries

The theory of fractional calculus is built upon definitions, theorems, and special notations. Moreover, fractional integrals and fractional derivatives are useful tools for designing and solving fractional differential equations. In this section, selected definitions, lemma, and theorems which are useful for the remaining of the work are given.

Definition 1. (see [22]). The Riemann–Liouville fractional integral of order of a function is defined asEquation (1) provided that the right-hand side of the integral is pointwise defined on ; moreover, is the well-known gamma function defined by .

Definition 2. (see [22]). The Caputo derivative of order of a function is defined as where , is the integer part of.

Definition 3. (see [20]). Let and and be an increasing function with , ; then, denotes the fractional integral of with respect to , and it is given byThis definition is known as -Caputo integral if the fractional integral is taken in the Caputo sense. It is used in sequel for formulating the solution to the -Caputo fractional model.

Definition 4. (see [20]). Let and , with be an increasing function and ; then, denotes the fractional derivative of with respect to (-Caputo if the fractional derivative is taken in the Caputo sense), and it is given by

Lemma 1. (see [20]). Let , a natural number, n, such that if , ; then,

Lemma 2. (see [23]). Error rate of estimate.

Let us consider an experimental data set represented by a vector of size n, obtained from a time- or space-dependent process such that the i^th observation is denoted by . If the experimental is fitted by a function , in a way that the i^th fitted data value is denoted by , the root-mean-squared deviation between the original data set and the fitted data set is computed by

On the contrary, the norm of the original data vector is obtained by averaging the sum of squared elements and then computing the square root of the results. This is expressed by the following formula:

Dividing equation (6) by equation (7) produces an interesting tool to evaluate the fitting error which is represented as

The rate ER evaluates the magnitude of the error that occurs when the estimated values are used in place of the original values. It has an advantage over commonly used root square error and mean square error because it standardizes the magnitude of the data set used.

3. Logistic Population Growth Model

The theory of the population growth modeling using logistic equation was introduced by an economist named Malthus [24–26]. In fact, given an initial population with growth and reproductivity capacity, the theoretical expectation would be that the population size will approach infinity as the time increases indefinitely. However, Malthus has proven in his theory that a population growth tends to stabilize at some point when the time approaches infinity. He then proposed a model of population growth, called the logistic growth model, which is defined aswhere r is the growth rate and K is the carrying capacity, which represents the maximum value that the population size may reach. Hence, the population size stabilizes when the carrying capacity is reached. A general solution of the classical logistic model defined by equation (9) iswhere N₀=N (0) is the initial size of the population at t = 0.

Since we aim to prove that the fractional differential equation is better than the classical differential equation, the fractional differential equation equivalent to problem defined by equation (9) is built as follows:or simply

The following lemma introduces the general formula of the solution to the fractional logistic equation with initial condition, defined by equation (11).

Lemma 3. Assume is an integrable function defined on ; then, the general solution of the fractional differential equation given by equation (12) is equivalent to the following integral equation:

It is important to highlight the use of the -Caputo fractional integral, equation (3) in Lemma 3.

Proof. Applying the operator to both sides of equation (12) leads to .

4. Main Result of Psi-Caputo Logistic Population Growth Model

In this section, the theoretical study is carried around equation (12) in order to prove existence and uniqueness of a solution.

Given , let denote the Banach space of all continuous functions from to , endowed with the norm defined by .

An operator associated with the problem defined by equation (12) can be defined as

The existence and uniqueness of solution for the problem given by equation (13) is discussed in sequel. Prior to establish proof of the main results, the following statements are established. These are mathematical hypotheses that are used in the proof sections.(A1) is continuous(A2)There exists such that (A3)There exists a function and a nondecreasing function such that (A4)There exists a constant such that

The following theorem establishes and proves existence of at least one solution to equation (12).

Theorem 1. Assume that (A1), (A3), and (A4) hold. Then, equation (12) has at least one solution in the interval .

Proof. The proof will be split into several steps. In the first step, we show that the operator maps bounded sets into bounded sets of . Let be a bounded set in ; then,Applying the norm, on both sides of equation (15) leads toThe next step in the proof is to show that the operator maps bounded sets into equi-continuous sets of .
Let with and . The following relation holds because of the above assumptions:where the right-hand side of equation (17) tends to zero as . That is, as .
Note that the right-hand side of equation (17) is independent of ; therefore, by Arzela–Ascoli theorem, we conclude that is completely continuous.
The last step to complete the assumptions of Leray–Schauder nonlinear alternative theorem is to show the boundedness of the set of all solution to equation . Assume that is a solution equation (12); then, it follows from equation (13) thatInverting both sides of equation (18) and dividing them by leads to the following relation:Recalling (A4), there exists a constant , which is indeed such that . Moreover, let us construct the set . It is obvious that the operator is continuous and completely continuous. Based on the constructed , such that for some . Consequently, by the nonlinear alternative of Leray–Schauder type, we deduce that has a fixed point which is a solution to the problem defined by equation (13).

Theorem 2. Let us assume that (A1) and (A2) hold. If , then the problem defined in equation (12) has a unique solution on .

Proof. Consider the operator defined in equation (14) and define a ball
withwhere .
First, let us show that . It is done as follows. For any and using equation (14), we have the following relation:On the contrary, triangular inequality is applied on and produced in the next relation:Substituting the appropriate fragment of equation in equation (21) by equation (22) leads to a new relation defined below:Equation (23) is enough to conclude that .
The next step of the proof is to show that the operator is a contraction. , and the following relation holds:From equation (24), is said to be a contraction. By the Banach contraction mapping theorem, the fractional differential equation, equation (12), has a unique solution on .

5. Population Growth with Carrying Capacity Approaching Infinity

The carrying capacity of a population is computed with accuracy using nonlinear optimization routine on the existing sample data over a given time interval. From equation (9), it is obvious that the carrying capacity has a high impact on the model. Some researchers focused in the past on the population growth behavior based on the carrying capacity dynamic (see [24, 27, 28]). Considering the human population growth modeling, we will observe that there exist some countries that try to control their population growth. China is a country with large population for which the concept of “large carrying capacity,” and is applicable. Birth control in China has reduced the population growth rate. In this regard, the true population size is usually smaller than what is theoretically expected using forecasting tools. This leads to a population size which is far away smaller than the carrying capacity. Mathematically, this is denoted as

That is,

Under the conditions defined by equation (25) and equation (26), the initial logistic model equation (9) becomes an exponential model:

A general solution of the exponential model defined by equation (27) is

The solution to the exponential model defined by equation (28) shows that the population size will increase infinitely as the time approaches infinity. Such result means that the Chinese population, for instance, will grow infinitely. This coincides with neither the Malthusian theory nor what is practically known about the Chinese population. In what follows, an alternative approach is proposed for the modeling of large size population with carrying capacity . Recalling the -derivative from Definition 4, the model defined by equation (27) can be written in fractional sense as

Applying the fractional integral operator to both sides of equation (29) and using Theorem 7.2 from [29] is enough to find the solution of the model defined by equation (29) aswhere is the Mittag-Leffler function.

In simulations, a choice of the Kernels is made based on their shape in a way to best fit the data.

6. Simulation Results

In this section, we proved based on the numerical study with the Chinese population that the proposed method in Section 5 is the most appropriate for fitting a population with large carrying capacity than the usual logistic and exponential models. Historical data of the Chinese population was retrieved from the World Bank website [30]. A nonlinear optimization routine “lsqcurvefit” from Matlab was used to estimate the parameter that would best fit the classical logistic model equation (9), the exponential model equation (27), and the psi-Caputo fractional model equation (30). Recalling Lemma 2, in which error metric is proposed, the following results were obtained. Physical meaning of parameters used in this section is given prior to presenting all the results’ scenarios. (i)Carrying capacity, K is the maximum size that a population can reach without causing any degradation(ii)The growth rate r of a population is the percentage change in size of a given population in unit period (in this case per year)(iii)Constant α is the fractional order of derivative involved in the model design

Figure 1 depicts the original data alongside fitting line of the estimation using the classical logistic approach. The best parameters values for this model were found as the growth rate of r = 0.0145, which means over the period from 1960 to 2020 the logistic model revealed that the overall growth rate of Chinese population was less than 2%. Moreover, the carrying capacity, K = 3.3198 × 10²³, is theoretically the maximum size that Chinese population could reach. Using the classical logistic model, the original data set would be fitted while making the total error rate of ER = 6.27%.

Figure 2 represents the solution obtained fitting original data set with classical exponential. Fitted model and original data are both shown in the figure. This approach has produced similar result as the classical logistic approach. In fact, the best estimated parameters are r = 0.0145 and K = 3.3198 × 10²³ for a minimum error rate of ER = 6.27%. At this point, a partial result concerning large carrying capacity is obtained. In fact, for a large carrying capacity, the logistic growth model and exponential growth model coincide.

Figure 3 represents the result obtained when using the -Caputo approach method with the kernel function being to fit the original data. The best parameters of the model were estimated as growth rate r = 0.0950, K = 3.3198 × 10²³ and fractional order of derivative α = 3.0096. Fitting true data using those parameters produced a total error of ER = 1.57%. This model is nonrealistic and irrelevant in this study since one expects the order of derivative, .

Figure 4 shows true data alongside fitted data using the -Caputo model with as the kernel function. The best parameters estimated for this model are the growth rate r = 0.0484 and the fractional order of derivative = 1.6455. The use of those parameters produced a total error of ER = 3.20%.

7. Conclusion

In this work, we studied population growth modeling using fractional calculus. We designed a fractional model and proved its existence and uniqueness, and then, we simulated the results using Chinese population data. At the end of the study, we observed the following. It can be cumbersome to give a practical or physical meaning to a fractional order of derivative associated to a model’s solution. However, it is still important to provide an explanation for what is obtained in a modeling process. Let us, for instance, state that, in a more general physical sense, the first derivative of a trajectory equation is the speed or velocity. For the population growth model, the first derivative represents the speed by which the said population grows. It is also important to give a meaning to a fractional order of derivative when used in fractional differential equations. In this regard, one can say that, given a modeling problem for which there exists a solution based on classical approach, if this solution involves a first-order derivative, then one would expect from its fractional counterpart to be optimal with a fractional derivative order value such that . In fact, if the fractional approach to a problem involving a first-order differential equation is optimal when , this coincides with the classical solution. Hence, the fractional approach in such case does not improve anything compared to the classical solution. However, if the solution is such that , then one can proudly say that the fractional differential approach is better compared to the classical approach for this problem.

In this study, when the carrying capacity of population is very large and that the said population is subjected to a restricted growth, then the logistic approach and the exponential approach coincide. However, they might not be suitable for modeling the growth of such population. The -Caputo method in both cases minimized the error. Moreover, the kernel function is the most appropriate in this study because it best fits the data with an , which is in line with the solution obtained using the classical differential equation approach, for which the solution is based on a first-order derivative.

An expansion of this work in future could consist of investigating the results using other types of fractional derivatives and integrals such as Hadamard and -Hilfer and compare the results with those obtained so far.

Data Availability

The numerical data used to support the findings of this study have been deposited in the World Bank repository (https://data.worldbank.org/indicator/SP.POP.TOTL?locations=CN).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Each of the authors, M.A, Y.Y.Y, and K.A. contributed to each part of this work equally and read and approved the final version of the manuscript.

Acknowledgments

This study was funded by the Deanship of Scientific Research (DSR), King Faisal University (KFU), Ahsa, Saudi Arabia. The authors, therefore, acknowledge technical and financial support under NASHER track with Grant no. 206151 of DSR at KFU.

References

L. Debnath, “A brief historical introduction to fractional calculus,” International Journal of Mathematical Education in Science and Technology, vol. 35, no. 4, pp. 487–501, 2004.
View at: Publisher Site | Google Scholar
N. I. Mahmudov, M. Awadalla, and K. Abuassba, “Nonlinear sequential fractional differential equations with nonlocal boundary conditions,” Advances in Difference Equations, vol. 2017, no. 1, p. 319, 2017.
View at: Publisher Site | Google Scholar
N. I. Mahmudov, M. Awadalla, and K. Abuassba, “Hadamard and caputo-hadamard FDE's with three point integral boundary conditions,” Nonlinear Analysis and Differential Equations, vol. 5, pp. 271–282, 2017.
View at: Publisher Site | Google Scholar
C. Li and S. Sarwar, “Existence and contiuation of solution for Caputo type fractional differential equations,” Electronic Journal of Differential Equations, vol. 2016, pp. 1–14, 2016.
View at: Google Scholar
A. Babakhani and V. Daftardar-Gejji, “Existence of positive solutions of nonlinear fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 278, no. 2, pp. 434–442, 2003.
View at: Publisher Site | Google Scholar
X. Zhang, X. Zhang, Z. Liu, W. Ding, H. Cao, and T. Shu, “On the general solution of impulsive systems with hadamard fractional derivatives,” Mathematical Problems in Engineering, vol. 2016, Article ID 2814310, 12 pages, 2016.
View at: Publisher Site | Google Scholar
M. A. Dokuyucu, “Analysis of a novel finance chaotic model via ABC fractional derivative,” Numerical Methods for Partial Differential Equations, Article ID 22598, 2020.
View at: Publisher Site | Google Scholar
M. A. Dokuyucu, E. Celik, H. Bulut, and H. M. Baskonus, “Cancer treatment model with the Caputo-Fabrizio fractional derivative,” The European Physical Journal Plus, vol. 133, no. 3, pp. 1–6, 2018.
View at: Publisher Site | Google Scholar
K. O. C. A. İlknur, E. Akçetin, and P. Yaprakdal, “Numerical approximation for the spread of SIQR model with caputo fractional order derivative,” Turkish Journal of Science, vol. 5, no. 2, pp. 124–139, 2020.
View at: Google Scholar
H. Sun, Y. Zhang, D. Baleanu, W. Chen, and Y. Chen, “A new collection of real world applications of fractional calculus in science and engineering,” Communications in Nonlinear Science and Numerical Simulation, vol. 64, pp. 213–231, 2018.
View at: Publisher Site | Google Scholar
I. Podlubny, “Fractional differential equations,” in Mathematics in Science and Engineering, vol. 198, Academic Press, San Diego, CA, USA, 1999.
View at: Google Scholar
I. Area and H. LosadaJ. Fi et al., “On a fractional order Ebola epidemic model,” Advances in Difference Equations, vol. 2015, no. 1, p. 278, 2015.
View at: Publisher Site | Google Scholar
Y. Sofuoglu and N. Ozalp, “A fractional order model on bilingualism,” Communications Faculty of Sciences University of Ankara Series A1, vol. 2014, pp. 81–89, 2014.
View at: Google Scholar
Y. Sofuoglu and N. Ozalp, “Fractional order bilingualism model without conversion from dominant unilingual group to bilingual group,” Differential Equations and Dynamical Systems, vol. 25, no. 1, pp. 1–9, 2015.
View at: Publisher Site | Google Scholar
R. Hilfer, Applications of Fractional Calculus in Physics, World Scientiﬁc, Hackensack, NJ, USA, 2001.
F. Mainardi, “Fractional relaxation-oscillation and fractional diffusion-wave phenomena,” Chaos, Solitons & Fractals, vol. 7, no. 9, pp. 1461–1477, 1996.
View at: Publisher Site | Google Scholar
R. Almeida, N. R. O. Bastos, and M. T. T. Monteiro, “Modeling some real phenomena by fractional differential equationsfferential. Equations,” Mathematical Methods in the Applied Sciences, vol. 39, no. 16, pp. 4846–4855, 2015.
View at: Publisher Site | Google Scholar
I. Area and J. J. Nieto, “Power series solution of the fractional logistic equation,” Physica A: Statistical Mechanics and Its Applications, vol. 573, Article ID 125947, 2021.
View at: Publisher Site | Google Scholar
L. N. Kaharuddin, C. Phang, and S. S. Jamaian, “Solution to the fractional logistic equation by modified Eulerian numbers,” The European Physical Journal Plus, vol. 135, no. 2, p. 229, 2020.
View at: Publisher Site | Google Scholar
R. Almeida, “A caputo fractional derivative of a function with respect to another function,” Communications in Nonlinear Science and Numerical Simulation, vol. 44, pp. 460–481, 2017.
View at: Publisher Site | Google Scholar
R. Almeida, “Functional differential equations involving the ψ-caputo fractional derivative,” Fractal and Fractional, vol. 4, no. 2, p. 29, 2020.
View at: Publisher Site | Google Scholar
S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives. Theory and Applications, Gordon and Breach, Yverdon, Switzerland, 1993.
Y. Y. Y. Noupoue, Y. Tandoğdu, and M. Awadalla, “On numerical techniques for solving the fractional logistic differential equation,” Advances in Difference Equations, no. 1, pp. 1–13, 2019.
View at: Publisher Site | Google Scholar
T. R. Malthus and J. Pullen, TR Malthus: Principles of Political Economy, vol. 2, Cambridge University Press, Cambridge, CA, USA, 1989.
T. R. Malthus, An Essay on the Principle of Population, Reeves & Turner, London, UK, 1872.
T. R. Malthus, Population: The First Essay, University of Michigan Press, Michigan, MI, USA, 1959.
J. A. Brander and M. S. Taylor, “The simple economics of Easter Island: a Ricardo-Malthus model of renewable resource use,” American Economic Review, vol. 88, pp. 119–138, 1998.
View at: Google Scholar
C. Bianca and L. Guerrini, “On the Dalgaard-Strulik model with logistic population growth rate and delayed-carrying capacity,” Acta Applicandae Mathematicae, vol. 128, no. 1, pp. 39–48, 2013.
View at: Publisher Site | Google Scholar
K. Diethelm, The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type, Springer Science & Business Media, Berlin, Germany, 2010.
Population, total. (n.d.). Retrieved from https://data.worldbank.org/indicator/SP.POP.TOTL?locations=CN.

Copyright

Copyright © 2021 Muath Awadalla 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.

PDF Download Citation

Download other formats

Order printed copies