Next Article in Journal
Fractional-Order Modeling and Nonlinear Dynamic Analysis of Forward Converter
Previous Article in Journal
Fractal and Multifractal Analysis of Microscopic Pore Structure of UHPC Matrix Modified with Nano Silica
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of a Fractional-Order Time-Delayed Solow Growth Model with Environmental Pollution

1
School of Applied Science, Beijing Information Science and Technology University, Beijing 100192, China
2
School of Statistics and Mathematics, Central University of Finance and Economics, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(6), 361; https://doi.org/10.3390/fractalfract8060361
Submission received: 6 May 2024 / Revised: 12 June 2024 / Accepted: 14 June 2024 / Published: 18 June 2024
(This article belongs to the Section General Mathematics, Analysis)

Abstract

:
Economic growth is resulting in serious environmental problems. Effectively establishing an economic growth model that considers environmental pollution is an important topic. To analyze the interplay between economic growth and environmental pollution, a fractional-order time-delayed economic growth model with environmental purification is proposed in this paper. The established model considers not only the environment and economic production but also the labor force population and total factor productivity. Furthermore, the asymptotic stability conditions and parameter stability interval are provided. Finally, in numerical experiments, the correctness of the theory is verified.

1. Introduction

Fractional calculus theory analysis and applications have become topics of great interest in current scientific research. The wider application background of fractional calculus has attracted the attention of many scholars from various fields, resulting in abundant research results [1,2,3,4]. In the field of finance, because of the “memory effect” of fractional calculus, fractional order equations can describe the long-term logarithmic prices of some financial assets well [5]. Compared with integer calculus, the main advantage of fractional calculus is its memory, and it has been proven to be a very suitable tool for describing the memory and genetic characteristics of various materials and processes. Financial and economic variables have a longer memory, and therefore, it is more appropriate to use fractional differential equation models to portray the dynamic behavior of economic systems, such as, exchange rates, gross domestic product (GDP), interest rates, production, and stock market prices, which are changing in terms of the financial and economic system. This provides a scientific approach to predict economic growth.
The study of the complex dynamics of economic systems has become a prominent issue in economics and macroeconomics in recent years. Several nonlinear continuous models have been proposed to explicate the core features of economic data based on the dynamic behavior of the system [6,7,8,9,10,11,12,13,14,15]. The results of investigating the dynamics of an economic system with chaotic behavior and a suboptimal control proposal to suppress the chaotic behavior are presented [8]. A fractional order economic quantity model with time-varying holding cost is discussed in detail with the help of numerical computations [9]. Based on the definition of Atangana–Baleanu–Caputo fractional derivative, the integer-order financial chaotic system with nonconstant demand elasticity is extended to a fractional-order system, and its nonlinear dynamic properties are analyzed [10]. The chaotic complexity of a financial mathematical model in terms of a new generalized Caputo fractional derivative is analyzed [11]. The conditions for the structural stability of a fractional order IS-LM-AS dynamic model with adaptive expectations are given [12]. A dynamic fractional-order discrete gray model for forecasting China’s total renewable energy capacity is proposed [13]. A fault-tolerant prescribed performance control approach for fractional-order economic and supply chain systems is presented [14].
Predicting economic growth is an important subject in economics, and accurate prediction can facilitate integrated economic planning and the development of rational economic policies to promote healthy economic growth [16,17,18]. In previous studies, scholars mainly focused on the chaotic motion of the financial system, and domestic scholars have achieved certain research results on economic growth by introducing delays. The Solow growth model [19] provides a theoretical foundation for breakthroughs in economic growth and a research framework that can be applied in subsequent work. And the Solow growth model demonstrates how saving rates, capital stock, the labor force growth rate, technological progress, and capital depreciation influence a country’s total output. The economy tends towards a stable state and emphasizes that technological progress is the ultimate driving force for long-term growth based on Solow’s theory. Because most economic processes are not only influenced by current states but also greatly rely on past relevant factors and indicators, mathematical models with time delay are more suitable for describing economic phenomena. Recently, the global attractivity of the quasi-periodicity of a new class of delayed classical growth models are proved [20]. The fractional order models serve to forecast the economic growth of Group of Twenty countries [21]. Based on the Solow model, a fractional-order time-delayed economic growth model is established to effectively capture memory characteristics in the economic growth path and explore the underlying growth factors [22]. Many results have been achieved using the Solow economic model [22,23,24,25], and it still has significant theoretical and practical value, particularly when considering fractional calculus theory, which is expected to yield meaningful results.
The economic system is an organic whole composed of interconnected and interactive economic elements. When addressing the issues related to the economic system, it is necessary to consider not only its economic benefits but also the impact of such benefits on the ecological environment [26]. The environment system has a certain self-regulation and self-recovering capacity. However, excessive pollution beyond its self-regulating capacity can cause irreversible damage. Therefore, the stability of the economic system still depends on the capacity of the environment system, external material exchange, and energy flow [27,28]. Furthermore, environmental pollution will inevitably have an impact on the economy, such as water pollution affecting crops [29,30]. Serious economic losses caused by nuclear pollution in Fukushima, Japan have affected economic development [31,32]. To better characterize the dynamic laws of the operation of economic and environmental systems, economic and environment systems are integrated to form an environmental economic growth system. Economic growth has always been an issue of great interest in macroeconomic research. However, with rapid economic growth, environmental pollution has posed a serious threat to human social development. Using the interplay between economic growth and environmental quality to analyze that negative impact has a certain time lag; that is, it is not instantaneous [26]. Economic growth, environmental pollution, and studies on the interactive influence between economic growth and environmental pollution are considered in Wuhan [33]. The importance and significance of the fractional order derivatives in the nonlinear environmental and economic model are provided [34].
Therefore, in this paper, environmental factors are incorporated into the classic Solow economic growth model, with the consideration that the main indicators in the economic environment system have the characteristic of “memory”, and a novel fractional-order time-delayed economic growth model with environmental purification is established. According to the Solow model, a fractional time-delayed economic model with environmental purification is provided to characterize the relationship between economic growth and environmental factors. Based on the strict assumption that technological progress is completely exogenous, the traditional Solow economic growth model is inconsistent with practical experience. Furthermore, environmental pollution will inevitably have an impact on the economy. Accordingly, with the consideration of environmental purification, a fractional economic growth model related to both capital and population is established in which technological progress is endogenous. Furthermore, a detailed stability analysis of the proposed model based on environmental purification is performed. The asymptotic stability condition of the equilibrium point is obtained and a stable parameter interval is provided. The influence of parameter variations on the stability of the established model is investigated.
The paper is structured as follows. In Section 2, some preliminaries and model descriptions are presented. The asymptotic stability conditions and parameter stability interval of the fractional-order time-delayed growth model with environmental pollution are provided in Section 3. Finally, a numerical simulation and discussion are given in Section 4 and Section 5.

2. Model Description

2.1. Definitions and Lemmas

The definition of the Caputo derivative has more advantages in actual applications, because its initial value has a measurable physical meaning [1,2,3]. So the definition of the Caputo derivative is adopted in this paper. In addition, according to its good explanation of the memory characteristics of economic variables [35], it can well capture economic processes.
Definition 1
([2]). The Caputo fractional derivative is defined as
D ι q f ( ι ) = 1 Γ ( n q ) 0 ι f ( n ) ( τ ) ( ι τ ) q n + 1 d τ ,
where n 1 < q < n , n Z + , Γ ( n q ) = 0 x n q 1 e x d x is the Gamma function [36].
Definition 2
([2,37]). The Laplace transform of the Caputo derivative is
{ D ι q f ( ι ) ; ς } = s q F ( ς ) k = 0 n 1 ς q k 1 f ( k ) ( 0 ) , n 1 < q < n ,
where F ( ς ) is the Laplace transform of f ( ι ) , i.e., F ( ς ) = 0 f ( ι ) e ς ι d ι .
The following linear fractional-order time-delayed system is considered:
D ι q X ( ι ) = A X ( ι ) + K X ( ι τ ) ,
where A = ( α i j ) n × n , K = ( κ i j ) n × n , X ( ι ) = ( x 1 ( ι ) , x 2 ( ι ) , , x n ( ι ) ) T , X ( ι τ ) = ( x 1 ( ι τ 1 ) , x 2 ( ι τ 2 ) , , x n ( ι τ n ) ) T .
The Laplace transform is applied to both sides of Equation (1), which yields the characteristic matrix [38]
Δ ( ς ) = ς q κ 11 e ς τ 11 α 11 κ 12 e ς τ 12 α 12 κ 1 n e ς τ 1 n α 1 n κ 21 e ς τ 21 α 21 ς q κ 22 e ς τ 22 α 22 κ 2 n e ς τ 2 n α 2 n κ n 1 e ς τ n 1 α n 1 κ n 2 e ς τ n 2 α n 2 s q κ n n e ς τ n n α n n .
According to the distribution of the eigenvalues of det ( Δ ( ς ) ) = 0 , the stability of system (1) is totally determined. If τ i j = 0 , Equation (1) can be rewritten as
D ι q X ( ι ) = A X ( ι ) + K X ( ι ) = Λ X ( ι ) ,
where the coefficient matrix Λ = A + K . Then, a useful Lemma about the stability of system (1) is introduced as:
Lemma 1
([38]). When q ( 0 , 1 ) , if all the eigenvalues of Λ satisfy | arg ( λ ) | > π 2 and the characteristic equation det ( Δ ( ς ) ) = 0 has no pure imaginary roots for any τ i j > 0 , i , j = 1 , , n , then the zero solution of Equation (1) is asymptotically stable by Lyapunov.

2.2. Model Description

In this subsection, a fractional-order time-delayed economic growth model with environmental pollution based on the Solow model is set up. According to [22,39], the following fractional-order generalization of the Solow-type equation with time delay is given as
D ι q L ι = r ( 1 L ι / L m ) L ι , D ι q K ι = s A ι K ι α L ι 1 α ( δ + g + n ) K ( ι τ ) , D ι q A ι = p A ι + w ¯ L ι + h K ( ι τ ) ,
where q ( 0 , 1 ) is the fractional order; L ι is the working population; r ( 0 , 1 ) denotes the natural growth rate of the working population; K ι is the capital stock; τ represents time delay; s ( 0 , 1 ) is the constant saving rate; δ ( 0 , 1 ) is the capital depreciation rate; g ( 0 , 1 ) is the growth rate of technology; n ( 0 , 1 ) is the population growth rate; Y ι , K ι are defined as Cobb–Douglas production function (gross domestic product) and capital stock at time ι . A ι is the index of total factor productivity (TFP) and α ( 0 , 1 ) is the output elasticity of capital stock.
The economic–environmental system consists of two subsystems: the economic subsystem and environmental subsystem. These two subsystems operate through the exchange of materials and energy. During the process of capital accumulation, the discharge of pollutants generated by the economic system into the environmental system increases the stock of environmental pollutants. By contrast, the accumulated pollutants in the environmental system cause economic losses for the economic system, thereby inhibiting the accumulation of capital. The economic losses caused by environmental pollution stimulate investors in the economic system to invest in pollution control, thereby removing a portion of the pollution and reducing the stock of environmental pollutants in the environmental system. Consequently, this reduces the economic losses caused by the environmental system on the economic system. An environmental purification factor is introduced into the classical Solow model.
Investment in pollution control, denoted by E, is directly responsive to the economic losses caused by environmental pollution, represented by G. The larger the economic losses from pollution, the greater the corresponding investment in pollution control. The relationship between the two can be expressed as
E = ρ G ,
where ρ is the degree of pollution control investment. The amount of waste pollution, denoted by Z, is directly related to gross domestic product Y ι . According to environmental Kuznets curve (EKC) theory [33,40], the relationship between waste pollution and gross domestic product is described as follows:
Z ( ι ) = ϵ Y ι e λ Y ι ,
where ϵ > 0 is referred to as the environmental pollution index, which indicates the severity of environmental pollution and λ > 0 is a parameter. The pollution amount, denoted by R, is related to pollution control investment E and waste pollution amount Z.
To define a pollution control function R ( Z , E ) that describes the amount of pollution removed, this process has a saturation effect, and its functional form is assumed as follows:
R ( Z , E ) = Z E / ( E + ω Z ) ,
where ω is a control parameter.
The function ψ ( P ι ) represents the natural purification capacity for pollution and exhibits strong nonlinear characteristics. The environmental purification capacity ψ ( P ι ) is related to the stock of environmental pollutants P ι and can be expressed as follows:
ψ ( P ι ) = σ P ι d β d β + P ι β ,
where σ , β , and d are parameters, and β > 1 .
The economic losses caused by pollution, denoted by G, are related to the environmental purification capacity. It is assumed that G is proportional to the difference between the nonlinear environmental purification capacity and the linear purification capacity. The pollution economic loss function is expressed as
G ( Y ι , P ι ) = l Y ι P t β + 1 ( d β + P ι β ) ,
where l is a parameter that represents the sensitivity of economic losses to the level of pollution, and β and d are parameters introduced in the environmental purification function.
Remark 1.
Note that the ψ ( P ι ) and G ( Y ι , P ι ) are the case of the Hill function [41,42]. The exponent β determines the steepness of the switch occurring around P ι . Considering that the higher the β values are, the steeper the Hill functions, and stronger the hysteresis they create; then, β > 1 is chosen here.
The change in the capital stock as K ι includes an increase in investment s Y ( ι ) , economic losses caused by pollution G, capital depreciation. Based on system (4), the evolution equation for the capital stock K ι is given by
D ι q K t = s A ι K ι α L ι 1 α G ( Y ι , P ι ) ( δ + g + n ) K ( ι τ ) .
The change in environmental pollutant stock P ( ι ) includes the input of waste pollution Z, the removal of pollution due to pollution control investment R, and the natural purification capacity ψ ( P ι ) . The evolution equation for the environmental system is given by
D ι q P ι = Z ι R ( Z ι , P ι ) ψ ( P ι ) .
From the above analysis, the following fractional-order time-delayed economic growth model with environmental purification can be provided:
D ι q L ι = r ( 1 L ι / L m ) L ι , D ι q K ι = s A ι K ι α L ι 1 α G ( Y ι , P ι ) ( δ + g + n ) K ( ι τ ) , D ι q A ι = p A ι + w ¯ L ι + h K ( ι τ ) , D ι q P ι = Z ι R ( Z ι , P ι ) ψ ( P ι ) .
The functional relationship for system (5) is as follows:
Y ι = A ι K ι α L ι 1 α , E ι = ρ G ι , Z ι = ϵ Y ι e λ Y ι , R ( Z ι , E ι ) = Z ι E ι / ( E ι + ω Z ι ) , ψ ( P ι ) = σ P ι d β d β + P ι β , G ( Y ι , P ι ) = l Y ι P ι β + 1 ( d β + P ι β ) .
A summary of the definitions of the parameters and variables in system (5) and Equation (6) are shown in Table 1.

3. Main Results

In this section, the focus is the local stability of system (5). Lemma 1 is a local stability theorem for system (5). Local properties can be analyzed for stability using the eigenvalue distribution of Δ ( s ) . Thus, simplification (part of linearization) is performed in this step. Note that λ > 0 and Y ι > 0 , and then e λ Y ι 1 ; hence, Z ι ϵ Y ι is considered. The equation Z ι is approximated to first order using the Taylor expansion, which meets the approximation requirement. According to Equation (6), R ( Z ι , E ι ) = Z ι E ι / ( E ι + ω Z ι ) . If E ι ω Z ι , then R ( ι ) E ι / ω . Therefore, R ( Z ι , E ι ) 1 ω E ι , that is, it is approximated by a linear function. The ability of environmental purification is related to the level of environmental pollutants. When the level of pollutants is low, the self-purification ability of the environmental system is strong. However, when the level of pollutants reaches a certain upper limit, the self-purification ability of the environment gradually weakens. This characteristic can be described using Hill functions. Hence,
ψ ( P ι ) = σ P ι d β d β + P ι β a n d G ( Y ι , P ι ) = l Y ι P ι β + 1 ( d β + P ι β ) .
To obtain the stability conditions of system (5), P ι β d β can be chosen. Then,
ψ ( P ι ) σ d β P ι 1 β a n d G ( Y ι , P ι ) l Y ι P ι .
The simplification of system (5) can be obtained:
D ι q L ι = r ( 1 L ι / L m ) L ι , D ι q K ι = s A ι K ι α L ι 1 α G ( Y ι , P ι ) ( δ + g + n ) K ( ι τ ) , D ι q A ι = p A ι + w ¯ L ι + h K ( ι τ ) , D ι q P ι = Z ι R ( Z ι , P ι ) ψ ( P ι ) .
The functional relationship in Equation (7) is as follows:
Y ι = A ι K ι α L ι 1 α , E ι = ρ G ι , Z ι = ϵ Y ι ι e λ Y ι , R ( Z ι , E ι ) 1 ω E ι , ψ ( P ι ) σ d β P ι 1 β , G ( Y ι , P ι ) l Y ι P ι .
Remark 2.
Actually, the functional relationship in system (7) is
Y ι = A ι K ι α L ι 1 α , E ι = ρ G ( Y ι , P ι ) , Z ι = ϵ Y ι e λ Y ι , R ( Z ι , E ι ) = Z ι E ι / ( E ι + ω Z ι ) , ψ ( P ι ) = σ P ι d β d β + P ι β , G ( Y ι , P ι ) = l Y ι P ι β + 1 ( d β + P ι β ) .
In order to obtain the local stability of system (5), the Equation (6) is reduced as Equation (8). If Equation (6) is used, the equilibrium point equation is difficult to solve, and the Jacobian determinant at the equilibrium point is also quite complex. Hence, it is very hard to obtain stability conditions and the parameter stability interval of system (5). However, according to the facts mentioned, in [38,43], the consideration here is locality, which means that linear equations conform to the linear form of the model. Therefore, when fractional-order systems have the same linear form, their stability can be studied through their linear equations, regardless of the complexity of the original equations. Thus, the simplification Equation (8) is reasonable. Based on this fact, to obtain stability conditions, the parameter σ = 0 is considered.
Theorem 1.
When q ( 0 , 1 ) and σ = 0 , if λ = 0 , r > 0 , p + b 22 + c 22 < 0 , p ( b 22 + c 22 ) h b 23 > 0 and b 44 < 0 , b ¯ 2 4 χ 2 c ¯ + b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 or b ¯ 2 4 χ 2 c ¯ b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 , the positive equilibrium point of system (7) is locally asymptotically stable by Lyapunov, where χ is any real root of the equation 8 χ 3 4 c ¯ χ 2 + 2 ( b ¯ d ¯ 8 e ¯ ) χ + e ¯ ( 4 c ¯ b ¯ 2 ) d ¯ 2 = 0 , and the coefficients of the equation are given as
b ¯ = 2 ( p + b 22 ) ( cos q π cos q π 2 + sin ( ± q π ) sin ( ± q π 2 ) ) , c ¯ = ( p + b 22 ) 2 c 22 2 + 2 p b 22 cos q π , d ¯ = 2 [ p ( c 22 2 b 22 2 b 22 ) h c 22 b 23 ] cos q π 2 , e ¯ = p 2 b 22 2 ( p c 22 h b 23 ) 2 ,
where b 22 , b 23 , b 44 and c 22 satisfy
b 22 = α ( s l P ¯ ) A ¯ K ¯ α 1 L ¯ 1 α , b 23 = ( s l P ¯ ) K ¯ α L ¯ 1 α , b 44 = ρ l ω A ¯ K ¯ α L ¯ 1 α , c 22 = ( δ + g + n ) ,
and the ( L ¯ , K ¯ , A ¯ , P ¯ ) is a solution of the equation
r ( 1 L ¯ / L m ) L ¯ = 0 , s A ¯ K ¯ α L ¯ 1 α l P ¯ A ¯ K ¯ α L ¯ 1 α ( δ + g + n ) K ¯ = 0 , p A ¯ + w ¯ L ¯ + h K ¯ = 0 , ( ϵ ρ l ω P ¯ ) A ¯ K ¯ α L ¯ 1 α = 0 .
Proof. 
If λ = 0 , σ = 0 and r > 0 , from the Equation (7), we can obtain the equilibrium solution equation as follows:
r ( 1 L ¯ / L m ) L ¯ = 0 , s A ¯ K ¯ α L ¯ 1 α l P ¯ A ¯ K ¯ α L ¯ 1 α ( δ + g + n ) K ¯ = 0 , p A ¯ + w ¯ L ¯ + h K ¯ = 0 , ( ϵ ρ l ω P ¯ ) A ¯ K ¯ α L ¯ 1 α = 0 .
Assume a positive equilibrium point ( L ¯ , K ¯ , A ¯ , P ¯ ) is derived from Equation (11), then the linear centralized system (7) at the point ( L ¯ , K ¯ , A ¯ , P ¯ ) can be written as:
D ι q Φ ( ι ) = A Φ ( ι ) + B ¯ Φ ( ι ) + C ¯ Φ ( ι τ ) ,
where
Φ ( ι ) = ( L ι , K ι , A ι , P ι ) T , A = r 0 0 0 0 0 0 0 0 0 p 0 0 0 0 0 ,
and B ¯ and C ¯ are the Jacobian matrices at equilibrium point ( L ¯ , K ¯ , A ¯ , P ¯ ) as
B ¯ = b 11 0 0 0 b 21 b 22 b 23 b 24 0 0 0 0 0 0 0 b 44 , C ¯ = 0 0 0 0 0 c 22 0 0 0 h 0 0 0 0 0 0 = 0 0 0 0 0 ( δ + g + n ) 0 0 0 h 0 0 0 0 0 0 , where
b 11 = 2 r , b 21 = ( 1 α ) ( s l P ¯ ) K ¯ α L ¯ α , b 22 = α ( s l P ¯ ) A ¯ K ¯ α 1 L ¯ 1 α , b 23 = ( s l P ¯ ) K ¯ α L ¯ 1 α , b 24 = l A ¯ K ¯ α L ¯ 1 α , b 44 = ρ l ω A ¯ K ¯ α L ¯ 1 α .
The Laplace transform is applied to both sides of system (12), and the characteristic matrix is given as:
Δ ( s ) = s q + r 0 0 0 b 21 s q + ( δ + g + n ) e s τ b 22 b 23 b 23 0 h e s τ s q p 0 0 0 0 s q b 44 .
The d e t ( Δ ( s ) ) = 0 is calculated as
( s q + r ) ( s q b 44 ) ( s 2 q a s q + ( δ + g + n ) s q e s τ + b e s τ + p b 22 ) = 0 ,
where
a = p + b 22 , b = p ( δ + g + n ) h b 23 .
The Equation (14) has no pure imaginary roots for any τ > 0 which is verified for the next part. We testify the fact by contradiction. When r > 0 , there are obviously no pure imaginary roots in the equation s q + r = 0 . Hence, the following equation is considered:
s 2 q a s q + ( δ + g + n ) s q e s τ + b e s τ + p b 22 = 0 .
Assume that s = ϖ i = | ϖ | ( cos ( π 2 ) + i sin ( ± π 2 ) ) is one of the pure imaginary roots of Equation (15), where ϖ is a positive real number. Taking s into Equation (15) yields
| ϖ | 2 q ( cos q π + i sin ( ± q π ) ) c 22 | ϖ | q ( c o s q π 2 + i sin ( ± q π 2 ) ) ( cos ϖ τ i sin ϖ τ ) a | ϖ | q ( cos q π 2 + i sin ( ± q π 2 ) ) + b ( cos ϖ τ i sin ϖ τ ) + p b 22 = 0 .
Both the real part and imaginary part of Equation (16) are zero, so we have
| ϖ | 2 q cos q π a | ϖ | q cos q π 2 + p b 22 = c 22 | ϖ | q sin ( ± q π 2 ) sin ϖ τ + ( c 22 | ϖ | q cos q π 2 b ) cos ϖ τ
and
| ϖ | 2 q sin ( ± q π ) a | ϖ | q sin ( ± q π 2 ) = ( b c 22 | ϖ | q cos q π 2 ) sin ϖ τ + c 22 | ϖ | q sin ( ± q π 2 ) cos ϖ τ .
The following is considered:
β 1 = | ϖ | 2 q cos q π a | ϖ | q cos q π 2 + p b 22 , β 2 = | ϖ | 2 q sin ( ± q π ) a | ϖ | q sin ( ± q π 2 ) , ρ 1 = c 22 | ϖ | q cos q π 2 b , ρ 2 = c 22 | ϖ | q sin ( ± q π 2 ) .
Thus, from Equations (17) and (18), we can obtain
ρ 1 cos ϖ τ + ρ 2 sin ϖ τ = β 1 , ρ 1 sin ϖ τ + ρ 2 cos ϖ τ = β 2 .
Sum the squares of two equations in Equation (20) on both sides:
ρ 1 2 + ρ 2 2 = β 1 2 + β 2 2 .
Expressions β 1 , β 2 , ρ 1 and ρ 2 are substituted, which yields
| ϖ | 4 q + b ¯ | ϖ | 3 q + c ¯ | ϖ | 2 q + d ¯ | ϖ | q + e ¯ = 0 ,
where
b ¯ = 2 a ( cos q π cos q π 2 + sin ( ± q π ) sin ( ± q π 2 ) ) , c ¯ = a 2 c 22 2 + 2 p b 22 cos q π , d ¯ = 2 ( a p b 22 + c 22 b ) cos q π 2 , e ¯ = p 2 b 22 2 b 2 .
If 0 < q < 1 and χ satisfy the inequality b ¯ 2 4 χ 2 c ¯ + b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 or b ¯ 2 4 χ 2 c ¯ b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 , there exists no real roots in Equation (22), where χ is any real root of the equation 8 χ 3 4 c ¯ χ 2 + 2 ( b ¯ d ¯ 8 e ¯ ) χ + e ¯ ( 4 c ¯ b ¯ 2 ) d ¯ 2 = 0 ; so, for any τ > 0 , the characteristic equation d e t ( Δ ( s ) ) = 0 has no pure imaginary roots.
If τ = 0 , the matrix Λ of Equation (12) is
Λ = r 0 0 0 b 12 b 22 + c 22 b 23 b 24 0 h p 0 0 0 0 b 44
and its characteristic equation is
f ( λ ^ ) = ( λ ^ + r ) ( λ ^ b 44 ) [ λ ^ 2 ( b 22 + c 22 + p ) λ ^ + p ( b 22 + c 22 ) h b 23 ] .
Its eigenvalues λ ^ 1 , 2 , 3 , 4 can be obtained. When 0 < p < ( 1 α ) ( δ + g + n ) , then
λ ^ 1 = r < 0 , λ ^ 2 = b 44 < 0 , λ ^ 3 + λ ^ 4 = p + b 22 + c 22 < 0 , λ ^ 3 λ ^ 4 = p ( b 22 + c 22 ) h b 23 > 0 .
Based on (25), the four eigenvalues of the matrix Λ have negative real parts, and all the eigenvalues of Λ satisfy | arg ( λ ^ ) | > π 2 . Based on Lemma 1, the positive equilibrium point ( L ¯ , K ¯ , A ¯ , P ¯ ) of system (7) is asymptotically stable by Lyapunov. This completes the proof. □
Remark 3.
In this paper, a fractional-order time-delayed economic growth model with environmental purification is proposed. The established model considers not only the environment and economic production but also the labor force population and total factor productivity. The delayed fractional-order economic growth model without pollution is given in [22]. So, our results obtained in this paper are further extended results than [22] on the analysis between the environment and economic production.
According to Equation (11),
L ¯ = L m , P ¯ = ϵ ω ρ l , s A ¯ K ¯ α L ¯ 1 α ( δ + g + n ) K ¯ l P ¯ A ¯ K ¯ α L ¯ 1 α = 0 , p A ¯ + h K ¯ + w ¯ L ¯ = 0 .
Based on Equation (26), it is noteworthy that the equilibrium points K ¯ and A ¯ are not easy to obtain. Hence, certain conditions are provided to facilitate the determination of the equilibrium points. The corresponding stability conditions are provided as two corollaries.
Corollary 1.
If q ( 0 , 1 ) , when w ¯ = 0 , λ = 0 , σ = 0 , h < 0 , r > 0 , 0 < p < ( 1 α ) ( δ + g + n ) , l ω 1 ϵ ω s ρ < 0 , and b ¯ 2 4 χ 2 c ¯ + b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 or b ¯ 2 4 χ 2 c ¯ b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 , then the positive equilibrium point ( L ¯ , K ¯ , A ¯ , P ¯ ) of system (7) is locally asymptotically stable by Lyapunov, where χ is any real root of the equation 8 χ 3 4 c ¯ χ 2 + 2 ( b ¯ d ¯ 8 e ¯ ) χ + e ¯ ( 4 c ¯ b ¯ 2 ) d ¯ 2 = 0 , and the coefficients of the equation are given as
b ¯ = 2 ( p + α ( δ + g + n ) ) ( cos q π cos q π 2 + sin ( ± q π ) sin ( ± q π 2 ) ) , c ¯ = p 2 + 2 α p ( δ + g + n ) ( 1 + cos q π ) + ( α 2 1 ) ( δ + g + n ) 2 , d ¯ = 2 α p ( δ + g + n ) [ α ( δ + g + n ) 1 ] cos q π 2 , e ¯ = p 2 α 2 ( δ + g + n ) 2 .
Proof. 
If w ¯ = 0 , as a result of taking τ = 0 , then
L ¯ = L m , P ¯ = ϵ ω ρ l , s A ¯ K ¯ α L ¯ 1 α l P ¯ A ¯ K ¯ α L ¯ 1 α ( δ + g + n ) K ¯ = 0 , p A ¯ + h K ¯ = 0 .
Solving the Equation (28), we can obtain the positive solution ( L ¯ , K ¯ , A ¯ , P ¯ ) of system (7) as
L ¯ = L m , K ¯ = ρ p ( δ + g + n ) ( ϵ ω s ρ ) h L m ( 1 α ) 1 α , A ¯ = h p K ¯ , P ¯ = ϵ ω ρ l .
Based on Equation (29), the coefficients can be computed as
b 11 = 2 r , b 21 = ( 1 α ) ( δ + g + n ) ρ p ( δ + g + n ) ( ϵ ω s ρ ) h L m ( 1 α ) 1 α , b 22 = ( δ + g + n ) α , b 23 = p h ( δ + g + n ) , b 24 = ρ l ( δ + g + n ) ( ϵ ω s ρ ) ρ p ( δ + g + n ) ( ϵ ω s ρ ) h L m ( 1 α ) 1 α , b 44 = ρ 2 l ( δ + g + n ) ω ( ϵ ω s ρ ) ρ p ( δ + g + n ) ( ϵ ω s ρ ) h L m ( 1 α ) 1 α .
According to the stability conditions in Theorem 1, if b 44 < 0 , it can be concluded that l ω 1 ϵ ω s ρ < 0 . According to p + b 22 + c 22 < 0 and p ( b 22 + c 22 ) h b 23 > 0 , it can be concluded that 0 < p < ( 1 α ) ( δ + g + n ) .
Furthermore, according to Equation (9),
b ¯ = 2 ( p + α ( δ + g + n ) ) ( cos q π cos q π 2 + sin ( ± q π ) sin ( ± q π 2 ) ) , c ¯ = p 2 + 2 α p ( δ + g + n ) [ ( 1 + cos q π ) + ( α 2 1 ) ( δ + g + n ) ] , d ¯ = 2 α p [ α ( δ + g + n ) 2 ( δ + g + n ) ] cos q π 2 , e ¯ = p 2 α 2 ( δ + g + n ) 2 .
This completes the proof. □
Corollary 2.
If q ( 0 , 1 ) , when λ = 0 , p = 0 , σ = 0 , h < 0 , w ¯ > 0 , r > 0 , l ω ( s ρ ϵ ω ) w ¯ h < 0 , s ρ ϵ ω ρ > 0 , and b ¯ 2 4 χ 2 c ¯ + b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 or b ¯ 2 4 χ 2 c ¯ b ¯ 3 4 b ¯ c ¯ + 8 d ¯ 8 χ + b ¯ 2 4 c ¯ < 0 , then the positive equilibrium point ( L ¯ , K ¯ , A ¯ , P ¯ ) of system (7) is locally asymptotically stable by Lyapunov, where χ is any real root of the equation 8 χ 3 4 c ¯ χ 2 + 2 ( b ¯ d ¯ 8 e ¯ ) χ + e ¯ ( 4 c ¯ b ¯ 2 ) d ¯ 2 = 0 , and the coefficients of the equation are given by
b ¯ = 2 α ( δ + g + n ) ( cos q π cos q π 2 + sin ( ± q π ) sin ( ± q π 2 ) ) , c ¯ = ( α 2 1 ) ( δ + g + n ) 2 , d ¯ = 2 h s ρ ϵ ω ρ ( w ¯ h ) α ( δ + g + n ) L m cos q π 2 , e ¯ = [ h ( s ρ ϵ ω ) ρ ( w ¯ h ) α L m ] 2 .
Proof. 
If p = 0 , h < 0 and w ¯ > 0 , taking τ = 0 , then we can obtain
L ¯ = L m , P ¯ = ϵ ω ρ l , s A ¯ K ¯ α L ¯ 1 α l P ¯ A ¯ K ¯ α L ¯ 1 α ( δ + g + n ) K ¯ = 0 , w ¯ L ¯ + h K ¯ = 0 .
Solving the equations, and getting the unique positive solution ( L ¯ , K ¯ , A ¯ , P ¯ ) of system (7)
L ¯ = L m , K ¯ = w ¯ h L m , A ¯ = ρ ( δ + g + n ) ( s ρ ϵ ω ) ( w ¯ h ) 1 α , P ¯ = ϵ ω ρ l .
Based on Equation (32), the coefficients can be computed as
b 11 = 2 r , b 21 = ( 1 α ) w ¯ h ( δ + g + n ) , b 22 = α ( δ + g + n ) , b 23 = ( ϵ ω s ρ ) ρ ( w ¯ h ) α L m , b 24 = ρ l ( δ + g + n ) ( ϵ ω s ρ ) w ¯ h L m , b 44 = ρ 2 l ( δ + g + n ) ω ( s ρ ϵ ω ) w ¯ h L m .
According to the stability conditions in Theorem 1, when b 44 < 0 , it can be concluded that l ω ( s ρ ϵ ω ) ω ¯ h < 0 , if p = 0 , h < 0 , w ¯ > 0 , s ρ ϵ ω ρ > 0 ,
b 22 + c 22 = ( δ + g + n ) ( α 1 ) < 0 ,
h b 23 = h ( s ρ ϵ ω ) ρ ( w ¯ h ) α L m > 0
are obtained.
According to Equation (9),
b ¯ = 2 α ( δ + g + n ) ( cos q π cos q π 2 + sin ( ± q π ) sin ( ± q π 2 ) ) , c ¯ = ( α 2 1 ) ( δ + g + n ) 2 , d ¯ = 2 h s ρ ϵ ω ρ ( w ¯ h ) α ( δ + g + n ) L m cos q π 2 , e ¯ = [ h ( s ρ ϵ ω ) ρ ( w ¯ h ) α L m ] 2 .
This completes the proof. □
Remark 4.
Note that, if coefficient l = 0 , this means that economic losses and production capital caused by pollution are not considered in the economic growth model. In [22], the prediction of China’s economic growth based on the delayed fractional-order economic growth model without pollution is discussed and potential economic growth factors are explored. Based on this model, by considering the fractional orders as parameters and optimizing them, an appropriate fractional order based on the economic data of China from 1978 to 2020 is found and China’s GDP in the next 30 years is predicted using the fractional-order delayed economic growth model. The factors that drive short-term high-speed economic growth are also found. The results indicate that China has a declining population dividend and capital accumulation deceleration. Therefore, the TFP is increasing along with technological progress and innovation. Based on the fractional-order economic growth model in [22], the environmental purification factor is considered in this paper.

4. Numerical Analysis

In this section, the effectiveness of the theoretical results is demonstrated through three numerical examples, and the impact of system parameters is further investigated.
The ABM predictor–corrector algorithm [44] and the computed step h = 0.01 are used to solve the fractional-order time-delayed economic growth model with environmental purification. The specific values of some parameters of system (7) are shown in Table 2.
Remark 5.
The time delay τ is chosen as τ = 2 . According to [22], the fitting result with time delay τ = 2 matches the original K t best. Analysis shows that when τ = 0 , which represents the case without time delay, capital stock will be slightly overestimated because historical states are ignored. By contrast, when τ = 4 , because of the excessive emphasis on the role and impact of economic variables in the capital accumulation that results from considering previous historical values, the historical data of K t will be underestimated. Therefore, the time delay τ = 2 is chosen properly.
Inputting these coefficients into model (7) yields a more specific model as follows:
D ι q L ι = 0.05 ( 1 L ι / 8 ) L ι , D ι q K ι = 0.4 Y ι l Y ι P ι 0.2 K ( ι 2 ) , D ι q A ι = p A ι + w ¯ L ι + h K ( ι 2 ) , D ι q P ι = ϵ Y ι e λ Y ι 18 l ω Y ι P ι σ d β P ι 1 β , Y ι = A ι K ι α L ι 1 α .
Example 1.
p = 0.04 , α = 0.4 , w ¯ = 0.02 , h = 0.01 , ϵ = 1 , ω = 1 , λ = 0 , σ = 0 , l = 1 / 18 are used to verify Theorem 1.
The initial value is chosen as E 0 = ( 4.0152 , 0.1383 , 1.7788 , 0.1 ) . From Equation (11), we can calculate the positive equilibrium point of system (33) and obtain E ¯ = ( L ¯ , K ¯ , A ¯ , P ¯ ) = ( 8.000 , 19.978 , 0.994 , 1 ) .Those can be obtained through further calculation:
b 22 = α ( s l P ¯ ) A ¯ K ¯ α 1 L ¯ 1 α = 0.079 , b 23 = ( s l P ¯ ) K ¯ α L ¯ 1 α = 3.973 , b 44 = ρ l ω A ¯ K ¯ α L ¯ 1 α = 11.467 , c 22 = ( δ + g + n ) = 0.2 .
It is verified that the corresponding conditions in Theorem 1 are fulfilled. Based on Theorem 1, the positive equilibrium point E ¯ is asymptotically stable. The convergence behaviors of the solution curve of system (33) about fractional order q are shown in Figure 1. From Figure 1, the smaller the fractional order, the slower the convergence speed.
Example 2.
p = 0.02 , α = 0.5 , h = 0.01 , w ¯ = 0 , ϵ = 1 , ω = 1 , λ = 0 , σ = 0 , l = 1 / 18 are used to verify Corollary 1. The initial value E 0 = ( 4.0152 , 0.1383 , 1.7788 , 0.1 ) is chosen. From Equation (11), the positive equilibrium point of system (33) is obtained as E ¯ = ( L ¯ , K ¯ , A ¯ , P ¯ ) =   ( 8.000 , 0.1686 , 0.0843 , 1 ) .
The corresponding conditions are satisfied in Corollary 1. Then, the positive equilibrium point E ¯ is asymptotically stable. The convergence behaviors of the solution curve of system (33) about the fractional order q are shown in Figure 2. From Figure 2, the smaller the fractional order, the slower the convergence speed. Similar numerical results for Corollary 2 can be obtained. Hence, verification is omitted.
Remark 6.
According to Theorem 1, the condition s ρ ϵ ω > 0 is satisfied. This condition also can be given by Corollaries 1 and 2. Furthermore, s l P ¯ > 0 is obtained. According to [16,45], countries that have high savings/investment rates tend to be richer. If stable economic growth is to be maintained, saving rates need to be higher than the pollution rate, otherwise it will be difficult to achieve economic growth. Hence, this condition is perfectly logical and reasonable. Next, if this condition is not satisfied, is the system (33) still stable?
Example 3.
The initial value E 0 = ( 4.0152 , 0.1383 , 1.7788 , 0.1 ) is chosen. The parameters of the system (33) are chosen as
q = 0.9 , α = 0.4 , p = 0.04 , h = 0.01 , w ¯ = 0.02 , λ = 0 , σ = 0 , l = 1 / 18 .
ϵ = 5 , ω = 3 are chosen. s ρ ϵ ω = 7.8 < 0 are obtained. Then, system (33) is unstable, which is shown as Figure 3a. A locally enlarged view of t [ 0 , 20 ] is also shown as Figure 3a. ϵ = 7.2 , ω = 1 are chosen, s ρ ϵ ω = 0 is obtained. Then, system (33) is unstable, as shown in Figure 3b. A locally enlarged view of t [ 0 , 20 ] is also shown in Figure 3b. However, when ϵ = 6.2 , ω = 1 are chosen, s ρ ϵ ω = 1 > 0 is obtained. Then, system (33) is also unstable, as shown in Figure 3c. A locally enlarged view of t [ 0 , 20 ] is also shown in Figure 3c. If ϵ = 3 , ω = 2 are chosen, s ρ ϵ ω = 1.2 > 0 is obtained. Then, system (33) is stable. The equilibrium point E ¯ = ( L ¯ , K ¯ , A ¯ , P ¯ ) = ( 8 , 53.1003 , 9.4197 , 6.0000 ) is asymptotically stable by Lyapunov, as shown in Figure 3d. Hence, the condition s ρ ϵ ω > 0 is a sufficient condition.

5. Discussion

Note that the conditions λ = 0 , σ = 0 in the stability analysis of model (33), when λ = 0 , approximate waste pollution Z ( ι ) as a linear function. When σ = 0 , the function ψ ( P ι ) = 0 ; hence, the natural purification capacity is not considered in model (33). In this section, the effect of parameters λ , σ in the model (33) is mainly considered.
Parameters from Example 1 are considered:
p = 0.04 , α = 0.4 , w ¯ = 0.02 , ϵ = 1 , h = 0.01 , ω = 1 , l = 1 / 18 .
Then, system (33) becomes
D ι q L ι = 0.05 ( 1 L ι / 8 ) L ι , D ι q K ι = 0.4 Y ι 1 / 18 Y ι P ι 0.2 K ( ι 2 ) , D ι q A ι = 0.04 A ι + 0.02 L ι 0.01 K ( ι 2 ) , D ι q P ι = Y ι e λ Y ι Y ι P ι σ P ι 1 β , Y ι = A ι K ι α L ι 1 α .
When λ = 0 , σ = 0 , according to Example 1, system (34) is stable. If λ = 0 and σ [ 0 , 3.55 ] , system (34) has a non-negative equilibrium point. When σ = 1 and λ [ 0 , 0.7 ] , system (34) has a non-negative real equilibrium point. When λ = 0 , σ = 0.2 , 1 , 2 , 3 , 3.55 are chosen and system (34) is asymptotically stable. Additionally, the convergence behavior of the solution curve of system (34) is given in Figure 4. The influence of σ on convergent behavior about A ι is shown in Figure 4a. K ι is shown in Figure 4b and P ι is shown in Figure 4c. From Figure 4, A ι , K ι , and P ι decrease along with σ increasing. This conclusion is drawn from system (34) ( 0.4 Y ι 1 / 18 Y ι P 0.2 K = 0 , 0.04 A + 0.02 L 0.01 K = 0 , Y ι Y ι P ι σ P ι 1 β = 0 ). However, L ι remains unchanged along with σ increasing, as shown in Figure 4d.
When σ = 0 and λ [ 0 , 0.7 ] , system (34) is asymptotically stable by Lyapunov. This numerical result is similar to λ = 0 , σ [ 0 , 3.55 ] . When λ = 0 , 0.1 , 0.2 , 0.4 , 0.6 , 0.7 are chosen, Figure 5 shows the convergence behavior of the solution curve of the system (34). Taking λ = 0 and λ > 0 ( λ = 0.1 , 0.2 , 0.4 , 0.6 , 0.7 ) , the convergence behaviors of system (34) are similar. So, using λ = 0 is reasonable for getting the asymptotic stability conditions.

