Abstract
Let be a nonnegative integer and a real number greater than or equal to 1. We present qualitative global behavior of solutions to a rational nonlinear higher-order difference equation with positive initial values , and show the global asymptotic stability of its positive equilibrium solution.
1. Introduction
Difference equations have wide applications in biology, computer science, digital signal processing, and economics. A general solution structure exists for linear difference equations [1]. However, in various situations of nonlinear higher-order difference equations, solution properties can only be observed by numerical simulation, and it is often exceedingly difficult to give a full mathematical proof for the properties predicted by numerical simulation and the conclusions formed on the basis of guesswork [2]. It is, therefore, important to make qualitative analysis on nonlinear higher-order difference equations, which is the topic of the current study. There have been some related studies on rational nonlinear difference equations in the literature (see, e.g., [3–8]). Global asymptotic properties of spectral functions are also crucial in determining algebro-geometric solutions of soliton equations (see, e.g., [9, 10]) and scattering data in matrix spectral problems (see, e.g., [11]).
An iterative algorithm to approximate a zero of a given function reads and an application of this to a quadratic function , gives Let be a nonnegative integer and a real number greater than or equal to 1. We would like to consider a more general rational nonlinear higher-order difference equation with positive initial values , which engender positive solutions. We take a transformation and then obtain another difference equation Obviously, the equilibrium solution of the rational nonlinear difference equation (3), , becomes the equilibrium solution of the transformed difference equation (5), .
If we further take , then we obtain the nonlinear difference equation discussed in [6, 7]: Introducing into (6) yields where . When , this gives the nonlinear difference equation in (2). The equation (7) in the case of was studied in [5] and its closed-form solution was presented in [6]. In the general case of , the asymptotic stability of the positive equilibrium solution of the equation (7) was proved in [7].
It is direct to see that the rational nonlinear higher-order difference equation, defined by (3), possesses three equilibria: . In this article, we would like to explore global behavior of solutions to the rational nonlinear higher-order difference equation (3), show the global asymptotic stability of its positive equilibrium solution , and present two illustrative examples of positive solutions.
2. Global Behavior
2.1. Classification of Solutions
First of all, based on the rational difference equation (3), one can have
Further from (9) and (10), we can easily derive the following solution properties.
Theorem 1. If is a solution to the rational nonlinear difference equation (3), then one has where .
If , the rational difference equation (3) becomes a first-order difference equation Then for , one has , because . For , one has , since , due to . Therefore, every solution decays to , when .
Generally, there are three types of solutions to the rational nonlinear higher-order difference equation (3).
Theorem 2. Let . If is a solution to the rational nonlinear higher-order difference equation (3), and then
(a) it is eventually equal to , more precisely , which occurs when for some ;
(b) it is eventually less than , more precisely , , which occurs when for some ; or
(c) it oscillates about , possessing at most consecutive increasing terms less than and at most consecutive decreasing terms greater than .
Proof. Equality (8) and property (11) directly tell that we have three types of solutions to the rational nonlinear higher-order difference equation (3).
The decreasing and increasing characteristics of oscillatory solutions in the third solution situation (c) can be proved as follows.
Suppose that are two integers satisfying . We express where can be written as by (9).
If for , then each term in is less than zero, and so , due to (14). If for , then each term in is greater than zero, and so , due to (14). This completes the proof.
Note that based on (8), we can see that there is no solution situation that a solution of (3) is eventually greater than .
2.2. Global Asymptotic Stability
When , the equilibrium solution of the first-order rational difference equation (13) is globally asymptotically stable, since it is a globally attractive equilibrium solution of a first-order difference equation (see [12] for a general theory).
For a general , we can show the same global asymptotic stability of the positive equilibrium solution of the rational nonlinear difference equation (3), by establishing the local asymptotic stability and the global attractivity, which imply the global asymptotic stability [2]. Instead, we establish a strong negative feedback property [13] to guarantee the global asymptotic stability of (see [14] for details on the strong negative feedback property).
Theorem 3. The positive equilibrium solution of the rational nonlinear higher-order difference equation (3) is globally asymptotically stable.
Proof. Based on the rational nonlinear difference equation (3), one can have From this equality and the equality in (10), we can obtain which leads to a strong negative feedback property: with equality for all if and only if . It, therefore, follows from a stability theorem (Corollary 3 of [14]) that the equilibrium solution of the rational nonlinear difference equation (3) is globally asymptotically stable. Thus, the proof is finished.
The above theorems with gives the results in [4] (), [5] () and [7] (a general ). There have also been similar studies on polynomial difference equations (see, for example, [15]) and other studies on positive rational function solutions, called lumps, to partial differential equations (see, e.g., [16]).
2.3. Illustrative Examples and an Open Question
To illustrate the global properties stated in Theorems 2 and 3, here we present two illustrative examples associated with two special cases:in Figure 1. From the plot pictures, we see that the convergence is achieved very fast in both cases.
Finally, let . For an oscillatory solution of the rational nonlinear difference equation (3), we define Since is oscillatory, Theorem 2 guarantees that both and have infinitely many numbers. A basic open question that we are very interested in is if is decreasing on and increasing on . We point out that through the above two examples, we failed to find any counterexample to this statement, but found that two cases could occur: either but or but for some .
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was in part supported by NNSFC under Grants 11371326, 11371086, and 11571079, NSF under Grant DMS-1664561, and the distinguished professorships of King Abdulaziz University, Saudi Arabia and North-West University, South Africa.