Fract. Calc. Appl. Anal. manuscript No.

(will be inserted by the editor)

## Asymptotic cycles in fractional generalizations of multidimensional maps

Mark Edelman 1

Received: ...... / Revised: ...... / Accepted: ......

Abstract In regular dynamics, discrete maps are model presentations of discrete dynamical systems, and they may approximate continuous dynamical systems. Maps are used to investigate general properties of dynamical systems and to model various natural and socioeconomic systems. They are also used in engineering. Many natural and almost all socioeconomic systems possess memory which, in many cases, is power-law-like memory. Generalized fractional maps, in which memory is not exactly the power-law memory but the asymptotically power-law-like memory, are used to model and investigate general properties of these systems.

In this paper we extend the definition of the notion of generalized fractional maps of arbitrary positive orders that previously was defined only for maps which, in the case of integer orders, converge to area/volume-preserving maps. Fractional generalizations of H´non and Lozi maps belong to the newly defined e class of generalized fractional maps. We derive the equations which define periodic points in generalized fractional maps. We consider applications of our results to the fractional and fractional difference H´non and Lozi maps. e

Keywords fractional maps (primary) · periodic points · power-law memory ·

bifurcations

Mathematics Subject Classification (2010) 26A33 (primary) · 47H99 · 34A99 · 37G15 · 70K50 · 39A70

Mark Edelman 1 , ∗

Department of Physics, Stern College at Yeshiva University 245 Lexington Ave, New York, NY 10016, USA Courant Institute of Mathematical Sciences, New York University 251 Mercer St., New York, NY 10012, USA E-mail: edelman@cims.nyu.edu ∗ corresponding author

## 1 Introduction

Many natural and socioeconomic systems possess memory. Systems with memory are also widely used in engineering. In many cases, this memory is a powerlaw-like memory. Maps with memory may be used to describe discrete systems with memory, to model continuous systems with memory, and all numerical schemes to solve fractional differential equations may be considered as maps with memory. A review of the economic applications of maps with memory may be found in [1] (see also [2,3]). Fractional maps allow modeling of the distributions of lifespans of living species consistent with the Gompertz distribution and a limited lifespan [4]. There is an abundance of papers in which fractional maps are used to model memristors (see, e.g, [5]), to encrypt images and signals (see, e.g. [6,7]), to model spread of infections (see, e.g., [8]), to control systems (see, e.g., [9,10]), etc.

The first maps with memory introduced and investigated in scientific publications were not fractional - the memory was directly introduced in the forms which were relevant to the physical formulation of the corresponding problems (see [11,12,13,14,15,16]). In the one-dimensional case, maps with memory may be written as

$$x _ { n + 1 } = \sum _ { k = 0 } ^ { n } V _ { \alpha } ( n, k ) G _ { K } ( x _ { k } ), & & ( 1. 1 )$$

where V α ( n, k ) and α characterize memory effects and K is a parameter. In many cases, maps are convolutions with V α ( n, k ) = V α ( n -k ). The initial applications were related to the case of exponentially decaying memory. It has been shown in [17] that directly introduced maps with power-law memory ( V α ( n, k ) = ( n -k ) α -1 coincide with the fractional maps introduced in [18] as solutions of fractional equations describing kicked systems. In the case of the fractional maps with time step h , corresponding fractional differential equations are recovered in the limit h → ∞ (see [17]). In 2006, Stanislavsky investigated the logistic map with memory, whose kernel was taken from an algorithm of the numerical fractional integration [19]. Fractional maps, introduced in [18], were initially investigated in [20,21,22,23,24,25].

Discrete fractional calculus was proposed in [26,27] (see also more recent reviews [9,10,28,29,30,31]). Solutions of fractional difference equations were initially investigated in [32,33,34,35,36]. In [37,38,39,40], the authors showed that solutions of fractional difference equations may be written as fractional (we will call them fractional difference) maps with the falling factorial memory.

It soon became obvious that fractional and fractional difference maps may be considered as particular forms of a wider class of maps which possess similar properties [4,30,31,41,42,43,44]. The notion of the generalized fractional maps (GFM) was introduced in [45], in which the equations defining asymptotically periodic points were derived. The definition of GFM was further considered in [46,47,48], in which the authors investigated maps' asymptotically periodic points, stability, and bifurcations.

GFM are particular forms of the Volterra difference equations of convolution type (see, e.g. Chapter 6 in [49]) with kernels, whose differences are absolutely summable, but the series of kernels are diverging. Proved in [50] Theorem 1.1 is the most important property of the GFM's kernels, which allows to investigate the asymptotic behavior of GFM.

GFM introduced in the above cited papers are maps which, in the case of integer orders, produce area/volume preserving maps. In this paper we extend the notion of GFM to include multi-dimensional maps, like H´ enon and Lozi maps, which are not area/volume preserving.

## 2 Preliminaries

As we mentioned in Section 1, Theorem 1.1 of [50] is the most important property of the GFM's kernels. It is formulated as the following:

## Theorem 1 If

$$\lim _ { k \to \infty } x _ { k } = X, \ \lim _ { N \to \infty } \sum _ { k = 0 } ^ { N } a _ { k } = S,$$

and the series is converging absolutely ( { a k } ∈ l 1 ), then

$$\Sigma = \lim _ { N \to \infty } \sum _ { k = 0 } ^ { N } x _ { N - k } a _ { k } = X S,$$

where k, N ∈ N 0 and x , a k k , X, S, Σ ∈ R . If S = ∞ , and Σ is finite, then X = 0 .

## 2.1 Regular H´ enon, Lozy, and p -dimensional maps

The H´ enon map [51], written as a difference equation, is:

$$x _ { n + 1 } - x _ { n } & = 1 - x _ { n } - a x _ { n } ^ { 2 } + y _ { n }, & ( 2. 3 ) \\ y _ { n + 1 } - y _ { n } & = b x _ { n } - y _ { n }. & ( 2. 4 )$$

In the classical H´non map, the parameter values are: e a = 1 4 and . b = 0 3. . The Lozi map [52], written as a difference equation, is:

$$x _ { n + 1 } - x _ { n } = 1 - x _ { n } - a | x _ { n } | + y _ { n },$$

$$y _ { n + 1 } - y _ { n } = b x _ { n } - y _ { n }.$$

In the general case of p -dimensional maps, the equations can be written as

$$\Delta x _ { i } ( t ) = - G _ { i, K } ( x _ { 1 } ( t ), x _ { 2 } ( t ), \dots, x _ { p } ( t ) ), \ 1 \leq i \leq p, \ t \in \mathbb { N }, \quad ( 2. 7 )$$

where the forward difference operator ∆ is defined as

$$\Delta x ( n ) = x ( n + 1 ) - x ( n ),$$

with the initial conditions x i (0) = x i, 0 . These equations, with x i,n = x i ( n ) for t = n , can be written as

$$\Delta x _ { i, n } = - G _ { i, K } ( x _ { 1, n }, x _ { 2, n }, \dots, x _ { p, n } ), \ 1 \leq i \leq p, \ n \geq 0, \quad ( 2. 9 )$$

which are equivalent to the map:

$$x _ { i, n + 1 } = F _ { i, K } ( x _ { 1, n }, x _ { 2, n }, \dots, x _ { p, n } ), \ 1 \leq i \leq p, \ n \geq 0 \quad \ ( 2. 1 0 )$$

with

$$G _ { i, K } ( x _ { 1, n }, x _ { 2, n }, \dots, x _ { p, n } ) = x _ { i, n } - F _ { i, K } ( x _ { 1, n }, x _ { 2, n }, \dots, x _ { p, n } ), \quad ( 2. 1 1 )$$

where K stands for a set of parameters.

## 2.2 Definition of fractional and fractional difference maps

The system Eq. (2.7) turns into a system of fractional Caputo delta h -difference equations when the first differences ∆x i ( ) are substituted by the orders 0 t &lt; α &lt; i 1 Caputo forward h -differences (see review [30] and references therein)

$$& ( _ { 0 } \Delta _ { h, * } ^ { \alpha _ { i } } x _ { i } ) ( t ) = \\ & - G _ { i, K } ( x _ { 1 } ( t + ( \alpha _ { i } - 1 ) h ), x _ { 2 } ( t + ( \alpha _ { i } - 1 ) h, \dots, x _ { p } ( t + ( \alpha _ { i } - 1 ) h ), ( 2. 1 2 )$$

where t ∈ { h N 1 } , with the initial conditions

$$x _ { i } ( 0 ) = x _ { i, 0 }, \ 1 \leq i \leq p,$$

which is equivalent to the system of Volterra difference equations of convolution type (see Theorem 9 in [30]):

$$& x _ { i, n + 1 } = x _ { i, 0 } \\ & - \frac { h ^ { \alpha _ { i } } } { \Gamma ( \alpha _ { i } ) } \sum _ { s = 0 } ^ { n } ( n - s - 1 + \alpha _ { i } ) ^ { ( \alpha _ { i } - 1 ) } G _ { i, K } ( x _ { 1, s }, x _ { 2, s }, \dots, x _ { p, s } ), \quad ( 2. 1 4 )$$

where x i,k = x i ( kh ). The definition of the falling factorial t ( α ) is

/negationslash

$$t ^ { ( \alpha ) } = \frac { \Gamma ( t + 1 ) } { \Gamma ( t + 1 - \alpha ) }, \ t \neq - 1, - 2, - 3 \dots. \quad \quad ( 2. 1 5 )$$

The falling factorial is asymptotically a power function:

$$\lim _ { t \to \infty } \frac { \Gamma ( t + 1 ) } { \Gamma ( t + 1 - \alpha ) t ^ { \alpha } } = 1, \quad \alpha \in \mathbb { R }. \quad \quad ( 2. 1 6 )$$

The h -falling factorial t ( α ) h is defined as

$$& t _ { h } ^ { ( \alpha ) } = h ^ { \alpha } \frac { \Gamma ( \frac { t } { h } + 1 ) } { \Gamma ( \frac { t } { h } + 1 - \alpha ) } = h ^ { \alpha } \left ( \frac { t } { h } \right ) ^ { ( \alpha ) }, \\ & \frac { t } { \iota } \neq - 1, - 2, - 3, \dots.$$

$$\frac { \iota } { h } \neq - 1, - 2, - 3, \dots$$

/negationslash

After substitution

$$G _ { i } ^ { 0 } ( x _ { 1, s }, x _ { 2, s }, \dots, x _ { p, s } ) = h _ { i } ^ { \alpha } G _ { i, K } ( x _ { 1, s }, x _ { 2, s }, \dots, x _ { p, s } ) / \Gamma ( \alpha _ { i } ), \quad ( 2. 1 8 )$$

assuming

$$U _ { \alpha } ( n ) & = ( n + \alpha - 2 ) ^ { ( \alpha - 1 ) }, \\ U _ { \alpha } ( 1 ) & = ( \alpha - 1 ) ^ { ( \alpha - 1 ) } = \Gamma ( \alpha ),$$

the system Eq. (2.14) may be written as

$$x _ { i, n } = x _ { i, 0 } - \sum _ { k = 0 } ^ { n - 1 } G _ { i } ^ { 0 } ( x _ { 1, k }, x _ { 2, k }, \dots, x _ { p, k } ) U _ { \alpha _ { i } } ( n - k ). \quad \ ( 2. 2 0 )$$

With the kernel Eq. (2.19), this system turns into a fractional difference p -dimensional map. When the kernel is (see [47])

$$U _ { \alpha } ( n ) = n ^ { \alpha - 1 }, \quad U _ { \alpha } ( 1 ) = 1,$$

the system may be called a fractional p -dimensional map.

We define the map as a generalized fractional p -dimensional map when the kernel U α ( n ) belongs to the space of functions D 0 ( N 1 ) defined as (see [46])

$$\mathbb { D } ^ { i } ( \mathbb { N } _ { 1 } ) \ = \ \{ f \colon \left | \sum _ { k = 1 } ^ { \infty } \Delta ^ { i } f ( k ) \right | > N, \\ \forall N, \ \ N \in \mathbb { N }, \sum _ { \real } ^ { \infty } \left | \Delta ^ { i + 1 } f ( k ) \right | = C, \ \ C \in \mathbb { R }.$$

$$\forall N, \ N \in \mathbb { N }, \sum _ { k = 1 } ^ { \infty } \left | \Delta ^ { i + 1 } f ( k ) \right | = C, \ C \in \mathbb { R } _ { + } \}. \quad \quad ( 2. 2 2 ) \\, \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

We also assume that U α ( n ) = 0 for n &lt; 1. For any real number a , N a = { a, a +1 , a +2 , a +3 , . . . } , and N = N 0 .

## 3 Periodic Points

Now we may transform Eq. (2.20), following the same steps as in Eq. (10) from [47]. Easily verifiable result for any 0 &lt; m &lt; l is

$$& x _ { i, l N + m + 1 } - x _ { i, l N + m } \\ & = \sum _ { k = 0 } ^ { N - 1 \, l - 1 } \sum _ { j = 1 } ^ { N - 1 } G _ { i } ^ { 0 } ( x _ { 1, l N + m - l k - j }, x _ { 2, l N + m - l k - j }, \,.., x _ { p, l N + m - l k - j } ) \Big [ U _ { \alpha _ { i } } ( l k + j )$$

$$& \coprod ^ { N - 1 } _ { \alpha _ { i } } ( l k + j + 1 ) \Big { ] } + \sum _ { k = 1 } ^ { N - 1 } G _ { i } ^ { 0 } ( x _ { 1, l N + m - l k }, x _ { 2, l N + m - l k }, \dots, x _ { p, l N + m - l k } ) \Big { [ } U _ { \alpha _ { i } } ( l k ) \\ & - U _ { \alpha _ { i } } ( l k + 1 ) \Big { ] } - G _ { i } ^ { 0 } ( x _ { 1, l N + m }, x _ { 2, l N + m }, \dots, x _ { p, l N + m } ) U _ { \alpha _ { i } } ( 1 ) \\ & + G _ { i } ^ { 0 } ( x _ { 1, m }, x _ { 2, m }, \dots, x _ { p, m } ) U _ { \alpha } ( l N ) + S _ { i, 0 },$$

where S i, 0 is a sum of a finite number of elements, which tends to zero as N → ∞ . Let's assume that in the limit N → ∞ the system converges to a period l sink ( -cycle) l

$$x _ { i, l, m } = \lim _ { N \to \infty } x _ { i, N l + m }, \ \ 0 < m < l + 1, \quad \ \ ( 3. 2 )$$

and consider the limit of Eq. (3.1) as N → ∞ . Following the same steps as in Eq. (12) - Eq. (18) from [47] and using the Theorem 1, we come to the following system of ( l -1) × p equations:

$$x _ { i, l, m + 1 } - x _ { i, l, m } = \sum _ { j = 0 } ^ { m - 1 } S _ { i, j + 1, l } G _ { i } ^ { 0 } ( x _ { 1, l, m - j }, x _ { 2, l, m - j }, \dots, x _ { p, l, m - j } )$$

$$& \nu - \nu \\ + & \sum _ { j = m } ^ { l - 1 } S _ { i, j + 1, l } G _ { i } ^ { 0 } ( x _ { 1, l, m - j + l }, x _ { 2, l, m - j + l }, \dots, x _ { p, l, m - j + l } ), \\ 0 & < m < l, \ \ 0 < i \leq p,$$

$$( 3. 3 )$$

where

$$S _ { i, j + 1, l } = \sum _ { k = 0 } ^ { \infty } \left [ U _ { \alpha _ { i } } ( l k + j ) - U _ { \alpha _ { i } } ( l k + j + 1 ) \right ], \ \ 0 \leq j < l. \quad ( 3. 4 )$$

It is easy to see that

$$\sum _ { j = 1 } ^ { l } S _ { i, j, l } = 0, \ \ 0 < i \leq p.$$

As in [47] (Eq. (20) - Eq. (23)), the limits when N →∞ of the totals of the projections of the periodic points

$$& \text{$p$voices on the previous points} \\ & \sum _ { m = 1 } ^ { l } x _ { i, l N + m } = l x _ { i, 0 } \\ & - \sum _ { m = 1 } ^ { l } \sum _ { j = 1 } ^ { l } \sum _ { k = 0 } ^ { N - 1 } G _ { i } ^ { 0 } ( x _ { 1, l N + m - l k - j }, x _ { 2, l N + m - l k - j }, \dots, x _ { p, l N + m - l k - j } ) U _ { \alpha _ { i } } ( l k + j ) \\ & - \sum _ { m = 1 } ^ { l } \sum _ { j = 1 } ^ { m } G _ { i } ^ { 0 } ( x _ { 1, m - j }, x _ { 2, m - j }, \dots, x _ { p, m - j } ) U _ { \alpha _ { i } } ( l N + j ),$$

after application of the second statement of Theorem 1, give additional p equations:

$$\sum _ { j = 1 } ^ { l } G _ { i } ^ { 0 } ( x _ { 1, l, j }, x _ { 2, l, j }, \dots, x _ { p, l, j } ) = 0, \ 0 < i \leq p,$$

which together with Eq. (3.3) define l (0 &lt; m ≤ l ) p -dimensional periodic points of the l -cycle x i,l,m .

## 4 Generalized fractional H´non map e

In the Generalized fractional H´non map, the functions generating the map e are:

$$G _ { 1 H } ^ { 0 } & = - \frac { h ^ { \alpha _ { 1 } } } { \Gamma ( \alpha _ { 1 } ) } ( 1 - x _ { 1, n } - a x _ { 1, n } ^ { 2 } + x _ { 2, n } ), \\ G _ { 2 H } ^ { 0 } & = - \frac { h ^ { \alpha _ { 2 } } } { \Gamma ( \alpha _ { 2 } ) } ( b x _ { 1, n } - x _ { 2, n } ). \\ \hat { \ } \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

The fixed points of the generalized fractional H´non map are the same as the e fixed points of the classical H´non map: e

$$x _ { 1 f } & = \frac { b - 1 \pm \sqrt { ( b - 1 ) ^ { 2 } + 4 a } } { 2 a }, \\ x _ { 2 f } & = b x _ { 1 f }.$$

## 4.1 The fractional H´non map e

Dissipative 2D maps may be introduced as solutions of differential equations with dissipation and kicks. Corresponding fractional maps may be obtained as solutions of the dissipative differential equations with fractional derivatives (see [53,54]). The generalized fractional dissipative map is obtained in [54] as a solution of the equation

$$_ { 0 } D _ { t } ^ { \alpha } x ( t ) - q _ { 0 } D _ { t } ^ { \alpha - 1 } x ( t ) = K G ( x ( t ) ) \sum _ { n = 0 } ^ { \infty } \delta ( t - n T ), \ \ 1 < \alpha \leq 2. \quad ( 4. 3 )$$

Then, the fractional H´non map is obtained by assuming e

$$G ( x ) = - \frac { q } { 1 + b } [ 1 + ( 1 + b ) x + a x ^ { 2 } ]$$

and integrating Eq. (4.3).

We should note that the technique to derive fractional maps from differential equations of the orders α ≤ 1 with kicks was developed in 2013 (see [23, 24]). So, the way in which Tarasov [54] generalized the H´non map was the e most appropriate in 2010. The fractional H´non map, which is a particular e form of the generalized fractional map proposed in this paper, is a solution

of the system of two Caputo fractional differential equations of the orders 0 &lt; α i ≤ 1, i = 1 2 with kicks (see Section 2.3 in [24] and Section 3 in [30]): ,

$$\overset { C } { _ { 0 } } D _ { t } ^ { \alpha _ { i } } x _ { i } ( t ) + G _ { i, K } ( x _ { 1 } ( t - \Delta ), x _ { 2 } ( t - \Delta ) ) \sum _ { n = - \infty } ^ { \infty } \delta \left ( \frac { t } { h } - ( n + \varepsilon ) \right ) = 0, \quad ( 4. 5 )$$

where ε &gt; ∆ &gt; 0, ε → 0, 0 &lt; α i ≤ 1, α i ∈ R , index K stands for the set of parameters, and the initial conditions

$$x _ { i } ( 0 + ) = b _ { i, 0 }, \ i = 1, 2.$$

The left-sided Caputo fractional derivative C 0 D x t α t ( ) is defined for t &gt; 0 [55, 56,57] as

$$\mathcal { C } _ { 0 } D _ { t } ^ { \alpha } x ( t ) = _ { 0 } I _ { t } ^ { n - \alpha } \ D _ { t } ^ { n } x ( t ) = \frac { 1 } { \Gamma ( n - \alpha ) } \int _ { 0 } ^ { t } \frac { D _ { \tau } ^ { n } x ( \tau ) d \tau } { ( t - \tau ) ^ { \alpha - n + 1 } }, \quad ( 4. 7 )$$

where n -1 &lt; α ≤ n . The system of equations Eq. (4.5) is equivalent to the 2-D map

$$x _ { i, n + 1 } = x _ { i, 0 } - \frac { h ^ { \alpha _ { i } } } { \Gamma ( \alpha _ { i } ) } \sum _ { k = 0 } ^ { n } G _ { i, K } ( x _ { 1, k }, x _ { 2, K } ) ( n - k + 1 ) ^ { \alpha _ { i } - 1 }, \quad ( 4. 8 )$$

The fractional H´non map is obtained when e

$$G _ { 1, K } & = - ( 1 - x _ { 1, n } - a x _ { 1, n } ^ { 2 } + x _ { 2, n } ), \\ G _ { 2, K } & = - ( b x _ { 1, n } - x _ { 2, n } ),$$

which corresponds to the definition of the generalized fractional H´non map e Eq. (4.1) with

$$G _ { i H } ^ { 0 } = - \frac { h ^ { \alpha _ { 1 } } } { \Gamma ( \alpha _ { 1 } ) } G _ { i, K }, \ i = 1, 2.$$

The numerical investigation of the fractional H´ enon map is beyond the scope of this paper and will be performed elsewhere.

## 4.2 Period two points of the fractional difference H´non map e

In this section, we will assume α 1 = α 2 = α . Then (see Eq. (31) from [48]),

$$S _ { 1, 2, 2 } & = S _ { 2, 2, 2 } = \Gamma ( \alpha ) 2 ^ { - \alpha }, \\ S _ { 1, 1, 2 } & = S _ { 2, 1, 2 } = - \Gamma ( \alpha ) 2 ^ { - \alpha },$$

and the equations defining T = 2 points { ( x 1 2 1 , , , x 2 2 1 , , ) , ( x 1 2 2 , , , x 2 2 2 , , ) } may be written as

$$& x _ { 1, 2, 2 } - x _ { 1, 2, 1 } = ( h / 2 ) ^ { \alpha } \Big \{ ( x _ { 1, 2, 2 } - x _ { 1, 2, 1 } ) [ 1 + a ( x _ { 1, 2, 2 } + x _ { 1, 2, 1 } ) ] \\ & - ( x _ { 2, 2, 2 } - x _ { 2, 2, 1 } ) \Big \},$$

$$& x _ { 2, 2, 2 } - x _ { 2, 2, 1 } = ( h / 2 ) ^ { \alpha } [ ( x _ { 2, 2, 2 } - x _ { 2, 2, 1 } ) - b ( x _ { 1, 2, 2 } - x _ { 1, 2, 1 } ) ], \\ & 2 - ( x _ { 1, 2, 2 } + x _ { 1, 2, 1 } ) - a ( x _ { 1, 2, 2 } ^ { 2 } + x _ { 1, 2, 1 } ^ { 2 } ) + ( x _ { 2, 2, 2 } + x _ { 2, 2, 1 } ) = 0, \\ & b ( x _ { 1, 2, 1 } + x _ { 1, 2, 2 } ) - ( x _ { 2, 2, 1 } + x _ { 2, 2, 2 } ) = 0. & ( 4. 1 2 )$$

If we exclude the fixed points, then the solution of these equations defining two T = 2 points can be written as

$$\max \cos \nu \text{ www} \, \cos \\ x _ { 1, 2, 1 } & = \frac { C _ { 1 } \pm \sqrt { 2 C _ { 2 } - C _ { 1 } ^ { 2 } } } { 2 }, \\ x _ { 2, 2, 1 } & = \frac { b } { 2 } C _ { 1 } + \frac { b } { C _ { 0 } } ( x _ { 1, 2, 1 } - \frac { C _ { 1 } } { 2 } ), \\ x _ { 1, 2, 2 } & = C _ { 1 } - x _ { 1, 2, 1 }, \\ x _ { 2, 2, 2 } & = b C _ { 1 } - x _ { 2, 2, 1 },$$

where

$$C _ { 0 } & = 1 - ( 2 / h ) ^ { \alpha }, \\ C _ { 1 } & = \frac { b - C _ { 0 } ^ { 2 } } { a C _ { 0 } }, \\ C _ { 2 } & = \frac { 1 } { a } [ 2 + ( b - 1 ) C _ { 1 } ]. \\ \wedge & - 1 \quad \omega \quad \text{www.tu.com} \, T = 9 \text{ in to of the localization } \mathbb { m } { \text{$\pi$} }$$

When h = 1 and α = 1, we recover two T = 2 points of the classical H´non e map:

$$x _ { 1, 2 } & = \frac { 1 - b \pm \sqrt { 4 a - 3 ( 1 - b ) ^ { 2 } } } { 2 a }, \\ y _ { 1, 2 } & = b x _ { 2, 1 }.$$

The fixed point T = 2 point bifurcation (we assume x of a bifurcation diagram with the fixed b ) occurs when

$$\dots \dots \dots \dots \dots \dots \\ x _ { 1, 1 b i f } = \frac { b - C _ { 0 } ^ { 2 } } { 2 a C _ { 0 } }, \\ x _ { 2, 1 b i f } = b x _ { 1, 1 b i f }, \\ a _ { b i f } = \frac { b - C _ { 0 } ^ { 2 } } { 4 C _ { 0 } } \left [ \frac { b } { C _ { 0 } } - C _ { 0 } - 2 ( b - 1 ) \right ]. \\ \text{host of our knowledge the only manner in which bitfraction diagrams}$$

To the best of our knowledge, the only paper in which bifurcation diagrams of the fractional H´non map were drawn is the article [58]. The diagrams ( e x 1 as a function of a ) were obtained for four different values of parameters b and α including the values b = 0 05 and . α = 0 6 for which the fixed and . T = 2 points (stable and unstable) are depicted in Fig. 1. The value of a at the bifurcation point a bif = 0 243 in Fig. 1 is the same value as the value calculated from . Eq. (4.16). The two fixed points in Fig. 1 are drawn only for values a &lt; a bif because, according to the numerical calculations form [58], they are unstable for those a . More than that, it is well known that the solution of Eq. (4.2)

10

M. Edelman

Fig. 1 Fixed and T = 2 points (stable and unstable) for the fractional difference H´non e map with b = 0 05 and . α = 0 6. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00134v2-figure-20e64731-d0a3-4e98-98a9-5671b077e1fa.jpg)

with the negative sign, which does not bifurcate, is always unstable. From the comparison of the figure from [58] and Fig. 1, it is obvious that the accuracy of calculations using formulae obtained in this paper is much better than from the direct map simulations.

## 5 Fractional difference Lozi map

In the fractional difference Lozi map, the functions generating the map are:

$$G _ { 1 H } ^ { 0 } = - \frac { h ^ { \alpha _ { 1 } } } { \underset { \dots } { \Gamma ( \alpha _ { 1 } ) } } ( 1 - x _ { 1, n } - a | x _ { 1, n } | + x _ { 2, n } ),$$

$$G _ { 2 H } ^ { 0 } = - \frac { \stackrel { \theta } { h ^ { \alpha _ { 2 } } } ^ { \theta ^ { \prime } } } { \Gamma ( \alpha _ { 2 } ) } ( b x _ { 1, n } - x _ { 2, n } ).$$

The fixed points of the generalized fractional Lozi map are the same as the fixed points of the classical Lozi map:

$$x _ { 1 f } & = \frac { 1 } { 1 + a - b }, \\ x _ { 2 f } & = b x _ { 1 f }$$

if x 1 f &gt; 0 and

$$x _ { 1 f } & = \frac { 1 } { 1 - a - b }, \\ x _ { 2 f } & = b x _ { 1 f }$$

if x 1 f ≤ 0. The equations defining period two points, in the case α 1 = α 2 = α may be obtained the same way as they were obtained in Section 4.2 for the H´ enon map:

$$\text{$row$ map}. \\ x _ { 1, 2, 2 } - x _ { 1, 2, 1 } = ( h / 2 ) ^ { \alpha } \left \{ x _ { 1, 2, 2 } - x _ { 1, 2, 1 } + a ( | x _ { 1, 2, 2 } | - | x _ { 1, 2, 1 } | ) \\ - ( x _ { 2, 2, 2 } - x _ { 2, 2, 1 } ) \right \}, \\ x _ { 2, 2, 2 } - x _ { 2, 2, 1 } = ( h / 2 ) ^ { \alpha } [ ( x _ { 2, 2, 2 } - x _ { 2, 2, 1 } ) - b ( x _ { 1, 2, 2 } - x _ { 1, 2, 1 } ) ], \\ 2 - ( x _ { 1, 2, 2 } + x _ { 1, 2, 1 } ) - a ( | x _ { 1, 2, 2 } | + | x _ { 1, 2, 1 } | ) + ( x _ { 2, 2, 2 } + x _ { 2, 2, 1 } ) = 0, \\ b ( x _ { 1, 2, 1 } + x _ { 1, 2, 2 } ) - ( x _ { 2, 2, 1 } + x _ { 2, 2, 2 } ) = 0. \\ \text{$c++\cdots\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{cases}$$

Stability and bifurcations in the regular Lozi map were investigated in [59], and some numerical results for the fractional Lozi map were obtained in [60]. We will not calculate the Lozi's map periodic points and analyze their stability here - it will be done in one of the following publications.

## 6 Conclusion

In this paper we defined fractional generalization of multidimensional maps and derived equations defining their periodic points. This allows calculations of periodic points in a large set of maps, which are already used in publications, including the H´non and Lozi maps. The problems of stability and bifurcations e in these maps will be investigated in publications that will follow.

Acknowledgements The author acknowledges continuing support from Yeshiva University and expresses his gratitude to the administration of Courant Institute of Mathematical Sciences at NYU for the opportunity to perform computations at Courant.

## Conflict of interest

The author declares that he has no conflict of interest.

## References

- 1. Tarasov, V.E., Tarasova, V.V.: Economic Dynamics with Memory: Fractional Calculus Approach. De Gruyter, Berlin, Boston (2021). DOI: 10.1515/9783110627459

- 2. Tarasov, V.E., Tarasova, V.V.: Long and short memory in economics: fractional-order difference and differentiation. Int. J. Management Social Sciences 5 , 327-334 (2016). DOI: 10.21013/jmss.v5.n2.p10
- 3. Tarasova, V.V., Tarasov, V.E.: Logistic map with memory from economic model. Chaos, Solitons and Fractals 95 , 84-91 (2017)
- 4. Edelman, M.: Evolution of systems with power-law memory: Do we have to die? (Dedicated to the Memory of Valentin Afraimovich). In: Skiadas, C.H., Skiadas C. (eds.) Demography of Population Health, Aging and Health Expenditures, pp. 65-85. Springer, eBook (2020)
- 5. He, S., Sun, K., Peng, Y., Wang, L.: Wang Modeling of discrete fracmemristor and its application. AIP Advances 10 , 015332 (2020)
- 6. Huang, L.-L., Baleanu, D, Wu, G.-C., Zeng, S.-D.: A new application of the fractional logistic map. Romanian Journal of Physics 61 , 1172-1179 (2016)
- 7. Ding, D., Wang, J., Wang, M. et al. Controllable multistability of fractional-order memristive coupled chaotic map and its application in medical image encryption. Eur. Phys. J. Plus 138 , 908 (2023).
- 8. Selvam1, A.G.M., Vianny, D.A.: Discrete Fractional Order SIR Epidemic Model of Childhood Diseases with Constant Vaccination and it's Stability. Int. J. Tech. Innova. Mod. Engin. Sci. 4 , 405-410 (2018)
- 9. Ostalczyk P.: Discrete Fractional Calculus: Applications in Control and Image Processing. World Scientific, River Edge, NJ (2016)
- 10. Ortigueira, M.: Discrete-time fractional difference calculus: origins, evolutions, and new formalisms. Fractal Fract. 7 , 502 (2023)
- 11. Fulinski, A., Kleczkowski, A.S.: Nonlinear maps with memory. Physica Scripta 335 , 119-122 (1987)
- 12. Fick, E, Fick, M., Hausmann, G.: Logistic equation with memory,. Physical Review A 44 , 2469-2473 (1991) Simulating memory effects with discrete
- 13. Giona, M.: Dynamics and relaxation properties of complex systems with memory. 1Nonlinearity 14 , 911-925 (1991)
- 14. Hartwich, K., Fick, E.: Hopf bifurcations in the logistic map with oscillating memory. Physics Letters A 177 , 305-310 (1993)
- 15. Gallas, J.A.C.: dynamical systems. Physica A 195 , 417-430 (1993)
- 16. Gallas, J.A.C.: Simulating memory effects with discrete dynamical systems. Physica A 198 , 339-339 (1993) (erratum)
- 17. Edelman, M.: On the fractional Eulerian numbers and equivalence of maps with long term power-law memory (integral Volterra equations of the second kind) to Gr¨nwald-Letnikov fractional difference (differential) equau tions. Chaos 25 , 073103 (2015)
- 18. Tarasov, V.E., Zaslavsky, G.M.: Fractional equations of kicked systems and discrete maps. J. Phys. A 41 , 435101 (2008)
- 19. Stanislavsky, A.A.: Long-term memory contribution as applied to the motion of discrete dynamical system. Chaos 16 , 043105 (2006)
- 20. Edelman, M., Tarasov, V.E.: Fractional standard map. Phys. Let. A 374 , 279-285 (2009)
- 21. Edelman, M.: Fractional standard map: Riemann-Liouville vs. Caputo. Commun. Nonlin. Sci. Numer. Simul. 16 , 4573-4580 (2011)

- 22. Edelman, M., Taieb, L.A.: New types of solutions of non-linear fractional differential equations. In: Almeida, A., Castro, L, Speck F.-O. (eds.) Advances in Harmonic Analysis and Operator Theory; Series: Operator Theory: Advances and Applications 229 , 139-155. Springer, Basel (2013)
- 23. Edelman, M.: Universal fractional map and cascade of bifurcations type attractors. Chaos 23 , 033127 (2013)
- 24. Edelman, M.: Fractional maps as maps with power-law memory. In: Afraimovich, V., Luo, A.C.J.,Fu, X. (eds.) Nonlinear Dynamics and Complexity, 79-120. Springer, New York (2014)
- 25. Xiao, H., Ma, Y., Li, C.: Chaotic vibration in fractional maps. Journal of Vibration and Control 20 , 964-972 2014. DOI:10.1177/1077546312473769
- 26. Gray, H.L., Zhang, N.F.: On a new definition of the fractional difference. Mathematics of Computation 50 , 513-529 (1988)
- 27. Miller, K.S., Ross, B.: Fractional difference calculus. In: Srivastava H.M., Owa S. (eds.) Univalent Functions, Fractional Calculus, and Their Applications, 139-151. Ellis Howard, Chichester (1989)
- 28. Agarwal, R.P.: Difference equations and inequalities. Marcel Dekker, New York (2000)
- 29. Goodrich, C., Peterson, A.: Discrete Fractional Calculus. Springer, New York (2015)
- 30. Edelman, M: Maps with power-law memory: direct introduction and Eulerian numbers, fractional maps, and fractional difference maps. In: Kochubei, A. and Luchko, Yu (eds.) Handbookp of Fractional Calculus With Applications: Theory, vol. 2, pp. 47-64. De Gruyter, Berlin (2019)
- 31. Edelman, M.: Dynamics of nonlinear systems with power-law memory. In: Tarasov, V.E. (ed.) Handbook of Fractional Calculus with Applications: Applications in Physics, vol. 4, pp. 103-132. De Gruyter, Berlin (2019)
- 32. Atici, F.M., Eloe, P.W.: Initial value problems in discrete fractional calculus, Proc. Am. Math. Soc. 137 , 981-989 (2009)
- 33. Anastassiou, G.A.: Discrete Fractional Calculus and Inequalities, http://arxiv.org/abs/0911.3370 (2009)
- 34. Atici, F. M., Eloe, P. W.: Discrete fractional calculus with the nabla operator. Electron. J. Qual. Theory Differ. Equ., Spec. Ed. I. 2009, No. 3, 1-12 (2009)
- 35. Anastassiou, G.A.: Nabla discrete fractional calculus and nabla inequalities. Mathematical and Computer Modelling 51 , 562-571 (2010). DOI:10.1016/j.mcm.2009.11.006
- 36. Bastos, N. R. O., Ferreira, R. A. C., Torres, D. F. M.: Discrete-time fractional variational problems, Signal Process. 91 , 513-524 (2011)
- 37. Chen, F., Luo, X., Zhou, Y.: Existence Results for Nonlinear Fractional Difference Equation. Adv. Differ. Equ. 2011 , 713201 (2011)
- 38. Wu, G.-C., Baleanu, D., Zeng, S.-D.: Discrete chaos in fractional sine and standard maps. Phys. Lett. A 378 , 484-487 (2014)
- 39. Wu, G.-C., Baleanu, D.: (2014), Discrete fractional logistic map and its chaos. Nonlin. Dyn. 75 , 283-287 (2014)

- 40. Peng, Y., Sun, K., He, S., Wang, L.: Comments on Discrete fractional logistic map and its chaos. Nonlinear Dyn. 75, 283-287 (2014). Nonlin. Dyn. 97 , 897-901 (2019)
- 41. Edelman, M.: Caputo standard α -family of maps: Fractional difference vs. fractional. Chaos 24 , 023137 (2014)
- 42. Edelman, M.: Fractional maps and fractional attractors. Part ii: Fractional difference α -families of maps. Discontinuity, Nonlinearity, and Complexity 4 , 391-402 (2015)
- 43. Edelman, M.: On stability of fixed points and chaos in fractional systems. Chaos 28 , 023112 (2018)
- 44. Edelman, M.: Universality in systems with power-law memory and fractional dynamics. In: Edelman, M., Macau, E., Sanjuan, M.A.F. (eds.) Chaotic, Fractional, and Complex Dynamics: New Insights and Perspectives. Ser.: Understanding Complex Systems, 147-171. Springer, eBook (2018)
- 45. Edelman, M.: Cycles in asymptotically stable and chaotic fractional maps. Nonlinear Dynamics 104 , 2829-2841 (2021)
- 46. Edelman, M., Helman, A.B.: Asymptotic cycles in fractional maps of arbitrary positive orders. Fract. Calc. Appl. Anal. 25 , 181-206 (2022). DOI: 10.1007/s13540-021-00008-w
- 47. Edelman, M., Helman, A. B., Smidtaite, R.: Bifurcations and transition to chaos in generalized fractional maps of the orders 0 &lt; α &lt; 1. Chaos 33 , 063123 (2023). DOI: 10.1063/5.0151812
- 48. Edelman, M.: Stability of fixed points in generalized fractional maps of the orders 0 &lt; α &lt; 1. Nonlinear Dynamics, 111 , 10247-10254 (2023). DOI: 10.1007/s11071-023-08359-0
- 49. Elaydi, S.: An Introduction to Difference Equations. Undergraduate Texts in Mathematics. Springer, New York (2005)
- 50. Edelman, M.: Comments on A note on stability of fractional logistic maps. Appl. Math. Lett. 129 , 107892 (2022)
- 51. H´ enon, M.: A two-dimensional mapping with a strange attractor. Communications in Mathematical Physics 50 , 69-77 (1976). DOI:10.1007/BF01608556
- 52. Lozi, R.: Un attracteur estrange (?) du type attracteur de H'enon. J. Physique (Paris) 39 , (Coll. C5), no. 8, 9-10 (1978)
- 53. Edelman, M., Tarasov, V.E.: Fractional dissipative standard map. Chaos 20 , 023127 (2010). https://doi.org/10.1063/1.3443235
- 54. Tarasov, V.E.: Fractional Zaslavsky and H´non Discrete Maps. In: Luo, e A.C.J., Afraimovich, V. (eds.) Long-range Interaction, Stochasticity and Fractional Dynamics. Ser.: Nonlinear Physical Science, 1-26. Springer, Berlin, Heidelberg (2010). https://doi.org/10.1007/978-3-642-12343-6 1
- 55. Kilbas, A.A., Srivastava, H.M., Trujillo, J.J.: Theory and Application of Fractional Differential Equations. Elsevier, Amsterdam (2006)
- 56. Podlubny I.: Fractional differential equations. Academic Press, San Diego (1999)
- 57. Samko, S.G., Kilbas, A.A., Marichev, O.I.: Fractional Integrals and Derivatives Theory and Applications. Gordon and Breach, New York (1993)

- 58. Hu, T.C.: Discrete chaos in fractional Henon map. Applied Mathematics 5 , 2243-2248 (2014). http://dx.doi.org/10.4236/am.2014.515218
- 59. Botella-Soler, V. et al.: Bifurcations in the Lozi map. Journal of Physics A: Mathematical and Theoretical 44 , 305101 (2011). DOI:10.1088/17518113/44/30/305101
- 60. Khennaoui, A.-A. et al.: On fractional-order discrete-time systems: Chaos, stabilization and synchronization. Chaos Solitons Fractals 119 , 150-162 (2019). https://doi.org/10.1016/j.chaos.2018.12.019

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Mark Edelman
Asymptotic cycles in fractional generalizations of multidimensional maps
In regular dynamics, discrete maps are model presentations of discrete
dynamical systems, and they may approximate continuous dynamical systems. Maps
are used to investigate general properties of dynamical systems and to model
various natural and socioeconomic systems. They are also used in engineering.
Many natural and almost all socioeconomic systems possess memory which, in many
cases, is power-law-like memory. Generalized fractional maps, in which memory
is not exactly the power-law memory but the asymptotically power-law-like
memory, are used to model and investigate general properties of these systems.
In this paper we extend the definition of the notion of generalized fractional
maps of arbitrary positive orders that previously was defined only for maps
which, in the case of integer orders, converge to area/volume-preserving maps.
Fractional generalizations of H'enon and Lozi maps belong to the newly defined
class of generalized fractional maps. We derive the equations which define
periodic points in generalized fractional maps. We consider applications of our
results to the fractional and fractional difference H'enon and Lozi maps.