6. Conclusions

A fractional-order time-delayed economic growth model with environmental purification is proposed in this paper to analyze the interplay between economic growth and environmental pollution. Time delay is considered in capital stock to describe the lag effect and memory features in economic operations. The established model is proposed in the form of a fractional-order differential equation. The stability conditions of the established model are obtained, and the parameter stability interval are provided. The theoretical results are verified in the simulation. The convergence behaviors of the solution curve of the proposed model about the fractional order q are further discussed. The impacts of parameter variation on the stability of the proposed model are analyzed.
Some potential research directions of fractional-order economic growth model will be explored based on the proposed model. Note that the theoretical analysis of the proposed model is provided in this paper. In the future, the economic data and pollution data will be collected to analyze the relationship between economic growth and environmental pollution based on system (5) and Equation (6). Furthermore, a country’s economic development cannot be separated from its energy consumption, and energy consumption directly leads to a large number of carbon emissions. It is of great significance to study the relationship between carbon dioxide emissions and economic growth for the implementation of energy conservation, emission reduction, and the development of low-carbon economy in cities. Based on the proposed model in this paper, some carbon dioxide emissions models (such as Kaya model [46], carbon emission simultaneous model [47]) can be used to study the relationship between carbon dioxide emissions and economic growth.

Author Contributions

Conceptualization, formal analysis, writing—review and editing, Y.G.; funding acquisition, software, writing—original draft, H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the National Natural Science Foundation of China under Grant No. 61903386, the Disciplinary Foundation of Central University of Finance and Economics, the Research Fund of Beijing Information Science & Technology University under Grant No. 2021XJJ64 and the Qin Xin Talents Cultivation Program Beijing Information Science & Technology University under Grant No. QXTCP C202119.

Data Availability Statement

No new data were created or analyzed in this study.

Acknowledgments

The authors would like to thank the respected reviewers for their kind comments and the editorial office for their advice.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley-Interscience: New York, NY, USA, 1993. [Google Scholar]
  2. Podlubny, I. Fractional Differential Equations; Academic Press: London, UK, 1999. [Google Scholar]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: New York, NY, USA, 2006. [Google Scholar]
  4. Concepcion, M.A.; Chen, Y.; Vinagre, M.; Xue, D. Fractional-Order Systems and Controls; Springer: London, UK, 2010. [Google Scholar]
  5. Meerschaert, M.M.; Scalas, E. Coupled continuous time random walks in finance. Physical A 2006, 370, 114–118. [Google Scholar] [CrossRef]
  6. Chen, L.; Chai, Y.; Wu, R. Control and synchronization of fractional-order financial system based on linear control. Discrete Dyn. Nat. Soc. 2011, 2011, 958393. [Google Scholar] [CrossRef]
  7. Fallahgoul, H.; Focardi, S.; Fabozzi, F. Fractional Calculus and Fractional Processes with Applications to Financial Economics Theory and Application; Elsevier: London, UK, 2017. [Google Scholar]
  8. Tusset, A.M.; Fuziki, M.E.; Balthazar, J.M.; Andrade, D.I.; Lenzi, G.G. Dynamic analysis and control of a financial system with chaotic behavior including fractional order. Fractal Fract. 2023, 7, 535. [Google Scholar] [CrossRef]
  9. Pakhira, R.; Mond, B.; Ghosh, U.; Sarkar, S. An EOQ model with fractional order rate of change of inventory level and time-varying holding cost. Soft Comput. 2024, 28, 3859–3877. [Google Scholar] [CrossRef]
  10. Lin, X.; Wang, Y.; Wang, J.; Zeng, W. Dynamic analysis and adaptive modified projective synchronization for systems with Atangana-Baleanu-Caputo derivative: A financial model with nonconstant demand elasticity. Chaos Soliton Fract. 2022, 160, 112269. [Google Scholar] [CrossRef]
  11. Alzaid, S.S.; Kumar, A.; Kumar, S.; Alkahtani, B.S.T. Chaotic behavior of financial dynamical system with generalized fractional operator. Fractals 2023, 31, 2340056. [Google Scholar] [CrossRef]
  12. Navarro, C.E.B.; Tomé, R.M.B. Qualitative behavior in a fractional order IS-LM-AS macroeconomic model with stability analysis. Math. Comput. Simulat. 2024, 217, 425–443. [Google Scholar] [CrossRef]
  13. Xia, L.; Ren, Y.; Wang, Y. Forecasting China’s total renewable energy capacity using a novel dynamic fractional order discrete grey model. Expert. Syst. Appl. 2024, 239, 122019. [Google Scholar] [CrossRef]
  14. Alsaadi, F.E.; Bekiros, S.; Yao, Q.; Liu, J.; Jahanshahi, H. Achieving resilient chaos suppression and synchronization of fractional-order supply chains with fault-tolerant control. Chaos Soliton Fract. 2023, 174, 113878. [Google Scholar] [CrossRef]
  15. Rogosin, S.; Karpiyenya, M. Fractional models for analysis of economic risks. Fract. Calc. Appl. Anal. 2023, 26, 2602–2617. [Google Scholar] [CrossRef]
  16. Robert, J.B.; Xavier, S. Economic Growth; The MIT Press: Cambridge, UK, 1998. [Google Scholar]
  17. Robert, J.B. Macroeconomics: A Modern Approach; South-Western College Pub: La Jolla, CA, USA, 2007. [Google Scholar]
  18. Novales, A.; Fernández, E.; Ruiz, J. Economic Growth: Theory and Numerical Solution Methods; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  19. Solow, R.M. Neoclassical growth theory. Handb. Macroecon. 1999, 1, 637–667. [Google Scholar]
  20. Duan, L.; Huang, C. Existence and global attractivity of almost periodic solutions for a delayed differential neoclassical growth model. Math. Method Appl. Sci. 2017, 40, 814–822. [Google Scholar] [CrossRef]
  21. Tejado, I.; Pérez, E.; Valério, D. Fractional derivatives for economic growth modelling of the group of twenty: Application to prediction. Mathematics 2020, 8, 50. [Google Scholar] [CrossRef]
  22. Lin, Z.; Wang, H. Modeling and application of fractional-order economic growth model with time delay. Fractal Fract. 2021, 5, 74. [Google Scholar] [CrossRef]
  23. Mankiw, N.G.; Romer, D.; Weil, D.N. A contribution to the empirics of economic growth. Q. J. Econ. 1992, 107, 407–437. [Google Scholar] [CrossRef]
  24. Machado, J.T.; Mata, M.E. Pseudo phase plane and fractional calculus modeling of western global economic downturn. Commun. Nonlinear Sci. 2015, 22, 396–406. [Google Scholar] [CrossRef]
  25. Tarasova, V.V.; Tarasov, V.E. Elasticity for economic processes with memory: Fractional differential calculus approach. Fract. Differ. Calc. 2016, 6, 219–232. [Google Scholar] [CrossRef]
  26. Massimiliano, F.; Gori, L.; Guerrini, L.; Sodini, M. A continuous time economic growth model with time delays in environmental degradation. J. Inform. Optim. Sci. 2019, 40, 185–201. [Google Scholar]
  27. Chang, C. A multivariate causality test of carbon dioxide emissions, energy consumption and economic growth in China. Appl. Energ. 2010, 87, 3533–3537. [Google Scholar] [CrossRef]
  28. Aghavee, V.M.; Aloo, A.S.; Shirazi, J.K. Energy, environment, and economy interactions in Iran with cointegrated and ECM simultaneous model. Procedia Econ. Financ. 2016, 36, 414–424. [Google Scholar] [CrossRef]
  29. Dolan, F.; Lamontagne, J.; Link, R.; Hejazi, M.; Reed, P.; Edmonds, J. Evaluating the economic impact of water scarcity in a changing world. Nat. Commun. 2021, 12, 1915. [Google Scholar] [CrossRef] [PubMed]
  30. Boretti, A.; Rosa, L. Reassessing the projections of the world water development report. NPJ Clean Water 2019, 2, 15. [Google Scholar] [CrossRef]
  31. Hirose, K.; Povinec, P.P. Ten years of investigations of Fukushima radionuclides in the environment: A review on process studies in environmental compartments. J. Environ. Radioactiv. 2022, 251, 106929. [Google Scholar] [CrossRef] [PubMed]
  32. Lu, Y.; Yuan, J.; Du, D.; Sun, B.; Yi, X. Monitoring long-term ecological impacts from release of Fukushima radiation water into ocean. Geogr. Sustain. 2021, 2, 95–98. [Google Scholar] [CrossRef]
  33. Rao, C.; Yan, B. Study on the interactive influence between economic growth and environmental pollution. Environ. Sci. Poll. Res. 2020, 27, 39442–39465. [Google Scholar] [CrossRef]
  34. Weera, W.; Zamart, C.; Sabir, Z.; Raja, M.; Alwabli, A.S.; Mahmoud, S.R.; Wongaree, S.; Botmart, T. Fractional order environmental and economic model investigations using artificial neural network. CMC-Comput. Mater. Con. 2023, 74, 1735–1748. [Google Scholar] [CrossRef]
  35. Wang, J.; Li, H. Surpassing the fractional derivative: Concept of the memory-dependent derivative. Comput. Math. Appl. 2011, 62, 1562–1567. [Google Scholar] [CrossRef]
  36. Artin, E. The Gamma Function; Courier Dover Publications: New York, NY, USA, 2015. [Google Scholar]
  37. Krantz, S.G.; Kress, S.; Kress, R. Handbook of Complex Variables; Springer: Boston, MA, USA, 1999. [Google Scholar]
  38. Wang, H.; Yu, Y.; Wen, G. Stability analysis of fractional-order Hopfield neural networks with time delays. Neural Netw. 2014, 55, 98–109. [Google Scholar] [CrossRef]
  39. Szydłowski, M.; Krawiec, A.; Toboła, J. Nonlinear oscillations in business cycle model with time lags. Chaos Soliton Fract. 2001, 12, 505–517. [Google Scholar] [CrossRef]
  40. Pasche, M. Technical progress, structural change and the environmental kuznets curve. Ecol. Econ. 2002, 42, 381–389. [Google Scholar] [CrossRef]
  41. Scheffer, M.; Carpenter, S.; Foley, J.A.; Folke, C.; Walker, B. Catastrophic shifts in ecosystems. Nature 2001, 413, 591–596. [Google Scholar] [CrossRef] [PubMed]
  42. Carpenter, S.R.; Ludwig, D.; Brock, W.A. Management of eutrophication for lake subject to potentially irreversible change. Ecol. Appl. 1999, 9, 751–771. [Google Scholar] [CrossRef]
  43. Wang, H.; Yu, Y.; Wen, G.; Zhang, S. Stability analysis of fractional-order neural networks with time delay. Neural Process. Lett. 2015, 42, 479–500. [Google Scholar] [CrossRef]
  44. Bhalekar, S.A.C.H.I.N.; Daftardar-Gejji, V.A.R.S.H.A. A predictor-corrector scheme for solving nonlinear delay differential equations of fractional order. J. Fract. Calc. Appl. 2011, 1, 1–9. [Google Scholar]
  45. Charles, I.J.; Dietrich, V. Introduction to Economic Growth, 3rd ed.; W. W. Norton: New York, NY, USA, 2013. [Google Scholar]
  46. Kaya, Y. Impact of Carbon Dioxide Emission Control on GNP Growth: Interpretation of Proposed Scenarios; IPCC Response Strategies Working Group Memorandum: Geneva, Switzerland, 1989. [Google Scholar]
  47. Alhassan, H.; Kwakwa, P.A.; Donkoh, S.A. The interrelationships among financial development, economic growth and environmental sustainability: Evidence from Ghana. Environ. Sci. Pollut. 2022, 29, 37057–37070. [Google Scholar] [CrossRef]
