Research article

Existence, continuous dependence and finite time stability for Riemann-Liouville fractional differential equations with a constant delay

  • Received: 24 November 2019 Accepted: 31 March 2020 Published: 21 April 2020
  • MSC : 34A08, 34A37

  • Non-linear scalar Riemann-Liouville fractional differential equation with a constant delay is studied on a finite interval. An initial value problem is set up in appropriate way combining the idea of the initial time interval in ordinary differential equations with delays and the properties of Riemann-Liouville fractional derivatives. The mild solution of the studied initial value problem is defined. The existence, uniqueness, continuous dependence on the initial functions, finite time stability of the mild solutions are investigated.

    Citation: Snezhana Hristova, Antonia Dobreva. Existence, continuous dependence and finite time stability for Riemann-Liouville fractional differential equations with a constant delay[J]. AIMS Mathematics, 2020, 5(4): 3809-3824. doi: 10.3934/math.2020247

    Related Papers:

    [1] Yanting Xiao, Yifan Shi . Robust estimation for varying-coefficient partially nonlinear model with nonignorable missing response. AIMS Mathematics, 2023, 8(12): 29849-29871. doi: 10.3934/math.20231526
    [2] Heng Liu, Xia Cui . Adaptive estimation for spatially varying coefficient models. AIMS Mathematics, 2023, 8(6): 13923-13942. doi: 10.3934/math.2023713
    [3] Yanping Liu, Juliang Yin . B-spline estimation in varying coefficient models with correlated errors. AIMS Mathematics, 2022, 7(3): 3509-3523. doi: 10.3934/math.2022195
    [4] Anum Iftikhar, Hongbo Shi, Saddam Hussain, Ather Qayyum, M. El-Morshedy, Sanaa Al-Marzouki . Estimation of finite population mean in presence of maximum and minimum values under systematic sampling scheme. AIMS Mathematics, 2022, 7(6): 9825-9834. doi: 10.3934/math.2022547
    [5] Sanaa Al-Marzouki, Christophe Chesneau, Sohail Akhtar, Jamal Abdul Nasir, Sohaib Ahmad, Sardar Hussain, Farrukh Jamal, Mohammed Elgarhy, M. El-Morshedy . Estimation of finite population mean under PPS in presence of maximum and minimum values. AIMS Mathematics, 2021, 6(5): 5397-5409. doi: 10.3934/math.2021318
    [6] Jieqiong Lu, Peixin Zhao, Xiaoshuang Zhou . Orthogonality based modal empirical likelihood inferences for partially nonlinear models. AIMS Mathematics, 2024, 9(7): 18117-18133. doi: 10.3934/math.2024884
    [7] Zawar Hussain, Atif Akbar, Mohammed M. A. Almazah, A. Y. Al-Rezami, Fuad S. Al-Duais . Diagnostic power of some graphical methods in geometric regression model addressing cervical cancer data. AIMS Mathematics, 2024, 9(2): 4057-4075. doi: 10.3934/math.2024198
    [8] Peng Lai, Wenxin Tian, Yanqiu Zhou . Semi-supervised estimation for the varying coefficient regression model. AIMS Mathematics, 2024, 9(1): 55-72. doi: 10.3934/math.2024004
    [9] Dayang Dai, Dabuxilatu Wang . A generalized Liu-type estimator for logistic partial linear regression model with multicollinearity. AIMS Mathematics, 2023, 8(5): 11851-11874. doi: 10.3934/math.2023600
    [10] Abdullah Mohammed Alomair, Weineng Zhu, Usman Shahzad, Fawaz Khaled Alarfaj . Non-parametric calibration estimation of distribution function under stratified random sampling. AIMS Mathematics, 2025, 10(2): 4457-4472. doi: 10.3934/math.2025205
  • Non-linear scalar Riemann-Liouville fractional differential equation with a constant delay is studied on a finite interval. An initial value problem is set up in appropriate way combining the idea of the initial time interval in ordinary differential equations with delays and the properties of Riemann-Liouville fractional derivatives. The mild solution of the studied initial value problem is defined. The existence, uniqueness, continuous dependence on the initial functions, finite time stability of the mild solutions are investigated.


    Recently, neural networks such as the Hopfield neural network, cellular neural network, Cohen-Grossberg neural network and bidirectional associative neural network have attracted much attention due to their significant promise for many practical applications[1,2,3,4,5,6]. For example, in the field of signal processing, based on the theoretical results of discrete-time high-order switched neural networks synchronization, the encryption and decryption scheme of multi-channel audio signal design has good security [3]. In the field of fault diagnosis, by using deep parameter-free reconstruction-classification networks with parameter-free adaptively rectified linear units, the fault characteristics of vibration signals under the same fault state under different operating conditions can be better captured [5]. In the field of image processing, the theoretical results of random synchronization control framework based on semi-Markov switching quaternion-valued neural networks can be effectively applied to image encryption [6]. In particular, stability properties of the neural networks play a significant role in their design for solving practical problems. There are some different factors to influence the stability of neural networks. In particular, time delay is often unavoidable. Moreover, for the neural networks with time delay, there exists many techniques to reduce the conservatism of stability conditions, such as the improved bounding technique [7], free-weighting matrix theory [8], integral inequality technique [9] and so on.

    It is well known that many dynamical systems may switch in different unpredictable modes, such as random failures [10]. As an effective tool, the Markovian jump process can be used to model these switching systems. Until now, there existed many literatures to investigate the stability, stabilization and observation for Markovian jump systems. For example, the authors investigated the adaptive sliding mode control problem of nonlinear Markovian jump systems with partly unknown transition probabilities, and obtained some conditions to guarantee the stochastic stability of the closed-loop system in [11]. The authors in [12] discussed the realization of H finite-time control for a class of uncertain stochastic time-delay systems with unmeasured states through sliding mode control, and obtained some conditions to guarantee the system state was stabilized within a limited time interval. In [13], the authors established some passivity analysis criteria for Markovian jump singularly perturbed systems with partially unknown probabilities by using the hidden Markov model, and presented a unified controller design method to ensure the passivity of the system. Much more literatures can be found in [14,15,16].

    In recent years, sliding mode control has become an effective robust control method because it is insensitive to model uncertainties, parameter variations and external disturbances. So, the sliding mode control has been used for lots of physical systems, such as robot manipulators, automotive engines and power systems. There are some existing results for linear or nonlinear systems[17,18,19,20,21,22]. For example, the authors in [17] investigated the adaptive sliding mode control issue for switched nonlinear systems with matched and mismatched uncertainties, designed the switched adaptive sliding mode control law and estimated the upper bound parameters of the matched uncertainty. The authors in [18] considered the synchronization of delayed chaotic neural networks with unknown disturbance via the observer-based sliding mode control, where the sliding surface involves an integral structure and a discontinuous controller.

    In addition, synchronization has received considerable attention from various research fields, such as secure communication [3], engineering [23] and chemistry [24]. Synchronization phenomenons can be observed in many real systems such as neural systems, lasers and electronic circuits [25,26,27,28]. For example, the authors in [25] dealt with chaos synchronization for master-slave piecewise linear systems and provided some new sufficient conditions by using a Lyapunov approach and the so-called S-procedure. In [26], the authors designed a proportional-derivative (PD) controller for the master-slave synchronization of chaotic Lurie systems and presented a new synchronization criterion based on Lyapunov functions with a quadratic form of states and nonlinear functions of the systems. In [27], the authors studied the finite-time lag synchronization issue of master-slave complex networks with unknown signal propagation delays by the linear and adaptive error state feedback approaches and discovered that the setting time was related to initial values.

    Motivated by the above discussion, we will consider the master-slave synchronization for uncertain neural networks with time-delay by using the sliding mode control method. The main contributions of our paper are as follows: (1) The master-slave synchronization for uncertain Markov jump neural networks with time-delay is considered by using the sliding mode control method, (2) the effect of time-delay and uncertainty for the neural networks stability are considered, and the uncertain parts in the neural networks only need be bounded other than any structure condition and (3) the provided sliding mode controller is very general.

    The rest of this letter is organized as follows. In Section 2, the considered neural network model and some preliminaries are given. In Section 3, some sufficient conditions are proposed by the sliding mode control. In Section 4, a numerical example is provided to illustrate the effectiveness of the method proposed. In the last section, conclusions are presented.

    Notation:

    Rn: the n-dimensional Euclidean space; Rn×m: the set of all n×m real matrices; ||||: the two-norm of a vector; ||||1: its one-norm; In: the n order unit matrix; λmin(H): the minimum eigenvalues of matrix H; the notation X>Y, where X,Y are symmetric matrices, meaning that XY is a positive definite symmetric matrix. For a given matrix ARn×n, AT denotes its transpose. sign() is the sign function. in a symmetric matrix denotes the symmetric terms.

    Let (Ω,F,Ft,P) be a probability space related to an increasing family {Ft}t0 of the σ-algebras FtF, where Ω is the sample space. F is the σ-algebras of the sample space and P is the probability measure defined on F.

    Consider the uncertain time delay neural networks with the Markovian jump defined on the probability space (Ω,F,Ft,P) as follows:

    M:{˙xm(t)=A(r(t))xm(t)+[B(r(t))+ΔB(r(t))]f(xm(t))+[C(r(t))+ΔC(r(t))]g(xm(td(t)))+J,xm(t)=φ(t),t[d,0], (2.1)

    where xm(t)=(xm1(t),xm2(t),...,xmn(t))TRn is the state of neuron networks and A(r(t))Rn×n and B(r(t))Rn×m and C(r(t))Rn×m are coefficient matrices. ΔB(r(t)) and ΔC(r(t))Rn×m denote the system's uncertain parts and satisfy

    [ΔB(r(t))ΔC(r(t))]=M(r(t))W(t)[N1(r(t))N2(r(t))],

    where matrix W(t) satisfies WT(t)W(t)I, M(r(t)),N1(r(t)),N2(r(t)) are some known matrices with appropriate dimensions. d(t) represents time delay and satisfies 0d(t)d and ˙d(t)μ1. f(xm(t)) and g(xm(td(t))) are the neuron activation functions, and JRn is a constant vector. {r(t),t0} is a finite state Markov jumping process and represents the switching process among different modes, which takes values in a state space L={1,2,...,l}, and l is the number of modes. Let =[πij]l×l denote the transition rate matrix, where the mode transition probabilities are

    Pr{r(t+Δt)=j|r(t)=i}={πijΔt+o(Δt),ij,1+πiiΔt+o(Δt),i=j, (2.2)

    where Δt>0 and limΔt0o(Δt)Δt=0, πij satisfies πij>0 with ij and πii=lj=1,jiπij for each mode i.

    In order to be notional convenience, for the i-th mode, system (2.1) can be rewritten as

    {˙xm(t)=Aixm(t)+(Bi+ΔBi)f(xm(t))+(Ci+ΔCi)g(xm(td(t)))+J,xm(t)=φ(t),t[d,0]. (2.3)

    Let (2.3) be the master system, then the slave system is

    S:{˙xs(t)=A(r(t))xs(t)+B(r(t))f(xs(t))+C(r(t))g(xs(td(t)))+J+u(t),xs(t)=ψ(t),t[d,0]. (2.4)

    Definition 1. [29] Master system (2.3) and slave system (2.4) are said to be asymptotic synchronization if

    limt||xm(t)xs(t)||=0

    for any initial conditions.

    In this paper, our objective is to design a suit controller u(t) such that the master system (2.3) and slave system (2.4) are in synchronization by using the sliding mode control method.

    To the end, writing the state error e(t)=xs(t)xm(t), and the corresponding state error system can be described by

    {˙e(t)=Aie(t)+BiF(e(t))+CiG(e(td(t)))ΔBif(xm(t))ΔCig(xm(td(t)))+u(t),e(t)=ψ(t)φ(t),t[d,0], (2.5)

    where F(e(t))=f(xs(t))f(xm(t)) and G(e(td(t)))=g(xs(td(t)))g(xm(td(t))).

    First, we take the sliding mode surface as

    σi(t)=e(t)+t0[(Ai+Ki)e(θ)BiF(e(θ))CiG(e(θd(θ)))]dθ. (2.6)

    Thus, the derivative of σi(t) is

    ˙σi(t)=˙e(t)+(Ai+Ki)e(t)BiF(e(t))CiG(e(td(t)))=Kie(t)ΔBif(xm(t))ΔCig(xm(td(t)))+u(t), (2.7)

    where Ki(1il) are some unknown matrices to be determined later. When the state trajectories reach the sliding mode surface, then ˙σi(t)=0 and σi(t)=0. So, we obtain the equivalent controller

    ueq(t)=Kie(t)+ΔBif(xm(t))+ΔCig(xm(td(t))). (2.8)

    Substituting (2.8) into (2.5), we have

    {˙e(t)=(Ai+Ki)e(t)+BiF(e(t))+CiG(e(td(t))),e(t)=ψ(t)φ(t),t[d,0]. (2.9)

    In this paper, the following assumptions for the neuron activation functions are needed.

    Assumption 1. [30] Assume that each component of the nonlinear function f() and g() are continuous and bounded and satisfy

    ηkfk(z1)fk(z2)z1z2η+k,
    θkgk(z1)gk(z2)z1z2θ+k,k=1,2,...,n,

    for any z1,z2R, where ηk>0,η+k>0,θk>0,θ+k>0 are some known positive constants.

    From Assumption 1, it is easy to get the following inequalities

    [eT(t)FT(e(t))][Θ1HΘ2HH][e(t)F(e(t))]0 (2.10)

    and

    [eT(t)GT(e(td(t)))][Θ3HΘ4HH][e(t)G(e(td(t)))]0, (2.11)

    where H=diag{h1,h2,...,hn} is a positive definite diagonal matrix,

    Θ1=diag{η1η+1,η2η+2,...,ηnη+n},Θ2=diag{η1+η+12,η2+η+22,...,ηn+η+n2},
    Θ3=diag{θ1θ+1,θ2θ+2,...,θnθ+n},Θ4=diag{θ1+θ+12,θ2+θ+22,...,θn+θ+n2}.

    Assumption 2.[31] Assume that each component of the nonlinear functions f() and g() are bounded, which means that there exists positive scalars Bf and Bg such that

    ||fk()||Bf,||gk()||Bg

    for k=1,2,...,n.

    Remark 1. In fact, the activation functions of neural networks are usually bounded. For example, the Logistic Sigmoid function h1(x)={1+eax}1 and the threshold value function

    h2(x)={1,x0,1,x<0,

    and so on.

    Throughout the paper, we need the following lemmas.

    Lemma 1.[32] For any positive definite symmetric matrix WRn×n and scalar τ>0, there is

    ttτxT(s)dsWttτx(s)dsτttτxT(s)Wx(s)ds.

    Lemma 2.[33] The linear matrix inequality

    [S11S12ST12S22]<0

    is equivalent to the following condition

    S22<0,S11S12S122ST12<0,

    where S11 and S22 are symmetric matrices.

    Now, we will analyze the synchronization condition and construct the sliding mode controller.

    Theorem 1. Under Assumption 1, if there exists positive definite symmetric matrices Pi,Wi,RiRn×n such that

    Φ1i=[Φ1i,1101dRiPiBiΘ2HPiCiΘ4HATiRiKTiRi(1μ)Wi00001dRi000H0BTiRiHCTiRi1dRi]<0 (3.1)

    for i=1,2,...,l, then the master system (2.1) and slave system (2.4) are in synchronization, where

    Φ1i,11=PiAiATiPiPiKiKTiPi+lj=1,jiπij(PjPi)+Wi1dRiΘ1HΘ3H.

    Proof. Constructing the following Lyapunov function

    V(t)=eT(t)Pie(t)+ttd(t)eT(θ)Wie(θ)dθ+0dtt+s˙eT(θ)Ri˙e(θ)dθds.

    The derivative of V(t) along with the trajectories of system (2.9) is

    ˙V(t)=2eT(t)Pi˙e(t)+eT(t)lj=1πijPje(t)+eT(t)Wie(t)(1˙d(t))eT(td(t))Wie(td(t))+d˙eT(t)Ri˙e(t)ttd˙eT(θ)Ri˙e(θ)dθ2eT(t)Pi[(Ai+Ki)e(t)+BiF(e(t))+CiG(e(td(t)))]+eT(t)lj=1,jiπij(PjPi)e(t)+eT(t)Wie(t)(1μ)eT(td(t))Wie(td(t))+d[(Ai+Ki)e(t)+BiF(e(t))+CiG(e(td(t)))]TRi[(Ai+Ki)e(t)+BiF(e(t))+CiG(e(td(t)))]1d[e(t)e(td)]TRi[e(t)e(td)].

    Letting ξ(t)=(eT(t),eT(td(t)),eT(td),FT(e(t))andGT(e(td(t)))T, one yields

    ˙V(t)ξT(t)[Φ2i+dΨTRiΨ]ξ(t), (3.2)

    where

    Φ2i=[Φ2i,1101dRiPiBiPiCi(1μ)Wi0001dRi00000],
    Ψ=[AiKi00BiCi],
    Φ2i,11=PiAiATiPiPiKiKTiPi+lj=1,jiπij(PjPi)+Wi1dRi.

    It is noted that

    ξT(t)[Φ2i+dΨTRiΨ]ξ(t)[eT(t)FT(e(t))][Θ1HΘ2HH][e(t)F(e(t))][eT(t)GT(e(td(t)))][Θ3HΘ4HH][e(t)G(e(td(t)))]=ξT(t)[Φ3i+dΨTRiΨ]ξ(t), (3.3)

    where

    Φ3i=[Φ1i,1101dRiPiBiΘ2HPiCiΘ4H(1μ)Wi0001dRi00H0H].

    By using Lemma 2, we know that Φ3i+dΨTRiΨ<0 is equivalent to Φ1i<0. Thus, one obtains ˙V(t)<0 and system (2.9) is asymptotically stable, which shows that the master system (2.1) and slave system (2.4) are in synchronization. The proof is completed.

    Theorem 2. Under Assumption 2 and the action of controller

    u(t)=Kie(t)ζ(t)sign(σi(t)), (3.4)

    master system (2.1) and slave system (2.4) are in synchronization, where

    ζ(t)=Bf||Mi||||N1i||+Bg||Mi||||N2i||+α,

    and α>0 is a positive scalar.

    Proof. Constructing the following Lyapunov function

    U(t)=12σTi(t)σi(t),

    then

    ˙U(t)=σTi(t)˙σi(t)=σTi(t)[Kie(t)ΔBif(xm(t))ΔCig(xm(td(t)))Kie(t)ζ(t)sign(σi(t))]||σi(t)||[||ΔBi||||f(xm(t))||+||ΔCi||||g(xm(td(t)))||]ζ(t)||σi(t)||. (3.5)

    Because of

    ||ΔBi||||f(xm(t))||+||ΔCi||||g(xm(td(t)))||||Mi||[Bf||N1i||+Bg||N2i||],

    then

    ˙U(t)α||σi(t)||=α2U(t).

    Thus, the state trajectories can attain the sliding mode surface in the finite time interval [0,T], where T2U(0)2α. The proof is completed.

    Remark 2. Compared with the uncertain system in [15], we especially considered the effect of time delay. From Theorems 1 and 2, we see that the designed sliding mode controller can realize the master-slave synchronization of the Markov jump neural networks.

    Remark 3. The control gain Ki can be obtained from matrix inequality (3.1). However, (3.1) is not a linear matrix inequality. In order to solve it, we can take Pi=Ri.

    From Theorem 1, we can obtain the following useful corollary as r(t)=1, which means that system (2.1) only has one mode.

    Corollary 1. Under Assumption 1, if there exist positive definite symmetric matrices P,W,RRn×n such that

    Φ1=[Φ1,1101dRPBΘ2HPCΘ4HATRKTR(1μ)W00001dR000H0BTRHCTR1dR]<0, (3.6)

    then the master system

    {˙xm(t)=Axm(t)+(B+ΔB)f(xm(t))+(C+ΔC)g(xm(td(t)))+J,xm(t)=φ(t),t[d,0] (3.7)

    and slave system

    {˙xs(t)=Axs(t)+Bf(xs(t))+Cg(xs(td(t)))+J+u(t),xs(t)=ψ(t),t[d,0] (3.8)

    are in synchronization, where

    Φ1,11=PAATPPKKTP+W1dRΘ1HΘ3H.

    Consider master system (2.1) and slave system (2.4) with the following parameters

    A1=[421052026],A2=[512070018],B1=[0.20.60.20.20.200.20.10.1],
    B2=[0.20.50.30.20.400.30.10.2],C1=[0.50.40.30.50.10.10.50.10.5],C2=[0.80.30.30.10.20.20.60.30.2],
    Π=[2211],M1=[0.20.50.7],M2=[0.30.80.4],
    N11=[0.10.20.1],N12=[0.050.10.03],
    N21=[0.40.10.3],N22=[0.10.020.07],μ=0.5,d=2,
    f(x(t))=[tanh(0.5x1(t)),tanh(0.4x2(t)),tanh(0.6x3(t))]T,
    g(x(t))=[tanh(0.3x1(t)),tanh(0.4x2(t)),tanh(0.2x3(t))]T,
    J1=0.5I3,J2=0.3I3,Θ1=Θ3=I3,Θ2=Θ4=I3,w(t)=sin(t).

    By using the linear matrix inequality (LMI) toolbox in the MATLAB, we obtain the following solutions of inequality (3.1):

    P1=[23.44057.37377.05257.373765.865820.94927.052520.949240.9154],P2=[21.21286.29397.45876.293961.528219.50707.458719.507038.4827],
    W1=[10.44523.21871.47953.218724.30057.12561.47957.125615.5240],W2=[5.79131.22002.45431.220016.47494.94762.45434.947610.7346],
    K1=[0.13660.00930.01880.00930.05500.02660.01880.02660.0902],K2=[0.13660.00930.01880.00930.05500.02660.01880.02660.0902],
    H=[76.637400059.824500045.2654].

    For the initial values xm(0)=(25,5,25)T and xs(0)=(30,5,25)T, Figure 1 is the Markovian jump process in different modes. Figure 2 is the state trajectories of error system (2.5), which shows that the error system is convergent and the master system synchronizes with the slave system. Figure 3 is the curve of the sliding mode surface, which shows that the state trajectories can arrive at the surface in a finite time interval.

    Figure 1.  Markovian jump process r(t).
    Figure 2.  The state trajectories of error system (2.5).
    Figure 3.  The curve of sliding mode surface (2.6).

    In this paper, the master-slave synchronization for uncertain neural networks with time delay by using the sliding mode control method has been studied. An integral sliding mode surface and sliding mode controller was designed. Moreover, the state trajectories of the neural networks can reach the sliding mode surface in finite time under the action of the controller. Sufficient conditions in terms of linear matrix inequalities were presented to guarantee the neural networks asymptotical stability. Finally, an example was provided to illustrate the validity of the proposed design method. In the future, we will consider how to solve some physical problems by applying the obtained theoretical results.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work is supported by the National Natural Science Foundation (NNSF) of China under Grant 11602134.

    The authors declare that they have no conflicts of interest.



    [1] H. Zhang, R. Ye, S. Liu, et al. LMI-based approach to stability analysis for fractional-order neural networks with discrete and distributed delays, Int. J. Syst. Sci., 49 (2018), 537-545. doi: 10.1080/00207721.2017.1412534
    [2] W. Zhang, J. Cao, R. Wu, et al. Lag projective synchronization of fractional-order delayed chaotic systems, J. Franklin Institute, 356 (2019), 1522-1534. doi: 10.1016/j.jfranklin.2018.10.024
    [3] W. Zhang, H. Zhang, J. Cao, et al. Synchronization in uncertain fractional-order memristive complex-valued neural networks with multiple time delays, Neural Networks, 110 (2019), 186-198. doi: 10.1016/j.neunet.2018.12.004
    [4] P. Dorato, Short time stability in linear time-varying systems, Proc. IRE Int. Convention Record, 4 (1961), 83-87.
    [5] D. F. Luo, Z. G. Luo, Uniqueness and novel finite-time stability of solutions for a class of nonlinear fractional delay difference systems, Discr. Dynam, Nature Soc., 2018 (2018), 1-7.
    [6] G. C. Wu, D. Baleanu, S. D. Zeng, Finite-time stability of discrete fractional delay systems: Gronwall inequality and stability criterion, Commun. Nonl. Sci. Numer. Simul., 57 (2018), 299-308. doi: 10.1016/j.cnsns.2017.09.001
    [7] V. N. Phat, N. T. Thanh, New criteria for finite-time stability of nonlinear fractional-order delay systems: A Gronwall inequality approach, Appl. Math. Lett., 83 (2018), 169-175. doi: 10.1016/j.aml.2018.03.023
    [8] D. F. Luo, Z. G. Luo, Existence and finite-time stability of solutions for a class of nonlinear fractional differential equations with time-varying delays and non-instantaneous impulses, Adv. Diff. Eq., 155, 2019.
    [9] D. Qian, C. Li, R. P. Agarwal, et al. Stability analysis of fractional differential system with Riemann-Liouville derivative, Math. Comput. Modell., 52 (2010), 862-874. doi: 10.1016/j.mcm.2010.05.016
    [10] M. Li, J. Wang, Finite time stability of fractional delay differential equations, Appl. Math. Lett., 64 (2017), 170-176. doi: 10.1016/j.aml.2016.09.004
    [11] M. Li, J. R. Wang, Exploring delayed Mittag-Leffler type matrix functions to study finite time stability of fractional delay differential equations, Appl. Math. Comput., 324 (2018), 254-265.
    [12] M. Li, J. R. Wang, Finite time stability and relative controllability of Riemann-Liouville fractional delay differential equations, Math. Meth. Appl. Sci., 2019 (2019), 1-17.
    [13] K. Diethelm, The Analysis of Fractional Differential Equations, Springer-Verlag, Berlin, Heidelberg, 2010.
    [14] I. Podlubny, Fractional Differential Equations, Academic Press: San Diego, 1999.
    [15] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, 2006.
    [16] H. Ye, J. Gao, Y. Ding, A generalized Gronwall inequality and its application to a fracnal differential equation, J. Math. Anal. Appl., 328 (2007), 1075-1081. doi: 10.1016/j.jmaa.2006.05.061
    [17] R. Agarwal, S. Hristova, D. O'Regan, Explicit solutions of initial value problems for linear scalar Riemann-Liouville fractional differential equations with a constant delay, Mathematics, 8 (2020), 1-14.
  • Reader Comments
  • © 2020 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)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4219) PDF downloads(402) Cited by(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog