\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2013 (2013), No. 56, pp. 1--9.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2013 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2013/56\hfil 
 Existence of positive almost-periodic solutions]
{Existence of positive almost-periodic solutions for a
Nicholson's blowflies model}

\author[Q.-L. Liu, H.-S. Ding \hfil EJDE-2013/56\hfilneg]
{Qing-Long Liu, Hui-Sheng Ding}  

\address{Qing-Long Liu \newline
College of Mathematics and Information Science,
Jiangxi Normal University,
Nanchang, Jiangxi 330022, China}
\email{758155543@qq.com}

\address{Hui-Sheng Ding \newline
College of Mathematics and Information Science,
Jiangxi Normal University,
Nanchang, Jiangxi 330022, China}
\email{dinghs@mail.ustc.edu.cn}

\thanks{Submitted October 15, 2012. Published February 21, 2013.}
\subjclass[2000]{34C27, 34K14}
\keywords{Nicholson's blowflies model; almost periodic solution; 
\hfill\break\indent exponential stability}

\begin{abstract}
 This article concerns the existence of solutions to a first-order 
 differential equation with time-varying delay, which is known as a
 Nicholson's blowflies model.
 By using fixed a point theorem and the Lyapunov functional method,
 we establish  the existence and locally exponential stability
 of almost periodic solutions for the model.
 Also, we apply our results to two examples, for which some
 earlier results can not be applied.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{example}[theorem]{Example}
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks

\section{Introduction and preliminaries}


Nicholson \cite{nicholson} and Gurney \cite{Gurney} proposed the
following delay differential equation as a population model
\begin{equation}\label{5}
x'(t)=-\delta x(t)+Px(t-\tau)e^{-ax(t-\tau)},
\end{equation}
where $x(t)$ is the size of the population at time $t$, $P$ is the
maximum per capita daily egg production, $1/a$ is the size
at which the population reproduces at its maximum rate, $\delta$ is
the per capita daily adult death rate, $\tau$ is the generation
time.


Nicholson's blowflies model and its analogous equations have
attracted much attention. Especially, equation \eqref{5} and its
variants have been of great interest for many mathematicians. As
pointed out in \cite{Fink}, compared with periodic effects, almost
periodic effects are more frequent in many biological dynamic
systems. Therefore, recently, the existence and stability of almost
periodic solutions for  \eqref{5} and its variants have
attracted more and more attention. We refer the reader to
\cite{Alzabut1,Alzabut2,Berezansky,liu1,liu2,liu3,long} and
references therein for some recent development on such topic.

Motivated by the above works, we study further this topic, by investigate the 
Nicholson's blowflies model
\begin{equation}\label{modeleq}
x'(t)=-\alpha(t)x(t)+\sum^{m}_{j=1}\beta_{j}(t) x(t-\tau_j(t))e^{-\gamma_{j}(t) x(t-\tau_j(t))},
\end{equation}
where $m$ is a given positive integer, and $\alpha, \beta_{j},
\gamma_{j}, \tau_{j} : \mathbb{R}\to (0,+\infty)$ are almost
periodic functions, $j=1,2,\dots,m$.

In this article, for a bounded continuous
function $g$ on $\mathbb{R}$, we denote
$$
g^{+}=\sup_{t\in{\mathbb{R}}}g(t),\quad 
g^{-}=\inf_{t\in\mathbb{R}}g(t).
$$
In addition, for $j=1,2,\dots,m$, we assume that
$$
\alpha^{-} >0,\quad  \beta_{j}^{-} >0, \quad \gamma_{j}^{-} >0,\quad
 \tau_{j}^{-} >0.
$$
Next, we recall some notation and basic results about almost
periodic functions. For more details, we refer the reader to
\cite{Fink}.


\begin{definition} \label{def1.1} \rm
A continuous function $f$: $\mathbb{R}\to \mathbb{R}$ is
called almost periodic if for each $\varepsilon > 0$, there exists
$l(\varepsilon) > 0$ such that every interval $I$ of length
$l(\varepsilon)$ contains a number $\tau$ with the property that
$$
\sup_{t\in\mathbb{R}}|f(t+\tau)-f(t)|<\varepsilon.
$$
We denote the set of all such functions by $AP(\mathbb{R})$.
\end{definition}

\begin{definition} \label{def} \rm
Let $x(\cdot)$ and $Q(\cdot)$ be $n\times n$ continuous matrix defined on
 $\mathbb{R}$. The linear system
\begin{equation}\label{6}
x'(t)=Q(t)x(t)
\end{equation}
is said to admit an exponential dichotomy on $\mathbb{R}$ if there
exist positive constants $k,\alpha$ and a projection $P$ such that
\begin{gather*}
\|X(t)PX^{-1}(s)\|\leq ke^{-\alpha (t-s)},\quad t\geq s,\\
\|X(t)(I-P)X^{-1}(s)\|\leq ke^{-\alpha (s-t)},\quad t\leq s,
\end{gather*}
for a fundamental solution matrix $X(t)$ of \eqref{6}.
\end{definition}

\begin{lemma}\label{lemma}
If the linear system \eqref{6} admits an exponential dichotomy, then
the almost periodic system
$$
x'(t)=Q(t)x(t)+g(t)
$$
has a unique almost periodic solution $x(t)$ given by
$$
x(t)=\int^{t}_{-\infty}X(t)PX^{-1}(s)g(s)ds
-\int^{+\infty}_{t}X(t)(I-P)X^{-1}(s)g(s)ds.
$$
\end{lemma}

\section{Main results}

\subsection{Existence of almost periodic solutions}
For the next theorem we use the following two assumptions:
\begin{itemize}
\item[(H1)] $n_{2}\geq n_{1}\geq \frac{1}{\min_{1\leq j\leq m}\{\gamma_{j}^{+}\}}$,
 where
$$
n_2=\sum_{j=1}^{m}\frac{\beta^+_j}{\alpha^-\gamma^-_je},\quad 
n_1=\sum_{j=1}^{m}\frac{n_2\beta^-_je^{-\gamma^+_jn_2}}{\alpha^+}.
$$

\item[(H2)] $\sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\alpha^{-}}
 <\min\{\frac{e^k}{|1-k|},e^2\}$,