Figure 1. Convergent behavior of system (33) about fractional order q. (a) Convergent behavior about fractional order q = 0.6 . (b) Convergent behavior about fractional order q = 0.7 . (c) Convergent behavior about fractional order q = 0.8 . (d) Convergent behavior about fractional order q = 0.9 .
Figure 1. Convergent behavior of system (33) about fractional order q. (a) Convergent behavior about fractional order q = 0.6 . (b) Convergent behavior about fractional order q = 0.7 . (c) Convergent behavior about fractional order q = 0.8 . (d) Convergent behavior about fractional order q = 0.9 .
Fractalfract 08 00361 g001
Figure 2. Convergent behavior of the system (33) about the fractional order q. (a) Convergent behavior about the fractional order q = 0.6 . (b) Convergent behavior about fractional order q = 0.7 . (c) Convergent behavior about the fractional order q = 0.8 . (d) Convergent behavior about the fractional order q = 0.9 .
Figure 2. Convergent behavior of the system (33) about the fractional order q. (a) Convergent behavior about the fractional order q = 0.6 . (b) Convergent behavior about fractional order q = 0.7 . (c) Convergent behavior about the fractional order q = 0.8 . (d) Convergent behavior about the fractional order q = 0.9 .
Fractalfract 08 00361 g002
Figure 3. Convergent behavior of system (33) about the parameters ϵ , ω . (a) Convergent behavior about ϵ = 5 , ω = 3 . (b) Convergent behavior about ϵ = 7.2 , ω = 1 . (c) Convergent behavior about ϵ = 6.2 , ω = 1 . (d) Convergent behavior about ϵ = 3 , ω = 2 .
Figure 3. Convergent behavior of system (33) about the parameters ϵ , ω . (a) Convergent behavior about ϵ = 5 , ω = 3 . (b) Convergent behavior about ϵ = 7.2 , ω = 1 . (c) Convergent behavior about ϵ = 6.2 , ω = 1 . (d) Convergent behavior about ϵ = 3 , ω = 2 .
Fractalfract 08 00361 g003
Figure 4. Convergent behavior of system (34) about parameters σ . (a) Influence of σ on the convergent behavior about A ι . (b) Influence of σ on convergent behavior of K ι . (c) Influence of σ on the convergent behavior of P ι . (d) Influence of σ on the convergent behavior of L ι .
Figure 4. Convergent behavior of system (34) about parameters σ . (a) Influence of σ on the convergent behavior about A ι . (b) Influence of σ on convergent behavior of K ι . (c) Influence of σ on the convergent behavior of P ι . (d) Influence of σ on the convergent behavior of L ι .
Fractalfract 08 00361 g004
Figure 5. Convergent behavior of system (34) about parameters λ . (a) Influence of σ on the convergent behavior of A ι . (b) Influence of σ on the convergent behavior of K ι . (c) Influence of σ on the convergent behavior of P ι . (d) Influence of σ on the convergent behavior of L ι .
Figure 5. Convergent behavior of system (34) about parameters λ . (a) Influence of σ on the convergent behavior of A ι . (b) Influence of σ on the convergent behavior of K ι . (c) Influence of σ on the convergent behavior of P ι . (d) Influence of σ on the convergent behavior of L ι .
Fractalfract 08 00361 g005
Table 1. Values of some parameters.
Table 1. Values of some parameters.
VariablesRepresentationsParametersRepresentations
Y ι Cobb–Douglas production function (gross domestic product)sConstant saving rate
A ι Total factor productivity (TFP) δ Capital depreciation rate
P ι Natural purification capacitygGrowth rate of technology
R ι Economic losses causedrNatural growth rate
ψ ι Environmental purification capacity ϵ Environmental pollution index
L ι Working populationqFractional order
K ι Capital stock L m Maximum number of labor force
E ι Investment in pollution controlnpopulation growth rate
Z ι Waste pollution τ time delay
G ι Pollution economic loss functionpTFP growth rate
ρ Degree of pollution control investment
α Output elasticity 0 < α < 1
w ¯ , l , h Scale factors
β Parameter and β > 1
d , ω , σ , λ Parameters
ι time
Table 2. Values of some parameters.
Table 2. Values of some parameters.
ParametersValuesParametersValues
Constant saving rate s s = 0.4 Maximum number of labor force L m L m = 8
Capital depreciation rate δ δ = 0.05 Population growth rate n n = 0.05
Growth rate of technology g g = 0.1 Time delay τ τ = 2
Natural growth rate r r = 0.05 Degree of pollution control investment ρ ρ = 18
Parameter β β = 1.6 Parameter d d = 1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gu, Y.; Wang, H. Stability Analysis of a Fractional-Order Time-Delayed Solow Growth Model with Environmental Pollution. Fractal Fract. 2024, 8, 361. https://doi.org/10.3390/fractalfract8060361

AMA Style

Gu Y, Wang H. Stability Analysis of a Fractional-Order Time-Delayed Solow Growth Model with Environmental Pollution. Fractal and Fractional. 2024; 8(6):361. https://doi.org/10.3390/fractalfract8060361

Chicago/Turabian Style

Gu, Yajuan, and Hu Wang. 2024. "Stability Analysis of a Fractional-Order Time-Delayed Solow Growth Model with Environmental Pollution" Fractal and Fractional 8, no. 6: 361. https://doi.org/10.3390/fractalfract8060361

APA Style

Gu, Y., & Wang, H. (2024). Stability Analysis of a Fractional-Order Time-Delayed Solow Growth Model with Environmental Pollution. Fractal and Fractional, 8(6), 361. https://doi.org/10.3390/fractalfract8060361

Article Metrics

Back to TopTop