Predator and n-classes-of-prey model incorporating extended Holling type Ⅱ functional response for n different prey species

1.   Introduction

Fractional differential equations with various types of fractional derivatives arise in modeling some dynamical processes (see, for example, ^[15] for the globally projective synchronization of Caputo fractional-order quaternion-valued neural networks with discrete and distributed delays, ^[18] for the quasi-uniform synchronization issue for fractional-order neural networks with leakage and discrete delays and ^[11] for Mittag-Leffler stability and adaptive impulsive synchronization of fractional order neural networks in quaternion field). In contrast to the classical derivative the fractional derivative is nonlocal and it depends significantly on its lower limit. As it is mentioned in ^[13], this leads to some obstacles for studying impulsive fractional differential equations.

Since many phenomena are characterized by abrupt changes at certain moments it is important to consider differential equations with impulses. In the literature there are two main approaches used to introduce impulses to fractional equations:

(i) With a fixed lower limit of the fractional derivative at the initial time- the fractional derivative of the unknown function has a lower limit equal to the initial time point over the whole interval of study;

(ii) With a changeable lower limit of the fractional derivative at each time of impulse- the fractional derivative on each interval between two consecutive impulses is changed because the lower limit of the fractional derivative is equal to the time of impulse.

Both interpretations of impulses are based on corresponding interpretations of impulses in ordinary differential equation, which coincide in the case of integer derivatives. However this is not the case for fractional derivatives. In the literature both types of interpretations are discussed and studied for Caputo fractional differential equations of order ${\alpha}\in(0, 1)$. We refer the reader to the papers ^{[6,7,12,13,16]} as well as the monograph ^[3].

We note in the case of the Caputo fractional derivative there is a similarity between both the initial conditions and the impulsive condition between fractional equations and ordinary equations (see, for example, ^[10] concerning the impulsive control law for the Caputo delay fractional-order neural network model). However for Riemann-Liouville fractional differential equations both the initial condition and impulsive conditions have to be appropriately given (which is different in the case of ordinary derivatives as well as the case of Caputo fractional derivatives). Riemann-Liouville fractional differential equations are considered, for example, in ^[1,2] for integral presentation of solutions in the case of the fractional order ${\alpha}\in(0, 1)$, ^[5,8,17] for the case of the fractional order ${\alpha}\in(1, 2)$.

In ^[14] the authors studied the following coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives of the form

where ${\alpha}, {\beta}\in (1, 2]$, $I = [0, T]$, $\phi_1, \phi_2:I\times {\mathbb{R}} \times {\mathbb{R}}\to {\mathbb{R}}$ are continuous functions, $\Delta u(t_j) = u(t_j^+)-u(t_j^-), \ \ \Delta u'(t_j) = u'(t_j^+)-u'(t_j^-)$, $\Delta v(t_k) = v(t_k^+)-v(t_k^-), \ \ \Delta v'(t_j) = v'(t_k^+)-v'(t_k^-)$, where $u(t_j^+), v(t_k^+)$ and $v(t_j^-), v(t_k^-)$ are the right limits and left limits respectively, $\mathcal{E}_j, \mathcal{E}^*_j, \mathcal{E}_k, \mathcal{E}^*_k:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, and $\mathcal{D}^{{\alpha}}, \mathcal{I}^{{\alpha}}$ are the ${\alpha}$-order Riemann-Liouville fractional derivative and integral operators respectively and $\mathcal{D}^{{\beta}-2} = \mathcal{I}^{2-{\beta}}$.

Since fractional integrals and derivatives have memories, and their lower limits are very important we will use the notations $^{RL}D^{{\alpha}}_{a, t}$ and $I^{{\beta}}_{a, t}$, respectively, instead of $\mathcal{D}^{{\alpha}}$ and $\mathcal{I}^{{\beta}}$, i.e. the Riemann-Liouville fractional derivative is defined by (see, for example ^[9])

and the Riemann-Liouville fractional integral $I^{{\beta}}_{a, t}$ of order ${\alpha} > 0$ is defined by (see, for example, ^[9])

where $a\geq 0$, ${\beta} > 0$ are given numbers.

Note there are some unclear parts in the statement of coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives (1.1), such as:

a). The presence of two different integers $p$ and $q$ in (1.1) leads to different domains of both the unknown functions $u$ and $v$. For example, in Corollary 1 ^[14] the solutions $u(t)$ and $v(t)$ are defined on $[0, t_{p+1}]$ and $[0, t_{q+1}]$ respectively, which causes some problems in the definitions of formulas (3.7) or (3.8) (^[14]);

b). The impulsive functions $\mathcal{E}_j, \mathcal{E}_k, \ j = 1, 2, \dots, p, \ k = 1, 2, \dots, q$ are assumed different but they are not (it is clear for example, for $j = k = 1$). The same is about the functions $\mathcal{E}^*_j, \mathcal{E}^*_k, \ j = 1, 2, \dots, p, \ k = 1, 2, \dots, q$.

In this paper we sort out the above mentioned points by setting up the cleared statement of the boundary value problem with the Riemann-Liouville (RL) fractional integral for the impulsive Riemann-Liouville fractional differential equation studied in ^[14], and we prove a new the integral presentations of the solutions. To be more general, we study two different interpretations for the presence of impulses in fractional differential equations. The first one is the case of the fixed lower limit of the RL fractional derivative at the initial time 0 and the second one is the case of the changed lower limit of the fractional RL derivative at any point of impulse. In both cases the integral presentation of the solution is provided.

2.   RL-fractional derivatives with fixed lower limit

2.1.   Statement of the problem

Define two different sequences of points of impulses

where $p, q$ are given natural numbers.

We will consider the following nonlinear boundary value problem for the coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives with a lower limit at $0$

where ${\alpha}, {\beta}\in (1, 2]$, $I = [0, T]$, $\phi_1, \phi_2:I\times{\mathbb{R}}\times{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $\Delta u(t_j) = u(t_j^+)-u(t_j^-), \ \ \Delta u'(t_j) = u'(t_j^+)-u'(t_j^-)$, $\Delta v({\tau}_k) = v({\tau}_k^+)-v({\tau}_k^-), \ \ \Delta v'({\tau}_k) = v'({\tau}_k^+)-v'({\tau}_k^-)$, where $u(t_j^+), v({\tau}_k^+), u'(t_j^+), v'({\tau}_k^+)$ and $u(t_j^-)$, $v(t_k^-)$, $u'(t_j^-)$, $v'(t_k^-)$ are the right limits and left limits respectively, $\mathcal{E}_j, \mathcal{E}^*_j, \mathcal{S}_k, \mathcal{S}^*_k:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, and $^{RL}D^{{\beta}}_{0, t}, \mathcal{I}^{{\alpha}}_{0, t}$ are the ${\alpha}$-order Riemann-Liouville fractional derivative and integral operators, respectively, $\mu_i, \nu_i, u_k, v_k, \ i = 1, 2,$ are given constants.

In the statement of the problem (2.1) some parts of (1.1) are cleared: there are two different points of impulses; the lower limits of the fractional integrals and fractional integrals are written; different functions at different points of impulses are used.

In our proofs we will use the following well known properties for fractional integrals (see, for example ^[9]).

We will apply the following auxiliary result which is a generalization of the result in ^[4] for an arbitrary lower limit of the fractional derivative:

Lemma 1. (^[4]). The general solution of the Riemann-Liouville fractional differential equation

is given by

where $c_0, c_1, a\geq 0$ are arbitrary real constants.

2.2.   Integral presentation of the linear problem

We will consider an appropriate boundary value problem for a scalar impulsive linear equation, we will prove a formula for its solution and later we will apply it to obatin the main result.

Consider the following boundary value problem for the linear impulsive fractional differential equation with Riemann-Liouville derivatives of the form

where $f:\ [0, T]\to {\mathbb{R}}$ is a continuous function, $\mathcal{E}_j, \mathcal{E}^*_j:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $u_1, u_2\in{\mathbb{R}}$.

Lemma 2. The solution of (2.5) satisfies the integral equation

where

Proof. We will use induction.

For $t\in (0, t_1]$ we apply Lemma 1 with $a = 0$ and we get

and

From the initial condition $I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1$ and equalities (2.2), (2.7) we get

i.e. $c_1 = \frac{u_1}{\Gamma({\alpha}-1)}.$

For $t\in (t_1, t_2]$ by Lemma 1 with $a = 0$ we get

and

From the impulsive condition $u(t_1+0)-u(t_1-0) = \mathcal{E}_1(u(t_1))$ we obtain

and from the impulsive condition $u'(t_1+0)- u'(t_1) = \mathcal{E}_1^*(u(t_1))$ we get

Thus we get the linear system w.r.t. $b_0$ and $b_1$

or

Therefore,

Assume the integral presentation of $u(t)$ is correct on $(t_{j-1}, t_j]$, i.e

Denote

and

Let $t\in (t_j, t_{j+1}], \ j = 2, \dots, p,$. By Lemma 1 with $a = 0$ we get

and

From the impulsive conditions and the equality (2.14) we obtain the linear system w.r.t. $k_0$ and $k_1$

or

Therefore,

From the boundary condition $\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = u_2$ we get

and

From (2.18) we have (2.6).

We will give an example to illustrate the claim of Lemma 2.

Example 1. Consider the following boundary value problem for the scalar RL fractional differential equation with an impulse at $t = 1$

The solution of (2.19) satisfies the integral equation

where

Consider the boundary value problem for the nonlinear impulsive fractional integro-differential equations with Riemann-Liouville derivatives of the form

where $\phi_1:\ [0, T]\times {\mathbb{R}}\to {\mathbb{R}}$ is a continuous function, $\mathcal{E}_j, \mathcal{E}^*_j:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $u_1, u_2\in{\mathbb{R}}$.

Corollary 1. The solution of (2.21) satisfies the integral equation

where

The proof of Corollary 1 is similar to the one of Lemma 2 with $f(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t))$ and we omit it.

Remark 1. Corollary 1 and the integral presentation (2.22) correct Theorem 3.1 and the formula (3.2) ^[14]. The main mistake in the proof of formula (3.2) ^[14] is the incorrect application of Lemma 1 with $a = 0$ on $(t_1, t_2]$ and taking the lower limit of the integral in (3.5) ^[14] incorrectly at $t_1 (\sigma_1)$ instead of 0. A similar comment applies to all the other intervals $(t_j, t_{j+1}]$.

2.3.   Integral presentation of the problem (2.1)

Following the proof of Lemma 2 and the integral presentation (2.22) of problem (2.21), we have the following result:

Theorem 1. The solution of (2.1) satisfies the integral equations

and

where

The proof is similar to the one of Lemma 2 applied twice to each of the both components $u$ and $v$ of the coupled system (2.1) for impulsive points $t_i, \ i = 1, 2, \dots, p$ and ${\tau}_i, \ i = 1, 2, \dots, p$ and the functions $f(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ and $f(t) = \phi_2(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ respectively.

Remark 2. Note the integral presentation (2.24), (2.25) of the solutions of the coupled system is the correction of Corollary 1 and integral presentation (3.7), (3.8) in ^[14].

3.   RL fractional derivatives with changeable lower limits

3.1.   Statement of the problem with changeable lower limit of the RL fractional derivatives

Consider the following nonlinear boundary value problem for the coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives with lower limits at impulsive points $t_i, \ i = 0, 1, 2, \dots, p-1$ and ${\tau}_k, \ k = 0, 1, 2, \dots, q-1$, respectively,

where ${\alpha}, {\beta}\in (1, 2]$, $\phi_1, \phi_2:I\times{\mathbb{R}}\times{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $\mathcal{P}_j, \mathcal{Q}_j, \ j = 1, 2, \dots, p,$ and $\mathcal{P}_j^*, \mathcal{Q}_j^*, \ j = 1, 2, \dots, q,$ are real numbers, $^{RL}D^{{\alpha}}_{t_i, t}$ and $^{RL}D^{{\beta}}_{{\tau}_k, t}$ are the ${\alpha}$-order Riemann-Liouville fractional derivatives with lower limits at $t_i$ and ${\tau}_k$, respectively, $\mu_k, \nu_k, u_k, v_k, \ i = 1, 2,$ are given constants.

Remark 3. Note problem (3.1) differs from problem (2.1):

$-$ The lower limits of the RL fractional derivatives $^{RL}D^{{\alpha}}_{t_j, t}$ and $^{RL}D^{{\beta}}_{{\tau}_k, t}$ in (3.1) are changed at any time of impulse $t_j$ and ${\tau}_k$, respectively.

$-$ The impulsive conditions are changed in (3.1). This is because the values of the unknown functions after the impulse, $u(t_j+0)$ and $v({\tau}_k+0)$, respectively, are considered as initial values at that point. But the RL fractional derivative has a singularity at its lower limit. It requires the chang of the impulsive conditions for the unknown functions.

3.2.   Integral presentation of the linear problem

Consider the following boundary value problem for the scalar linear impulsive fractional equation with Riemann-Liouville derivatives of the form

where the function $f:\ [0, T]\to {\mathbb{R}}$ is a continuous function, $\mathcal{P}_j, \mathcal{Q}_j, \ j = 1, 2, \dots, p$ are real numbers, $u_1, u_2\in {\mathbb{R}}$.

Now we will provide an integral presentation of the solution of (3.2).

Lemma 3. The solution of (3.2) satisfies the integral equation

where

Proof. We will use induction.

For $t\in (0, t_1]$ similar to Lemma 2 we get

where $c_1 = \frac{u_1}{\Gamma({\alpha}-1)}.$

For $t\in (t_1, t_2]$ by Lemma 1 with $a = t_1$ we get

and

From the impulsive condition $I^{2-{\alpha}}_{t_1, t}u'(t)|_{t = t_1} = \mathcal{P}_1u(t_1)+\mathcal{Q}_1$, equalities (2.2) and $I^{{\alpha}}_{t_1, t}f(t))|_{t = t_1} = 0$ we obtain

and

Since the integral $I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3}$ does not converge, it follows that $b_1 = 0$ and

Thus,

Similarly, for $t\in (t_j, t_{j+1}], \ j = 1, 2, \dots, p,$ by Lemma 1 with $a = t_j$ we get

and

From the impulsive conditions we obtain $k_1 = 0$ and

From the boundary condition $\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = u_2$ we get

and

and we obtain (3.3).

Example 2. Consider the equation

Note the Eq (3.13) is similar to (2.19) but the lower limit of the fractional derivative is changed at the point of the impulse. The solution of (3.13) satisfies the integral equation

It is clear the change of the lower limits of the fractional derivatives has a huge influence on the solution of the equation.

3.3.   Integral presentation of the problem (3.1)

Based on the integral presentation of the linear problem (3.2) and Lemma 3, we obtain the following result:

Theorem 2. The solution of (3.1) satisfies the integral equations

and

where

The proof is similar to the one of Lemma 3 applied twice to each of the both components $u$ and $v$ of the coupled system (3.1) for impulsive points $t_i, \ i = 1, 2, \dots, p$ and ${\tau}_i, \ i = 1, 2, \dots, p$ and the functions $f(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ and $f(t) = \phi_2(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ respectively.

4.   Conclusions

In this paper we set up and study a scalar nonlinear integro-differential equation with Riemann-Liouville fractional derivative and impulses. We consider a boundary value problem for the studied equation with Riemann-Liouville fractional derivative of order in $(1, 2)$. Note for Riemann-Liouville fractional differential equations both the initial condition and impulsive conditions have to be appropriately given (which is different in the case of ordinary derivatives as well as the case of Caputo fractional derivatives). We consider both interpretations in the literature on the presence of impulses in fractional differential equations: With fixed lower limit of the fractional derivative at the initial time point and with lower limits changeable at each impulsive time point. In both cases we set up in an appropriate way impulsive conditions which are dependent on the Riemann-Liouville fractional derivative. We obtain integral presentations of the solutions in both cases. These presentations could be successfully used for furure studies of existence, stability and other qualitative properties of the solutions of the integro-differential equations with Riemann-Liouville fractional derivative of order in $(1, 2)$ and impulses.

Acknowledgments

S. H. is partially supported by the Bulgarian National Science Fund under Project KP-06-N32/7.

Conflict of interest

The authors declare that they have no competing interests.