where $k=n_1 \min_{1\leq j\leq m}\{\gamma_{j}^{-}\}$.
Moreover, if $k=1$, we denote
$$
\min\big\{\frac{e^k}{|1-k|},e^2\big\}=e^2.
$$
\end{itemize}

\begin{theorem}\label{theorem1}
Under assumptions {\rm (H1)--(H2)}, Equation \eqref{modeleq} has a unique
almost periodic solution in
$$
\Omega=\{\varphi\in AP(\mathbb{R}) : n_1\leq \varphi(t)\leq n_2, 
\forall t\in\mathbb{R} \}.
$$
\end{theorem}

\begin{proof} 
For each $\varphi\in AP(\mathbb{R})$, we consider the  almost periodic 
differential equation
\begin{equation}\label{eq1}
 x'(t)=-\alpha(t)x(t)+\sum^{m}_{j=1}\beta_{j}(t)\varphi(t-\tau_j(t))
e^{-\gamma_{j}(t) \varphi(t-\tau_j(t))}.
\end{equation}
Since $\alpha^->0$, we know that
$$
x'(t)=-\alpha(t)x(t)
$$
admits an exponential dichotomy on $\mathbb{R}$ with $P=I$.
Combining this with Lemma \ref{lemma}, we conclude that
\eqref{eq1} has a unique almost periodic solution:
\begin{equation}\label{eq2}
x^\varphi(t)=\int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}(s)\varphi(s-\tau_j(s))
e^{-\gamma_j(s)\varphi(s-\tau_j(s))}\Big)ds.
\end{equation}
Now, we define a mapping $T$ on $AP(\mathbb{R})$ by
$$
(T\varphi)(t)=x^{\varphi}(t),\quad t\in \mathbb{R}.
 $$

Next, we show that $T(\Omega)\subset \Omega$. It suffices to
prove that
$$
n_1\leq (T\varphi)(t)\leq n_2
$$
for all $t\in\mathbb{R}$ and $\varphi\in\Omega$. Noticing that
\begin{equation}\label{1}
\sup_{x\geq 0}x e^{-\gamma^-_{j}x}=\frac{1}{\gamma^-_{j}e},\quad 1\leq j\leq
m,
\end{equation}
 for all $t\in\mathbb{R}$ and $\varphi\in\Omega$, we have
\begin{align*}
(T\varphi)(t)=x^\varphi(t)
&= \int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}(s)\varphi(s-\tau_j(s))
e^{-\gamma_j(s)\varphi(s-\tau_j(s))}\Big)ds
\\
&\leq  \int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}^{+}\varphi(s-\tau_j(s))
 e^{-\gamma_j^{-}\varphi(s-\tau_j(s))}\Big)ds
