# A Mathematical Model for the Dynamics of a Fish Algae Consumption Model with Impulsive Control Strategy

Authors: Jin Yang; Min Zhao
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

A dynamic mathematical model of fish algae consumption with an impulsive control strategy is proposed and analyzed in detail. It is shown that the system has a globally asymptotically stable algae-eradication periodic solution which can be obtained using the Floquet theory of impulsive differential equations and small-amplitude perturbation techniques. The conditions for the permanence of the system can also be determined. Numerical results for impulsive perturbations show the rich dynamic behavior of the system. All these results may be useful in controlling eutrophication.


---

## Body

## 1. Introduction

One of the most common ecological and environmental problems of inland water bodies is eutrophication, which diminishes water quality by spurring excessive growth of algae and increasing amounts of suspended organic material [1–3]. Because of climate change and discharges of industrial wastewater into rivers, the self-purification capacity of the water body decreases. This provides favorable conditions for algae growth, leading to repeated outbreaks of abnormal proliferations of algae blooms in tributaries; recently, this phenomenon has been increasing [4]. Therefore, research on how to control the amount of algae in the environment is of great theoretical importance and practical significance.There are many ways to control algae populations. The most widely used method is chemical control; this is not discussed in detail here because it has many negative impacts. More recently, another important method, called biological control, has been introduced. Biological control is the practice of using natural enemies (fish) to suppress algae, as has already been done for pest control [5–8]. These natural enemies can prevent algae populations from increasing to levels that cause eutrophication.Mathematical models of ecological population dynamics have not only to account for growth and interactions, but also to improve the understanding of the functioning of food chains and webs and their dependence on some environmental conditions [9–12]. Impulsive differential equations have been widely used to study the mathematical properties of impulsive predator-prey models or three-species food-chain models [13–17]. Moreover, the theory of impulsive differential equations is being recognized, not only to be richer than the corresponding theory of differential equations without impulses, but also to represent a more natural framework for the mathematical modeling of real-world phenomena [18–20].Recently, because of eutrophication, the nuisance algal (cyanobacteria and green algae) blooms frequently come forth in the Zeya reservoir which is located in Wenzhou. The algal blooms can cause clogging and blocking of the filtration system and result in millions of people without drinking water. In order to apply the biological principle to control algal blooms in the Zeya reservoir, we carry out good biological control agents for the algal blooms. Firstly, we bring up two species of filter-feeding fish, silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis), which are thought to be suited to control algal biomass directly in freshwater reservoir, especially cyanobacteria and green algae. Secondly, the decision maker must give cost-benefit and maintenance considerations; the fish should be harvested at an appropriate time. Finally, in order to prevent the rapid growth of the cyanobacteria and green algae, we must release a certain amount of filter-feeding fish at a fixed time. These processes are subject to short-term perturbations. Consequently, it is obvious to assume that these perturbations act instantaneously in the form of impulses. Hence, a fish algae consumption model can be described by the following differential equations:(1.1)dx(t)dt=r1x(t)(G0-x(t)G1-x(t))-c1x2(t)-u1x(t)z(t)b1+x(t),dy(t)dt=r2y(t)(1-(y(t)/ym)1-(y(t)/yn))-c2y2(t)-u1ay(t)z(t)x(t)+ay(t)+b2,t≠nT,dz(t)dt=-mz(t)+u2ay(t)z(t)x(t)+ay(t)+b2+u3x(t)z(t)b1+x(t),Δx(t)=0,Δy(t)=0,t=nT,Δz(t)=-δz(t)+p,
where x(t), y(t), z(t) are the densities of the cyanobacteria population, green algae population, and fish population at time t, Δx(t)=x(t+)-x(t);Δy(t)=y(t+)-y(t), and Δz(t)=z(t+)-z(t); r1,r2 is a growth parameter which is related to the biological characteristics of populations and the rationalization of environmental resources; r1G0  (0≤G0≤G1) is the carrying capacity of the cyanobacteria population x; G1 is the limiting value of available resources; ym(0≤(ym/yn)≤1) is the maximum density of green algae population y (i.e., the environmental carrying capacity); yn is a nutritional parameter which is related to the resource conditions of the environment; u1 is the cropping rate of fish; u2 and u3 are the average growth rate of fish that depends on the cyanobacteria and green algae for food; a is the capture rate (as a proportion); m is the average mortality rate for fish; b1 and b2 are the half-saturation constants; c1and c2 are the intraspecific competition rates of the cyanobacteria and green algae; T is the period of the impulsive effect; N is the set of all nonnegative integers; δ(δ∈[0,1]) is the proportion of fish harvested at fixed times t=nT; and p>0 is the number of fish released at times t=nT.The rest of this paper is arranged as follows: in the next section, some useful notations and basic properties of a system with impulsive effect are given, and then the main theorems and their proofs are developed using the Floquet theory. In Section3, numerical results will be presented which exhibit the rich dynamic behaviors of system (1.1), and these results will be briefly discussed. Conclusions are given in the last section.


## 2. Preliminaries and Mathematical Analysis

LetR+=[0,∞),R3+={X∈R3∣X>0}. Denote by f=(f1,f2,f3) the mapping defined by the right-hand sides of the first, second, and third equations of system (1.1). Let V:R+×R+3→R+; then Vis said to belong to class V0 if(1)
V is continuous in (nT,(n+1)T]×R+3, and for each X∈R+3,n∈N, lim(t,y)→(nT+,X)V(t,y)=V(nT+,X) exists.(2)
V is locally Lipschitzian in X.Definition 2.1 (see [18]).
LetV∈V0; then for (t,x)∈(nT,(n+1)T]×R+3, the upper right derivative of V(t) with respect to the impulsive differential system (1.1) is defined as
(2.1)D+V(t,X)=lim⁡h→0+sup⁡1h[V(t+h,X+hf(t,X))-V(t,X)].Remark 2.2.
(1) The solution of system (1.1) is a piecewise continuous function with X:R+→R+3; X(t) is continuous on (nT,(n+1)T], and X(nT+)=lim(t→nT+)X(nT) exists.
(2) The smoothness properties off guarantee the global existence and uniqueness of solutions of system (1.1) (for details, see [18, 19]).
Becausex′(t)=y′(t)=z′(t)=0 whenever x(t)=y(t)=z(t)=0   (t≠nT) and z(nT+)=(1-δ)z(nT)+p, the following lemma can be stated.Lemma 2.3.
AssumeX(t) is a solution of system (1.1) such that(1)
ifX(0+)≥0, then X(t)≥0  for all t≥0,(2)
ifX(0+)>0, then X(t)>0  for all t>0.
Next, an important comparison theorem on impulsive differential equations [18] will be stated.Lemma 2.4 (comparison theorem).
Assume thatV∈V0 and that
(2.2)D+V(t,X)≤g(t,V(t,X)),t≠nT,V(t,X(t+))≤φn(V(t,X)),t=nT,
where g:R+×R+→R is continuous on (nT,(n+1)T]×R+3 and for u∈R+, n∈N, lim(t,y)→(nT+,u)g(t,v)=g(nT+,u) exists, with φn:R+→R+ nondecreasing. Let r(t) be the maximal solution of the scalar impulsive differential equation:
(2.3)du(t)dt=g(t,u(t)),t≠nT,u(t+)=φn(u(t)),t=nT,u(0+)=u0,
existing on [0,∞). Then V(0+,X0)≤u0, which implies that V(t,X(t))≤r(t), t≥0, where X(t) is any solution of system (1.1).
If the cyanobacteria and green algae are eradicated, then system (1.1) will reduce to the following system:
(2.4)dz(t)dt=-mz(t),t≠nT,z(t+)=(1-δ)z(t)+p,t=nT,z(0+)=z0.Therefore, the following result can be easily obtained.Lemma 2.5.
z*(t)=pexp(-m(t-nT))/(1-(1-δ)exp(-mT)) is a positive periodic solution of (2.4) with initial value z*(0+)=p/(1-(1-δ)exp(-mT)),t∈(nT,(n+1)T],n∈N.
z(t)=(z(0+)-p/(1-(1-δ)exp(-mT)))exp(-mt)+z*(t) is a general solution of (2.4) with initial value z0=z(0+)≥0, where t∈(nT,(n+1)T], n∈N.
For a positive periodic solutionz*(t) of (2.4) and every solution z(t) of (2.4) with z0≥0, |z(t)-z*(t)|→0, t→∞.From Lemma2.5, the general solution z(t) converges to the positive periodic solution z*(t) of (2.4) when t→∞, and the complete expression for the cyanobacteria and green algae-eradication periodic solution (0,0,z*(t)) of system (1.1) becomes(2.5)(0,0,z*(t))=(0,0,pexp⁡(-m(t-nT))1-(1-δ)exp⁡(-mT)),
for t∈(nT,(n+1)T].Now, the Floquet theorem for impulsive equations [20] will be introduced:(2.6)dx(t)dt=A(t)x(t),t≠tk,x(t+)=x(t)+Bkx(t),t=tk,
with conditions as follows(1)
A(·)∈PC(R,Cn×n) and A(t+T)=A(t)(t∈R), where PC(R,Cn×n) is the set of all piecewise continuous matrix functions and Cn×n is the set of n×n matrices.(2)
Bk∈(Cn×n) and det(E+Bk)≠0, tk≤Tk+1.(3)
There exists aq∈N such that Bk+q=Bk, tk+q=tk+T.LetH(t) be a fundamental matrix of system (2.6), then there exists a unique nonsingular matrix M∈Cn×n such that H(t+T)=H(t)M. M is similar to the matrix of system (2.6) and has the same eigenvalues u1,u2,…,un. These eigenvalues u1,u2,…,un of system (2.6) are called Floquet multipliers.Lemma 2.6 (Floquet theorem [20]).
Assume that (1)–(3) hold, then system (2.6) is(1)
stable if and only ifui≤1(i=1,…,n), and to those ui=1, there correspond simple elementary divisors;(2)
asymptotically stable if and only if allui<1(i=1,…,n);(3)
unstable ifui>1 for some i=1,…,n.After these preliminaries, the main theorems of this paper will now be presented, then it will be proven that each solution of system (1.1) is ultimately bounded.Theorem 2.7.
There exists a positive constantM such that x(t)≤M, y(t)≤M, and z(t)≤M for each solution of system (1.1) with all t large enough.Proof.
DefineV(t,x(t)) such that
(2.7)V(t,x(t))=x(t)+u2y(t)+u1z(t).
Because(dx(t)/dt)≤r1x(t)-c1x2(t), then x(t)≤(r1/c1), and the upper right derivative of V(t,x(t)) can be calculated along a solution of system (1.1) when t≠nT:
(2.8)D+(V(t))+LV(t)≤(L+r1)x(t)-c1x2(t)+(Lu2+u2r2)y(t)-c2u2y2(t)+(Lu1-u1m+u1r1u3b1c1)z(t).
Let0<L<m-(r1u3/b1c1), then (L+r1)x(t)-c1x2(t), (Lu2+u2r2)y(t)-c2u2y2(t) are both bounded, so D+(V(t))+LV(t) is bounded, and L1 and L2 can be chosen such that
(2.9)dD+(V(t))dt≤-L1V(t)+L2,t≠nT,V(t+)≤V(t)+u1p,t=nT,
where L1 and L2 are positive constants. According to Lemma 2.4,
(2.10)V(t)≤(V(0+)-L2L1)exp⁡(-L1t)+u1p(1-exp⁡(-L1nT))1-exp⁡(-L1T)exp⁡(-L1(t-nT))+L2L1,
and hence
(2.11)V(t)≤L2L1+u1pexp⁡(L1T)exp⁡(L1T)-1,
as t→∞, so V(t,x(t)) is ultimately bounded. Therefore, each positive solution of system (1.1) is ultimately bounded, and the proof is completed.Theorem 2.8.
If(2.12)r1G0G1T<u1p(1-exp⁡(-mT))mb1(1-(1-δ)exp⁡(-mT)),r2T<u1ap(1-exp⁡(-mT))mb2(1-(1-δ)exp⁡(-mT)),
then the periodic solution (0,0,z*(t)) is locally asymptotically stable.Proof.
The local stability of the periodic solution(0,0,z*(t)) may be determined by considering the behavior of small-amplitude perturbations of the solution.
Definex(t)=u(t),y(t)=v(t),z(t)=w(t)+z*(t).
Then the linearization of system (1.1) becomes
(2.13)du(t)dt=(r1G0G1-u1z*(t)b1)u(t),dv(t)dt=(r2-u1az*(t)b2)v(t),t≠nT,dw(t)dt=u3z*(t)b1u(t)+u2az*(t)b2v(t)-mw(t),Δu(t)=0,Δv(t)=-δv(t),t=nT,Δw(t)=0,
and, as a result,
(2.14)(u(t)v(t)w(t))=Φ(t)(u(0)v(0)w(0)),0≤t≤T,
where Φ(t) satisfies
(2.15)dΦ(t)dt=(r1G0G1-u1z*(t)b1000r2-u1az*(t)b20u3z*(t)b1u2az*(t)b2-m)Φ(t),Φ(0)=I, the identity matrix, and
(2.16)(u(nT+)v(nT+)w(nT+))=(100010001-δ)(u(nT)v(nT)w(nT)).
Therefore, the stability of the periodic solution(0,0,z*(t)) is determined by the eigenvalues
(2.17)θ=(100010001-δ)Φ(T).
Therefore, all eigenvalues ofθ are given by
(2.18)λ1=exp⁡(∫0T(r1G0G1-u1z*(t)b1)dt),λ2=exp⁡(∫0T(r2-u1az*(t)b2)dt),λ3=(1-δ)exp⁡(-mT)<1.
According to Lemma2.6, (0,0,z*(t)) is locally asymptotically stable if  λ1<1,λ2<1, that is to say
(2.19)r1G0G1T<u1p(1-exp⁡(-mT))mb1(1-(1-δ)exp⁡(-mT)),r2T<u1ap(1-exp⁡(-mT))mb2(1-(1-δ)exp⁡(-mT)),
which completes the proof.Theorem 2.9.
If(2.20)r1G0G1T<u1p(1-exp⁡(-mT))mb1(1-(1-δ)exp⁡(-mT)),r2T<u1ap(1-exp⁡(-mT))mb2(1-(1-δ)exp⁡(-mT)),r1-u1b1+M(pexp⁡(-mT)1-(1-δ)exp⁡(-mT))<0,r2-u1aM+aM+b2(pexp⁡(-mT)1-(1-δ)exp⁡(-mT))<0,
then the periodic solution (0,0,z*(t)) is globally asymptotically stable.Proof.
From Theorem2.8, the periodic solution (0,0,z*(t)) is locally asymptotically stable. Now it will be proved to be a global attractor. Let V(t)=x(t)+y(t), then by Theorem 2.7, there exists a constant M>0 such that x(t)≤M, y(t)≤M, and z(t)≤M. Obviously,
(2.21)V′|(1.1)≤x(t)(r1-u1z(t)b1+M)-c1x2(t)+y(t)(r2-u1az(t)M+aM+b2)-c2y2(t).
However, for system (1.1),
(2.22)dz(t)dt≥-mz(t),t≠nT,z(t+)=(1-δ)z(t)+p,t=nT.
Therefore, there exists aτ>0, and an ε>0 small enough, such that z(t)≥z1*(t)-ε for all t>τ. This leads to
(2.23)z(t)≥z1*(t)-ε=pexp⁡(-mT)1-(1-δ)exp⁡(-mT)-ε.
Letγ=Δpexp(-mT)/(1-(1-δ)exp(-mT))-ε; it is well known that r1-(u1γ/(b1+M))<0 and r2-(u1aγ/(M+aM+b2))<0. Therefore, when t>τ,
(2.24)V′|(1.1)≤x(t)(r1-u1z(t)b1+M)-c1x2(t)+y(t)(r2-u1az(t)M+aM+b2)-c2y2(t)<0.
Therefore,V(t)→0, x(t)→0, y(t)→0 as t→∞, and it follows that the periodic solution (0,0,z*(t)) is a global attractor. This completes the proof.Theorem 2.10.
The system (1.1) is permanent, if
(2.25)r1G0G1T>u1p(1-exp⁡(-mT))mb1(1-(1-δ)exp⁡(-mT)),r2T>u1ap(1-exp⁡(-mT))mb2(1-(1-δ)exp⁡(-mT)),r1-u1b1+M(pexp⁡(-mT)1-(1-δ)exp⁡(-mT))>0,r2-u1aM+aM+b2(pexp⁡(-mT)1-(1-δ)exp⁡(-mT))>0.Proof.
Assume thatX(t)=(x(t),y(t),z(t)) is any solution of system (1.1) with X(0)>0. From Theorem 2.7, x(t)≤M, y(t)≤M, and z(t)≤M with t≥0. From Lemma 2.5, z(t)≥z*(t)-ε for all t large enough and some t with z(t)>pexp(-mT)/(1-(1-δ)exp(-mT))-ε=ξ1. Therefore, it is necessary to find ξ2 and ξ3 such that x(t)≥ξ2 and y(t)≥ξ3 for t large enough. Select an ε1 small enough so that
(2.26)φ1=exp⁡(∫nT(n+1)T(r1G0G1-r1MG1-M-c1M-u1b1(v*(t)+ε1))dt)>1.
Next, it is necessary to prove thatx(t)≥ξ2 for t large enough, where ξ2 is a positive constant.
It is easy to prove that there exists at1∈(0,∞) such that x(t)≥ξ2. Otherwise, if x(t)<ξ2 for all t>0,
(2.27)dz(t)dt≤(u2aMb2+u3Mb1-m)z(t),t≠nT,z(t+)=(1-δ)z(t)+p,t=nT,z(0+)=z0,
and therefore z(t)≤v(t) and v(t)→v*(t), where v(t) is the solution of
(2.28)dv(t)dt=(u2aMb2+u3Mb1-m)v(t),t≠nT,v(t+)=(1-δ)v(t)+p,t=nT,v(0+)=v0,
with
(2.29)v*(t)=pexp⁡(-(m-(u2aM/b2)-(u3M/b1))(t-nT))1-(1-δ)exp⁡(-(m-(u2aM/b2)-(u3M/b1))T),t∈(nT,(n+1)T].
Therefore, there exists aT1>0 such that z(t)≤v(t)≤z*(t)+ε1, and
(2.30)dx(t)dt≥x(t)[r1G0G1-r1MG1-M-c1M-u1b1(v*(t)+ε1)],
for t>T1. Let N1∈N and N1T≥T2≥T1. Integrating (2.30) over (nT,(n+1)T]n>N1 yields
(2.31)x((n+1)T)≥x(nT+)exp⁡(∫nT(n+1)T(r1G0G1-r1MG1-M-c1M-u1b1(v*(t)+ε1))dt)=x(nT+),exp⁡(∫nT(n+1)T(r1G0G1-r1MG1-M-c1M-u1b1(v*(t)+ε1))dt)=x(nT+)φ1.
Therefore,x((N1+n)T)≥x(N1T)φ1n→∞ as n→∞, but this leads to a contradiction because x(t) is bounded. Therefore, there exists a t1>0 such that x(t1)≥ξ2.
Second, ifx(t1)≥ξ2 for all t≥t1, then the proof is complete. Hence, it is only necessary to consider those solutions which leave the region R={x(t):x(t)<ξ2} and reenter it. Let t*(t)=inft≥t1{x(t)<ξ2}, then x(t)≥ξ2, t∈(t,t*), and t*∈(n1T,(n1+1)T). However, obviously x(t*)=ξ2 because X(t)=(x(t),y(t),z(t)) is continuous. It is claimed here that there must exist a t2∈((n1+1)T,(n1+1)T+T1) such that x(t2)≥ξ2. Otherwise, x(t)<ξ2, for t∈((n1+1)T,(n1+1)T+T1),T1=n2T+n3T, n2,n3∈N, and then
(2.32)(n2-1)T>ln⁡(ε1/(M+p))-((u2aM/b2)+(u3M/b1)-m),exp⁡(φ2(n2+1)T)φ1n3>1.
Consider (2.28) with v(t*+)=z(t*+); this leads to
(2.33)v(t)=(v((n1+1)T+)-p1-(1-δ)exp⁡(-((u2aM/b2)+(u3M/b1)-m)T))×exp⁡(-(u2aMb2+u3Mb1-m)(t-(n1+1)T))+v*(t),
for t∈(nT,(n+1)T], n1+1<n<n1+n2+n3+1, and then
(2.34)|v(t)-v*(t)|<(M+p)exp⁡(-(u2aMb2+u3Mb1-m)(t-(n1+1)T))<ε1,z(t)≤v(t)≤v*(t)+ε1,
for t∈[(n1+n2+1)T,(n1+1)T+T1], which implies (2.30). Integrating (2.30) over t∈[(n1+n2+1)T,(n1+1)T+T1] yields
(2.35)x((n1+n2+n3+1)T)≥x((n1+n2+1)T)φ1n3.
For everyt∈(t*,(n1+1)T), there are only two possible cases: one, if x(t)<ξ2 for t∈(t*,(n1+1)T), when t∈(t*,(n1+n2+1)T), then x(t)<ξ2, and hence
(2.36)dx(t)dt≥x(t)[r1G0G1-r1ξ2G1-M-c1M-u1b1M]=φ2x(t).
Integrating (2.36) over t∈(t*,(n1+n2+1)T),x((n1+n2+1)T)≥x(t*)exp(φ2(n2+1)T) is obvious.
Therefore,x((n1+n2+n3+1)T)≥ξ2exp(φ2(n2+1)T)φ1n3>ξ2, but this leads to a contradiction. Hence, let t3=inft>t*{x(t)≥ξ2}, then x(t3)=ξ2, and (2.36) holds. Integrating (2.36) over [t*,t3),
(2.37)x(t)≥x(t*)exp⁡(φ2(t-t*))≥ξ2exp⁡(φ2(n1+n3+1)T)=ξ5.
x(t3)≥ξ2 is also true for t>t3. Hence, x(t)≥ξ5, t>t3. Two, there exists a t5∈(t*,(n1+1)T] such that x(t5)≥ξ2. Therefore, let t4=inft>t*{x(t)≥ξ2}, then x(t)<ξ2 for t∈[t*,t4) and x(t4)=ξ2. For t∈[t*,t4), (2.36) holds, and integrating it over t∈[t*,t4), x(t)≥x(t*)exp(φ2(t-t*))>ξ5. This process can be continued because x(t4)≥ξ2 and x(t)≥ξ5, t>t4. This yields the two possible cases, x(t)≥ξ5, t≥t1, and y(t)≥ξ4, t≥t2. Let Ω={(x,y,z):x(t)≥ξ2,y(t)≥ξ3,z(t)≥ξ1,x(t)+y(t)+z(t)≤3M}. Then the set Ω is a global attractor. Moreover, every solution X(t)=(x(t),y(t),z(t)) of system (1.1) will eventually enter and remain in the set Ω. Therefore, system (1.1) is permanent, and the proof is completed.


## 3. Numerical Analysis

### 3.1. Bifurcation Analysis

In Section2, it has been proved that the cyanobacteria and green algae-eradication periodic solution is locally asymptotically stable, and conditions have been established for the permanence of system (1.1). Now, to substantiate these theoretical results, the following parameters will be considered: u1=0.5, u2=0.6, u3= 0.4, r1=0.5, r2=0.5, G0=10, G1=15, ym=15, yn=20, c1=0.1, c2=0.12, a=0.75, b1=0.45, b2=0.6, δ=0.55, T=40 with initial values x0=0.5, y0=0.6, z0=0.7.From Theorem2.8, it is known that the periodic solution (0,0,z*(t)) is locally asymptotically stable if p>pmax=8.15; this periodic solution (0,0,z*(t)) is shown in Figure 1. It is clear that the variable fish z oscillates in a stable cycle, but the cyanobacteria population x and green algae population y rapidly decrease to zero when p>pmax. When p is smaller than pmax, then the cyanobacteria and green algae-eradication periodic solution will become unstable, and it is possible that the cyanobacteria population, green algae population, and fish population can coexist.Dynamic behavior of system (1.1). When p>pmax=8.15, the cyanobacteria and green algae will become extinct. Time series evolving according to biological control system (1.1): (a) the cyanobacteria population x, (b) the green algae population y, and (c) the fish population z.
(a)(b)(c)Figure2 shows typical bifurcation diagrams for system (1.1) with respect to p in the range p∈[0,10]. As p increases, system (1.1) exhibits rich dynamic behavior, such as chaotic bands with periodic windows, crises, period-halving bifurcations, quasi-periodic oscillations, and narrow or wide periodic windows. When p<0.89, the cyanobacteria population x enters a chaotic band. As p increases beyond 0.89, the cyanobacteria population x enters a chaotic band with periodic windows, crisis phenomena appear, and the chaotic attractor suddenly changes into a periodic attractor. Then the periodic attractor loses its stability, and system (1.1) again enters a chaotic band. Finally, the chaotic attractor changes into a periodic attractor again; details of these results are shown in Figure 3. When 1.15<p<2.7, there is a cascade of period-halving bifurcation which leads system (1.1) to a T-periodic solution (Figure 4), where the period-halving bifurcation is the opposite of the bifurcation observed earlier. When p>2.75, the green algae population y decreases to zero because of the increasing number of fish population z and the competition between the cyanobacteria population x and green algae population y. Then the cyanobacteria population x and fish population z can coexist in a stable cycle, but the numbers of the cyanobacteria population x will decrease because of short supply of resources. When p>pmax= 8.15 the cyanobacteria population x will be eradicated because the number of fish population z released is continually increasing, then the cyanobacteria and green algae-eradication periodic solution occurs. Once the cyanobacteria populationx and green algae population y have decreased to zero, fish population z will be eradicated after a period of time because of lack of food. Therefore, it is apparent that the value of p can effectively control the size of the cyanobacteria population and green algae population. All these numerical simulations are consistent with the theoretical proofs presented earlier.Bifurcation diagrams for system (1.1) showing the effect of p. (a) x versus p, (b) y versus p, and (c) z versus p.
(a)(b)(c)Crisis. (a) Chaotic attractor whenp=0.88, (b) phase diagram of periodic solution whenp=0.9, (c) chaotic attractor when p=1, and (d) phase diagram of periodic solution whenp=1.1.
(a)(b)(c)(d)Period-halving bifurcation. (a) Chaotic attractor whenp=1.3, (b) phase diagram of a 4T-periodic solution when p=1.55, (c) phase diagram of a 2T-periodic solution whenp=2, and (d) a period-T attractor whenp=2.5.
(a)(b)(c)(d)


### 3.2. Strange Attractors and Power Spectra

Now the commonly used method called power spectra will be used to study the qualitative nature of strange attractors [21]. By calculating the largest Lyapunov exponent for strange attractor (a) (Figure 5(a)), this value is determined to be 0.18219. Obviously, the strange attractor is a chaotic attractor. Furthermore, the spectrum of chaotic attractor (a) (Figure 5(b)) is composed of dense wave bands and sparse wave peaks with low wave troughs. These results agree with the observation that the chaotic attractor comes into being because some cycles lose weak stability.Strange attractor and power spectrum. (a) Strange attractor whenp=1.25, (b) power spectrum of attractor (a).
(a)(b)


### 3.3. The Largest Lyapunov Exponent

Convincing evidence for deterministic chaos has come from several recent experiments [22–24]. From these results, it is clear that chaos plays a very important role in these studies, and therefore detecting and exploring chaos is very important [25–31]. Therefore, the largest Lyapunov exponent is considered to be the most useful diagnostic tool for chaotic systems. The largest Lyapunov exponent λ must be positive for a chaotic attractor, otherwise, if λ is negative, the system will enter a stable state or become a periodic attractor. The largest Lyapunov exponent for system (1.1) was calculated for various values of p, and Figure 6 shows the results for p from 0 to 2.Figure 6
The largest Lyapunov exponents (LLE) for system (1.1) with p varying between 0 and 2.


## 3.1. Bifurcation Analysis

In Section2, it has been proved that the cyanobacteria and green algae-eradication periodic solution is locally asymptotically stable, and conditions have been established for the permanence of system (1.1). Now, to substantiate these theoretical results, the following parameters will be considered: u1=0.5, u2=0.6, u3= 0.4, r1=0.5, r2=0.5, G0=10, G1=15, ym=15, yn=20, c1=0.1, c2=0.12, a=0.75, b1=0.45, b2=0.6, δ=0.55, T=40 with initial values x0=0.5, y0=0.6, z0=0.7.From Theorem2.8, it is known that the periodic solution (0,0,z*(t)) is locally asymptotically stable if p>pmax=8.15; this periodic solution (0,0,z*(t)) is shown in Figure 1. It is clear that the variable fish z oscillates in a stable cycle, but the cyanobacteria population x and green algae population y rapidly decrease to zero when p>pmax. When p is smaller than pmax, then the cyanobacteria and green algae-eradication periodic solution will become unstable, and it is possible that the cyanobacteria population, green algae population, and fish population can coexist.Dynamic behavior of system (1.1). When p>pmax=8.15, the cyanobacteria and green algae will become extinct. Time series evolving according to biological control system (1.1): (a) the cyanobacteria population x, (b) the green algae population y, and (c) the fish population z.
(a)(b)(c)Figure2 shows typical bifurcation diagrams for system (1.1) with respect to p in the range p∈[0,10]. As p increases, system (1.1) exhibits rich dynamic behavior, such as chaotic bands with periodic windows, crises, period-halving bifurcations, quasi-periodic oscillations, and narrow or wide periodic windows. When p<0.89, the cyanobacteria population x enters a chaotic band. As p increases beyond 0.89, the cyanobacteria population x enters a chaotic band with periodic windows, crisis phenomena appear, and the chaotic attractor suddenly changes into a periodic attractor. Then the periodic attractor loses its stability, and system (1.1) again enters a chaotic band. Finally, the chaotic attractor changes into a periodic attractor again; details of these results are shown in Figure 3. When 1.15<p<2.7, there is a cascade of period-halving bifurcation which leads system (1.1) to a T-periodic solution (Figure 4), where the period-halving bifurcation is the opposite of the bifurcation observed earlier. When p>2.75, the green algae population y decreases to zero because of the increasing number of fish population z and the competition between the cyanobacteria population x and green algae population y. Then the cyanobacteria population x and fish population z can coexist in a stable cycle, but the numbers of the cyanobacteria population x will decrease because of short supply of resources. When p>pmax= 8.15 the cyanobacteria population x will be eradicated because the number of fish population z released is continually increasing, then the cyanobacteria and green algae-eradication periodic solution occurs. Once the cyanobacteria populationx and green algae population y have decreased to zero, fish population z will be eradicated after a period of time because of lack of food. Therefore, it is apparent that the value of p can effectively control the size of the cyanobacteria population and green algae population. All these numerical simulations are consistent with the theoretical proofs presented earlier.Bifurcation diagrams for system (1.1) showing the effect of p. (a) x versus p, (b) y versus p, and (c) z versus p.
(a)(b)(c)Crisis. (a) Chaotic attractor whenp=0.88, (b) phase diagram of periodic solution whenp=0.9, (c) chaotic attractor when p=1, and (d) phase diagram of periodic solution whenp=1.1.
(a)(b)(c)(d)Period-halving bifurcation. (a) Chaotic attractor whenp=1.3, (b) phase diagram of a 4T-periodic solution when p=1.55, (c) phase diagram of a 2T-periodic solution whenp=2, and (d) a period-T attractor whenp=2.5.
(a)(b)(c)(d)


## 3.2. Strange Attractors and Power Spectra

Now the commonly used method called power spectra will be used to study the qualitative nature of strange attractors [21]. By calculating the largest Lyapunov exponent for strange attractor (a) (Figure 5(a)), this value is determined to be 0.18219. Obviously, the strange attractor is a chaotic attractor. Furthermore, the spectrum of chaotic attractor (a) (Figure 5(b)) is composed of dense wave bands and sparse wave peaks with low wave troughs. These results agree with the observation that the chaotic attractor comes into being because some cycles lose weak stability.Strange attractor and power spectrum. (a) Strange attractor whenp=1.25, (b) power spectrum of attractor (a).
(a)(b)


## 3.3. The Largest Lyapunov Exponent

Convincing evidence for deterministic chaos has come from several recent experiments [22–24]. From these results, it is clear that chaos plays a very important role in these studies, and therefore detecting and exploring chaos is very important [25–31]. Therefore, the largest Lyapunov exponent is considered to be the most useful diagnostic tool for chaotic systems. The largest Lyapunov exponent λ must be positive for a chaotic attractor, otherwise, if λ is negative, the system will enter a stable state or become a periodic attractor. The largest Lyapunov exponent for system (1.1) was calculated for various values of p, and Figure 6 shows the results for p from 0 to 2.Figure 6
The largest Lyapunov exponents (LLE) for system (1.1) with p varying between 0 and 2.


## 4. Conclusions

In this paper, the effects of impulsive perturbations on a fish algae model have been investigated. It has been proved by means of the Floquet theory of impulsive differential equations and small-amplitude perturbation techniques that system (1.1) has a stable algae-eradication periodic solution. Furthermore, the conditions for system (1.1) to be permanent have been determined using the comparison theorem. Typical bifurcation diagrams have been analyzed in detail, revealing that the system exhibits very rich dynamics. From Theorem 2.8, it is clear that the cyanobacteria and green algae-eradication periodic solution (0,0,z*(t)) is locally asymptotically stable if p>pmax=8.15, but that the cyanobacteria population x and green algae population y rapidly decrease to zero. When the value of p>2.75, green algae population y decreases to zero, but the cyanobacteria population x and fish population z can coexist in a cycle; when 2.25<p<2.75, the cyanobacteria population, green algae population and fish population can coexist in a cycle. When 0<p<2.25, system (1.1) undergoes complex dynamics. From the above analysis, it is apparent that the numbers of both algae species can be controlled effectively using an impulsive control strategy.



---

*Source: 452789-2012-04-17.xml*