The global asymptotic behavior of a nonautonomous competitor-competitor-mutualist model is investigated, where all the coefficients are time-dependent and asymptotically approach periodic functions, respectively. Under certain conditions, it is shown that the limit periodic system of this asymptotically periodic model admits two positive periodic solutions (u1T,u2T,u3T),(u1T,u2T,u3T) such that uiT≤uiT (i=1,2,3), and the sector {(u1,u2,u3):uiT≤ui≤uiT,i=1,2,3} is a global attractor of the asymptotically periodic model. In particular, we derive sufficient conditions that guarantee the existence of a positive periodic solution which is globally asymptotically stable.
1. Introduction
In this paper, we investigate the global asymptotic behavior of solutions for the following competitor-competitor-mutualist diffusion model:
(1)u1t-d1Δu1=g1u1(1-u1a1-a2u21+a3u3)inΩ×(0,∞),u2t-d2Δu2=g2u2(1-b1u1-u2b2)inΩ×(0,∞),u3t-d3Δu3=g3u3(1-u3c1+c2u1)inΩ×(0,∞),(2)∂ui∂n=0on∂Ω×(0,∞),i=1,2,3,(3)ui(x,0)=ui0(x)onΩ,i=1,2,3,
where u1(x,t), u2(x,t), and u3(x,t) are the densities of a mutualist-competitor, a competitor, and a mutualist population, respectively. Ω⊂RN is a bounded smooth domain, ∂/∂n is an outward normal derivative on ∂Ω.
In 1983, Rai et al. [1] firstly presented and studied a general competitor-competitor-mutualist ordinary differential equation (ODE) model. Zheng [2] studied the problem (1)–(3) in the case where all coefficients are positive constants. He proved the local stability of the unique positive constant steady-state solution under suitable condition on the reaction rates by the method of spectral analysis for linearized operator. Xu [3] investigated the global asymptotic stability of the unique positive constant steady-state solution under some assumptions by the iteration method. Pao [4] considered the model with time delays, and, under a very simple condition on the reaction rates, proved that the time-dependent solution with any nontrivial initial function converges to the positive steady-state solution by the method of upper and lower solutions. Chen and Peng [5] proved some existence results concerning nonconstant positive steady-states for the model with cross-diffusion and demonstrated that the cross-diffusion can create patterns when the corresponding model without cross-diffusion fails. Li et al. [6] proved that this model with cross-diffusion possesses at least one coexistence state if cross-diffusions and cross-reactions are weak by the Schauder fixed point theory and the method of upper and lower solutions and its associated monotone iterations. Fu et al. [7] investigated the global asymptotic behavior and the global existence of time-dependent solutions for the model with cross-diffusion when the space dimension is at most 5. Very recently, Tian and Ling [8] proved that, under some conditions, a corresponding predator-prey-mutualist model with cross-diffusion admits at least a nonhomogeneous stationary solution by the stability analysis for the positive uniform solution and the Leray-Schauder degree theory and carried out numerical simulations for a Turing pattern.
For the model (1) with T-periodic coefficients, Tineo [9] studied the asymptotic behavior of positive solutions by the method of upper and lower solutions. Du [10] investigated the existence of positive T-periodic solutions by using the degree and bifurcation theories. Pao [11] proved the existence of maximal and minimal T-periodic solutions by the method of upper and lower solutions. Wang et al. [12] considered the local asymptotic behavior of the time-dependent solutions and the existence of periodic solutions to the model in an unbounded domain. Zhou and Fu [13] investigated the global asymptotic behavior of the time-dependent solutions and the existence of periodic solutions for the model with discrete delays. Very recently, replacing the usual -Δu term by a degenerate elliptic operator as -Δum, Wang and Yin [14] proved the existence of maximal and minimal T-periodic solutions to the model with time delays by the Schauder fixed point theorem. It is important to note that the uniqueness of positive periodic solution is not considered in the previous references.
When a3=0, (1) reduces to the competition diffusion system
(4)ut-d1Δu=u(a-bu-cv),vt-d2Δv=v(d-eu-fv),
where a=g1, b=g1/a1, c=g1/a2, d=g2, e=g2b1, and f=g2/b2. The system (4) is a diffusion extension of the well-known Lotka-Volterra system
(5)dudt=u(a-bu-cv),dvdt=v(d-eu-fv).
In the case that a, b, c, d, e, and f are positive T-periodic functions, the existence and asymptotic stability of periodic solutions for (5) was studied by Gopalsamy [15], Alvarez and Lazer [16], and Ahmad [17] in the 1980’s. The global asymptotic behavior of (5) was studied by Ahmad and Lazer [18] and Tineo [19]. Denote fL=infx∈Xf(x) and fM=supx∈Xf(x) for any function f:X→R. If a, b, c, d, e, and f are positive asymptotically T-periodic functions on R, Peng and Chen [20] proved that if the conditions
(6)aL-ε0dM+ε0>cM+ε0fL-ε0,dL-ε0aM+ε0>eM+ε0bL-ε0
are satisfied for a certain sufficiently small ε0>0, then any positive solutions of (5) asymptotically approach the unique positive periodic solution for the limit periodic system of (5).
It is well known that periodic reaction diffusion equations are of particular interests since they can take into account seasonal fluctuations occurring in the phenomena appearing in the models, and they have been extensively studied by many researchers (see, e.g., [9–14, 19, 21]). However, so far, the research work on asymptotically periodic systems is much fewer than on the periodic ones. In fact, asymptotically periodic systems describe our world more realistically and more accurately than periodic ones to some extent. Therefore, for asymptotically periodic systems, studying the dynamics behavior is important and necessary (see, e.g., [22–27]).
In this paper, we study the global asymptotic behavior of positive solutions for the asymptotically periodic system (1). Under some conditions, it is shown that any positive solutions of the models asymptotically approach the unique strictly positive periodic solutions of the corresponding periodic system. This means that the results in Tineo [9] and the results for ODE model in Peng and Chen [20] can be extended to the asymptotically periodic reaction diffusion system and the 3-species diffusion system, respectively. Furthermore, using the method of the present paper, we note that the corresponding conclusions hold for the time-dependent n-species Lotka-Volterra systems. More specifically, we provide a way of how to use the method of upper and lower solutions to study asymptotic behavior of solutions for asymptotically periodic reaction diffusion systems. As one can see, the optimal bounds and uniqueness of positive periodic solutions will play an important role in the study of the global asymptotic behavior of periodic solutions.
2. Permanence and Extinction
For the sake of convenience, we introduce the two signs ~ and ≺ for functions u,v:Ω¯×R→[0,∞). u is said to approach v asymptotically in notation, u~v, if limt→∞|u(x,t)-v(x,t)|=0 uniformly for x in Ω¯. Furthermore, if (φ1,φ2,…,φn) and (ψ1,ψ2,…,ψn) are vector functions, then (φ1,φ2,…,φn)~(ψ1,ψ2,…,ψn) if and only if φi~ψi(i=1,2,…,n). We say that u(x,t) is asymptotically smaller than v(x,t) and write u(x,t)≺v(x,t) if lim¯t→∞(u(x,t)-v(x,t))≤0 uniformly for x∈Ω¯. It is clear that u(x,t)≺v(x,t) if and only, if for any ε>0, there exists a corresponding t1>0 such that u(x,t)<v(x,t)+ε on Ω×[t1,∞).
Assume the following.
di,Ai,Bi,Ci, and Gi are positive smooth and T-periodic functions on Ω¯×R.
ai,bi,ci, and gi are positive smooth functions on Ω¯×R, and
(7)(ai,bi,ci,gi)~(Ai,Bi,Ci,Gi).
By (H2), the limit periodic system of (1), (2) is given as follows:
(8)u1t-d1Δu1=G1u1(1-u1A1-A2u21+A3u3)inΩ×(0,∞),u2t-d2Δu2=G2u2(1-B1u1-u2B2)inΩ×(0,∞),u3t-d3Δu3=G3u3(1-u3C1+C2u1)inΩ×(0,∞),∂ui∂n=0on∂Ω×(0,∞)(i=1,2,3).
As a complement, we state the following main result which comes from [9, Theorem 0.3].
Theorem 1.
Assume that (H1) holds, and
(9)A2MB2M<1+A3LC1L,(10)A1MB1M<1.
Then (8) has the periodic solutions (u1T,u2T,u3T) and (u1T,u2T,u3T) such that uiT≥ui≥uiT>0(i=1,2,3) for any positive T-periodic solution (u1,u2,u3) of (8). Moreover, given that ε>0 and a solution (u1,u2,u3) of (8) with ui(x,0)≥(≢)0, there exists t1>0 such that uiT(x,t)-ε<ui(x,t)<uiT(x,t)+ε on Ω×(t1,∞).
In order to get the conditions for the permanence of (1)–(3), we need to make the following optimal bounds.
Lemma 2.
If (9) and (10) hold and (u1,u2,u3) is a positive smooth T-periodic solution of (8), then
(11)εui≤ui≤δui(i=1,2,3),
where δu1 is the unique positive root of p1x2+q1x+r1=0 and
(12)p1=A3MC2M,r1=A1M(A2LB2L-1-A3MC1M),q1=1+A3MC1M-A1MA3MC2M-A1MA2LB1MB2L.εu1 is the unique positive root of p2x2+q2x+r2=0, and
(13)p2=A3LC2L,r2=A1L(A2MB2M-1-A3LC1L),q2=1+A3LC1L-A1LA3LC2L-A1LA2MB1LB2M,δu2=B2M-B1LB2Mεu1,εu2=B2L-B1MB2Lδu1,δu3=C1M+C2Mδu1,εu3=C1L+C2Lεu1.
Proof.
By the maximum principle (see Lemma 1.2 of [18]), we have
(14)1-u1MA1M-A2Lu2L1+A3Mu3M≥0,1-u1LA1L-A2Mu2M1+A3Lu3L≤0,1-B1Lu1L-u2MB2M≥0,1-B1Mu1M-u2LB2L≤0,1-u3MC1M+C2Mu1M≥0,1-u3LC1L+C2Lu1L≤0.
Hence, p1u1M2+q1u1M+r1≤0. Since p1>0 and r1<0 (by (9)), we can see immediately that u1M≤δu1. Similarly, if εu1 is the unique positive root of p2x2+q2x+r2=0, then u1L≥εu1. So,
(15)εu2≤u2L≤u2M≤δu2,εu3≤u3L≤u3M≤δu3.
Evidently, εu3>0. By (9) and (10), we have
(16)p1B1M2+q1B1M+r1>0,from which it follows that εu2>0. This completes the proof.
Corollary 3.
Assume that (9) and (10) hold. If Ai, Bi, and Ci are positive constants, then (9) has the unique positive periodic solution (r,B2(1-B1r),C1+C2r), where r is the unique positive root of p1x2+q1x+r1=0.
The main results in this section are the following theorems.
Theorem 4 (permanence).
Assume that (H1), (H2), (9), and (10) hold. Then (8) has the positive T-periodic solutions (u1T,u2T,u3T) and (u1T,u2T,u3T) such that uiT≤uiT(i=1,2,3). Moreover, if (u1,u2,u3) is the solution of (1)–(3) with smooth initial values ui0(x)≥(≢)0, then
(17)uiT≺ui≺uiT(i=1,2,3).
Remark 5.
Under the assumptions of Theorem 4, the system (1), (2) is permanent, the sector 〈uT,uT〉={u∈C(Ω¯×R):uT≤u≤uT} is a global periodic attractor of (1), (2), and its trivial and semitrivial periodic solutions are unstable. Furthermore, if Ai, Bi, and Ci are positive constants, then (r,B2(1-B1r),C1+C2r) is the unique globally asymptotically stable solution of (8).
Theorem 6.
Assume that (H1) and (H2) hold. Then one has the following conclusions.
(Extinction of u2) Assume that (9) holds and that A1LB1L≥1. Then (8) has a T-periodic solution (U1,U2,U3) such that U1>0, U2=0, U3>0, and
(18)limt→∞|ui(x,t)-Ui(x,t)|=0(i=1,2,3)
uniformly on Ω¯, for any positive solution (u1,u2,u3) of (1)–(3).
(Extinction of u1) Assume that (10) holds and that
(19)A2LB2L≥1+A3MC1M+A3MA1MC2M.
Then (8) has a T-periodic solution (U1,U2,U3) with U1=0, U2>0, and U3>0 satisfying (18), where (u1,u2,u3) is any positive solution of (1)–(3).
Proof of Theorem 4.
By (9) and (10), there exists a sufficiently small ε0>0 such that, for δ∈(0,ε0),
(20)[(G1+δ)(A2+δ)G1-δ]M[(G2+δ)(B2+δ)G2-δ]M<1+(A3L-δ)[(G3-δ)(C1-δ)G3+δ]L,(21)[(G1+δ)(A1+δ)G1-δ]M[(G2+δ)(B1+δ)G2-δ]M<1.
Consider two auxiliary systems as follows:
(22)u1t-d1Δu1=u1[(G1+δ)-(G1-δ)u1A1+δ-(G1-δ)(A2-δ)u21+(A3+δ)u3]inΩ×(0,∞),u2t-d2Δu2=u2[(G2+δ)u2B2-δ(G2-δ)-(G2+δ)(B1+δ)u1-(G2+δ)u2B2-δ]inΩ×(0,∞),u3t-d3Δu3=u3[(G3-δ)u3(C1+δ)+(C2+δ)u1(G3+δ)-(G3-δ)u3(C1+δ)+(C2+δ)u1]inΩ×(0,∞),∂ui∂n=0on∂Ω×(0,∞)(i=1,2,3),(23)u1t-d1Δu1=u1[(G1-δ)-(G1+δ)u1A1-δ-(G1+δ)(A2+δ)u21+(A3-δ)u3]inΩ×(0,∞),u2t-d2Δu2=u2[G2-δu2B2+δ(G2+δ)-(G2-δ)(B1-δ)u1-(G2-δ)u2B2+δ]inΩ×(0,∞),u3t-d3Δu3=u3[(G3-δ)-(G3+δ)u3(C1-δ)+(C2-δ)u1]inΩ×(0,∞),∂ui∂n=0on∂Ω×(0,∞)(i=1,2,3).
By (20), (21), and Theorem 1, (22) has the positive T-periodic solutions (U1δ,u2δ,U3δ) and (U1δ,u2δ,U3δ) such that Uiδ≤ui≤Uiδ(i=1,3) and u2δ≤u2≤u2δ, for any positive T-periodic solution (u1,u2,u3) of (22). Moreover, if (u1,u2,u3) is a solution of (22) with nontrivial nonnegative initial values, then, for any ε>0, there exists tε>0 such that
(24)Uiδ(x,t)-ε<ui(x,t)<Uiδ(x,t)+ε(i=1,3),u2δ(x,t)-ε<u2(x,t)<u2δ(x,t)+ε,
for all x∈Ω¯ and t>tε. Similarly, (23) has the positive T-periodic solutions (u1δ,U2δ,u3δ) and (u1δ,U2δ,u3δ) such that uiδ≤ui≤uiδ(i=1,3) and U2δ≤u2≤U2δ, for any positive T-periodic solution (u1,u2,u3) of (23). Furthermore, if (u1,u2,u3) is a solution of (23) with nontrivial nonnegative initial values, then, for the previous ε>0, there exists tε′>0 such that, for all x∈Ω¯ and t>tε′,
(25)uiδ(x,t)-ε<ui(x,t)<uiδ(x,t)+ε(i=1,3),U2δ(x,t)-ε<u2(x,t)<U2δ(x,t)+ε.
Now we prove
(26)uiδ≤uiT≤Uiδ,uiδ≤uiT≤Uiδ,
where (u1T,u2T,u3T) and (u1T,u2T,u3T) are positive T-periodic solutions of (8) (see Theorem 1). Let (U1,u2,U3), (p1,p2,p3), and (u1,U2,u3) be the solutions of (22), (8), and (23), respectively, which all satisfy the same initial conditions. It is easily testified that (U1,U2,U3), (u1,u2,u3) are the upper and lower solutions of (8) and (3), respectively. So from [28, Corollary 5.2.10], we see that
(27)ui≤pi≤Ui(i=1,2,3).
For sufficiently small m>0 and δ>0, define
(28)r1=[(G1+δ)(A1+δ)G1-δ]M,r2=[(G2+δ)(B2+δ)G2-δ]M,r3=[(G3+δ)(C1+δ)G3-δ]M+[(G3+δ)(C2+δ)G3-δ]Mr1,s0=[(G3-δ)(C1-δ)G3+δ]L+[(G3-δ)(C2-δ)G3+δ]Lm.
Choose (u10(x),u20(x),u30)=(r1,m,r2). Then (r1,r2,r3) and (m,m,s) are the ordered upper and lower solutions of (22) and (3) (also of (23) and (3) and of (8) and (3)). Applying the same technique from [18, Theorem 4.1], we can prove that the solution (U1,u2,U3) of (22) and (3), the solution (u1,U2,u3) of (23) and (3), and the solution (p1,p2,p3) of (8) and (3) satisfy, respectively,
(29)limn→∞(U1(x,t+nT),u2(x,t+nT),U3(x,t+nt))=(U1δ(x,t),u2δ(x,t),U3δ(x,t)),limn→∞(u1(x,t+nt),U2(x,t+nT),u3(x,t+nT))=(u1δ(x,t),U2δ(x,t),u3δ(x,t)),limn→∞(p1(x,t+nT),p2(x,t+nT),p3(x,t+nT))=(u1T(x,t),u2T(x,t),u3T(x,t)).
It follows from (27) that uiδ≤uiT≤Uiδ(i=1,3) and that u2δ≤u2T≤U2δ. Similarly, choose (u10(x),u20(x),u30(x))=(m,r2,s); we can prove the other inequalities of (26).
Denote by ε1ui,δ1ui,ε2ui, and δ2ui the optimal bounds for the positive periodic solutions of problems (22) and (23), respectively, (see Lemma 2). Then
(30)δ2ui≤δui≤δ1ui,ε2ui≤εui≤ε1ui(i=1,3),δ1u2≤δu2≤δ2u2,ε1u2≤εu2≤ε2u2.
Moreover,
(31)ε2ui≤uiδ≤Uiδ≤δ1ui(i=1,3),ε1u2≤u2δ≤U2δ≤δ2u2,uiδ1<uiδ2,Uiδ1>Uiδ2
for 0<δ2<δ1<ε0. By using the dominated convergence theorem and the bootstrap arguments (see [18]), we have
(32)limδ→0+(U1δ,u2δ,U3δ)=(u1s,u2s,u3s),limδ→0+(u1δ,U2δ,u3δ)=(u1s,u2s,u3s)
uniformly for (x,t) on Ω¯×R, and (u1s,u2s,u3s) and (u1s,u2s,u3s) are the positive T-periodic solutions of (8).
From Theorem 1, we see that uiT≤uis and uis≤uiT(i=1,2,3). It follows from (26) and (32) that uis≤uiT,uiT≤uis(i=1,2,3). So,
(33)uis=uiT,uis=uiT(i=1,2,3).
Therefore, given that ε>0, by (32) and (33), there exists δ0∈(0,ε0) such that
(34)uiT-ε2<uiδ0≤Uiδ0<uiT+ε2(i=1,2,3).
Since (ai,bi,ci,gi)~(Ai,Bi,Ci,Gi), for the previous δ0, there exists Tδ0>0 such that, for t>Tδ0,
(35)Ai-δ0<ai<Ai+δ0,…,Gi-δ0<gi<Gi+δ0.
Denote by (u1,u2,u3) the solution of (1)–(3), and denote by (u1*,u2*,u3*), (U1,U2,U3), and (u1*,u2*,u3*) the solutions of problems (22), (8), and (23), respectively, which all satisfy the same initial conditions ui(x,Tδ0+1)=ui0(x)≥(≢)0(i=1,2,3). Analogous to (27), we have ui*≤ui≤ui*(i=1,2,3) for t>Tδ0. By (24), there exists T1>Tδ0+1 such that
(36)Uiδ0-ε2<ui*<Uiδ0+ε2(i=1,3),u2δ0-ε2<u2*<u2δ0+ε2
for t>T1. Similarly, by (25), there exists T2>Tδ0+1 such that
(37)uiδ0-ε2<ui*<uiδ0+ε2(i=1,3),U2δ0-ε2<u2*<U2δ0+ε2
for t>T2. Hence, by (34)–(37), we have
(38)uiT-ε<ui*≤ui≤ui*<uiT+ε
for t>max{T1,T2}. This completes the proof.
3. Global Stability
In order to get conditions of global stability for (1)–(3), we need the following result.
Lemma 7.
Let (u,v,w) be a T-periodic solution for the linear problem
(39)ut-d1Δu-Σaiuxi=M1(-Au+Bv-Cw),vt-d2Δv-Σbivxi=M2(Du-Ev),wt-d3Δw-Σciwxi=M3(Fu-Gw),∂u∂n=∂v∂n=∂w∂n=0on∂Ω×R,
where di, Mi, A, B, C, D, E, F, and G are positive smooth T-periodic functions on Ω¯×R and where ai, bi, and ci are smooth T-periodic functions. If
(40)B(DE)M+C(FG)M<A,
then u=v=w=0.
Proof.
Let (u,v,w) be a smooth T-periodic solution of (39), and let positive constants ε, m, k, and l be chosen so that
(41)M1L(A-Bk-Cl)≥ε,M2L(Ek-D)≥εk,M3L(Gl-F)≥εl,uM≤m,vM≤mk,wM≤ml.
Such choices are obviously possible because (40) holds.
It is easy to verify that m(1,k,l)e-εt and -m(1,k,l)e-εt are a pair of ordered upper and lower solutions of (39). Thus, -me-εt≤u(x,t)≤me-εt. This implies that u(x,t)=0 because u(x,t) is T-periodic in t.
Similarly, v=0=w. This completes the proof.
Lemma 8 (uniqueness).
Assume that (9) and (10) hold. If
(42)B1MB2Mεu2+A3MC2Mδu32(1+A3Lεu3)εu32<εu1(1+A3Lεu3)A1MA2Mδu1δu2,
then the problem (8) has a unique positive T-periodic solution.
Proof.
Let α=u1T/u1T-1, β=1-u2T/u2T and γ=u3T/u3T-1. By Theorem 1, we have
(43)αt-d1Δα-2d1U1TΣu1Txiαxi=g1u1Tu1T[-u1TαA1+A2u2Tβ1+A3u3T-A2A3u2Tu3Tγ(1+A3u3T)(1+A3u3T)],βt-d2Δβ-2d2U2TΣu2xiTβxi=g2u2Tu2T[B1u1Tα-u2TβB2],γt-d3Δγ-2d3U3TΣu3Txiγxi=g3u3T(C1+C2u1T)u3T[C2u1Tu3TαC1+C2u1T-u3Tγ],∂α∂n=∂β∂n=∂γ∂n=0on∂Ω×(0,∞).
It follows from Lemmas 2 and 7 and the conditions (9), (10), and (42) that α=β=γ=0. This completes the proof.
Theorem 9.
If all conditions of Theorem 4 and (42) are satisfied, then
(44)(u1T,u2T,u3T)=(u1T,u2T,u3T)~(u1,u2,u3)
for any positive solution (u1,u2,u3) of (1)–(3).
Proof.
By some elementary calculations, we know that Theorem 9 is an immediate corollary of Theorem 4 and Lemma 8.
Example 10.
Consider the following asymptotically periodic system:
(45)u1t-(2+sint)u1xx=(1+sin2(t+x)e-t2)u1(1-u1-u21+u3),u2t-(2-sint)u2xx=u2(1-(38+124cos2(t+x))u1-u2)in(0,1)×R,u3t-u3xx=u3(1-u35/12+(7/12)sin2(t-x)+u1),uix(0,t)=uix(1,t)=0onR(i=1,2,3).
It is not hard to verify that all conditions of Theorem 9 are satisfies. Thus, any positive solution of (45) asymptotically approach the unique positive periodic solution of the limit periodic system of (45).
4. Case a3=0
The following results are natural generalizations of the main results in [20] which can be proved in the similar way as to prove Theorems 4 and 9.
Theorem 11.
Assume the following.
d1, d2, A, B, C, D, E, and F are positive smooth T-periodic functions on Ω¯×R.
a, b, c, d, e, and f are positive smooth functions on Ω¯×R:
(46)(a,b,c,d,e,f)~(A,B,C,D,E,F),(FD)L>(CA)M,(BA)L>(ED)M.
Then the limit periodic system of (4)
(47)Ut-k1ΔU=U(A-BU-CV)inΩ×R,Vt-k2ΔV=V(D-EU-FV)inΩ×R,∂U∂n=0=∂V∂non∂Ω×R
has the positive T-periodic maxmini solution (UT,VT) and minimax solution (UT,VT). Moreover, if (u,v) is any positive solution of (4) with smooth initial value (u0,v0), then UT≺u≺UT and VT≺v≺VT. In addition, if
(48)(CA)M(ED)M[(FD)M-(CA)L][(BA)M-(ED)L]<(BA)L(FD)L[(FD)L-(CA)M][(BA)L-(ED)M],
then (47) has the unique positive T-periodic solution (U,V) and
(49)(u(·,t),v(·,t))~(U(·,t),V(·,t)).
Acknowledgments
The authors would like to thank the referees for their helpful comments. This work is supported by the China National Natural Science Foundation (Grants nos. 11061031 and 11261053) and the Fundamental Research Funds for the Gansu University.
RaiB.FreedmanH. I.AddicottJ. F.Analysis of three-species models of mutualism in predator-prey and competitive systems1983651135010.1016/0025-5564(83)90069-XMR714289ZBL0532.92025ZhengS. N.A reaction-diffusion system of a competitor-competitor-mutualist model1987124125428010.1016/0022-247X(87)90038-2MR883526ZBL0658.35053XuS.Global stability of a reaction-diffusion system of a competitor-competitor-mutualist model201115416171627MR2848977ZBL1233.35028PaoC. V.The global attractor of a competitor-competitor-mutualist reaction-diffusion system with time delays20076792623263110.1016/j.na.2006.09.027MR2343991ZBL1121.35021ChenW.PengR.Stationary patterns created by cross-diffusion for the competitor-competitor-mutualist model2004291255056410.1016/j.jmaa.2003.11.015MR2039069ZBL1060.35146LiM.LinZ.LiuJ.Coexistence in a competitor-competitor-mutualist model201034113400340710.1016/j.apm.2010.02.029MR2651775ZBL1201.35097FuS. M.GaoH. Y.CuiS. B.Global solutions for the competitor-competitor-mutualist model with cross-diffusion2008511153164MR2407998ZBL1164.35047TianC.LingZ.Turing pattern formation in a predator-prey-mutualist system20111263224323710.1016/j.nonrwa.2011.05.022MR2832962ZBL1231.35275TineoA.Asymptotic behavior of solutions of a periodic reaction-diffusion system of a competitor-competitor-mutualist model1994108232634110.1006/jdeq.1994.1037MR1270583ZBL0806.35095DuY.Positive periodic solutions of a competitor-competitor-mutualist model19969510431066MR1392094ZBL0858.35057PaoC. V.Periodic solutions of parabolic systems with nonlinear boundary conditions1999234269571610.1006/jmaa.1999.6412MR1695336ZBL0932.35111WangR.-N.XiaoT.-J.LiangJ.Asymptotic behavior of solutions for systems of periodic reaction-diffusion equations in unbounded domains200513281298MR2257700ZBL1122.35059ZhouL.FuY.Existence and stability of periodic quasisolutions in nonlinear parabolic systems with discrete delays2000250113916110.1006/jmaa.2000.6986MR1893882ZBL0970.35004WangY.YinJ.Periodic solutions of a class of degenerate parabolic system with delays20113801576810.1016/j.jmaa.2011.02.044MR2786183ZBL1217.35014GopalsamyK.Exchange of equilibria in two-species Lotka-Volterra competition models198224216017010.1017/S0334270000003659MR671810AlvarezC.LazerA. C.An application of topological degree to the periodic competing species problem198628220221910.1017/S0334270000005300MR862570AhmadS.Convergence and ultimate bounds of solutions of the nonautonomous Volterra-Lotka competition equations1987127237738710.1016/0022-247X(87)90116-8MR915064ZBL0648.34037AhmadS.LazerA. C.Asymptotic behaviour of solutions of periodic competition diffusion system198913326328410.1016/0362-546X(89)90054-0MR986448ZBL0686.35060TineoA.Existence of global coexistence state for periodic competition diffusion systems199219433534410.1016/0362-546X(92)90178-HMR1178407ZBL0779.35058PengQ. L.ChenL. S.Asymptotic behavior of the nonautonomous two-species Lotka-Volterra competition models19942712536010.1016/0898-1221(94)90085-XMR1284130ZBL0798.92023FuS.CuiS.Persistence in a periodic competitor-competitor-mutualist diffusion system2001263123424510.1006/jmaa.2001.7612MR1865278ZBL0995.35008DimbourW.N'GuérékataG. M.S-asymptotically ω-periodic solutions to some classes of partial evolution equations2012218147622762810.1016/j.amc.2012.01.029MR2892728GaoH.WangK.WeiF.DingX.Massera-type theorem and asymptotically periodic logistic equations2006751268128310.1016/j.nonrwa.2005.11.008MR2260914ZBL1162.34325HenríquezH. R.Asymptotically periodic solutions of abstract differential equations20138013514910.1016/j.na.2012.10.010MR3010761ZBL1266.34100WangJ.FeckanM.ZhouY.Nonexistence of periodic solutions and asymptotically periodic solutions for fractional differential equations201318224625610.1016/j.cnsns.2012.07.004MR2971825ZBL1253.35204WeiF.WangK.Asymptotically periodic solution of N-species cooperation system with time delay20067459159610.1016/j.nonrwa.2005.03.019MR2235220WeiF.WangK.Uniform persistence of asymptotically periodic multispecies competition predator-prey systems with Holling III type functional response2005170299499810.1016/j.amc.2004.12.040MR2175261ZBL1084.92036YeQ. X.LiZ. Y.1990Beijing, ChinaScience PressMR1113267