\\
&\leq  \int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\gamma^-_je}\Big)ds
\\
&\leq  \sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\gamma^-_je}
 \int^t_{-\infty}e^{-\int^{t}_{s}\alpha^-du}ds
\\
& = \sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\gamma^-_je}
 \int^t_{-\infty}e^{\alpha^-(s-t)}ds
\\
& = \sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\alpha^-\gamma^-_je}
 = n_2.
\end{align*}
On the other hand, by (H1), we know that
$n_1\geq\frac{1}{\min_{1\leq j\leq m}\{\gamma_j^+\}}$. Thus,
we have 
\begin{equation}\label{2}
\inf_{n_1\leq x\leq n_2}x
e^{-\gamma^+_jx}= n_2e^{-\gamma^+_jn_2},\quad 1\leq j\leq m,
\end{equation}
 which yields 
\begin{align*}
(T\varphi)(t)=x^\varphi(t)
&= \int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}(s)\varphi(s-\tau_j(s))
e^{-\gamma_j(s)\varphi(s-\tau_j(s))}\Big)ds
\\
&\geq  \int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}^{-}\varphi(s-\tau_j(s))
e^{-\gamma_j^{+}\varphi(s-\tau_j(s))}\Big)ds
\\
&\geq  \int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}^{-}n_2e^{-\gamma_j^{+}n_2}\Big)ds
\\
&\geq  \sum^{m}_{j=1}\beta_{j}^{-}n_2e^{-\gamma_j^{+}n_2}
 \int^t_{-\infty}e^{-\int^{t}_{s}\alpha^+du}ds
\\
&\geq  \sum^{m}_{j=1}\beta_{j}^{-}n_2e^{-\gamma_j^{+}n_2}
 \int^t_{-\infty}e^{\alpha^+(s-t)}ds
\\
& =  \sum^{m}_{j=1}\frac{n_2\beta_{j}^{-}e^{-\gamma_j^{+}n_2}}{\alpha^+}
 =  n_1
\end{align*}
for all $t\in\mathbb{R}$ and $\varphi\in\Omega$.

By a direct calculation, one can obtain
\begin{equation}\label{3}
|xe^{-x}-ye^{-y}|\leq
\max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}|x-y|,\quad
x,y\geq k.
\end{equation} 
Denoting
$$
\Gamma_j(s)=\gamma_j(s)\varphi(s-\tau_j(s))e^{-\gamma_j(s)\varphi(s-\tau_j(s))}
-\gamma_j(s)\psi(s-\tau_j(s))e^{-\gamma_j(s)\psi(s-\tau_j(s))},
$$
for $\varphi,\psi\in \Omega$, we have
\begin{align*}
&\|T\varphi-T\psi\|\\
& = \sup_{t\in\mathbb{R}}|(T\varphi)(t)-(T\psi)(t)|\\
&= \sup_{t\in\mathbb{R}}\Big|\int^t_{-\infty}
 e^{-\int^{t}_{s}\alpha(u)du}\sum^{m}_{j=1}
 \frac{\beta_{j}(s)}{\gamma_j(s)}\Gamma_j(s)ds\Big|\\
&\leq \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}
 \sup_{t\in\mathbb{R}}\int^t_{-\infty}e^{-\int^{t}_{s}\alpha^-du}
 \sum^{m}_{j=1}\beta_{j}^+ |\varphi(s-\tau_j(s))-\psi(s-\tau_j(s))|ds\\
&\leq  \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}
 \sum^{m}_{j=1}\beta_{j}^+\|\varphi-\psi\|\cdot\sup_{t\in\mathbb{R}}
 \int^t_{-\infty}e^{\alpha^-(s-t)}ds\\
&=\max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}
 \sum^{m}_{j=1}\frac{\beta_{j}^+}{\alpha^-}\|\varphi-\psi\|\\
&=\frac{\sum^{m}_{j=1}\frac{\beta_{j}^+}{\alpha^-}}{\min\{\frac{e^k}{|1-k|},e^2\}}
 \|\varphi-\psi\|.
\end{align*}
Now, by (H2), $T$ has a unique fixed point in $\Omega$, i.e., Equation
\eqref{modeleq} has a unique almost periodic solution in $\Omega$.
\end{proof}

\begin{remark}\label{rmk2.2} \rm
Compared with some earlier results (cf. \cite{liu1,long}), in
Theorem \ref{theorem1}, we removed the following two restrictive
conditions:
$$\gamma^-_j>1,\quad n_1\geq  \frac{1}{\min_{1\leq j\leq
m}\{\gamma_j^-\}}.
$$
\end{remark}

For the next corollary we use the assumption
\begin{itemize}
\item[(H3)] there exists $\lambda>0$ such that
$$
\frac{\sum^{m}_{j=1}
\|\beta_{j}\|_{S^\lambda}}{1-e^{-\lambda\alpha^-}}
<\min\{\frac{e^k}{|1-k|},e^2\},
$$
where 
$$
k=n_1\cdot \min_{1\leq j\leq
m}\{\gamma_{j}^{-}\},\quad  \mbox{and} \quad
\|\beta_{j}\|_{S^\lambda}=\sup_{t\in\mathbb{R}}\int^{t+\lambda}_t
\beta_{j}(s)ds.
$$
\end{itemize}

\begin{corollary}\label{theorem1-corollary}
Under assumptions {\rm (H1), (H3)}, Equation \eqref{modeleq} has a unique
almost periodic solution in $\Omega$.
\end{corollary}

\begin{proof}
Let $T$ be the self mapping on $\Omega$ in Theorem
\ref{theorem1}, i.e.,
 $$
(T\varphi)(t)=
\int^t_{-\infty}e^{-\int^{t}_{s}\alpha(u)du}
\Big(\sum^{m}_{j=1}\beta_{j}(s)\varphi(s-\tau_j(s))
 e^{-\gamma_j(s)\varphi(s-\tau_j(s))}\Big)ds
$$
for $\varphi\in \Omega$ and $t\in \mathbb{R}$. For all $t\in\mathbb{R}$ and
$\varphi,\psi\in \Omega$, by \eqref{3}, we obtain
\begin{align*}
&|(T\varphi)(t)-(T\psi)(t)|\\
&\leq \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}
\int^t_{-\infty}e^{-\alpha^-(t-s)}\sum^{m}_{j=1}\beta_{j}(s)
 |\varphi(s-\tau_j(s))-\psi(s-\tau_j(s))|ds
\\
&\leq \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}
 \|\varphi-\psi\|\cdot\sum^{m}_{j=1}\int^t_{-\infty}
e^{-\alpha^-(t-s)}\beta_{j}(s)ds
\\
&= \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}\|\varphi-\psi\|
\sum^{m}_{j=1}\int_0^{+\infty}e^{-\alpha^-s}\beta_{j}(t-s)ds
\\
&= \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}
 \|\varphi-\psi\|\cdot\sum^{m}_{j=1}
\sum^{+\infty}_{k=0}\int_{k\lambda}^{(k+1)\lambda}e^{-\alpha^-s}\beta_{j}(t-s)ds
\\
&\leq \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}\|\varphi-\psi\|
\sum^{m}_{j=1} \sum^{+\infty}_{k=0}e^{-\alpha^-k\lambda}
 \int_{k\lambda}^{(k+1)\lambda}\beta_{j}(t-s)ds
\\
&= \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}\|\varphi-\psi\|
\sum^{m}_{j=1} \sum^{+\infty}_{k=0}e^{-\alpha^-k\lambda}
\int_{t-k\lambda-\lambda}^{t-k\lambda}\beta_{j}(s)ds
\\
&\leq \max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}\|\varphi-\psi\|
\sum^{m}_{j=1} \sum^{+\infty}_{k=0}e^{-\alpha^-k\lambda}
\|\beta_{j}\|_{S^{\lambda}}
\\
&=\max\big\{\frac{|1-k|}{e^k},\frac{1}{e^2}\big\}\|\varphi-\psi\|
\frac{\sum^{m}_{j=1} \|\beta_{j}\|_{S^{\lambda}}}{1-e^{-\lambda\alpha^-}}.
\end{align*}
Then, by (H3), $T$ is a contraction, and thus \eqref{modeleq}
has a unique almost periodic solution in $\Omega$.
\end{proof}

\begin{remark} \label{rmk2.4}\rm
Note that in some cases, (H3) is weaker than (H2)
(see for example \ref{exam2}).
\end{remark}

\subsection{Locally exponential stability of the almost periodic solution}

\begin{theorem}\label{theorem2}
Suppose that {\rm (H1), (H2)} are
satisfied, $x(t)$ is the unique almost periodic solution of 
\eqref{modeleq} in $\Omega$, and $y(t)$ is an arbitrary solution of
\eqref{modeleq} with 
\begin{equation}\label{4}
y(t)\geq n_1,\quad \forall t\in[-r,+\infty),
\end{equation} 
where $r=\max_{1 \leq j \leq m}\{\tau_j^+\}$. Then, there exists 
a constant $\lambda>0$
such that
$$
|x(t)-y(t)|\leq Me^{-\lambda t},\quad \forall t\in[-r,+\infty), 
$$
where $M=\max_{-r\leq t \leq 0}|x(t)-y(t)|$.
\end{theorem}

\begin{proof}
Let $\widetilde{k}=\max\{\frac{|1-k|}{e^k},\frac{1}{e^2}\}$. Then,
by (H2), we have
$$
\sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\alpha^{-}}<\frac{1}{\widetilde{k}}\,.
$$
Thus, there exists $\lambda >0$ such that
\begin{equation}\label{maodun}
\lambda -\alpha^-+\widetilde{k}\sum^{m}_{j=1}\beta^+_je^{\lambda
r}<0.
\end{equation}
Next, setting $z(t)=x(t)-y(t)$, we consider Lyapunov functional
$$
V(t)=|x(t)-y(t)|e^{\lambda t}=|z(t)|e^{\lambda t}.
$$
It can see easily  that
$$
V(t)\leq M,\quad  t \in [-r,0].
$$
Now, we claim that
\begin{equation}\label{claim}
V(t)\leq M, \quad  t \in (0,+\infty).
\end{equation}
Otherwise, the set 
$$
S=\{t>0: V(t)>M\}\neq\emptyset.
$$
 We denote
$t_1=\inf S$. It is not difficult to prove that
$$
V(t_1)=M,\quad  V(t)\leq M,\quad \forall t\in(0,t_1).
$$
In addition, for each $\delta>0$, there exists
$t'\in(t_1,t_1+\delta)$ such that $V(t')\geq M$. It follows that
 $$
D^+V(t_1)=\inf_{\delta>0}\sup_{t\in(t_1,t_1+\delta)}{\frac{V(t)-V(t_1)}{t-t_1}}
\geq 0 .
$$
On the other hand, by using \eqref{3}, we obtain
\begin{align*}
D^+V(t_1) & =  D^+(|z(t_1)|e^{\lambda t_1})\\ 
&\leq  e^{\lambda t_1}D^+|z(t_1)|+\lambda |z(t_1)|e^{\lambda t_1}
\\ 
&\leq  e^{\lambda t_1}sgn(z(t_1))z'(t_1)+\lambda|z(t_1)|e^{\lambda t_1}
\\ 
&\leq  e^{\lambda t_1}\big\{-\alpha(t_1)z(t_1)sgn(z(t_1))+\sum^{m}_{j=1}\frac{\beta_{j}(t_1)}{\gamma_{j}(t_1)}|\Gamma_j(t_1)|\big\}+\lambda|z(t_1)|e^{\lambda t_1}
\\ 
&\leq  -\alpha(t_1)e^{\lambda t_1}|z(t_1)|+\widetilde{k}\sum^{m}_{j=1}\beta_{j}^+|z(t_1-\tau_j(t_1))|e^{\lambda t_1}+\lambda|z(t_1)|e^{\lambda t_1}
\\ 
&\leq  -\alpha^-M+\widetilde{k}\sum^{m}_{j=1}\beta_{j}^+Me^{\lambda\tau_j(t_1)}+\lambda M
\\
&\leq (\lambda-\alpha^- + \widetilde{k}\sum^{m}_{j=1}\beta_{j}^+e^{\lambda r})M,
\end{align*}
where
$$
\Gamma_j(t_1)=\gamma_j(t_1)x(t_1-\tau_j(t_1))
e^{-\gamma_j(t_1)x(t_1-\tau_j(t_1))}-\gamma_j(t_1)y(s-\tau_j(t_1))
e^{-\gamma_j(t_1)y(t_1-\tau_j(t_1))}.
$$
Noting that $D^+V(t_1)\geq 0$, we obtain
$$
\lambda-\alpha^- + \widetilde{k}\sum^{m}_{j=1}\beta_{j}^+e^{\lambda
r}\geq 0,
$$ 
which contradicts \eqref{maodun}. Hence,
\eqref{claim} holds; i.e.,
$$
|x(t)-y(t)|\leq Me^{-\lambda t},\quad t\in[-r,+\infty).
$$
This completes the proof.
\end{proof}


