Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access November 4, 2022

Existence of positive periodic solutions for first-order nonlinear differential equations with multiple time-varying delays

  • Xiaoling Han and Ceyu Lei EMAIL logo
From the journal Open Mathematics

Abstract

This study elucidates the sufficient conditions for the first-order nonlinear differential equations with periodic coefficients and time-varying delays to have positive periodic solutions. Our results are proved using the Krasnosel’skii fixed point theorem. In this article, we have identified two sets Δ and and proved that at least one positive periodic solution exists in the interval between the point belonging to Δ and the point belonging to . We propose simple conditions that guarantee the existence of sets Δ and . In addition, we obtain the necessary conditions for the existence of positive periodic solutions of the first-order nonlinear differential equations when the periodic coefficients satisfy certain conditions. Finally, examples and numerical simulations are used to illustrate the validity of our results.

MSC 2010: 34A34; 34K13; 47H10; 92D25

1 Introduction

Consider the first-order nonlinear equation

(1.1) x ( t ) = a ( t ) x ( t ) + i = 1 n b i ( t ) f ( x ( t τ i ( t ) ) ) μ i = 1 n H i ( t ) x ( t τ i ( t ) ) ,

where t R . Let ω > 0 and μ [ 0 , 1 ] .

We assume that

(H1) a : R [ 0 , 1 ] and b i : R [ 0 , ) ( 1 i n ) are ω -periodic continuous functions satisfying a ( t ) = a ( t + ω ) and b i ( t ) = b i ( t + ω ) for t R ;

(H2) τ i : R [ 0 , ) ( 1 i n ) and H i : R [ 0 , ) ( 1 i n ) are ω -periodic continuous functions satisfying τ i ( t ) = τ i ( t + ω ) and H i ( t ) = H i ( t + ω ) for t R ;

(H3) f : [ 0 , ) [ 0 , ) is a Lipschitz continuous function, i.e., there is a non-negative constant L , for any x , y , the following inequality holds:

f ( x ) f ( y ) L x y .

Equation (1.1) is a delay differential equation. Delay differential equations are mainly used to describe dynamic systems that rely on current and past historical states. The time-delay phenomenon is widely used in the fields of population dynamics [1,2, 3,4] and infectious diseases [5,6]. Scholars have shown that delay differential equations can more accurately reflect the changing laws of objective things than differential equations without time delay. In the past few decades, delay differential equations have received the attention of many scholars and achieved certain results.

For equation (1.1), when the nonlinear term f ( x ) has different expressions, model (1.1) also has different application backgrounds. For example, when the nonlinear term f ( N ) = N e α N , then equation (1.1) will degenerate into the following Nicholson’s blowflies model [7] through certain assumptions:

N ( t ) = δ N ( t ) + B N ( t τ ) e γ N ( t τ ) ,

where N ( t ) is the size of the population at time t , B is the maximum per capita daily egg production rate, δ is per capita daily adult death rate, τ denotes the approximate time of the life cycle, and 1 γ is the size at which the population reproduces at its maximum rate. Nicholson’s blowflies model has received wide attention due to its extensive practical significance, and its theoretical achievements have made remarkable progress in the past few decades, see [8,9,10, 11,12]. For example, Li and Du [13] studied the existence of positive periodic solutions of the generalized Nicholson’s blowflies model,

x ( t ) = δ ( t ) x ( t ) + i = 1 m p i ( t ) x ( t τ i ( t ) ) e q i ( t ) x ( t τ i ( t ) ) ,

where m is a positive integer, p i , δ , q i C ( R + , ( 0 , ) ) , and τ i C ( R + , R + ) are T -periodic functions for i = 1 , 2 , 3 , , m .

When the nonlinear term f ( x ) = x m 1 + x l , then model (1.1) represents the hematopoiesis model. The term i = 1 n b i ( t ) x m ( t τ i ( t ) ) 1 + x l ( t τ i ( t ) ) represents the current production (density) of blood cells that is affected by the past blood cell count (density). The hematopoiesis model was proposed by Mackey and Glass in 1977 [14]. Since then, hematopoiesis models have been studied by many scholars [15,16,17, 18,19]. For example, Liu et al. [20] used the fixed point theorem to study the existence of the positive periodic solution of the following hematopoiesis model and gave the necessary conditions for the existence of a unique positive periodic solution,

x ( t ) = a ( t ) x ( t ) + i = 1 n b i ( t ) 1 + x m ( t τ i ( t ) ) ,

where m > 0 is a constant, a , b i , τ i are positive continuous ω -periodic functions.

In recent years, many authors have studied the existence of positive ω -periodic solutions for periodic ecological models with multiple delays. See [21,22, 23,24] and references therein. However, in this article, we will study the existence of positive periodic solutions of the first-order nonlinear differential equation with multiple time-varying delays. When nonlinear term f takes some special forms, problem (1.1) can degrade into some well-known models. Therefore, it is more universal than literature [7,13,20]. First, we give the existence condition of positive periodic solutions of equation (1.1), then give the uniqueness condition of positive periodic solution by using the contraction mapping principle, and finally give the oscillation condition of periodic solutions.

Here, we introduce the notations required to describe our main results. For each i = 1 , 2 , , n , let

b ̲ i = min 0 t ω b i ( t ) , b ¯ i = max 0 t ω b i ( t ) ,

H ̲ i = min 0 t ω H i ( t ) , H ¯ i = max 0 t ω H i ( t ) .

Define

α = 1 e 0 ω a ( r ) d r 1 and β = e 0 ω a ( r ) d r e 0 ω a ( r ) d r 1 .

Let

Δ = x R : α ω i = 1 n b ̲ i min ρ x u x f ( u ) α μ ω i = 1 n H ¯ i + 1 x > 0 and = x R : β ω i = 1 n b ¯ i max ρ x u x f ( u ) β ρ μ ω i = 1 n H ̲ i + 1 x ,

where

ρ = α β = 1 e 0 ω a ( r ) d r ( 0 , 1 ) .

In this article, we also assume that (H4) holds:

(H4) β ω i = 1 n b ¯ i > β ρ ω μ i = 1 n H ̲ i + 1 .

For simplicity, we refer to points belonging to Δ (resp., ) as Δ -points (resp., -points) of (1.1). Note that there are no positive real numbers belonging to both Δ and .

2 Auxiliary lemmas and preparations

For convenience, in this section, we would like to introduce some notations, definitions, lemmas, and assumptions which are used in what follows.

Definition 1

[25] Let M be a real Banach space. A nonempty, closed, convex set P M is a cone if it satisfies the following two conditions:

  1. x P , λ 0 imply λ x P ;

  2. x P , x P imply x = θ , where θ is the zero element of P .

Definition 2

[25] An operator K : M M is completely continuous if it is continuous and maps bounded sets into relatively compact set.

The following is the well-known Krasnoselskii’s fixed point theorem in a cone.

Lemma 1

[25] Let M be a Banach space, and let P M be a cone. Assume that Ω 1 and Ω 2 are open subset of M with θ Ω 1 Ω ¯ 1 Ω 2 , and let

K : P ( Ω ¯ 2 Ω 1 ) P

be a completely continuous operator such that

  1. K u u , u P Ω 1 , and K u u , u P Ω 2 ; or

  2. K u u , u P Ω 1 , and K u u , u P Ω 2 .

Then K has a fixed point in P ( Ω ¯ 2 Ω 1 ) .

Let

M = { x C ( R , R ) : x ( t + ω ) = x ( t ) for t R }

be the Banach space of ω -periodic continuous functions equipped with the norm

x = max t R x ( t ) = max t [ 0 , ω ] x ( t ) .

Define a subset in M by

P = { x M : x ( t ) ρ x , t [ 0 , ω ] } .

It is easy to see that P is a cone in M .

Lemma 2

The positive ω -periodic function x is a solution of (1.1) if and only if x M satisfies that

(2.1) x ( t ) = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t , s ) d s ,

where

G ( t , s ) = e t s a ( r ) d r e 0 ω a ( r ) d r 1 f o r t s t + ω .

Proof

“Only if” part: Let x M be a solution of equation (1.1). We have

( x ( t ) + a ( t ) x ( t ) ) e 0 t a ( r ) d r = i = 1 n b i ( t ) f ( x ( t τ i ( t ) ) ) μ i = 1 n H i ( t ) x ( t τ i ( t ) ) e 0 t a ( r ) d r ,

which is equivalent to

d d t x ( t ) e 0 t a ( r ) d r = i = 1 n b i ( t ) f ( x ( t τ i ( t ) ) ) μ i = 1 n H i ( t ) x ( t τ i ( t ) ) e 0 t a ( r ) d r .

The integration from t to t + ω gives

x ( t + ω ) e 0 t + ω a ( r ) d r x ( t ) e 0 t a ( r ) d r = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) e 0 s a ( r ) d r d s .

By the periodic properties, we obtain that

x ( t ) e 0 t a ( r ) d r e t t + ω a ( r ) d r 1 = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) e 0 s a ( r ) d r d s .

Thus,

x ( t ) = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) e t s a ( r ) d r e 0 ω a ( r ) d r 1 d s

for t R . Hence, we obtain the expression (2.1).

“If” part: Take the derivative of (2.1) with respect to t to obtain (1.1).

Thus, the proof of Lemma 2 is complete.□

Define a functional G : R × R [ 0 , ) by

G ( t , s ) = e t s a ( r ) d r e 0 ω a ( r ) d r 1 for t s t + ω .

It is clear that G is strictly increasing with respect to the second component s [ t , t + ω ] . Also, we see that

G ( t , t ) = 1 e 0 ω a ( r ) d r 1 = α

and

G ( t , t + ω ) = e 0 ω a ( r ) d r e 0 ω a ( r ) d r 1 = β .

Hence, when s [ t , t + ω ] and t R , we have

(2.2) α = G ( t , t ) G ( t , s ) G ( t , t + ω ) = β .

Moreover, G has the relationship

G ( t + ω , s + ω ) = e t + ω s + ω a ( r ) d r e 0 ω a ( r ) d r 1 = e t s a ( r ) d r e 0 ω a ( r ) d r 1 = G ( t , s ) .

Define an operator K : M M by:

(2.3) ( K x ) ( t ) = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t , s ) d s .

The operator K has the following properties.

Lemma 3

The operator K defined by (2.3) maps P into P.

Proof

Any element x of P is ω -periodic because P M . According to conditions (H1), (H2), and (2.3), we have

( K x ) ( t + ω ) = t + ω t + 2 ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t + ω , s ) d s = t t + ω i = 1 n b i ( s + ω ) f ( x ( s + ω τ i ( s + ω ) ) ) μ i = 1 n H i ( s + ω ) x ( s + ω τ i ( s + ω ) ) G ( t + ω , s + ω ) d s = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t , s ) d s = ( K x ) ( t ) .

This means that K x M . By (2.2), we have

(2.4) ( K x ) ( t ) β t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) d s

and

(2.5) ( K x ) ( t ) α t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) d s .

Form (2.4) and (2.5), we see that

( K x ) ( t ) α β K x = ρ K x .

Hence,

K P P .

The proof of Lemma 3 is complete.□

Lemma 4

The operator K : P P is completely continuous.

Proof

We need to verify the following two points:

  1. K is continuous;

  2. K maps any bounded subset of P into a relatively compact subset of P .

Point (i): Let { x n } converge to x in P . From (2.1), we obtain

( K x n ) ( t ) ( K x ) ( t ) = t t + ω i = 1 n b i ( s ) f ( x n ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x n ( s τ i ( s ) ) G ( t , s ) d s t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t , s ) d s β t t + ω i = 1 n b i ( s ) f ( x n ( s τ i ( s ) ) ) f ( x ( s τ i ( s ) ) ) d s + μ β t t + ω i = 1 n H i ( s ) x n ( s τ i ( s ) ) x ( s τ i ( s ) ) d s L β ω i = 1 n b ¯ i + μ β ω i = 1 n H ¯ i x n x ,

that is,

( K x n ) ( t ) ( K x ) ( t ) 0 as n .

Therefore, operator K is continuous.

Point (ii): For any x P , t R , we have

( K x ) ( t ) = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t , s ) d s β t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) d s + t t + ω μ i = 1 n H i ( s ) x ( s τ i ( s ) ) d s β ω i = 1 n ( b ¯ i f + μ H ¯ i x i ) < + ,

where

f = max x P f ( x ) , x i = max t [ 0 , ω ] x ( t τ i ( t ) ) .

Therefore, K P is uniformly bounded.

On the other hand, let x P and t 1 , t 2 [ 0 , ω ] with t 1 < t 2 . We have

( K x ) ( t 2 ) ( K x ) ( t 1 ) = t 2 t 2 + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t 2 , s ) d s t 1 t 1 + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t 1 , s ) d s = t 2 t 1 i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t 2 , s ) d s

+ t 1 t 1 + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t 2 , s ) d s + t 1 + ω t 2 + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t 2 , s ) d s t 1 t 1 + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t 1 , s ) d s 2 β i = 1 n ( b ¯ i f + μ H ¯ i x i ) t 2 t 1 + β i = 1 n ( b ¯ i f + μ H ¯ i x i ) t 1 t 1 + ω G ( t 2 , s ) G ( t 1 , s ) d s .

We can write

G ( t 2 , s ) G ( t 1 , s ) = 1 e 0 ω a ( r ) d r 1 e t 2 s a ( r ) d r e t 1 s a ( r ) d r e t 2 s a ( r ) d r e 0 ω a ( r ) d r 1 1 e t 1 t 2 a ( r ) d r .

This immediately implies that

t 1 t 1 + ω G ( t 2 , s ) G ( t 1 , s ) d s t 1 t 1 + ω e t 2 s a ( r ) d r e 0 ω a ( r ) d r 1 1 e t 1 t 2 a ( r ) d r d s ω β t 2 t 1 .

Therefore, for any x P , we obtain

( K x ) ( t 2 ) ( K x ) ( t 1 ) 2 β i = 1 n ( b ¯ i f + μ H ¯ i x i ) t 2 t 1 + β i = 1 n ( b ¯ i f + μ H ¯ i x i ) ω β t 2 t 1 2 β i = 1 n ( b ¯ i f + μ H ¯ i x i ) + β i = 1 n ( b ¯ i f + μ H ¯ i x i ) ω β t 2 t 1 β i = 1 n ( b ¯ i f + μ H ¯ i x i ) ( 2 + ω β ) t 2 t 1 ,

which implies that the operator K is equicontinuous. By using the Ascoli-Arzela theorem, K : P P is relatively compact. Hence, K : P P is completely continuous.□

3 Simple conditions to ensure the existence of -point and Δ -point

We can prove that (1.1) has sufficient conditions for -point and Δ -point. For -point, we can choose ε such that

β ω i = 1 n b ¯ i β ρ ω μ i = 1 n H ̲ i + 1 < 1 ε .

If function f satisfies

(3.1) lim x 0 f ( x ) x = 0 ,

then we can find a sufficiently small number x > 0 such that f ( x ) < ε x for 0 x x . Since f is continuous on [ 0 , ) , we can see that f ( x ) < ε x . Hence, from (3.1), we have

β ω i = 1 n b ¯ i max ρ x x x f ( x ) β ρ ω μ i = 1 n H ̲ i + 1 < 1 ε max ρ x x x f ( x ) x .

Hence, we see that x .

For Δ -point, we can choose G such that

α ω i = 1 n b ̲ i α ω μ i = 1 n H ¯ i + 1 > 1 G .

If for given G , there exists M > 0 , so that when x > M = ρ x , there is f ( x ) > G x . Then, we have

α ω i = 1 n b ̲ i min ρ x x x f ( x ) α ω μ i = 1 n H ¯ i + 1 > 1 G min ρ x x x f ( x ) x .

Hence, we see that x Δ .

4 Main results

Theorem 1

Under conditions (H1)–(H3), if there exist a Δ -point x and a -point x of (1.1), then equation (1.1) has at least one positive ω -periodic solution x ^ satisfying min { x , x } x ^ max { x , x } .

Proof

To apply Lemma 1, we have to only find open bounded subsets Ω 1 and Ω 2 of M that satisfy either conditions (i) and (ii). Here, we define the subsets Ω 1 and Ω 2 of M by

Ω 1 = { x M : x < min { x , x } }

and

Ω 2 = { x M : x < max { x , x } } ,

respectively. Any -point x and any Δ -point x never have the same value. Hence, the inclusion relation θ Ω 1 Ω ¯ 1 Ω 2 holds. And from Lemma 4, K : P ( Ω ¯ 2 Ω 1 ) P is completely continuous. There are two cases to consider depending on which is larger, x or x . We present the proof only of the case x < x because the proofs of both cases are essentially the same.

Any ω -periodic function x of P Ω 1 M satisfies x ( t ) ρ x = ρ x for t R . Hence, we have

(4.1) ρ x x ( s τ i ( s ) ) x for s R and i = 1 , 2 , , n .

By (2.2) and (4.1), we obtain

( K x ) ( t ) α t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) d s α ω i = 1 n b ̲ i min ρ x x x f ( x ) α ω μ i = 1 n H ¯ i x x

for t R . From these inequalities, we can see that K x x = x for x P Ω 1 .

Any ω -periodic function x of P Ω 2 M satisfies x ( t ) ρ x = ρ x for t R . Hence, we have

(4.2) ρ x x ( s τ i ( s ) ) x for s R and i = 1 , 2 , , n .

By (2.2) and (4.2), we obtain

( K x ) ( t ) β t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) d s β ω i = 1 n b ¯ i max ρ x x x f ( x ) β ρ ω μ i = 1 n H ̲ i x x

for t R . From these inequalities, we can see that K x x = x for x P Ω 2 .

Thus, we have confirmed that assumption (ii) of Lemma 1 is satisfied in the case x x (the assumption (i) of Lemma 1 holds in the case x x ). From Lemma 1, we can conclude that the operator K has a fixed point x ^ in P ( Ω ¯ 2 Ω 1 ) . Hence, the fixed point x ^ is a positive ω -periodic solution that satisfied the properties

x ^ ( t ) ρ x ^ for t R and x x ^ x .

Then, x ^ is a positive ω -periodic solution to equation (1.1). The proof is complete.□

Theorem 2

If L β ω i = 1 n b ¯ i + μ β ω i = 1 n H ¯ i < 1 holds, then equation (1.1) has a unique positive ω -periodic solution x.

Proof

From Lemma 4, we know that K is a continuous operator from P to P . For any x , y P , we have

( K x ) ( t ) ( K y ) ( t ) = t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) x ( s τ i ( s ) ) G ( t , s ) d s t t + ω i = 1 n b i ( s ) f ( y ( s τ i ( s ) ) ) μ i = 1 n H i ( s ) y ( s τ i ( s ) ) G ( t , s ) d s β t t + ω i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) f ( y ( s τ i ( s ) ) ) d s + μ β t t + ω i = 1 n H i ( s ) x ( s τ i ( s ) ) y ( s τ i ( s ) ) d s L β ω i = 1 n b ¯ i + μ β ω i = 1 n H ¯ i x y .

According to condition L β ω i = 1 n b ¯ i + μ β ω i = 1 n H ¯ i < 1 , K is a compressed map. Thus, from the contraction mapping principle, equation (1.1) has a unique positive ω -periodic solution x . The proof of Theorem 2 is complete.□

5 Necessary condition

For convenience, let us make an assumption:

H inf ( t ) = lim x 0 + inf i I n μ i = 1 n H i ( t ) x ( t τ i ( t ) ) x ( t ) , H sup ( t ) = lim x 0 + sup i I n μ i = 1 n H i ( t ) x ( t τ i ( t ) ) x ( t ) ,

where I n = { 1 , 2 , 3 , , n } .

Theorem 3

Assume that (H1)–(H4) hold and that

(5.1) a ( t ) + H sup ( t ) a ( t ) + H inf ( t ) i = 1 n b i ( t )

for all t. Then every positive solution of equation (1.1) tends to zero as t .

Proof

Let x ( t ) be any positive solution of equation (1.1). Equation (1.1) can be changed into

d d t x ( t ) e 0 t a ( r ) + μ i = 1 n H i ( r ) x ( r τ i ( r ) ) x ( r ) d r = i = 1 n b i ( t ) f ( x ( t τ i ( t ) ) ) e 0 t [ a ( r ) + μ i = 1 n H i ( r ) x ( r τ i ( r ) ) x ( r ) ] d r .

Integrating the above from t 0 > 0 to t > t 0 , we have

x ( t ) = x ( t 0 ) e t 0 t a ( r ) + μ i = 1 n H i ( r ) x ( r τ i ( r ) ) x ( r ) d r + t 0 t i = 1 n b i ( s ) f ( x ( s τ i ( s ) ) ) e t s a ( r ) + μ i = 1 n H i ( r ) x ( r τ i ( r ) ) x ( r ) d r d s .

From (3.1) and (5.1),

(5.2) x ( t ) x ( t 0 ) e t 0 t a ( r ) + μ i = 1 n H i ( r ) x ( r τ i ( r ) ) x ( r ) d r + ε t 0 t i = 1 n b i ( s ) x ( s τ i ( s ) ) e t s [ a ( r ) + H sup ( r ) ] d r d s x ( t 0 ) e t 0 t [ a ( r ) + H inf ( r ) ] d r + ε t 0 t [ a ( s ) + H sup ( s ) ] x ( s τ i ( s ) ) e t s [ a ( r ) + H sup ( r ) ] d r d s .

Let ζ = lim t sup x ( t ) , then 0 ζ < . Below we prove that ζ = 0 . We divide it into three cases.

Case 1. When x ( t ) > 0 . Choose t 0 > 0 such that x ( t ) > 0 for t > t 0 . Then 0 < x ( t 0 τ i ( t 0 ) ) < x ( t τ i ( t ) ) < x ( t ) for t > t 0 . From (1.1),

0 < a ( t ) x ( t ) + i = 1 n b i ( t ) f ( x ( t τ i ( t ) ) ) μ i = 1 n H i ( t ) x ( t τ i ( t ) ) x ( t ) a ( t ) + μ i = 1 n H i ( t ) x ( t τ i ( t ) ) x ( t ) + i = 1 n b i ( t ) f ( x ( t τ i ( t ) ) ) x ( t ) [ a ( t ) + H inf ( t ) ] + ε i = 1 n b i ( t ) x ( t τ i ( t ) ) x ( t ) [ a ( t ) + H inf ( t ) ] + ε i = 1 n b i ( t ) x ( t ) x ( t ) [ ε i = 1 n b i ( t ) a ( t ) H inf ( t ) ] 0 .

This contradiction shows that Case 1 is impossible.

Case 2. When x ( t ) < 0 . Choose t 0 > 0 such that x ( t ) < 0 for t > t 0 . Then ζ < x ( t τ i ( t ) ) < x ( t 0 τ i ( t 0 ) ) for t > t 0 . From (5.1) and (5.2), we have

(5.3) x ( t ) x ( t 0 ) e t 0 t [ a ( r ) + H inf ( r ) ] d r + ε t 0 t [ a ( s ) + H sup ( s ) ] x ( s τ i ( s ) ) e t s [ a ( r ) + H sup ( r ) ] d r d s x ( t 0 ) e t 0 t [ a ( r ) + H inf ( r ) ] d r + ε max i I n x ( t 0 τ i ( t 0 ) ) 1 e t 0 t [ a ( r ) + H sup ( r ) ] d r .

Let t in (5.3), we obtain

ζ ε max i I n x ( t 0 τ i ( t 0 ) ) .

Again let t 0 in the above, we have that ζ ε ζ . Because ε ( 0 , 1 ) , which implies that ζ = 0 .

Case 3. When x ( t ) is oscillatory. In this case, there is t n with t n as n such that

x ( t n ) = 0 for n = 1 , 2 , , and lim n x ( t n ) = ζ .

From (1.1) and (3.1), we have

a ( t n ) x ( t n ) + μ i = 1 n H i ( t n ) x ( t n τ i ( t n ) ) = i = 1 n b i ( t n ) f ( x ( t n τ i ( t n ) ) ) ε i = 1 n b i ( t n ) x ( t n τ i ( t n ) ) .

Transforming the above formula, we have

(5.4) x ( t n ) [ a ( t n ) + H inf ( t n ) ] ε i = 1 n b i ( t n ) x ( t n τ i ( t n ) ) .

Choose l = l ( n ) I n such that

(5.5) x ( t n ) τ l ( t n ) = max i I n x ( t n τ i ( t n ) ) .

From (5.4) and (5.5), we have

(5.6) x ( t n ) [ a ( t n ) + H inf ( t n ) ] ε x ( t n τ l ( t n ) ) i = 1 n b i ( t n ) .

Set ξ = lim n sup x ( t n τ l ( t n ) ) , then ξ ζ . Thus, (5.6) becomes

ζ [ a ( t n ) + H inf ( t n ) ] ε ξ i = 1 n b i ( t n ) ε ζ i = 1 n b i ( t n ) ,

which is a contradiction. So we have that ζ = 0 . The proof is complete.□

From Theorem 3, we have the following results immediately.

Corollary 1

Let (H1)–(H4) and (5.1) hold. Then equation (1.1) has no positive ω -periodic solution.

Corollary 2

Let (H1)–(H4) hold, and let a ( t ) + H sup ( t ) < i = 1 n b i ( t ) . Then equation (1.1) has at least one positive ω -periodic solution.

6 Example and numerical simulation

In this section, we give an example to illustrate the correctness of our main results.

Example 1

Consider the delayed periodic Nicholson’s blowflies models with a time-varying delay:

(6.1) x ( t ) = 1 2 + 1 3 sin t x ( t ) + ( 5 + 2 sin t ) x ( t τ ( t ) ) e 1 5 ( 2 + cos t ) x ( t τ ( t ) ) ( 1 + sin t ) x ( t τ ( t ) ) .

Proof

Note that a ( t ) = 1 2 + 1 3 sin t , b ( t ) = 5 + 2 sin t , τ ( t ) = 1 cos t , H ( t ) = 2 + 2 sin t , μ = 1 2 , and ω = 2 π . There are

b ¯ ( t ) = 7 , b ̲ ( t ) = 3 , H ¯ ( t ) = 4 , H ̲ ( t ) = 0 .

Then

α =