Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Dissipativity and contractivity of the second-order averaged L1 method for fractional Volterra functional differential equations

  • Received: 16 December 2022 Revised: 10 February 2023 Accepted: 16 February 2023 Published: 03 March 2023
  • This paper focuses on the dissipativity and contractivity of a second-order numerical method for fractional Volterra functional differential equations (F-VFDEs). Firstly, an averaged L1 method for the initial value problem of F-VFDEs is presented based on the averaged L1 approximation for Caputo fractional derivative together with an appropriate piecewise interpolation operator for the functional term. Then the averaged L1 method is proved to be dissipative with an absorbing set and contractive with an algebraic decay rate. Finally, the numerical experiments further confirm the theoretical results.

    Citation: Yin Yang, Aiguo Xiao. Dissipativity and contractivity of the second-order averaged L1 method for fractional Volterra functional differential equations[J]. Networks and Heterogeneous Media, 2023, 18(2): 753-774. doi: 10.3934/nhm.2023032

    Related Papers:

    [1] Wenrui Guan, Xun Wang . A Dual-channel Progressive Graph Convolutional Network via subgraph sampling. Electronic Research Archive, 2024, 32(7): 4398-4415. doi: 10.3934/era.2024198
    [2] Yongmei Zhang, Zhirong Du, Lei Hu . A construction method of urban road risky vehicles based on dynamic knowledge graph. Electronic Research Archive, 2023, 31(7): 3776-3790. doi: 10.3934/era.2023192
    [3] Chuang Ma, Helong Xia . A one-step graph clustering method on heterogeneous graphs via variational graph embedding. Electronic Research Archive, 2024, 32(4): 2772-2788. doi: 10.3934/era.2024125
    [4] Xinzheng Xu, Xiaoyang Zhao, Meng Wei, Zhongnian Li . A comprehensive review of graph convolutional networks: approaches and applications. Electronic Research Archive, 2023, 31(7): 4185-4215. doi: 10.3934/era.2023213
    [5] Zhenyu Yang, Lei Wu, Peian Wen, Peng Chen . Visual Question Answering reasoning with external knowledge based on bimodal graph neural network. Electronic Research Archive, 2023, 31(4): 1948-1965. doi: 10.3934/era.2023100
    [6] Suhua Wang, Zhen Huang, Hongjie Ji, Huinan Zhao, Guoyan Zhou, Xiaoxin Sun . PM2.5 hourly concentration prediction based on graph capsule networks. Electronic Research Archive, 2023, 31(1): 509-529. doi: 10.3934/era.2023025
    [7] Bingjie Zhang, Junchao Yu, Zhe Kang, Tianyu Wei, Xiaoyu Liu, Suhua Wang . An adaptive preference retention collaborative filtering algorithm based on graph convolutional method. Electronic Research Archive, 2023, 31(2): 793-811. doi: 10.3934/era.2023040
    [8] Xia Liu, Liwan Wu . FAGRec: Alleviating data sparsity in POI recommendations via the feature-aware graph learning. Electronic Research Archive, 2024, 32(4): 2728-2744. doi: 10.3934/era.2024123
    [9] Hui Xu, Jun Kong, Mengyao Liang, Hui Sun, Miao Qi . Video behavior recognition based on actional-structural graph convolution and temporal extension module. Electronic Research Archive, 2022, 30(11): 4157-4177. doi: 10.3934/era.2022210
    [10] Ruxin Xue, Jinggui Huang, Zaitang Huang, Bingyan Li . Reconstructed graph spatio-temporal stochastic controlled differential equation for traffic flow forecasting. Electronic Research Archive, 2025, 33(4): 2543-2566. doi: 10.3934/era.2025113
  • This paper focuses on the dissipativity and contractivity of a second-order numerical method for fractional Volterra functional differential equations (F-VFDEs). Firstly, an averaged L1 method for the initial value problem of F-VFDEs is presented based on the averaged L1 approximation for Caputo fractional derivative together with an appropriate piecewise interpolation operator for the functional term. Then the averaged L1 method is proved to be dissipative with an absorbing set and contractive with an algebraic decay rate. Finally, the numerical experiments further confirm the theoretical results.



    1. Introduction

    Time delay systems represent one of the most popular class of systems. Time delay, whether occurs in the system state, the control input, or the measurement, is often inevitable in practical systems and can be a source of instability and poor performance [21,19,1,61,41,20,94,95,79,72]. The future evolution of the system state of a time delay system depends not only on its current value, but also on its past values [124,15,88,87,55,121,111,81,34,14,117,46]. Many processes have time delay characteristics in their dynamics. Since time delays often appear in engineering, biological and economical systems, and sometimes they may poorly affect the performance of a system. The problem of stability of IDSs and impulsive stabilization of delay systems have been extensively investigated [119,108,104,96,36,98,18,24]. For example, [104] studied the stability of a class of nonlinear impulsive switching systems with time-varying delays. Based on the common Lyapunov function method and Razumikhin technique, several stability criteria are established for nonlinear impulsive switching systems with time-varying delays. In [98], by structuring hybrid impulsive and feedback controllers, synchronization problem of the memristive delayed neural networks is proposed. Then, based on differential inclusions, several synchronization criteria for the memristive delayed neural networks are obtained by impulsive control theories, special inequalities and the Lyapunov-type functional. In literatures, the research results concerning time delay systems can be classified into two types. One is delay-independent conditions, the other is delay-dependent conditions. Delay-dependent conditions are less conservative compared with the delay-independent conditions because they incorporate the information of time delays. Various techniques have been developed in literatures to derive delay-dependent conditions, such as the Lyapunov-Razumikhin method and the Lyapunov-Krasovskii functional method [18]. In [24], the authors proposed the equivalence between stability conditions for switched systems and the Lyapunov-Krasovskii functional stability conditions for discrete-time delay systems. This provides us another method to investigate time delay systems.

    Impulsive effects exist widely in the world. As we know, the state of systems are often subject to instantaneous disturbances and experience abrupt changes at certain instants, which might be caused by frequent changes or other suddenly noises. These systems are called impulsive systems, which are governed by impulsive differential equations or impulsive difference equations [42,118,59,57,2]. In the past decades, there has been a growing interest in the theory of impulsive dynamical systems because of their applications to various problems arising in communications, control technology, impact mechanics, electrical engineering, medicine, biology, and so on [54,28,22,17,33,32,106,3,112,74,100,75,83,31]. For example, [54] investigated the pth moment exponential stability of impulsive stochastic functional differential equations. Some sufficient conditions are obtained to ensure the pth moment exponential stability of the equilibrium solution by the Razumikhin method and Lyapunov functions. In [100], the authors proposed a IDS model for insulin therapy for both type 1 and type 2 diabetes mellitus with time delay in insulin production. It is shown that impulsive exogenous insulin infusions can mimic natural pancreatic insulin production. From the viewpoint of impulsive effects, the stability analysis for impulsive time-delay systems can be generally classified into two groups: impulsive stabilization and impulsive perturbation. The case where a given system without impulses is unstable or stable and can be turned into uniformly stable, uniformly asymptotically stable, and even exponentially stable under proper impulsive control, it is regarded as impulsive stabilization problem. Now it has been shown that impulsive stabilization problems can be widely applied to many fields such as orbital transfer of satellite, dosage supply in pharmacokinetics, ecosystems management, and synchronization in chaotic secure communication systems [35,116,9]. Alternatively, the case where a given system without impulses is stable and can remain stable under certain impulsive interference, it is regarded as impulsive perturbation problem. Up to now, many interesting results dealing with impulsive perturbations of time-delay systems have been reported [89,91,13].

    By analyzing related literatures, this paper provides a comprehensive and intuitive overview for IDSs, which include the basic theory of IDSs, stability analysis with impulsive control and impulsive perturbation, and delayed impulses. Essentially, it provides an overview on the progress of stability and stabilization problem of IDSs. The rest of this paper is organized as follows. In Section 2, some notations and definitions of stability are presented. Section 3 covers the effects of impulses for IDSs. Section 4 considers the impulsive control problem. Section 5 considers the impulsive perturbation problem. Section 6 covers the delayed impulses. Section 7 concludes the paper and discuss the future research direction on this topic.


    2. Impulsive delayed systems

    Notations. Let R denote the set of real numbers, R+ the set of positive numbers, Z+ the set of positive integer, N the set of nonnegative integer. Rn is the n-dimensional real spaces equipped with the Euclidean norm ||||. Rn×m denotes the n×m-dimensional real spaces. I denotes the identity matrix with appropriate dimensions. For any interval JR, set SRk(1kn), C(J,S)={φ:JS is continuous} and C1(J,S)={φ:JS is continuously differentiable}. PC(J,S)={φ:JS is continuous everywhere except at finite number of points t, at which φ(t+), φ(t) exist and φ(t+)=φ(t)}. In particular, for given τ>0, set Cτ=C([τ,0],Rn) and PCτ=PC([τ,0], Rn) with the norm ||||τ defined by ||φ||τ=sup{||φ(s)||: s[τ,0]}. Λ={1,2,,n}. S(ρ)={xRn: |x|<ρ}. A function α:[0,)[0,) is of class K if α is continuous, strictly increasing, and α(0)=0. In addition, if α is unbounded, then it is of class K.

    In recent years, there are many results on impulsive delayed systems (IDSs). Roughly speaking, an impulsive dynamical system consists of three elements: a continuous-time dynamical equation, which governs the evolution of the system between reset (impulsive) events; a difference equation, which describes the way the system states are instantaneously changed; and finally a criterion for determining when the states of the system are to be reset. In addition, it is well known that time-delays phenomena frequently appear in many practical problems, such as biological systems, mechanical, transmissions, fluid transmissions, networked control systems [122,47,123,86,80,25,45]. Therefore, it is not surprising that IDSs with time delays have become an attractive research field. In the following, consider the impulsive functional differential equation

    {˙x(t)=f(t, xt), t>t0, ttk,Δx=Ik(t, x(t)), t=tk, kN, (1)

    where f,Ik:R+×PCτRn, Δx=x(tk)x(tk). For each tt0, xtPCτ is defined by xt(s)=x(t+s),s[τ,0]. The impulse times {tk} satisfy 0t0<<tk<, limktk=. With Eq. (1), one associates an initial condition of the form

    xt0=ϕ, (2)

    where t0R+ and ϕPCτ.

    Definition 2.1. [85]. A function x(t) is called a solution of the initial value problem (1) and (2) if x:[t0τ, β)Rn, for some β (t0β<), is continuous for t[t0τ,β){tk, kΛ}, x(tk) and x(t+k) exist, and x(t+k)=x(tk) for tk[t0τ, β), and satisfies (1) and (2).

    Under the following hypotheses (H1)(H4), the initial value problem (1) and (2) exists with a unique solution which will be written in the form x(t,t0,ϕ), see [85] for detailed information.

    (H1) f is continuous on [tk1, tk)×PCτ for each kZ+, and for all φPCτ and kZ+, the limit lim(t,ϕ)(tk,φ)f(t,ϕ)=f(tk,φ) exists.

    (H2) f is locally Lipschitzian in ϕ in each compact set in PCτ.

    (H3) For each kZ+, Ik(t,x)C(R+×S(ρ),Rn).

    (H4) There exists a ρ1>0 (ρ1<ρ) such that xS(ρ1) implies that x+Ik(tk,x)S(ρ) for all kZ+.

    Assume that conditions (H1)(H4) hold and moreover, f(t,0)0 and Ik(tk,0)0, then x(t)0 is the solution of (1), which is called the trivial solution.

    Definition 2.2. [56]. The trivial solution of system (1) is said to be

    (S1) stable if for any ε>0 and t0R+, there exists some δ=δ(ε, t0)>0, such that ||ϕ||τ<δ implies ||x(t, t0, ϕ)||<ε, for all t>t0;

    (S2) uniformly stable if the δ in (S1) is independent of t0;

    (S3) asymptotically stable if (S1) hods and for any t0R+, there exists some η=η(t0)>0 such that if ϕPCτ with ||ϕ||τ<η, then limtx(t,t0,ϕ)=0;

    (S4) uniformly asymptotically stable if (S2) holds and there exists some η>0 such that for every γ>0, there exists some T=T(η,γ)>0 such that if ϕPCτ with ||ϕ||τ<η, then ||x(t, t0, ϕ)||<γ for tt0+T;

    (S5) globally stable (globally exponentially stable) if (S1) holds and for any δ, there exist λ>0 and k(δ)>0, when ||ϕ||τ<δ we have

    ||x(t, t0, ϕ)||k(δ)eλ(tt0), tt0.

    Definition 2.3. [38]. A map x: R+Rn is said to be an ωperiodic solution of the system (1) provided:

    (ⅰ) x satisfies (1) and is a piecewise continuous map with first-class discontinuity points;

    (ⅱ) x satisfies x(t+ω)=x(t) for tτk and x(τk+ω+)=x(τ+k) for kN.

    Definition 2.4. [38]. Let x=x(t, t0, ϕ) be an ω-periodic solution of the system (1) with initial condition (t0, ϕ). Then x is said to be globally attractive if for any solution x=x(, t0, ϕ) of the system (1) through (t0, ϕ), |xx|0 as t.

    Definition 2.5. [109]. The function V: [T0τ,)×RnR+ belongs to class v0 if

    (1) V is continuous on each of the sets [tk1,tk)×R+ and lim(t,u)(tk,v)V(t,u)=V(tk,v);

    (2) V(t,x) is locally Lipschitzian in x and V(t,0)0, tt0.

    Definition 2.6. [109]. Let Vv0, the upper right-hand derivative of V with respect to system (1) is defined by

    D+V(t,ψ(0))=limh0+sup1h[V(t+h,ψ(0)+hf(t,ψ))V(t,ψ(0))],

    for (t,ψ)[tk1,tk)×PCτ.


    3. Effects of impulses

    Generally speaking, existing results on stability for IDSs can be classified into two groups: impulsive stabilization and impulsive perturbation. In the case where a given equation without impulses is unstable or stable, it can be tended to uniformly stable, uniformly asymptotically stable even exponentially stable under proper impulsive control. Such case is regarded as impulsive stabilization. In the case where a given equation without impulses is stable, and it can remain the stability behavior under certain impulsive interference, it is regarded as impulsive perturbation. At each discontinuous point tk, suppose that V(tk)μV(tk), where the constant μ represents the impulsive strength. There are three kinds of impulsive strength μ for IDSs:

    (1) When |μ|<1, the impulses are beneficial for the stability of the system. Hence we call them beneficial impulses when |μ|<1, which can be categorized as impulsive controllers that would enhance the stability of the system.

    (2) When |μ|=1, the impulses are neither harmful nor beneficial for the stability of IDSs since the absolute value of the differences between the nodes' states remains unchanged. Hence, the impulses with |μ|=1 are named as inactive impulses.

    (3) When |μ|>1, the impulses might destroy the stability of the system. Thus we call them harmful impulses if |μ|>1, which can be regarded as the perturbations that would suppress the stability of the system.

    In the following, two examples are given to illustrate the effects of impulses for the IDSs.

    Example 3.1. Consider a simple IDS:

    {˙x(t)=0.5x(t)+0.6x(t1),t0,x(k)=μx(k),kN, (3)

    where μ=0.4. Figs. 3.1 (a) and (b) show the time response of states x(t) of (3) with or without impulsive control, respectively. Obviously, it can be seen from the figures that the system (3) is unstable when there is no impulsive control and becomes stable under proper impulsive control, which can be obtained by the results in [48]. It implies that impulses contribute to system dynamics.

    Figure 3.1. (a) State trajectory of system (3) without impulse control. (b) State trajectory of system (3) with impulsive stabilization.

    Example 3.2. Consider another IDS:

    {˙x(t)=0.7x(t)+0.35x(t1),t0,x(k)=μx(k),kN, (4)

    where μ=1.3. Figs. 3.2 (a) and (b) show the time response of states x(t) of (4) with or without impulsive perturbation, respectively. It can be seen from the figures that the system (4) is stable when there is no impulsive perturbation and remains stable under proper impulsive perturbation.

    Figure 3.2. (a) State trajectory of system (4) without impulsive perturbation. (b) State trajectory of system (4) with impulsive perturbation.

    The above two examples fully illustrate the different effects of impulse for stability on the IDSs. In resent years, there are many researches on stability analysis for IDSs. For example, in [63], criteria on uniform asymptotic stability were established for impulsive delay differential equations by using Lyapunov functions and Raxumikhin techniques. [44] presented some sufficient conditions for global exponential stability for a class of delay difference equations with impulses by means of constructing an extended impulsive delay difference inequality. In [39], authors addressed the impulsive systems with unbounded time-varying delay and introduced a new impulsive delay inequality that involves unbounded and non-differentiable time-varying delay. Some sufficient conditions ensuring stability and stabilization of impulsive time-invarying and time-varying systems are derived, respectively. [115] investigated the synchronization problem of coupled switched neural networks (SNNs) with mode-dependent impulsive effects and time delays. The impulses considered here include those that suppress synchronization or enhance synchronization. Based on switching analysis techniques and the comparison principle, the exponential synchronization criteria are derived for coupled delayed SNNs with mode-dependent impulsive effects. In addition, the concept of ''average impulsive interval" was introduced in [68] by referring to the concept of average dwell time [23] to characterize how often or how seldom impulses occur.

    Definition 3.3. [68]. The average impulsive interval (AII) of the impulsive sequence ξ={t1,t2,} is equal to τα, if there exist positive integer N0 and positive number τα such that

    TtταN0Nξ(T,t)Ttτα+N0,  Tt0, (5)

    where Nξ(T,t) denotes the number of impulsive times of the impulsive sequence ξ on the interval [t,T). The idea behind it is that there may exist some consecutive impulse signals separated by less than or greater than τα, but the average interval between consecutive impulse signals is τα.

    For most impulsive signals, the occurrence of impulses is not uniformly distributed. Fig.3.3 presents a specific form of a non-uniformly distributed impulsive sequence. One may observe from Fig. 3.3 that the impulses seldom occur in some time intervals, but frequently occur in some other intervals. For such impulsive signals, it is possible that the lower bound of the impulsive intervals is small or the upper bound is quite large. Hence, many previous results cannot be effectively applied to dynamical systems with the impulsive signal shown in Fig. 3.3.

    Figure 3.3. Impulsive sequence.

    In addition, note that inequality (5) can be rewritten as follows:

    Nξ(T,t)TtταN0, Tt0 (6)
    Nξ(T,t)Ttτα+N0, Tt0. (7)

    When the original system without impulsive perturbation is stable, and the impulsive effects are harmful, in order to guarantee the stability, the impulses should not occur frequently. That is, there should always exist a requirement that tktk1>h1, where the quantity h1 can be regarded as a measure to ensure that the harmful impulses do not occur too frequently. Thus, condition (7) enforces an upper bound on the number of impulses. Conversely, when the original system without impulsive effects is unstable and the impulsive effects are beneficial, in order to ensure the stability of the system, it is usually assumed that the frequency of impulses should not be too low. Therefore, there should always exist a requirement that tktk1<h2, where the quantity h2 is chosen to guarantee that there will be no overly long impulsive intervals. Thus condition (6) enforces a lower bound on the number of impulses.

    Hence, the concept of AII is suitable for characterizing a wide range of impulsive signals. In recent years, there are many results using the concept of AII [110,7,71,68,32]. For instance, the global exponential synchronization of delayed complex dynamical networks with nonidentical nodes and stochastic perturbations was studied in [110]. By combining adaptive control and impulsive control schemes, the considered network can be synchronized onto any given goal dynamics. With respect to impulsive control, the concept named AII with ''elasticity number'' of impulsive sequence is utilized to get less conservative synchronization criterion. [7] investigated the problems of impulsive stabilization and impulsive synchronization of discrete-time delayed neural networks (DDNNs), where two types of DDNNs with stabilizing impulses were studied. [71] established finite time stability (FTS) criteria for the nonlinear impulsive systems, where by using AII method, less conservative conditions were obtained for the FTS problem on the impulsive systems.


    4. Impulsive stabilization on IDSs

    Impulsive control is to change the state of a system by discontinuous control input at certain time instances. From the control point of view, impulsive control is of distinctive advantage, since control gains are only needed at discrete instances. Thus, there are many interesting results on impulsive stabilization of IDSs. Considering the IDSs (1), authors in [63] established some criteria on uniform asymptotic stability by using Lyapunov functions and Razumikhin techniques, which is given as follows:

    Theorem 4.1. [63]. Assume that there exist functions a, b, cK, pPC(R+,R+), g,ˆgK, and Vv0, such that the following conditions hold:

    (ⅰ) b(||x||)V(t,x)a(||x||), for all (t,x)[t0τ,+)×Rn;

    (ⅱ) D+V(t,ψ(0))p(t)c(V(t,ψ(0))), for all ttk in R+ and ψPCτ whenever V(t, ψ(0))g(V(t+s,ψ(s))) for s[τ,0];

    (ⅲ) V(tk,ψ(0)+Ik(tk,ψ))g(V(tk,ψ(0))), for all (tk,ψ)R+×PCτ for which ψ(0)=ψ(0);

    (ⅳ) τ0=supkZ{tktk1}<, M1=supt0t+ιtp(s)ds<, M2=infq>0qg(q)dsc(s)>M1.

    Then the trivial solution of (1) is uniformly asymptotically stable.

    It follows from the definitions M1 and M2 in (iv) that the value of M2 is greater than M1, which indicates that there may be some increasing on V between impulses, but the decreasing on V at impulses are needed to guarantee the stability. Notice from (iv) that there is a restriction on upper bound of impulse intervals due to the effect of stabilizing impulses. In other words, the time interval without impulsive control cannot be too long. Similar results on stability without time delay can be found in [62].

    Based on the ideas given in [63], authors in [105] and [103] further investigated exponential stability and global exponential stability of solutions for IDSs (1), respectively, which play important effects on exponential stability analysis of impulsive time-delay systems.

    Theorem 4.2. [103]. Assume that there exist a function Vv0 and constants p,c,c1,c2>0 and α>τ, λ>c such that

    (ⅰ) c1||x||pV(t,x)c2||x||p, for any tR+ xRn;

    (ⅱ) D+V(t,ψ(0))cV(t,ψ(0)), for all t[tk1,tk),kN, whenever qV(t, ψ(0))V(t+s,ψ(s)) for s[τ,0], where q>e2λα is a constant;

    (ⅲ) V(tk,ψ(0)+Ik(tk,ψ))dkV(tk,ψ(0)), where dk>0, kN, are constants;

    (ⅳ) τtktk1α and ln(dk)+λα<λ(tk+1tk).

    Then the trivial solution of (1) is globally exponentially stable and the convergence rate is λp.

    For all the above studies [63,105,103], authors have investigated for the uniform asymptotical stability and global exponential stability of IDSs under the assumption that τ<tktk1, where τ is the finite delay or some positive constant. From this point of view, author in [36] studied the globally exponential stabilization of impulsive functional differential equations with infinite delays or finite delays by using Lyapunov functions and improved Razumikhin technique, where there is no any restriction on the lower bound of the impulse interval. Hence, the obtained results in [36] have wider applications.

    Consider the following impulsive functional differential equations:

    {˙x(t)=F(t, x()), t>t0, ttk,Δx=Ik(t, x(t)), t=tk, kN,xt0=ϕ(s), αs0. (8)

    Some detailed information can be found in [36]; here we omit it.

    Theorem 4.3. [36]. Assume that there exist a function Vv0 and constants p>0, q>1, c1>0, c2>0, m>0 and γ>0 such that

    (ⅰ) c1||x||mV(t,x)c2||x||m, for any (t,x)[0,+)×Rn;

    (ⅱ) For any t0R+ and ψPCτ, if eγθV(t+θ,ψ(θ))qV(t,ψ(0)), θ[α,0], ttk, then D+V(t,ψ(0))pV(t,ψ(0));

    (ⅲ) For all (tk,ψ)R+×PCτ, V(tk,ψ(0)+Ik(tk,ψ))1/qV(tk,ψ(0));

    (ⅳ) supkZ+tktk1<lnq/p.

    Then the trivial solution of (8) is globally exponentially stable.

    It should be noted that in [36], for tt00>α, f(t,xt) in IDSs (1) was replaced by f(t,x(s)), where s[t+α,t] or f(t,x()) be a Volterra type functional. In [36], let α=. Moreover, it follows from [36] that Theorem 4.3 can be applied to systems with finite or/and infinite delays since α[,0]. Based on this point, [49] considered the existence, uniqueness, and global stability of periodic solutions for a class of recurrent neural networks with discrete and continuously distributed delays. By using contraction mapping theorem, some new sufficient conditions ensuring the existence, uniqueness, and global stability of periodic solutions were obtained. As we know, since the periodic oscillations can be presented by systems model with periodic coefficients, the properties of periodic solution are very important to study the dynamical behaviors of systems. With the development of impulsive control theory, it is known that an equilibrium point can be viewed as a special periodic solution of impulsive control systems with arbitrary period. Hence, the analysis of periodic solutions of impulsive control systems is more general than that of an equilibrium point.

    Hence, investigation of periodic solution for system is indispensable for practical design and engineering applications of models. The periodic solution problem of IDSs has found many applications such as associative memories, pattern recognition, machine learning, robot motion control, and so on [76,97,82,27]. For example, in [49], a class of recurrent neural networks with discrete and continuously distributed delays was considered. Sufficient conditions for the existence, uniqueness, and global exponential stability of a periodic solution were obtained by using contraction mapping theorem and stability theory on impulsive functional differential equations. [38] dealt with the periodic solutions problem for impulsive differential equations. By using Lyapunov's second method and the contraction mapping principle, some conditions ensuring the existence and global attractiveness of unique periodic solutions were derived, and so on [37,84,101,26].

    In addition, the stability analysis is much more complicated because of the existence of impulsive effects and stochastic effects at the same time. In [12], based on the Razumikhin techniques and Lyapunov functions, several criteria on the global exponential stability and instability of impulsive stochastic functional differential systems were obtained. The results show that impulses make contribution to the exponential stability of stochastic differential systems with any time delay even they are originally unstable. In [77], authors investigated the pth moment and almost sure exponential stability of impulsive stochastic functional differential equations with finite delay by using Lyapunov method. The obtained results do not need the strong condition of impulsive gain |dk|<1, which is more general than those given in [12]. By using Lyapunov functions method [77] and stochastic analysis approaches, stability theorems were derived in terms of linear matrix inequality (LMI), which can overcome the effect of time-delay and impulses. Those results can guarantee the neural networks to be robustly exponentially stable in the mean square. [78] studied robust stability for a class of uncertain stochastic neural networks, where the 'uncertainties' mean the uncertain parameters, which take values in some intervals. [6] was devoted to prove some sufficient conditions ensuring practical pth moment exponential and almost sure exponential stability of solutions to impulsive stochastic functional differential equations.

    In what follows, an example is given to illustrate the existence and attractiveness of the periodic solution for impulsive control system.

    Example 4.4. [49]. Consider the following neural networks with discrete delays:

    {˙x1(t)=[0.15+0.05sin2πωt]x1(t)+2j=1[0.30.01cos2πω(t+j)])×fj(xj(t)+2j=1[0.1+0.04sin2πω(t+j)]fj(xj(t1))+cos2πωt, ttk,˙x2(t)=[0.2+0.1sin2πωt]x2(t)+2j=1[0.17+0.02sin2πω(t+j)]×fj(xj(t))+2j=1[0.1+0.01sin2πω(t+j)]fj(xj(t1))+sin2πωt, ttk, (9)

    subject to impulses:

    xi(tk)=1ρxi(tk),kZ+, i=1,2,

    where f1=f2=tanh(x),ω>0 and ρ>1 are some real constants.

    In the simulations, one may observe that system (9) with ω=2 or 4 has no periodic solution which is globally exponentially stable when there is no impulsive effect, which can be seen from Fig.4. 2.(a-d). However, via the impulsive control strategies that we established, system (9) may admit a unique periodic solution which is globally exponentially stable. For instance, when ω=2 or 4, one may choose ρ=1.5 and μ=0.5 such that the conditions in [49] hold and so the periodic solution can be guaranteed. The corresponding numerical simulations are shown in Fig.4. 3.(a-d).

    Figure 4.2. (a) State trajectories of system (9) with ω = 2 without impulsive control. (b) Phase portraits of system (9) with ω = 2 without impulsive control. (c) State trajectories of system (9) with ω = 8 without impulsive control. (d) Phase portraits of system (9) with ω = 8 without impulsive control.
    Figure 4.3. (a) State trajectory of system (9) with ω=2,ρ=1.5,μ=0.5. (b) Phase portraits of system (9) with ω=2,ρ=1.5,μ=0.5. (c) State trajectory of system (9) with ω=8,ρ=1.5,μ=0.5. (d) Phase portraits of system (9) with ω=8,ρ=1.5,μ=0.5.

    5. Impulsive perturbation on IDSs

    In this case where a given equation without impulses is stable, and it can remain the stability behavior under certain impulsive interference, it is regarded as impulsive perturbation problem. Considering the IDSs (1), authors in [63] established a criteria on uniform asymptotic stability, where there exists impulsive perturbation.

    Theorem 5.1. [63]. Assume that there exist functions a, b, cK, pPC(R+,R+), g,ˆgK, and Vv0, such that the following conditions hold:

    (ⅰ) b(||x||)V(t,x)a(||x||), for all (t,x)[t0τ,+)×Rn;

    (ⅱ) D+V(t,ψ(0))p(t)c(V(t,ψ(0))), for all ttk in R+ and ψPCτ whenever g(V(t, ψ(0)))V(t+s,ψ(s)) for s[τ,0];

    (ⅲ) V(tk,ψ(0)+Ik(tk,ψ))ˆg(V(tk,ψ(0))), for all (tk,ψ)R+×PCτ for which ψ(0)=ψ(0);

    (ⅳ) μ=infkZ{tktk1}>0, M2=supq>0g(q)q(ds)/(c(s)), M1=inft0t+μtp(s)ds>M2,

    where s<ˆg(s)<g(s). Then the trivial solution of (1) is uniformly asymptotically stable.

    Theorem 5.1 is in some ways the opposite of Theorem 4.1. Here the derivative of V is always non-positive, which implies that it is nonincreasing along solutions between impulses. In the absence of impulses, the trivial solution of system (1) is uniformly asymptotically stable. Theorem 5.1 allows for significant increases in V, and correspondingly the solutions themselves, at impulse times but only as long as these are balanced sufficiently by the decrease of V between impulses. These techniques are based in part on earlier work in the study of boundedness properties of solutions of impulsive differential equations without delay (see [4,5]). Years later, Liu et al. employed the method of Lyapunov functionals for the study of exponential stability of impulsive systems with time delay in [66]. Several exponential stability criteria were established as follows:

    Theorem 5.2. [66]. Assume that hypotheses (H1)(H3) in [66] are satisfied and there exist V1v0, V2v0(), p1, p2>0 with p1p2 and constants α, l, c, c1, c2, c3>0, dk0, kN, such that

    (ⅰ) c1||x||p1V1(t,x)c2||x||p1, 0V2(t,ψ)c3||ψ||p2τ, for any tR+ xRn, ψPCτ;

    (ⅱ) V1(tk,x+Ik(tk,x))dkV(tk,x), where xRn, kN;

    (ⅲ) for V(t,ψ)=V1(t,ψ(0))+V2(t,ψ), D+V(t,ψ)cV(t,ψ), for t[tk1,tk),ψPCτ,kN,

    (ⅳ) For any kN, τtktk1l, and ln(dk+c3c1e(p2p11)ckl)(α+c)l.

    Then the trivial solution of (1) is exponentially stable.

    Based on the idea in [66], ref. [60] investigated input-to-state stability (ISS) and integral input-to-state stability (iISS) of impulsive and switching hybrid systems with time-delay, using the method of multiple Lyapunov-Krasovskii functionals. It is also shown that the results in the present paper can be applied to systems with arbitrarily large delays and, therefore, improve the results in [66]. In addition, [68] presented a new concept AII, which can be used to describe impulsive signals with a wider range of impulsive interval. Based on this concept, a unified synchronization criterion was obtained for impulsive directed dynamical networks with desynchronizing impulses or synchronizing impulses. However, there is no time delay in it. [114] addressed the stability problem of a class of delayed neural networks with time-varying impulses. Different from the results in Lu et al. [68], the impulses considered here are time-varying, which can cover the results in Lu et al. [68]. From the obtained results in [114], we can see that even if destabilizing impulsive effects occur frequently, the delayed neural networks can also be stable if stabilizing impulses can prevail over the influence of destabilizing impulsive effects.

    To data, impulsive perturbation problem has been widely investigated [89,91,13,113,92,40,65,58,38,43,73,67,102,90,69,70,120]. For example, in [113], the global exponential stability of complex-valued impulsive systems was addressed. Some new sufficient conditions were obtained to guarantee the global exponential stability by the Lyapunov-Razumikhin theory. Authors in [92] studied the problem of impulsive effects on global exponential stability for a class of impulsive n-dimensional neural networks with unbounded delays and supremums. The robust exponential stability of nonlinear impulsive switched delayed systems was investigated in [67]. In [48], the stability problem of impulsive functional differential equations (IFDEs) was considered. Several criteria ensuring the uniform stability of IFDEs with finite or infinite delay were derived by establishing some new Razumikhin conditions. authors investigated the stability problem of time-delay systems with persistent impulses and focused on the discussion of systems with unbounded time-varying delays in [40]. The above results have been discussed the stability for impulsive perturbation systems, where the impulsive perturbation exists in system but not destroy the stability property.

    Sometimes, the impulsive perturbation can change the dynamics of a time-delay system. For example, a stable equilibrium may becomes a periodic attractor under impulsive perturbation. To show this observation, we illustrate the following example:

    Example 5.3. Consider the 2D impulsive delayed neural networks

    {˙x(t)=Ax(t)+Bf(x(tτ)),tt0,x(k)=Kx(tk),kN, (10)

    where

    A=(0.9000.9),  B=(0.5110.4),  f(x(tτ))=tanh(x(tτ)), τ=1.

    One may observe that when there is no impulsive perturbation on the system (10), i.e., K=I, the system (10) admits a stable equilibrium, see 5.3 (a). When there is impulsive perturbation such as K=1.2I and tk=0.8, kN, then the system (10) admits a periodic attractor, see Fig. 5.3 (b). Hence, such example shows that impulse can change the dynamics of time-delay systems and leads to complex dynamics.

    Figure 5.3. (a) State trajectories of (10) without impulse. (b) State trajectories of system (10) with impulsive perturbation.

    6. Delayed impulse

    Over the past few decades, many stability criteria for IDSs have been proposed. However, most existing results on IDSs do not consider the effect of delayed impulses. Of current interest is the delayed impulses of dynamical systems arising in such applications as automatic control, secure communication and population dynamics, [64,29,10,8,30]. Delayed impulse describes a phenomenon where impulsive transients depend on not only their current but also historical states of the system, see [16,107,11,52,51]. For instance, in communication security systems based on impulsive synchronization [8,30], there exist transmission and sampling delays during the information transmission process, where the sampling delay created from sampling the impulses at some discrete instances causes the impulsive transients depend on their historical states. Another example, in population dynamics such as fishing industry [52,51,99,50], effective impulsive control such as harvesting and re- leasing can keep the balance of fishing, and the quantities of every impulsive harvesting or releasing are not only measured by the current numbers of fish but also depend on the numbers in recent history due to the fact that the immature fish need some time to grow.

    In the previous literature concerning impulsive systems, the impulses are usually assumed to take the form:

    Δx(tk)=x(tk)x(tk)=Bkx(tk),

    which indicates the state 'jump' at the impulse times tk only depends on the current state. When considering the delay effect in impulse, the impulses are updated as follows:

    Δx(tk)=x(tk)x(tk)=Bkx(tkτ),

    where τ is the possible time delay in impulses. It indicates that the state 'jump' at the impulse times tk depends on not only the current state but also the state in history. Such impulses are regarded as a better way to model many practical problem. In fact, many practical systems can be modeled as differential systems with delayed impulses [64,29,10,8,30,16,107,11,52,51,99,50,53]. For example, in [64,29], the stability problem of impulsive systems with stabilizing delayed impulses was studied. More exactly, the asymptotic stability in [29] was investigated for a class of delay-free autonomous systems with linear delayed impulses of the form

    x(t)=x(t)+Bkx((tdk)),  t=tk, kN,

    where {tk} is an impulsive time sequence and dk are impulse input delays. A sufficient condition for asymptotic stability involving the sizes of impulse input delays was derived. [10] studied the problem of exponential stability of nonlinear time-delay systems with more general delayed impulses, which includes the linear delayed impulses of the form

    x(t)=C0kx(t)+C1kx((tdk)),  t=tk, kN,

    as a special case. The results in [10] dealt with both destabilizing delayed impulses and stabilizing delayed impulses, and derived the corresponding Lyapunov-type sufficient conditions for exponential stability. In the application of networked control systems, due to the finite speed of computation, a type of delayed impulses which are called sensor-to-controller delay and controller-to-actuator delay do exist in a working network [107,11]. Authors in [52] studied the delayed impulsive control of nonlinear differential systems, where the impulsive control involves the delayed state of the system for which the delay is state-dependent. [51] focused on stability problem of nonlinear differential systems with impulses involving state-dependent delay based on Lyapunov methods. Some general and applicable results for uniform stability, uniform asymptotic stability and exponential stability of the systems were derived in [51] by using the impulsive control theory and some comparison arguments. It shows how restrictions on the change rates of states and impulses should be imposed to achieve system's stability, in comparison with general impulsive delay differential systems with state-dependent delay in the nonlinearity, or the differential systems with constant delays.

    With the development of impulsive control theory, some recent works have focused on input-to-state stability (ISS) property of time-delay control system under the delayed impulsive control. For example, [116] addressed the ISS and integral input-to-state stability (iISS) of nonlinear systems with distributed delayed impulses. [53] studied the ISS property of nonlinear systems with delayed impulses and external input affecting both the continuous dynamics and the state impulse map. However, it seems that there have been few results that consider the effect of delayed impulses on ISS property for nonlinear systems, which still remains as an important direction in research fields.


    7. Conclusion and future work

    IDS is a very important research area with wide applications. Stability analysis is one of the fundamental problems for IDSs. This paper has overviewed the research area of IDSs with emphasis on the following topics:

    (ⅰ) We described the general IDSs and presented the existence and uniqueness of solutions for IDSs. Moreover, we introduced the effects of impulses on stability for IDSs, which includes impulsive stabilization and impulsive perturbation. Examples were given to illustrate the effects of impulses. In addition, the concept of AII was introduced to characterize how often or how seldom seldom impulses occur.

    (ⅱ) We presented sufficient conditions for IDSs, where the impulses contribute to system dynamics. In this sense, an example was illustrated to show the existence of periodic attractor under impulsive control, where there is originally no periodic attractor.

    (ⅲ) We presented sufficient conditions for IDSs, where the impulsive effects are harmful. In this sense, we illustrated an example to show the effects of impulsive perturbation for periodic solutions. It indicates that the dynamics of time-delay systems can be changed under impulsive perturbations.

    (ⅳ) We introduced the delayed impulses. Some interesting results on stability or ISS properties involving delayed impulse have been presented.

    Although IDSs and their control theory have been developed for many years, there are still some shortcomings and problems to be solved:

    (1) Many results on delayed impulses have been derived. However, most of them only considered the negative effect of time delay which exists in impulses. How to study the positive effect for such time delay is still a difficult problem.

    (2) In recent years, the design and optimization of impulsive controller has been a hot and frontier problem for impulsive control theory, but there has been not much research progress in impulsive systems involving time delays. In particular, how to design an optimal impulsive controller under the constraints of engineering background is a key scientific problem that needs to be solved.

    (3) Since impulsive control strategy usually has simple structure in which only discrete control are needed to achieve the desired performance, event-triggered impulsive control deserves increasing attention and some related control strategies have been proposed [93]. However, all the previous works have focused on the design of event-triggered impulsive control strategies for some specific systems, and there is no unified research method for general nonlinear systems, which leads to that the derived results have limitations in applications.


    Acknowledgments

    This work was supported by National Natural Science Foundation of China (11301308, 61673247, 11601269), and the Research Fund for Distinguished Young Scholars and Excellent Young Scholars of Shandong Province (JQ201719, ZR2016JL024). The paper has not been presented at any conference.




    [1] S. Abbas, Existence of solutions to fractional order ordinary and delay differential equations and applications, Electron. J. Differ. Equ., 2011 (2011), 1–11.
    [2] J. Alzabut, A. G. M. Selvam, R. A. El-Nabulsi, V. Dhakshinamoorthy, M. E. Samei, Asymptotic stability of nonlinear discrete fractional pantograph equations with non-local initial conditions, Symmetry, 13 (2021), 473. https://doi.org/10.3390/sym13030473 doi: 10.3390/sym13030473
    [3] J. A. D. Applelby, I. Győri, D. W. Reynolds, On exact convergence rates for solutions of linear systems of Volterra difference equations, J. Difference Equ. Appl., 12 (2006), 1257–1275. https://doi.org/10.1080/10236190600986594 doi: 10.1080/10236190600986594
    [4] K. Balachandran, S. Kiruthika, J. J. Trujillo, Existence of solutions of nonlinear fractional pantograph equations, Acta Math. Sci., 33 (2013), 712–720. https://doi.org/10.1016/S0252-9602(13)60032-6 doi: 10.1016/S0252-9602(13)60032-6
    [5] M. Benchohra, J. Henderson, S. K. Ntouyas, A. Ouahab, Existence results for fractional order functional differential equations with infinite delay, J. Math. Anal. Appl., 338 (2008), 1340–1350. https://doi.org/10.1016/j.jmaa.2007.06.021 doi: 10.1016/j.jmaa.2007.06.021
    [6] N. D. Cong, H. T. Tuan, Existence, uniqueness, and exponential boundedness of global solutions to delay fractional differential equations, Mediterr. J. Math., 14 (2017), 1–12. https://doi.org/10.1007/s00009-017-0997-4 doi: 10.1007/s00009-017-0997-4
    [7] Y. Jalilian, R. Jalilian, Existence of solution for delay fractional differential equations, Mediterr. J. Math., 10 (2013), 1731–1747. https://doi.org/10.1007/s00009-013-0281-1 doi: 10.1007/s00009-013-0281-1
    [8] B. Ji, H. Liao, Y. Gong, L. Zhang, Adaptive second-order Crank-Nicolson time-stepping schemes for time-fractional molecular beam epitaxial growth models, SIAM J. Sci. Comput., 42 (2020), B738–B760. https://doi.org/10.1137/19M1259675 doi: 10.1137/19M1259675
    [9] E. Kaslik, S. Sivasundaram, Analytical and numerical methods for the stability analysis of linear fractional delay differential equations, J. Comput. Appl. Math., 236 (2012), 4027–4041. https://doi.org/10.1016/j.cam.2012.03.010 doi: 10.1016/j.cam.2012.03.010
    [10] K. Krol, Asymptotic properties of fractional delay differential equations, Appl. Math. Comput., 218 (2011), 1515–1532. https://doi.org/10.1016/j.amc.2011.04.059 doi: 10.1016/j.amc.2011.04.059
    [11] V. Lakshmikantham, Theory of fractional functional differential equations, Nonlinear Anal., 69 (2008), 3337–3343. https://doi.org/10.1016/j.na.2007.09.025 doi: 10.1016/j.na.2007.09.025
    [12] D. Li, C. Zhang, Long time numerical behaviors of fractional pantograph equations, Math. Comput. Simul., 172 (2020), 244–257. https://doi.org/10.1016/j.matcom.2019.12.004 doi: 10.1016/j.matcom.2019.12.004
    [13] S. Li, A review of theoretical and numerical analysis for nonlinear stiff Volterra functional differential equations, Front. Math. China, 4 (2009), 23–48. https://doi.org/10.1007/s11464-009-0003-y doi: 10.1007/s11464-009-0003-y
    [14] S. Li, High order contractive Runge-Kutta methods for Volterra functional differential equations, SIAM J. Numer. Anal., 47 (2010), 4290–4325. https://doi.org/10.1137/080741148 doi: 10.1137/080741148
    [15] S. Li, Numerical Analysis for Stiff Ordinary and Functional Differential Equations, Xiangtan: Xiangtan University Press, 2018.
    [16] C. Lubich, On the stability of linear multistep methods for Volterra convolution equations, IMA J. Numer. Anal., 3 (1983), 439–465. https://doi.org/10.1093/imanum/3.4.439 doi: 10.1093/imanum/3.4.439
    [17] K. Mustapha, An L1 approximation for a fractional reaction-diffusion equation, a second-order error analysis over time-graded meshes, SIAM J. Numer. Anal., 58 (2020), 1319–1338. https://doi.org/10.1137/19M1260475 doi: 10.1137/19M1260475
    [18] J. Robinson, Infinite-dimensional Dynamical Systems, Cambridge: Cambridge University Press, 2001.
    [19] J. Shen, F. Zeng, M. Stynes, Second-order error analysis of the averaged L1 scheme ¯L1 for time-fractional initial-value and subdiffusion problems, Sci. China Math., 2022. https://doi.org/10.1007/s11425-022-2078-4
    [20] A. Stuart, A. Humphries, Dynamical Systems and Numerical Analysis, Cambridge: Cambridge Unversity Press, 1996.
    [21] M. Stynes, A survey of the L1 scheme in the discretisation of time-fractional problems, Numer. Math. Theor. Meth. Appl., 15 (2022), 1173–1192. https://doi.org/10.4208/nmtma.OA-2022-0009s doi: 10.4208/nmtma.OA-2022-0009s
    [22] R. Temam, Infinite-dimensional Dynamical System in Mechanics and Physics, Berlin: Springer, 1997.
    [23] D. Wang, A. Xiao, Dissipativity and contractivity for fractional-order systems, Nonlinear Dyn., 80 (2015), 287–294. https://doi.org/10.1007/s11071-014-1868-1 doi: 10.1007/s11071-014-1868-1
    [24] D. Wang, A. Xiao, H. Liu, Dissipativity and stability analysis for fractional functional differential equations, Fract. Calc. Appl. Anal., 18 (2015), 1399–1422. https://doi.org/10.1515/fca-2015-0081 doi: 10.1515/fca-2015-0081
    [25] D. Wang, A. Xiao, S. Sun, Asymptotic behavior of solutions to time fractional neutral functional differential equations, J. Comput. Appl. Math., 382 (2021), 113086. https://doi.org/10.1016/j.cam.2020.113086 doi: 10.1016/j.cam.2020.113086
    [26] D. Wang, A. Xiao, J. Zou, Long-time behavior of numerical solutions to nonlinear fractional ODEs, ESAIM Math. Model. Numer. Anal., 54 (2020), 335–358. https://doi.org/10.1051/m2an/2019055 doi: 10.1051/m2an/2019055
    [27] D. Wang, J. Zou, Dissipativity and contractivity analysis for fractional functional differential equations and their numerical approximations, SIAM J. Numer. Anal., 57 (2019), 1445–1470. https://doi.org/10.1137/17M1121354 doi: 10.1137/17M1121354
    [28] F. Wang, D. Chen, X. Zhang, Y. Wu, The existence and uniqueness theorem of the solution to a class of nonlinear fractional order system with time delay, Appl. Math. Lett., 53 (2016), 45–51. https://doi.org/10.1016/j.aml.2015.10.001 doi: 10.1016/j.aml.2015.10.001
    [29] W. Wang, C. Zhang, Dissipativity of variable-stepsize Runge- Kutta methods for nonlinear functional differential equations with application to Nicholson's blowflies models, Commun. Nonlinear Sci. Numer. Simulat., 97 (2021), 105723. https://doi.org/10.1016/j.cnsns.2021.105723 doi: 10.1016/j.cnsns.2021.105723
    [30] Z. Wang, X. Huang, G. Shi, Analysis of nonlinear dynamics and chaos in a fractional order financial system with time delay, Comput. Math. Appl., 62 (2011), 1531–1539. http://dx.doi.org/10.1016/j.camwa.2011.04.057 doi: 10.1016/j.camwa.2011.04.057
    [31] L. Wen, S. Li, Dissipativity of Volterra functional differential equations, J. Math. Anal. Appl., 324 (2006), 696–706. https://doi.org/10.1016/j.jmaa.2005.12.031 doi: 10.1016/j.jmaa.2005.12.031
    [32] L. Wen, Y. Yu, S. Li, Dissipativity of Runge- Kutta methods for Volterra functional differential equations, Appl. Numer. Math., 61 (2011), 368–381. https://doi.org/10.1016/j.apnum.2010.11.002 doi: 10.1016/j.apnum.2010.11.002
    [33] L. Wen, Y. Yu, W. Wang, Generalized Halanay inequalities for dissipativity of Volterra functional differential equations, J. Math. Anal. Appl., 347 (2008), 169–178. https://doi.org/10.1016/j.jmaa.2008.05.007 doi: 10.1016/j.jmaa.2008.05.007
    [34] Y. Yan, C. Kou, Stability analysis for a fractional differential model of HIV infection of CD4+ T-cells with time delay, Math. Comput. Simul., 82 (2012), 1572–1585. http://dx.doi.org/10.1016/j.matcom.2012.01.004 doi: 10.1016/j.matcom.2012.01.004
    [35] Y. Zhou, M. Stynes, Optimal convergence rates in time-fractional discretisations: the L1, ¯L1 and Alikhanov schemes, East Asian J. Appl. Math., 12 (2022), 503–520. https://doi.org/10.4208/eajam.290621.220921 doi: 10.4208/eajam.290621.220921
  • Reader Comments
  • © 2023 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(1331) PDF downloads(37) Cited by(0)

Article outline

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog