On Hyers–Ulam Mittag-Leffler stability of discrete fractional Duffing equation with application on inverted pendulum

A human being standing upright with his feet as the pivot is the most popular example of the stabilized inverted pendulum. Achieving stability of the inverted pendulum has become common challenge for engineers. In this paper, we consider an initial value discrete fractional Duffing equation with forcing term. We establish the existence, Hyers–Ulam stability, and Hyers–Ulam Mittag-Leffler stability of solutions for the equation. We consider the inverted pendulum modeled by Duffing equation as an example. The values are tabulated and simulated to show the consistency with theoretical findings.


Introduction
The understanding of the real-world problems by replicating into mathematical models proves to be an effective tool. Analyzing the developed model thus provides a wide insight into the considered phenomena. In [1], Rayleigh introduced a nonlinear damping function in a second-order oscillator equation. Extensive study of this equation using vacuum-tube circuits for analogue simulation was carried out by Van der Pol [2]. Besides, a model of heart beat was constructed using an electrical circuit with coupled relaxation oscillators and simulations of normal heart beat and of certain disorders were convincingly obtained by Van der Pol and Van der Mark [3]. Certain damped and driven oscillators are modeled by the Duffing equation, a second-order differential equation with cubic nonlinearity named after Georg Duffing [4]. The motion of a damped oscillator described by the equation has more complex potential than simple harmonic motion. This equation is used to illustrate the motion of a mass attached to a nonlinear spring and a linear damper. The Duffing equation is given bÿ x + θẋ + δx + ηx 3 = γ sin(ωt), (1.1) in modeling hard spring oscillators (iron core inductor), soft spring oscillators (nonlinear LRC circuit), negative stiffener, and nonharmonic oscillator [5]. The inverted pendulum equation, which is framed from the Duffing equation, is used in rocket propeller, segway, and robotics.
Fractional calculus, which is regarded as 21st century calculus, has its origin during the same period as that of the ordinary calculus. Continuous fractional-order equations prove to be significant in modeling nuclear reactor dynamics, chaotic dynamical systems, chemical kinetics, population dynamics, and circuit theory [6]. Qualitative analysis of the solutions of fractional-order equations representing real-life phenomena plays a predominant role in understanding the nature and behavior of the models [7,8]. Intensive interest shown by researchers during this decade toward discrete fractional calculus demands the need for the development of the methods equivalent to the fractional differential equations. This opinion is very much strengthened by the increase in number of researchers involved in the development of the methodology for discrete fractional calculus [9][10][11][12].
The works by Anastassiou [13], Atici and Eloe [14][15][16][17], Goodrich [18], and Miller and Ross [19] have laid the foundations for the field of discrete fractional calculus. Time-scale calculus unifies the theory of difference equations with that of differential equations, and qualitative properties such as oscillation and nonoscillation of the dynamic equations on discrete time scales were studied in [20][21][22]. Chen [23,24] was the first author who studied the stability results of the nonlinear fractional difference equations. The response given by Hyers to the question put forth by Ulam during his talk [25,26] was marked as the origin for the study on stability of functional equations. The Ulam stability of integerand arbitrary-order differential equations were established in [27][28][29][30]. Recently, boundary value impulsive integrodifferential equations and coupled systems of Hilfer-Hadamardtype fractional equations are considered for discussion of stability in the Hyers-Ulam sense [31,32]. Ulam stability analysis of nabla fractional equations was carried out in [33][34][35], and in [36], the Ulam-Hyers stability of discrete fractional boundary value problems was studied. Here we consider the discrete-time forced fractional-order Duffing equation without damping.
Denote Q := [j + 2, j + T] ∩ N j+2 , where T ∈ N and N j = {j, j + 1, . . .}, j ∈ R. Our main equation has the form where υ * is the Caputo like difference operator, δ and η control the linear stiffness and nonlinearity in restoring force, p : Q → R is the driving force with A, B ∈ R + . The restoring force represented by δψ + η(ψ) 3 is vital in determining the nature of the spring to be used in the model. The positive real values of δ and η describe the hardening spring, and η < 0 (δ > 0) denotes soft spring.
The choice of the operator plays a crucial role in developing models arising in physics. Here the construction of the physical model using the Caputo difference operator is preferred over the Riemann-Liouville operator to overcome some limitations of the Riemann-Liouville operator in modeling real-life problems. One of the limitations concern the initial conditions defined for the physical problems. Initial conditions of a Caputotype fractional difference operator are traditional integer-order ( k , k ∈ N 0 ) conditions, whereas for a Riemann-Liouville type operator, they are defined in terms of a fractional sum or difference terms ( α , α ∈ R), which fail to provide physical interpretation for the model.
The plan of the paper is as follows. Section 2 imparts some necessary definitions, lemmas, and an existence result. Section 3 presents the Hyers-Ulam stability followed with Hyers-Ulam Mittag-Leffler stability in Sect. 4. Appropriate examples accompanied with simulation are provided in Sect. 5.

