# Finite-Time Stability Analysis and Control for a Class of Stochastic Singular Biological Economic Systems Based on T-S Fuzzy Model

Authors: Shuangyun Xing; Qingling Zhang; Yi Zhang
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

This paper studies the problem of finite-time stability and control for a class of stochastic singular biological economic systems. It shows that such systems exhibit the distinct dynamic behavior when the economic profit is a variable rather than a constant. Firstly, the stochastic singular biological economic systems are established as fuzzy models based on T-S fuzzy control approach. These models are described by stochastic singular T-S fuzzy systems. Then, novel sufficient conditions of finite-time stability are obtained for the stochastic singular biological economic systems, and the state feedback controller is designed so that the population (state of the systems) can be driven to the bounded range by the management of the open resource. Finally, by using Matlab software, numerical examples are given to illustrate the effectiveness of the obtained results.


---

## Body

## 1. Introduction

At present, mankind is facing the problems of shortage of resources and worsening environment. So there has been rapidly growing interest in the analysis and modeling of biological systems. From the view of human needs, the exploitation of biological resources and harvest of population are commonly practiced in the fields of fishery, wildlife, and forestry management, which has already been tackled in the classical literature (see [1–5]). In these works, the main objective was to maximize a utility function representing the sustainable economic rent (i.e., net economic income). The results may not be realistic due to assumptions under which the models are studied. Indeed, the validity of the assumption stipulating that price is constant seems to be limited to a small-scale fishery, wildlife, or forestry. Also, the assumption of free access to valuable resources must be discarded because it inevitably leads to an overexploitation. For instance, the Moroccan octopus crisis arising between 1996 and 1997 generated a biological (overexploitation of the resource) and economic instability (raising the resource price).To avoid these situations, Jerry and Raïssi [6] have studied the biological and economic stabilizability by adding a dynamic prevailing on the market in order to take into account resource price variations based on a standard open-access fisheries model (see [1, 7–10]); the model was given by
(1)dx(t)=[rx(t)(1-x(t)k)-Ex(t)]dt,dp(t)=[s(a-p(t)-Ex(t))]dt,
where x(t) is the population densities at time t, p(t) is the unit price of the population at time t. The constants r>0, k>0 are the intrinsic growth rate and the carrying capacity of the population, respectively. E>0 is constant representing harvesting effort. a is a positive constant parameter, representing the market capacity. s is the price speed adjustment, and it reflects market competition (i.e., perfect elasticity of price).In daily life, economic profit is a very important factor for governments, merchants, and even every citizen, so it is necessary to investigate biological economic systems, which are always described by differential-algebraic equations. At present, most of differential-algebraic equations can be found in the general power systems [11, 12], economic administration [13], robot system [14], mechanical engineering [15], and so on. However, to the best of our knowledge, the reports of such systems are few in biology. Considering the economic theory of fishery resource [16] proposed by Gordon in 1954, this paper studies a class of singular biological economic model as follows:
(2)dx(t)=[rx(t)(1-x(t)k)-Ex(t)]dt,dp(t)=[s(a-p(t)-Ex(t))]dt,0=Ep(t)x(t)-cEx(t)-m(t),
where c>0 is harvesting cost per unit harvesting effort for the population, m(t)≥0 is the economic profit per unit harvesting effort at time t. From the previous modeling view, the economic profit is variable. In order to facilitate the analysis, many authors assume that the economic profit is constant (see [17–20]). In fact, the economic profit is not always constant and may vary with numerous factors, such as seasonality, revenue, market demand, and harvesting cost and so on. Hence, it is more reasonable that the economic profit is a variable from real world point of view, and model (2) is more suitable for the analysis of the dynamic behavior of biological economic systems. The economic profit is affected by the unit price of the population, and the population densities of our discussed model and other factors are omitted. However, owing to difficulty arising in analysis, few results are concerned with such systems.It is well known that the white noise always exists and we cannot omit the influence of the white noise to many dynamical systems. In reality, due to continuous fluctuations in the environment (e.g., variation in intensity of sunlight, temperature, water level, etc.), parameters involved in models are not absolute constants, but they always fluctuate around some average value. As a result, the population density never attains a fixed value with the advancement of time but rather exhibits continuous oscillation around some average values. Based upon these factors, stochastic population models have received more and more attention [21–23]. Many authors studied the effect of the stochastic perturbation to the biological economic system with different functional responses, such as [24–27]. However, only a small amount of work has been studied for the stochastic perturbation biological economic systems with differential-algebraic equations. In this paper, we introduce stochastic perturbation to system (2) and obtain the following stochastic singular system:
(3)dx(t)=[rx(t)(1-x(t)k)-Ex(t)]dt+σx(t)dW(t),dp(t)=s(a-p(t)-Ex(t)),0=Ep(t)x(t)-cEx(t)-m(t),
where W(t) is one-dimensional Brownian motion with W(0)=0 and σ2 is intensities of the white noise.The proposed system (3) has a more widely practical background in this paper. For example, we show that there is a stationary distribution of system (3) if the white noise is small. While if the white noise is large, we prove that the prey population will either extinct or its distribution converges to a probability measure. So, we may study the economic profit through the changes of population density and the variation of market price in the stochastic environment. For this reason, we say the stochastic model is more realistic than the deterministic model.The aim of this paper is to discuss the finite-time stability and control problem of the system (2) and (3). It is now worth pointing out that the stability, stabilization, and control performance mentioned previously concern the desired behavior of the controlled dynamic systems over an infinite-time interval, which always deals with the asymptotic property of system trajectories [28–30]. But in some practical processes, the asymptotically stable system over an infinite-time interval does not mean that it has good transient characteristics, for instance, biochemistry reaction system, robot control system, communication network system, and so forth. Therefore, it is necessary to study the transient behavior over a finite-time interval.On the other hand, the literature on finite-time stability (FTS) (or short-time stability) of systems has attracted particular interests of researchers. Some earlier results on finite-time stability date back to the 1960s [31]. Recently, the concept of finite-time stability has been revisited in the light of linear matrix inequality theory. Several sufficient conditions for finite-time stability and stabilization of continuous time systems or discrete time systems have been presented [32, 33]. References [34, 35] have extended the definition of finite-time stability to linear systems with impulsive effects or singular systems with impulsive effects, respectively, and derived some sufficient conditions for finite-time stability and stabilization problems. Very recently, [36] generalized the results of [32] to linear stochastic systems. However, to date and to the best of our knowledge, no work so far has been done on the problem of finite-time stability for a class of stochastic singular biological economic systems. The problems are important in many practice applications, which motivate the main purpose of our study.T-S fuzzy model, proposed by Takagi and Sugeno [37] in 1985, has been shown to be a powerful tool for modeling complex nonlinear systems. It is wellknown that, by means of the T-S fuzzy model, a nonlinear system can be represented by a weighted sum of some simple linear subsystems and then can be stabilized by a model-based fuzzy control. Therefore, many issues related to the stability analysis of the complex nonlinear systems can be discussed using the T-S fuzzy model.In this paper, we investigate the problem of finite-time stability and control for stochastic singular biological economic systems. Our results are totally different from those previous results. Firstly, the system (2) and (3) are established for fuzzy models based on T-S fuzzy control approach. Then, sufficient conditions of finite-time stability are obtained for the system (2) and (3). One of the contributions of this paper is that the economic profit, which is a variable rather than a constant of the biological economic systems, is discussed. We study the economic profit fluctuations through the unit price range and the population changes. The other contribution is to study the finite-time stability for the biological economic system with stochastic perturbation, and sufficient criterions are presented for the solvability of the problem, which can be reduced to a feasibility problem involving restricted linear matrix, inequalities with a fixed parameter. Based on this, the state feedback controller is designed so that the biological population can be controlled in the bounded range by the management of the open resource. The corresponding sufficient conditions are given as well. Finally, some numerical simulations on feasible region are given to illustrate the effectiveness of our approach.The rest of the paper is organized as follows. In Section2, fuzzy modeling and preliminaries are stated in this paper. In Section 3, the finite-time stability of stochastic biological economic systems is studied, which provides some sufficient conditions for systems to be FTS. Meanwhile, theorems of sufficient conditions for FTS are established. In Section 4, finite-time controllers are designed for the stochastic biological economic systems. Numerical examples are provided in Section 5. Conclusions are given in Section 6.Notations. In this paper, the superscript “T” stands for matrix transposition. The notation X>0 means X is a positive definite matrix. ε(x) denotes the expectation of stochastic variable x. deg(·) denotes the degree of the determinant. In addition, the symbol “*” denotes the transposed elements in the symmetric positions of a matrix and diag(·) stands for a block-diagonal matrix. λmin(P) and λmax(P) denote the smallest and the largest eigenvalue of matrix P, respectively.


## 2. Fuzzy Modeling and Preliminaries

Considering the model (2), for simplicity, we use the following coordinate transformation
(4)ξ(t)=p(t)-a.
Then, the model (2) can be rewritten as follows:
(5)dx(t)=[rx(t)(1-x(t)k)-Ex(t)]dt,dξ(t)=[-sEx(t)-sξ(t)]dt,0=(a-c)Ex(t)+Ex(t)ξ(t)-m(t).
Let X(t)=[x(t)ξ(t)m(t)]Τ; the previous model can be expressed by fuzzy model, which consists of the following two rules:
(6)R1:Ifx(t)isF1,thenΞX˙(t)=A1X(t)R2:Ifx(t)isF2,thenΞX˙(t)=A2X(t),whereΞ=[100010000], A1=[2r-E00-sE-s0aE-cE-kE-1], A2=[-E00-sE-s0aE-cEkE-1].F
i
(
i
=
1,2
) is the fuzzy set, hi(x(t))≥0 is a membership function of Fi, and h1(x(t))=(1/2)(1-(x(t)/k)), h2(x(t))=(1/2)(1+(x(t)/k)), and ∑i=12hi(x(t))=1.By using a standard fuzzy singleton inference method, a singleton fuzzifier to produce a fuzzy inference, and a weighted center-average defuzzifier, we can obtain the global fuzzy system as follows:(7)ΞX˙(t)=AhX(t),
where Ah=∑i=12hi(x(t))Ai.As for the model (3), using the similar method, we can obtain
(8)ΞdX(t)=AhX(t)dt+CX(t)dW(t),
where C=[σ00000000].Definition 1.
(i) The pair(Ξ,Ah) is said to be regular if det(λΞ-Ah) is not identically zero.
(ii) The pair(Ξ,Ah) is said to be impulse free if deg(det(λΞ-Ah))=rank(Ξ).Remark 2.
In this paper, through the calculation, it is known that the pair(Ξ,Ah) is regular.Definition 3.
Given two positive real numbersc1, c2 with c1≤c2, a positive definite matrix R, system (7) is said to be finite-time stable with respect to (c1,c2,T,R) if
(9)XΤ(0)ΞΤRΞX(0)≤c1⟹XΤ(t)ΞΤRΞX(t)<c2,23222222222222222222222222∀t∈[0,T].Definition 4.
Given two positive real numbersc1, c2 with c1≤c2, a positive definite matrix R, system (8) is said to be stochastic finite-time stable with respect to (c1,c2,T,R) if
(10)ε{XΤ(0)ΞΤRΞX(0)}≤c1⟹ε{XΤ(t)ΞΤRΞX(t)}<c2,22222222222222222222222222222222∀t∈[0,T].Remark 5.
Asymptotic stability and FTS are independent concepts of each other as it was described in [31]. In certain cases, asymptotic stability can be also considered as an additional requirement while restricting our attention over a finite-time interval. The definition of FTS can be interpreted in terms of ellipsoidal domains. The set defined by {XΤ(0)ΞΤRΞX(0)≤c1} contains all the admissible initial states. Instead, the inequality {XΤ(t)ΞΤRΞX(t)<c2} defines a time-varying ellipsoid which bounds the state trajectory over the finite-time interval t∈[0,T]. In addition, we can fix the upper bound beforehand in finite-time control, and in general, the initial and ultimates bound are not the same. Thus, this concept will be more meaningful in practice, such as aircraft control, networked control, and robot operating.Definition 6 (see [38]).
LetV(x(t),t) be the stochastic Lyapunov function of the stochastic system, where x(t) satisfies the following stochastic differential equation:
(11)dx(t)=f(t)dt+g(t)dW(t).
Define a weak infinitesimal operator ℒ of random process {x(t),t>0}:
(12)ℒV(x(t),t)=Vt(x(t),t)+Vx(x(t),t)f(t)+12tr[gΤ(t)Vxx(x(t),t)g(t)].Lemma 7 (see [39], Gronwalls inequality).
Letv(t) be a nonnegative function such that
(13)v(t)≤a+b∫0tv(s)ds,0≤t≤T,
for some constants a,b≥0, then we have
(14)v(t)≤aexp(bt),0≤t≤T.Lemma 8 (Schur’s complement).
For real matricesN,MT=M, and R=RT>0, the following three conditions are equivalent:(1)
M+NR-1NT<0,(2)
(MNNT-R)<0,(3)
(M-N-NT-R)<0.Lemma 9 (see [40]).
The following items are true.(i)
Assume thatrank(Ξ)=r, there exist two orthogonal matrices U and V such that Ξ has the decomposition as
(15)Ξ=U[Σr0*0]VT=U[Ir0*0]vT,
where Σr=diag{δ1,δ2,…,δr} with δk>0 for all k=1,2,…,r. Partition U=[U1U2], V=[V1V2] and v=[V1ΣrV2] with ΞV2=0 and U2ΤΞ=0.(ii)
IfP satisfies
(16)ΞPΤ=PΞΤ≥0.
Then P~=UΤPv with U and v satisfying (8) if and only if
(17)P~=[P11P120P22]
with P11≥0∈ℝr×r. In addition, when P is nonsingular, we have P11>0 and det(P22)≠0. Furthermore, P satisfying (16) can be parameterized as
(18)P=Ξv-ΤYv-1+UZV2Τ,
where Y=diag{P11,Ψ}, Z=[P12ΤP22Τ]Τ, and Ψ∈ℝ(n-r)×(n-r) is an arbitrary parameter matrix.(iii)
If P is a nonsingular matrix,R and Ψ are two positive definite matrices, P and Ξ satisfy (16), Y is a diagonal matrix from (18), and the following equality holds:
(19)P-1Ξ=ΞΤR1/2QR1/2Ξ.
Then the positive definite matrix Q=R-1/2UY-1UΤR-1/2 is a solution of (19).


## 3. FTS of the Stochastic Singular Biological Economic System

This section provides finite-time stability analysis results for the stochastic singular biological economic system, some novel sufficient conditions for finite-time stability of the stochastic singular system are given. Through the previous analysis, we know that the system (2) and (3) are equivalent to the system (7) and (8), respectively. Because the system (7) and (8) are regular, the following theorems are given.


### 3.1. FTS of the System (7)

Theorem 10.
For a given time constantT>0 and a scalar α>0, the system (7) is finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P and a symmetric positive definite matrix Q such that the following conditions hold:
(20)PΞΤ=ΞPΤ≥0,(21a)
A
1
P
Τ
+
P
A
1
Τ
-
α
Ξ
P
Τ
<
0
,
(21b)
A
2
P
Τ
+
P
A
2
Τ
-
α
Ξ
P
Τ
<
0
,(22)
P
-
1
Ξ
=
Ξ
Τ
R
1
/
2
Q
R
1
/
2
Ξ
,
(23)
λ
max
(
Q
)
c
1
e
α
T
-
c
2
λ
min
(
Q
)
<
0
.Proof.
Consider the following Lyapunov function:(24)V(X(t),t)=XΤ(t)P-1ΞX(t).
We first prove that the following condition (25) is equivalent to (21a) and (21b):
(25)V˙(X(t),t)<αV(X(t),t).
When t∈[0,T], we calculate the time derivative of V(X(t),t) along the solution of system (7) and obtain
(26)V˙(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah)X(t)
which leads to
(27)V˙(X(t),t)-αV(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah-αP-1Ξ)X(t).
So, (27) is equivalent to the following inequality
(28)AhΤP-Τ+P-1Ah-αP-1Ξ<0.
Pre- and postmultiplying (28) by P and PΤ, we get
(29)AhPΤ+PAhΤ-αΞPΤ<0.
When Ah=∑i=12hi(x(t))Ai, through the matrix decomposition, we can see that (29) is equivalent to (21a) and (21b).
Thus, dividing byV(X(t),t) on both sides of (25) and integrating it from 0 to t with t∈[0,T], one gets
(30)V(X(t),t)<V(X(0),0)eαt.
By condition (22), one can deduce that
(31)V(X(t),t)=XΤ(t)P-1ΞX(t)=XΤ(t)ΞΤR1/2QR1/2ΞX(t)≥λmin(Q)XΤ(t)ΞΤRΞX(t),V(x(0),0)eαt=XΤ(0)ΞΤR1/2QR1/2ΞX(0)eαt≤λmax(Q)XΤ(0)ΞΤRΞX(0)eαt≤λmax(Q)c1eαT.
From (31), it is easy to obtain
(32)XΤ(t)ΞΤRΞX(t)≤λmax(Q)λmin(Q)c1eαT.
Considering condition (23) and inequality (32), it is easy to get XΤ(t)ΞΤRΞX(t)<c2 for all t∈[0,T]. This completes the proof of the theorem.


### 3.2. FTS of the System (8)

Theorem 11.
For a given time constantT>0 and a scalar α>0, the system (8) is stochastic finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P and a symmetric positive definite matrix Q such that (20), (22), and (23) hold and the following matrix inequalities are satisfied:(33a)
[
A
1
P
Τ
+
P
A
1
Τ
-
α
Ξ
P
Τ
(
Ξ
C
P
Τ
)
Τ
Ξ
C
P
Τ
-
Q
~
]
<
0
,
(33b)
[
A
2
P
Τ
+
P
A
2
Τ
-
α
Ξ
P
Τ
(
Ξ
C
P
Τ
)
Τ
Ξ
C
P
Τ
-
Q
~
]
<
0
,whereQ~=R-1/2Q-1R-1/2.Proof.
Let us consider the quadratic Lyapunov function(34)V(X(t),t)=XΤ(t)P-1ΞX(t).
Let ℒ be the infinitesimal generator (also called the averaged derivative). We first show the following condition:
(35)ℒV(X(t),t)<αV(X(t),t).
Now we will prove that condition (35) is equivalent to (33a) and (33b). Applying Ito’s formula, it follows that
(36)ℒV(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah+CΤP-1ΞC)X(t)
which leads to
(37)ℒV(X(t),t)-αV(X(t),t)=XΤ(t)(AhΤP-Τ2222222222+P-1Ah+CΤP-1ΞC-αP-1Ξ)X(t).
Combining condition (22) with Q~=R-1/2Q-1R-1/2, we can obtain
(38)ℒV(X(t),t)-αV(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah222222222+CΤΞΤQ~-1ΞC-αP-1Ξ)X(t).
So, (37) is equivalent to
(39)AhΤP-Τ+P-1Ah+CΤΞΤQ~-1ΞC-αP-1Ξ<0.
Pre- and postmultiplying (39) by P and PΤ, we get
(40)PAhΤ+AhPΤ+PCΤΞΤQ~-1ΞCPΤ-αΞPΤ<0.
When Ah=∑i=12hi(x(t))Ai, through the matrix decomposition and by Lemma 8, we can see that (40) is equivalent to (33a) and (33b).
Hence, integrating both sides of (35) from 0 to t with t∈[0,T] and then taking the expectation, it yields
(41)ε{V(X(t),t)}<V(X(0),0)+α∫0tε{V(X(s),s)}ds.
By Lemma 7, we can obtain
(42)ε{V(X(t),t)}<V(X(0),0)eαt.
From (34), it follows that
(43)ε{V(X(t),t)}=ε{XΤ(t)ΞΤR1/2QR1/2ΞX(t)}≥λmin(Q)ε{XΤ(t)ΞΤRΞX(t)},V(X(0),0)eαt=XΤ(0)ΞΤR1/2QR1/2ΞX(0)eαt≤λmax(Q)XΤ(0)ΞΤRΞX(0)eαt≤λmax(Q)c1eαT.
From (43), it is easy to obtain
(44)ε{XΤ(t)ΞΤRΞX(t)}<λmax(Q)λmin(Q)c1eαT.
Considering condition (23) and inequality (44), it is easy to get ε{XΤ(t)ΞΤRΞX(t)}<c2 for all t∈[0,T]. This completes the proof of the theorem.


## 3.1. FTS of the System (7)

Theorem 10.
For a given time constantT>0 and a scalar α>0, the system (7) is finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P and a symmetric positive definite matrix Q such that the following conditions hold:
(20)PΞΤ=ΞPΤ≥0,(21a)
A
1
P
Τ
+
P
A
1
Τ
-
α
Ξ
P
Τ
<
0
,
(21b)
A
2
P
Τ
+
P
A
2
Τ
-
α
Ξ
P
Τ
<
0
,(22)
P
-
1
Ξ
=
Ξ
Τ
R
1
/
2
Q
R
1
/
2
Ξ
,
(23)
λ
max
(
Q
)
c
1
e
α
T
-
c
2
λ
min
(
Q
)
<
0
.Proof.
Consider the following Lyapunov function:(24)V(X(t),t)=XΤ(t)P-1ΞX(t).
We first prove that the following condition (25) is equivalent to (21a) and (21b):
(25)V˙(X(t),t)<αV(X(t),t).
When t∈[0,T], we calculate the time derivative of V(X(t),t) along the solution of system (7) and obtain
(26)V˙(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah)X(t)
which leads to
(27)V˙(X(t),t)-αV(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah-αP-1Ξ)X(t).
So, (27) is equivalent to the following inequality
(28)AhΤP-Τ+P-1Ah-αP-1Ξ<0.
Pre- and postmultiplying (28) by P and PΤ, we get
(29)AhPΤ+PAhΤ-αΞPΤ<0.
When Ah=∑i=12hi(x(t))Ai, through the matrix decomposition, we can see that (29) is equivalent to (21a) and (21b).
Thus, dividing byV(X(t),t) on both sides of (25) and integrating it from 0 to t with t∈[0,T], one gets
(30)V(X(t),t)<V(X(0),0)eαt.
By condition (22), one can deduce that
(31)V(X(t),t)=XΤ(t)P-1ΞX(t)=XΤ(t)ΞΤR1/2QR1/2ΞX(t)≥λmin(Q)XΤ(t)ΞΤRΞX(t),V(x(0),0)eαt=XΤ(0)ΞΤR1/2QR1/2ΞX(0)eαt≤λmax(Q)XΤ(0)ΞΤRΞX(0)eαt≤λmax(Q)c1eαT.
From (31), it is easy to obtain
(32)XΤ(t)ΞΤRΞX(t)≤λmax(Q)λmin(Q)c1eαT.
Considering condition (23) and inequality (32), it is easy to get XΤ(t)ΞΤRΞX(t)<c2 for all t∈[0,T]. This completes the proof of the theorem.


## 3.2. FTS of the System (8)

Theorem 11.
For a given time constantT>0 and a scalar α>0, the system (8) is stochastic finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P and a symmetric positive definite matrix Q such that (20), (22), and (23) hold and the following matrix inequalities are satisfied:(33a)
[
A
1
P
Τ
+
P
A
1
Τ
-
α
Ξ
P
Τ
(
Ξ
C
P
Τ
)
Τ
Ξ
C
P
Τ
-
Q
~
]
<
0
,
(33b)
[
A
2
P
Τ
+
P
A
2
Τ
-
α
Ξ
P
Τ
(
Ξ
C
P
Τ
)
Τ
Ξ
C
P
Τ
-
Q
~
]
<
0
,whereQ~=R-1/2Q-1R-1/2.Proof.
Let us consider the quadratic Lyapunov function(34)V(X(t),t)=XΤ(t)P-1ΞX(t).
Let ℒ be the infinitesimal generator (also called the averaged derivative). We first show the following condition:
(35)ℒV(X(t),t)<αV(X(t),t).
Now we will prove that condition (35) is equivalent to (33a) and (33b). Applying Ito’s formula, it follows that
(36)ℒV(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah+CΤP-1ΞC)X(t)
which leads to
(37)ℒV(X(t),t)-αV(X(t),t)=XΤ(t)(AhΤP-Τ2222222222+P-1Ah+CΤP-1ΞC-αP-1Ξ)X(t).
Combining condition (22) with Q~=R-1/2Q-1R-1/2, we can obtain
(38)ℒV(X(t),t)-αV(X(t),t)=XΤ(t)(AhΤP-Τ+P-1Ah222222222+CΤΞΤQ~-1ΞC-αP-1Ξ)X(t).
So, (37) is equivalent to
(39)AhΤP-Τ+P-1Ah+CΤΞΤQ~-1ΞC-αP-1Ξ<0.
Pre- and postmultiplying (39) by P and PΤ, we get
(40)PAhΤ+AhPΤ+PCΤΞΤQ~-1ΞCPΤ-αΞPΤ<0.
When Ah=∑i=12hi(x(t))Ai, through the matrix decomposition and by Lemma 8, we can see that (40) is equivalent to (33a) and (33b).
Hence, integrating both sides of (35) from 0 to t with t∈[0,T] and then taking the expectation, it yields
(41)ε{V(X(t),t)}<V(X(0),0)+α∫0tε{V(X(s),s)}ds.
By Lemma 7, we can obtain
(42)ε{V(X(t),t)}<V(X(0),0)eαt.
From (34), it follows that
(43)ε{V(X(t),t)}=ε{XΤ(t)ΞΤR1/2QR1/2ΞX(t)}≥λmin(Q)ε{XΤ(t)ΞΤRΞX(t)},V(X(0),0)eαt=XΤ(0)ΞΤR1/2QR1/2ΞX(0)eαt≤λmax(Q)XΤ(0)ΞΤRΞX(0)eαt≤λmax(Q)c1eαT.
From (43), it is easy to obtain
(44)ε{XΤ(t)ΞΤRΞX(t)}<λmax(Q)λmin(Q)c1eαT.
Considering condition (23) and inequality (44), it is easy to get ε{XΤ(t)ΞΤRΞX(t)}<c2 for all t∈[0,T]. This completes the proof of the theorem.


## 4. Finite-Time Control of the Stochastic Singular Biological Economic System

Unstable fluctuations have always been regarded as unfavorable ones from the ecological managers’ point of view. Meanwhile, price variation depends on the gap between the demand function and supply function. In order to plan harvesting strategies and maintain the sustainable development of system, it is necessary to take action to stabilize biological population and market directed economy. This paper proposes the following state feedback control method for the system (2) and (3):
(45)dx(t)=[rx(t)(1-x(t)k)-Ex(t)+u(t)]dt,dp(t)=s(a-p(t)-Ex(t))dt,0=Ep(t)x(t)-cEx(t)-m(t),(46)dx(t)=[rx(t)(1-x(t)k)-Ex(t)+u(t)]dt+σx(t)dW(t),dp(t)=s(a-p(t)-Ex(t))dt,0=Ep(t)x(t)-cEx(t)-m(t),
where u(t) is control variable, representing the management of the open resource.Then, we use the similar T-S fuzzy method shown in Section2 and consider the following state feedback fuzzy controller:
(47)u(t)=∑i=12hi(x(t))LiX(t),
where Li is the state feedback gain to be designed.The model (45) and (46) can be expressed by the following fuzzy models:
(48)ΞX˙(t)=∑i=12hi(x(t))∑j=12hj(x(t))(Ai+BLj)X(t),(49)ΞdX(t)=∑i=12hi(x(t))∑j=12hj(x(t))×[(Ai+BLj)X(t)dt+CX(t)dW(t)],
where B=[100]Τ.In the sequence, we provide finite-time stabilization results for above the mentioned closed-loop system (48) and (49). At the same time, some novel sufficient conditions for finite-time stabilization of the closed-loop singular biological economic system are given.


### 4.1. Finite-Time Control of the System (45)

