Journal of Optimization

Volume 2019, Article ID 9512879, 12 pages

https://doi.org/10.1155/2019/9512879

## Theoretical Analysis of an Imprecise Prey-Predator Model with Harvesting and Optimal Control

Department of Applied Mathematics with Oceanology and Computer Programming, Vidyasagar University, Midnapore, West Bengal-721101, India

Correspondence should be addressed to Anjana Das; moc.liamg@ajna.8891

Received 16 May 2018; Accepted 10 December 2018; Published 21 January 2019

Academic Editor: K. F. C. Yiu

Copyright © 2019 Anjana Das and M. Pal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In our present paper, we formulate and study a prey-predator system with imprecise values for the parameters. We also consider harvesting for both the prey and predator species. Then we describe the complex dynamics of the proposed model system including positivity and uniform boundedness of the system, and existence and stability criteria of various equilibrium points. Also the existence of bionomic equilibrium and optimal harvesting policy are thoroughly investigated. Some numerical simulations have been presented in support of theoretical works. Further the requirement of considering imprecise values for the set of model parameters is also highlighted.

#### 1. Introduction

The eternal relationship between prey and predators is one of the major topics to be discussed in recent science. Scientists from various fields are currently engaged in finding out different interactions among prey populations and predator populations. With the help of mathematical modeling, one can describe the strong and competitive relationship between these two types of creatures. However the discussion on this fascinating topic, with the help of some mathematical tool, was started during the first quarter of the twentieth century, thanks to the age-breaking works of Lotka [1] and Volterra [2]. Influenced by those works, researchers are still engaged in theoretical study of the ecological system with the help of mathematical modeling. In their book, Kot [3], Britton (2003) described some ecological phenomena on various ecological interactions including prey-predator interactions. Smith [4] has demonstrated various aspects in theoretical ecology with the help of some basic mathematical models. May [5] has also considered and analyzed some other types of ecological systems including prey-predator dynamics, with the help of some sophisticated mathematical models, which are comparatively complex in nature. Some other research works on theoretical ecology including predator-prey dynamics can be found, like Cushing [6], Hadeler and Freedman [7], Chen et al. [8], Kar [9, 10], Kar et al. [11], Chakraborty et al. [12], and references therein.

Harvesting is however a common and quite natural phenomenon. In fishery harvesting is used frequently as the biological resources are mostly renewable resources. On an exploited fishery system with interacting prey and predator species, researchers are considering harvesting on either prey species or predator species or harvesting on both prey and predator species. Martin and Ruan [13] discussed the dynamics of harvesting of prey populations whereas, in his article, Kar [9] describes phenomenon on selective harvesting on a prey-predator system. Further harvesting in predator species or both of prey and predator species can be found in literature also (Kar and Pahari [14], Zhang and Zhang [15], Jana et al. [16–18], Pal et al. [19], Walters et al. [20], Liu and Zhang [21], etc.). From a bioeconomic point of view, harvesting on a species should be in a balance to both keep the resource live and keep fishermen in profitable mode. In his two books Clark [22, 23] describes different harvesting policies in some realistic ecological systems with optimal outcome.

In this regard, in our present article, we consider a prey-predator type ecological system with harvesting on both the species. However, till now most of the models are proposed by considering only precise set of parameters but the natural world may not be precise every time. In many situations at experimental field like birth and death rate of different individuals of the same species, interaction between two different species, etc., may be imprecise. For this purpose introduction of fuzzy sets (Zadeh [24]) is now considered as a revolutionary work. However, consideration of interval-valued parameters was due to the broad application of fuzzy sets. The imprecise parameters set may not always belong to the interval but they may belong to any interval of positive number. Hence interval-valued parameter set of an imprecise mathematical model would be regardlessly better from a realistic point of view.

The rest of the paper is organized in the following manner: In Section 2, we put some preliminaries on interval value numbers. In Section 3, on the basis of some realistic assumptions we formulate our predator-prey system and then convert it to an imprecise parametric system whose dynamical behavior is thoroughly discussed in Section 4. Section 5 is devoted to discussing the existence of bionomic equilibrium, and optimal harvesting policies are studied in Section 6 keeping harvesting parameter as the control variable. In Section 7, we validate our theoretical results through some numerical simulation works, and in the last section we present some key findings.

#### 2. Preliminaries

Here we give the definition of interval numbers with some operations. We use interval-valued function in lieu of interval number.

*Definition 1 (interval number). *We denote interval number as and define it as where is called the set of all real numbers and , are the lower and upper limits of the interval number, respectively.

A real number can also be used in form of interval number as .

The basic operations between any two interval numbers are as follows:(i).(ii).(iii) where is a real number.(iv), .(v).

*Definition 2 (interval-valued function). *For an interval the interval-valued function can be created as for .

#### 3. Predator-Prey Model with Harvesting with Imprecise Parameter

##### 3.1. Crisp Model

In this article, we consider only two species, namely, prey species and predator species. Let denote the prey biomass and the predator class at any time . Let the prey population grow logistically with intrinsic growth rate and environmental carrying capacity . Also let the predator attack the prey in the predation rate following the mass action law. Thus the differential equation for the prey population becomes Let be the conversion factor from the prey population to the matured predator population, be the natural death rate, and be the intraspecific competition rate for the predator populations (due to Ruan et al. [25]). Then the differential equation of the predator population reduces toHere are all positive parameters.

Next if we consider that both the species are harvested, this is carried out on assuming the demand in the market of both species (prey and predator). Taking as the harvesting effort for both species and & as the catchability coefficient of the prey species and predator species, respectively, then our system reduces tosubject to the initial conditions

##### 3.2. Fuzzy Model

The environment and other factors including temperature and food habits caused the parameters to be imprecise. So they should be taken as interval number rather than a single value. Let be the corresponding interval numbers for , respectively. Then the prey-predator model with combined harvesting effort becomeswhere

For the interval number we consider the interval-valued function for . Similarly taking the other interval numbers in the same way as function form, we get the model assubject to the initial conditionsHere , where the value of depends on the underlying environment.

#### 4. Dynamical Behavior

In this section, we describe a thorough dynamical behavior of the proposed model system. To do so, we first check the positivity of the solutions of crisp system and uniform boundedness of the solution of the same system. Now, it can also be concluded that uniform boundedness and positivity in the solutions also hold for the corresponding fuzzy systems, if these things hold in crisp system.

##### 4.1. Positivity

First we consider the corresponding crisp system in following form.Now, on integration, we have, from above system of equations,andHence from above, two expressions related to two state variables will always be positive. Thus the solution of corresponding crisp problem will be nonnegative and so the solution of the corresponding fuzzy system will also be nonnegative.

##### 4.2. Uniform Boundedness

In this section we now study the uniform boundedness of the proposed imprecise system. Now from the first expression of system (5), we haveNow, by simple mathematics, it can be concluded that has maximum value , which is obtained for Thus from above, we haveNow Integrating both sides of the above inequality and then applying the theory of differential inequality due to ( see Birkhoff and Rota [26]), we haveNow on letting , we haveHence the biomass density of prey population is uniformly bounded with an upper and lower limit and , respectively.

Next we are targeting to show that the biomass of predator population is uniformly bounded. In this regard from (5), we haveThus similarly to the above, it is to be claimed that the right hand side of the above expression has maximum value at and this maximum value is .

Thus proceeding in the same way as prey populations and with the help of Birkhoff and Rota [26], we can writeNow on letting , we haveSo the biomass density of predator populations is also uniformly bounded with lower and upper bound, respectively, and

Hence the biomass density of both the population species is uniformly bounded.

##### 4.3. Existence of Equilibria

The equilibrium points of this system are given below.

(1) Trivial equilibrium: .

(2) Axial equilibrium: [where ] exists if .

(3) Interior equilibrium: where

The interior equilibrium exists

ifandhold ifwhereand

##### 4.4. Local Asymptotic Stability

In this section we state and prove the local asymptotic stability criteria at different equilibrium points. Also the corresponding conditions for which the system is stable at different equilibria are given below.

*Case 1. *For trivial equilibrium the variational matrix at is given by the following.Therefore, the eigenvalues are given by .

Here ; then is asymptotically stable if , i.e., if which implies .

In the next theorem, we state the stability criteria of trivial equilibrium point

Theorem 3. *Trivial equilibrium point of the system is locally asymptotically stable if holds.*

*Case 2. *At axial equilibrium the variational matrix iswhereand . Then the eigenvalues of the characteristic equation of are and . The first one of them is negative since . Now is asymptotically stable if the second one is negative, i.e.,which implies

In the next theorem we will state the local asymptotic stability criteria of the axial equilibrium or the predator free equilibrium

Theorem 4. *The axial equilibrium is locally asymptotically stable if*

In this condition the trivial equilibrium becomes unstable.

*Case 3. *The variational matrix for interior equilibrium is written below.The characteristic equation of is given by whereand Here and since and . Then the values of are negative.

Therefore, The system is locally asymptotically stable at and we state this criteria in the following theorem.

Theorem 5. *The interior equilibrium of the system exists and is locally asymptotically stable ifwhereand*

##### 4.5. Global Stability

Here we will discuss the global asymptotic stability criteria of the system around its interior equilibrium point. In next theorem we study the criteria.

Theorem 6. *The interior equilibrium of the system is globally asymptotically stable provided it is locally asymptotically stable there.*

*Proof. *A Lyapunov function is constructed here as follows where is suitable positive constant to be determined in the subsequent steps.

Taking derivative with respect to along the solutions of the system, we haveNowThenIf we consider , then reduces to the following.It is seen from the above that .

That is, the system is globally stable around its interior equilibrium .

#### 5. Bionomic Equilibrium

In this section we study the bionomic equilibrium of the competitive predator-prey model. Here we consider the following parameters: (1) : fishing cost per unit effort, (2) : price per unit biomass of the prey, (3) : price per unit biomass of the predator. The net revenue at any time is given byThe interior equilibrium point of the system is on the line given belowThis biological equilibrium line meets x-axis at and y-axis at , whereand It is seen that always but is feasible if

The ‘zero-profit line’ is given byEquation (6) together with the above condition represents the bionomic equilibrium of prey-predator harvesting system.

For the points on the equilibrium line where , the fishery becomes useless. Because it cannot produce any positive economic revenue.

These three cases may arise in bionomic equilibrium.

*Case 1. *When fishing or harvesting of predator species is not possible, then gives that .

*Case 2. *When harvesting of prey is not possible, then gives that with .

*Case 3. *When the bionomic equilibrium is at a point where both and , then the fishing of prey and predator is possible. Here where Since and , then the following two conditions hold.where Then we conclude the bionomic equilibrium shorty in the following theorem.

Theorem 7. *The bionomic equilibrium () always exists, exists when , and exists when conditions (49), (50) and (52), (53) hold simultaneously.*

#### 6. Optimal Harvesting Policy

Here both prey and predator populations are considered as fish populations. The optimal net profit is obtained from fishing. We discuss in this section the optimal harvesting policy. We consider the profit gained from harvesting taking the cost as a quadratic function and focusing on the conservation of fish population. The price assumed here is inversely proportional to the available biomass of fish (prey and predator); i.e., if the biomass increases, the price decreases (see Chakraborty et al. (2011)). Let be the constant harvesting cost per unit effort and and be, respectively, the constant price per unit biomass of the prey and predator. Now our target is to get the maximum net revenues from fishery. Then the optimal control problem can be created in the following way:subject to the system of differential equations (6) and the initial conditions (7). and are economic constants and is the instantaneous discount rate.

Here the control is bounded in and our object is to find an optimal control such that where is the control set defined byHere the convexity of the objective functional with respect to the control variable along with the compactness of the range values of the state variables can be combined to give the existence of the optimal control Now the optimal control can be found by using Pontryagin’s maximum principle (Pontryagin et al. (1962)). To optimize the objective functional , we construct the Hamiltonian of the system as follows:Here the variables and are adjoint variables and the transversality conditions are as follows.First we use the optimality condition to obtain the optimal effort which is as follows: The adjoint equations are and, therefore, we have the following theorem regarding the optimal value of the harvesting effort.

Theorem 8. *There exists an optimal control , corresponding to the optimal solutions for the state variables as and such that this control optimizes the objective functional over the region Moreover, there exist adjoint variables and satisfying the first order differential equations given in (60) with the transversality conditions given in (58), where, at the optimal harvesting level, the values of the state variables and are, respectively, and *

#### 7. Numerical Simulation

In this section, we analyze our mathematical model through some simulation works. The main difference of our proposed model compared to other models of the same type is the consideration of interval-valued parameters instead of fixed-valued parameters. Inclusion of the parameter assumes the value corresponding to the parameters of the system as an interval. In this regard we first analyze the importance of considering the parameter in Figure 1. For simulation purpose we consider the following parametric values: , , , , , , , For different parametric values of , we have obtained various types of dynamical behavior of the proposed prey-predator system. From Figure 1, it can be said that lower values of make the system unstable at the interior equilibrium point whereas the higher values of gradually make the system locally asymptotically stable around the interior equilibrium point. As the numerical value of increases, the instability solutions slowly become stable (unstable branches at are less than the number of unstable branches at , but still at the system is unstable but asymptotically stable at a higher value ()).