Abstract
The Beverton–Holt model is a classical population model which has been considered in the literature for the discrete-time case. Its continuous-time analogue is the well-known logistic model. In this paper, we consider a quantum calculus analogue of the Beverton–Holt equation. We use a recently introduced concept of periodic functions in quantum calculus in order to study the existence of periodic solutions of the Beverton–Holt q-difference equation. Moreover, we present proofs of quantum calculus versions of two so-called Cushing–Henson conjectures.
1. Introduction
The Beverton–Holt difference equation has wide applications in population growth [Citation1] and is given by
where ν>1, K(t)>0 for all t∈ℕ0, and x(0)>0. We call the sequence K the carrying capacity and ν the inherent growth rate [Citation8]. The periodically forced Beverton–Holt equation, which is obtained by letting the carrying capacity be a periodic positive sequence K(t) with period ω∈ℕ, i.e., K(t+ω)=K(t) for all t∈ℕ0, has been treated with the methods found in [Citation8–10]. For the Beverton–Holt dynamic equation on time scales, one article has been presented by Bohner and Warth [Citation7]. In [Citation7], a general Beverton–Holt equation is given, which reduces to Equation (1) in the discrete case and to the well-known logistic equation in the continuous case. The approach given in [Citation7] opened a new path to the study of the discrete Beverton–Holt equation, which was pursued by Bohner et al. in [Citation6]. The crucial idea in [Citation6,Citation7] was to rewrite Equation (1) as
and thus identify the continuous version of the discrete Beverton–Holt equation (2) as
which turned out to be the usual logistic equation.Footnote1 This approach was generalized to any so-called dynamic equation of the form
hence accommodating both the continuous and discrete equations (2) and (3). However, the restriction on the time scale
was that it should be periodic. Hence, ℕ0 and ℝ (and also hℕ0 with h>0) were allowed, but
for q>1 was not.
In this paper, we are filling this gap by studying a quantum calculus version of the Beverton–Holt equation, namely, a Beverton–Holt q-difference equation. This became possible by using a new definition of periodic functions in quantum calculus which was introduced by the authors in [Citation3, Definition 3.1] (see also [Citation4]). Using this concept, we are interested in seeking ω-periodic solutions of the Beverton–Holt q-difference equation given by
where a is 1-periodic and K is ω-periodic, and
Using this notation and also our Assumptions (7) below, we can easily rewrite Equation (5) as
One can now observe the similarity of the discrete (additive) recursion (1) and the quantum (multiplicative) recursion (6).
The set-up of this paper is as follows. Section 2 contains some preliminaries on quantum calculus. We approach the periodic solutions of the Beverton–Holt q-difference equation (5) by some strategies presented in Section 3. In Sections 4 and 5, we formulate and prove the first and the second Cushing–Henson conjectures for the q-difference equations case, respectively.
2. Some auxiliary results
Definition 2.1
We say that a function is regressive provided
The set of all regressive functions will be denoted by
. Moreover,
is called positively regressive and we write
provided
Definition 2.2 Exponential function
Let and
. The exponential function ep(·, t0) on
is defined by
Remark 2.3
See [[Citation5], Theorem 2.44] If then ep(t, t0)>0 for all t≥t0,
.
Definition 2.4 See [3, Definition 3.1]
A function is called ω-periodic if
Theorem 2.5 See [5, Theorem 2.36]
If then
(i) e0(t, s)=1 and ep(t, t)=1;
(ii) ep(t, s)=1/ep(s, t);
(iii) ep(t, s)ep(s, r)=ep(t, r);
(iv) ep(σ(t), s)=(1+μ(t)p(t))ep(t, s);
(v) (1/ep(·, s))Δ(t)=−p(t)/ep(σ(t), s).
The integral on is defined as follows.
Definition 2.6
Let m, n∈ℕ0 with m<n. For we define
Theorem 2.7 Integration by parts, see [5, Theorem 1.77]
For and
we have
and
Theorem 2.8 Jensen's inequality, see [12, Theorem 2.2]
Let and c, d∈ℝ. Suppose
and
. If F∈C((c, d), ℝ) is convex, then
If F is strictly convex, then ‘≤’ can be replaced by ‘<’.
3. The Beverton–Holt equation
Throughout this paper, we use the following assumptions and notation:
Note that Assumption (7) implies that 0<λ<1,
, and
Note also that β is well defined as q/λ ¬∈{−1, 1} since λ ¬∈{−q, q} due to 0<λ<1.
In the q-difference equation (5), we assume x(t)>0 for all and substitute
Then, using the quotient rule [Citation5, Theorem 1.20 (v)], Equation (5) becomes
The general solution of Equation (9) is given by applying variation of parameters [Citation5, Theorem 2.77] twice as
where
. Now, we require an ω-periodic solution x¯ of Equation (5). This means that x¯ satisfies x¯(t)=qωx¯(qω t) for all
. This in turn means that a solution ū=1/x¯ of Equation (9) satisfies
Lemma 3.1
Assume Assumption (7). If Equation (9) has a solution ū satisfying Equation (12), then
Proof
Assume Equation (9) has a solution ū satisfying Equation (12). Then,
Thus, ū satisfies the required initial condition. ▪
4. The first Cushing–Henson conjecture
Now we state and prove the first Cushing–Henson conjecture for the Beverton–Holt q-difference equation (5).
conjecture 4.1 First Cushing–Henson conjecture
The Beverton–Holt q-difference model (5) with an ω-periodic carrying capacity K has a unique ω-periodic solution x¯ that globally attracts all solutions.
Using Equation (10) and Lemma 3.1, the solution ū of Equation (9) can be written as
Theorem 4.2
Assume Assumption (7) and let ū be given by Equation (13). Then, x¯:=1/ū is an ω-periodic solution of the Beverton–Holt q-difference equation (5).
Proof
By Equation (11), we have
so that
since by putting t0=qm and t=qn, we have
Hence, ū satisfies Equation (12) and thus x¯ is indeed ω-periodic. ▪
Now we are ready to prove the validity of the first Cushing–Henson conjecture.
Theorem 4.3
Assume Assumption (7). The solution x¯ of Equation (5) given in Theorem 4.2 is globally attractive.
Proof
First note that K is bounded. Indeed, define
For any m∈ℕ0, there exist ℓ∈ℕ0 and 0≤k≤ω−1 such that m=ℓ ω+k, and thus
Now let x be any solution of Equation (5) with x(t)>0 for all
. We have
which due to [Citation2, Theorem 2] tends to zero as t→∞. ▪
5. The second Cushing–Henson conjecture
Now we state and prove the second Cushing–Henson conjecture for the Beverton–Holt q-difference equation (5).
conjecture 5.1 Second Cushing–Henson conjecture
The average of the ω-periodic solution x¯ of Equation (5) is strictly less than the average of the ω-periodic carrying capacity K times the constant (q−λ)/(1−λ).
In order to prove the second Cushing–Henson conjecture, we use the following series of auxiliary results.
Lemma 5.2
Assume Assumption (7). Then, for any we have
Proof
Using Theorems 2.5 and 2.7, we get
which shows Equation (14). ▪
Lemma 5.3
Assume Assumption (7). Then, for any we have
Proof
Using Theorems 2.5 and 2.7, we get
which shows Equation (15). ▪
Now note that Equation (13) implies that for t0≤t<qω t0, we have
where
Lemma 5.4
Assume Assumptions (7) and (17). Then, for t0≤s<qω t0, we have
Proof
Using Lemma 5.3 and β qωλ−ω−β−1=0, we obtain
which shows Equation (18). ▪
Lemma 5.5
Assume Assumptions (7) and (17). Then, for t0≤t<qω t0, we have
Proof
Using Lemma 5.2 and β qωλ−ω−β−1=0, we obtain
which shows Equation (19). ▪
Now we are ready to prove the validity of the second Cushing–Henson conjecture.
Theorem 5.6
Let x¯ be the unique ω-periodic solution of Equation (5). If ω≠1, then
Proof
Since K is ω-periodic with ω≠1, tK(t) cannot be a constant. In addition, F(x)=1/x is strictly convex. Thus, we may use Jensen's inequality (Theorem 2.8) for the single inequality in the forthcoming calculation to obtain
which shows Inequality (20). The proof is complete. ▪
Theorem 5.7
If K is 1-periodic, then Inequality (20) becomes an equality, i.e.,
Proof
Since K is 1-periodic, we have
Now it is easy to check that x¯ given by
is 1-periodic and satisfies
Hence, x¯ is the unique 1-periodic solution of Equation (5). Thus, (21) holds. ▪
Remark 5.8 Note that the factor
is not present in the statements of the Cushing–Henson conjectures for the continuous and the discrete cases. However, in q-calculus, the presence of such quantities is common, and when replacing q by 1 (i.e., letting q→1+) in these terms, the continuous case is usually recovered. Note that replacing q by 1 in the factor (22) yields the corresponding ‘continuous’ (and also ‘discrete’) factor 1.
Notes
Note that, as Jim Cushing points out, the terminology ‘discrete logistic equation’ differs, due to Robert May [Citation11], in a slight but essential way, from the discrete model (2).
References
- R.J.H. Beverton and S.J. Holt, On the Dynamics of Exploited Fish Population, Fishery Investigations (Great Britain, Ministry of Agriculture, Fisheries, and Food), Vol. 19, H. M. Stationery Office, London, 1957.
- M. Bohner, Some oscillation criteria for first order delay dynamic equations, Far East J. Appl. Math. 18(3) (2005), pp. 289–304
- M. Bohner and R. Chieochan, Floquet theory for q-difference equations, Sarajevo J. Math. 8(21) (2) (2012), pp. 1–12
- M. Bohner and R. Chieochan, Positive periodic solutions for higher-order functional q-difference equations, J. Appl. Funct. Anal. 8(1) (2013), pp. 14–22
- M. Bohner and A. Peterson, Dynamic Equations on Time Scales: An Introduction with Applications, Birkhäuser Boston, Boston, MA 2001.
- M. Bohner, S. Stević, and H. Warth, The Beverton–Holt difference equation, in Discrete Dynamics and Difference Equations, S. Elaydi, H. Oliveira, J. Ferreira, and J. Alves, eds., World Scientific, Hackensack, NJ, 2010, pp. 189–193
- M. Bohner and H. Warth, The Beverton–Holt dynamic equation, Appl. Anal. 86(8) (2007), pp. 1007–1015. doi: 10.1080/00036810701474140
- J. M. Cushing and S. M. Henson, A periodically forced Beverton–Holt equation, J. Difference Equ. Appl. 8(12) (2002), pp. 1119–1120. doi: 10.1080/1023619021000053980
- S. Elaydi and R.J. Sacker, Global stability of periodic orbits of nonautonomous difference equations in population biology and the Cushing–Henson conjectures, in Proceedings of the Eighth International Conference on Difference Equations and Applications, S. Elaydi, G. Ladas, B. Aulbach, and O. Dosly, eds., Chapman & Hall/CRC, Boca Raton, FL, 2005, pp. 113–126.
- S. Elaydi and R. J. Sacker, Nonautonomous Beverton–Holt equations and the Cushing-Henson conjectures, J. Difference Equ. Appl. 11(4–5) (2005), pp. 337–346. doi: 10.1080/10236190412331335418
- R. M. May, Simple mathematical models with very complicated dynamics, Nature 261 (1976), pp. 459–467. doi: 10.1038/261459a0
- F.-H. Wong, C.-C. Yeh and W.-C. Lian, An extension of Jensen's inequality on time scales, Adv. Dyn. Syst. Appl. 1(1) (2006), pp. 113–120