An improved atomic search algorithm for optimization and application in ML DOA estimation of vector hydrophone array

Peng Wang; Weijia He; Fan Guo; Xuefang He; Jiajun Huang; Peng Wang; Weijia He; Fan Guo; Xuefang He; Jiajun Huang

doi:10.3934/math.2022308

AIMS Mathematics

2022, Volume 7, Issue 4: 5563-5593. doi: 10.3934/math.2022308

Previous Article Next Article

Research article

An improved atomic search algorithm for optimization and application in ML DOA estimation of vector hydrophone array

Department of Mathematics, North University of China, Taiyuan, Shanxi 030051, China

Received: 30 September 2021 Revised: 22 December 2021 Accepted: 03 January 2022 Published: 10 January 2022
MSC : 65K05, 65K10

The atom search optimization (ASO) algorithm has the characteristics of fewer parameters and better performance than the traditional intelligent optimization algorithms, but it is found that ASO may easily fall into local optimum and its accuracy is not higher. Therefore, based on the idea of speed update in particle swarm optimization (PSO), an improved atomic search optimization (IASO) algorithm is proposed in this paper. Compared with traditional ASO, IASO has a faster convergence speed and higher precision for 23 benchmark functions. IASO algorithm has been successfully applied to maximum likelihood (ML) estimator for the direction of arrival (DOA), under the conditions of the different number of signal sources, different signal-to-noise ratio (SNR) and different population size, the simulation results show that ML estimator with IASO algorithum has faster convergence speed, fewer iterations and lower root mean square error (RMSE) than ML estimator with ASO, sine cosine algorithm (SCA), genetic algorithm (GA) and particle swarm optimization (PSO). Therefore, the proposed algorithm holds great potential for not only guaranteeing the estimation accuracy but also greatly reducing the computational complexity of multidimensional nonlinear optimization of ML estimator.

Keywords:

Citation: Peng Wang, Weijia He, Fan Guo, Xuefang He, Jiajun Huang. An improved atomic search algorithm for optimization and application in ML DOA estimation of vector hydrophone array[J]. AIMS Mathematics, 2022, 7(4): 5563-5593. doi: 10.3934/math.2022308

Related Papers:

[1]	Sitalakshmi Venkatraman, Ramanathan Venkatraman . Big data security challenges and strategies. AIMS Mathematics, 2019, 4(3): 860-879. doi: 10.3934/math.2019.3.860
[2]	Atiqe Ur Rahman, Muhammad Saeed, Mazin Abed Mohammed, Alaa S Al-Waisy, Seifedine Kadry, Jungeun Kim . An innovative fuzzy parameterized MADM approach to site selection for dam construction based on sv-complex neutrosophic hypersoft set. AIMS Mathematics, 2023, 8(2): 4907-4929. doi: 10.3934/math.2023245
[3]	Pengyu Liu, Jie Jian . Effects of network topology and trait distribution on collective decision making. AIMS Mathematics, 2023, 8(5): 12287-12320. doi: 10.3934/math.2023619
[4]	Shahryar Sorooshian, Seyedh Mahboobeh Jamali, Nader Ale Ebrahim . Performance of the decision-making trial and evaluation laboratory. AIMS Mathematics, 2023, 8(3): 7490-7514. doi: 10.3934/math.2023376
[5]	Erdal Bayram, Gülşah Çelik, Mustafa Gezek . An advanced encryption system based on soft sets. AIMS Mathematics, 2024, 9(11): 32232-32256. doi: 10.3934/math.20241547
[6]	Esmail Hassan Abdullatif Al-Sabri, Muhammad Rahim, Fazli Amin, Rashad Ismail, Salma Khan, Agaeb Mahal Alanzi, Hamiden Abd El-Wahed Khalifa . Multi-criteria decision-making based on Pythagorean cubic fuzzy Einstein aggregation operators for investment management. AIMS Mathematics, 2023, 8(7): 16961-16988. doi: 10.3934/math.2023866
[7]	Muhammad Arshad, Muhammad Saeed, Khuram Ali Khan, Nehad Ali Shah, Wajaree Weera, Jae Dong Chung . A robust MADM-approach to recruitment-based pattern recognition by using similarity measures of interval-valued fuzzy hypersoft set. AIMS Mathematics, 2023, 8(5): 12321-12341. doi: 10.3934/math.2023620
[8]	Muhammad Ihsan, Muhammad Saeed, Atiqe Ur Rahman, Hüseyin Kamacı, Nehad Ali Shah, Wajaree Weera . An MADM-based fuzzy parameterized framework for solar panels evaluation in a fuzzy hypersoft expert set environment. AIMS Mathematics, 2023, 8(2): 3403-3427. doi: 10.3934/math.2023175
[9]	Esra Korkmaz . An innovative algorithm based on weighted fuzzy soft multisets and its application in selecting optimal construction materials. AIMS Mathematics, 2024, 9(10): 27512-27534. doi: 10.3934/math.20241336
[10]	Jia-Bao Liu, Rashad Ismail, Muhammad Kamran, Esmail Hassan Abdullatif Al-Sabri, Shahzaib Ashraf, Ismail Naci Cangul . An optimization strategy with SV-neutrosophic quaternion information and probabilistic hesitant fuzzy rough Einstein aggregation operator. AIMS Mathematics, 2023, 8(9): 20612-20653. doi: 10.3934/math.20231051

Abstract

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

$\begin{align} \left\{ \begin{array}{ll} \left\{ \begin{array}{ll} \mathcal{D}^{{\alpha}}u(t)-\phi_1(t, \mathcal{I}^{\alpha} u(t), \mathcal{I}^{\beta} v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = t_i, \ \ i = 1, 2, \dots, p, \\ \Delta u(t_j)-\mathcal{E}_j(u(t_j)) = 0, \ \ \ \ \ \ \Delta u'(t_j)-\mathcal{E}_j^*(u(t_j)) = 0, \ \ j = 1, 2, \dots, p, \\ \nu_1\mathcal{D}^{{\alpha}-2}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_2 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = u_2, \end{array} \right.\\ \left\{ \begin{array}{ll} \mathcal{D}^{{\beta}}v(t)-\phi_2(t, \mathcal{I}^{\alpha} u(t), \mathcal{I}^{\beta} v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = t_k, \ \ k = 1, 2, \dots, q, \\ \Delta v(t_k)-\mathcal{E}_k(v(t_k)) = 0, \ \ \ \ \ \ \Delta v'(t_k)-\mathcal{E}_j^*(v(t_k)) = 0, \ \ k = 1, 2, \dots, q, \\ \nu_3\mathcal{D}^{{\beta}-2}v(t)|_{t = 0} = v_1, \ \ \ \ \ \mu_2 v(t)|_{t = T}+\nu_4 \mathcal{I}^{{\beta}-1} v(t)|_{t = T} = v_2, \end{array} \right.\end{array} \right. \end{align}$

(1.1)

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])

$\begin{equation} ^{RL}D^{{\alpha}}_{a, t}u(t) = \frac{1}{\Gamma(2-{\alpha})}\Big(\frac{d}{dt}\Big)^2\int_a^t\frac{u(s)}{(t-s)^{{\alpha}-1}}ds, \ \ \ t > a, \ \ {\alpha}\in(1, 2), \end{equation}$

(1.2)

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

$\begin{equation} I^{{\beta}}_{a, t}u(t) = \frac{1}{\Gamma({\beta})}\int _a^t(t-s)^{{\beta}-1}u(s)ds, t > a, \end{equation}$

(1.3)

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

$0 = t_0 < t_1 < t_2 < \dots < t_p < t_{p+1} = T\ \ \mbox{and}\ \ \ 0 = {\tau}_0 < {\tau}_1 < {\tau}_2 < \dots < {\tau}_q < {\tau}_{q+1} = T,$

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$

$\begin{align} \left\{ \begin{array}{ll} \left\{ \begin{array}{ll} ^{RL}D^{{\alpha}}_{0, t}u(t)-\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = t_i, \ \ i = 1, 2, \dots, p, \\ \Delta u(t_j)-\mathcal{E}_j(u(t_j)) = 0, \ \ \ \ \ \ \Delta u'(t_j)-\mathcal{E}_j^*(u(t_j)) = 0, \ \ j = 1, 2, \dots, p, \\ - \mathcal{I}^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1}_{0, t} u(t)|_{t = T} = u_2, \end{array} \right.\\ \left\{ \begin{array}{ll} ^{RL}D^{{\beta}}_{0, t}v(t)-\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t}v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = {\tau}_k, \ \ k = 1, 2, \dots, q, \\ \Delta v({\tau}_k)-\mathcal{S}_k(v({\tau}_k)) = 0, \ \ \ \ \ \ \Delta v'({\tau}_k)-\mathcal{S}_k^*(v({\tau}_k)) = 0, \ \ k = 1, 2, \dots, q, \\ \mathcal{I}^{2-{\beta}}_{0, t}v(t)|_{t = 0} = v_1, \ \ \ \ \ \mu_2 v(t)|_{t = T}+\nu_2 \mathcal{I}^{{\beta}-1}_{0, t} v(t)|_{t = T} = v_2, \end{array} \right.\end{array} \right. \end{align}$

(2.1)

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]).

$\begin{equation} \label{{7}}\begin{split} &I^{{\alpha}}_{a, t}I^{{\beta}}_{a, t}u(t) = I^{{\alpha}+{\beta}}_{a, t}u(t), \ \ \ {\alpha}, {\beta} > 0, \\ &I^{{\alpha}}_{a, t}(t-a)^q = \frac{\Gamma(q+1)}{\Gamma(q+{\alpha}+1)}(t-a)^{q+{\alpha}}, \ \ {\alpha} > 0, \ q > -1.\ \ t > a.\end{split} \end{equation}$

(2.2)

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

$\begin{equation} \label{{9}}^{RL}D^{{\alpha}}_{a, t}w(t) = g(t), \ \ \ t\in (a, T], \ \ \ {\alpha}\in(1, 2) \end{equation}$

(2.3)

is given by

$\begin{equation} \label{{400}}\begin{split} w(t)& = c_1(t-a)^{{\alpha}-1}+c_0(t-a)^{{\alpha}-2}+I^{{\alpha}}_{a, t}\ g(t)\\ & = c_1(t-a)^{{\alpha}-1}+c_0(t-a)^{{\alpha}-2}+ \frac{1}{\Gamma({\alpha})}\int_a^t(t-s)^{{\alpha}-1}g(s)ds, \ t\in(a, T], \end{split} \end{equation}$

(2.4)

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

$\begin{equation} \label{{999}}\begin{split}&^{RL}D^{{\alpha}}_{0, t}u(t) = f(t), \ \ \ t\in (0, T], \ \ t\not = t_j, \ j = 1, 2, \dots p, \ \ {\alpha}\in(1, 2), \\ &\Delta u(t_j) = \mathcal{E}_j(u(t_j)), \ \ \ \ \ \ \Delta u'(t_j) = \mathcal{E}_j^*(u(t_j))\ \ j = 1, 2, \dots, p, \\ &I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t}u(t)|_{t = T} = u_2, \end{split} \end{equation}$

(2.5)

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

$u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}f(s)ds, & t\in (0, t_1], \\ c_0t^{{\alpha}-1}+\Big(\sum _{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum _{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum _{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum _{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, & t\in (t_j, t_{j+1}], \\ \end{array} \right.$

where

$\begin{equation} \label{{302}} \begin{split} c_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}f(s)ds\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds+\frac{u_2}{\mu_1}.\end{split} \end{equation}$

(2.6)

Proof. We will use induction.

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

$\begin{equation} \label{{200}} u(t) = c_0t^{{\alpha}-1}+c_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(2.7)

and

$\begin{equation} \label{{209}} u'(t) = c_0({\alpha}-1)t^{{\alpha}-2}+c_1({\alpha}-2)t^{{\alpha}-3}+\frac{{\alpha}-1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(2.8)

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

$I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = c_0\frac{\Gamma({\alpha})}{\Gamma(1)}t|_{t = 0}+c_1\frac{\Gamma({\alpha}-1)}{\Gamma(0)}|_{t = 0}+I^{2}_{0, t}f(t)|_{t = 0} = c_1\Gamma({\alpha}-1),$

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

$\begin{equation} \label{{210}} u(t) = b_0t^{{\alpha}-1}+b_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(2.9)

and

$\begin{equation} \label{{219}} u'(t) = b_0({\alpha}-1)t^{{\alpha}-2}+b_1({\alpha}-2)t^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{0}^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(2.10)

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

$\begin{equation} \label{{300}} \begin{split} u(t_1^+) = &b_0t_1^{{\alpha}-1}+b_1t_1^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{0}^{t_1}(t_1-s)^{{\alpha}-1}f(s)ds\\ &-c_0t_1^{{\alpha}-1}-\frac{u_1}{\Gamma({\alpha}-1)}t_1^{{\alpha}-2}-\frac{1}{\Gamma({\alpha})}\int_0^{t_1}(t_1-s)^{{\alpha}-1}f(s)ds\\ = &(b_0-c_0)t_1^{{\alpha}-1}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})t_1^{{\alpha}-2} = \mathcal{E}_1(u(t_1)). \end{split} \end{equation}$

(2.11)

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

$\begin{equation} \label{{310}} \begin{split} u'(t_1^+) = &b_0({\alpha}-1)t_1^{{\alpha}-2}+b_1({\alpha}-2)t_1^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{0}^{t_1}(t_1-s)^{{\alpha}-2}f(s)ds\\ &-c_0({\alpha}-1)t_1^{{\alpha}-2}-\frac{u_1}{\Gamma({\alpha}-1)}({\alpha}-2)t_1^{{\alpha}-3}-\frac{1}{\Gamma({\alpha}-1)}\int_0^{t_1}(t_1-s)^{{\alpha}-2}f(s)ds\\ = &(b_0-c_0)({\alpha}-1)t_1^{{\alpha}-2}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})({\alpha}-2)t_1^{{\alpha}-3} = \mathcal{E}_1^*(u(t_1)).\end{split} \end{equation}$

(2.12)

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

$\begin{equation} \label{{8}}\begin{split} &(b_0-c_0)({\alpha}-1)t_1^{{\alpha}-2}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})({\alpha}-2)t_1^{{\alpha}-3} = \mathcal{E}_1^*(u(t_1))\\ &(b_0-c_0)t_1^{{\alpha}-1}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})t_1^{{\alpha}-2} = \mathcal{E}_1(u(t_1))\nonumber\end{split} \end{equation}$

