Discrete Dynamics in Nature and Society

Volume 2016 (2016), Article ID 5036305, 10 pages

http://dx.doi.org/10.1155/2016/5036305

## The Modeling and Control of a Singular Biological Economic System with Time Delay in a Polluted Environment

^{1}School of Automation, Nanjing University of Science and Technology, Nanjing, Jiangsu 210094, China^{2}School of Science, Shenyang University of Technology, Shenyang, Liaoning 110870, China^{3}State Key Laboratory of Synthetical Automation for Process Industries, Northeastern University, Shenyang, Liaoning 110819, China

Received 16 July 2016; Accepted 20 September 2016

Academic Editor: Daniele Fournier-Prunaret

Copyright © 2016 Yi Zhang et al. 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

This paper brings up the idea of a biological economic system with time delay in a polluted environment. Firstly, by proper linear transformation and parametric method, the singular time-delay systems are transformed to differential time-delay systems. Then, using center manifold theory and Poincare normal form method, the direction of Hopf bifurcation and the stability and period of its periodic orbits are analysed. At last, we have performed numerical simulation to support the analytical results.

#### 1. Introduction

Environmental pollution has been increasingly influencing the biological systems. In order to investigate the development and dynamics of population of the biological systems, it is necessary to consider the factor of pollution when establishing a mathematical model. In addition, delay is also a kind of common phenomenon in reality and it has great influence on the dynamic behavior of system. Therefore, the delay differential equations are needed to describe the system when the influence of time delay is considered. Time delay can lead to the imbalance of the system and the emergence of a variety of bifurcations, among which Hopf bifurcation is the most common. The properties of Hopf bifurcation consist of the stability of the periodic solutions, the direction of bifurcations, the period, and so forth. In recent years, the theory of delay system has gradually been generalised to many important fields by domestic and foreign scholars, including the applications in circuit communiment-system [1], electrodynamics [2], optical [3], ecological-system [4], and economics [5]. Many research findings on biological applications also emerged, such as the analysis of the stability of a class of stochastic system with time delay [6], investigation on nonautonomous competitive Lotka-Volterra systems with infinite delay [7], researching on dynamic behavior of a class of prey-predator model with time delay in a polluted environment [8], the analysis and control of a class of singular prey-predator model with discrete delay [9] which studies the prey-predator system with commercial harvesting and double time delays, and the dynamic behavior analysis and optimal control of a class of economic model with stage structure and pregnancy delay [10].

There are many kinds of research methods for delay differential systems. Among them, the most commonly used ones [11] are the center manifold method and the Poincare normal form method. Being always an important mathematical means to investigate the bifurcation problems with parameter and the qualitative theory of differential equations, more attention has been paid to the Poincare normal form method for a long time, home and abroad. In [12], the author lays the foundation of the center manifold standard method by combining the normal form theory and the center manifold theorem, and the method was used on the investigation of Hopf bifurcation. When it comes to related properties of the Hopf bifurcation, the center manifold standard method is usually used to reduce the dimension of high-dimension system, which isolates the asymptotic behaviors of complex systems, so that we can investigate the original system in a center manifold of low dimension, which is much simpler. This paper takes a singular biological economy system with time-delay in a polluted environment and analyses it using the stability theory of singular system, the theory of economic system, the theory of the Hopf bifurcation of delay differential system, and so forth.

#### 2. Model Formulation

A single-creature model with stage structure is investigated in [5, 13]:where and are the densities of the immature number of creatures and mature number of creatures at time , denotes the birth rate of immature creatures, , are the death rates of immature creatures and mature creatures, denotes the conversion rate from immature creatures into mature creatures, and denotes the intraspecific effect coefficient. All coefficients are positive.

A single-creature model in the polluted environment is investigated in [5]:where is the creatures density, is the concentration of environment pollutants, denotes the intrinsic growth rate when there is no pollution, denotes the capacity of the environment, can be interpreted as the measuring response function of the reduction of creatures because of the pollution factor, denotes the amount of pollutants that are inputted by the outside, and can be interpreted as the reduction of pollutant concentration because of other factors. Assume that endotoxin excretion rate and purification rate are relatively small in an organism body, and thus it can be neglected.

Considering the need of a period of time when the immature creatures change into mature creatures, based on system (1) and system (2), the following system is proposed:where is the capture capability of mature creatures at time , denotes the unit price, denotes the unit cost, and denotes the economic profit. is the total revenue, and is the total cost. All the parameters are positive [14, 15].

#### 3. Stability Analysis

Theorem 1. *The positive equilibrium of system (3) is locally asymptotically stable when ; the positive equilibrium of system (3) is unstable when ; system (3) undergoes a Hopf bifurcation at the positive equilibrium when .*

*Proof. *System (3) will be investigated in this chapter [16].

Considering system (3), letFrom (4), the following root can be obtained: where is decided by the following equation: The discriminant of the equation isWhen , the parameters of system (3) meet conditions , and then is the unique positive equilibrium point of system (1), whereWhen , system (3) has two equilibrium points , , whereThe stability of the unique positive equilibrium point of system is investigated as an example; the stability of the other equilibrium points can be investigated in the same way.

In order to make the research more convenient, system (3) can be written as [17]where .

In order to investigate the local stability of the positive equilibrium point, make the following transformation on system (10):Then,By generating system (3), the following system can be obtained:In order to derive the formula determining the properties of the positive equilibrium of system (13), we consider local parametric of the fourth equation of system (13) as literatures [18], which is given as follows:where, and is a continuous map.By generating system (13), the following system can be obtained:The linearized system of parametric system (17) at can be given as follows:The characteristic equation of system (18) is as follows:where , , , , and let to be a root of (19), , in (19); thenSeparate the imaginary part and real part from (20) and the following equations can be obtained:Then, the following equation can be obtained by calculatingLet ; then (22) can be changed as the following form:When , , (22) characteristic has no real root; when , , (22) characteristic has one real root , in whichWhen , , (22) characteristic has one real root , in whichWhen , , (22) characteristic have two real roots , in whichAccording to (21), if , , can be written as the following form: Then, (19) has a pair of pure imaginary roots at ; assuming that (19) has a solution , when , , in (19), taking the derivative of with respect to in (19),whenWhen , , ( is the minimum of )DenoteThen, we can getAccording to the above analysis, we can prove the theorem.

#### 4. Direction and Stability of the Hopf Bifurcation

Theorem 2. *The properties of Hopf bifurcation are determined by (70), the detailed contents are as follows:*(1)*The direction of Hopf bifurcation: the Hopf bifurcation is supercritical (resp., subcritical) when resp., and the bifurcation periodic solutions exist when *(2)*The stability of the bifurcating periodic solutions: the bifurcation periodic solutions are stable (resp., unstable) if (resp., ).*(3)*The period of the bifurcating periodic solutions: the period increases (resp., decreases) if (resp., ).*

*Proof. *Considering system (19), let , , ; then system (3) can be written asin which , , .By Riesz representation theorem, there exist in interval andin which is a bounded variation function. Letfor , definingThen, (33) is equivalent toin which , .

For ,It is easy to calculate that and are the eigenvalues of and ; then we can obtain the feature vectors of and of and ; let ; according to (34) and (37), we can getThen,By the same way, let ; we can getBy (37), we can getdue to ; then,Next, the coordinates to describe the center manifold at will be calculated. Assuming that is the root of (33) at , defineEquation (47) satisfied the following function in center manifold :in which , are the local coordinates of center manifold that, in the directions of , if is a real root of (39), when , , the following equation can be obtained:The equation can be written asin whichdue to , ; thenAccording to (51), we can get thatcompared with the parameter of (51); then we get In order to calculate , we need to get and by (39) and (41):in whichBy comparing the parameters, thenAccording to (51) and (55), when ,Considering the parameters of (56), thenIn the same way, Calculate vectors and and let , by (51), (56), and (57) and we can getBy the definition of (57), then Put (60) and (61) to (65), and we can get Let ; according to (66), we can get In the same way, put (61) and (63) to (65), we can getLet ; the root of (67) can be obtained According to , , (60) and (61), can be easily obtained. Furthermore, we can see that each in (54) is determined by parameters and delays. Thus, we can compute the following quantities.

We have the following results:This completes the proof.

*5. Numerical Simulation*

*By selecting some related data from China environment protection database and doing the appropriate treatment, the following parameters can be obtained [19, 20]: , , , , , , , , , , , and .*

*Then, system (3) can be written as follows: is the unique positive equilibrium of system (70). , , , and is the bifurcation parameter of system (70). When , the dynamical responses of system (70) are shown by Figure 1, system (70) is stable at , and the population and economic profits develop sustainably in this case; when , the dynamical responses of system (70) are shown by Figure 2, system (70) is unstable, and the population and economic profits cannot develop sustainably in this case; Figure 3 shows that the dynamic behavior of the population changes with time delay and the Hopf bifurcation exists when .*