\subsection{Examples}


In this section, we give two examples to illustrate our results.


\begin{example}\label{exam1}\rm
We consider the following Nicholson's blowflies model with multiple
time-varying delays:
\begin{equation}
\label{lizi2}
x'(t)=-\alpha(t)x(t)+\sum^{m}_{j=1}\beta_{j}(t) x(t-\tau_j(t))
e^{-\gamma_{j}(t) x(t-\tau_j(t))},
\end{equation}
where
\begin{gather*}
m=2, \quad \alpha(t)=18+\frac{|\sin \sqrt{2}t+\sin \sqrt{3}t|}{2},\\
\beta_1(t)=e^{e-1}(10+0.005|\sin\sqrt{3}t+\sin \sqrt{2}{t}|),\\
\beta_2(t)=e^{e-1}(10+0.005|\sin \sqrt{5}{t}+\sin \sqrt{3}{t}|),\\
\tau_1(t)=e^{|\cos\sqrt{2}t+\cos t|},\ \tau_2(t)
 =e^{|\cos \sqrt{2} {t}+\cos\sqrt{3}t|},\\
\gamma_1(t)=\gamma_2(t)=0.25+0.025|\sin t+\sin\sqrt{3}t|.
\end{gather*}
By a direct calculation, we can obtain
\begin{gather*}
\alpha^-=18,\quad \alpha^+=19,\quad
\beta^-_1=\beta^-_2=10e^{e-1},\\
\beta^+_1=\beta^+_2=10.01e^{e-1},\quad
\gamma^-_1=\gamma^-_2=0.25,\quad \gamma^+_1=\gamma^+_2=0.3,\\
n_2=\sum^{2}_{j=1}\frac{\beta^+_j}{\alpha^-\gamma^-_je}
 =\frac{80.08e^{e-2}}{18}>9,\\
3.4<n_1=\sum_{j=1}^{2}\frac{n_2\beta^-_je^{-\gamma^+_jn_2}}{\alpha^+}<3.5.
\end{gather*}
It is easy to see that 
$$
n_{2}\geq n_{1}\geq \frac{1}{\min_{1\leq j\leq
2}\{\gamma_{j}^{+}\}}=\frac{10}{3}.
$$ 
So (H1) holds. In addition, we have
\begin{gather*}
\sum^{2}_{j=1}\frac{\beta^{+}_{j}}{\alpha^{-}}=\frac{20.02e^{e-1}}{18}<e^2,\\
\frac{\exp(0.25n_1)}{1-0.25n_1}>e^2.
\end{gather*}
Thus, (H2) holds. By Theorem \ref{theorem1} and Theorem
\ref{theorem2}, Equation \eqref{lizi2} has a unique almost periodic
solution $x(t)$ in $\Omega$, and every solution $y(t)$ satisfying
\eqref{4} converges exponentially to $x(t)$ as $t\to +\infty$.
\end{example}

\begin{remark} \label{rmk2.7}\rm
In Example \ref{lizi2}, $\gamma_1^-=\gamma_2^-< 1$. Thus, the
results in \cite{liu1} cannot be applied to Example \ref{lizi2}. In
addition, here we have
$$
n_1< \frac{1}{\min_{1\leq j\leq m}\{\gamma_j^-\}}=4.
$$ 
So the results in \cite{long} can not be
applied to Example \ref{exam1}.
\end{remark}

\begin{example}\label{exam2}\rm
Let $m=1$, $\alpha(t)\equiv 1$, $\gamma_1(t)\equiv 1$, and
$\tau_1(t)=1+|\sin t+\sin \pi t|$ for all $t\in \mathbb{R}$. Moreover,
$\beta_1$ is a continuous $\frac{1}{2}$-periodic function, which is
defined on $[0,\frac{1}{2}]$ by
$$
\beta_1(t)=\begin{cases}
e^2,& t\in [\frac{1}{4}-\varepsilon,\frac{1}{4}+\varepsilon],\\
e^{e-1}, & t\in [0,\frac{1}{4}-2\varepsilon]\cup 
 [\frac{1}{4}+2\varepsilon,\frac{1}{2}],\\
\mbox{line segments}& t\in [\frac{1}{4}-2\varepsilon,\frac{1}{4}-\varepsilon]
\cup[\frac{1}{4}+\varepsilon,\frac{1}{4}+2\varepsilon],
\end{cases}
$$
where $\varepsilon\in (0,\frac{1}{8})$ is a fixed constant.
It is easy to check that
$$
\alpha^+=\alpha_-=\gamma^+=\gamma_-=1,\quad n_2=e,\quad n_1=1,
$$
and 
$$
\|\beta_1\|_{S^{1/2}}=\int^{1/2}_0 \beta_1(s)ds\leq
4\varepsilon e^2+(\frac{1}{2}-4\varepsilon)e^{e-1}.
$$ 
Then we conclude that (H1) holds and (H3) holds for sufficiently small
$\varepsilon$ and $\lambda=\frac{1}{2}$ since
$$
\lim_{\varepsilon\to 0 }4\varepsilon
e^2+(\frac{1}{2}-4\varepsilon)e^{e-1}=\frac{1}{2}e^{e-1}<e^2(1-e^{-1/2}).
$$
Then, by Corollary \ref{theorem1-corollary}, we know that
$$
x'(t)=-x(t)+\beta_1(t)x(t-\tau(t))e^{-x(t-\tau(t))}
$$
has an almost periodic solution provided that $\varepsilon$ is
sufficiently small.
\end{example}

\begin{remark} \label{rmk2.9} \rm
In Example \ref{exam2}, we have
$$
\sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\alpha^{-}}=\beta_1^+=e^2.
$$
However, in \cite{liu1,long},
$$
\sum^{m}_{j=1}\frac{\beta^{+}_{j}}{\alpha^{-}}<e^2
$$
is a key assumption. So the results in \cite{liu1,long} can not be
applied to Example \ref{exam2}.
\end{remark}

\subsection*{Acknowledgements}
This work was supported by the following grants:
11101192 from the NSF of China,
211090  from Key Project of Chinese Ministry ofEducation,
20114BAB211002 from the NSF of Jiangxi Province,
and GJJ12173 from the Jiangxi Provincial Education Department.



\begin{thebibliography}{00}

\bibitem{Alzabut1} Jehad O. Alzabut;
\emph{Almost periodic solutions for an impulsive delay Nicholson's 
blowflies model}, J. Comput. Appl. Math. 234 (2010) 233--239.

\bibitem{Alzabut2} Jehad O. Alzabut, Y. Bolat, T. Abdeljawad;
\emph{Almost periodic dynamic of a discrete Nicholson's blowflies model
involving a linear harvesting term}, Advances in Difference Equations
2012, 2012:158.

\bibitem{Berezansky} L. Berezansky, E. Braverman, L. Idels;
\emph{Nicholson's blowflies differential equations revisited: main results 
and open problems}, Appl. Math. Modelling
34 (2010) 1405--1417.

\bibitem{liu1} W. Chen, B. W. Liu;
\emph{Positive almost periodic solution for a class of Nicholson's 
blowflies model with multiple time-varying delays}, J. Comput. Appl. Math.
235 (2011) 2090--2097.

\bibitem{Gurney} W. S. Gurney, S. P. Blythe, R. M. Nisbet;
\emph{Nicholsons blowflies} (revisited), Nature 287 (1980) 17--21.

\bibitem{Fink} A. M. Fink;
\emph{Almost Periodic Differential Equations}, 
in: Lecture Notes in Mathematics, vol. 377, Springer, Berlin, 1974.

\bibitem{liu2} B. W. Liu;
\emph{Global stability of a class of Nicholson's blowflies model with 
patch structure and multiple time-varying delays}, Nonlinear Anal. Real
World Appl. 11 (2010), 2557--2562

\bibitem{liu3} B. W. Liu;
\emph{The existence and uniqueness of positive periodic
solutions of Nicholson-type delay systems}, Nonlinear Anal. Real
World Appl. 12 (2011), 3145--3151.


\bibitem{long} F. Long;
\emph{Positive almost periodic solution for a class of Nicholson's blowflies
model with a linear harvesting term}, Nonlinear Anal. Real World
Appl. 13 (2012) 686--693.

\bibitem{nicholson} A. J. Nicholson;
\emph{An outline of the dynamics of animal populations},
 Aust. J. Zool. 2 (1954) 9--65.

\end{thebibliography}
\end{document}
