A New Fractional Discrete Memristive Map with Variable Order and Hidden Dynamics

: This paper introduces and explores the dynamics of a novel three-dimensional (3D) fractional map with hidden dynamics. The map is constructed through the integration of a discrete sinusoidal memristive into a discrete Duffing map. Moreover, a mathematical operator, namely, a fractional variable-order Caputo-like difference operator, is employed to establish the fractional form of the map with short memory. The numerical simulation results highlight its excellent dynamical behavior, revealing that the addition of the piecewise fractional order makes the memristive-based Duffing map even more chaotic. It is characterized by distinct features, including the absence of an equilibrium point and the presence of multiple hidden chaotic attractors.


Introduction
Fractional calculus deals with fractional derivatives (for continuous-time systems) as well as fractional difference operators (for discrete-time systems) [1].Starting from 1974, attention has been focused on the potential of discrete fractional calculus, so that several non-integer-order difference operators have been proposed [2,3].
Unlike integer-order derivatives, fractional-order derivatives offer increased accuracy, serving as effective tools for characterizing memory effects across diverse materials and processes.Recently, there has been a notable surge of efforts directed towards exploring discrete fractional calculus [3].Greater attention has been dedicated to studying chaotic behaviors in fractional discrete systems nonlinear maps described by non-integer-order difference equations, showcasing heightened sensitivity to initial conditions [4][5][6][7].For instance, in reference [4], the emergence of chaos within a fractional logistic map is revealed, while reference [5] examines the delayed version of the previously studied map.Furthermore, in [6], the chaotic behaviors of non-integer-order higher-dimensional multicavity chaotic maps have been analyzed.In recent research, Wu et al. [8] developed a method to create fractional difference equations with variable orders designed to exhibit short memory behavior, which are useful in characterizing the memory properties of discrete systems.However, limited exploration has been conducted into variable-order fractional maps where the fractional order changes over discrete time.For instance, reference [9] delved into exploring the chaotic behavior of a variable-order fractional Tinkerbell map through bifurcations and phase portraits.
Very recently, researchers have explored the manifestation of chaos in discrete systems involving fractional discrete memristors [10][11][12][13][14][15].These non-integer-order maps originate from the concept of memristors, which represent nonlinear circuit elements relating electric charge and magnetic flux linkage [13].For example, reference [11] introduces a novel chaotic fractional memristive-based map with concealed attractors, revealing diverse dynamical behaviors, including coexisting hidden dynamics and initial offset boosting.In [12], mathematical models of fractional discrete memristors are illustrated, although further details on the specific aspects depicted could enhance clarity.Additionally, diverse chaotic behaviors were observed, indicating that the complexity of the systems is influenced by the memristive parameters [12].
Drawing from previous considerations regarding systems with variable fractional orders, this paper aims to contribute to the field of discrete memristive systems.The primary contribution lies in the introduction of a novel 3D fractional discrete memristive map characterized by a variable order and a short memory effect.This achievement is realized by integrating a discrete memristive with the Duffing map.Building upon the framework outlined in reference [8], we propose a fractional-order discrete memristive-based Duffing map with short memory, thereby introducing greater flexibility to the utilized models.The resulting system exhibits distinctive features, including the coexistence of multiple hidden chaotic attractors.The structure of the manuscript is as follows: Section 2 presents the fractional discrete memristive-based Duffing map through the integration of a discrete sinusoidal memristor into a discrete Duffing map.Moreover, a mathematical operator, namely, the fractional variable-order Caputo-like difference operator, is employed to establish the fractional form of the map with short memory.In Section 3, the chaotic dynamics of the novel fractional map are analyzed under piecewise fractional-order functions over three intervals.Specifically, the proposed fractional map shows the coexistence of different types of hidden attractors.The multistability phenomena of hidden attractors are found and investigated in detail.Finally, 0-1 is computed as a measurement of the degree of complexity of the dynamic behavior of the conceived fractional discrete memristive-based Duffing map.This computation aims to demonstrate the advantage of using short memory in modeling fractional memristive-based maps.

Preliminaries
Among the various difference operators established in discrete fractional calculus, this paper specifically focuses on the ν-Caputo-like difference operator.The ν-Caputo fractional difference operator for a function g(s) is defined as [16] : where s ∈ N a+m−ν , m = ⌈ν⌉ + 1 and ν ̸ ∈ N. (s − τ − 1) (m−ν−1) with the fractional sum of order ν for a function g on a time scale N a = {a, a + 1, a + 2, . . .} is defined as with s ∈ N a+ν .

Model Description
Recently, the addition of a discrete memristor to a two-dimensional Duffing map has resulted in the development of a 3D discrete memristor-based Duffing map, rendering this study's coupled plane intriguing [17].According to the research conducted by Bao et al. [17], the proposed model exhibited multistable behavior.In this work, we present an improved three-dimensional fractional memristive Duffing map with hidden dynamics and provide a detailed analysis of its dynamics, including the examination of equilibrium points and their stability.Through introducing the Caputo's difference operator, a new type of discrete memristive-based Duffing map is established, which can be represented as follows: ( where the ν-Caputo fractional difference operator is denoted by C ∆ ν a .Although we have devised a novel approach utilizing a sinusoidal discrete memristive to preserve lengthy memories beginning from their initial states, this technique is ineffective due to its large storage capacity.From a practical perspective, certain states may not require storage, indicating that only a short memory model is necessary.Recently, Wu et al. [8] proposed a technique for constructing fractional difference equations with variable order.Their study demonstrated how fractional variable order can describe the short-memory property of discrete systems.Inspired by this study, we now propose the following variable-order analogue of the fractional map (3): where h ∈ R + denotes the discretization step size and a is the initial state.Let L be a positive integer and Beside each subinterval, the fractional order ν(s) is an important parameter to describe memory effects.Let ν(s) be a piecewise constant function defined by: Notice that the initial conditions are initialized at each interval.As pointed out by Wu in [8], the fractional difference equation C ∆ ν(s) x(t) reduces to C ∆ ν 0 x(t) with initial condition s 0 = a; while when s ∈ (T 1 , T 2 ], the initial condition is considered to s 1 and the fractional difference equation C ∆ ν(s) x(t) reduces to C ∆ ν 1 x(t).Therefore, by using the ν-fractional sum defined above, the numerical formula can be obtained for n = 0, . . ., L − 1: for n = L, . . ., 2L − 1 . . .This is a new class of fractional variable-order memristive-based Duffing map which holds "short-memory effects" to the past L states.As one can see from Equation (4), when n = 0, . . ., L − 1, the new fractional discrete memristive-based Duffing map (5) state depends on the past L variable, i.e, x(0), x(1), . . ., x(L − 1); while when n = L, . . ., 2L − 1 x(i) depends on x 1 (L), . . ., x 1 (2L − 1).Hence, generally when n = kL, . . ., kL + L − 1, the state depends on x(kL), . . ., x(L(k + 1) − 1).Therefore, we can investigate the chaotic behavior on each sub-domain.This idea of short memory provides more space and chaos in the fractional discrete systems and makes the numerical computation more practical.

Analysis of Equilibrium Point
Examining equilibrium points is a fundamental method for evaluating system stability.We may find the equilibrium points of the suggested fractional memristive-based Duffing map (FMD) (4) with piecewise function variable order and short memory by solving these equations: Upon examining Equation ( 8), it becomes evident that the system parameter γ acts as a pivotal controller, influencing the nature of equilibria and functioning as a bifurcation parameter.Notably, when γ is not equal to zero, Equation ( 8) has no solution, meaning that the novel FMD map with short memory has no equilibrium point.According to recent classification, studies of chaotic attractors can be divided into two classes: self-excited attractors and hidden attractors.In particular, a chaotic attractor is called hidden if its basin of attraction does not intersect with a small neighborhood of the equilibria, while a self-attractor has a basin of attraction that exists from an unstable equilibrium point [18].Therefore, one can confirm that the proposed map (4) with no equilibrium point belongs to the family of hidden attractors.

Dynamical Analysis
To demonstrate the effectiveness of the fractional variable order and short memory in the discrete memristive-based Duffing map (4), we set L = 4000 and m = 3.Our analysis focuses on three intervals: [T 0 , T 1 ], (T 1 , T 2 ], and (T 2 , T 3 ], each corresponding to fractional variable orders ν 0 , ν 1 , and ν 2 , respectively.Numerical simulations include bifurcation diagrams of the state variable z, maximum Lyapunov exponents, phase diagrams, and time-state representations.The specific parameter configurations and simulation results are detailed below.

Coexisting Multiple Hidden Attractors
Multistability in a nonlinear system exhibits a rich diversity of stable states, augmenting the system's adaptability.Specifically, as the number of coexisting attractors in a dynamical system grows to infinity, the coexistence of infinitely many attractors linked to the initial condition of a particular state variable is believed to epitomize extreme multistability.For the new fractional memristive Duffing map, varied initial states can lead to the coexistence of multiple attractors within specific regions of the control parameter α 1 .The bifurcation diagram is generated using constant parameters: α 2 = 0.244, K = −2, and γ = −2.Figure 1 illustrates the bifurcation diagrams for the variable z, specifically for the initial condition (0.3, −0.2, 3), represented by the red curve.The blue curve is generated with the initial condition (0.3, −0.2, 1), and the black curve represents the initial condition (0.3, −0.2, 0.1).For the fractional-order values, ν 0 = 0.99 in [T 0 , T 1 ], ν 1 = 0.98 in (T 1 , T 2 ], and ν 2 = 0.97 in (T 2 , T 3 ].The plot of the maximum Lyapunov exponent, linked to the bifurcation diagram of initial condition (0.3, −0.2, 0.1), is generated by varying the control parameter α 1 across the range of 1.5 to 1.95.From the bifurcation diagram presented in Figure 1 and its corresponding graph of the maximum LEs, it is evident that the analyzed model is capable of exhibiting the phenomena of multistability.Upon closer examination of Figure 1, it becomes apparent that all bifurcation diagrams share a consistent structure as the IC z 0 variation.This observation implies that the suggested discrete fractional FMD map displays extremely complicated dynamics behavior, including chaos, periodicity, period doubling bifurcation, and homogeneous multiple attractors.It is noteworthy that the system exhibits regular dynamics when the value of α 1 is less than or equal to 1.75.However, beyond this threshold (α 1 ≥ 1.75), the system undergoes a transition into chaotic behavior through period-doubling bifurcation.Furthermore, Figure 2 showcases distinct hidden attractors, emphasizing the diverse dynamic characteristics of the novel FMD map with variable order and short memory as α 1 varies.To reveal the phenomena of homogeneous multistability in the FMD map with fractional and variable order, the previous parameters are maintained unchanged, with α 1 set to 1.8, and the initial condition chosen as (1, −0.2, z 0 ).When the memristive initial value z 0 varies in the range [−4, 6], the bifurcation diagram of the state variable z max in the short interval (T 2 , T 3 ] for the memory L = 3000 and the corresponding Lyapunov exponents are shown in Figure 3.We observe that the bifurcation diagram contains numerous elongated strip-shaped regions, each with the same shape but positioned differently.These elongated strips have the potential to transform into a single structure attractor.Moreover, Figure 3b reveals that the FMD with short memory and variable order possesses an unchanging Lyapunov exponent.This suggests its ability to consistently generate a persistent shape and structure for the hidden attractor.In Figure 4, the coexistence of hidden attractors is demonstrated for three distinct initial conditions: z 0 = 0.1, z 0 = 2, and z 0 = 3. Figure 4a displays three periodic hidden attractors, all sharing a consistent structure but with different positions.
Similarly, set α 1 to 1.8 and leave the bifurcation parameters unaltered.The unique FMD with variable order and short memory exhibits multistability, as demonstrated by the occurrence of three coexisting chaotic attractors.This coexisting of the homogeneous attractors property enhances the suitability of the fractional Duffing map for various chaos-based industrial applications.Moreover, compared with the constant order results shown in Figure 2, it can be seen that introduced variable order has a huge impact on the multistability phenomena of the fractional memristive-based Duffing map.

The Invariance Property
We consider the parameters γ and α 1 , with α 1 set to 1.8 and γ set to 0.09.The system parameter α 2 of the fractional memristive-based Duffing map (4) is assigned as a bifurcation parameter.For the piecewise fractional-order function ν(s): To illustrate the invariance property clearly, the bifurcation diagrams are plotted for two distinct initial conditions and values of K.The red diagram is plotted for (x 0 , y 0 , z 0 , K) = (−0.3,−0.2, 0.1, 2), whereas the blue diagram is essentially plotted for (x 0 , y 0 , z 0 , K) = (0.3, 0.2, 0.1, −2).As α 2 increases from 0.1 to 0.6 in Figure 5, the system displays a range of dynamics, such as coexisting behaviors, periodic behavior, and chaos.The coexisting behaviors are more observable during the periodic window.3.3.Effect of the Fractional Variable Order (ν 1 , ν 2 , ν 3 ) The fractional order ν affects the hidden dynamics of the FMD with short memory and variable order but has no effect on the existence of the equilibrium point.The fixed initial condition (1, −0.2, 0.1), α 1 = 1.77, α 2 = 0.24, K = −2, and γ = 0.09.We examine the dynamics of the piecewise fractional-order function over three intervals: [T 0 , T 1 ], (T 1 , T 2 ], and (T 2 , T 3 ]: [8001, 12,000]. The 0-1 test for chaos is a recent method used to distinguish between regularity and chaos.The concept was initially proposed by Gottwald and Melbourne in [19] and subsequently extended to both continuous and discrete models for analyzing actual data.This test aids in identifying periodic or chaotic attractors.In practical terms, the system is characterized as chaotic when the translation components p(r) = ∑ r i=1 g(i) cos(ic) and q(r) = ∑ r i=1 g(i) sin(ic), where r ∈ {1, 2, . . ., N} and c is chosen randomly in (0, π), exhibit Brownian-like trajectories in the p − q plane and when the asymptotic growth rate K approaches 1.However, the system becomes regular as K approaches 0 and p and q show bounded-like trajectories.Here, we use the 0-1 test based on the generated sequence (x i ) i=1:N to analyze the presence of chaos of the FMD with short memory.The bifurca- tion diagram and the asymptotic growth rate K of the variable-order map are shown in Figure 6. Figure 6 illustrates how, for the short memory L = [8000, 12,000], the state of the map (6) shifts from period states to chaos once more as ν 2 decreases for the short memory L = [8000, 12,000].For ν 2 ∈ (0.82, 1), the novel map exhibits chaos.The suggested FVDM (4) with piecewise function variable order is shown to have both periodic and chaotic hidden attractor behavior.It should be noted that the states evolutions for the FVDM (4) variable order and short memory differ from those of the fractional discrete memory-based Duffing map with constant order.Specifically, as ν decreases below 0.84, the fractional memristive-based Duffing map with constant order ν = ν 0 = ν 1 = . . .ν m enters an unbounded attractor, while the system (4) suggested herein reveals hidden chaotic behaviors, suggesting that the adoption of the short memory principle has increased the chaotic behavior of the memristive-based Duffing map.
We have chosen to graph the chaotic attractors for the incommensurate orders ν 0 , ν 1 , and ν 2 in the three distinct intervals [0, 4000], [4000, 8000], [8000,12,000], respectively, to illustrate the varying dynamics produced by the variable order.The results are shown in Figure 7.The type of attractor in the first integral is distinctly different, and the map exhibits periodic behavior.Furthermore, when the piecewise order shifts to the second interval, chaos emerges, demonstrating that the piecewise function amplifies the chaotic dynamics of the fractional memristive-based Duffing map.Remark 1.As a result of the previous numerical findings, we can conclude that there is a significant improvement in the chaotic area of the two-dimensional Duffing map when applying the sinusoidal discrete memristive and piecewise fractional-order function.Moreover, compared with the constant fractional-order Duffing map shown in [20], it can be seen that introducing the variable order has a large impact on the multistability phenomena of the fractional Duffing map.It is worth noting that all the applications described above regarding the use of variable order are new and have not been applied before in discrete memristor maps.The only paper in the literature that deals with chaos in variable order is [9].First of all, it should be noted that in reference [9], the variable order maps are achieved using very complex variable-order functions and the variable order does not provide a short memory effect, whereas the piecewise variable-order technique here involves a short memory effect, implying the effective application of the novel system.

Conclusions
In this study, fractional variable order and a discrete sinusoidal memristive are added to the Duffing map.In essence, we suggest a Duffing map with short memory that is based on a variable-order fractional discrete memristive.We have theoretically evaluated the stability of this map, which suggests that it can produce hidden attractors.The new fractional map demonstrated the presence of many attractors and chaos, according to the numerical results.Furthermore, our results have demonstrated that the fractional Duffing map becomes more complex as a result of the addition of the memristive and variable order.Because of its rich and complex behavior, this research can support theoretical understanding and advancements in encryption and secure communication.Utilizing variable order as confidential information can enhance the range of security choices, leading to the development of a novel picture encryption technique based on chaos.