Next Article in Journal
Generalized Bayes Estimation Based on a Joint Type-II Censored Sample from K-Exponential Populations
Previous Article in Journal
From Cell–Cell Interaction to Stochastic and Deterministic Descriptions of a Cancer–Immune System Competition Model
Previous Article in Special Issue
The Dynamical Behavior of a Three-Dimensional System of Exponential Difference Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Stability of Traveling Waves for the Lotka–Volterra Competition System with Three Species

1
School of Mathematics and Physics, University of South China, Hengyang 421001, China
2
Basic Department, Hunan Institute of Traffic Engineering, Hengyang 421001, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(9), 2189; https://doi.org/10.3390/math11092189
Submission received: 1 March 2023 / Revised: 21 April 2023 / Accepted: 28 April 2023 / Published: 6 May 2023
(This article belongs to the Special Issue Difference and Differential Equations and Applications)

Abstract

:
The stability of traveling waves for the Lotka–Volterra competition system with three species is investigated in this paper. Specifically, we first show the asymptotic behavior of traveling wave solutions and then establish the local stability and the global stability under the weighted functional space. For local stability, the spectrum approach is used, while for global stability, the comparison principle and squeezing theorem are combined.

1. Introduction

The aim of this paper is to study the stability of traveling waves for the Lotka–Volterra competition system with three species as follows:
u ˜ t = u ˜ x x + u ˜ 1 u ˜ b 12 v ˜ b 13 w ˜ , v ˜ t = d 1 v ˜ x x + α v ˜ 1 b 21 u ˜ v ˜ , w ˜ t = d 2 w ˜ x x + β w ˜ 1 b 31 u ˜ w ˜ .
To proceed, we first transform the variables so that u ˜ = u , v ˜ = 1 v , w ˜ = 1 w and the system (1) is converted into the following cooperative system:
u t = u x x + u ( 1 u b 12 + b 12 v b 13 + b 13 w ) , v t = d 1 v x x + α ( 1 v ) ( b 21 u v ) , w t = d 2 w x x + β ( 1 w ) ( b 31 u w ) ,
with the initial value u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x ) and w ( x , 0 ) = w 0 ( x ) for x R . In this system, u , v and w are the population densities of three species, respectively; d i ( i = 1 , 2 ) is the diffusion coefficient of species i; b 1 j and b j 1 ( j = 2 , 3 ) denote the competition coefficients between the other two species j and the first species; and α and β stand for the growth rates of the two species of v , w , respectively. All the coefficients are positive. Further, we can understand that there are three species u , v and w living together, and species u is a predator, while species v , w are both prey. However, v , w do not directly affect each other, and the predator u acts as a mediator for v and w.
The Lotka–Volterra model is well-known for better describing changes in biological populations, and many mathematicians are interested in its dynamics. In particular, many studies on the existence, stability, and invasion speed of traveling wave solutions have been generated on the two species competitive model, see [1,2,3,4,5,6,7,8,9]. For the three-species competition model, the studies on the dynamical behaviors are also receiving increased attention. The existence of traveling wave solutions for the three-species system has been extensively studied in [10,11,12,13,14]. In addition, many scholars [15,16,17,18] investigated the speed selection, and for more studies on other aspects of the three-species system, please see [19,20,21]. Among them, Pan et al. [15] converted the competitive system into a cooperative system and investigated the determinism of the invasion velocity by the upper and lower solution method. We shall directly employ some results in [15] for this study.
For a competitive system, understanding the conditions under which a species survives or dies is always an important and interesting topic in dynamics, and traveling wave solutions can be used to help us answer this question. By a simple calculation, we can find that (2) admits at least five equilibrium points in the range { ( u , v , w ) | 0 u 1 , 0 v 1 , 0 w 1 } , i.e., (0, 0, 0), (0, 1, 0), (0, 0, 1), (0, 1, 1) and (1, 1, 1). Then, this paper focuses on the traveling waves connecting the equilibrium points e 0 = ( 0 , 0 , 0 ) and e 1 = ( 1 , 1 , 1 ) in the form
( u , v , w ) ( x , t ) = ( U ¯ , V ¯ , W ¯ ) ( z ) , z = x c t ,
where c is called the wave speed and ( U ¯ , V ¯ , W ¯ ) is called the wave profile. For the convenience of discussion, we always assume that
1 2 < b 12 + b 13 < 1 , b 21 > 1 , b 31 > 1 .
A Similar assumption has been made in other papers studying the three-species model, such as [15,17,22], and the assumption is essential for stability properties in this paper. This condition means that v , w are weak competitors of u and it makes the point ( 0 , 0 , 0 ) unstable and the point ( 1 , 1 , 1 ) is stable. By substituting (3) into (2), we have
U ¯ z z + c U ¯ z + U ¯ 1 b 12 b 13 U ¯ + b 12 V ¯ + b 13 W = 0 , d 1 V ¯ z z + c V ¯ z + α ( 1 V ¯ ) b 21 U ¯ V ¯ = 0 , d 2 W ¯ z z + c W ¯ z + β ( 1 W ¯ ) b 31 U ¯ W ¯ = 0 , ( U ¯ , V ¯ , W ¯ ) ( ) = e 1 , ( U ¯ , V ¯ , W ¯ ) ( + ) = e 0 .
The existence of the traveling wave has been given in other related literature. Pan et al. [15] gave the existence of the traveling wave when c c * and the minimal wave speed is linearly determined for c * = c 0 = 2 1 b 12 b 13 . Apart from that, the asymptotic behavior of ( U ¯ , V ¯ , W ¯ ) near the equilibrium point ( 0 , 0 , 0 ) is also given in [15], see the following lemma.
Lemma 1. 
For any c > c 0 and constants C 1 > 0 , C 3 > 0 , C 4 > 0 , or C 1 = 0 with C 2 > 0 , C 3 > 0 , C 4 > 0 , when z , ( U ¯ , V ¯ , W ¯ ) has the following asymptotic behavior:
U ¯ ( z ) V ¯ ( z ) W ¯ ( z ) C 1 1 α b 21 Γ 1 μ 1 β b 31 Γ 2 μ 1 e μ 1 z + C 2 1 α b 21 Γ 1 μ 2 β b 31 Γ 2 μ 2 e μ 2 z + C 3 0 1 0 e μ 3 z + C 4 0 0 1 e μ 4 z ,
where Γ 1 ( μ ) = d 1 μ 2 c μ α , Γ 2 ( μ ) = d 2 μ 2 c μ β and
μ 1 ( c ) = 1 2 [ c c 2 4 ( 1 b 12 b 13 ) ] , μ 2 ( c ) = 1 2 [ c + c 2 4 ( 1 b 12 b 13 ) ] , μ 3 ( c ) = 1 2 d 1 ( c + c 2 + 4 d 1 α ) , μ 4 ( c ) = 1 2 d 2 ( c + c 2 + 4 d 2 β ) .
Throughout this article, for better determining the weight function later, we always assume that μ 1 is the minimum between μ i ( i = 1 , 2 , 3 , 4 ) . To make the assumption true, we summarize the required parameter conditions and we can find that restrictions are only proposed for c and d i ( i = 1 , 2 ) . It is not contrary with other assumptions in our paper.
In order to study the stability of the traveling wave, we need to determine the solution with ( u 0 ( x ) , v 0 ( x ) , w 0 ( x ) ) as the initial value whether converges to ( U ¯ , V ¯ , W ¯ ) . Hence, a change of variables ( u , v , w ) ( x , t ) = ( U , V , W ) ( z , t ) further transforms (2) into a partial differential model
U t = U z z + c U z + U 1 b 12 b 13 U + b 12 V + b 13 W , V t = d 1 V z z + c V z + α ( 1 V ) b 21 U V , W t = d 2 W z z + c W z + β ( 1 W ) b 31 U W , U ( z , 0 ) = u 0 ( z ) , V ( z , 0 ) = v 0 ( z ) , W ( z , 0 ) = w 0 ( z ) , z R .
We know that ( U ¯ , V ¯ , W ¯ ) is the steady-state to the above new system. We need to add the following extra assumption about the steady-state in order to obtain global stability:
U ¯ max b 12 V ¯ + b 13 W ¯ , 1 b 21 V ¯ , 1 b 31 W ¯ .
It is not difficult to find that we can demonstrate that the condition (9) is not empty by using the linear selection condition in Theorem 4.1 in [15]. By choosing V ¯ = b 21 U ¯ , W ¯ = b 31 U ¯ and combining the condition (4), we have b 12 V ¯ + b 13 W ¯ = ( b 12 b 21 + b 13 b 31 ) U ¯ U ¯ because of the linear selection condition 2 ( 1 b 12 b 13 ) + b 12 b 21 + b 13 b 31 0 .
The attention, which focused on the stability of traveling waves, has increased and various methods have been shed light on, where the weighted energy method and the spectral analysis were widely used. In terms of local stability, Hou and Li [23] demonstrated the local stability of traveling waves of nonlinear reaction-diffusion equations in different weighted Banach spaces by employing a new method to analyze the location of the spectra. To investigate the stability of the traveling wave solutions with non-critical wave speeds, Leung et al. [24] similarly analyzed the spectrum of the linearization operator in the exponentially weighted Banach space. In terms of global stability, Wu and Xing [25] proved that traveling front solutions with critical speeds are globally exponentially stable in some exponentially weighted spaces. By using a combination of the weighted energy method and the Green’s function technique, the global stability of monostable traveling waves for nonlocal time-delayed reaction-diffusion equations was given in [26]. For additional research on stability by using the weighted energy method, see also [27,28,29,30].
More specifically, in recent years, there have been numerous investigations on the stability of the Lotka–Volterra diffusion model. Chen et al. [31] applied the weighted energy method to study the nonlinear stability of a discrete three-species Lotka–Volterra competitive diffusion system with monostable traveling wavefronts. The global asymptotic stability of a diffuse multispecies Lotka–Volterra interaction model for the non-homogeneous coexistence equilibrium state was established by using the Lyapunov function method in [32]. Ma and Guo [33] combined the monotonic dynamical systems theory, the sub-super solutions method, master spectrum theory to study the global asymptotic stability of the coexisting steady state of a competitive Lotka–Volterra reaction-diffusion model with an advection term arising. Alhasanat and Ou [2] showed the global stability of the traveling waves of the Lotka–Volterra diffusion model by using the upper–lower solution method together with the squeezing technique. Further reading on the stability of the Lotka–Volterra diffusion model may be found at [28,34,35,36,37,38].
Research on the existence of traveling waves and the choice of linear and nonlinear minimal wave speeds for the three-species competition model has been successful. However, the stability of traveling waves has received less attention. In light of this, we investigate both the local and global stability of the steady-state ( U ¯ , V ¯ , W ¯ ) under the weighted functional space in this research.
Theorem 1. 
For any c > c * and the weight function w ( z )
w ( z ) = e a z z 0 , z > z 0 , 1 , z z 0
with some constants z 0 , a ( μ 1 , m i n { μ 2 , μ 3 , μ 4 } ) , the traveling wave solution ( U ¯ , V ¯ , W ¯ ) ( z ) is locally stable in the weighted functional space L w p , which is defined in Definition 1.
Theorem 2. 
Suppose c > c * , conditions (4)–(9) hold true and the initial data of the solution ( U , V , W ) ( z , t ) to (8) are
U 0 ( z ) = U ( z , 0 ) , V 0 ( z ) = V ( z , 0 ) , W 0 ( z ) = W ( z , 0 ) ,
which satisfy
( 0 , 0 , 0 ) ( U 0 , V 0 , W 0 ) ( z ) ( 1 , 1 , 1 ) , z R , lim inf z ( U 0 , V 0 , W 0 ) ( z ) > ( 0 , 0 , 0 )
and
U 0 ( z ) U ¯ ( z ) L w ( R ) , V 0 ( z ) V ¯ ( z ) L w ( R ) , W 0 ( z ) W ¯ ( z ) L w ( R ) ,
then the traveling wave solution exists globally with
( 0 , 0 , 0 ) ( U , V , W ) ( z , t ) ( 1 , 1 , 1 ) , ( z , t ) R × R + ,
and for positive constants k and η, there are
sup z R | U ( z , t ) U ¯ ( z ) | k e η t , t > 0 , sup z R | V ( z , t ) V ¯ ( z ) | k e η t , t > 0 , sup z R | W ( z , t ) W ¯ ( z ) | k e η t , t > 0 ,
i.e., any solution satisfying the conditions converges exponentially to the equilibrium solution ( U ¯ , V ¯ , W ¯ ) ( z ) .
Despite the fact that the local stability of the Lotka–Volterra competition system with three species has been demonstrated in [28], we refer to its methodology for the verification of global stability before introducing a new weighted functional space to prove our Theorem 1. The spectral problem is explored in the weighted functional space to determine the sign of the real part of the eigenvalues and further obtain the result of local stability. For global stability, to prove our Theorem 2, we construct the upper solution based on the assumptions (4)–(9), and then the comparison principle is utilized for global stability.
The rest of this paper is organized as follows. In Section 2, we linearize the model and perform a spectral analysis on it in the suitable weighted functional space, which led to the conclusion of local stability. Then, also under the weighted functional space, the global stability is proved by combining the upper-lower solution method and the squeezing theorem in Section 3. Conclusions are shown in Section 4.

