## Stability and Bifurcation Analysis of Discrete Dynamical Systems 2020

View this Special IssueResearch Article | Open Access

Mirela Garić-Demirović, Samra Moranjkić, Mehmed Nurkanović, Zehra Nurkanović, "Stability, Neimark–Sacker Bifurcation, and Approximation of the Invariant Curve of Certain Homogeneous Second-Order Fractional Difference Equation", *Discrete Dynamics in Nature and Society*, vol. 2020, Article ID 6254013, 12 pages, 2020. https://doi.org/10.1155/2020/6254013

# Stability, Neimark–Sacker Bifurcation, and Approximation of the Invariant Curve of Certain Homogeneous Second-Order Fractional Difference Equation

**Academic Editor:**Abdul Qadeer Khan

#### Abstract

We investigate the local and global character of the unique equilibrium point and boundedness of the solutions of certain homogeneous fractional difference equation with quadratic terms. Also, we consider Neimark–Sacker bifurcations and give the asymptotic approximation of the invariant curve.

#### 1. Introduction and Preliminaries

In this paper, the subject of our consideration is the following difference equation:with positive parameters , and where initial conditions are positive numbers. By substituting , equation (1) reduces to the following equation:

Equations (1) and (2) are the special cases of the following homogeneous rational difference equation:

Equation (3) is a general homogeneous rational difference equation with quadratic terms and it is very complicated for investigation in many special cases. The function associated with the right side of equation (3) is of the form , with the following monotonicity property: it is either monotonically decreasing in the first variable and monotonically increasing in the second one or monotonically increasing in the first variable and monotonically decreasing in the second variable. Using the theory of monotone maps, it was possible to investigate the global dynamics of some special cases of equation (3) [1–4] in the situations when corresponding function is monotonically decreasing in the first variable and monotonically increasing in the second variable.

In [4], the authors investigated the local and global character of the unique equilibrium point and the existence of Neimark–Sacker and period-doubling bifurcations of equation (3). They have also studied the local and global stability of the minimal period-two solutions for some special cases of the parameters. The special case when and , i.e.,was considered in [2]. It was shown that equation (4) is characterized by three types of global behavior with respect to the existence of a unique positive equilibrium and existence of one or two minimal period-two solutions, one of which is locally asymptotically stable and the other is a saddle point. An important feature of this equation is the coexistence of an equilibrium and the minimal period-two solution which are both locally asymptotically stable. Also, the basins of attraction of these solutions are described in detail.

In [3], a special case of equation (3) was considered when , i.e.,

It was shown that (5) exhibits period-two bifurcation and that stable manifolds of the minimal period-two solutions represent the boundaries of the basin of attraction of locally stable equilibrium point and the basins of attraction of points and . Further, in the situations when equilibrium is a saddle point, corresponding stable manifold separates the basins of attraction of the points and .

The investigation of the special cases of equation (3) when corresponding function is monotonically increasing in the first variable and monotonically decreasing in the second variable is significantly harder, since by now, there is no any general result for this type of monotonicity. Some initial steps, regarding this, were taken in [5], where the considered case was , , , i.e.,

In [5], the authors have successfully used the embedding method to demonstrate the boundedness of the solutions, and then they determined the invariant interval of equation (6). That was a crucial idea for proving that local asymptotic stability (which holds when ) implies global asymptotic stability of the unique positive equilibrium point when . Also, the existence of Neimark–Sacker bifurcation is shown and asymptotic approximation of the invariant curve is computed.

The investigation of the dynamics of equation (4) and its special cases has been the subject of many research studies for the last ten years. Some of these cases, as we have seen, have been successfully realized. However, only the ideas from [5] finally made the problem of investigating the behavior of equation (1) or equation (2) solvable. Namely, note that equation (2) has the formwhere , , and . Therefore, we will consider equation (7) instead of equation (1). However, the application of the embedding method on equation (7) is significantly harder compared with application on equation (6) (because the form is more complicated). We will investigate local and global stability of a unique equilibrium point and boundedness of the solutions of equation (7) and examine the existence of Neimark–Sacker bifurcation. Also, in the situation when Neimark–Sacker bifurcation appears, we will give the asymptotic approximation of the invariant curve.