Theorem 12.
For a given time constantT>0 and a scalar α>0, there exists a state feedback controller (47) such that the closed-loop system (48) is finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P, a symmetric positive definite matrix Q, and matrix Yj, j=1,2, such that (22), (23) hold, and the following matrix inequalities are satisfied:
(50)PΞΤ=ΞPΤ≥0,(51a)Φii<0,i=j,i,j=1,2,(51b)Φij+Φji<0,i<j,i,j=1,2,
whereΦij=PAiΤ+YjΤBΤ+AiPΤ+BYj-αΞPΤ. Then, the desired state feedback gain can be chosen as Lj=YjP-Τ.Proof.
Considering Theorem10 and the closed-loop system (48), we can obtain the following matrix inequality
(52)∑i=12∑j=12hi(x(t))hj(x(t))Πij<0,
where Πij=(Ai+BLj)PΤ+P(Ai+BLj)Τ-αΞPΤ.
LetYj=LjPΤ. Then,
(53)∑i=12∑j=12hi(x(t))hj(x(t))×(PAiΤ+YjΤBT+AiPΤ+BYj-αΞPΤ)<0.
Obviously, (53) is equivalent to
(54)∑i=12hi2(x(t))Φii+∑i=12∑j=12hi(x(t))hj(x(t))(Φij+Φji)<0.
The above inequality (54) is equivalent to (51a) and (51b).
Then, by a similar argument in the proof of Theorem10 and using (22)-(23), we can easily get the above results. This completes the proof of the theorem.


### 4.2. Finite-Time Control of the System (46)

Theorem 13.
For a given time constantT>0 and a scalar α>0, there exists a state feedback controller (47) such that the closed-loop system (49) is stochastic finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P, a symmetric positive definite matrix Q, and matrix Yj, j=1,2, such that (50), (22), and (23) hold, and the following matrix inequalities are satisfied:(55a)[Ψii(ΞCPΤ)ΤΞCPΤ-Q~]<0i=j,i,j=1,2,(55b)[Ψij+Ψji(ΞCPΤ)ΤΞCPΤ-Q~]<0i<j,i,j=1,2,
whereΨij=PAiΤ+YjΤBΤ+AiPΤ+BYj-αΞPΤ and Q~=R-1/2Q-1R-1/2. Then, the desired state feedback gain can be chosen as Lj=YjP-Τ.Proof.
Considering Theorem11 and the closed-loop system (49), this proof follows a similar way with the proof in Theorem 12, and it is easy to obtain the previous results. This completes the proof of the theorem.Remark 14.
The state feedback control method can be implemented by the management of the open resource so that the biological population can be controlled at the bounded range, and those unfavorable phenomena will be eliminated. In our daily life, managers can take action, such as adjusting revenue, drawing out favorable policy to encourage fishery, and abating pollution, in order to control the population density and economic profit.


## 4.1. Finite-Time Control of the System (45)

Theorem 12.
For a given time constantT>0 and a scalar α>0, there exists a state feedback controller (47) such that the closed-loop system (48) is finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P, a symmetric positive definite matrix Q, and matrix Yj, j=1,2, such that (22), (23) hold, and the following matrix inequalities are satisfied:
(50)PΞΤ=ΞPΤ≥0,(51a)Φii<0,i=j,i,j=1,2,(51b)Φij+Φji<0,i<j,i,j=1,2,
whereΦij=PAiΤ+YjΤBΤ+AiPΤ+BYj-αΞPΤ. Then, the desired state feedback gain can be chosen as Lj=YjP-Τ.Proof.
Considering Theorem10 and the closed-loop system (48), we can obtain the following matrix inequality
(52)∑i=12∑j=12hi(x(t))hj(x(t))Πij<0,
where Πij=(Ai+BLj)PΤ+P(Ai+BLj)Τ-αΞPΤ.
LetYj=LjPΤ. Then,
(53)∑i=12∑j=12hi(x(t))hj(x(t))×(PAiΤ+YjΤBT+AiPΤ+BYj-αΞPΤ)<0.
Obviously, (53) is equivalent to
(54)∑i=12hi2(x(t))Φii+∑i=12∑j=12hi(x(t))hj(x(t))(Φij+Φji)<0.
The above inequality (54) is equivalent to (51a) and (51b).
Then, by a similar argument in the proof of Theorem10 and using (22)-(23), we can easily get the above results. This completes the proof of the theorem.


## 4.2. Finite-Time Control of the System (46)

Theorem 13.
For a given time constantT>0 and a scalar α>0, there exists a state feedback controller (47) such that the closed-loop system (49) is stochastic finite-time stable with respect to (c1,c2,T,R) if there exist a nonsingular matrix P, a symmetric positive definite matrix Q, and matrix Yj, j=1,2, such that (50), (22), and (23) hold, and the following matrix inequalities are satisfied:(55a)[Ψii(ΞCPΤ)ΤΞCPΤ-Q~]<0i=j,i,j=1,2,(55b)[Ψij+Ψji(ΞCPΤ)ΤΞCPΤ-Q~]<0i<j,i,j=1,2,
whereΨij=PAiΤ+YjΤBΤ+AiPΤ+BYj-αΞPΤ and Q~=R-1/2Q-1R-1/2. Then, the desired state feedback gain can be chosen as Lj=YjP-Τ.Proof.
Considering Theorem11 and the closed-loop system (49), this proof follows a similar way with the proof in Theorem 12, and it is easy to obtain the previous results. This completes the proof of the theorem.Remark 14.
The state feedback control method can be implemented by the management of the open resource so that the biological population can be controlled at the bounded range, and those unfavorable phenomena will be eliminated. In our daily life, managers can take action, such as adjusting revenue, drawing out favorable policy to encourage fishery, and abating pollution, in order to control the population density and economic profit.


## 5. Numerical Examples

To illustrate the results obtained, let us consider the following particular cases. We choose the ecological parameters according to biological sardine data of the Moroccan western coastline area, and these parameters are given in [6].Example 1.
Consider finite-time stabilization problem of the biological economic model, and select the ecological parameters asr=0.015, k=1702.8, E=0.008, s=0.005, and c=200, a=800 in appropriate units. Then, the closed-loop system (48) can be described with the following data:
(56)Ξ=[100010000],A1=[0.02200.00000.0000-0.0000-0.00500.00004.8000-13.6224-1.0000],A2=[-0.00800.00000.0000-0.0000-0.00500.00004.800013.6224-1.0000],B=[100]Τ.
Let α=0.001, c1=10000, c2=10000000, T=1000, R=I. By Theorem 12, we can calculate that
(57)P=105×[0.01450.00030.04630.00030.00020.0016001.8983],Q=[0.0007-0.00140-0.00140.05800000.0294],L1=[-36.699169.2352-0.0098],L2=[-36.653869.2506-0.0147].
So, we claim that XΤ(t)ΞΤRΞX(t)<10000000, ∀t∈[0,1000]. The trajectory is shown in Figure 2.
Figure1 shows the state trajectory of the open-loop biological economic system when t∈[0,1000]. Figure 2 shows the state trajectory of the closed-loop biological economic system with state feedback controller (47) when t∈[0,1000].Figure 1
Trajectory of the open-loop biological economic system.Figure 2
Trajectory of the closed-loop biological economic system.Example 2.
Consider finite-time stabilization problem of the biological economic model with the white noise in Example1. σ2 is intensities of the white noise, and let σ2=0.01. Select the ecological parameters as r=0.015, k=1702.8, E=0.008, s=0.005, c=200, and a=800 in appropriate units. The closed-loop system (49) can be formulated with the following parameters:
(58)Ξ=[100010000],A1=[0.02200.00000.0000-0.0000-0.00500.00004.8000-13.6224-1.0000],B=[100]Τ,A2=[-0.00800.00000.0000-0.0000-0.00500.00004.800013.6224-1.0000],C=[0.100000000].
Let α=0.001, c1=10000, c2=10000000, T=1000, and R=I. By Theorem 13, we can calculate that
(59)P=105×[0.08370.00200.26790.00200.00040.0098003.7772],Q=[0.0001-0.00070-0.00070.03020000.0049],L1=[-8.657947.2591-0.0345],L2=[-8.622147.2785-0.0364].
So, we claim that ε{XΤ(t)ΞΤRΞX(t)}<10000000, for all t∈[0,1000]. The trajectory is shown in Figure 4.
Figure3 shows the state trajectory of the open-loop stochastic singular biological economic system with the white noise σ2. Taking σ2=0.01, we know that the population density and the unit price are not stable over the finite-time interval [0,1000]. The economic profit is random fluctuation in Figure 3, because the economic profit is affected by the unit price of the population and the population density of our proposed biological economic model in the random disturbance environment.
Figure4 shows the state trajectory of the closed-loop stochastic singular biological economic system with the white noise σ2=0.01 through the state feedback controller (47). From Figure 4, it can be seen that the economic profit tends to be stable over the finite-time interval [0,1000].Figure 3
Trajectory of the open-loop stochastic singular system.Figure 4
Trajectory of the closed-loop stochastic singular system.


## 6. Conclusions

In this paper, we study the problems of finite-time stability and control for a class of stochastic singular biological economic systems based on T-S fuzzy model. By the theorems we presented, some novel sufficient conditions are derived which ensure that the systems are stable in finite time. And then, the corresponding controller design methods are obtained. Finally, numerical examples are presented to illustrate the effectiveness of the proposed results. Our proposed results in this paper can also be applied to the study of other classes of systems.



---

*Source: 946491-2013-06-10.xml*