2. The Local Stability

We first introduce a weighted functional space L w p different from [28] before studying the local stability for the subsequent proof of global stability.
Definition 1. 
L p ( R ) is the well-known Lebesgue space of integrable functions. Define a weighted functional space L w p as follows:
L w p = f ( z ) : w ( z ) f ( z ) L p ( R ) , p 1 , z R .
The norm is
f ( z ) L w p = w ( z ) | f ( z ) | p d z 1 / p ,
and the weight function is
w ( z ) = 1 w 1 ( z ) , 1 w 2 ( z ) , 1 w 3 ( z ) ,
where
w 1 ( z ) = e p ¯ z z 0 , z > z 0 , 1 , z z 0 , w 2 ( z ) = e q ¯ z z 0 , z > z 0 , 1 , z z 0 , w 3 ( z ) = e r ¯ z z 0 , z > z 0 , 1 , z z 0 ,
with some constants z 0 , p ¯ , q ¯ and r ¯ , where p ¯ , q ¯ , r ¯ are positive.
In this paper, we study the local stability in the presence of perturbations. By analyzing the behavior of the traveling waves under this small perturbation over a long period of time, the solution can be considered as locally stable if it converges to the steady-state solution.
Let
U ( z , t ) = U ¯ ( z ) + δ ϕ 1 ( z ) e λ t , V ( z , t ) = V ¯ ( z ) + δ ϕ 2 ( z ) e λ t , W ( z , t ) = W ¯ ( z ) + δ ϕ 3 ( z ) e λ t ,
where δ 1 , ϕ 1 ( z ) , ϕ 2 ( z ) , ϕ 3 ( z ) are real functions and λ is a parameter.
Let Φ = ( ϕ 1 , ϕ 2 , ϕ 3 ) T and in order to facilitate the exploration of the spectrum of the operator L on the space L w p , we write Φ in the following form:
Φ = ϕ 1 ϕ 2 ϕ 3 = w 1 ψ 1 w 2 ψ 2 w 3 ψ 3 ,
where ψ i ( i = 1 , 2 , 3 ) belong to L p .
By substituting (20) into (8) and linearizing it at ( U ¯ , V ¯ , W ¯ ) , we can obtain the following spectral problem:
λ Φ = L Φ : = D Φ + c Φ + J ( z ) Φ ,
where
D = 1 0 0 0 d 1 0 0 0 d 2 ,
J ( z ) = 1 b 12 b 13 2 U ¯ + b 12 V ¯ + b 13 W ¯ b 12 U ¯ b 13 U ¯ α b 21 ( 1 V ¯ ) α ( 1 b 21 U ¯ + 2 V ¯ ) 0 β b 31 ( 1 W ¯ ) 0 β ( 1 b 31 U ¯ + 2 W ¯ ) .
By examining the maximum real part sign of the spectrum λ of the operator L , we can now evaluate the local stability of the traveling wave solution.
Combine (19) and substitute (21) into (22) to obtain
λ Ψ = L w Ψ : = D Ψ + M ( z ) Ψ + N ( z ) Ψ ,
where Ψ = ( ψ 1 , ψ 2 , ψ 3 ) T ,
M ( z ) = c + 2 w 1 w 1 0 0 0 c + 2 d 1 w 2 w 2 0 0 0 c + 2 d 2 w 3 w 3 ,
and
N ( z ) = w 1 w 1 + c w 1 w 1 0 0 0 d 1 w 2 w 2 + c w 2 w 2 0 0 0 d 2 w 3 w 3 + c w 3 w 3 + 1 b 12 b 13 2 U ¯ + b 12 V ¯ + b 13 W ¯ b 12 U ¯ w 2 w 1 b 13 U ¯ w 3 w 1 α b 21 ( 1 V ¯ ) w 1 w 2 α 1 b 21 U ¯ + 2 V ¯ 0 β b 31 ( 1 W ¯ ) w 1 w 3 0 β 1 b 31 U ¯ + 2 W ¯ .
Then, we can use the following details from Theorem A.2 in [39] to determine the essential spectrum of the operator L w . After choosing the weight function to compel the essential spectrum to locate in the left-half complex plane, we may determine the sign of the maximum real part of the point spectrum in the weighted space. We choose
p ¯ ( μ 1 , μ 2 ) , q ¯ ( 0 , μ 3 ) , r ¯ ( 0 , μ 4 )
such that
p ¯ μ 1 < q ¯ p ¯ , p ¯ μ 1 < r ¯ p ¯ ,
where μ i ( i = 1 , 2 , 3 , 4 ) are defined in (7). M ( z ) and N ( z ) are bound by the preconditions mentioned above. Therefore, we define
lim z ± M ( z ) = M ± , lim z ± N ( z ) = N ± ,
and an algebraic curves S ± ,
S ± : = { λ det ( τ 2 D + i τ M ± + N ± λ I ) = 0 , < τ < } .
The union of areas within or on the curves S + and S contains the essential spectrum of the operator L w . If we prove that max ( R e ( λ ) ) < 0 for z ± , respectively, then S ± are on the left-half complex plane, which implies that the essential spectrum of L w lies on the left-half complex plane, for further details, see [28].
Because μ 1 is the smallest parameter, we choose p ¯ = q ¯ = r ¯ = a , where a is a constant and a ( μ 1 , m i n { μ 2 , μ 3 , μ 4 } ) . Then, the weight function is as follows:
w ( z ) = e a z z 0 , z > z 0 , 1 , z z 0 ,
where z 0 is defined in (19).
In order to obtain the local stability for (22), we next determine the sign of the major eigenvalue in the point spectrum.
Lemma 2. 
For Φ L w p , the real part of the eigenvalue λ of (22) is all negative.
Proof. 
Consider an associated linear partial differential system
f t = D f z z + c f z + J ( z ) f ,
where f ( z , t ) = ( f 1 ( z , t ) , f 2 ( z , t ) , f 3 ( z , t ) ) . By comparing with (22), we know that e λ t Φ is a solution of the above system with the same λ and Φ as in (22). According to the well-known Krein–Rutman theorem in [40], a compact linear operator which is strongly positive has a simple principal eigenvalue with a strongly positive eigenvector. For each given initial data set ϕ L p , let R t = f ( t , z , ϕ ) indicate the solution semiflow of (33) and we can find that R t satisfies the requirements of the theorem. So, we have
| e λ t | < e λ m a x t ,
where λ m a x is the simple principal eigenvalue. To proceed, we can prove λ < 0 by contradiction for two cases.
Case 1. λ = 0 .
For any c > c * , obviously we have U ¯ ( z ) C 1 e μ 1 z , C 1 > 0 , as z and V ¯ ( z ) , W ¯ ( z ) are the same asymptotic behaviors as U ¯ ( z ) . So, the operator L defined in (22) has an eigenvalue λ = 0 with the one-sign eigenvector ( U ¯ , V ¯ , W ¯ ) ( z ) , which is strongly positive. Because of (32) and a > μ 1 , we can check that ( U ¯ , V ¯ , W ¯ ) ( z ) is not inside the weighted functional space L w p .
Case 2. λ > 0 .
For z + and Φ L w p , there is obviously that Φ ¯ ( z ) = ( U ¯ , V ¯ , W ¯ ) ( z ) > Φ ( z ) except for sets of measure zero.
For z , assume that the eigenfunction of (22) possesses the asymptotic behavior as k e μ z for some positive constants k , μ . Thus, we can obtain the corresponding characteristic equation as follows:
μ 2 + c μ 1 λ b 12 b 13 0 d 1 μ 2 + c μ + α ( 1 b 21 ) λ 0 0 0 d 2 μ 2 + c μ + β ( 1 b 31 ) λ = 0 ,
which exist three positive roots,
μ ^ 1 = c + c 2 + 4 ( 1 + λ ) 2 , μ ^ 2 = c + c 2 + 4 d 1 ( λ + α ( b 21 1 ) ) 2 , μ ^ 3 = c + c 2 + 4 d 2 ( λ + β ( b 31 1 ) ) 2 .
These statements demonstrate that μ increases with λ , which also implies that Φ ¯ ( z ) > Φ ( z ) .
Hence, we can choose a sufficiently large k ¯ such that k ¯ Φ ¯ | Φ | . Both Φ ¯ and | Φ | e λ t are solutions of (33), so by using the comparison principle, we have k ¯ Φ ¯ | Φ | e λ t , which is not correct for a sufficiently large t. Thus, we complete the proof. □

3. The Global Stability

This section will analyze the global stability of the equilibrium solution ( U ¯ , V ¯ , W ¯ ) ( z ) in the weighted functional space L w p with p = , where the norm is defined as f L w = e s s s u p z R | w ( z ) f ( z ) | with a special weight function w ( z ) . Based on Theorem 1, we choose the weight function w ( z ) , which is defined in (32), and let
a = μ 1 + ϵ ,
where ϵ is a sufficiently small positive number.
First of all, assume the initial data of the solution of (8) as
U 0 ( z ) = U ( z , 0 ) , V 0 ( z ) = V ( z , 0 ) , W 0 ( z ) = W ( z , 0 ) ,
which satisfy
( 0 , 0 , 0 ) ( U 0 , V 0 , W 0 ) ( z ) ( 1 , 1 , 1 ) , z R , lim inf z ( U 0 , V 0 , W 0 ) ( z ) > ( 0 , 0 , 0 )
and
U 0 ( z ) U ¯ ( z ) L w ( R ) , V 0 ( z ) V ¯ ( z ) L w ( R ) , W 0 ( z ) W ¯ ( z ) L w ( R ) .
Based on the above conditions, for z R , we define
U 0 + ( z ) = max U 0 ( z ) , U ¯ ( z ) , U 0 ( z ) = min U 0 ( z ) , U ¯ ( z ) , V 0 + ( z ) = max V 0 ( z ) , V ¯ ( z ) , V 0 ( z ) = min V 0 ( z ) , V ¯ ( z ) , W 0 + ( z ) = max W 0 ( z ) , W ¯ ( z ) , W 0 ( z ) = min W 0 ( z ) , W ¯ ( z )
which can be viewed as the initial value of the solutions ( U + , V + , W + ) and ( U , V , W ) for (8). That is to say that ( U + , V + , W + ) and ( U , V , W ) satisfy
U t ± = U z z ± + c U z ± + U ± ( 1 b 12 b 13 U ± + b 12 V ± + b 13 W ± ) , V t ± = d 1 V z z ± + c V z ± + α ( 1 V ± ) ( b 21 U ± V ± ) , W t ± = d 2 W z z ± + c W z ± + β ( 1 W ± ) ( b 31 U ± W ± ) , ( U 0 ± , V 0 ± , W 0 ± ) ( z ) = ( U ± , V ± , W ± ) ( z , 0 ) .
By using the comparison principle, we have
( 0 , 0 , 0 ) ( U , V , W ) ( z , t ) ( U , V , W ) ( z , t ) ( U + , V + , W + ) ( z , t ) ( 1 , 1 , 1 ) , ( z , t ) R × R + , ( 0 , 0 , 0 ) ( U , V , W ) ( z , t ) ( U ¯ , V ¯ , W ¯ ) ( z ) ( U + , V + , W + ) ( z , t ) ( 1 , 1 , 1 ) , ( z , t ) R × R + .
Next, we need to demonstrate the convergence of ( U ± , V ± , W ± ) ( z , t ) to the wavefront ( U ¯ , V ¯ , W ¯ ) ( z ) in the subsequent lemmas, respectively.
Lemma 3. 
Under the conditions (38)–(42), ( U + , V + , W + ) ( z , t ) converges to ( U ¯ , V ¯ , W ¯ ) ( z ) .
Proof. 
We define
F ( z , t ) = U + ( z , t ) U ¯ ( z ) , G ( z , t ) = V + ( z , t ) V ¯ ( z ) , H ( z , t ) = W + ( z , t ) W ¯ ( z ) ,
with the initial value
F ( z , 0 ) = U 0 + ( z ) U ¯ ( z ) , G ( z , 0 ) = V 0 + ( z ) V ¯ ( z ) , H ( z , 0 ) = W 0 + ( z ) W ¯ ( z ) .
It is simple to see from inequality (43) that
( 0 , 0 , 0 ) ( F , G , H ) ( z , t ) ( 1 , 1 , 1 ) , z R ,
for t 0 . Afterwards, by combining (5) and (42) and performing some transformations, we obtain
F t = F z z + c F z + ( 1 b 12 b 13 ) F + ( F + U ¯ ) ( F + b 12 G + b 13 H ) + ( U ¯ + b 12 V ¯ + b 13 W ¯ ) F , G t = d 1 G z z + c G z + α ( b 21 F G ) + α ( G + V ¯ ) ( b 21 F + G ) + α ( b 21 U ¯ + V ¯ ) G , H t = d 2 H z z + c H z + β ( b 31 F H ) + β ( H + W ¯ ) ( b 31 F + H ) + β ( b 31 U ¯ + W ¯ ) H .
Let
F G H ( z , t ) = e a ( z z 0 ) F ¯ G ¯ H ¯ ( z , t ) , ( z , t ) ( R , R + ) ,
where F ¯ , G ¯ , H ¯ L ( R ) and z 0 is defined in (32). We will then demonstrate this lemma in two scenarios.
Case 1. Assume z z 0 , + for any fixed z 0 .
Substituting (48) into (47), we have
F ¯ G ¯ H ¯ t = D F ¯ G ¯ H ¯ z z + Q F ¯ G ¯ H ¯ z + A ( a ) F ¯ G ¯ H ¯ + ( U ¯ + b 12 V ¯ + b 13 W ¯ ) F ¯ α ( b 21 U ¯ + V ¯ ) G ¯ β ( b 31 U ¯ + W ¯ ) H ¯ + ( e a ( z z 0 ) F ¯ + U ¯ ) ( F ¯ + b 12 G ¯ + b 13 H ¯ ) α ( e a ( z z 0 ) G ¯ + V ¯ ) ( b 21 F ¯ + G ¯ ) β ( e a ( z z 0 ) H ¯ + W ¯ ) ( b 31 F ¯ + H ¯ ) : = L 1 ( F ¯ , G ¯ , H ¯ ) L 2 ( F ¯ , G ¯ , H ¯ ) L 3 ( F ¯ , G ¯ , H ¯ ) ,
where D is defined in (23),
Q = c 2 a 0 0 0 c 2 d 1 a 0 0 0 c 2 d 2 a
and
A ( a ) = Γ 3 ( a ) 0 0 α b 21 Γ 1 ( a ) 0 β b 31 0 Γ 2 ( a ) ,
where Γ 3 ( a ) = a 2 c a + 1 b 12 b 13 and Γ 1 , Γ 2 are given in Lemma 1. Assume that ( ξ 1 , ξ 2 , ξ 3 ) = ( ξ 1 ( a ) , ξ 2 ( a ) , ξ 3 ( a ) ) is the eigenvector of the matrix A ( a ) at eigenvalue a 2 c a + 1 b 12 b 13 and a direct calculation gives
ξ 1 = Γ 3 ( a ) Γ 1 ( a ) = ( 1 d 1 ) ( μ 1 2 + ϵ ) + 1 b 12 b 13 + α , ξ 2 = α b 21 , ξ 3 = β b 31 ( Γ 3 ( a ) Γ 1 ( a ) ) Γ 3 ( a ) Γ 2 ( a ) = β b 31 [ ( 1 d 1 ) ( μ 1 2 + ϵ ) + 1 b 12 b 13 + α ] ( 1 d 2 ) ( μ 1 2 + ϵ ) + 1 b 12 b 13 + β .
Then, we also define
F ¯ 1 ( z , t ) = k 1 ξ 1 e η 1 t , G ¯ 1 ( z , t ) = k 1 ξ 2 e η 1 t , H ¯ 1 ( z , t ) = k 1 ξ 3 e η 1 t , ( z , t ) ( R , R + ) ,
where k 1 , η 1 are positive. Since F ¯ ( z , 0 ) , G ¯ ( z , 0 ) , H ¯ ( z , 0 ) L w , thus we can choose
k 1 m a x z R { F ¯ ( z , 0 ) ξ 1 , G ¯ ( z , 0 ) ξ 2 , H ¯ ( z , 0 ) ξ 3 } .
For z + , by using (9), substituting (54) into the right side of (49) and performing the calculation, we find
L 1 ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) < 0 , L 2 ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) < 0 , L 3 ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) < 0 .
This means that we can find a suitable η 1 such that the inequality
F ¯ 1 G ¯ 1 H ¯ 1 t = η 1 k 1 ξ 1 ξ 2 ξ 2 e η 1 t L 1 ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) L 2 ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) L 3 ( F ¯ 1 , G ¯ 1 , H ¯ 1 )
holds.
Hence, ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) is equivalent to an upper solution. Then, by using the comparison principle on an unbounded domain, see [41], we have
( F , G , H ) ( z , t ) = ( F ¯ , G ¯ , H ¯ ) e a ( z z 0 ) ( F ¯ 1 , G ¯ 1 , H ¯ 1 ) e a ( z z 0 ) = k 1 ( ξ 1 , ξ 2 , ξ 3 ) e a ( z z 0 ) η 1 t , ( z , t ) z 0 , + × R + .
Now, we also need to verify the convergence of ( F , G , H ) to (0,0,0) at z , z 0 .
Case 2. Assume z , z 0 for any fixed z 0 .
System (47) can be represented in another form:
F G H t = D F G H z z + c F G H z + J ( z ) F G H + ( F + b 12 G + b 13 H ) F α ( b 21 F + G ) G β ( b 31 F + H ) H
where J ( z ) is defined in (24) and we write J ( z ) as J ( z ) = ( J i j ) 3 × 3 . Now, we present a new 3 × 3 matrix B ϵ 1 ,
B ϵ 1 = 1 + ϵ 1 b 12 + ϵ 1 b 13 + ϵ 1 ϵ 1 α ( 1 b 21 ) + ϵ 1 0 ϵ 1 0 β ( 1 b 31 ) + ϵ 1 = ( B i j ) 3 × 3 ,
for some given small ϵ 1 > 0 . When z , z 0 , due to the fact that ( U ¯ , V ¯ , W ¯ ) is nearing (1, 1, 1) for any z in this range, the inequality J i j < B i j ( i , j = 1 , 2 , 3 ) holds.
If we build an autonomous system related to B ϵ 1 with ( F ^ , G ^ , H ^ ) ( t ) as the solution:
F ^ G ^ H ^ t = B F ^ G ^ H ^ + ( F ^ + b 12 G ^ + b 13 H ^ ) F ^ α ( b 21 F ^ + G ^ ) G ^ β ( b 31 F ^ + H ^ ) H ^ ,
and the initial value satisfies
F ^ ( 0 ) F ¯ ( z , 0 ) , G ^ ( 0 ) G ¯ ( z , 0 ) , H ^ ( 0 ) H ¯ ( z , 0 ) , z R ,
then we can verify that ( F ^ , G ^ , H ^ ) ( t ) is an upper solution to the system (58).
We must now determine if ( F ^ , G ^ , H ^ ) ( t ) converges to ( 0 , 0 , 0 ) as t . We can use the Jacobi matrix J ( 0 , 0 , 0 ) to examine the behavior close to ( 0 , 0 , 0 ) , which is one of its fixed points. By using (60), the equation J ( 0 , 0 , 0 ) = B ϵ 1 has three eigenvalues denoted as λ ^ 3 < λ ^ 2 < λ ^ 1 < 0 . As a result, the point at ( 0 , 0 , 0 ) is stable, meaning that the flow in the F ^ G ^ H ^ -space converges to the origin for every ( F ^ , G ^ , H ^ ) ( 0 ) in the range [ 0 , 1 ] × [ 0 , δ 1 ] × [ 0 , δ 2 ] with 0 < δ i 1 ( i = 1 , 2 ) . The maximum possible value of δ i ( i = 1 , 2 ) depends on the position of the nonconstant fixed point to the system (3.24) near or inside the box [ 0 , 1 ] × [ 0 , 1 ] × [ 0 , 1 ] . If the point is far away from the box, then δ i ( i = 1 , 2 ) can be 1; If the point is near the boundary of the box, then the maximum possible value of δ 1 in ( b 21 1 ϵ 1 α , 1 ) and δ 2 in ( b 31 1 ϵ 1 β , 1 ) ; if the point is inside the box, then δ 1 is close to b 21 1 ϵ 1 α and δ 2 is close to b 31 1 ϵ 1 β . Thus, we find that
( F ^ , G ^ , H ^ ) = k ^ 1 ( C ^ 1 , C ^ 2 , C ^ 3 ) e λ ^ 1 t , t .
Here, k ^ 1 > 0 and ( C ^ 1 , C ^ 2 , C ^ 3 ) is the eigenvector of B ϵ 1 with the eigenvalue λ ^ 1 .
Finally, we have
( F , G , H ) ( z 0 , t ) k 1 ( ξ 1 , ξ 2 , ξ 3 ) e η 1 t k ^ 1 ( ξ 1 , ξ 2 , ξ 3 ) e λ ¯ 1 t
at z = z 0 by choosing a large enough k ^ 1 and λ ¯ 1 = min { η 1 , λ ^ 1 } . And by comparison on the domain , z 0 × 0 , , see [42], we find that
( F , G , H ) ( z , t ) k ^ 1 ( ξ 1 , ξ 2 , ξ 3 ) e λ ¯ 1 t , ( z , t ) , z 0 × R + .
Up to here, the proof is complete. □
Lemma 4. 
Under the above conditions (38)–(42), ( U , V , W ) ( z , t ) converges to ( U ¯ , V ¯ , W ¯ ) ( z ) .
Proof. 
We define
I ( z , t ) = U ¯ ( z ) U ( z , t ) , K ( z , t ) = V ¯ ( z ) V ( z , t ) , S ( z , t ) = W ¯ ( z ) W ( z , t ) ,
with the initial value
I ( z , 0 ) = U ¯ ( z ) U 0 ( z ) , K ( z , 0 ) = V ¯ ( z ) V 0 ( z ) , S ( z , 0 ) = V ¯ ( z ) W 0 ( z ) .
By inequalities (43), it is easy to see that
( 0 , 0 , 0 ) ( I , K , S ) ( z , t ) ( 1 , 1 , 1 ) , z R , t 0 .
Then, repeat the steps above, and I, K and S satisfy the system
I K S t = D I K S z z + c I K S z + J ( z ) I K S ( I + b 12 K + b 13 S ) I α ( b 21 I + K ) K β ( b 31 I + S ) S ,
where J ( z ) is defined in (24). Similarly, we analyze it in two cases.
Case 1. Let ( z , t ) z 0 , + × R + .
By using an approach similar to the proof of Lemma 3 with (9) and the facts I < U ¯ , K < V ¯ , S < W ¯ . There exist η 2 > 0 and
k 2 e a ( z z 0 ) m a x z R { I ( z , 0 ) ξ 1 , K ( z , 0 ) ξ 2 , S ( z , 0 ) ξ 3 }
such that
( I , K , S ) ( z , t ) k 2 ( ξ 1 , ξ 2 , ξ 3 ) e η 2 t , ( z , t ) z 0 , + × R + .
Case 2. Let ( z , t ) , z 0 × R + .
Now, we need to introduce w ( z ) defined in (32) with a = μ 1 + ϵ to study the stability under the weighted functional space L w p . Defined
I ^ K ^ S ^ t = B ϵ 1 I ^ K ^ S ^ 1 w ( z ) ( I ^ + b 12 K ^ + b 13 S ^ ) I ^ α ( b 21 I ^ + K ^ ) K ^ β ( b 31 I ^ + S ^ ) S ^ ,
and the initial date satisfies
I ^ ( 0 ) I ( z , 0 ) , K ^ ( 0 ) K ( z , 0 ) , S ^ ( 0 ) S ( z , 0 ) , z R .
We can check that ( I ^ , K ^ , S ^ ) ( t ) is an upper solution to the system (68). As in Lemma 3, ( I ^ , K ^ , S ^ ) also converges to ( 0 , 0 , 0 ) when all initial value on the space [ 0 , 1 ] × [ 0 , 1 ] × [ 0 , 1 ] except ( 1 , 1 , 1 ) by analyzing the phase plane. Finally, for some k ^ 2 , λ ¯ 2 > 0 , we have
( I , K , S ) ( z , t ) k ^ 2 ( ξ 1 , ξ 2 , ξ 3 ) e λ ¯ 2 t , ( z , t ) , z 0 × R + .
This completes the proof. □
In the end, we can prove Theorem 2 on the global stability.
Proof of Theorem 2. 
From (43), we have
| I ( z , t ) | | U ( z , t ) U ¯ ( z ) | | F ( z , t ) | , | K ( z , t ) | | V ( z , t ) V ¯ ( z ) | | G ( z , t ) | , | S ( z , t ) | | W ( z , t ) W ¯ ( z ) | | H ( z , t ) | ,
for ( z , t ) R × R + . Combining Lemmas 3 and 4 and the squeezing theorem, it is easy to find that, for all ( z , t ) R × R + ,
| U ( z , t ) U ¯ ( z ) | k e η t , t > 0 , | V ( z , t ) V ¯ ( z ) | k e η t , t > 0 , | W ( z , t ) W ¯ ( z ) | k e η t , t > 0 ,
where k , η > 0 . Hence, the proof is done. □

