Qualitative behavior of a second-order fuzzy difference equation

Abstract

In this paper, we study the qualitative behavior of the positive solutions of a second-order rational fuzzy difference equations with initial conditions and parameters are positive fuzzy numbers. More precisely, we investigate existence of positive solutions, boundedness, persistence and stability analysis of a second-order fuzzy rational difference equation. Some numerical examples are given to verify our theoretical results.

Keywords

Fuzzy difference equations existence of solutions boundedness character local and global behavior

1 Introduction

The theory of discrete dynamical systems and difference equations developed greatly during the last twenty-five years of the twentieth century. Applications of discrete dynamical systems and difference equations have appeared recently in many areas. The theory of difference equations occupies a central position in applied analysis. There is no doubt that the theory of difference equations continue to play an important role in mathematics as a whole. Nonlinear difference equations of order greater than one are of paramount importance in applications. Such equations also appear naturally as discrete analogues and as numerical solutions of differential and delay differential equations which model various diverse phenomena in biology, ecology, physiology physics, engineering, economics, probability theory, genetics, psychology and resource management. It is very interesting to investigate the behavior of solutions of higher-order rational difference equations and to discuss the local asymptotic stability of its equilibrium points. Rational difference equations have been studied by several authors. Especially there has been a great interest in the study of the attractivity of the solutions of such equations. Rational difference equations are special type of nonlinear difference equations. When we divide two difference equations (may be linear or nonlinear), we obtain a rational difference equation (which is always nonlinear). The solutions of rational difference equations can be found rarely in compact or closed forms. For basic theory and applications of nonlinear difference equations we refer the interested reader to [1, 2]. For applications of nonlinear difference equations in population dynamics we refer to [3–7 , 18–21].

In case of fuzzy difference equations, parameters, initial conditions and solutions of these equations are fuzzy numbers, while in case of difference equations all these are usually real numbers. Hence study of fuzzy difference equations is more interesting as well as complicated. Using difference equation approach we can investigate dynamical aspects of advanced models in epidemiology, finance and other branches of science. In difference equations assuming the transmission coefficients and other parameter as positive fuzzy numbers, we can explore several dynamics of humans health related issues. These issues include cardiology, hydrology, health services and vital status ascertainment in clinical trials. The readers are referred to see the monographs [22 , 29]. For basic theory of fuzzy difference equations see [8]. In [9], embedding problem of non-compact fuzzy number space is given. It is very interesting to investigate the dynamics of nonlinear fuzzy difference equations [10 –15]. Bajo and Liz [16] investigated the global behavior of difference equation: $x_{n + 1} = \frac{x_{n - 1}}{a + {bx}_{n - 1} x_{n}},$ for all values of real parameters a, b.

Our aim in this paper is to investigate the qualitative behavior of following second-order fuzzy rational difference equation $x_{n + 1} = \frac{x_{n - 1}}{A + {Bx}_{n - 1} x_{n}},$ (1) where A, B and initial conditions x_-1, x₀ are positive fuzzy numbers.

2 Preliminaries

Definition 1. Let X be a normed space. We will denote by $𝔼 (X)$ the family of fuzzy sets u : X → [0, 1] which has the following properties:

u is normal, i . e ., there exist an x in X such that u (x) = 1,

u is fuzzy convex function, i . e ., min { u (x₁) , u (x₂) } ≤ u (tx₁ + (1 - t) x₂) for all t ∈ [0, 1], and x₁, x₂ ∈ X,

u is upper semi-continuous,

The support of u, $supp u = \bar{⋃_{α \in [0, 1]} {[u]}_{α}} = \bar{{x : u (x) > α}}$ is compact.

An element $u \in 𝔼 (X)$ is called a fuzzy vector. If $X = ℝ$ , then an element $u \in 𝔼 (ℝ)$ is called fuzzy number.

Remark 1. In particular, if { [a_α, b_α] } _α∈(0,1] is a family of compact intervals in $ℝ$ such that

[a_β, b_β] ⊂ [a_α, b_α] for all 0 < α ≤ β

$[lim_{m \to \infty} a_{α_{m}}, lim_{m \to \infty} b_{α_{m}}] = [a_{α}, b_{α}]$

whenever {α_m } ⊂ (0, 1] in an increasing sequence such that

lim_{m \to \infty} α_{m} = α

for m ≥ 1 then, the family { [a_α, b_α] } _α∈(0,1] represent the α-level sets for the fuzzy number

u : ℝ \to [0, 1]

, defined by

u (x) = \sup {α \in [0, 1] : x \in [a_{α}, b_{α}]} .

Moreover, α-cuts of u are given by ${[u]}_{α} = {x \in ℝ : u (x) \geq α}$ for α ∈ [0, 1]. We say a fuzzy number u is positive if supp u ⊂ (0, ∞). Let u, v be positive fuzzy numbers with [u] _α = [u_l,α, u_r,α] and [v] _α = [v_l,α, v_r,α] for α ∈ (0, 1], then we consider the following metric: $D (u, v) = sup_{α \in (0, 1]} \max {| u_{l, α} - v_{l, α} |, | u_{r, α} - v_{r, α} |} .$

A sequence (x_n) of fuzzy number is said to be bounded and persists if there exist positive real numbers β, γ such that supp x_n ⊂ [β, γ] for n = 1, 2, ⋯. Moreover, if (x_n) be a sequence of positive fuzzy numbers and x be any positive fuzzy number such that [x_n] _α = [L_l,α, R_r,α] and [x] _α = [L_α, R_α] for α ∈ (0, 1] and n = 0, 1, 2, ⋯ then (x_n) converges to x with respect to metric D, if $lim_{n \to \infty} D (x_{n}, x) = 0$ .

Definition 2. Let (x_n) be a positive solution of the Equation (1) and $\bar{x}$ be its equilibrium point, i . e ., $\bar{x} = \frac{\bar{x}}{A + B {\bar{x}}^{2}}$ , then

$\bar{x}$ is stable, if for given ɛ > 0 there exist a δ = δ (ɛ) > 0 such that if D (x_-i, x) ≤ δ for i∈ { - 1, 0 }, then $D (x_{n}, \bar{x}) \leq ɛ$ for all n > 0.

$\bar{x}$ is said to be unstable if it is not stable.

$\bar{x}$ is locally asymptotically stable, if it is stable and $lim_{n \to \infty} D (x_{n}, \bar{x}) = 0$ .

$\bar{x}$ is said to be an attractor if $lim_{n \to \infty} D (x_{n}, \bar{x}) = 0$ .

$\bar{x}$ is globally asymptotically stable, if it is an attractor and locally asymptotically stable.

3 Main results

This section is devoted to the main results, which are associated with the proposed problem. Arguing as in [15], we have the following result.

Theorem 1. Let A, B are positive fuzzy numbers, then Equation (1) has unique positive solution {x_n} with initial conditions x_-1, x₀ are positive fuzzy numbers.

Proof. Let {x_n} be a positive solution of Equation (1). The α-cuts of positive fuzzy numbers x_n, A, B are given by $\begin{matrix} {[x_{n}]}_{α} = [L_{l, α}, R_{r, α}], {[A]}_{α} = [A_{l, α}, A_{r, α}], \\ {[B]}_{α} = [B_{l, α}, B_{r, α}] . \end{matrix}$

Then from Equation (1), one has $\begin{matrix} {[x_{n + 1}]}_{α} & = & [L_{l + 1, α}, R_{r + 1, α}] \\ = & {[\frac{x_{n - 1}}{A + {Bx}_{n - 1} x_{n}}]}_{α} = \frac{[x_{n - 1}]_{α}}{[A]_{α} + [B]_{α} [x_{n - 1}]_{α} [x_{n}]_{α}} \\ = & \frac{[L_{n - 1, α}, R_{n - 1, α}]}{[A_{l, α}, A_{r, α}] + [B_{l, α}, B_{r, α}] [L_{n - 1, α}, R_{n - 1, α}] [L_{n, α}, R_{n, α}]} \\ = & [\frac{L_{n - 1, α}}{A_{r, α} + B_{r, α} R_{n - 1, α} R_{n, α}}, \frac{R_{n - 1, α}}{A_{l, α} + B_{l, α} L_{n - 1, α} L_{n, α}}] . \end{matrix}$

Hence, we obtain the following system

$\begin{matrix} L_{n + 1, α} & = & \frac{L_{n - 1, α}}{A_{r, α} + B_{r, α} R_{n - 1, α} R_{n, α}} \\ R_{n + 1, α} & = & \frac{R_{n - 1, α}}{A_{l, α} + B_{l, α} L_{n - 1, α} L_{n, α}} \end{matrix}$ (2) for n = 0, 1, 2, ⋯ , and α ∈ (0, 1]. Let (L_n,α, R_n,α) be a solution of the System (2) with initial conditions L_-1,α, L_0,α, R_-1,α, R_0,α for every α ∈ (0, 1]. Then, it is sufficient to prove that (x_n) is a solution of original System (1) with initial conditions x_-1, x₀ and [x_n] _α = [L_n,α, R_n,α] for every n = 0, 1, 2, ⋯ and α ∈ (0, 1]. Moreover, for any fuzzy number $u \in 𝔼 (ℝ)$ with ${[u]}_{α} = [\underline{u} (α), \bar{u} (α)]$ , one has $\underline{u} (α)$ is non-decreasing and left continuous, while $\bar{u} (α)$ is non-increasing and left continuous for 0 < α ≤ 1. As A, B, x_-1, x₀ are positive fuzzy numbers, so for α₁, α₂ ∈ (0, 1] with α₁ ≤ α₂, one has $\begin{matrix} 0 < A_{l, α_{1}} \leq A_{l, α_{2}} \leq A_{r, α_{2}} \leq A_{r, α_{1}}, \\ 0 < B_{l, α_{1}} \leq B_{l, α_{2}} \leq B_{r, α_{2}} \leq B_{r, α_{1}}, \\ 0 < L_{- 1, α_{1}} \leq L_{- 1, α_{2}} \leq R_{- 1, α_{2}} \leq R_{- 1, α_{1}}, \\ 0 < L_{0, α_{1}} \leq L_{0, α_{2}} \leq R_{0, α_{2}} \leq R_{0, α_{1}} . \end{matrix}$

Now, it is easy to show by induction on n that 0 < L_n,α₁ ≤ L_n,α₂ ≤ R_n,α₂ ≤ R_n,α₁, n = 0, 1, 2, ⋯.

Similarly, by induction on n, one can show that L_n,α, R_n,α are left continuous for all n = 0, 1, 2, ⋯, and 0 < α ≤ 1. Furthermore, it can be inductively proved that $\bar{⋃_{α \in [0, 1]} [L_{n, α}, R_{n, α}]} \subset (0, \infty)$ , $i . e ., \bar{\cup_{α \in [0, 1]} [L_{n, α}, R_{n, α}]}$ is compact. Hence [L_n,α, R_n,α] determines a sequence (x_n) of positive fuzzy numbers such that [x_n] _α = [L_n,α, R_n,α] for every n = 0, 1, 2, ⋯, and α ∈ (0, 1] . ■

In order to study the further dynamics of the fuzzy difference Equation (1) we will use the results concerning the behavior of the solutions of the corresponding system of two parametric ordinary difference equations given by:

$u_{n + 1} = \frac{u_{n - 1}}{a_{1} + b_{1} v_{n - 1} v_{n}}, v_{n + 1} = \frac{v_{n - 1}}{a_{2} + b_{2} u_{n - 1} u_{n}},$ (3) where the parameters a₁, a₂, b₁, b₂ as well as the initial conditions u_-1, u₀, v_-1, v₀ are positive real numbers.

Let a₁ < 1, a₂ < 1, then O = (0, 0) and $E = (\sqrt{\frac{1 - a_{2}}{b_{2}}}, \sqrt{\frac{1 - a_{1}}{b_{1}}})$ be equilibrium points of the System (3).

Theorem 2. Let ${(u_{n}, v_{n})}_{n = - 1}^{\infty}$ be any solution of the System (3), then the following statements are true:

For every m ≥ 0 the following results hold: $\begin{matrix} 0 \leq u_{n} \leq {(\frac{1}{a_{1}})}^{m + 1} u_{- 1}, if n = 2 m + 1, \\ 0 \leq u_{n} \leq {(\frac{1}{a_{1}})}^{m + 1} u_{0}, if n = 2 m + 2, \\ 0 \leq v_{n} \leq {(\frac{1}{a_{2}})}^{m + 1} v_{- 1}, if n = 2 m + 1, \\ 0 \leq v_{n} \leq {(\frac{1}{a_{2}})}^{m + 1} v_{0}, if n = 2 m + 2 . \end{matrix}$

If 1 < a₁ and 1 < a₂, then every solution ${(u_{n}, v_{n})}_{n = - 1}^{\infty}$ of the System (3) is bounded.

For the equilibrium point O = (0, 0) of the System (3) the following results hold:

Let 1 < a₁ and 1 < a₂, then equilibrium point O = (0, 0) of the System (3) is locally asymptotically stable.

If 1 > a₁ or 1 > a₂, then equilibrium point O = (0, 0) of the System (3) is unstable.

If 1 > a₁ and 1 > a₂, then equilibrium point $E = (\sqrt{\frac{1 - a_{2}}{b_{2}}}, \sqrt{\frac{1 - a_{1}}{b_{1}}})$ of the System (3) is unstable.

Let 1 > a₁ and 1 > a₂. Then for k = -1, 0 following statements are true:

If $(u_{k}, v_{k}) \in (0, {(\frac{1 - a_{2}}{b_{1}})}^{\frac{1}{2}}) \times ({(\frac{1 - a_{1}}{b_{1}})}^{\frac{1}{2}}, \infty)$ then $(u_{n}, v_{n}) \in (0, {(\frac{1 - a_{2}}{b_{1}})}^{\frac{1}{2}}) \times ({(\frac{1 - a_{1}}{b_{1}})}^{\frac{1}{2}}, \infty)$ .

If $(u_{k}, v_{k}) \in ({(\frac{1 - a_{2}}{b_{2}})}^{\frac{1}{2}}, \infty) \times (0, {(\frac{1 - a_{1}}{b_{1}})}^{\frac{1}{2}})$ , then $(u_{n}, v_{n}) \in ({(\frac{1 - a_{2}}{b_{2}})}^{\frac{1}{2}}, \infty) \times (0, {(\frac{1 - a_{1}}{b_{1}})}^{\frac{1}{2}})$ .

Let 1 < a₁ and 1 < a₂, then the equilibrium point O = (0, 0) of the System (3) is globally asymptotically stable.

Proof. (1) The proof follows from induction.

(2) Let α = max {u_-1, u₀} and β = max{ v_-1, v₀ }, then it is easy to see from (1) that 0 ≤ u_n < α, 0 ≤ v_n < β for all n = 0, 1, 2, ⋯.

Hence, every solution ${(u_{n}, v_{n})}_{n = - 1}^{\infty}$ of the System (3) is bounded.

(4) To construct corresponding linearized form of System (3), we consider the following transformations: $(u_{n}, u_{n - 1}, v_{n}, v_{n - 1}) \to (f, f_{1}, g, g_{1})$ where f = u_n+1, g = v_n+1, f₁ = u_n, g₁ = v_n. The linearized system of (3) about the equilibrium point (0, 0) is given by: $X_{n + 1} = F_{J} (0, 0) X_{n},$ Where $(\begin{matrix} u_{n} \\ u_{n - 1} \\ \begin{matrix} v_{n} \\ v_{n - 1} \end{matrix} \end{matrix})$ and $F_{J} (0, 0) = (\begin{matrix} 0 & \frac{1}{a_{1}} & \begin{matrix} 0 & 0 \end{matrix} \\ 1 & 0 & \begin{matrix} 0 & 0 \end{matrix} \\ \begin{matrix} 0 \\ 0 \end{matrix} & \begin{matrix} 0 \\ 0 \end{matrix} & \begin{matrix} \begin{matrix} 0 \\ 1 \end{matrix} & \begin{matrix} \frac{1}{a_{2}} \\ 0 \end{matrix} \end{matrix} \end{matrix})$ .

The characteristic polynomial of F_J (0, 0) is given by $P (λ) = λ^{4} - (\frac{1}{a_{1}} + \frac{1}{a_{2}}) λ^{2} - \frac{1}{a_{1} a_{2}} .$ (4)

Since all eigenvalues of Jacobian matrix F_J (0, 0) about (0, 0) lie in open unit disk |λ| < 1, if 1 < a₁ and 1 < a₂. Hence, the equilibrium point (0, 0) is locally asymptotically stable.

Similarly, all eigen values of Jacobian matrix F_J (0, 0) about (0, 0) lie in open unit disk |λ| < 1, if 1 > a₁ and 1 > a₂. Hence, the equilibrium point (0, 0) is unstable in this case.

(4) The linearized system of (2) about the equilibrium point $E = (\sqrt{\frac{1 - a_{2}}{b_{2}}}, \sqrt{\frac{1 - a_{1}}{b_{1}}})$ is given by:

X_n+1 = F_J (E) X_n, where $X_{n} = (\begin{matrix} u_{n} \\ u_{n - 1} \\ \begin{matrix} v_{n} \\ v_{n - 1} \end{matrix} \end{matrix})$ and $F_{J} (E) (\begin{matrix} 0 & 1 & - \begin{matrix} \frac{\sqrt{1 - a_{1}} \sqrt{1 - a_{2}} \sqrt{b_{1}}}{\sqrt{b_{2}}} & - \frac{\sqrt{1 - a_{1}} \sqrt{1 - a_{2}} \sqrt{b_{1}}}{\sqrt{b_{2}}} \end{matrix} \\ 1 & 0 & \begin{matrix} 0 & 0 \end{matrix} \\ \begin{matrix} - \frac{\sqrt{1 - a_{1}} \sqrt{1 - a_{2}} \sqrt{b_{2}}}{\sqrt{b_{1}}} \\ 0 \end{matrix} & \begin{matrix} - \frac{\sqrt{1 - a_{1}} \sqrt{1 - a_{2}} \sqrt{b_{2}}}{\sqrt{b_{1}}} \\ 0 \end{matrix} & \begin{matrix} \begin{matrix} 0 \\ 1 \end{matrix} & \begin{matrix} 0 \\ 0 \end{matrix} \end{matrix} \end{matrix})$ .

The characteristic polynomial of F_J (E) is given by

$\begin{matrix} P (λ) & = & λ^{4} + (a_{1} + a_{2} - a_{1} a_{2} - 3) λ^{2} \\ + 2 (a_{1} + a_{2} - a_{1} a_{2} - 1) λ \\ + a_{1} + a_{2} - a_{1} a_{2} . \end{matrix}$ (5)

The roots of characteristic polynomial P (λ) are given by $- 1, 1 \pm \sqrt{a_{1} a_{2} - a_{1} - a_{2} + 1} .$

It is sufficient to prove that any one of these roots has absolute value greater than one. For this consider $\begin{matrix} | 1 + \sqrt{a_{1} a_{2} - a_{1} - a_{2} + 1} | \\ = | 1 + \sqrt{(1 - a_{1}) (1 - a_{2})} | > 1 . \end{matrix}$

(5) The proof of 5(i) and 5(ii) follow from induction.

(6) For 1 < a₁ and 1 < a₂, from (i) of (3) of Theorem 2 (0, 0) is locally asymptotically stable. Moreover, from (2) the solution { (u_n, v_n)} of the System (3) is bounded if 1 < a₁ and 1 < a₂. Now, it is sufficient to prove that { (u_n, v_n)} is decreasing. From the System (3) one has $\begin{matrix} u_{n + 1} & = & \frac{u_{n - 1}}{a_{1} + b_{1} v_{n - 1} v_{n}} \\ \geq & \frac{u_{n - 1}}{a_{1}} < u_{n - 1} . \end{matrix}$

This implies that u_2n+1 < u_2n-1 and u_2n+3 < u_2n+1. Hence the subsequences {u_2n+1} , { u_2n+2 } are decreasing, i . e ., the sequence {u_n} is decreasing. Similarly, one has $\begin{matrix} v_{n + 1} & = & \frac{v_{n - 1}}{a_{2} + b_{2} u_{n} u_{n - 1}} \\ \leq & \frac{v_{n - 1}}{a_{2}} < v_{n - 1} . \end{matrix}$

This implies that v_2n+1 < v_2n-1 and v_2n+3 < v_2n+1. Hence the subsequences {v_2n+1} , { v_2n+2 } are decreasing, i . e ., the sequence {v_n} is decreasing. Hence $\lim_{n \to \infty} x_{n} = \lim_{n \to \infty} y_{n} = 0 ■$

The following results give the rate of convergence of solutions of a system of difference equations $X_{n + 1} = (A + B (n)) X_{n},$ (6) where X_n is an m-dimensional vector, A ∈ C^m×m is a constant matrix, and B : Z⁺ → C^m×m is a matrix function satisfying $∥ B (n) ∥ \to 0$ (7) as n→ ∞, where ∥· ∥ denotes any matrix norm which is associated with the vector norm $∥ B (n) ∥ = \sqrt{x^{2} + y^{2}} .$

Proposition 1. (Perron’s theorem) [17] Suppose that condition (8) holds. If X_n is a solution of (7), then either X_n = 0 for all large n or $ρ = lim_{n \to \infty} {(∥ X_{n} ∥)}^{\frac{1}{n}}$ (8) exist and is equal to the modulus of one of the eigen values of matrix A.

Proposition 2. [17] Suppose that condition (8) holds. If X_n is a solution of (7), then either X_n = 0 for all large n or $ρ = lim_{n \to \infty} \frac{∥ X_{n + 1} ∥}{∥ X_{n} ∥}$ (9) exist and is equal to the modulus of one of the eigen values of matrix A.

Let {(u_n, v_n)} be any solution of the System (3) such that $\lim_{n \to \infty} u_{n} = \bar{u}$ , $\lim_{n \to \infty} v_{n} = \bar{v}$ , where $(\bar{u}, \bar{v}) = (0, 0) .$ To find the error terms, from the System (3) we have $\begin{matrix} u_{n + 1} - \bar{u} & = & \frac{u_{n - 1}}{a_{1} + b_{1} v_{n - 1} v_{n}} - \frac{\bar{u}}{a_{1} + b_{1} {\bar{v}}^{2}} \\ = & \frac{(a_{1} + b_{1} v_{n} \bar{v}) (u_{n - 1} - \bar{u}) - b_{1} u_{n - 1} \bar{v} (v_{n} - \bar{v}) - b_{1} v_{n} \bar{u} (v_{n - 1} - \bar{v})}{(a_{1} + b_{1} v_{n - 1} v_{n}) (a_{1} + b_{1} {\bar{v}}^{2})} \\ = & a_{n} (u_{n} - \bar{u}) + b_{n} (u_{n - 1} - \bar{u}) \\ + c_{n} (v_{n} - \bar{v}) + d_{n} (v_{n - 1} - \bar{v}) \end{matrix}$ and $\begin{matrix} v_{n + 1} - \bar{v} & = & \frac{v_{n - 1}}{a_{2} + b_{2} u_{n} u_{n - 1}} - \frac{\bar{v}}{a_{2} + b_{2} {\bar{u}}^{2}} \\ = & \frac{(a_{2} + b_{2} u_{n} \bar{u}) (v_{n - 1} - \bar{v}) - b_{2} v_{n - 1} \bar{u} (u_{n} - \bar{u}) - b_{2} u_{n} \bar{v} (u_{n - 1} - \bar{u})}{(a_{2} + b_{2} u_{n} u_{n - 1}) (a_{2} + b_{2} {\bar{u}}^{2})} \\ = & a_{n}^{1} (u_{n} - \bar{u}) + b_{n}^{1} (u_{n - 1} - \bar{u}) + c_{n}^{1} (v_{n} - v) \\ + d_{n}^{1} (v_{n - 1} - \bar{v}) \end{matrix}$

Let $e_{n}^{1} = u_{n} - \bar{u}$ and $e_{n}^{2} = v_{n} - \bar{v}$ , one has $e_{n + 1}^{1} = a_{n} e_{n}^{1} + b_{n} e_{n - 1}^{1} + c_{n} e_{n}^{2} + d_{n} e_{n - 1}^{2},$ and $e_{n + 1}^{2} = a_{n}^{1} e_{n}^{1} + b_{n}^{1} e_{n - 1}^{1} + c_{n}^{1} e_{n}^{2} + d_{n}^{1} e_{n - 1}^{2},$ where $\begin{matrix} a_{n} = 0, b_{n} & = & \frac{a_{1} + b_{1} v_{n} \bar{v}}{(a_{1} + b_{1} v_{n - 1} v_{n}) (a_{1} + b_{1} {\bar{v}}^{2})}, \\ c_{n} & = & - \frac{b_{1} u_{n - 1} \bar{v}}{(a_{1} + b_{1} v_{n - 1} v_{n}) (a_{1} + b_{1} {\bar{v}}^{2})} \\ dn & = & - \frac{b_{1} v_{n} \bar{u}}{(a_{1} + b_{1} v_{n - 1} v_{n}) (a_{1} + b_{1} {\bar{v}}^{2})} \\ a_{n}^{1} & = & - \frac{b_{2} u_{n} \bar{v}}{(a_{2} + b_{2} u_{n} u_{n - 1}) (a_{2} + b_{2} {\bar{u}}^{2})} \\ b_{n}^{1} & = & - \frac{b_{2} u_{n} \bar{v}}{(a_{2} + b_{2} u_{n} u_{n - 1}) (a_{2} + b_{2} {\bar{u}}^{2})} \\ c_{n}^{1} & = & 0 \\ d_{n}^{1} & = & \frac{a_{2} + b_{2} u_{n} \bar{u}}{(a_{2} + b_{2} u_{n} u_{n - 1}) (a_{2} + b_{2} {\bar{u}}^{2})} . \end{matrix}$

Moreover, $\begin{matrix} \lim_{n \to \infty} a_{n} = \lim_{n \to \infty} c_{n} = \lim_{n \to \infty} d_{n} = 0, \lim_{n \to \infty} b_{n} = \frac{1}{a_{1}}, \\ \lim_{n \to \infty} a_{n}^{1} = \lim_{n \to \infty} b_{n}^{1} = \lim_{n \to \infty} c_{n}^{1} = 0 \lim_{n \to \infty} d_{n}^{1} = \frac{1}{a_{2}} . \end{matrix}$

Now the limiting system of error terms can be written as $(\begin{matrix} e_{n + 1}^{1} \\ e_{n}^{1} \\ \begin{matrix} e_{n + 1}^{2} \\ e_{n}^{2} \end{matrix} \end{matrix}) = (\begin{matrix} 0 & \frac{1}{a_{1}} & \begin{matrix} 0 & 0 \end{matrix} \\ 1 & 0 & \begin{matrix} 0 & 0 \end{matrix} \\ \begin{matrix} 0 \\ 0 \end{matrix} & \begin{matrix} 0 \\ 0 \end{matrix} & \begin{matrix} \begin{matrix} 0 \\ 1 \end{matrix} & \begin{matrix} \frac{1}{a_{2}} \\ 0 \end{matrix} \end{matrix} \end{matrix}) (\begin{matrix} e_{n}^{1} \\ e_{n - 1}^{1} \\ \begin{matrix} e_{n}^{2} \\ e_{n - 1}^{2} \end{matrix} \end{matrix}),$ which is similar to linearized system of (3) about the equilibrium point (0, 0).

Using Proposition 1, one has following result.

Theorem 3. Assume that { (u_n, v_n)} be a positive solution of the System (2) such that $\lim_{n \to \infty} u_{n} = 0$ and $\lim_{n \to \infty} v_{n} = 0$ . Then, the error term $E_{n} = (\begin{matrix} e_{n}^{1} \\ e_{n - 1}^{1} \\ \begin{matrix} e_{n}^{2} \\ e_{n - 1}^{2} \end{matrix} \end{matrix})$ of every solution of (3) satisfies both of the following asymptotic relations, $\begin{matrix} \lim_{n \to \infty} {(∥ E_{n} ∥)}^{\frac{1}{n}} = | λ F_{J} (0, 0) | \\ \lim_{n \to \infty} \frac{∥ E_{n + 1} ∥}{∥ E_{n} ∥} = | λ F_{J} (0, 0) |, \end{matrix}$ where λF_J (0, 0) are the characteristic roots of Jacobian matrix F_J (0, 0).

Theorem 4. Assume that (x_n) be a positive solution of (1) and 1 < A_l,α, then (x_n) is bounded and persists. Moreover, (x_n) converges to 0.

Proof. Consider the following system $s_{n + 1} = \frac{s_{n - 1}}{γ_{A} + γ_{B} t_{n - 1} t_{n}}, t_{n + 1} = \frac{t_{n - 1}}{β_{A} + β_{B} s_{n - 1} s_{n}}$ (10)

Where $[A_{l, α}, A_{r, α}] \subset \bar{\cup_{α \in [0, 1]} [A_{l, α}, A_{r, α}]} \subset [β_{A}, γ_{A}],$ and $[B_{l, α}, B_{r, α}] \subset \bar{\cup_{α \in [0, 1]} [B_{l, α}, B_{r, α}]} \subset [β_{B}, γ_{B}] .$

Let (s_n, t_n) be a solution of the System (11) with initial conditions s_-1 = β_-1, s₀ = β₀, t_-1 = γ_-1, t₀ = γ₀, where β_i, γ_i for i∈ { - 1, 0 } are given as $\begin{matrix} [L_{- 1, α}, R_{- 1, α}] & \subset & \bar{\cup_{α \in [0, 1]} [L_{- 1, α}, R_{- 1, α}]} \\ \subset & [β_{- 1}, γ_{- 1}], \end{matrix}$ and $[L_{0, α}, R_{0, α}] \subset \bar{\cup_{α \in [0, 1]} [L_{0, α}, R_{0, α}]} \subset [β_{0}, γ_{0}] .$

Then, it follows that $\begin{matrix} s_{1} & = & \frac{s_{- 1}}{γ_{A} + γ_{B} t_{- 1} t_{0}} \leq L_{1, α}, \\ t_{1} & = & \frac{t_{- 1}}{β_{A} + β_{B} s_{- 1} s_{0}} \geq R_{1, α} . \end{matrix}$

Hence by induction one has s_n ≤ L_n,α, t_n ≥ R_n,α for all n = 1, 2, ⋯. Assume that 1 < A_l,α, then it follows that 1 < γ_A, 1 < β_A. Hence the solution (s_n, t_n) of the System (11) is bounded and persist, which is the solution (x_n) of the Equation (1). Next, from (6) of the Theorem 2, it is easy to see that $lim_{n \to \infty} L_{n, α} = lim_{n \to \infty} R_{n, α} = 0$ . Then it follows that $lim_{n \to \infty} D (x_{n}, 0) = \lim_{n \to \infty} \sup_{α \in (0, 1]} {\max {| L_{n, α} - 0 |, | R_{r, α} - 0 |}} = 0 .$ ■

Theorem 5. Assume that A, B are positive real numbers and A < 1, then unique positive equilibrium point $\bar{x} = \frac{\bar{x}}{A + B {\bar{x}}^{2}}$ of the Equation (1) is unstable.

Proof. Let (x_n) be a positive solution of the Equation (1) and A, B are positive real numbers with A < 1. Then, from (5)(i) and (5)(ii) of Theorem 2, either $lim_{n \to \infty} R_{n, α} = 0$ and $lim_{n \to \infty} L_{n, α} = \infty$ , or $lim_{n \to \infty} R_{n, α} = \infty$ and $lim_{n \to \infty} L_{n, α} = 0$ . Hence, there are no positive numbers β, γ such that ∪_α∈[0,1] [L_n,α, R_n,α] ⊂ [β, γ]. This completes the proof. ■

4 Examples

Example 1. Consider the System (3) with initial conditions u₀ = 0.05, v₀ = 0.1, u_-1 = 0.07, v_-1 = 0.06. Moreover, choosing the parameters a₁ = a₂ = 1.006, b₁ = b₂ = 2.8.

Then, the System (3) can be written as

$\begin{matrix} u_{n + 1} & = & \frac{u_{n - 1}}{1.006 + 2.8 v_{n - 1} v_{n}}, \\ v_{n + 1} & = & \frac{v_{n - 1}}{1.006 + 2.8 u_{n - 1} u_{n}}, \end{matrix}$ (11) with initial condition u₀ = 0.05, v₀ = 0.1, u_-1 = 0.07, v_-1 = 0.06. The plot of solutions for the System (12) is shown in Fig. 1 and its global attractor is shown in Fig. 2.

Fig.1

Plot of solutions for the System (12).

Fig.2

An attractor for the System (12).

Example 2. Consider System (3) with initial conditions u₀ = 0.2, v₀ = 0.4, u_-1 = 0.5, v_-1 = 0.05. Moreover, choose the parameters a₁ = 1.05, a₂ = 1.03, b₁ = 10, b₂ = 2.

Then, the System (3) can be written as

$\begin{matrix} u_{n + 1} & = & \frac{u_{n - 1}}{1.05 + 10 v_{n - 1} v_{n}}, \\ v_{n + 1} & = & \frac{v_{n - 1}}{1.03 + 2 u_{n - 1} u_{n}}, \end{matrix}$ (12) with initial condition u₀ = 0.2, v₀ = 0.4, u_-1 = 0.5, v_-1 = 0.05. The plot of u_n for the System (13) is shown in Fig. 3, the plot of v_n is shown in Fig. 4 and an attractor of the System (3) is shown in Fig. 5.

Fig.3

Plot of u_n for the System (13).

Fig.4

Plot of v_n for the System (13).

Fig.5

An attractor of System (13).

5 Conclusion

The present work is related to the qualitative nature of the positive solutions of a second-order rational fuzzy difference equations with initial conditions as well as the parameters are positive fuzzy numbers. Furthermore, the existence of positive solutions, boundedness, persistence, and the parametric conditions for stability analysis are investigated. Due to the diverse applications of fuzzy difference equations in the applied sciences, the researchers are more interested to overcome the problems faced in this numerical technique. Since it is well known that, in fuzzy difference equations parameters are positive fuzzy numbers and problems due to these parameters may be handle through a strong and precise parameterized fuzzy mathematical tools. Therefore, well-known parameterized mathematical tools known as fuzzy soft set [25] and soft fuzzy rough sets (tool used for reduction of parameters) [24 , 26–28] can be used to overcome the upcoming difficulties faced in the aforementioned numerical technique.

Footnotes

Acknowledgments

This work was partially supported by the Higher Education Commission of Pakistan.

References

Grove

E.A.

and Ladas

, Periodicities in Nonlinear Difference Equations, Chapman and Hall/CRC Press, Boca Raton, 2004.

Sedaghat

, Nonlinear Difference Equations: Theory with Applications to Social Science Models, Springer, Netherland, 2003.

Din

, Dynamics of a discrete-Lotka-Volterra model, Advances in Difference Equations (2013), 95.

Din

and Donchev

, Global character of a host-parasite model, Chaos, Solitons & Fractals 54 (2013), 1–7.

Din

, Global stability of a population model, Chaos, Solitons & Fractals 59 (2014), 119–128.

Din

, Ibrahim

T.F.

and Khan

K.A.

, Behavior of a competitive system of second-order difference equations, The Scientific World Journal 2014 (2014), 9. Article ID 283982.

Din

and Elsayed

E.M.

, Stability analysis of a discrete ecological model, Computational Ecology and Software 4(2) (2014), 89–103.

Lakshmikantham

and Vatsala

A.S.

, Basic theory of fuzzy difference equations, Journal of Difference Equations & Applications 8(11) (2002), 957–968.

and Zhang

, Embedding problem of non-compact fuzzy number space

\tilde{E}

, Fuzzy Sets & Systems 1 (2000), 110, 135–142.

10.

Bassanezi

R.C.

, Barros

L.C.D.

and Tonelli

P.A.

, Attractors and asymptotic stability for fuzzy dynamical systems, Fuzzy Sets & Systems 113(3) (2000), 473–483.

11.

Papaschinopoulos

and Papadopoulos

B.K.

, On the fuzzy difference equation

x_{n + 1} = A + \frac{B}{x_{n}}

, Soft Computing 6 (2002), 456–461.

12.

Papaschinopoulos

and Papadopoulos

B.K.

, On the fuzzy difference equation

x_{n + 1} = A + \frac{x_{n}}{x_{n - m}}

, Fuzzy Sets & Systems 129 (2002), 73–81.

13.

Papaschinopoulos

and Stefanidou

, Boundedness and asymptotic behavior of the solutions of a fuzzy difference equation, Fuzzy Sets & Systems 140 (2003), 523–539.

14.

Stefanidou

and Papaschinopoulos

, A fuzzy difference equation of rational form, Journal of Nonlinear Mathematical Physics 12(2) (2005), 300–315.

15.

Zhang

Q.H.

, Yang

L.H.

and Liao

D.X.

, Behavior of solutions to a fuzzy nonlinear difference equation, Iranian Journal of Fuzzy Systems 9(2) (2012), 1–12.

16.

Bajo

and Liz

, Global behavior of a second-order nonlinear difference equation, Journal of Difference Equations & Applications 17(10) (2011), 1471–1486.

17.

Pituk

, More on Poincare’s and Perron’s theorems for difference equations, Journal of Difference Equations & Applications 8 (2002), 201–216.

18.

Din

, Asymptotic behavior of an anti-competitive system of second-order difference equations, Journal of the Egyptian Mathematical Society 24 (2016), 37–43.

19.

Din

, Stability analysis of a biological network, Network Biology 4(3) (2014), 123–129.

20.

Din

, Khan

K.A.

and Nosheen

, Stability analysis of a system of exponential difference equations, Discrete Dynamics in Nature and Society 2014 (2014), 11. Article ID 375890.

21.

Din

, On a system of a fourth-order rational difference equations, Acta Universitatis Apulensis 39 (2014), 137–150.

22.

Moye

L.A.

and Kapadia

A.S.

23.

Umekken

S.A.

, Can

and Barak

M.A.

, Fuzzy difference equations in finance, International Journal of Scientific and Innovative Mathematical Research 2(8) (2014), 729–735.

24.

Zhan

, Irfan Ali

and Mehmood

, On a novel uncertain soft set model: Z-soft fuzzy rough set model and corresponding decision making method, Applied Soft Computing 56 (2017), 446–457.

25.

, Liu

and Zhan

, A survey of decision making methods based on certain hybrid soft set models, Artificial Intelligence Review 47 (2017), 507–530.

26.

Zhan

, Liu

and Herawan

, A novel soft rough set: Soft rough hemirings and its multicriteria group decision making, Applied Soft Computing 54 (2017), 393–402.

27.

Zhan

and Zhu

, A novel soft rough fuzzy set: Z-soft rough fuzzy ideals of hemirings and corresponding decision making, Soft Computing 21 (2017), 1923–1936.

28.

Zhan

, Liu

and Zhu

, Another approach to rough soft hemirings and corresponding decision making, Soft Computing 21 (2017), 3769–3780.

29.

Zadeh

L.A.

, Fuzzy sets, Information and Control 8 (1965), 338–353.