Mathematical background
In this section, we provide some fundamental definitions and lemmas and state an existence result used throughout this work.

Definition 2.1 ([16])
The Fractional sum of order ν > 0 for a function p is given by where p is defined for r = j mod (1), and -ν ψ is defined for n = (j + ν) mod (1) and In particular, if 1 < γ < 2 and j = 0, then this relation becomes where p is defined on N 2 .
where n ∈ Q.
Proof Let ψ be a solution of (1.2). Then from (2.2) we have which satisfies (2.3). On the other hand, if ψ is a solution of (2.3), then by comparison of (2.2) and (2.3) we get which takes the form Letting n = 1, 2, 3, . . . yields It is evident that ψ satisfies (1.2). The proof is complete.
To ensure the existence of solution, we consider To prove the existence, we apply the Krasnoselskii fixed point theorem. Let W be a nonempty, closed, bounded, and convex subset of a Banach space (S; · ). Suppose that 3 .

Theorem 2.6 Problem (1.2) has a solution in the set W , provided that
where p = sup n∈Q |p(n)|. Moreover, the solution is unique if Proof First, we can easily check that |F 1 ψ(n)| ≤ A + BT. On the other hand, we compute We can choose the constants δ, η, T such that A + BT + W * ≤ K . For such a choice, |F 1 ψ + The continuity of F 1 is easy to check as it is just a function of n. Besides, the set F 1 W is bounded. Thus F 1 W is contained in a compact subset of W .
On the other hand, we can easily see that for y 1 , y 2 ∈ W and F = F 1 + F 2 , we get We can see that for T ≥ 2, K * is nonnegative. Hence under the condition K * = 1 Γ (υ) (T + υ -1) υ-1 (T -2)(δ + η3L 2 ) ≤ 1, the mapping F is a contraction. Applying the Banach fixed point theorem, the existence and uniqueness of solution is ensured. We can see that we do not need any additional assumptions for uniqueness.

Hyers-Ulam stability
This section provides results on the Hyers-Ulam stability of (1.2).
We make the following assumptions before proving the stability of (1.2).

Theorem 3.4
Assume that (H2) holds. Let φ ∈ R solve (3.1) for some > 0, and let ψ ∈ R be the solution of
Proof If φ solves (4.1), then using Remark (4.2) and (2.2), we have that the solution to (B2) satisfies This completes the proof.

Applications
Springs are elastic in nature, and thus the original shape is regained after it is subject to some stress. They follow Newton's third law of motion: the harder you pull, the harder it hits back. The rapid increase in restoring force of a spring than suggested by Hooke's law remains the main criterion for a spring to be nonlinear and hard. The motion of such hard springs plays a significant role in study and understanding of nonlinear physics. For a hard spring oscillator, the increase in amplitude results in decrease in period.
In this section, we consider a pendulum with its center of mass above its pivot point. One of the common challenges for engineers and researchers is achieving the stability of an inverted pendulum. The applications of the inverted pendulum varies from personal transporters to electronic unicycles. The motion of an inverted pendulum with forcing term (τ sin(ωt)) described by Duffing equation is [5] m 2ẍ + bẋ + (kmg )x + 1 6 mg x 3 + τ sin(ωt) = 0, (5.1) where bẋ is the damping term, k is the stiffness constant, m is the mass provided at the pendulum top by two or more strong magnets, is the effective length of the pendulum, g is the acceleration due to gravity, the maximum torque is denoted by τ , and ω is the driving frequency.   where n ∈ [0, 13] ∩ N 0.5 . We will prove the Hyers-Ulam stability of (5.3). Straightforward calculations show that E n + 1.5, ψ(n + 1.5) = 0.01ψ(n + 1.5) + 0.2 cos(n + 1.5) is Lipschitz continuous with L = 0.01. The value of ξ in Theorem 3.4 is 0.3224 < 1. We now ensure that inequality (3.1) holds. Let = 0.72 and ψ(n) = n (2) 20 for n ∈ [0, 13] ∩ N 0 . Then 1.5 * ψ(n) + 0.01ψ(n + 1.5) + 0.2 cos(n + 1.5) = -0.5 2 (n) (2) 20 + 0.0979 + 0.2 cos(n + 1.5) ≤ 0.7008 < .

Conclusion
Following the trend in investigating equations of fractional order, we consider a discrete fractional form of Duffing equation with forcing term. We accommodate the newly established discrete fractional calculus to determine sufficient conditions for the existence, Hyers-Ulam, stability, and Hyers-Ulam Mittag-Leffler stability for the addressed equation. We analyze practical examples describing the undamped inverted pendulum and forced simple harmonic case as applications of the theoretical results. Stability conditions are obtained numerically for different lengths of the pendulum and the values are thus tabulated and represented graphically. We believe that results of this paper are of utmost importance for audience engaged in studying stability of mathematical models describing real physio-electrical phenomena.