Global Stabilization of a Single-Species Ecosystem with Markovian Jumping under Neumann Boundary Value via Laplacian Semigroup

: By applying impulsive control, this work investigated the global stabilization of a single-species ecosystem with Markovian jumping, a time delay and a Neumann boundary condition. Variational methods, a ﬁxed-point theorem, and Laplacian semigroup theory were employed to derive the unique existence of the global stable equilibrium point, which is a positive number. Numerical examples illuminate the feasibility of the proposed methods.


Introduction
As pointed out in [1], the following logistic system has been widely discussed and studied due to its importance in the development of ecology: Here, Z(t) is the population's quantity or density at the time t, and K and R > 0 are the intrinsic growth rate of the environmental capacity and the population. Because all the solutions of nonlinear ecosystems are difficult to provide accurately, people pay more attention to the long-term dynamic trend of an ecosystem, i.e., the long-term trend of population density (see, for example, [1][2][3][4][5]). People especially want to know whether the population will tend to a positive constant after a long time, which is related to the final long-term existence of the population. For example, the long-term behavior of the following random single-species ecosystem was studied in [2]: Animal populations will inevitably spread because of climate, foraging and random walking. Hence, the reaction-diffusion ecological models well simulate a real ecosystem. Particularly, reaction-diffusion ecosystems were studied in [6,7]. For example, in [8], a single-species Markovian jumping ecosystem with diffusion and delayed feedback under a Dirichlet boundary value was investigated: Markov systems often occur in various engineering technologies (see, for example, [9][10][11]). Particularly, a Markovian jumping delayed feedback model well reflects the influence of stochastic factors on time delays in the changes of populations, such as weather, temperature, humidity, ventilation status, etc. However, the case of a single-species ecosystem with a Neumann boundary value is seldom researched. In fact, a Neumann zero boundary value model well simulates the biosphere boundary without population migration. For example, freshwater fish do not enter the sea through rivers. Inspired by some ideas or methods of the related literature , we are willing to study the global stabilization of a Markovian jumping delayed feedback diffusion ecosystem with a single species equipped with the Neumann zero boundary value.
The main contributions are as follows.
• The uniqueness proof of the positive equilibrium solution is presented in this paper, while it was given in previous work that only involved the existence of the positive equilibrium solution. • In the case of a single-species ecosystem with impulses, it is the first study using a Laplacian semigroup to globally stabilize the ecosystem. • A numerical example is designed to illuminate the advantages of Theorem 2 against [22] (Theorem 4.2), as a result of reducing the algorithm's conservatism.

Preparatory Work
Consider the reaction-diffusion delayed ecosystem: Here, u(t, x) is the population density at time t and space point x. a > 0 and b > 0 are described similarly as those of [8]. τ(t) ∈ [0, τ], and Γ(s, x) is the bounded initial value function on [−τ, 0] × Ω. For convenience, c(r(t)) is denoted simply by c r for r(t) = r ∈ S.
In addition, due to the limited resources of nature, the population density should have an upper limit. At the same time, the population density should also have a lower limit because a low population density does not allow male and female animals to meet in the vast biosphere and reproduce. Hypothesis 1 (H1). There exist two positive constants N 1 and N 2 and a decimal k 0 ∈ (0, 1) with Remark 1. The boundedness assumption in (H1) brings innovations. It ensures that the initial value maintains a certain distance from the upper and lower bounds, so that an impulse with an appropriate frequency and intensity can ensure that the dynamic behavior of the system with such an initial value will not exceed the bounds.
, and u * (x) satisfies (H1), and Now, it is easy to conclude from Definition 1 that u * ≡ a b is a stationary solution of the system (5). Moreover, letting Here, the positive solution u * ≡ a b of the ecosystem (5) corresponds to the zero solution of the system (8). Thus, the stabilization of the above-mentioned zero solution will be investigated below. Furthermore, employing an impulse control on the natural ecosystem (5) or (8) results in whose zero solution corresponds to the stationary solution u * ≡ a b of the following system: Here, we assume that 0 < t 1 < t 2 < · · · , and each t k (k ∈ N) represents a fixed impulsive instant. Additionally, lim
In this paper, the following condition is also required: Hypothesis 2 (H2). There are two constants M > 0 and γ > 0 such that e t∆ e t∆ w (see [23]). Lemma 1. (see, for example, [14]) . Ω ⊂ R m is a bounded domain, and its smooth boundary ∂Ω is of class where λ 1 is the least positive eigenvalue of the following Neumann boundary problem: where ν is the external normal direction of ∂Ω.

Lemma 2. ([36]
) If f is a contraction mapping on a complete metric space H, there must exist a unique point u ∈ H, satisfying f (u) = u.

Main Result
Firstly, the unique existence of the stationary solution of the ecosystem (5) should be considered. Moreover, the unique stationary solution of the ecosystem should be positive. Based on these two points, we present the following unique existence theorem: If, in addition, the following condition is satisfied: then the positive solution u * is the unique stationary solution of the system (5).
Thus, Definition 1 yields that u * > 0 defined in Theorem 1 is the unique stationary solution of the system (5).
Below, we claim that u * is the unique stationary solution of the ecosystem (5). Indeed, if u * and v * (x) are two different stationary solutions of the system (5), then Poincare inequality and the boundary condition yield The condition (12), Definition 2 and the continuity of u * and v * lead to which contradicts the inequality (13). This completes the proof.

Remark 3.
As far as our knowledge, Theorem 1 is the first theorem to provide the unique existence of the stationary solution of a single-species ecosystem under a Neumann boundary value.

Remark 4.
This paper provides the unique existence of a stationary solution of a reaction-diffusion system. However, there are many previous articles related to reaction-diffusion systems that only involve the existence of the equilibrium point. For example, in [14], only the unique existence of the constant equilibrium point of a reaction-diffusion system with a Neumann boundary value was given, but the stationary solutions of a reaction-diffusion system may include the non-constant stationary solutions. Because the solution u(t, x) of a reaction-diffusion system involves not only the time variable t but also the space variable x, its stationary solution should be u * (x), independent of the time variable t. Obviously, u * (x) must not be a constant equilibrium point, which may be dependent upon the space variable x. Thereby, it is not inappropriate to prove that the equilibrium point is the unique constant equilibrium point in [14]. It must be proved that it is the unique stationary solution, just like that of this paper. A similar example can also be found in [12]. Note that the system (10) has the same elliptic equation as that of the system (5), and hence, each stationary solution of the system (5) is that of the system (10), and vice versa. i.e., Theorem 1 shows that u * ≡ a b is also the unique stationary solution of the system (10). Next, the global stability of the stationary solution u * ≡ a b should be investigated.
Theorem 2. Set p 1. Suppose all the conditions of Theorem 1 hold. Assume, in addition, the condition (H2) holds, and 0 < r < 1, r(t) = r ∈ S; then, u * ≡ a b of the system (10) is globally exponentially stable in the pth moment; equivalently, the null solution of the system (9) is globally exponentially stable in the pth moment, where Proof. Let the normed space H be the functions space consisting of functions g(t, x) : , where g satisfies: (A1) g is pth moment continuous at t 0 with t = t k (k ∈ N) ; (A2) for any given x ∈ Ω, lim   (A3) g(s, x) = ξ(s, x), ∀ s ∈ [−τ, 0], x ∈ Ω ; (A4) e αt g(t) p → 0 as t → +∞, where α is a positive scalar with 0 < α < qγ. It is easy to verify that H is a complete metric space equipped with the following distance: Construct an operator Θ such that, for any given U ∈ H, Below, we want to prove that Θ : H → H, and it takes four steps to achieve the goal.
Step 1. It is claimed that, for U ∈ H, Θ(U) must be pth moment continuous at t 0 with ] means the boundedness of U, and let δ be a scalar small enough; a routine proof yields that if which proves the claim. Then, (A1) is verified.
Step 4. Verifying (A4), i.e., for U ∈ H, verifying Indeed, e q(t−s)∆ c r U(s, x)ds p Moreover, The Holder inequality yields On the other hard, e αt U i (t) p → 0 means that, for any ε > 0, there exists t * > 0 such that all e αt U i (t) p < ε. Therefore, which, together with the arbitrariness of ε, means that Combining (20)-(28) yields (19). It follows from the above four steps that Finally, we claim that Θ is a contractive mapping on H. Indeed, for any U, V ∈ H, the Holder inequality and (H2) yield Similarly, Assume t j−1 < t t j ; then, the definition of the Riemann integral b a e s ds yields It follows from (30) Now, the definition of r implies that Θ : H → H is contractive such that there must exist the fixed point U of Θ in H, which means that U is the solution of the system (9), satisfying e αt U p → 0, t → +∞ so that e αt u − u * p → 0, t → +∞ . Therefore, the zero solution of the system (9) is the globally exponential stability in the pth moment; equivalently, u * ≡ a b of the system (10) is the globally exponential stability in the pth moment.

Remark 5.
To the best of our knowledge, this is the first paper to employ impulsive control and the Laplacian semigroup to globally stabilize a single-species ecosystem. Remark 6. This paper reports the global stability of a single-species ecosystem, while the stability in [3] did not involve the global one. This means that the stability in [3] depends heavily on the choice of initial value, while the global stability does not need such a choice. On the other hand, Equation (5) involves the space state, while the models in [3] did not involve the spatial location. In fact, population migration has a great impact on population stability, so the spatial state should be considered in the ecosystem model.

Example 2.
This example uses all the data provided in Example 1. Assume, in addition, p = 1.005, M k ≡ 1.02, µ = 5; obviously, the condition (H2) holds in Example 1. Moreover, we can obtain, by direct calculations, that and, hence, the condition (14) is satisfied.
Thereby, Theorem 2 yields that the null solution of the system (9) is globally exponentially stable in the pth moment; equivalently, u * ≡ 0.423 of the system (10) is globally exponentially stable in the pth moment (see Figures 1 and 2).

Remark 7.
Obviously, Example 2 illuminates that Theorem 2 is better than [22] (Theorem 4.2) due to reducing the conservatism of the algorithm, because the pulse intensity in Theorem 2 does not require the pulse intensity to be less than 1 but allows it to be greater than 1, while the latter requires that the pulse intensity is less than 1 (see [22] (Theorem 4.2)).

Remark 8.
In Example 1 and Example 2, it follows from 0.3 = N 1 u N 2 = 0.523 and u * = 0.423 that 0 |U| 0.123. A computer simulation of the dynamics of the state U(t, x) of the system (9) illuminated the feasibility of Theorems 1 and 2 (see Figures 1 and 2). In addition, (H1) is the only common condition in much related literature (see, for example, [37] (Definition 1)).

Conclusions
In this paper, there are some improvements on mathematical methods, because it is the first paper to employ fixed-point theory, variational methods and a Laplacian semigroup to obtain the unique existence of the globally stable positive equilibrium point of a singlespecies Markovian jumping delayed ecosystem. Numerical examples are provided to show the feasibility of the artificial management of nature by way of impulse control.