4. Conclusions

We examined if traveling waves in the Lotka–Volterra competition model with three species (2) display both local and global stability under the condition (4). Theorem 1 demonstrates, utilizing linearization and the crucial spectrum analysis, that the traveling wave solution is locally stable in a weighted functional space. Additionally, Theorem 2 demonstrates that all solutions converge to the wavefront solution using the upper-and-lower solution method and the squeezing theorem under the added constraint (9).

Author Contributions

S.H., C.P. and L.W. analyzed the method and revised the manuscript text together, S.H. and C.P. wrote and prepared the original draft, while L.W. supervised the writing. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Teaching Research and Reform Project of South China University of Technology grant (C9213136).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alhasanat, A.; Ou, C. Minimal-speed selection oftraveling waves to the Lotka-Volterra competition model. J. Differ. Equ. 2019, 266, 7357–7378. [Google Scholar] [CrossRef]
  2. Alhasanat, A.; Ou, C. Stability of Traveling Waves to the Lotka-Volterra Competition Model. Complexcity 2019, 2019, 6569520. [Google Scholar] [CrossRef]
  3. Alhasanat, A.; Ou, C. On the conjecture for the pushed wavefront to the diffusive Lotka-Volterra competition model. J. Math. Biol. 2020, 80, 1413–1422. [Google Scholar] [CrossRef]
  4. Bao, X.; Wang, Z.C. Existence and stability of time periodic traveling waves for a periodic bistable Lotka-Volterra competition system. J. Differ. Equ. 2013, 255, 2402–2435. [Google Scholar] [CrossRef]
  5. Li, X.; Jiang, D.; Mao, X. Population dynamical behavior of non-autonomous Lotka-Volterra competitive system with random perturbation. Discret. Contin. Dyn. Syst. Ser. A 2009, 24, 523–545. [Google Scholar] [CrossRef]
  6. Huang, W. Problem on minimum wave speed for Lotka-Volterra reaction-diffusion competition model. J. Dyn. Differ. Equ. 2010, 22, 285–297. [Google Scholar] [CrossRef]
  7. Wang, J.; Liu, Q.; Luo, Y. The numerical analysis of the long time asymptotic behavior for Lotka-Volterra competition model with diffusion. Numer. Funct. Anal. Optim. 2019, 40, 1–21. [Google Scholar] [CrossRef]
  8. Wang, H.; Pan, C.; Ou, C. Existence, uniqueness and stability of forced waves to the Lotka-Volterra competition system in a shifting environment. Stud. Appl. Math. 2022, 148, 186–218. [Google Scholar] [CrossRef]
  9. Chen, S.; Shi, J. Global dynamics of the diffusive Lotka-Volterra competition model with stage structure. Calc. Var. Partial Differ. Equ. 2020, 59, 2–19. [Google Scholar] [CrossRef]
  10. Wong, J. The analysis of a finite element method for the three-species Lotka-Volterra competition-diffusion with Dirichlet boundary conditions. J. Comput. Appl. Math. 2009, 223, 421–437. [Google Scholar] [CrossRef]
  11. Hung, L. Traveling wave solutions of competitive-cooperative Lotka-Volterra systems of three species. Nonlinear Anal.-Real World Appl. 2011, 12, 3691–3700. [Google Scholar] [CrossRef]
  12. Hou, X.; Li, Y. Traveling waves in a three species competition-cooperation system. Commun. Pure Appl. Anal. 2017, 16, 1103–1120. [Google Scholar] [CrossRef]
  13. Mimura, M.; Tohma, M. Dynamic coexistence in a three-species competition-diffusion system. Ecol. Complex. 2015, 21, 215–232. [Google Scholar] [CrossRef]
  14. Chen, C.; Hung, L.; Mimura, M. Exact travelling wave solutions of three-species competition-diffusion systems. Discret. Contin. Dyn. Syst. Ser. B 2012, 17, 2653–2669. [Google Scholar] [CrossRef]
  15. Pan, C.; Wang, H.; Ou, C. Invasive speed for a competition-diffusion system with three species. Discret. Contin. Dyn. Syst. Ser. B 2022, 27, 3515–3532. [Google Scholar] [CrossRef]
  16. Liu, Q.; Liu, S.; Lam, K. Stacked invasion waves in a competition-diffusion model with three species. J. Differ. Equ. 2021, 271, 665–718. [Google Scholar] [CrossRef]
  17. Guo, J.; Wang, Y.; Wu, C. The minimal speed of traveling wave solutions for a diffusive three species competition system. Taiwan. J. Math. 2015, 19, 1805–1829. [Google Scholar] [CrossRef]
  18. Wu, Q.; Pan, C.; Wang, H. Speed determinacy of the traveling waves for a three species time-periodic Lotka-Volterra competition system. Math. Methods Appl. Sci. 2022, 45, 6080–6095. [Google Scholar] [CrossRef]
  19. Chen, C.; Hung, L.; Tohma, M. Semi-exact equilibrium solutions for three-species competition-diffusion systems. Hiroshima Math. J. 2013, 43, 176–206. [Google Scholar] [CrossRef]
  20. Kan-on, Y.; Mimura, M. Singular perturbation approach to a 3-component reaction-diffusion system arising in population dynamics. SIAM J. Math. Anal. 1998, 29, 1519–1536. [Google Scholar] [CrossRef]
  21. Chen, G.; Wu, S. Invasion Entire Solutions for a Three Species Competition-diffusion System. Taiwan. J. Math. 2018, 22, 859–880. [Google Scholar] [CrossRef]
  22. Guo, J.; Wu, C. The existence of traveling wave solutions for a bistable three-component lattice dynamical system. J. Differ. Equ. 2016, 260, 1445–1455. [Google Scholar] [CrossRef]
  23. Hou, X.; Li, Y. Local stability of traveling-wave solutions of nonlinear reaction-diffusion equations. Discret. Contin. Dyn. Syst. 2006, 15, 681–701. [Google Scholar] [CrossRef]
  24. Leung, W.; Hou, X.; Feng, W. Traveling wave solutions for Lotka-Volterra system re-visited. Discret. Contin. Dyn. Syst. Ser. B 2011, 15, 171–196. [Google Scholar] [CrossRef]
  25. Wu, Y.; Xing, X. Stability of traveling waves with critical speeds for p-degree Fisher-type equations. Discret. Contin. Dyn. Syst. 2008, 20, 1123–1139. [Google Scholar] [CrossRef]
  26. Mei, M.; Ou, C.; Zhao, X. Global Stability of Monostable Traveling Waves For Nonlocal Time-Delayed Reaction-Diffusion Equations. SIAM J. Math. Anal. 2010, 42, 2762–2790. [Google Scholar] [CrossRef]
  27. Lv, G.; Wang, M. Nonlinear stability of travelling wave fronts for delayed reaction diffusion equations. Nonlinearity 2010, 23, 845–873. [Google Scholar] [CrossRef]
  28. Chang, C. The stability of traveling wave solutions for a diffusive competition system of three species. J. Math. Anal. Appl. 2018, 459, 564–576. [Google Scholar] [CrossRef]
  29. Li, Y.; Li, W.; Yang, Y. Stability of traveling waves of a diffusive susceptible-infective-removed (SIR) epidemic model. J. Math. Phys. 2016, 57, 041504. [Google Scholar] [CrossRef]
  30. Guo, Y. Exponential stability analysis of travelling waves solutions for nonlinear delayed cellular neural networks. Dynam. Syst. 2017, 32, 490–503. [Google Scholar] [CrossRef]
  31. Chen, G.; Wu, S.; Hsu, C. Stability of traveling wavefronts for a discrete diffusive competition system with three species. J. Math. Anal. Appl. 2019, 474, 909–930. [Google Scholar] [CrossRef]
  32. Vargas-De-León, C. Global stability of nonhomogeneous coexisting equilibrium state for the multispecies Lotka-Volterra mutualism models with diffusion. Math. Methods Appl. Sci. 2021, 45, 2123–2131. [Google Scholar] [CrossRef]
  33. Ma, L.; Guo, S. Positive solutions in the competitive Lotka-Volterra reaction-diffusion model with advection terms. Proc. Am. Math. Soc. 2021, 143, 3013–3019. [Google Scholar] [CrossRef]
  34. Yuan, L. Global dynamics of a Lotka-Volterra competition-diffusion-advection system in heterogeneous environments. J. Math. Pures Appl. 2019, 121, 47–82. [Google Scholar]
  35. Lei, Y.; Liu, Z.; Zhou, L. Large Time Behavior in a Fractional Chemotaxis-Navier-Stokes System with Competitive Kinetics. Acta Appl. Math. 2022, 179, 3. [Google Scholar] [CrossRef]
  36. Xu, M.; Guo, S. Dynamics of a delayed Lotka-Volterra model with two predators competing for one prey. Discret. Contin. Dyn. Syst.-Ser. B 2022, 27, 5573. [Google Scholar] [CrossRef]
  37. Wang, Y.; Li, H.; Li, X. Travelling wave fronts of Lotka-Volterra reaction-diffusion system in the weak competition case. Proc. R. Soc. Edinb. Sect. A Math. 2022, 152, 912–938. [Google Scholar] [CrossRef]
  38. Li, Z. Stability and Hopf bifurcation analysis in a Lotka-Volterra competition-diffusion-advection model with time delay effect. Nonlinearity 2021, 34, 3271–3313. [Google Scholar] [CrossRef]
  39. Henry, D. Geometric Teory of Semilinear Parabolic Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  40. Hess, P. Periodic-Parabolic Boundary Value Problems and Positivity; Longman Scientifc Technical: Harlow, UK, 1991. [Google Scholar]
  41. Aronson, D.; Weinberger, H. Multidimensional nonlinear diffusion arising in population genetics. Adv. Math. 1978, 30, 33–76. [Google Scholar] [CrossRef]
  42. Tieme, H. Asymptotic estimates of the solutions of nonlinear integral equations and asymptotic speeds for the spread of populations. J. Reine Angew. Math. 1979, 306, 94–121. [Google Scholar]
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

Hu, S.; Pan, C.; Wang, L. Global Stability of Traveling Waves for the Lotka–Volterra Competition System with Three Species. Mathematics 2023, 11, 2189. https://doi.org/10.3390/math11092189

AMA Style

Hu S, Pan C, Wang L. Global Stability of Traveling Waves for the Lotka–Volterra Competition System with Three Species. Mathematics. 2023; 11(9):2189. https://doi.org/10.3390/math11092189

Chicago/Turabian Style

Hu, Shulin, Chaohong Pan, and Lin Wang. 2023. "Global Stability of Traveling Waves for the Lotka–Volterra Competition System with Three Species" Mathematics 11, no. 9: 2189. https://doi.org/10.3390/math11092189

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop