Global Analysis of Solutions of a New Class of Rational Difference Equation

The study suggests asymptotic behavior of the solution to a new class of difference equations: . where a, bi, α and β are positive real numbers for i = 0, 1, · · · , k , and the initial conditions ψ-j, ψ-j+1, · · ·, ψ0 are randomly positive real numbers where j = 2k + 1. Accordingly, we consider the stability, boundedness and periodicity of the solutions of this recursive sequence. Indeed, we give some interesting counter examples in order to verify our strong results.

1. Introduction

The global asymptotic behavior of the solutions and oscillation of solution are two such qualitative properties which are very important for applications in many areas such as control theory, mathematical biology, neural networks, etc. It is impossible to use computer based (numerical) techniques to study the oscillation or the asymptotic behavior of all solutions of a given equation due to the global nature of these properties. Therefore, these properties have received the attention of several mathematicians and engineers.

Currently, much attention has given to study the properties of the solutions of the recursive sequences from scientists in various disciplines. Specifically, the topics dealt with include the following:

- Finding equilibrium points for the recursive sequences;

- Investigating the local stability of the solutions of the recursive sequences;

- Finding conditions which insure that the solutions of the recursive sequences are bounded;

- Investigating the global asymptotic stability of the solutions of the recursive sequences;

- Finding conditions which insure that the solutions of equation are periodic with positive prime period two or more;

- Finding conditions for oscillation of solutions.

Closely related global convergence results were well-gained from these articles  -  . Khuong in  studied the dynamics the recursive sequences

${\psi }_{\eta +1}=a+{\left(\frac{{\psi }_{\eta -k}}{{\psi }_{\eta -r}}\right)}^{p}.$

For further related and special cases of this difference equations see    ,    .

Elsayed  studied the periodicity, the boundedness of the positive solution of the recursive sequences

${\psi }_{\eta +1}=a+\frac{b{\psi }_{\eta -l}+c{\psi }_{\eta -k}}{c{\psi }_{\eta -l}+d{\psi }_{\eta -k}}.$

Abdelrahman  considered analytical investigation of the solution of the recursive sequence

${\psi }_{\eta +1}=a{\psi }_{\eta -k}+\delta {\left(\frac{b{\psi }_{\eta -l}}{c{\psi }_{\eta -k}+d{\psi }_{\eta -l}}\right)}^{\alpha }.$

By new method, Elsayed  investigated the periodic solution of the equation

${\psi }_{\eta +1}=\alpha +\beta \frac{{\psi }_{\eta }}{{\psi }_{\eta -1}}+\gamma \frac{{\psi }_{\eta -1}}{{\psi }_{\eta }}.$

Also, Moaaz  completed the results of  .

In this work, we deal with some qualitative behaviour of the solutions of the recursive sequence

${\psi }_{\eta +1}=a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{{\psi }_{\eta -\left(2i+1\right)}}{\alpha {\psi }_{\eta -2i}+\beta {\psi }_{\eta -\left(2i+1\right)}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\eta =0,1,2,\cdots$ (1.1)

where $a,{b}_{i},\alpha$ and $\beta$ are positive real numbers for $i=0,1,\cdots ,k$ , and the the initial conditions ${\psi }_{-j},{\psi }_{-j+1},\cdots ,{\psi }_{0}$ are arbitrary positive real numbers where $j=2k+1$ .

In the next, we will and to many of the basic concepts. Before anything, the concept of equilibrium point is essential in the study of the dynamics of any physical system. A point $\stackrel{¯}{\psi }$ in the domain of the function $\Phi$ is called an equilibrium point of the equation

${\psi }_{\eta +1}=\Phi \left({\psi }_{\eta },{\psi }_{\eta -1},\cdots ,{\psi }_{\eta -k}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\eta =0,1,2,\cdots$ (1.2)

if $\stackrel{¯}{\psi }$ is a fixed point of $\Phi$ [ $\Phi \left(\stackrel{¯}{\psi },\stackrel{¯}{\psi },\cdots ,\stackrel{¯}{\psi }\right)=\stackrel{¯}{\psi }$ ]. For a stability of equilibrium point, equilibrium point $\stackrel{¯}{\psi }$ of equation (2) is said to be locally stable if for all $\epsilon >0$ there exists $\delta >0$ such that, if ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ with ${\sum }_{i=0}^{k}|{\psi }_{-i}-\stackrel{¯}{\psi }|<\delta$ . As well, $\stackrel{¯}{\psi }$ is said to be locally asymptotically stable if it is locally stable and there exists $\gamma >0$ such that, if ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ with ${\sum }_{i=0}^{k}|{\psi }_{-i}-\stackrel{¯}{\psi }|<\gamma$ , then ${\mathrm{lim}}_{\eta \to \infty }{\psi }_{\eta }=\stackrel{¯}{\psi }$ . Also, $\stackrel{¯}{\psi }$ is said to be a global attractor if used for every ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ , we have ${\mathrm{lim}}_{\eta \to \infty }{\psi }_{\eta }=\stackrel{¯}{\psi }$ . On the other hand, $\stackrel{¯}{\psi }$ is said to be unstable if it is not locally stable.

Finally, Equation (1.2) is called permanent and bounded if there exists numbers r and R with $0 such that for any initial conditions ${\psi }_{-\nu }\in \left(0,\infty \right)$ for $\nu =0,1,\cdots ,k$ there exists a positive integer N which depends on these initial conditions such as $r<{\psi }_{\eta } for all $\eta \ge N$ .

The linearized equation of Equation (1.1) about the equilibrium point $\stackrel{^}{\psi }$ is

${y}_{\eta +1}=\underset{i=0}{\overset{k}{\sum }}{p}_{i}{y}_{\eta -i}$ (1.3)

where

${p}_{i}=\frac{\partial F}{\partial {\psi }_{\eta -i}}\left(\stackrel{^}{\psi },\stackrel{^}{\psi },\cdots ,\stackrel{^}{\psi }\right).$

Theorem 1.1.  Assume that ${p}_{i}\in ℝ$ for $i=0,1,\cdots ,k$ . The $+ve$ equilibrium of (1.1) is locally asymptotically stable if

$|{p}_{0}|+|{p}_{1}|+\cdots +|{p}_{k}|<1.$ (1.4)

2. Local Stability of Equation (1.1)

The $+ve$ equilibrium point of Equation (1.1) is

$\stackrel{^}{\psi }=a+\underset{i=0}{\overset{k}{\sum }}\frac{{b}_{i}\stackrel{^}{\psi }}{\alpha \stackrel{^}{\psi }+\beta \stackrel{^}{\psi }},$

and so,

$\stackrel{^}{\psi }=a+\frac{B}{\alpha +\beta },$

where

$B=\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}.$

Let $f\in C\left({\left(0,\infty \right)}^{2k+2},\left(0,\infty \right)\right)$ defined by

$f\left({u}_{0},{u}_{1},\cdots ,{u}_{2k+1}\right)=a+\underset{i=0}{\overset{k}{\sum }}\frac{{b}_{i}{u}_{2i+1}}{\alpha {u}_{2i}+\beta {u}_{2i+1}}.$ (2.1)

Therefore it follows that

$\frac{\partial f}{\partial {u}_{2r}}=-\frac{{b}_{r}\alpha {u}_{2r+1}}{{\left(\alpha {u}_{2r}+\beta {u}_{2r+1}\right)}^{2}}$ (2.2)

and

$\frac{\partial f}{\partial {u}_{2r+1}}=\frac{{b}_{r}\alpha {u}_{2r}}{{\left(\alpha {u}_{2r}+\beta {u}_{2r+1}\right)}^{2}},$ (2.3)

for $r=0,1,\cdots ,k$ .

Theorem 2.1. Let $\stackrel{^}{\psi }$ be $+ve$ equilibrium of Equation (1.1). If

$\left(\alpha -\beta \right)B

than $\stackrel{^}{\psi }$ is locally stable.

Proof. From (2.2) to (2.3), we obtain

$\frac{\partial f}{\partial {u}_{2r}}\left(\stackrel{^}{\psi },\cdots ,\stackrel{^}{\psi }\right)=-\frac{{b}_{r}\alpha }{\left(\alpha +\beta \right)\left(a\left(\alpha +\beta \right)+B\right)}={p}_{2r},$

and

$\frac{\partial f}{\partial {u}_{2r+1}}\left(\stackrel{^}{\psi },\cdots ,\stackrel{^}{\psi }\right)=\frac{{b}_{r}\alpha }{\left(\alpha +\beta \right)\left(a\left(\alpha +\beta \right)+B\right)}={p}_{2r+1},$

for $r=0,1,\cdots ,k$ . Thus, the linearized equation of (1.1) is

${y}_{\eta +1}={p}_{0}{y}_{\eta }+{p}_{1}{y}_{\eta -1}+\cdots +{p}_{2k+1}{y}_{\eta -\left(2k+1\right)}.$

It follows by Theorem 1.1 that Equation (1.1) is locally stable if

$|\frac{{b}_{0}\alpha }{\left(\alpha +\beta \right)}|+|\frac{{b}_{0}\alpha }{\left(\alpha +\beta \right)\rho }|+\cdots +|\frac{{b}_{k}\alpha }{\left(\alpha +\beta \right)\rho }|+|\frac{{b}_{k}\alpha }{\left(\alpha +\beta \right)\rho }|<1,$

where $\rho =\left(a\left(\alpha +\beta \right)+B\right)$ , and hence,

$\frac{2\alpha }{\left(\alpha +\beta \right)\rho }B<1.$

Thus, we find

$2\alpha B<\left(\alpha +\beta \right)\rho ,$

and so,

$\left(\alpha -\beta \right)B

Hence, the proof is complete.

In order to verify and support our theoretical outcomes and discussions, in this concern, we investigate several interesting numerical examples.

Example 2.1. By Theorem 2.1, the $+ve$ equilibrium Equation (1.1) with $a=2$ , $k=2$ , ${b}_{i}=0.3$ , $\beta =0.1$ and $\alpha =1$ , is locally stable (see Figure 1).

3. Global Stability of Equation (1.1)

In the following theorem, we check into the global stability of the recursive sequence (1.1).

Theorem 3.1. The $+ve$ equilibrium $\stackrel{^}{\psi }$ of Equation (1.1) is global attractor if

$B=a\left(\alpha -\beta \right).$

Figure 1. The stable solution corresponding to difference Equation (1.1).

Proof. We consider the function as follow:

$f\left({u}_{0},{u}_{1},\cdots ,{u}_{2k+1}\right)=a+\underset{i=0}{\overset{k}{\sum }}\frac{{b}_{i}{u}_{2i+1}}{\alpha {u}_{2i}+\beta {u}_{2i+1}}.$

From (2.2) and (2.3), we note that $f$ is increasing in ${u}_{2r+1}$ and decreasing in ${u}_{2r}$ for all $r=0,1,\cdots ,k$ . Suppose that $\left(\lambda ,\mu \right)$ is a solution of the system

$\lambda =f\left(\mu ,\lambda ,\mu ,\lambda ,\cdots ,\lambda \right)$

$\mu =f\left(\lambda ,\mu ,\lambda ,\mu ,\cdots ,\mu \right).$

Then, we find

$\lambda =a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{\lambda }{\alpha \mu +\beta \lambda },$

and

$\mu =a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{\mu }{\alpha \lambda +\beta \mu }.$

Hence, we get

$\lambda =a+\frac{B\lambda }{\alpha \mu +\beta \lambda },$ (3.1)

and

$\mu =a+\frac{B\mu }{\alpha \lambda +\beta \mu }.$ (3.2)

By (3.1) and (3.2), we obtain

$\beta \left({\lambda }^{2}-{\mu }^{2}\right)-\left(a\beta +B-a\alpha \right)\left(\lambda -\mu \right)=0.$

Thus,

$\left(\lambda -\mu \right)\left(\beta \left(\lambda +\mu \right)-\left(a\beta +B-a\alpha \right)\right)=0.$

Since $B=a\left(\alpha -\beta \right)$ , we have that $\mu =\lambda$ . Hence, the proof of Theorem 3.1 is complete. □

4. Periodic Solutions

In this section, we enumerate some basic facts concerning the existence of two period solutions.

Theorem 4.1. Equation (1.1) has prime period-two solutions if

$\left(a\beta +B-a\alpha \right)\left(\alpha -\beta \right)>4a\alpha \beta .$ (4.1)

Proof. Assume that Equation (1.1) has a prime period-two solution

$\cdots ,\rho ,\sigma ,\rho ,\sigma ,\rho ,\sigma ,\cdots$

We shall prove that condition (4.1) holds. From Equation (1.1), we see that

${\psi }_{\eta }={\psi }_{\eta -2}=\cdots ={\psi }_{\eta -2k}=\sigma ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\psi }_{\eta +1}={\psi }_{\eta -1}=\cdots ={\psi }_{\eta -\left(2k+1\right)}=\rho ,$

and hence,

$\rho =a+\frac{B\rho }{\alpha \sigma +\beta \rho },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sigma =a+\frac{B\sigma }{\alpha \rho +\beta \sigma }.$

Thus, we get

$\alpha \rho \sigma +\beta {\rho }^{2}=a\alpha \sigma +a\beta \rho +B\rho ,$ (4.2)

and

$\alpha \rho \sigma +\beta {\rho }^{2}=a\alpha \sigma +a\beta \sigma +B\sigma .$ (4.3)

From (4.3) and (4.2), we have

$\beta \left({\rho }^{2}-{\sigma }^{2}\right)=a\alpha \left(\sigma -\rho \right)+a\beta \left(\rho -\sigma \right)+B\left(\rho -\sigma \right).$

Dividing $\left(\rho -\sigma \right)$ , then we find

$\left(\rho +\sigma \right)=\frac{a\beta +B-a\alpha }{\beta }$ (4.4)

By combining (4.2) and (4.3), we obtain

$2\alpha \rho \sigma +\beta \left({\rho }^{2}+{\sigma }^{2}\right)=\left(a\alpha +a\beta +B\right)\left(\rho +\sigma \right).$

Since ${\rho }^{2}+{\sigma }^{2}={\left(\rho +\sigma \right)}^{2}-2\rho \sigma$ , we get

$\rho \sigma =\frac{a\alpha \left(a\beta +B-a\alpha \right)}{\beta \left(\alpha -\beta \right)}.$ (4.5)

Now, evident is that (4.4) and (4.5) that $\rho$ and $\sigma$ are both two positive distinct roots of the quadratic equation

${u}^{2}+\left(\rho +\sigma \right)u+\rho \sigma =0.$ (4.6)

Hence, we obtain

$\frac{\left(a\beta +B-a\alpha \right)}{\beta }>\frac{4a\alpha }{\alpha -\beta },$

which has the same extent as

$\left(a\beta +B-a\alpha \right)\left(\alpha -\beta \right)>4a\alpha \beta .$

Hence, the proof is complete. □

The next numerical example is mimicry to enhance our results.

Example 4.1. By Theorem 4.1, Equation (1.1) with $\alpha =0.02$ , $\beta =0.01$ , $a=500$ , ${b}_{0}=2$ , ${b}_{1}=200$ and ${b}_{2}=20$ , has prime period two solution (see Figure 2)

5. Boundedness

Theorem 5.1. Every solution of Equation (1.1) is bounded and persists.

Proof. Let ${\left\{{\psi }_{\eta }\right\}}_{\eta =-k}^{\infty }$ be a Solution (1.1), we can conclude from (1.1) that

${\psi }_{\eta +1}=a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{{\psi }_{\eta -\left(2i+1\right)}}{\alpha {\psi }_{\eta -2i}+\beta {\psi }_{\eta -\left(2i+1\right)}}>a.$

Then

${\psi }_{\eta }>a\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{all}\text{\hspace{0.17em}}\eta >-\left(2k+1\right).$

Also, from Equation (1.1), we see that

$\begin{array}{c}{\psi }_{\eta +1}=a+\underset{i=0}{\overset{k}{\sum }}\text{ }\text{ }{b}_{i}\frac{{\psi }_{\eta -\left(2i+1\right)}}{\alpha {\psi }_{\eta -2i}+\beta {\psi }_{\eta -\left(2i+1\right)}}\\

then,

$a<{\psi }_{\eta }\le a+\frac{B}{\beta }\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{all}\text{\hspace{0.17em}}\eta >-\left(2k+1\right).$

Thus, the solution is bounded and persists and the proof is complete. □

Conclusion 1. In this paper, we study a asymptotic behavior of solutions of a general class of difference Equation (1.1). Our results extend and generalize to the earlier ones. Moreover, we obtain the next results:

- The $+ve$ equilibrium point $\stackrel{^}{\psi }$ of Equation (1.1) is local stable if $\left(\alpha -\beta \right)B . Also, if $B=a\left(\alpha -\beta \right)$ , then $\stackrel{^}{\psi }$ is global attractor.

Figure 2. Prime period two solution of Equation (1.1).

- Equation (1.1) has a prime period-two solutions if $\left(a\beta +B-a\alpha \right)\left(\alpha -\beta \right)>4a\alpha \beta$ .

- Every solution of (1.1) is bounded and persists.

Acknowledgements

The author is very grateful to the reviewers for their valuable suggestions and useful comments on this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

Moaaz, O. , Elabbasy, E. and Alsaeed, S. (2017) Global Analysis of Solutions of a New Class of Rational Difference Equation. American Journal of Computational Mathematics, 7, 495-503. doi: 10.4236/ajcm.2017.74036. 