The special case of equation (3), when and , was considered in [1]. In the region of parameters where , corresponding function is monotonically increasing by first variable and monotonically decreasing by second variable, and through very complicated calculations, using the so-called “M-m” theorems, it is shown that in some areas of parametric space of parameters and , unique positive equilibrium is globally stable. The existence of period-doubling bifurcation is proved in the case , when corresponding function is monotonically decreasing by first variable and monotonically increasing by second variable. Using the theory of monotone maps, the global stability of minimal period-two solution is shown for some special values of parameters.

Notice that equation (3) is a special case (and probably the most complicated equation of the form (3, 3)) of the following general second-order rational difference equation with quadratic termsthat has caught the attention of mathematical researchers over the last ten years ([1–12]).

The following lemma gives us the type of local stability of a unique positive equilibrium point of equation (7) depending on different values of parameters and .

Lemma 1. *Equation (7) has a unique equilibrium point , which is*(a)*locally asymptotically stable if *(b)*nonhyperbolic if *(c)*a repeller if *

*Proof. *Denote . The linearized equation associated with equation (7) about the equilibrium point has the formwhere and . NoticeSincethe conclusion follows.

This paper is organized as follows. In Section 2, using the embedding method [5, 13] and the so-called “M-m” theorems [14–17], we prove global asymptotic stability of a unique positive equilibrium for and conduct the semicycle analysis as well. In Section 3, using Neimark–Sacker theorem [15, 18–21], we give reduction to the normal form and perform computation of the coefficients of the aforementioned bifurcation, based on the computational algorithm developed in [22]. Furthermore, we determine the asymptotic approximation of the invariant curve and give a visual evidence.

#### 2. Global Asymptotic Stability

In this section, we will show that all solutions of equation (7) are bounded, and using the so-called “M-m” theorem, we will obtain sufficient conditions for unique positive equilibrium to be globally asymptotically stable. Similar to [5, 13], we apply the method of embedding. First, we substitutein equation (7) and obtain

Then, by substitutingin equation (13), we have

Notice that the solutions of equation (15) are bounded:that is,

Since every solution of equation (7) is also a solution of equation (15) with initial values , , and , we see that the solutions of equation (7) are also bounded.

From (7) and (15), we getwhich impliesthat is,

By replacing in (20), we obtain the following equation:that is,

*Remark 1. *Note that equation (22) has the unique equilibrium point , which is the same equilibrium as in equation (7). It follows fromFurthermore, notice that every solution of equation (7) is also a solution of equation (22) with initial values , and and that it is of the form , where

Lemma 2. *Every interval of the form , whereis an invariant interval for the function .*

*Proof. *As we know, for , interval is invariant for the function ifFor , we have that . If satisfiesthen we obtain . It further implies that for every , there exists such , which means that is invariant for the function , where . Since we can assume that and , is the invariant for .

Lemma 3. *Interval is an attracting interval for equation (22).*

*Proof. *It is clear that we need to show that every solution of equation (22) must enter interval . Note that for arbitrary initial conditions , , and , it holdsIf , , then for , by Lemma 2. Otherwise, if , let us prove that there must be some such that for all . Namely, suppose that for arbitrary initial conditions . Then,So, from equation (22), we obtainwhich impliesthat is,By induction, we concludeSince for , the right side in (33) is decreasing sequence converging toThen, from (33), we getSince the case is not possible (otherwise there would exist another positive equilibrium different from ), it implies that there is some such thatfor all , i.e., every solution of equation (22) must enter the interval .

Theorem 1. *If , then the equilibrium point of equation (7) is globally asymptotically stable.*

*Proof. *It is enough to prove that is an attractor of equation (22). Since there exists the invariant and attracting interval when , we need to check the conditions of Theorem A.0.5 in [14]:i.e.,By subtracting the second equation in (38) from the first, we getfrom which or for ,If , thenSince for , it implies that . Therefore,By substituting (42) into (38), we obtain the following quadratic equation:whose discriminant has the formIt is clear that if and that there are no real solutions for equation (43). If , then , which is a contradiction with assumption . When , equation (43) has two positive roots ; it implies that (for example, , for , and , where is in invariant interval ). Therefore, the conditions of Theorem A.0.5 in [14] are satisfied for and every solution of equation (22) converges to . By Remark 1, the unique equilibrium of equation (7) is an attractor. By using Lemma 1, we conclude that is globally asymptotically stable (see Figure 1).

**(a)**

**(b)**

Conjecture 1. *If , then the equilibrium point is globally asymptotically stable.*

#### 3. Neimark–Sacker Bifurcations

The following results are obtained by applying the algorithm from Theorem 1 and Corollary 1 in [5] (see also [22]). If we make a change of variable , we will shift the equilibrium point to the origin. Then, the transformed equation is given by

Setand write equation (1) in the equivalent form:

Let denote the corresponding map defined by

Then, the Jacobian matrix of is given by

The eigenvalues of are and where

Lemma 4. *If , then has equilibrium point at and eigenvalues of Jacobian matrix of at are and where**Moreover, satisfies the following:*(i)* for *(ii)* at *(iii)*Eigenvectors associated to the are**such that , , and , where .*

*Proof. *Let . Then, we obtainAfter straightforward calculation for , we obtain andfrom which it follows that for . Furthermore, we getIt is easy to see that and .

Let , where is a sufficiently small parameter. From Lemma 4, we can transform system (47) into the normal formand there are smooth functions and so that in polar coordinates, the function is given by

Now, we compute following the procedure in [22]. Notice that if and only if . First, we compute and defined in [22]. For , we havewhere

Hence, for , system (47) is equivalent to

Define the basis of by , where ; then, we can represent as

Let

Now, we have

Sincethen

By using , , and , we have

Finally, we get

If is fixed point of , then invariant curve can be approximated bywhere

Thus, we prove the following result.

Theorem 2. *Let ; then, there is a neighborhood of the equilibrium point and such that for and , -limit set of solution of equation (7) is the equilibrium point if and belongs to a closed invariant curve encircling the equilibrium point if . Furthermore, and invariant curve can be approximated by*

Since , a nondegenerate Neimark–Sacker bifurcation occurs at the critical value . We proved and , so it implies if , then an attracting closed curve exists, surrounding the unstable fixed point, when the parameter crosses the bifurcation value (supercritical Neimark–Sacker bifurcation). As increases, the attracting closed curve decreases in size and merges with the fixed point at , leaving a stable fixed point (subcritical Neimark–Sacker bifurcation). All orbits starting outside or inside the closed invariant curve, except at the origin, tend to the attracting closed curve.

The asymptotic approximation of the invariant curve is shown in Figure 2, and some orbits and trajectories in the case where the unique equilibrium point is stable or nonhyperbolic are given in Figures 3 and 4. Using parameter value and decreasing the dynamical parameter , the unique positive equilibrium point loses its stability via Neimark–Sacker bifurcation leading to chaos as depicted in Figure 5(a). Also, the Lyapunov exponents corresponding to Figure 5(a) are shown in Figure 5(b), which verifies the existence of chaos after the occurrence of Neimark–Sacker bifurcation.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

#### 4. Conclusion

The investigation of the dynamic stability of homogeneous difference equation (3) and all its special cases is very complicated. The corresponding function associated with the right side of equation (3) is monotonically increasing in the first variable and monotonically decreasing in the second variable or monotonically decreasing in the first variable and monotonically increasing in the second variable or it switches its type of monotonicity between the first and second case or vice versa, depending on the parameters which appear in the function.

The theory of monotone maps (or more precisely, the theory of competitive maps) was used for determining the dynamics of equation (3) or some special case of equation (3), in the scenario when corresponding function is monotonically decreasing in the first variable and monotonically increasing in the second variable (see [2, 3]).

The investigation of the special cases of equation (3) when corresponding function is monotonically increasing in the first variable and monotonically decreasing in the second variable is significantly harder because there is no general result for this type of monotonicity.

In every situation, when the corresponding equation does not possess minimal period-two solutions, the global stability of a unique equilibrium usually can be determined by applying the so-called “M-m” theorems and finding an invariant interval of the map before that, of course. However, it is very often impossible or extremely complicated to conduct, as we saw in [1]. That is the case with equation (1) which does not possess minimal period-two solution (since the corresponding function is monotonically increasing by first and monotonically decreasing by second variable). So, for that reason, we studied equation (7) instead of equation (1). By using the method of embedding, we were able to connect equation (7) with equation (15). Namely, we have shown that every solution of equation (7) is also a solution of equation (15) with initial conditions , , , and . Furthermore, we have shown that every solution of equation (15) is bounded, which implies that every solution of equation (7) is also bounded. After that, we linked equation (15) with equation (22) and showed that every solution of equation (7) is also a solution of equation (22) with initial conditions , , and . Additionally, we determined the invariant and attracting interval for the function that is associated with the right side of equation (22) and successfully applied “M-m” theorem to get conditions for parameters and under which the equilibrium of equation (22) and therefore of equation (7) is globally asymptotically stable. The area of the regions in plane where unique equilibrium is locally asymptotically stable and is not globally asymptotically stable is small (see Figure 1). We expect to prove Conjecture 1 in some of our future studies.

Finally, using Neimark–Sacker theorem [15, 18–21], we gave reduction to the normal form and performed computation of the coefficients of bifurcation, based on the computational algorithm developed in [22]. Furthermore, we determined the asymptotic approximation of the invariant curve and provided visual evidence (see Figures 2–5). Also, based on the computational algorithm in [23], we calculated the Lyapunov exponents corresponding to Figure 5(a) to confirm the existence of chaos after the occurrence of Neimark–Sacker bifurcation as in [24, 25] (see Figure 5(b)).

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This study was supported in part by the Fundamental Research Funds of Bosnia and Herzegovina (FMON no. 01-6211-1-IV/19).

#### References

- M. Garić-Demirović, S. Hrustić, and M. Nurkanović, “Stability and periodicity of certain homogeneous second-order fractional difference equation with quadratic terms,”
*Advances in Dynamical Systems and Applications*, vol. 14, no. 2, pp. 149–178, 2019. View at: Publisher Site | Google Scholar - M. Garić-Demirović, M. R. S. Kulenović, and M. Nurkanović, “Global dynamics of certain homogeneous second-order quadratic fractional difference equations,”
*The Scientific World Journal*, vol. 2013, Article ID 210846, p. 10, 2013. View at: Publisher Site | Google Scholar - M. Garić-Demirović, M. R. S. Kulenović, and M. Nurkanović, “Basins of attraction of certain homogeneous second order quadratic fractional difference equation,”
*Journal of Concrete and Applicable Mathematics*, vol. 13, no. 1-2, pp. 35–50, 2015. View at: Google Scholar - M. Garić-Demirović, M. Nurkanović, and Z. Nurkanović, “Stability, periodicity and Neimark-Sacker bifurcation of certain homogeneous fractional difference equation,”
*International Journal of Difference Equations*, vol. 12, no. 1, pp. 27–53, 2017. View at: Google Scholar - T. Khyat, M. R. S. Kulenović, and E. Pilav, “The Naimark-Sacker bifurcation and asymptotic approximation of the invariant curve of a certain difference equation,”
*Journal of Computational Analysis and Applications*, vol. 23, no. 8, pp. 1335–1346, 2017. View at: Google Scholar - S. Jašarević Hrustić, M. R. S. Kulenović, and M. Nurkanović, “Global dynamics and bifurcations of certain second order rational difference equation with quadratic terms,”
*Qualitative Theory of Dynamical Systems*, vol. 15, no. 1, pp. 283–307, 2016. View at: Publisher Site | Google Scholar - S. Kalabušić, M. Nurkanović, and Z. Nurkanović, “Global dynamics of certain mix monotone difference equation,”
*Mathematics*, vol. 6, no. 1, p. 10, 2018. View at: Publisher Site | Google Scholar - M. R. S. Kulenović, S. Moranjkić, M. Nurkanović, and Z. Nurkanović, “Global asymptotic stability and Neimark-Sacker bifurcation of certain mix monotone difference equation,”
*Discrete Dynamics in Nature and Society*, vol. 2018, Article ID 7052935, 22 pages, 2018. View at: Publisher Site | Google Scholar - M. R. S. Kulenovic, S. Moranjkić, and Z. Nurkanović, “Naimark-Sacker bifurcation of second order rational difference equation with quadratic terms,”
*The Journal of Nonlinear Sciences and Applications*, vol. 10, no. 7, pp. 3477–3489, 2017. View at: Publisher Site | Google Scholar - M. R. S. Kulenović, S. Moranjkić, and Z. Nurkanović, “Global dynamics and bifurcation of perturbed Sigmoid Beverton-Holt difference equation,”
*Mathematical Methods in the Applied Sciences*, vol. 39, pp. 2696–2715, 2016. View at: Publisher Site | Google Scholar - S. Moranjkić and Z. Nurkanović, “Local and global dynamics of certain second-order rational difference equations containing quadratic terms,”
*Advances in Dynamical Systems and Applications*, vol. 12, no. 2, pp. 123–157, 2017. View at: Google Scholar - H. Sedaghat, “Global behaviours of rational difference equations of orders two and three with quadratic terms,”
*Journal of Difference Equations and Applications*, vol. 15, no. 3, pp. 215–224, 2009. View at: Publisher Site | Google Scholar - E. J. Janowski and M. R. S. Kulenović, “Attractivity and global stability for linearizable difference equations,”
*The Journal of Computational and Applied Mathematics*, vol. 57, no. 9, pp. 1592–1607, 2009. View at: Publisher Site | Google Scholar - M. R. S. Kulenović and G. Ladas,
*Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures*, Chapman and Hall/CRC, Boca Raton, FL, USA, 2001. - M. R. S. Kulenović and O. Merino,
*Discrete Dynamical Systems and Difference Equations with Mathematica*, Chapman and Hall/CRC, Boca Raton, FL, USA, 2000. - M. R. S. Kulenović and M. Nurkanović, “Asymptotic behavior of a two dimensional linear fractional system of difference equations,”
*Radovi Matematički (Sarajevo Journal of Mathematics)*, vol. 11, no. 1, pp. 59–78, 2002. View at: Google Scholar - M. R. S. Kulenović and M. Nurkanović, “Asymptotic behavior of a system of linear fractional difference equations,”
*Journal of Inequalities and Applications*, vol. 2005, no. 2, pp. 127–143, 2005. View at: Publisher Site | Google Scholar - J. K. Hale and H. Kocak, “Dynamics and bifurcations,” in
*Text in Applied Mathematics*, vol. 3, Springer-Verlag, New York, NY, USA, 1991. View at: Google Scholar - Y. Kuznetsov,
*Elements of Applied Bifurcation Theory*, Springer, New York, NY, USA, 1998. - S. Robinson,
*Stability, Symbolic Dynamics and Chaos*, CRC Press, Boca Raton, FL, USA, 1995. - S. Wiggins,
*Introduction to Applied Nonlinear Dynamical Systems and Chaos, Text in Applied Mathematics*, vol. 2, Springer-Verlag, New York, NY, USA, Second edition, 2003. - K. Murakami, “The invariant curve caused by Neimark-Sacker bifurcation,”
*Dynamics of Continuous, Discrete and Impulsive Systems*, vol. 9, pp. 121–132, 2002. View at: Google Scholar - M. Sandri, “Numerical calculation of Lyapunov exponents,”
*The Mathematica Journal*, vol. 6, pp. 78–84, 1996. View at: Google Scholar - A. Al-khedhairi, A. E. Matouk, and S. S. Askar, “Bifurcations and chaos in a novel discrete economic system,”
*Advances in Mechanical Engineering*, vol. 11, no. 4, pp. 1–15, 2019. View at: Publisher Site | Google Scholar - A. E. Matouk, A. A. Elsadany, and B. Xin, “Neimark-Sacker bifurcation analysis and complex nonlinear dynamics in a heterogeneous quadropoly game with an isoelastic demand function,”
*Nonlinear Dynamics*, vol. 89, no. 4, pp. 2533–2552, 2017. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2020 Mirela Garić-Demirović et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.