Research article Special Issues

Barycentric rational collocation method for semi-infinite domain problems

  • Received: 25 December 2022 Revised: 01 February 2023 Accepted: 02 February 2023 Published: 07 February 2023
  • MSC : 65D32, 65D30, 65R20

  • The barycentric rational collocation method for solving semi-infinite domain problems is presented. Following the barycentric interpolation method of rational polynomial and Chebyshev polynomial, matrix equation is obtained from discrete semi-infinite domain problem. Truncation method and transformation method are presented to solve linear and nonlinear differential equation defined on the semi-infinite domain problems. At last, three numerical examples are presented to valid our theoretical analysis.

    Citation: Jin Li. Barycentric rational collocation method for semi-infinite domain problems[J]. AIMS Mathematics, 2023, 8(4): 8756-8771. doi: 10.3934/math.2023439

    Related Papers:

    [1] Mohammed Aldandani, Omar Naifar, Abdellatif Ben Makhlouf . Practical stability for nonlinear systems with generalized conformable derivative. AIMS Mathematics, 2023, 8(7): 15618-15632. doi: 10.3934/math.2023797
    [2] Yuanlin Ding . Existence and stability analysis of solutions for periodic conformable differential systems with non-instantaneous impulses. AIMS Mathematics, 2025, 10(2): 4040-4066. doi: 10.3934/math.2025188
    [3] Muhammad Bilal, Javed Iqbal, Ikram Ullah, Aditi Sharma, Hasim Khan, Sunil Kumar Sharma . Novel optical soliton solutions for the generalized integrable (2+1)- dimensional nonlinear Schrödinger system with conformable derivative. AIMS Mathematics, 2025, 10(5): 10943-10975. doi: 10.3934/math.2025497
    [4] Xuelian Jin . Exponential stability analysis and control design for nonlinear system with time-varying delay. AIMS Mathematics, 2021, 6(1): 102-113. doi: 10.3934/math.2021008
    [5] Aisha F. Fareed, Emad A. Mohamed, Mokhtar Aly, Mourad S. Semary . A novel numerical method for stochastic conformable fractional differential systems. AIMS Mathematics, 2025, 10(3): 7509-7525. doi: 10.3934/math.2025345
    [6] Awais Younus, Khizra Bukhsh, Manar A. Alqudah, Thabet Abdeljawad . Generalized exponential function and initial value problem for conformable dynamic equations. AIMS Mathematics, 2022, 7(7): 12050-12076. doi: 10.3934/math.2022670
    [7] Tan Zhang, Pianpian Yan . Asymmetric integral barrier function-based tracking control of constrained robots. AIMS Mathematics, 2024, 9(1): 319-339. doi: 10.3934/math.2024019
    [8] Noha M. Kamel, Hamdy M. Ahmed, Wafaa B. Rabie . Solitons unveilings and modulation instability analysis for sixth-order coupled nonlinear Schrödinger equations in fiber bragg gratings. AIMS Mathematics, 2025, 10(3): 6952-6980. doi: 10.3934/math.2025318
    [9] Narongrit Kaewbanjak, Watcharin Chartbupapan, Kamsing Nonlaopon, Kanit Mukdasai . The Lyapunov-Razumikhin theorem for the conformable fractional system with delay. AIMS Mathematics, 2022, 7(3): 4795-4802. doi: 10.3934/math.2022267
    [10] Erfeng Xu, Wenxing Xiao, Yonggang Chen . Local stabilization for a hyperchaotic finance system via time-delayed feedback based on discrete-time observations. AIMS Mathematics, 2023, 8(9): 20510-20529. doi: 10.3934/math.20231045
  • The barycentric rational collocation method for solving semi-infinite domain problems is presented. Following the barycentric interpolation method of rational polynomial and Chebyshev polynomial, matrix equation is obtained from discrete semi-infinite domain problem. Truncation method and transformation method are presented to solve linear and nonlinear differential equation defined on the semi-infinite domain problems. At last, three numerical examples are presented to valid our theoretical analysis.



    Since the second half of the 20th century, fractional calculus found several useful applications, as well as extensions of classical results, through the Riemann-Liouville integral and the Riemann-Liouville and Caputo derivatives; however, new definitions of fractional operators have been recently proposed, which aim to improve the modelling of many systems. Thus, several operators are now being considered in the development of the theory and the applications of fractional calculus [1,2,3,4,5,6,7,8,9,10]; some recent reviews and analysis of the most used definitions are found in [11,12].

    However, the main disadvantage of the original fractional operators is that they fail to accomplish with some rules existing in integer-order calculus, e.g. the Leibniz product rule, the quotient rule, and the chain rule. For these reasons, R. Khalil et al. proposed a new definition, called the conformable fractional derivative, which fulfilled these properties along with some other classical calculus results [13]. After that, T. Abdeljawad extended the study of this derivative to other properties in calculus and linear systems [14]. Thenceforth, several works have developed the theory and applications of this operator [15,16,17,18,19,20,21,22,23].

    On the other hand, U. Katugampola proposed other fractional operators, which also accomplish with the classical results [24]. These operators are also being studied and applied [25,26,27,28,29,30].

    Furthermore, based on the tools proposed by Khalil et al. and Katugampola, Akkurt et al. proposed the so-called Generalized fractional derivative [31], which generalizes both operators and also satisfies the classical properties and results from integer-order calculus. This derivative is general because of the freedom to choose its kernel; however, a better nomenclature for this operator may be "generalized conformable derivative", given that it encompasses the conformable-type derivatives formerly proposed, and thus in this work we will refer to it in this way. This operator has also been considered for further studies and applications [32,33,34,35].

    Due to the level of generality and its accomplishment of the classical calculus results, the idea of developing theoretical and applied results using the generalized conformable derivative is appealing, because this would encompass results using the conformable operators. In the work [36], another generalized conformable operator [37] has been assessed, and it was observed that the results obtained with this operator outperformed the ones obtained with integer and fractional operators, presenting even finite-time stability.

    On another topic, there exist several definitions and configurations of observers, which are useful for different applications. Consider e.g. the intermediate estimator [38], which is useful for observing states and faults simultaneously in systems with Lipschitzian nonlinearities without requiring to fulfill the observer matching condition from other adaptive observers. As another example, the interval observer in [39] is useful for estimation in discrete-time linear systems with unknown but bounded disturbance, which is less restrictive than other similar approaches. Therefore, in order to assess the operators studied in this paper, we will consider the design of a nonlinear quadratic regulator (NQR) [40] and a high-gain observer [41]. The NQR provides robust stability with a minimized energy-like performance index, besides being computationally efficient. Moreover, the high-gain observer is robust in the presence of external disturbances or uncertainties (like noise in the measurement), has a rapid exponential decay of the estimation error, and the speed of convergence increases by augmenting the gain of the observer. Besides, both observers are simple to implement and practical for many applications.

    As it was mentioned, the generalized conformable derivative proposed by Akkurt et al. generalizes other operators and satisfies the classical calculus rules, thus, it is of great interest to develop results in theory and applications, which would enclose existing or possible results using conformable-type derivatives. For example, it may be very useful to obtain generalized results in control theory, such as modelling, control algorithms, observer design, stability proofs, among others. Therefore, in this work we use the generalized conformable derivative to design a NQR-based estimator and a high-gain observer, for a class of nonlinear systems with the generalized conformable derivative in their dynamics. Then, with the aid of a generalized exponential function, some Lyapunov-like generalized stability theorems are proven. These stability results are considered in the design of the observers, and later some numerical simulations are carried out in two systems. Finally, for evaluating the performance of the operator studied, the results obtained are compared with those obtained with the integer and the fractional-order operators.

    The main contributions of this work are the following:

    ● A generalized conformable exponential function, which implies exponential stability for the class of systems studied.

    ● Some Lyapunov-like theorems to verify the stability of a class of nonlinear systems using the derivative in question, generalizing the work presented in [42].

    ● The design of stable nonlinear observers for these type of systems, which are validated by numerical simulations applied to the generalized models of the simple pendulum and the Van der Pol chaotic oscillator.

    The paper is structured as follows: In Section 2, the generalized conformable derivative is presented, with some of its properties and existing results. In Section 3, the generalized conformable exponential function is presented, in order to prove some Lyapunov-like theorems for systems with the generalized conformable derivative. Section 4 shows the design of nonlinear observers for systems with this derivative, which are proven to be exponentially stable in a generalized sense. In Section 5, numerical simulations are carried out in the models mentioned, in order to evaluate the performance of the observers; moreover, they are compared with their integer and fractional-order versions. Finally, some conclusions, potential applications and future work are given in Section 6.

    In this section, the definition of the generalized conformable derivative is presented, along with some of its properties. For the sake of comparison, consider first the original definition proposed by Khalil et al.

    Definition 1 ([13]) Given a function f:[0,)R, the conformable derivative of f of order α is defined by

    Tα(f)(t)=limϵ0f(t+ϵt1α)f(t)ϵ,

    t>0,α(0,1). If f is α-differentiable in some (0,a),a>0, and limt0+f(α)(t) exists, then define

    f(α)(0)=limt0+f(α)(t).

    Now, we present the definition of the generalized operator proposed by Akkurt et al.

    Definition 2 ([31]) Let k:[a,b]R be a strictly increasing, continuously differentiable nonnegative map such that k(t),k(t)>0, whenever t>a. Given a function f:[a,b]R and α(0,1) a real, we say that the generalized conformable derivative of f of order α is defined by

    Dαf(t):=limϵ0f(tk(t)+k(t)eϵkα(t)k(t))f(t)ϵ,

    given that the limit exists. We say then that f is α-differentiable.

    If f is α-differentiable in (0,a) and limt0+f(α)(t) exists, then define

    Dαf(0)=limt0+Dαf(t).

    Henceforth, the operator Dα will always refer to the generalized conformable derivative (GCD), and the α-differentiability property will refer to this operator.

    Theorem 1. ([31]) Let f:[a,b]R be a continuously differentiable function and t>a. Then, f is α-differentiable at t and

    Dαf(t)=k1α(t)k(t)f(t).

    Also, if f is continuous at t=a, then

    Dαf(a)=k1α(a)k(a)f(a).

    Remark 1. There exist α-differentiable functions which are not differentiable in the usual sense, and thus the order α and the kernel k(t) may help to obtain different results between the GCD and the integer-order derivative.

    Theorem 2. ([31]) If a function f:[a,b]R is α-differentiable at a>0,α(0,1], then f is continuous at a.

    Theorem 3. ([31]) Let α(0,1] and f,g be α-differentiable at a point t>0. Then

    1. Dα[cf(t)+dg(t)]=cDαf(t)+dDαg(t),c,dR (linearity).

    2. Dα[tn]=k1α(t)k(t)ntn1,nR.

    3. Dαc=0, for all constant functions f(t)=c.

    4. Dα[f(t)g(t)]=f(t)Dαg(t)+g(t)Dαf(t) (product rule).

    5. Dα[f(t)g(t)]=f(t)Dαg(t)g(t)Dαf(t)[g(t)]2 (quotient rule).

    6. Dα[f(t)g(t)]=k1α(t)k(t)f(g(t))g(t) (chain rule).

    Moreover, consider the following result.

    Lemma 1. Let f:[a,)R be α-differentiable on (a,). If Dαf(t)0 (respectively 0) t(a,), then f is an increasing function (respectively decreasing).

    Proof. This proof can be derived from the extension of the Mean Value Theorem for the operator in question [31].

    Remark 2. Let x:[a,)Rn be α-differentiable on (a,). Let P be a symmetric positive definite matrix. Then

    Dα[x(t)TPx(t)]=2x(t)TPDαx(t),t>a.

    On the other hand, consider the generalized conformable integral.

    Definition 3. ([31]) Let ta0. Also, let f be a function defined on (a,t] and αR. Let k:[a,b]R be a continuous nonnegative map such that k(t),k(t)0. Then, the α-generalized conformable integral of f is defined by

    Iα(f)(t)=tak(x)f(x)(k(x))1αdx,

    if the Riemann improper integral exists.

    Consider the following class of nonlinear systems

    Dαx=f(t,x),t>t0,x(t0)=x0, (3.1)

    where xRn,f:R+×RnRn is a given nonlinear function satisfying f(t,0)=0,t0, and α(0,1].

    Throughout this paper, the function |||| represents the Euclidean norm.

    Definition 4. The origin of system (3.1) is said to be

    i) stable, if for every ε>0 and t0R+,δ(ε,t0) such that for any x0Rn,||x0||<δ||x(t)||<ε,tt0.

    ii) attractive, if for any t00,c(t0)>0 such that for any x0Rn,||x0||<climtx(t)=0.

    iii) asymptotically stable, if it is stable and attractive.

    iv) globally asymptotically stable, if it is asymptotically stable for any x0Rn.

    The next function will be the main instrument used in the stability proofs for systems involving the operator considered.

    Definition 5. The generalized conformable exponential function Eαk(), with parameters α,γ and t0, and which depends on the kernel k(t), is defined as

    Eαk(γ,t,t0)=exp{γα[kα(t)kα(t0)]},tt0, (3.2)

    where α(0,1],γR and k(t) is as stated in Definition 2.

    Remark 3. It is not difficult to verify that Eαk(γ,t,t0) is α-differentiable, and that

    Dα[Eαk(γ,t,t0)]=γEαk(γ,t,t0).

    Now, the notion of generalized conformable exponential stability is introduced.

    Definition 6. The origin of system (3.1) is said to be generalized conformable exponentially stable (GCES) if

    ||x||C||x0||Eαk(γ,t,t0),

    with t>t0 and C,γ>0.

    Lemma 2. Let g:[t0,)R+ be an α-differentiable function on (t0,) such that

    Dαg(t)γg(t),

    where γ>0 and α(0,1]. Then

    g(t)g(t0)Eαk(γ,t,t0).

    Proof. Let h(t)=g(t)Eαk(γ,t,t0). From Theorem 3, we have

    Dαh(t)=g(t)Dα[Eαk(γ,t,t0)]+Dα[g(t)]Eαk(γ,t,t0)γg(t)Eαk(γ,t,t0)γg(t)Eαk(γ,t,t0).

    Since Dαh(t)0, from Lemma 1 h(t) is a decreasing function. Hence, h(t)h(t0), which gives the result.

    The results that follow may be considered as Lyapunov-type stability theorems for the class of systems considered in this work.

    Theorem 4. Let x=0 be an equilibrium point for system (3.1). Let the function V:R+×RnR be α-differentiable with respect to t and x, and consider ai(i=1,2,3) as arbitrary positive constants. If the following conditions are satisfied:

    (i) a1||x||2V(t,x)a2||x||2,

    (ii) DαV(t,x)a3||x||2,

    then the origin of system (3.1) is GCES.

    Proof. From conditions (i) and (ii) we have

    DαV(t,x)a3a2V(t,x).

    Applying Lemma 2 to the last inequality, we get

    V(t,x)V(t0,x0)Eαk(a3a2,t,t0).

    Also from (i), we have

    a1||x||2V(t0,x0)Eαk(a3a2,t,t0)a2Eαk(a3a2,t,t0)||x0||2.

    Thus

    ||x||2(a2a1)Eαk(a3a2,t,t0)||x0||2,

    and

    ||x||C||x0||Eαk(γ,t,t0),

    with C=(a2/a1)1/2 and γ=a3/(2a2). Therefore, the origin of system (3.1) is GCES.

    Remark 4. Note that GCES implies the classical exponential stability, taking α=1 and k(t)=t.

    Definition 7. A continuous function k:R+R+ is said to belong to class K if it is strictly increasing and k(0)=0. It belongs to class K if limt+k(t)=+.

    Theorem 5. Let x=0 be an equilibrium point for system (3.1). Let the function V:R+×RnR be α-differentiable with respect to t and x, and consider ki(i=1,2,3) as functions of class K. If the following conditions are satisfied:

    (i) k1(||x||)V(t,x)k2(||x||),

    (ii) DαV(t,x)k3(||x||),

    then the origin of system (3.1) is locally asymptotically stable.

    Furthermore, if kiK(i=1,2,3), then the origin of system (3.1) is globally asymptotically stable.

    Proof. First, from condition (ii) and the fact that k(t) is strictly increasing and nonnegative, one has that DαV(t,x)0 for almost all tt0. Thus, from Lemma 1 V(t,x) is monotonically decreasing, but given that V(t,x)0 from (i), this implies the existence of some real number L=limtV(t,x)0.

    Moreover, condition (i) is equivalent to k12(V(t,x))||x||k11(V(t,x)), since k1 and k2 are monotonically increasing functions. By the same argument, one has that k3(||x||)k3(k12(V(t,x))). Therefore, condition (ii) implies that:

    DαV(t,x)k3k12(V(t,x)).

    Now, suppose that L>0. Furthermore, define

    λ=k3k12(L)V(t0,x0);

    it is clear that λ>0 when L>0. Then,

    DαV(t,x)k3k12(V(t,x))k3k12(L)sinceLV(t,x),tt0k3k12(L)V(t0,x0)V(t0,x0)=λV(t0,x0).

    Finally, by considering that V(t,x)V(t0,x0) (from Lemma 1), one has that

    DαV(t,x)λV(t,x),

    hence, from Lemma 2 one has

    V(t,x)V(t0,x0)Eαk(λ,t,t0),

    and consequently, L=limtV(t,x)=0, but it was supposed that L>0, which is a contradiction. Thus, L=0 and from (i), one has that limtx(t)=0 and therefore, the origin of system (3.1) is locally asymptotically stable.

    Considering the case where kiK(i=1,2,3), it is straightforward that the origin of system (3.1) is globally asymptotically stable.

    Remark 5. Note that Theorems 4 and 5 involve the generalized conformable derivative. Then, these results may serve as meta-theorems for proving generalized exponential and asymptotic stability for a more general class of dynamical systems, which consider a larger family of differential operators.

    In this section, a pair of observers are proposed for a class of nonlinear systems with the GCD in their dynamics. Moreover, the stability results presented in the past section are used to prove that the estimation error of both observers is GCES.

    Consider the following class of systems

    Dαx=f(x,u),x(0)=x0,y=Cx, (4.1)

    where 0<α<1, xRn is the state vector, uRm is the input vector, yR is the output and f:Rn×RmRn is locally Lipschitz in x and uniformly bounded in u.

    Former system may be rewritten in the following canonical form [43]

    Dαx=Ax+Υ(x,u),x(0)=x0,y=Cx, (4.2)

    where A is an upper shift matrix (A:RnRn,Ai,j=δi+1,j, with δi,j the Kronecker delta), C=[1000], the pair (A,C) is observable, and Υ(x, u) satisfies

    ||Υ(x,u)Υ(ˆx,u)||φ||xˆx||, (4.3)

    in a region D. The observers will be designed for the system expressed in this canonical equivalent.

    First, consider the following generalized conformable estimator, based on the design of a nonlinear quadratic regulator (NQR)

    Dαx=Aˆx+Υ(ˆx,u)+mi=1Ki(yCˆx)2i1, (4.4)

    where ˆx,KiRn,1im. For this system, consider the following result.

    Lemma 3. ([44]) Given a stable n×m matrix ˆA and γ>0, there exists a positive definite, symmetric matrix P such that

    ˆATP+PˆA+γ2PP+I<0, (4.5)

    if and only if there exists another positive definite, symmetric matrix ˆP such that

    ˆAˆP+ˆPˆAT+γ2ˆPˆP+I<0.

    We have the following LMI, equivalent to (4.5)

    [ˆATPPˆAIγPγPI]>0. (4.6)

    Note that, for some ε>0,

    ˆATP+PˆA+γ2PP+I+εI=0. (4.7)

    Now, let υ=Υ(x,u)Υ(ˆx,u). Considering the Lipschitz condition (4.3), the estimation error e=xˆx and the solution P of (4.6), we have the following inequality [44]

    2eTPυφ2eTPPe+eTe. (4.8)

    Remark 6. A observer is said to be generalized conformable exponentially stable if its estimation error e is GCES.

    Theorem 6. Consider system (4.2) with (A,C) observable. If ˉA=AK1C is a stable matrix and KiC>0, then the estimator (4.4) is a generalized conformable exponentially stable observer for system (4.1).

    Proof. From (4.2) and (4.4), the fractional dynamics of the estimation error is

    Dαe=ˉAe+Fmi=2Ki(Ce)2i1.

    Since ˉA is stable and φ>0, from Lemma 3 there exists some P>0. Consider V=||e||2P=eTPe a candidate Lyapunov function that satisfies the Rayleigh-Ritz inequality

    λmin  (P)||e||2Vλmax  (P)||e||2. (4.9)

    Moreover, from Remark 2, (4.7) and (4.8), we have

    DαV=2eTPDαe=2eTP[ˉAe+υmi=2Ki(Ce)2i1]eT[ˆATP+PˆA+γ2PP+I]e2mi=2(Ce)2i2eTPKiCeε||e||2mi=2(Ce)2i2eTPKiCe.

    Given that, KiC is a positive real number for any i, one has that PKiC>0, and consequently,

    DαVε||e||2. (4.10)

    Finally, from Theorem 4, (4.9) and (4.10) it follows that e=0 is GCES and we have

    ||e||C||e0||Eαk(γ,t,t0),

    with C=λmax  (P)λmin  (P),γ=ε2λmax  (P)ande(t0)=e0.

    Now, consider the following generalized conformable high-gain observer (HGO)

    Dαˆx=Aˆx+Υ(ˆx,u)+F1CT(yCˆx), (4.11)

    where F=limtFθ(t), with Fθ(t) a positive definite matrix solution of

    DαFθ(t)=θFθ(t)ATFθ(t)Fθ(t)A+CTC,F0=Fθ(0).

    Remark 7. Given that F is a constant matrix, DαF=0 and thus it can be obtained from the following algebraic equation

    θFATFFA+CTC=0. (4.12)

    The coefficients of F are given by

    (F)i,j=αi,jθi+j1,

    where αi,j is a symmetric positive definite matrix independent of θ [41,45].

    Theorem 7. The high-gain observer (4.11) is a generalized conformable exponentially stable observer for system (4.1).

    Proof. Consider again the observation error e=xˆx. From (4.2) and (4.11), the fractional dynamics of the error is

    Dαe=(A+F1CTC)e+υ,

    with υ=Υ(x,u)Υ(ˆx,u). Consider V=||e||2F=eTFe a candidate Lyapunov function that satisfies the Rayleigh-Ritz inequality

    λmin  (F)||e||2Vλmax  (F)||e||2. (4.13)

    Moreover, from Remark 2 and (4.12) we have

    DαV=2eTFDαe=2eTF[(A+F1CTC)e+υ]θeTFe+2eTFυ.

    Thus, considering that eTFυφλmax  (F)||e||2, one has that

    DαV(θλmin  (F)2φλmax  (F))||e||2.

    Finally, for θ>2φλmax  (F)/λmin  (F), the system is asymptotically stable. In addition, according to Theorem 4, it follows that e=0 is GCES and we have

    ||e||C||e0||Eαk(γ,t,t0),

    with C=λmax  (F)λmin  (F),γ=θλmin  (F)2φλmax  (F)2λmax  (F)ande(t0)=e0.

    Remark 8. Note that the results of Theorems 4.1 and 4.2 involve the generalized conformable derivative, and given that this operator encompasses the conformable type derivatives, these results may serve as meta-theorems for designing exponentially stable estimators for the class of systems that consider this family of operators.

    In this part, numerical simulations of the observation scheme developed in the former section are presented. In order to assess the performance of the observers, the simulations were performed on two models. Also, these results will be compared to the ones obtained with the integer [40,41] and fractional-order [46,47] versions of the same observers. The simulations were performed using Simulink® from MATLAB®.

    For the fractional-order dynamics, the Caputo operator was considered, given to the need of having integer-order initial conditions; to implement this operator, the ninteger fractional derivative block from D. Valério was used. Furthermore, in order to implement the GCD, the function k(t)=2t2+2t+2 was used for both systems and observers.

    Consider the simple pendulum, which has been extended and studied in its fractional version [48,49,50,51]. The model of the simple pendulum with GCD is the following:

    Dαx1=x2,Dαx2=gLsin(x1),y=x1,

    where x1=θ (angular position), x2=ω (angular velocity), g=9.81m/s2 and L=1m. Rewrite the system to its canonical form:

    Dαx=Ax+Υ(x),x(0)=x0,y=Cx, (5.1)

    with A=[0100],Υ(x)=[0gLsin(x1)] with Lipschitz constant φ=gL, and C=[10].

    The fractional dynamics in the model of the pendulum, as has been seen with other mechanical systems, adds the effect of damping on the trajectories [50,51], which represents an additional parameter which value depends on the fractional order. In the case of the dynamics with GCD, due to the freedom of choice of the order α and the kernel k(t), the resulting model can be seen as a pendulum with some time-varying parameters (a nonautonomous system). From the graphs of the states shown in this section it can be observed, for example, a variation in the frequency of oscillation.

    From (4.4), taking m=3, the NQR-based estimator for system (5.1) is

    Dαˆx=Aˆx+Υ(ˆx)+K1C(xˆx)+K2[C(xˆx)]3+K3[C(xˆx)]5,y=ˆx,

    where Ki=[Ki1Ki2]T.

    For the simulation α=0.98 was selected, the gains were chosen as K1=[4.00763.1305]T, K2=[4.9054.905]T, K3=[52.4525]T, and the initial conditions are x1(0)=π/2, x2(0)=0, ˆx1(0)=π, ˆx2(0)=0.

    Figures 1 and 2 show the estimations of the states by the NQR-based observer with the GCD. Moreover, Figures 3 and 4 show the behaviour of the estimation errors for each state in the integer, fractional, and generalized conformable cases, for the same observer gains and initial conditions. It can be seen that, although the fractional observer performs slightly better that the integer-order one, the estimation errors obtained with the GCD converge to zero faster than the other.

    Figure 1.  Comparison between x1 and ˆx1 obtained from the NQR-based estimator with the GCD for the simple pendulum.
    Figure 2.  Comparison between x2 and ˆx2 obtained from the NQR-based estimator with the GCD for the simple pendulum.
    Figure 3.  Comparison between the e1 obtained from the NQR-based estimator with integer, fractional and GC derivatives for the simple pendulum.
    Figure 4.  Comparison between the e2 obtained from the NQR-based estimator with integer, fractional and GC derivatives for the simple pendulum.

    Now, consider the high-gain observer (4.11). For the pendulum system, the observer takes the structure

    Dαˆx1=ˆx2+2θ(x1ˆx1),Dαˆx2=gLsin(ˆx1)+θ2(x1ˆx1),y=ˆx.

    Simulations for this observer were performed using the same initial conditions and θ=25. Figures 5 and 6 show the estimations of the states with the GCD, and Figures 7 and 8 show the behaviour of the estimation errors for each state in the integer, fractional, and generalized conformable cases. In the three versions the same gain and initial conditions have been used. It can be seen, as in the former case, that the estimation errors obtained with the GCD converge to zero faster than with the integer and fractional-order versions, which again show a similar performance with a slight advantage for the fractional case.

    Figure 5.  Comparison between x1 and ˆx1 obtained from the high-gain observer with the GCD for the simple pendulum.
    Figure 6.  Comparison between x2 and ˆx2 obtained from the high-gain observer with the GCD for the simple pendulum.
    Figure 7.  Comparison between the e1 obtained from the high-gain observer with integer, fractional and GC derivatives for the simple pendulum.
    Figure 8.  Comparison between the e2 obtained from the high-gain observer with integer, fractional and GC derivatives for the simple pendulum.

    The subject of chaotic oscillators is a research area of great interest due to their theoretical results and applications. Particularly, the Van der Pol oscillator was proposed to study oscillations in vacuum tube circuits [52]. This model has also been extended and studied in its fractional-order counterpart [53,54,55,56]. Using the GCD, the model of the system is

    Dαx1=x2,Dαx2=x1ε(x211)x2,y=x1,

    where ε is the control parameter. The system is rewritten in its canonical form:

    Dαx=Ax+Υ(x),x(0)=x0,y=Cx, (5.2)

    with A=[0100],Υ(x)=[0x1ε(x211)x2] with Lipschitz constant φ=max{1+2ε|x1||x2|+ε+ε|x21|}, and C=[10].

    Hence, the Van der Pol system will perform as a master system, to which the observers designed for it will serve as slaves, looking to synchronize with it.

    From (4.4) with m=3, the NQR-based estimator for system (5.2) is

    Dαˆx=Aˆx+Υ(ˆx)+K1C(xˆx)+K2[C(xˆx)]3+K3[C(xˆx)]5,y=ˆx,

    where Ki=[Ki1Ki2]T.

    For the simulation, the gains were selected as K1=[2.30941.5166]T, K2=[1.151.15]T, K3=[30.575]T. The parameters are α=0.9, ε=0.1 and the initial conditions were chosen as x1(0)=0.25, x2(0)=1.2, ˆx1(0)=3, ˆx2(0)=1.8. With these values, the Lipschitz constant is set as φ=2.3.

    Figures 9 and 10 show the estimations of the states by the NQR-based observer with the GCD. Moreover, Figures 11 and 12 show the behaviour of the estimation errors for each state in the integer, fractional, and generalized conformable cases, for the same observer gains and initial conditions. It can be seen that although the fractional observer performs slightly better that the integer-order one, the estimation errors obtained with the GCD converge to zero faster than the other.

    Figure 9.  Comparison between x1 and ˆx1 obtained from the NQR-based estimator with the GCD for the Van der Pol oscillator.
    Figure 10.  Comparison between x2 and ˆx2 obtained from the NQR-based estimator with the GCD for the Van der Pol oscillator.
    Figure 11.  Comparison between the e1 obtained from the NQR-based estimator with integer, fractional and GC derivatives for the Van der Pol oscillator.
    Figure 12.  Comparison between the e2 obtained from the NQR-based estimator with integer, fractional and GC derivatives for the Van der Pol oscillator.

    Moreover, Figure 13 shows the phase portrait obtained with the estimated states ˆx1 and ˆx2 obtained from the NQR-based estimator with the GCD.

    Figure 13.  Phase portrait of the Van der Pol oscillator with the estimated states obtained from the NQR-based estimator with GCD.

    Now, consider the high-gain observer (4.11). For the Van der Pol oscillator, the observer takes the structure

    Dαˆx1=ˆx2+2θ(x1ˆx1),Dαˆx2=ˆx1ε(ˆx211)ˆx2+θ2(x1ˆx1),y=ˆx.

    Simulations for this observer were performed using the same initial conditions and θ=1. Figures 14 and 15 show the estimations of the states with the GCD, and Figures 16 and 17 show the behaviour of the estimation errors for each state in the integer, fractional, and generalized conformable cases. In the three versions the same gain and initial conditions have been used. It can be seen, as in the former case, that the estimation errors obtained with the GCD converge to zero faster than with the integer and fractional-order versions, which again show a similar performance with a slight advantage to the fractional case.

    Figure 14.  Comparison between x1 and ˆx1 obtained from the high-gain observer with the GCD for the Van der Pol oscillator.
    Figure 15.  Comparison between x2 and ˆx2 obtained from the high-gain observer with the GCD for the Van der Pol oscillator.
    Figure 16.  Comparison between the e1 obtained from the high-gain observer with integer, fractional and GC derivatives for the Van der Pol oscillator.
    Figure 17.  Comparison between the e2 obtained from the high-gain observer with integer, fractional and GC derivatives for the Van der Pol oscillator.

    Finally, Figure 18 shows the phase portrait obtained with the estimated states ˆx1 and ˆx2 obtained from the high-gain observer with the GCD.

    Figure 18.  Phase portrait of the Van der Pol oscillator with the estimated states obtained from the high-gain observer with GCD.

    Given the definition of the generalized conformable derivative, which encompasses some noninteger operators of conformable type, it may be very useful to obtain generalized results in control theory, such as modelling, control algorithms, observer design, stability proofs, among others. Thus, in this work an observation scheme, composed of a NQR-based estimator and a high-gain observer, was designed for a class of nonlinear systems which dynamics involve this derivative. With the aid of the proposed generalized conformable exponential function, it was proven that the observers were exponentially stable in a generalized way by means of some Lyapunov-like theorems. The estimators were applied to the generalized models of the simple pendulum and the Van der Pol oscillator, evaluating their performance by means of numerical simulations, and comparing these results with the ones obtained with their integer and fractional-order counterparts.

    According to the simulation results, the integer and fractional-order estimators had a similar performance regarding the behaviour of the estimation error, being slightly better the fractional-order case. However, both cases were outperformed by the one with the generalized conformable derivative, which exhibited a faster speed of convergence with similar or less overshoot, for both systems. However, these results may be improved with an adequate choice of gains and the selection of the kernel k(t). Furthermore, comparing the performance of each observer, both presented excellent outcomes, with slight differences that occurred due to the choice of their parameters.

    To conclude, the following may be considered as future directions for this work:

    ● One of the limitations of this work is that the performance of the generalized conformable derivative relies on the selection of the kernel k(t), which in this paper has been empirical. Thus, it may be possible to obtain optimal results with a structured selection method for the kernel.

    ● It may be interesting to apply this operator to other kind of estimators, like the reduced-order, algebraic, sliding-mode observers and the Kalman filter, and compare their advantages and theoretical considerations.

    ● Moreover, another theme of interest is the extension of the theory of the generalized conformable derivative to include time-varying and complex orders, and even real orders higher than one; besides, to consider also their discrete-time versions. It would be desired to study all these cases in order to compare the results obtained in applications, like the scheme proposed here, and in other existing systems and methodologies.

    ● Finally, some immediate potential applications of the generalized conformable derivative are related to the grey system. This model possesses several uses, such as the prediction of the production and consumption of combustibles and renewable energies, the calculation of wind turbine capacities, among others. Recently, the gray system has been extended to its conformable [57] and generalized conformable [58] versions; thus, there is a good amount of potential works related to these areas.

    Guillermo Fernández-Anaya would like to acknowledge the support and funding provided by División de Investigación y Posgrado (DINVP) de la Universidad Iberoamericana. Universidad Iberoamericana Ciudad de México, Departamento de Física y Matemáticas, Prolongación Paseo de la Reforma 880, Colonia Lomas de Santa Fe, Álvaro Obregón, Ciudad de México 01219, México.

    No potential conflict of interest was reported by the authors.



    [1] M. Maleki, I. Hashim, S. Abbasbandy, Analysis of IVPs and BVPs on semi-infinite domains via collocation methods, J. Appl. Math., 2012 (2012), 1–21. https://doi.org/10.1155/2012/696574 doi: 10.1155/2012/696574
    [2] S. A. Odejide, Y. A. S. Aregbesola, Applications of method of weighted residuals to problems with Semi-infinite domain, Rom. J. Phys., 56 (2011), 14–24.
    [3] F. Auteri, L. Quartapelle, Galerkin-Laguerre spectral solution of self-similar boundary layer problems, Commun. Comput. Phys., 12 (2012), 1329–1358. https://doi.org/10.4208/cicp.130411.230911a doi: 10.4208/cicp.130411.230911a
    [4] A. O. Adewumi, S. O. Akindeinde, A. A. Aderogba, Laplace-weighted residual method for problems with semi-infinite domain, J. Mod. Method Numer. Math., 7 (2016), 59–66.
    [5] H. F. Ismael, H. Bulut, H. M. Baskonus, W. Gao, Dynamical behaviors to the coupled schrdinger-boussinesq system with the beta derivative, AIMS Math., 6 (2021), 7909–7928. https://doi.org/10.3934/math.2021459 doi: 10.3934/math.2021459
    [6] H. Jafari, N. Kadkhoda, D. Baleanu, Lie group theory for nonlinear fractional K(m, n) type equation with variable coefficients, Meth. Math. Model. Comput. Complex Syst., 2021,207–227. https://doi.org/10.1007/978-3-030-77169-0_8
    [7] D. Baleanu, Z. B. Guvenc, J. Machado, New trends in nanotechnology and fractional calculus applications, Springer Netherlands, 2010. https://doi.org/10.1007/978-90-481-3293-5
    [8] A. Atangana, D. Baleanu, A. Alsaedi, New properties of conformable derivative, Open Math., 13 (2015), 1–10. https://doi.org/10.1515/math-2015-0081 doi: 10.1515/math-2015-0081
    [9] Gulnur Yel, C. Cattani, H. M. Baskonus, W. Gao, On the complex simulations with dark-bright to the hirota-maccari system, J. Comput. Nonlinear Dyn., 6 (2021), 16. https://doi.org/ 10.1115/1.4050677 doi: 10.1115/1.4050677
    [10] Y. M. Li, H. M. Baskonus, A. M. Khudhur, Investigations of the complex wave patterns to the generalized calogero-bogoyavlenskii-schiff equation, Soft Comput., 25 (2021), 6999–7008. https://doi.org/10.1007/s00500-021-05627-2 doi: 10.1007/s00500-021-05627-2
    [11] S. Rashid, S. Parveen, H. Ahmad, Y. M. Chu, New quantum integral inequalities for some new classes of generalized ϕ-convex functions and their scope in physical systems, Open Phys., 19 (2021), 35–50. https://doi.org/10.1515/phys-2021-0001 doi: 10.1515/phys-2021-0001
    [12] S. Rashid, D. Baleanu, Y. M. Chu, Some new extensions for fractional integral operator having exponential in the kernel and their applications in physical systems, Open Phys., 18 (2020), 478–491. https://doi.org/10.1515/phys-2020-0114 doi: 10.1515/phys-2020-0114
    [13] L. Xu, Y. M. Chu, S. Rashid, A. A. El-Deeb, K. S. Nisar, On new unified bounds for a family of functions via fractional-calculus theory, J. Funct. Space., 2020 (2020), 1–9. https://doi.org/10.1155/2020/4984612 doi: 10.1155/2020/4984612
    [14] S. Rashid, M. Can, D. Baleanu, M. C. Yu, Generation of new fractional inequalities via n polynomials s-type convexity with applications, Adv. Differential Equ., 2020 (2020), 1–20. https://doi.org/10.1186/S13662-020-02720-Y doi: 10.1186/S13662-020-02720-Y
    [15] S. Rashid, Z. Hammouch, D. Baleanu, M. C. Yu, New generalizations in the sense of the weighted non-singular fractional integral operator, Fractalsy, 28 (2020), 2040003. https://doi.org/10.1142/S0218348X20400034 doi: 10.1142/S0218348X20400034
    [16] S. Rashid, H. Kalsoom, Z. Hammouch, R. Ashraf, Y. M. Chu, New multi-parametrized estimates having pth-order differentiability in fractional calculus for predominating-convex functions in hilbert space, Symmetry, 12 (2020), 222. https://doi.org/10.1016/s0362-546x(01)00646-0 doi: 10.1016/s0362-546x(01)00646-0
    [17] Y. B. Yang, S. R. Kuo, H. H. Hung, Frequency-independent infinite elements for analysing semi-infinite problems, Int. J. Numer. Method Eng., 39 (1996), 3553–3569. https://doi.org/10.1002/(SICI)1097-0207(19961030)39:20<3553::AID-NME16>3.0.CO;2-6 doi: 10.1002/(SICI)1097-0207(19961030)39:20<3553::AID-NME16>3.0.CO;2-6
    [18] A. Akgül, A novel method for the solution of Blasius equation in semi-infinite domains, IJOCTA, 7 (2017), 225–233. https://doi.org/10.11121/ijocta.01.2017.00363 doi: 10.11121/ijocta.01.2017.00363
    [19] P. Berrut, M. S. Floater, G. Klein, Convergence rates of derivatives of a family of barycentric rational interpolants, Appl. Numer. Math., 61 (2011), 989–1000. https://doi.org/10.1016/j.apnum.2011.05.001 doi: 10.1016/j.apnum.2011.05.001
    [20] J. P. Berrut, S. A. Hosseini, G. Klein, The linear barycentric rational quadrature method for Volterra integral equations, SIAM J. Sci. Comput., 36 (2014), 105–123. https://doi.org/10.1137/120904020 doi: 10.1137/120904020
    [21] M. Floater, H. Kai, Barycentric rational interpolation with no poles and high rates of approximation, Numer. Math., 107 (2007), 315–331. https://doi.org/10.1007/s00211-007-0093-y doi: 10.1007/s00211-007-0093-y
    [22] G. Klein, J. Berrut, Linear rational finite differences from derivatives of barycentric rational interpolants, SIAM J. Numer. Anal., 50 (2012), 643–656. https://doi.org/10.1137/110827156 doi: 10.1137/110827156
    [23] G. Klein, J. Berrut, Linear barycentric rational quadrature, BIT Numer. Math., 52 (2012), 407–424. https://doi.org/10.1007/s10543-011-0357-x doi: 10.1007/s10543-011-0357-x
    [24] L. H. Wang, M. H. Hu, Z. Zhong, F. Yang, Stabilized lagrange interpolation collocation method: A meshfree method incorporating the advantages of finite element method, Comput. Method. Appl. M., 404 (2023), 115780. https://doi.org/10.1016/j.cma.2022.115780 doi: 10.1016/j.cma.2022.115780
    [25] Z. H. Qian, L. H. Wang, A meshfree stabilized collocation method (SCM) based on reproducing kernel approximation, Comput. Method. Appl. M., 371 (2020), 113303. https://doi.org/10.1016/j.cma.2020.113303 doi: 10.1016/j.cma.2020.113303
    [26] Z. H. Qian, L. H. Wang, Y. Gu, C. Z. Zhang, An efficient meshfree gradient smoothing collocation method (GSCM) using reproducing kernel approximation, Comput. Method. Appl. M., 374 (2021), 113573. https://doi.org/10.1016/j.cma.2020.113573 doi: 10.1016/j.cma.2020.113573
    [27] M. N. Rasoulizadeh, M. J. Ebadi, Z. Avazzadeh, O. Nikan, An efficient local meshless method for the equal width equation in fluid mechanics, Eng. Anal. Bound. Elem., 131 (2021), 258–268. https://doi.org/10.1016/j.enganabound.2021.07.001 doi: 10.1016/j.enganabound.2021.07.001
    [28] O. Nikan, Avazzadeh, An efficient localized meshless technique for approximating nonlinear sinh-Gordon equation arising in surface theory, Eng. Anal. Bound. Elem., 130 (2021), 268–285. https://doi.org/10.1016/j.enganabound.2021.05.019 doi: 10.1016/j.enganabound.2021.05.019
    [29] O. Nikan, Z. Avazzadeh, A locally stabilized radial basis function partition of unity technique for the sine-Gordon system in nonlinear optics, Math. Comput. Simul., 199 (2022), 394–413. https://doi.org/10.1016/j.matcom.2022.04.006 doi: 10.1016/j.matcom.2022.04.006
    [30] O. Nikan, Z. Avazzadeh, M. N. Rasoulizadeh, Soliton wave solutions of nonlinear mathematical models in elastic rods and bistable surfaces, Eng. Anal. Bound. Elem., 143 (2022), 14–27. https://doi.org/10.1016/j.enganabound.2022.05.026 doi: 10.1016/j.enganabound.2022.05.026
    [31] O. Nikan, Z. Avazzadeh, M. N. Rasoulizadeh, Soliton solutions of the nonlinear sine-Gordon model with Neumann boundary conditions arising in crystal dislocation theory, Nonlinear Dyn., 106 (2021), 783–813. https://doi.org/10.1007/s11071-021-06822-4 doi: 10.1007/s11071-021-06822-4
    [32] Z. Avazzadeh, O. Nikan, J. A. T. Machado, Solitary wave solutions of the generalized Rosenau-KdV-RLW equation, Mathematics, 8 (2020), 1601. https://doi.org/10.3390/math8091601 doi: 10.3390/math8091601
    [33] J. Li, Y. Cheng, Linear barycentric rational collocation method for solving heat conduction equation, Numer. Meth. Part. D. E., 37 (2021), 533–545. https://doi.org/10.1002/num.22539 doi: 10.1002/num.22539
    [34] J. Li, Y. Cheng, Linear barycentric rational collocation method for solving second-order Volterra integro-differential equation, Comput. Appl. Math., 39 (2020). https://doi.org/10.1007/s40314-020-1114-z
    [35] J. Li, Y. L. Cheng, Z. C. Li, Z. K. Tian, Linear barycentric rational collocation method for solving generalized Poisson equations, MBE, 20 (2023), 4782–4797. https://doi.org/10.3934/mbe.2023221 doi: 10.3934/mbe.2023221
    [36] J. Li, Y. Cheng, Barycentric rational method for solving biharmonic equation by depression of order, Numer. Meth. Part. D. E., 37 (2021), 1993–2007. https://doi.org/10.1002/num.22638 doi: 10.1002/num.22638
    [37] Z. Wang, S. Li, Barycentric interpolation collocation method for nonlinear problems, National Defense Industry Press, Beijing, 2015.
    [38] Z. Wang, Z. Xu, J. Li, Mixed barycentric interpolation collocation method of displacement-pressure for incompressible plane elastic problems, Chin. J. Appl. Mech., 35 (2018), 195–201.
    [39] Z. Wang, L. Zhang, Z. Xu, J. Li, Barycentric interpolation collocation method based on mixed displacement-stress formulation for solving plane elastic problems, Chin. J. Appl. Mech., 35 (2018), 304–309.
  • 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(1948) PDF downloads(77) Cited by(3)

Figures and Tables

Tables(15)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog