# Permanence, Extinction, and Almost Periodic Solution of  a Nicholson's Blowflies Model with Feedback Control and  Time Delay

Authors: Haihui Wu; Shengbin Yu
Journal: Discrete Dynamics in Nature and Society
 (2013)

---

## Abstract

A Nicholson's blowflies model with feedback control and time delay is studied. By applying the comparison
theorem of the differential equation and fluctuation lemma and constructing a suitable Lyapunov functional, sufficient conditions which guarantee the permanence, extinction, and existence of a unique globally attractive positive almost periodic solution of the system are obtained. It is proved that the feedback control variable and time delay have no influence on the permanence and extinction of the system.


---

## Body

## 1. Introduction

Letf(t) be any continuous bounded function defined on [0,+∞); we set
(1)fl=inft≥0f(t),fu=supt≥0f(t).In order to describe the dynamics of Nicholson’s blowflies, Gurney et al. [1] proposed the following mathematical model in 1980:
(2)N˙(t)=-δN(t)+PN(t-τ)e-aN(t-τ),
where N(t) is the size of the population at time t, P is the maximum per capita daily egg production rate, (1/a) is the size at which the population reproduces at its maximum rate, δ is the per capita daily adult death rate, and τ is the generation time. Kulenović and Ladas [2], Győri and Ladas [3], and Győri and Trofimchuk [4] investigated the oscillatory behaviors of the solutions of (2). For the attractivity, Kulenović et al. [5] and So and Yu [6] have shown that, when P>δ, every positive solution N(t) of (2) tends to a positive equilibrium N*=(1/a)ln(P/δ) as t→∞ if
(3)(eδτ-1)(Pδ-1)<1.Reference [5] further showed that, for P≤δ, every nonnegative solution of (2) tends to zero as t→∞, and for P>δ, (2) is uniformly persistent. Furthermore, Li and Fan [7] considered the following nonautonomous equation:(4)x˙(t)=x(t)(-δ(t)+p(t)exp{-α(t)x(t)}),
where α(t),δ(t), and p(t) are all positive ω-periodic functions. The authors show that (4) has a unique globally attractive ω-periodic positive solution if
(5)p(t)>δ(t)fort∈[0,ω].
Their results improved the results of Saker and Agarwal [8] who considered system (4) with α(t)=a  (a is a constant).Recently, Wang and Fan [9] proposed the following discrete Nicholson’s blowflies model with feedback control:
(6)x(n+1)=x(n)exp{-δ(n)+p(n)exp{-α(n)x(n)}-c(n)μ(n)},Δμ(n)=-a(n)μ(t)+b(n)x(n-m).
Sufficient conditions are established for the permanence and the extinction of the system (6). They show that the bounded feedback terms do not have any influence on the permanence or extinction of (6). The authors in [9] also proposed the following continuous model:
(7)x˙(t)=x(t)(-δ(t)+p(t)exp{-α(t)x(t)}-c(t)μ(t)),μ˙(t)=-a(t)μ(t)+b(t)x(t-τ);
however, they did not discuss the dynamic behaviors of the system (7). Considering that continuous models can excellently show the dynamic behaviors of those populations who have a long life cycle, overlapping generations, and large quantity, sufficient conditions for the permanence, global attractivity, and the existence of a unique, globally attractive, strictly positive almost periodic solution of the system (7) with τ=0 are obtained by Yu [10]. As pointed out by Nindjin et al. [11], time delay plays an important role in many biological dynamical systems, being particularly relevant in ecology and a model with time delay is a more realistic approach to the understanding of dynamics. Hence, it is necessary to study the model (7) which contains time delay.In the following discussion, we always assume thatδ(t),p(t),α(t),c(t),a(t),b(t) are all continuous, positive almost periodic functions. Also, from the viewpoint of mathematical biology, we consider (7) together with the following initial conditions:
(8)x(θ)=φ(θ)≥0,θ∈[-τ,0]φ(0)>0,μ(θ)=ψ(θ)≥0,θ∈[-τ,0]ψ(0)>0,
where φ(s) and ψ(s) are continuous on [-τ,0]. It is not difficult to see that solutions of (7) and (8) are well defined for all t≥0 and satisfy
(9)x(t)>0,μ(t)>0,fort≥0.The aim of this paper is, by constructing a suitable Lyapunov functional and applying the analysis technique of Feng et al. [12], to obtain sufficient conditions for the existence of a unique globally attractive positive almost periodic solution of the system (7) with initial condition (8).This paper is organized as follows. In Section2, by applying the analysis technique of [13, 14] and Fluctuation lemma [15, 16], we present the permanence and the extinction of model (7) and (8). In Section 3, by constructing a suitable Lyapunov functional, a sufficient conditions for the existence of a unique globally attractive positive almost periodic solution of the system (7) and (8). Examples together with their numeric simulations are stated in Section 4. For more works on almost periodic solutions of the ecosystem with feedback control, one could refer to [17–23] and the references cited therein.


## 2. Permanence and Extinction

Now let us state several lemmas which will be useful in proving the main result of this section.Lemma 1 (see [13]).
Assume thata>0, b(t)>0 is a boundedness continuous function and x(0)>0. Further suppose that (i)
(10)x˙(t)≤-ax(t)+b(t)
then for all t≥s,
(11)x(t)≤x(t-s)exp{-as}+∫t-st‍b(τ)exp{a(τ-t)}dτ.
Particularly, if b(t) is bounded above with respect to M, then
(12)limsupt→+∞x(t)≤Ma.(ii)
Also suppose that(13)x˙(t)≥-ax(t)+b(t);
then for all t≥s,
(14)x(t)≥x(t-s)exp{-as}+∫t-st‍b(τ)exp{a(τ-t)}dτ.
Particularly, if b(t) is bounded above with respect to m, then
(15)liminft→+∞x(t)≥ma.Lemma 2 (see [20]).
Ifa>0, b>0 and x˙≥x(b-ax), when t≥0 and x(0)>0, one has
(16)liminft→+∞x(t)≥ba.
If a>0, b>0, and x˙≤x(b-ax), when t≥0 and x(0)>0, one has
(17)limsupt→+∞x(t)≤ba.Lemma 3.
Let(x(t),μ(t))T be any solution of system (7) with initial condition (8); there exists positive numbers M1 and M2, which are independent of the solution of the system, such that
(18)limsupt→+∞x(t)≤M1,limsupt→+∞μ(t)≤M2.Proof.
Let(x(t),μ(t))T be any solution of system (7) satisfying initial condition (8). Since xexp{-αlx(t)}≤(1/αle) for x>0, according to the positivity of solution and the first equation of system (7), for t≥0,
(19)x˙(t)≤x(t)(-δl+puexp{-αlx(t)})≤-δlx(t)+puαle,
where e is the mathematical constant.
By applying Lemma1(i) to (19), we have
(20)limsupt→+∞x(t)≤puδlαle=ΔM1.
Hence, there exists T1>0 such that
(21)x(t)≤2M1,∀t≥T1.
Equation (21) together with the second equation of (7) leads to
(22)μ˙(t)≤-alμ(t)+2buM1,∀t≥T1+τ.
Using Lemma 1(i) again, one has
(23)limsupt→+∞μ(t)≤2buM1al=ΔM2.
Obviously, Mi(i=1,2) are independent of the solution of system (7). Equations (20) and (23) show that the conclusion of Lemma 3 holds. The proof is completed.Lemma 4.
Assume that(H1)
(24)p(t)>δ(t),t≥0,holds. Then there exists positive constants m1 and m2, which are independent of the solution of system (7), such that
(25)liminft→+∞x(t)≥m1,liminft→+∞μ(t)≥m2.Proof.
Let(x(t),μ(t))T be any solution of system (7) satisfying initial condition (8). From Lemma 3, there exists a T2>T1+τ such that for all t≥T2,x(t)≤M, μ(t)≤M, where M=2max{M1,M2}. According to the first equation of system (7) and the positivity of solution, for t≥T2,
(26)x˙(t)=x(t)(-δ(t)+p(t)exp{-α(t)x(t)}-c(t)μ(t))≥x(t)(-δ(t)+p(t)exp{-α(t)M}-c(t)M)=Δ-Q(t)x(t),
where Q(t)=δ(t)-p(t)exp{-α(t)M}+c(t)M.
Integrating both sides of (26) from η(η≤t) to t leads to
(27)x(t)x(η)≥exp{-∫ηt‍Q(s)ds},
or
(28)x(η)≤x(t)exp{∫ηt‍Q(s)ds}.
Particularly, taking η=t-τ, one can get
(29)x(t-τ)≤x(t)exp{∫t-τt‍Q(s)ds}.
Substituting (29) into the second equation of system (7) leads to
(30)μ˙(t)≤-alμ(t)+bux(t)exp{∫t-τt‍Q(s)ds}.
Applying Lemma 1(i) to the above differential inequality, for 0≤s≤t, one has
(31)μ(t)≤μ(t-s)exp{-als}+∫t-st‍bux(η)exp{∫η-τη‍Q(u)du}exp{al(η-t)}dη≤from(28)μ(t-s)exp{-als}+∫t-st‍bux(t)exp{∫ηt‍Q(u)du}exp{∫η-τη‍Q(u)du}×exp{al(η-t)}dη≤μ(t-s)exp{-als}+bux(t)∫t-st‍exp{∫ηt‍Q(u)du}×exp{∫η-τη‍Q(u)du}dη,
where we used the fact maxη∈[t-s,t]exp{al(η-t)}=exp{0}=1.
Note that there exists aK, such that 2cuMexp{-als}<(β/2), as s≥K, where β=inft≥0(p(t)-δ(t)). In fact, we can choose K>(1/al)ln(4cuM/β). And so, fixing K, combined with (31), we can obtain
(32)μ(t)≤Mexp{-alK}+bux(t)∫t-Kt‍exp{∫ηt‍Q(u)du}×exp{∫η-τη‍Q(u)du}dη≤Mexp{-alK}+Dx(t),
for all t>T2+K, where D=supt≥T3(bu∫t-Kt‍exp{∫ηt‍Q(u)du}exp{∫η-τη‍Q(u)du}dη)>0.
Considering thate-x≥1-x, for x>0, from the first equation of system (7) and the positivity of the solution, for t>T2+K, we can get
(33)x˙(t)=x(t)(-δ(t)+p(t)exp{-α(t)x(t)}-c(t)μ(t))≥x(t)(-δ(t)+p(t)-p(t)α(t)x(t)-2cuMexp{-alK}-cuDx(t))≥x(t)(β2-(puαu+cuD)x(t)).
By Lemma 2, we have
(34)liminft→+∞x(t)≥β2(puαu+cuD)=Δm1.
Thus, there exists T3>T2+K such that for all t>T3,
(35)x(t)≥m12.
Equation (35) together with the second equation of (7) leads to
(36)μ˙(t)≥-auμ(t)+blm12,∀t>T3.
By applying Lemma 1(ii) to the above differential inequality, we have
(37)liminft→+∞μ(t)≥blm12au=Δm2.
Obviously, mi(i=1,2) are independent of the solution of system (7). Equations (34) and (37) show that the conclusion of Lemma 4 holds. The proof is completed.From Lemmas3–4 and the definition of permanence, we can obtain the following conclusion.Theorem 5.
Assume that(H1) holds; then system (7) with initial condition (8) is permanent.As a direct corollary of Theorem  2 in [24], from Theorem 5, we have the following.Corollary 6.
Suppose that(H1) holds; then system (7) admits at least one positive ω-periodic solution if δ(t),p(t),α(t),c(t),a(t),b(t) are all continuous positive ω-periodic functions.Theorem 7.
Let(x(t),μ(t))T be any positive solution of the system (7) with initial condition (8). Assume that
(38)∫0∞‍(p(t)-δ(t))dt=-∞orp(t)≤δ(t),t≥0
holds; then
(39)limt→+∞x(t)=0,limt→+∞μ(t)=0.Proof.
Firstly, from the the first equation of (7),
(40)x˙(t)≤x(t)(p(t)-δ(t)).
If the former case of (38) holds, then
(41)0<x(t)≤x(0)exp[∫0t‍(p(t)-δ(t))dt]⟶0,ast⟶∞,
which shows that limt→+∞x(t)=0.
If the latter case of (38) holds, from (40) we have x˙(t)<0 or x(t) is decreasing; therefore, limt→+∞x(t)=q∈[0,+∞). Hence limsupt→+∞=liminft→+∞x(t)=q. We only need to show that q=0. Otherwise, if q>0, then there exists a T4>0, such that x(t)>(q/2)  for t≥T4. According to the Fluctuation lemma, there exists a sequence ξn→∞ as n→∞ such that x˙(ξn)→0, x(ξn)→limsupt→∞=q, as n→∞. We can choose a large enough number N such that ξn>T4 for n>N; hence, x(ξn)>(q/2) for all n>N.
Forn>N, p(t)≤δ(t) together with the first equation of (7) leads to(42)x˙(ξn)≤x(ξn)(-δ(ξn)+p(ξn)exp{-α(ξn)x(ξn)})≤x(ξn)(-δ(ξn)+δ(ξn)exp{-α(ξn)q2)})≤x(ξn)(-1+exp{-αlq2)})δl.
Let n→∞; we obtain that 0≤q(-1+exp{-αl(q/2)})δl or exp{-αl(q/2))}>1 which is impossible. Hence, q=0 or
(43)limt→+∞x(t)=0.
Now, we come to prove that(44)limt→+∞μ(t)=0.
For any ϵ>0, according to (43), there exists a T5>0, such that
(45)x(t)<ϵ∀t>T5.
Then, for t>T5+τ,
(46)μ˙(t)≤-alμ(t)+buϵ.
Thus, by applying Lemma 1(i) to the above differential inequality, we have
(47)0<μ(t)≤buϵal
which implies that
(48)limt→+∞μ(t)=0.
The proof is complete.


## 3.  Existence of a Unique Almost Periodic Solution

Now, we give the definition of the almost periodic function.Definition 8 (see [25, 26]).
A functionf(t,x), where f is an m-vector, t is a real scalar, and x is an n-vector, is said to be almost periodic in t uniformly with respect to x∈X⊂Rn, if f(t,x) is continuous in t∈R and x∈X, and if for any ε>0, there is a constant l(ε)>0, such that in any interval of length l(ε), there exists τ such that the inequality
(49)∥f(t+τ)-f(t)∥=∑i=1m‍|fi(t+τ,x)-fi(t,x)|<ε
is satisfied for all t∈(-∞,+∞),x∈X. The number τ is called an ε-translation number of f(t,x).Definition 9 (see [25, 26]).
A functionf:R→R is said to be an asymptotically almost periodic function if there exists an almost periodic function q(t) and a continuous function r(t) such that
(50)f(t)=q(t)+r(t),t∈R,r(t)⟶0ast⟶∞.
We denote byS(E) the set of all solutions z(t)=(x(t),μ(t))T of system (7) satisfying m1≤x(t)≤M1,m2≤μ(t)≤M2 for all t∈R.Lemma 10.
One hasS(E)≠∅.Proof.
Sinceδ(t),p(t),α(t),c(t),a(t),b(t) are almost periodic functions, there exists a sequence {tn},   tn→∞ as n→∞ such that
(51)δ(t+tn)⟶δ(t),p(t+tn)⟶δ(t),α(t+tn)⟶δ(t),c(t+tn)⟶δ(t),a(t+tn)⟶δ(t),b(t+tn)⟶δ(t),
as n→∞ uniformly on R. Suppose z(t)=(x(t),μ(t))T is a solution of (7) satisfying m1≤x(t)≤M1, m2≤μ(t)≤M2 for t>T. Obviously, the sequence (z(t+tn)) is uniformly bounded and equicontinuous on each bounded subset of R. Therefore, by the Ascoli-Arzela theorem, there exists a subsequence of {tn}, which we still denote by {tn}, such that x(t+tn)→m(t),μ(t+tn)→n(t),asn→∞ uniformly on each bounded subset of R. For any T1∈R, we may assume that tn+T1≥T for all n. For t≥0, we have(52)x(t+tn+T1)-x(tn+T1)=∫T1t+T1‍x(s+tn)(-δ(s+tn)+p(s+tn)×exp{-α(s+tn)x(s+tn)}-c(s+tn)μ(s))ds,μ(t+tn+T1)-μ(tn+T1)=∫T1t+T1‍(-a(s+tn)μ(s+tn)+b(s+tn)x(s+tn-τ))ds.
Applying Lebesgue’s dominated convergence theorem and letting n→∞ in to previous equations, we obtain
(53)m(t+T1)-m(T1)=∫T1t+T1‍m(s)(-δ(s)+p(s)×exp{-α(s)m(s)}-c(s)n(s))ds,n(t+T1)-n(T1)=∫T1t+T1‍(-a(s)n(s)+b(s)m(s-τ))ds,
for all t≥0. Since T1∈R is arbitrarily given, (m(t),n(t))T is a solution of system (7) on R. It is clear that m1≤m(t)≤M1, m2≤n(t)≤M2 for t∈R. That is to say, (m(t),n(t))T∈S(E). This completes the proof.Lemma 11 (see [27]).
Letf be a nonnegative function defined on [0,+∞) such that f is integrable on [0,+∞) and is uniformly continuous on [0,+∞). Then, limt→+∞f(t)=0.Theorem 12.
In addition to(H1), further suppose that (H2) there exists a h>0, such that
(54)plαlexp(αuM1)-bu>h,al-cu>h,where M1 is defined in (23); then system (7) with initial conditions (8) is globally attractive. That is to say, for any two positive solutions, one has
(55)limt→+∞|x(t)-x*(t)|=0,limt→+∞|μ(t)-μ*(t)|=0.Proof.
Let(x*(t),μ*(t))T and (x(t),μ(t))T be any two positive solutions of system (7)-(8). Theorem 5 implies there exist positive constants T,mi, and Mi(i=1,2) such that for t≥T(56)m1≤x(t)≤M1,m1≤x*(t)≤M1,m2≤μ(t)≤M2,m2≤μ*(t)≤M2,
where mi and Mi(i=1,2) are defined in Lemma 3 and Lemma 4. Set V(t)=V1(t)+V2(t), where
(57)V1(t)=|lnx(t)-lnx*(t)|,V2(t)=|μ(t)-μ*(t)|+bu∫t-τt‍|x(u)-x*(u)|du.
Calculating the upper right derivatives of V1(t) along the solution of (7) leads to
(58)D+V1(t)=sgn[x(t)-x*(t)]×(p(t)(exp{-α(t)x(t)}-exp{-α(t)x*(t)})+c(t)(μ*(t)-μ(t)))=sgn[x(t)-x*(t)]×(p(t)(-α(t)exp{-ξ(t)}(x(t)-x*(t)))+c(t)(μ*(t)-μ(t)))≤-p(t)α(t)(exp{-ξ(t)}|x(t)-x*(t)|)+c(t)|μ(t)-μ*(t)|,
where we used the elementary mean value theorem of differential calculus and ξ(t) lies between α(t)x(t) and α(t)x*(t). Then, for t≥T, we have
(59)αlm1≤ξ(t)≤αuM1.
Hence, by (58), we can have
(60)D+V1(t)≤-plαlexp(αuM1)|x(t)-x*(t)|+cu|μ(t)-μ*(t)|.
Calculating the upper right derivatives of V2(t) along the solution of (7), one has
(61)D+V2(t)=sgn[μ(t)-μ*(t)]×(-a(t)(μ(t)-μ*(t))+b(t)(x(t-τ)-x*(t-τ)))+bu(|x(t)-x*(t)|-|x(t-τ)-x*(t-τ)|)≤-al|μ(t)-μ*(t)|+bu|x(t)-x*(t)|.
According to (58), (66), and condition (H2), we can obtain
(62)D+V(t)≤(bu-plαlexp(αuM1))|x(t)-x*(t)|+(cu-al)|μ(t)-μ*(t)|<h[|μ(t)-μ*(t)|+|x(t)-x*(t)|].
Integrating both sides of (62) from T to t leads to
(63)V(t)+h∫Tt‍[|μ(s)-μ*(s)|+|x(s)-x*(s)|]ds<V(T)<+∞,t≥T.
Then,
(64)∫Tt‍[|μ(s)-μ*(s)|+|x(s)-x*(s)|]ds<V(T)h<+∞,t≥T.
Hence, |μ(t)-μ*(t)|+|x(t)-x*(t)|∈L1([T,+∞)). By system (7) and Theorem 5, we get μ(t),μ*(t),x(t),x*(t), and their derivatives are bounded on [T,+∞), which implies that |μ(t)-μ*(t)|+|x(t)-x*(t)| is uniformly continuous on [T,+∞). By Lemma 11, we obtain
(65)limt→+∞|x(t)-x*(t)|=0,limt→+∞|μ(t)-μ*(t)|=0.
The proof of Theorem 12 is complete.Theorem 13.
Suppose all conditions of Theorem12 hold; then there exists a unique almost periodic solution of systems (7) and (8).Proof.
According to Lemma10, there exists a bounded positive solution u(t)=(u1(t),u2(t))T of (7) with initial condition (8). Then there exists a sequence {tk′},{tk′}→∞ as k→∞, such that (u1(t+tk′),u2(t+tk′))T is a solution of the following system:
(66)x˙(t)=x(t)(-δ(t+tk′)+p(t+tk′)exp{-α(t+tk′)x(t)}-c(t+tk′)μ(t)),μ˙(t)=-a(t+tk′)μ(t)+b(t+tk′)x(t-τ).
According to Theorem 5 and the fact that δ(t),p(t),α(t),c(t),a(t),b(t) are all continuous, positive almost periodic functions, we know that both {ui(t+tk′)}(i=1,2) and its derivative function {u˙(t+tk′)}(i=1,2) are uniformly bounded; thus, {ui(t+tk′)}(i=1,2) are uniformly bounded and equi-continuous. By Ascoli’s theorem, there exists a uniformly convergent subsequence {ui(t+tk)}⊆{ui(t+tk′)} such that for any ε>0, there exists a K(ε)>0 with the property that if m,k≥K(ε), then
(67)|ui(t+tm)-ui(t+tk)|<ε,i=1,2.
That is to say, ui(t)(i=1,2) are asymptotically almost periodic functions. Hence there exists two almost periodic functions ri(t+tk)(i=1,2) and two continuous functions si(t+tk)(i=1,2) such that
(68)ui(t+tk)=ri(t+tk)+si(t+tk),i=1,2,
where
(69)limk→+∞ri(t+tk)=ri(t),limk→+∞si(t+tk)=0,i=1,2,ri(t)(i=1,2) are also almost periodic functions.
Therefore,(70)limk→+∞ui(t+tk)=ri(t),i=1,2.
On the other hand,
(71)limk→+∞u˙i(t+tk)=limk→+∞limh→0ui(t+tk+h)-ui(t+tk)h=limh→0limk→+∞ui(t+tk+h)-ui(t+tk)h=limh→0ri(t+h)-ri(t)h.
So r˙i(t)(i=1,2) exist. Moreover,
(72)r˙1(t)=limk→+∞u˙1(t+tk)=limk→+∞{u1(t+tk)×(-δ(t+tk)-c(t+tk)u2(t+tk)+p(t+tk)×exp{-α(t+tk)u1(t+tk)})}=r1(t)(-δ(t)+p(t)exp{-α(t)r1(t)}-c(t)r2(t)),r˙2(t)=limk→+∞u˙2(t+tk)=limk→+∞{-a(t+tk)r2(t+tk)+b(t+tk)r1(t+tk-τ)}=-a(t)r2(t)+b(t)r1(t-τ).
These show that (r1(t),r2(t))T satisfied system (7). Hence, (r1(t),r2(t))T is a positive almost periodic solution of (7). Then, it follows from Theorem 12 that system (7) has a unique positive almost periodic solution. The proof is completed.Without the feedback terms, that isa(t)=0,b(t)=0,c(t)=0, and (7) becomes the following equation:
(73)x˙(t)=x(t)(-δ(t)+p(t)exp{-α(t)x(t)}).
Equation (73) with periodic coefficients has been studied by Li and Fan [7] and Saker and Agarwal [8] with α(t)=a. Since the periodic case is a special case of almost periodic, hence, as a direct corollary of Theorem 13, we have the following.Corollary 14.
Supposeδ(t),p(t),α(t),c(t),a(t),b(t) are all continuous positive ω-periodic functions and p(t)>δ(t) for t∈[0,ω]; then (73) has a unique globally attractive ω-periodic positive solution.Remark 15.
Li and Fan in [7] show that (73) has a unique globally attractive ω-periodic positive solution if p(t)>δ(t)fort∈[0,ω], which is the same as Corollary 14. Thus, Theorem 13 supplements and generalizes results in [7, 8].


## 4.  Examples and Numeric Simulations

Now we give several examples together with their numeric simulations to show the feasibility of our main results.Example 16.
Consider the following example:(74)x˙(t)=x(t)(-5-sin(5t)+10exp{-(30+cos(11t))x(t)}-(1+0.5sin(7t))μ(t)),μ˙(t)=-(2.5+0.5cos(7t))μ(t)+(5.8+0.2sin(3t))x(t-2).
In this case, corresponding to system (7), we have δ(t)=5+sin(5t),p(t)=10,α(t)=30+cos(11t),a(t)=2.5+0.5cos(7t),b(t)=5.8+0.2sin(3t),c(t)=1+0.5sin(7t),τ=2. According to the proof of Lemmas 3 and 4, one has
(75)M1=3pu2δlαle=0.04757,M2=3buM12al=0.214066,m1=pl-δu-cuM22plαu=0.005934,m2=blm12au=0.0055384.
Hence,
(76)p(t)>δ(t),plαlexp(αuM1)-bu≅1261.193896>0,al-cu=0.5>0.
Thus, all the conditions of Theorem 13 are satisfied, and so, there exists a unique almost periodic solution of systems (74). Figure 1 shows this property.Figure 1
Dynamics of thex(t)) and μ(t) of system (74) with the initial values (x(0),μ(0))T=(0.02,0.05)T,(0.03,0.03)T,(0.04,0.08)T,(0.005,0.1)T; here t∈[0,100].Example 17.
Consider the following example:(77)x˙(t)=x(t)(-2-1(t+1)2+2exp{-(30+cos(11t))x(t)}-(1+0.5sin(7t))μ(t)1(t+1)2),μ˙(t)=-(2.5+0.5cos(7t))μ(t)+(5.8+0.2sin(3t))x(t-2).
In this case, we have(78)p(t)<δ(t).
Hence, By Theorem 7, we know that any positive solution of system (77) satisfies limt→+∞x(t)=0,limt→+∞μ(t)=0. Numerical simulation also confirms our result (see Figure 2).Figure 2
The dynamic behavior of system (77) with initial condition (x(0),μ(0))T=(0.03,0.05)T,(0.02,0.04)T,(0.005,0.01)T,(0.002,0.003)T; here t∈[0,100].


## 5. Conclusion

In this paper, we consider a Nicholson’s blowflies model with feedback control and time delay. It is shown that feedback control variable and time delay have no influence on the permanence and extinction of the system. Also, by constructing a suitable Lyapunov functional, a set of sufficient conditions which ensure the existence of a unique globally attractive positive almost periodic solution of the system is established. Moreover, compared with the main result of the relative discrete model (see [9]), we can see that the continuous and discrete models have similar results on permanence and the extinction of the Nicholson’s blowflies model with feedback control and time delay. At the end of this paper, two examples together with their numerical simulations show the verification of our main results. Our results supplement and generalize the results in [7, 8].



---

*Source: 798961-2013-05-07.xml*