$b_0 = c_0+t_1^{2-{\alpha}}\mathcal{E}_1^*(u(t_1))-({\alpha}-2)t_1^{1-{\alpha}}\mathcal{E}_1(u(t_1))$

$b_1 = ({\alpha}-1)t_1^{2-{\alpha}}\mathcal{E}_1(u(t_1))-t_1^{3-{\alpha}}\mathcal{E}_1^*(u(t_1))+\frac{u_1}{\Gamma({\alpha}-1)}.$

Therefore,

$\begin{equation} \label{{220}} \begin{split}u(t)& = c_0t^{{\alpha}-1}+\Big(t_1^{2-{\alpha}}\mathcal{E}_1^*(u(t_1))-({\alpha}-2)t_1^{1-{\alpha}}\mathcal{E}_1(u(t_1))\Big)t^{{\alpha}-1}\\ &\ \ \ \ \ +\Big(({\alpha}-1)t_1^{2-{\alpha}}\mathcal{E}_1(u(t_1))-t_1^{3-{\alpha}}\mathcal{E}_1^*(u(t_1))+\frac{u_1}{\Gamma({\alpha}-1)}\Big)t^{{\alpha}-2}\\ &\ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \ \ \ t\in(t_1, t_2].\end{split} \end{equation}$

(2.13)

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

$\begin{equation} \label{{2209}} \begin{split}u(t)& = c_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))-({\alpha}-2)\sum\limits_{k = 1}^{j-1}t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))\Big)t^{{\alpha}-1}\\ &\ \ \ \ \ +\Big(({\alpha}-1)\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{j-1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))+\frac{u_1}{\Gamma({\alpha}-1)}\Big)t^{{\alpha}-2}\\ &\ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \ \ \ t\in(t_{j-1}, t_{j}].\end{split} \end{equation}$

(2.14)

Denote

$m_0 = c_0+\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))-({\alpha}-2)\sum\limits_{k = 1}^{j-1}t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))$

and

$m_1 = ({\alpha}-1)\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{j-1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))+\frac{u_1}{\Gamma({\alpha}-1)}.$

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

$\begin{equation} \label{{1210}} u(t) = k_0t^{{\alpha}-1}+k_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \end{equation}$

(2.15)

and

$\begin{equation} \label{{1219}} u'(t) = k_0({\alpha}-1)t^{{\alpha}-2}+k_1({\alpha}-2)t^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{0}^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(2.16)

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

$\begin{equation} \label{{18}}\begin{split} &(k_0-m_0)({\alpha}-1)t_j^{{\alpha}-2}+(k_1-m_1)({\alpha}-2)t_j^{{\alpha}-3} = \mathcal{E}_2^*(u(t_2))\\ &(k_0-m_0)t_j^{{\alpha}-1}+(k_1-m_1)t_j^{{\alpha}-2} = \mathcal{E}_j(u(t_j))\nonumber\end{split} \end{equation}$

$k_0 = c_0+\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]$

$k_1 = \frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k)).$

Therefore,

$\begin{equation} \label{{1220}} \begin{split}u(t)& = c_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ &\ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ &\ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \ \ \ t\in(t_j, t_{j+1}], \ j = 1, 2, \dots, p.\end{split} \end{equation}$

(2.17)

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

$\mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = \frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds$

and

$\begin{equation} \label{{888}}\begin{split} &\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T}\\ & = \mu_1 c_0T^{{\alpha}-1}+\mu_1 \Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)T^{{\alpha}-1}\\ &\ \ \ \ \ +\mu_1 \Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{{\alpha}-2}\\ &\ \ \ \ \ +\mu_1 \frac{1}{\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ +\nu_1\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds\\ & = u_2. \end{split} \end{equation}$

(2.18)

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$

$\begin{equation} \label{{99987}}\begin{split}&^{RL}D^{1.5}_{0, t}u(t) = t, \ \ \ t\in (0, 1]\cup (1, 2], \\ &\Delta u(1) = 1, \ \ \ \ \ \ \Delta u'(1) = 0, \\ &I^{0.5}_{0, t}u(t)|_{t = 0} = 0, \ \ \ \ \ u(t)|_{t = 2}+ I^{0.5}_{0, t}u(t)|_{t = 2} = 1.\end{split} \end{equation}$

(2.19)

The solution of (2.19) satisfies the integral equation

$u(t) = \left\{ \begin{array}{ll} c_0t^{0.5}+\frac{t^{-0.5}}{\Gamma(0.5)}+\frac{0.266667 t^{2.5}}{\Gamma(1.5)}, & t\in (0, 1] \\ c_0t^{0.5}+0.5t^{0.5} +\Big(\frac{1}{\Gamma(0.5)}+0.5\Big)t^{-0.5}+\frac{0.266667 t^{2.5}}{\Gamma(1.5)}, & t\in (1, 2], \\ \end{array} \right.$

where

$\begin{equation} \label{{30298}} \begin{split} c_0 & = 0.25-\frac{1.50849}{2^{0.5}\Gamma(1.5)}-\frac{ 1}{2^{0.5}\Gamma(0.5)}\int _0^2(2-s)^{-0.5}u(s)ds.\end{split} \end{equation}$

(2.20)

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

$\begin{equation} \label{{88}}\begin{split}&^{RL}D^{{\alpha}}_{0, t}u(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t)), \ \ \ t\in (0, T], \ \ t\not = t_j, \ j = 1, 2, \dots p, \ \ {\alpha}\in(1, 2)\\ &\Delta u(t_j) = \mathcal{E}_j(u(t_j)), \ \ \ \ \ \ \Delta u'(t_j) = \mathcal{E}_j^*(u(t_j))\ \ j = 1, 2, \dots, p, \\ &I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t}u(t)|_{t = T} = u_2, \end{split} \end{equation}$

(2.21)

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

$\begin{equation} \label{{66}} \begin{split} u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), )ds , \ \ t\in (0, t_1] \\ c_0t^{{\alpha}-1}+\Big(\sum _{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum _{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum _{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum _{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s))ds, \ \ \ \ \ \ t\in (t_j, t_{j+1}], \\ \end{array} \right. \end{split} \end{equation}$

(2.22)

where

$\begin{equation} \label{{87}} \begin{split} c_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds+\frac{u_2}{\mu_1}.\end{split} \end{equation}$

(2.23)

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

$\begin{equation} \label{{44}}\begin{split} u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ t\in (0, t_1], \\ c_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ \ \ t\in (t_j, t_{j+1}], \\ \end{array} \right. \end{split} \end{equation}$

(2.24)

and

$\begin{equation} \label{{45}}\begin{split} v(t) = \left\{ \begin{array}{ll} b_0t^{{\alpha}-1}+\frac{v_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}\phi_2(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ t\in (0, t_1] \\ b_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{S}_k^*(v(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{S}_k(v(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{S}_k(v(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{S}_k^*(v(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}\phi_2(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ \ \ t\in (t_j, t_{j+1}], \\ \end{array} \right. \end{split} \end{equation}$

(2.25)

where

$\begin{equation} \label{{5302}} \begin{split} c_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds+\frac{u_2}{\mu_1}\\ b_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{S}_k^*(v(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{S}_k(v(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{v_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{S}_k(v(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{S}_k^*(v(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_2(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ -\frac{\nu_2}{\mu_2T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}v(s)ds+\frac{v_2}{\mu_2}.\end{split} \end{equation}$

(2.26)

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,

$\begin{align} \left\{ \begin{array}{ll} \left\{ \begin{array}{ll} ^{RL}D^{{\alpha}}_{t_i, t}u(t)-\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)) = 0\ \ \ \mbox{for}\ t\in (t_i, t_{i+1}]\ \ i = 0, 1, 2, \dots, p, \\ I^{2-{\alpha}}_{t_j, t}u'(t)|_{t = t_j} = \mathcal{P}_ju(t_j)+\mathcal{Q}_j, \ \ j = 1, 2, \dots, p, \\ I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t} u(t)|_{t = T} = u_2, \end{array} \right.\\ \left\{ \begin{array}{ll} ^{RL}D^{{\beta}}_{{\tau}_k, t}v(t)-\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t}v(t)) = 0\ \ \ \mbox{for}\ t\in ({\tau}_k, {\tau}_{k+1}]\ \ k = 0, 1, 2, \dots, q, \\ I^{2-{\alpha}}_{{\tau}_k, t}u'(t)|_{t = {\tau}_k} = \mathcal{P}_k^*v({\tau}_k)+\mathcal{Q}_jk^*, \ \ k = 1, 2, \dots, q, \\ I^{2-{\beta}}_{0, t}v(t)|_{t = 0} = v_1, \ \ \ \ \ \mu_2 v(t)|_{t = T}+\nu_2 I^{{\beta}-1}_{0, t} v(t)|_{t = T} = v_2, \end{array} \right.\end{array} \right. \end{align}$

(3.1)

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

$\begin{equation} \label{{1999}}\begin{split}&^{RL}D^{{\alpha}}_{t_i, t}u(t) = f(t), \ \ \ \mbox{for}\ t\in (t_i, t_{i+1}], \ \ i = 0, 1, 2, \dots, p, \ \ {\alpha}\in(1, 2)\\ &I^{2-{\alpha}}_{t_j, t}u'(t)|_{t = t_j} = \mathcal{P}_ju(t_j)+\mathcal{Q}_j, \ \ j = 1, 2, \dots, p, \\ &I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t}u(t)|_{t = T} = u_2, \end{split} \end{equation}$

(3.2)

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

$u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+I^{{\alpha}}_{0, t}f(t), \ \ \ t\in (0, t_1] \\ \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)(t-t_m)^{{\alpha}-1}\prod\limits_{k = 1}^m\frac{\mathcal{P}_k(t_k-t_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ \ \ \ \ +I^{{\alpha}}_{t_m, t}f(t)+(t-t_m)^{{\alpha}-1}\sum\limits_{k = 1}^m\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}f(t)|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}, \\ \ \ \ \ \ \mathit{\mbox{for}}\ \ \ \ t\in (t_m, t_{m+1}], m = 1, 2, \dots, p, \end{array} \right.$

where

$\begin{equation} \label{{1302}}\begin{split} c_0& = \frac{u_2}{\mu_1}(T-t_p)^{1-{\alpha}}\mathcal{M} -\frac{u_1}{t_1\Gamma({\alpha}-1)} -(T-t_p)^{1-{\alpha}}\mathcal{M} I^{{\alpha}}_{t_p, t}f(t)|_{t = T}\\ &-\mathcal{M}\sum\limits_{k = 1}^p\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}f(t)|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1} (T-t_p)^{1-{\alpha}}\mathcal{M}\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds, \\ \mathcal{M}& = ({\alpha}-1)^p\Gamma^p({\alpha}-1)\prod\limits_{k = 1}^p\mathcal{P}_k(t_k-t_{k-1})^{1-{\alpha}}. \end{split} \end{equation}$

(3.3)

Proof. We will use induction.

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

$\begin{equation} \label{{1200}} u(t) = c_0t^{{\alpha}-1}+c_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(3.4)

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

$\begin{equation} \label{{12105}} u(t) = b_0(t-t_1)^{{\alpha}-1}+b_1(t-t_1)^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{t_1}^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(3.5)

and

$\begin{equation} \label{{12195}} u'(t) = b_0({\alpha}-1)(t-t_1)^{{\alpha}-2}+b_1({\alpha}-2)(t-t_1)^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{t_1}^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(3.6)

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

$\begin{equation} \label{{1309}} \begin{split}& I^{2-{\alpha}}_{t_1, t}u'(t) = b_0({\alpha}-1)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-2}+b_1({\alpha}-2)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3} + I^{2-{\alpha}}_{t_1, t} I^{{\alpha}-1}_{t_1, t}f(t)\\ & = b_0({\alpha}-1)\Gamma({\alpha} -1)+b_1({\alpha}-2)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3} + I^{1}_{t_1, t}f(t)\end{split} \end{equation}$

(3.7)

and

$\begin{equation} \label{{6309}} \begin{split}& I^{2-{\alpha}}_{t_1, t}u'(t)|_{t = t_1} = b_0({\alpha}-1)\Gamma({\alpha} -1)+b_1({\alpha}-2)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3}|_{t = t_1} = \mathcal{P}_1u(t_1)+\mathcal{Q}_1 < \infty.\end{split} \end{equation}$

(3.8)

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

$b_0 = \frac{\mathcal{P}_1}{({\alpha}-1)\Gamma({\alpha}-1)}(c_0(t_1-t_0)^{{\alpha}-1}+\frac{u_1}{t_1\Gamma({\alpha}-1)}(t_1-t_0)^{{\alpha}-1}+I^{{\alpha}}_{t_0, t}f(t)|_{t = t_1})+\frac{\mathcal{Q}_1}{({\alpha}-1)\Gamma({\alpha}-1)}.$

Thus,

$\begin{equation} \label{{12196}} \begin{split}u(t)& = \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)\frac{\mathcal{P}_1(t_1-t_0)^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}(t-t_1)^{{\alpha}-1}\\ &\ \ \ +\frac{\mathcal{P}_1}{({\alpha}-1)\Gamma({\alpha}-1)}I^{{\alpha}}_{t_0, t}f(t)|_{t = t_1}(t-t_1)^{{\alpha}-1}\\ &\ \ \ +\frac{\mathcal{Q}_1}{({\alpha}-1)\Gamma({\alpha}-1)}(t-t_1)^{{\alpha}-1} +I^{{\alpha}}_{t_1, t}f(t).\end{split} \end{equation}$

(3.9)

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

$\begin{equation} \label{{12106}} u(t) = k_0(t-t_j)^{{\alpha}-1}+k_1(t-t_j)^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{t_j}^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds \end{equation}$

(3.10)

and

$\begin{equation} \label{{12198}} u'(t) = k_0({\alpha}-1)(t-t_j)^{{\alpha}-2}+k_1({\alpha}-2)(t-t_j)^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{t_j}^t(t-s)^{{\alpha}-2}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds. \end{equation}$

(3.11)

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

$k_0 = \frac{\mathcal{P}_j}{({\alpha}-1)\Gamma({\alpha}-1)}u(t_j)+\frac{\mathcal{Q}_j}{({\alpha}-1)\Gamma({\alpha}-1)}.$

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

$\mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = \frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds$

and

$\begin{equation} \label{{6888}}\begin{split} &\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T}\\ & = \mu_1 \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)(T-t_p)^{{\alpha}-1}\prod\limits_{k = 1}^p\frac{\mathcal{P}_k(t_k-t_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ +\mu_1 I^{{\alpha}}_{t_p, t}f(t)|_{t = T}+\mu_1(T-t_p)^{{\alpha}-1}\sum\limits_{k = 1}^p\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}f(t)|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ +\nu_1\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds\\ & = u_2 \end{split} \end{equation}$

(3.12)

and we obtain (3.3).

Example 2. Consider the equation

$\begin{equation} \label{{99988}}\begin{split}&^{RL}D^{1.5}_{0, t}u(t) = t, \ \ \ t\in (0, 1], \ \ \ \ \ \ \ ^{RL}D^{1.5}_{1, t}u(t) = t, \ t\in(1, 2], \\ &I^{0.5}_{1, t}u'(t)|_{t = 1} = 1, \\ &I^{0.5}_{0, t}u(t)|_{t = 0} = 0, \ \ \ \ \ u(t)|_{t = 2}+ I^{0.5}_{0, t}u(t)|_{t = 2} = 1.\end{split} \end{equation}$

(3.13)

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

$u(t) = \left\{ \begin{array}{ll} \frac{0.266667 t^{2.5}}{\Gamma(1.5)}, & t\in (0, 1], \\ \frac{1}{\Gamma(1.5)}(-0.4 (t-1)^{2.5} + \frac{2 (t-1)^{1.5} t}{3})+\frac{(t-1)^{0.5}}{(0.5)\Gamma(0.5)}, & t\in (1, 2]. \end{array} \right.$

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

$u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+I^{{\alpha}}_{0, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)), \ \ \ t\in (0, t_1] \\ \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)(t-t_m)^{{\alpha}-1}\prod _{k = 1}^m\frac{\mathcal{P}_k(t_k-t_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ \ \ \ \ +I^{{\alpha}}_{t_m, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))\\ \ \ \ \ +(t-t_m)^{{\alpha}-1}\sum _{k = 1}^m\frac{\mathcal{P}_{k}\ I^{{\alpha}}_{t_{k-1}, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod _{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}, \\ \ \ \ \ \ \mathit{\mbox{for}}\ \ \ \ t\in (t_m, t_{m+1}], m = 1, 2, \dots, p, \end{array} \right.$

and

$v(t) = \left\{ \begin{array}{ll} b_0t^{{\alpha}-1}+\frac{v_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+I^{{\alpha}}_{0, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)), \ \ \ t\in (0, {\tau}_1] \\ \Big(b_0+\frac{v_1}{{\tau}_1\Gamma({\alpha}-1)}\Big)(t-{\tau}_m)^{{\alpha}-1}\prod _{k = 1}^m\frac{\mathcal{P}_k^*({\tau}_k-{\tau}_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ \ \ \ \ +I^{{\alpha}}_{{\tau}_m, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))\\ \ \ \ \ +(t-{\tau}_m)^{{\alpha}-1}\sum _{k = 1}^m\frac{\mathcal{P}_{k}^*\ I^{{\alpha}}_{{\tau}_{k-1}, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = {\tau}_{k}}+\mathcal{Q}_{k}^*}{({\alpha}-1)\Gamma({\alpha}-1)}\prod _{j = k+1}^{m}\frac{\mathcal{P}_j^*({\tau}_j-{\tau}_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}, \\ \ \ \ \ \ \mathit{\mbox{for}}\ \ \ \ t\in ({\tau}_m, {\tau}_{m+1}], m = 1, 2, \dots, q, \end{array} \right.$

where

$\begin{equation} \label{{53029}} \begin{split} c_0& = \frac{u_2}{\mu_1}(T-t_p)^{1-{\alpha}}\mathcal{M} -\frac{u_1}{t_1\Gamma({\alpha}-1)} \\ &\ \ \ -(T-t_p)^{1-{\alpha}}\mathcal{M} I^{{\alpha}}_{t_q, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = T}\\ &\ \ \ -\mathcal{M}\sum\limits_{k = 1}^p\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{p}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1} (T-t_p)^{1-{\alpha}}\mathcal{M}\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds, \\ \mathcal{M}& = ({\alpha}-1)^p\Gamma^p({\alpha}-1)\prod\limits_{k = 1}^p\mathcal{P}_k(t_k-t_{k-1})^{1-{\alpha}}, \nonumber\end{split} \end{equation}$

$\begin{equation} \label{{53022}} \begin{split}b_0& = \frac{v_2}{\mu_2}(T-{\tau}_q)^{1-{\alpha}}\mathcal{C} -\frac{v_1}{{\tau}_1\Gamma({\alpha}-1)} \\ &\ \ \ -(T-{\tau}_q)^{1-{\alpha}}\mathcal{C} \ I^{{\alpha}}_{{\tau}_q, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))_{t = T}\\ &\ \ \ -\mathcal{C}\sum\limits_{k = 1}^q\frac{\mathcal{P}_{k}^*\ I^{{\alpha}}_{{\tau}_{k-1}, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = {\tau}_{k}}+\mathcal{Q}_{k}^*}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{q}\frac{\mathcal{P}_j^*({\tau}_j-{\tau}_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ -\frac{\nu_2}{\mu_2} (T-t_q)^{1-{\alpha}}\mathcal{C}\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds, \\ \mathcal{C}& = ({\alpha}-1)^q\Gamma^q({\alpha}-1)\prod\limits_{k = 1}^q\mathcal{P}_k^*({\tau}_k-{\tau}_{k-1})^{1-{\alpha}}.\nonumber\end{split} \end{equation}$

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.

References

[1]	B. Jiao, Z. Lian, X. Gu, A dynamic inertia weight particle swarm optimization algorithm, Chaos Soliton Fract, 37 (2008), 698–705. https://doi.org/10.1016/j.chaos.2006.09.063 doi: 10.1016/j.chaos.2006.09.063
[2]	D. E. Goldberg, Genetic algorithm in search optimization and machine learning, Addison Wesley, 8 (1989), 2104–2116. https://dl.acm.org/doi/book/10.5555/534133 doi: 10.5555/534133
[3]	S. Kirpatrick, C. D. Gelatt, M. P. Vecchi, Optimization by simulated annealing, Readings Computer Vision, 220 (1983), 671–680. https://doi.org/10.1126/science.220.4598.671 doi: 10.1126/science.220.4598.671
[4]	S. Mirjalili, SCA: A Sine Cosine Algorithm for solving optimization problems, Knowl-based Syst., 96 (2016), 120–133. https://doi.org/10.1016/j.knosys.2015.12.022 doi: 10.1016/j.knosys.2015.12.022
[5]	Z. Zhang, J. Lin, Y. Shi, Application of artificial bee colony algorithm to maximum likelihood DOA estimation, J. Bionic. Eng., 10 (2013), 100–109. https://doi.org/10.1016/S1672-6529(13)60204-8 doi: 10.1016/S1672-6529(13)60204-8
[6]	S. Feng, Z. Zhang, Y. Shi, Introduction of bat algorithm into maximum likelihood DOA estimation, Modern Electronics Technique, 39 (2016), 26–29. https://doi.org/10.16652/j.issn.1004-373x.2016.08.007 doi: 10.16652/j.issn.1004-373x.2016.08.007
[7]	X. Fan, L. Pang, P. Shi, G. Li, X. Zhang, Application of bee evolutionary genetic algorithm to maximum likelihood direction-of-arrival estimation, Math. Probl. Eng., 2019 (2019), 1–11. https://doi.org/10.1155/2019/6035870 doi: 10.1155/2019/6035870
[8]	M. Jain, V. Singh, A. Rani, A novel nature-inspired algorithm for optimization: Squirrel search algorithm, Swarm Evol. Comput., 44 (2018), 148–175. https://doi.org/10.1016/j.swevo.2018.02.013 doi: 10.1016/j.swevo.2018.02.013
[9]	W. Zhao, L. Wang, Z. Zhang, Atom search optimization and its application to solve a hydrogeologic parameter estimation problem, Knowl-based Syst., 163 (2018), 283–304. https://doi.org/10.1016/j.knosys.2018.08.030 doi: 10.1016/j.knosys.2018.08.030
[10]	H. C. Corben, P. Stehle, Classical Mechanics, Physics Today, 6 (1953). https://doi.org/10.1063/1.3061288 doi: 10.1063/1.3061288
[11]	J. P. Ryckaert, G. Ciccotti, H. J. C Berendsen, Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes, J. Comput. Phys., 23 (1977), 327–341. https://doi.org/10.1016/0021-9991(77)90098-5 doi: 10.1016/0021-9991(77)90098-5
[12]	A. Stone, The theory of intermolecular forces, Pure. Appl. Chem., 51 (1979), 1627–1636. https://doi.org/10.1351/pac197951081627 doi: 10.1351/pac197951081627
[13]	J. E. Jones, On the determination of molecular fields Ⅱ. From the equation of state of a gas, P. Roy. Soc. A-Math. Phy., 106 (1924), 463–477. https://doi.org/10.2307/94265 doi: 10.2307/94265
[14]	W. Zhao, L. Wang, Z. Zhang, A novel atom search optimization for dispersion coefficient estimation in groundwater, Future Gener Comp. Sy., 91 (2018), 601–610. https://doi.org/10.1016/j.future.2018.05.037 doi: 10.1016/j.future.2018.05.037
[15]	A. M. Agwa, A. A. El-Fergany, G. M. Sarhan, Steady-State modeling of fuel cells based on atom search optimizer, Energies, 12 (2019), 1884. https://doi.org/10.3390/en12101884 doi: 10.3390/en12101884
[16]	A. Almagboul Mohammed, F. Shu, Y. Qian, X. Zhou, J. Wang, J. Hu, Atom search optimization algorithm based hybrid antenna array receive beamforming to control sidelobe level and steering the null, Aeu-int J. Electron. C., 111 (2019), 152854. https://doi.org/10.1016/j.aeue.2019.152854 doi: 10.1016/j.aeue.2019.152854
[17]	S. Barshandeh, A new hybrid chaotic atom search optimization based on tree-seed algorithm and Levy flight for solving optimization problems, Eng. Comput-germany, 37 (2021), 3079–3122. https://doi.org/10.1007/s00366-020-00994-0 doi: 10.1007/s00366-020-00994-0
[18]	K. K. Ghosh, R. Guha, S. Ghosh, S. K. Bera, R. Sarkar, Atom Search Optimization with Simulated Annealing-a Hybrid Metaheuristic Approach for Feature Selection, arXiv preprint arXiv: 2005.08642, (2020). https://arXiv.org/pdf/2005.08642v1
[19]	M. A. Elaziz, N. Nabil, A. A. Ewees, S. Lu, Automatic data clustering based on hybrid atom search optimization and Sine-Cosine algorithm, 2019 IEEE Congress on Evolutionary Computation (CEC), (2019), 2315–2322. https://doi.org/10.1109/CEC.2019.8790361 doi: 10.1109/CEC.2019.8790361
[20]	P. Sun, H. Liu, Y. Zhang, L. Tu, Q. Meng, An intensify atom search optimization for engineering design problems, Appl. Math. Model., 89 (2021), 837–859. https://doi.org/10.1016/j.apm.2020.07.052 doi: 10.1016/j.apm.2020.07.052
[21]	L. Xu, J. Chen, Y. Gao, Off-Grid DOA estimation based on sparse representation and rife algorithm, Microelectron J., 59 (2017), 193–201. https://doi.org/10.2528/PIERM17070404 doi: 10.2528/PIERM17070404
[22]	A. Peyman, Z. Kordrostami, K. Hassanli, Design of a MEMS bionic vector hydrophone with piezo-gated MOSFET readout, Prog. Electromagn Res. M., 98 (2020), 104748. https://doi.org/10.1016/j.mejo.2020.104748 doi: 10.1016/j.mejo.2020.104748
[23]	H. Song, M. Diao, T. Tang, J. Qin, Vector-Sensor Array DOA Estimation Based on Spatial Time-Frequency Distribution, 2018 Eighth International Conference on Instrumentation & Measurement, Computer, Communication and Control (IMCCC), (2020), 1351–1356. https://doi.org/10.1109/IMCCC.2018.00280 doi: 10.1109/IMCCC.2018.00280
[24]	M. Cao, X. Mao, L. Huang, Elevation, azimuth, and polarization estimation with nested electromagnetic vector-sensor arrays via tensor modeling, Eurasip J. Wirel. Comm., 2020 (2020), 153. https://doi.org/10.1186/s13638-020-01764-8 doi: 10.1186/s13638-020-01764-8
[25]	V. Baron, A. Finez, S. Bouley, F. Fayet, J. I. Mars, B. Nicolas, Hydrophone array optimization, conception, and validation for localization of acoustic sources in deep-Sea mining, IEEE J. Oceanic. Eng., 46 (2021), 555–563. https://doi.org/10.1109/JOE.2020.3004018 doi: 10.1109/JOE.2020.3004018
[26]	W. Wand, Q. Zhang, W. Shi, J. Shi, X. Wang, Iterative sparse covariance matrix fitting direction of arrival estimation method based on vector hydrophone array, Xibei Gongye Daxue Xuebao, 38 (2020), 14–23. https://doi.org/10.1051/jnwpu/20203810014 doi: 10.1051/jnwpu/20203810014
[27]	K. Aghababaiyan, R. G.Zefreh, V. Shah-Mansouri, 3D-OMP and 3D-FOMP algorithms for DOA estimation, Phys. Commun-amst, 31 (2018), 87–95. https://doi.org/10.1016/j.phycom.2018.10.005 doi: 10.1016/j.phycom.2018.10.005
[28]	K. Aghababaiyan, V. Shah-Mansouri, B. Maham, High-precision OMP-based direction of arrival estimation scheme for hybrid non-uniform array, IEEE Commun. Lett., 24 (2019), 354–357. https://doi.org/10.1109/LCOMM.2019.2952595 doi: 10.1109/LCOMM.2019.2952595
[29]	A. Nehorai, E. Paldi, , Acoustic vector-sensor array processing, IEEE T. Signal. Proces., 42 (1994), 2481–2491. https://doi.org/10.1109/ACSSC.1992.269285 doi: 10.1109/ACSSC.1992.269285
[30]	K. T. Wong, M. D. Zoltowski, Root-MUSIC-based azimuth-elevation angle-of-arrival estimation with uniformly spaced but arbitrarily oriented velocity hydrophones, IEEE T. Signal. Proces., 47 (1999), 3250–3260. https://doi.org/10.1109/78.806070 doi: 10.1109/78.806070
[31]	K. T. Wong, M. D. Zoltowski, Uni-vector-sensor ESPRIT for multisource azimuth, elevation, and polarization estimatio, IEEE T. Antenn. Propag., 45 (1997), 1467–1474. https://doi.org/10.1109/8.633852 doi: 10.1109/8.633852
[32]	I. Ziskind, M. Wax, Maximum likelihood localization of multiple sources by alternating projection, IEEE Trans. Acoust. Speech Signal Process, 36 (1988), 1553–1560. https://doi.org/10.1109/29.7543 doi: 10.1109/29.7543
[33]	M. Feder, E. Weinstein, Parameter estimation of superimposed signals using the EM algorithm, IEEE Trans. Acoust. Speech Signal Process, 36 (1988), 477–489. https://doi.org/10.1109/29.1552 doi: 10.1109/29.1552
[34]	Y. Zheng, L. Liu, X. Yang, SPICE-ML Algorithm for Direction-of-Arrival Estimation, Sensors, 20 (2019), 119. https://doi.org/10.3390/s20010119 doi: 10.3390/s20010119
[35]	Y. Hu, J. Lu, X. Qiu, Direction of arrival estimation of multiple acoustic sources using a maximum likelihood method in the spherical harmonic domain, Appl. Acoust., 135 (2018), 85–90. https://doi.org/10.1016/j.apacoust.2018.02.005 doi: 10.1016/j.apacoust.2018.02.005
[36]	J. W. Paik, K. H. Lee, J. H. Lee, Asymptotic performance analysis of maximum likelihood algorithm for direction-of-arrival estimation: Explicit expression of estimation error and mean square error, Applied Sciences, 10 (2020), 2415. https://doi.org/10.3390/app10072415 doi: 10.3390/app10072415
[37]	S. Jesus, Efficient ML direction of arrival estimation assuming unknown sensor noise powers, arXiv preprint arXiv: 2001.01935, (2020), https: //arXiv: 2001.01935
[38]	Y. Yoon, Y. H. Kim, Optimizing taxon addition order and branch lengths in the construction of phylogenetic trees using maximum likelihood, J. Bioinf. Comput. Biol., 18 (2020), 837–859. https://doi.org/10.1142/S0219720020500031 doi: 10.1142/S0219720020500031
[39]	P. Vishnu, C. S. Ramalingam, An improved LSF-based algorithm for sinusoidal frequency estimation that achieves maximum likelihood performance, 2020 International Conference on Signal Processing and Communications (SPCOM), (2020), 1–5. https://doi.org/10.1109/SPCOM50965.2020.9179546
[40]	M. Li, Y. Lu, Genetic algorithm based maximum likelihood DOA estimation, RADAR 2002, 2002 (2002), 502–506. https://doi.org/10.1109/RADAR.2002.1174766 doi: 10.1109/RADAR.2002.1174766
[41]	A. Sharma, S. Mathur, Comparative analysis of ML-PSO DOA estimation with conventional techniques in varied multipath channel environment, Wireless Pers. Commun., 100 (2018), 803–817. https://doi.org/10.1007/s11277-018-5350-0 doi: 10.1007/s11277-018-5350-0
[42]	P. Wang, Y. Kong, X. He, M. Zhang, X. Tan, An improved squirrel search algorithm for maximum likelihood DOA estimation and application for MEMS vector hydrophone array, IEEE Access, 7 (2019), 118343–118358. https://doi.org/10.1109/ACCESS.2019.2936823 doi: 10.1109/ACCESS.2019.2936823
[43]	L. Cai, H. Tian, H. Chen, J. Hu, A random maximum likelihood algorithm based on limited PSO initial space, Computer Modernization, 282 (2019), 60–65. https://doi.org/10.3969/j.issn.1006-2475.2019.02.011 doi: 10.3969/j.issn.1006-2475.2019.02.011

This article has been cited by:

1.	Vitor Miguel Ribeiro, Pioneering paradigms: unraveling niche opportunities in green finance through bibliometric analysis of nation brands and brand culture, 2024, 6, 2643-1092, 287, 10.3934/GF.2024012
2.	Manuel Mojica, Pedro R. Palos-Sanchez, Edgar Cabanas, Is there innovation management of emotions or just the commodification of happiness? A sentiment analysis of happiness apps, 2024, 1460-1060, 10.1108/EJIM-11-2023-0963
3.	Pedro R. Palos-Sanchez, Raúl J. Chang-Tam, José A. Folgado-Fernández, The Role of Neobanks and FinTech in Sustainable Finance and Technology. The Customer/User Perspective for Entrepreneurs, 2025, 27730328, 100109, 10.1016/j.stae.2025.100109

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)