
In this paper, we investigate the qualitative behavior of the fuzzy difference equation
zn+1=Azn−sB+Cs∏i=0zn−i
where n∈N0=N∪{0},(zn) is a sequence of positive fuzzy numbers, A,B,C and the initial conditions z−j,j=0,1,...,s are positive fuzzy numbers and s is a positive integer. Moreover, two examples are given to verify the effectiveness of the results obtained.
Citation: İbrahim Yalçınkaya, Durhasan Turgut Tollu, Alireza Khastan, Hijaz Ahmad, Thongchai Botmart. Qualitative behavior of a higher-order fuzzy difference equation[J]. AIMS Mathematics, 2023, 8(3): 6309-6322. doi: 10.3934/math.2023319
[1] | Qianhong Zhang, Ouyang Miao, Fubiao Lin, Zhongni Zhang . On discrete-time laser model with fuzzy environment. AIMS Mathematics, 2021, 6(4): 3105-3120. doi: 10.3934/math.2021188 |
[2] | Iqbal Ahmad, Mohd Sarfaraz, Syed Shakaib Irfan . Common solutions to some extended system of fuzzy ordered variational inclusions and fixed point problems. AIMS Mathematics, 2023, 8(8): 18088-18110. doi: 10.3934/math.2023919 |
[3] | Omar Kahouli, Amina Turki, Mohamed Ksantini, Mohamed Ali Hammami, Ali Aloui . On the boundedness of solutions of some fuzzy dynamical control systems. AIMS Mathematics, 2024, 9(3): 5330-5348. doi: 10.3934/math.2024257 |
[4] | Sumbal Ahsan, Rashid Nawaz, Muhammad Akbar, Saleem Abdullah, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar . Numerical solution of system of fuzzy fractional order Volterra integro-differential equation using optimal homotopy asymptotic method. AIMS Mathematics, 2022, 7(7): 13169-13191. doi: 10.3934/math.2022726 |
[5] | Zena Talal Yassin, Waleed Al-Hayani, Ali F. Jameel, Ala Amourah, Nidal Anakira . Solving fuzzy system of Fredholm integro-differential equations of the second kind by using homotopy analysis method. AIMS Mathematics, 2025, 10(1): 1704-1740. doi: 10.3934/math.2025078 |
[6] | Tareq Eriqat, Rania Saadeh, Ahmad El-Ajou, Ahmad Qazza, Moa'ath N. Oqielat, Ahmad Ghazal . A new analytical algorithm for uncertain fractional differential equations in the fuzzy conformable sense. AIMS Mathematics, 2024, 9(4): 9641-9681. doi: 10.3934/math.2024472 |
[7] | Li Chen, Suyun Wang, Yongjun Li, Jinying Wei . New results for fractional ordinary differential equations in fuzzy metric space. AIMS Mathematics, 2024, 9(6): 13861-13873. doi: 10.3934/math.2024674 |
[8] | Muhammad Akram, Ghulam Muhammad, Tofigh Allahviranloo, Ghada Ali . A solving method for two-dimensional homogeneous system of fuzzy fractional differential equations. AIMS Mathematics, 2023, 8(1): 228-263. doi: 10.3934/math.2023011 |
[9] | Saima Rashid, Fahd Jarad, Khadijah M. Abualnaja . On fuzzy Volterra-Fredholm integrodifferential equation associated with Hilfer-generalized proportional fractional derivative. AIMS Mathematics, 2021, 6(10): 10920-10946. doi: 10.3934/math.2021635 |
[10] | Nurain Zulaikha Husin, Muhammad Zaini Ahmad . Hybridization of the shooting and Runge-Kutta Cash-Karp methods for solving Fuzzy Boundary Value Problems. AIMS Mathematics, 2024, 9(11): 31806-31847. doi: 10.3934/math.20241529 |
In this paper, we investigate the qualitative behavior of the fuzzy difference equation
zn+1=Azn−sB+Cs∏i=0zn−i
where n∈N0=N∪{0},(zn) is a sequence of positive fuzzy numbers, A,B,C and the initial conditions z−j,j=0,1,...,s are positive fuzzy numbers and s is a positive integer. Moreover, two examples are given to verify the effectiveness of the results obtained.
Difference equations play an important role in modeling problems that arise in biology, physics, engineering, finance and many other areas [1,2]. In many cases, the obtained models are restricted in their ability to describe phenomena due to the incomplete or uncertain information about the variables, parameters and initial conditions available. To take into account these uncertainties or lack of precision, we may use a fuzzy environment in variables, parameters and initial conditions, by turning general difference equations into fuzzy difference equations. The most striking applications of the fuzzy notion can be seen in some branches of engineering. See [3,4] for a few applications from mechanical engineering. The fuzzy notion also has applications in differential equations, which are closely related to difference equations. See [5]. Most numerical methods convert differential equations to difference equations, and their relationship is discussed in [6].
In [7], Bajo and Liz investigated the behavior of the ordinary difference equation
xn+1=xn−1a+bxn−1xn, | (1.1) |
where the parameters a,b and the initial conditions x−1,x0 are real numbers.
Moreover, in [8], Rahman et al. investigated the behavior of the fuzzy difference equation
xn+1=xn−1A+Bxn−1xn, | (1.2) |
where A,B and the initial conditions x−1,x0 are positive fuzzy numbers.
In [9], Shojaei et al. investigated the behavior of the ordinary difference equation
xn+1=αxn−2β+γxn−2xn−1xn, | (1.3) |
where the parameters α,β,γ and the initial conditions x−2,x−1,x0 are real numbers.
Also, in [10], Atak et al. investigated the behavior of the fuzzy difference equation
zn+1=zn−2C+zn−2zn−1zn , | (1.4) |
where (zn) is a sequence of positive fuzzy numbers, C and the initial conditions z−2,z−1,z0 are positive fuzzy numbers.
Motivated by above mentioned studies, in this paper, we investigate the qualitative behavior of the fuzzy difference equation
zn+1=Azn−sB+Cs∏i=0zn−i , | (1.5) |
where n∈N0=N∪{0},(zn) is a sequence of positive fuzzy numbers, A,B,C and the initial conditions z−j(j=0,1,...,s) are positive fuzzy numbers and s is a positive integer. Moreover, some examples are given to verify the effectiveness of the results obtained.
Recently some researchers have studied fuzzy difference equations and properties of their solutions in different approaches. See, for example, [11,12,13,14,15,16,17,18,19,20]. We also refer to [21,22] for some fundamentals of ordinary difference equations.
In this section, we give some definitions which will be used in this paper. For more details see [23,24,25].
Definition 2.1. Consider a fuzzy set A which is a function from the set of real numbers R into the interval [0,1]. We say that A is a fuzzy number if it satisfies the following properties
(a) A is normal, i.e., ∃x0∈R with A(x0)=1,
(b) A is fuzzy convex, i.e., A(tx1+(1−t)x2)≥min{A(x1),A(x2)},∀t∈[0,1] and x1,x2∈R,
(c) A is upper semicontinuous on R,
(d) A is compactly supported, i.e., ¯{x∈R:A(x)>0} is compact.
Let us denote by RF the space of all fuzzy numbers. For 0<α≤1 and A∈RF, we denote α-cuts of fuzzy number A by [A]α={x∈R,A(x)≥α} and [A]0=¯{x∈R,A(x)≥0}. We call [A]0, the support of fuzzy number A and denote it by supp(A).
The fuzzy number A is called positive if supp(A)⊂(0,∞). We denote by R+F, the space of all positive fuzzy numbers.
Definition 2.2. (a) Let A,B∈RF with [A]α=[Aαl,Aαr] and [B]α=[Bαl,Bαr] for α∈[0,1]. We define ||A|| on the space of fuzzy numbers as follow;
||A||=supmax{|Aαl|,|Aαr|}, |
where sup is taken for all α∈[0,1]. We recall the following metric
D(A,B)=sup{max{|Aαl−Bαl|,|Aαr−Bαr|}}, |
where sup is taken for all α∈[0,1].
(b) Let (xn) be a sequence of positive fuzzy numbers and x∈RF. Then, we say that
lim n→∞xn=xiff lim n→∞D(xn,x)=0. |
The following lemma and definition are given in [24].
Lemma 2.1. Let X,Y∈RF and [X]α=[Xαl,Xαr],[Y]α=[Yαl,Yαr] for α∈[0,1] be the α-cuts of X,Y, respectively. Let Z be a fuzzy number such that [Z]α=[Zαl,Zαr] for α∈[0,1]. Then, MIN{X,Y}=Z (resp. MAX{X,Y}=Z) if and only if min{Xαl,Yαl}=Zαl and min{Xαr,Yαr}=Zαr (resp. max{Xαl,Yαl}=Zαl and max{Xαr,Yαr}=Zαr).
Definition 2.3. (a) We say that a sequence of positive fuzzy numbers (xn) is bounded and persistent if there exist n0∈N and C,D∈R+F such that MIN{xn,C}=C and MAX{xn,D}=D for n≥n0.
(b) We say that (xn) for n∈N0 is an unbounded sequence if the ||xn|| for n∈N0 is an unbounded sequence.
We need the following lemma which is given in [18] in the next section.
Lemma 2.2. Let f be a continuous function from R+×R+×...×R+ into R+ and B0,B1,...,Bk∈RF. Then,
[f(B0,B1,...,Bk)]α=f([B0]α,[B1]α,...,[Bk]α)forα∈(0,1]. |
In this section, we study the existence and properties of the positive solutions to (1.5). If (zn) is a sequence of positive fuzzy numbers which satisfies (1.5), we say (zn) is a positive solution of (1.5).
Theorem 3.1. Consider (1.5) where A,B,C∈R+F. Then, for any positive fuzzy numbers z−j(j=0,1,...,s) there exists a unique positive solution (zn) of (1.5) with the initial conditions z−j (j=0,1,...,s).
Proof. Suppose that there exists a sequence of fuzzy numbers (zn) satisfying (1.5) with the initial conditions z−j (j=0,1,...,s). Consider the α-cuts
[zn]α=[Lαn,Rαn], [A]α=[Aαl,Aαr],[B]α=[Bαl,Bαr],[C]α=[Cαl,Cαr] | (3.1) |
for n=−s,−s+1,... and all α∈(0,1]. Then, from (1.5)-(3.1) and Lemma 2.2, it follows that
[zn+1]α=[Azn−sB+Cs∏i=0zn−i]α=[A]α[zn−s]α[B]α+[C]αs∏i=0[zn−i]α=[Aαl,Aαr][Lαn−s,Rαn−s][Bαl,Bαr]+[Cαl,Cαr]s∏i=0[Lαn−i,Rαn−i]=[AαlLαn−sBαr+Cαrs∏i=0Rαn−i,AαrRαn−sBαl+Cαls∏i=0Lαn−i] |
from which we have
Lαn+1=AαlLαn−sBαr+Cαrs∏i=0Rαn−i, Rαn+1=AαrRαn−sBαl+Cαls∏i=0Lαn−i, | (3.2) |
for n∈N0 and α∈(0,1]. Then, it is clear that for any (Lαj,Rαj),j=−s,−s+1,...,0, there exists a unique solution (Lαn,Rαn) with the initial conditions (Lαj,Rαj),j=−s,−s+1,...,0 for all α∈(0,1].
Now, we prove that [Lαn,Rαn] for all α∈(0,1], where (Lαn,Rαn) is the solution of the system (3.2) with the initial conditions (Lαj,Rαj),j=−s,−s+1,...,0, determines the solution (zn) of (1.5) with the initial conditions z−j (j=0,1,...,s) such that
[zn]α=[Lαn,Rαn], α∈(0,1], n=−s,−s+1,.... | (3.3) |
Since A,B,C,z−j(j=0,1,...,s)∈R+F for any α1,α2∈(0,1] and α1≤α2, we have
0<Aα1l≤Aα2l≤Aα2r≤Aα1r,0<Bα1l≤Bα2l≤Bα2r≤Bα1r,0<Cα1l≤Cα2l≤Cα2r≤Cα1r,0<Lα1j≤Lα2j≤Rα2j≤Rα1j, | (3.4) |
for j=−s,−s+1,...,0. We prove by the induction that
0<Lα1n≤Lα2n≤Rα2n≤Rα1n, n∈N0. | (3.5) |
From (3.4), we see that (3.5) holds for n=−s,−s+1,...,0. Suppose that (3.5) is valid for n≤k,k∈{1,2,...}. Then, from (3.2), (3.4) and (3.5) for n≤k, it follows that
Lα1k+1=Aα1lLα1k−sBα1r+Cα1rs∏i=0Rα1k−i≤Aα2lLα2k−sBα2r+Cα2rs∏i=0Rα2k−i=Lα2k+1≤Aα2rRα2k−sBα2l+Cα2ls∏i=0Lα2k−i=Rα2k+1≤Aα1rRα1k−sBα1l+Cα1ls∏i=0Lα1k−i=Rα1k+1. |
Therefore, (3.5) is satisfied. Moreover, from (3.2), we have
Lα1=AαlLα−sBαr+Cαrs∏i=0Rα−i, Rα1=AαrRα−sBαl+Cαls∏i=0Lα−i, for all α∈(0,1]. | (3.6) |
Then, since A,B,C,z−j (j=0,1,...,s)∈R+F, we have that Aαl,Aαr,Bαl,Bαr,Cαl,Cαr,Lα−j,Rα−j(j=0,1,...,s) are left continuous. So, from (3.6) we see that Lα1 and Rα1 are also left continuous. Working inductively we can easily obtain that Lαn and Rαn are left continuous for n∈N.
Now, we prove that ¯∪α∈(0,1][Lαn,Rαn] is compact. It is sufficient to prove that ∪α∈(0,1][Lαn,Rαn] is bounded. Let n=1, since A,B,C,z−j (j=0,1,...,s)∈R+F, there exist constants
MA,NA,MB,NB,MC,NC,M−j,N−j>0 |
for j=0,1,...,s such that
[Aαl,Aαr]⊂[MA,NA],[Bαl,Bαr]⊂[MB,NB],[Cαl,Cαr]⊂[MC,NC],[Lα−j,Rα−j]⊂[M−j,N−j] | (3.7) |
for j=0,1,...,s. Therefore, from (3.6) and (3.7) we can easily prove that
[Lα1,Rα1]⊂[MAM−sNB+Ncs∏i=0N−i,NAN−sMB+Mcs∏i=0M−i] | (3.8) |
from which it is clear that
∪α∈(0,1][Lα1,Rα1]⊂[MAM−sNB+Ncs∏i=0N−i,NAN−sMB+Mcs∏i=0M−i] for all α∈(0,1]. | (3.9) |
Also, (3.9) implies that ¯∪α∈(0,1][Lα1,Rα1] is compact and ¯∪α∈(0,1][Lα1,Rα1]⊂(0,∞). Working inductively, we can easily see that ¯∪α∈(0,1][Lαn,Rαn] is compact and
¯∪α∈(0,1][Lαn,Rαn]⊂(0,∞), for n∈N. | (3.10) |
Therefore, using (3.5), (3.10) and since Lαn,Rαn are left continuous, we see that [Lαn,Rαn] determines a sequence of positive fuzzy numbers (zn) such that (3.3) holds.
We prove now that (zn) is the solution of (1.5) with the initial conditions z−j (j=0,1,...,s). Since
[zn+1]α=[Lαn+1,Rαn+1]=[AαlLα−sBαr+Cαrs∏i=0Rα−i,AαrRαn−sBαl+Cαls∏i=0Lαn−i]=[Azn−sB+Cs∏i=0zn−i]α |
for all α∈(0,1], we have that (zn) is the solution of (1.5) with the initial conditions z−j (j=0,1,...,s).
Suppose that there exists another solution (˜zn) of (1.5) with the initial conditions z−j (j=0,1,...,s). Then, we can easily show that
[˜zn]α=[Lαn,Rαn] for α∈(0,1] and n∈N0. | (3.11) |
Then, from (3.3) and (3.11), [zn]α=[˜zn]α for α∈(0,1] and n=−s,−s+1,... from which we get zn=˜zn for n=−s,−s+1,.... Thus, the proof is completed.
In order to study the further dynamics of (1.5), we apply the results concerning the following system of ordinary difference equations
un+1=a1un−sb1+c1s∏i=0vn−i, vn+1=a2vn−sb2+c2s∏i=0un−i, n∈N0, | (3.12) |
where the parameters a1,b1,c1,a2,b2,c2 and initial conditions u−j,v−j(j=0,1,...,s) are positive real numbers and s is a positive integer. Note that system (3.12) can be written as
xn+1=xn−sq+s∏i=0yn−i, yn+1=yn−sp+s∏i=0xn−i, n∈N0, | (3.13) |
by the change of variables un=(a2c2)1s+1xn,vn=(a1c1)1s+1yn with q=b1a1 and p=b2a2. The equilibrium points of (3.13) are the solutions of the equations
¯x=¯xq+s∏i=0¯y, ¯y=¯yp+s∏i=0¯x. |
We need the following results concerning the behavior of the solutions of the system (3.13) which has been presented in [26].
If p<1,q<1, then (3.13) has equilibrium points (0,0) and (s+1√1−p,s+1√1−q). In addition, if p<1,q=1, then system (3.13) has an equilibrium (s+1√1−p,0) and if p=1,q<1, then system (3.13) has an equilibrium (0,s+1√1−q). Finally, if p>1,q>1, then system (3.13) has the unique equilibrium (0,0).
Theorem 3.2. Let (xn,yn) be any positive solution of system (3.13), then the following statements are true:
(1) For every m≥0, the following results hold:
0≤xn≤{(1q)m+1x−s, n=(s+1)m+1, (1q)m+1x−s+1,n=(s+1)m+2, ⋮⋮(1q)m+1x0, n=(s+1)m+s+1 | (3.14) |
and
0≤yn≤{(1p)m+1y−s, n=(s+1)m+1, (1p)m+1y−s+1,n=(s+1)m+2, ⋮⋮(1p)m+1y0, n=(s+1)m+p+1. | (3.15) |
(2) If p<1 and q<1, then the following statements are true for j=−s,−s+1,...,0.
(i) If (xj,yj)∈(0,s+1√1−p)×(s+√1−q,∞), then (xn,yn)∈(0,s+1√1−p)×(s+√1−q,∞).
(ii) If (xj,yj)∈(s+1√1−p,∞)×(0,s+1√1−q), then (s+1√1−p,∞)×(0,s+1√1−q).
(3) If p>1 and q>1, then every positive solution (xn,yn) of system (3.13) converges to (0,0) as n→∞.
The following corollary can be obtained from Theorem 3.2.
Corollary 3.1. Let (xn,yn) be positive solution to system (3.13), then the following statements are true.
(1) If p>1 and q>1, then every positive solution (xn,yn) of system (3.13) is bounded and persistent.
(2) If p<1 and q<1, then the system (3.13) has unbounded solutions:
(i) If xj∈(0,s+1√1−p)andyj∈(s+1√1−q,∞) for j=−s,−s+1,...,0, then limn→∞xn=0 and limn→∞yn=∞.
(ii) If xj∈(s+1√1−p,∞) and yj∈(0,s+1√1−q) for j=−s,−s+1,...,0, then limn→∞yn=0 and limn→∞xn=∞.
Proof. (1) Let β=max{x−s,x−s+1,...,x0} and γ=max{y−s,y−s+1,...,y0}, then it is easy to see from system (3.13) that
0≤xn≤β and 0≤yn≤γ |
for all n∈N0. Hence, every positive solution (xn,yn) of system (3.13) is bounded and persistent.
(2) We only prove the part (i), since the part (ii) can be proved similarly.
(i) Assume that (xn,yn) is a positive solution of system (3.13) such that xj∈(0,s+1√1−p)andyj∈(s+√1−q,∞) for j=−s,−s+1,...,0. Then, from system (3.13), we obtain the following inequalities
x1=x−sq+s∏i=0y−i<x−sq+s∏i=0s+1√1−q=x−s,x2=x−s+1q+s∏i=0y−i+1<x−s+1q+s∏i=0s+1√1−q=x−s+1,x3=x−s+2q+s∏i=0y−i+2<x−s+2q+s∏i=0s+1√1−q=x−s+2,⋮ |
and
y1=y−sp+s∏i=0x−i>y−sp+s∏i=0s+1√1−p=y−s,y2=y−s+1p+s∏i=0x−i+1>y−s+1p+s∏i=0s+1√1−p=y−s+1,y3=y−s+2p+s∏i=0x−i+2>y−s+2p+s∏i=0s+1√1−p=y−s+2,⋮ |
from which it follows that
limn→∞xn=0 andlimn→∞yn=∞. |
So, the proof is completed.
Theorem 3.3. Consider the fuzzy difference Eq (1.5). The the following statements are true.
(1) If Aαr<Bαl for all α∈(0,1], then every positive solution of (1.5) is bounded and persistent.
(2) If there exists an ¯α∈(0,1] such that B¯αr<A¯αl, then the Eq (1.5) has unbounded solutions.
Proof. (1) Consider the following system of ordinary difference equations
sn+1=γAsn−pβB+βCp∏i=0tn−i, tn+1=βAtn−pγB+γCp∏i=0sn−i, n∈N0 | (3.16) |
where
[A]α=[Aαl,Aαr]⊂¯∪α∈(0,1][Aαl,Aαr]⊂[γA,βA],[B]α=[Bαl,Bαr]⊂¯∪α∈(0,1][Bαl,Bαr]⊂[γB,βB],[C]α=[Cαl,Cαr]⊂¯∪α∈(0,1][Cαl,Cαr]⊂[γC,βC]. | (3.17) |
Let (sn,tn) be a solution of system (3.16) with the initial conditions (s−j,t−j)=(γ−j,β−j) for j=0,1,...,p where γ−j and β−j are given
[Lα−j,Rα−j]⊂¯∪α∈(0,1][Lα−j,Rα−j]⊂[γ−j,β−j] for j=−p,−p+1,...,0. | (3.18) |
Then, from (3.16) and (3.17) it folows that
s1=γAs−pβB+βCp∏i=0t−i≤AαlLα−pBαr+Cαrp∏i=0Rα−i=Lα1 | (3.19) |
and
t1=βAt−pγB+γCp∏i=0s−i≥AαrR−p,αBαl+Cαlp∏i=0Lα−i=Rα1. | (3.20) |
Hence, by induction, we get sn≤Lαn and Rαn≤tn for n∈</p><p>N</p><p>. Assume that Aαr<Bαl for all α∈(0,1], then it follows that γA<βB and βA<γB. From (2) of Theorem 3.2, the solution (sn,tn) of system (3.15) is bounded and persistent, which is the solution (zn) of (1.5). This completes the proof of (1).
(2) Suppose that there exists an ¯α∈(0,1] such that B¯αr<A¯αl. If A¯αl=a1,A¯αr=a2,B¯αl=b2,B¯αr=b1,L¯αn=xn and R¯αn=yn for n=−s,−s+1,..., then we can apply (i) of (2) in Corollary 3.1 to system (3.2) (We can use (ii) of (2) in Corollary 3.1, too). If there exists an ¯α∈(0,1] such that B¯αr<A¯αl and xj∈(0,s+1√1−p)andyj∈(s+1√1−q,∞) for j=−s,−s+1,...,0, then there exist solutions (xn,yn) of system (3.13) where ¯α=α with initial conditions (x−j,y−j) for j=0,1,...,s such that
lim n→∞xn=0 and lim n→∞yn=∞. | (3.21) |
Moreover, if x−j<y−j(j=0,1,...,s), we can find z−j∈R+F such that
[zj]α=[Lαj,Rαj] for α∈(0,1] | (3.22) |
and
[zj]¯α=[L¯αj,R¯αj]=[xj,yj] | (3.23) |
for j=−s,−s+1,...,0. Let (zn) be a positive solution of (1.5) with the initial conditions z−j(j=0,1,...,s) and [zn]α=[Lαn,Rαn] for α∈(0,1]. Since (3.22) and (3.23) hold and (Lαn,Rαn) satisfies system (3.2), we have
[zn]¯α=[L¯αn,R¯αn]=[xn,yn]. | (3.24) |
Therefore, from (3.21), (3.24) and since
||zn||=supα∈(0,1]max{|Lαn|,|Rαn|}≥max{|L¯αn|,|R¯αn|}=R¯αn |
where sup is taken for all α∈(0,1], it is clear that solution (zn) is unbounded. This completes the proof of (2).
Theorem 3.4. If Aαr<Bαl for all α∈(0,1], then every positive solution (zn) of (1.5) converges to 0 as n→∞.
Proof. Let (zn) be a positive solution of (1.5) such that (3.3) holds with Aαr<Bαl for all α∈(0,1]. Then, we can apply (3) of Theorem 3.2 to system (3.2). So, we get
lim n→∞Lαn=lim n→∞Rαn=0. | (3.25) |
Therefore, from (3.25), we get
lim n→∞D(zn,0)=lim n→∞(supα∈(0,1]{max{|Lαn−0|,|Rαn−0|}})=0. |
This completes the proof.
In this section, to verify obtained results, we give two numerical examples for s=3 with different values of A,B,C where the initial conditions z−3,z−2,z−1,z0 are satisfied
z−3(x)={4x−0.42,0.1≤x≤0.6,4.4−4x2,0.6≤x≤1.1,z−2(x)={5x−12,0.2≤x≤0.6,5−5x2,0.6≤x≤1,z−1(x)={4x−12,0.25≤x≤0.75,5−4x2,0.75≤x≤1.25,z0(x)={5x−2.52,0.5≤x≤0.9,6.5−5x2,0.9≤x≤1.3. | (4.1) |
From (4.1), we get
[z−3]α=[2α+0.44,4.4−2α4],[z−2]α=[2α+15,5−2α5],[z−1]α=[2α+14,5−2α4],[z0]α=[2α+2.505,6.50−2α5] |
for all α∈[0,1].
Example 4.1. Consider Eq (1.5) where the initial conditions are satisfied (4.1) and A,B,C are satisfied
A={4x−1,0.25≤x≤0.5,3−4x,0.5≤x≤0.75,B={x−1,1≤x≤2,3−x,2≤x≤3,C={2x−1,0.5≤x≤1,3−2x,1≤x≤1.5. | (4.2) |
Then, from (4.2), we get [A]α=[α+14,3−α4],[B]α=[α+1,3−α] and [C]α=[α+12,3−α2] for all α∈(0,1]. By Theorem 3.1, there exists a unique solution. Since Aαr<Bαl for all α∈[0,1], then by case (1) in Theorem 3.3, the positive solution (zn) of fuzzy difference Eq (1.5) is bounded and persistent and by Theorem 3.4, it converges to 0 as n→∞. For α1=0.2 and α2=0.8, the α−cuts of the solution Lαn and Rαn are depicted in Figures 1 and 2, respectively.
Example 4.2. Consider Eq (1.5) where the initial conditions are satisfied (4.1) and A,B,C are satisfied
A={x−2,2≤x≤3,4−x,3≤x≤4,B={4x−1,0.25≤x≤0.5,3−4x,0.5≤x≤0.75,C={x−1,1≤x≤2,3−x,2≤x≤3. | (4.3) |
Then, from (4.3), we get [A]α=[α+2,4−α],[B]α=[α+14,3−α4] and [C]α=[α+1,3−α] for all α∈(0,1]. By Theorem 3.1 there exists a unique positive solution. For any α∈[0,1], we have Aαl>Bαr. So, by case (2) in Theorem 3.3, the corresponding fuzzy difference equation has unbounded solutions. For α1=0.2 and α2=0.8, the α-cuts of the solution Lαn and Rαn are depicted in Figures 3 and 4, respectively.
In this study, we investigated behavior of the fuzzy difference equation zn+1=Azn−s/(B+Cs∏i=0zn−i). We have shown that, under certain conditions, the positive solutions of this equation converge to zero. We have also considered the case where the solutions are unbounded. Finally, we have supported our theoretical results via two numerical examples. This study extends the results in the references [8,10].
The authors declare that they have no conflicts of interest.
[1] | R. P. Agarwal, Difference equations and inequalities, New York, 1993. |
[2] | E. P. Popov, Automatic regulation and control, Moscow, 1966. |
[3] |
M. Bakır, S. Akan, E. Özdemir, Regional aircraft selection with fuzzy piprecia and fuzzy marcos: a case study of the turkish airline industry, Facta Univ. Ser. Mech., 19 (2021), 423–445. https://doi.org/10.22190/FUME210505053B doi: 10.22190/FUME210505053B
![]() |
[4] |
D. Božanić, A. Milić, D. Tešić, W. Salabun, D. Pamučar, D numbers-FUCOM-Fuzzy RAFSI model for selecting the group of construction machines for enabling mobility, Facta Univ. Ser. Mech., 19 (2021), 447–471. https://doi.org/10.22190/FUME210318047B doi: 10.22190/FUME210318047B
![]() |
[5] |
F. Rabiei, F. A. Hamid, M. Rashidi, Z. Ali, K. Shah, K. Hosseini, et al., Numerical simulation of fuzzy Volterra integro-differential equation using improved Runge-Kutta method, J. Appl. Comput. Mech., 9 (2023), 72–82. https://doi.org/10.22055/JACM.2021.38381.3212 doi: 10.22055/JACM.2021.38381.3212
![]() |
[6] |
J. H. He, F. Y. Ji, H. Mohammad-Sedighi, Difference equation vs differential equation on different scales, Int. J. Numer. Method. H., 31 (2021), 391–401. https://doi.org/10.1108/HFF-03-2020-0178 doi: 10.1108/HFF-03-2020-0178
![]() |
[7] |
I. Bajo, E. Liz, Global behavior of a second-order non-linear difference equation, J. Differ. Equ. Appl., 17 (2011), 1471–1486. https://doi.org/10.1080/10236191003639475 doi: 10.1080/10236191003639475
![]() |
[8] |
G. Rahman, Q. Din, F. Faizullah, F. M. Khan, Qualitative behavior of a second-order fuzzy difference equation, J. Intell. Fuzzy Syst., 34 (2018), 745–753. https://doi.org/10.3233/JIFS-17922 doi: 10.3233/JIFS-17922
![]() |
[9] |
M. Shojaei, R. Saadeti, H. Adibi, Stability and periodic character of a rational third-order difference equation, Chaos Soliton. Fract., 39 (2009), 1203–1209. https://doi.org/10.1016/j.chaos.2007.06.029 doi: 10.1016/j.chaos.2007.06.029
![]() |
[10] | I. Yalcinkaya, N. Atak, D. T. Tollu, On a third-order fuzzy difference equation, J. Prime Res. Math., 17 (2021), 59–69. |
[11] |
K. A. Chrysafis, B. K. Papadopoulos, G. Papaschinopoulos, On the fuzzy difference equations of finance, Fuzzy Set. Syst., 159 (2008), 3259–3270. https://doi.org/10.1016/j.fss.2008.06.007 doi: 10.1016/j.fss.2008.06.007
![]() |
[12] |
E. Deeba, A. De Korvin, Analysis by fuzzy difference equations of a model of CO2 level in blood, Appl. Math. Lett., 12 (1999), 33–40. https://doi.org/10.1016/S0893-9659(98)00168-2 doi: 10.1016/S0893-9659(98)00168-2
![]() |
[13] | E. Hatir, T. Mansour, I. Yalcinkaya, On a fuzzy difference equation, Utilitas Mathematica, 93 (2014), 135–151. |
[14] |
A. Khastan, Fuzzy logistic difference equation, Iran. J. Fuzzy Syst., 15 (2018), 55–66. https://doi.org/10.22111/IJFS.2018.4281 doi: 10.22111/IJFS.2018.4281
![]() |
[15] | A. Khastan, Z. Alijani, On the new solutions to the fuzzy difference equation xn+1=A+B/xn, Fuzzy Set. Syst., 358 (2019), 64–83. https://doi.org/10.1016/j.fss.2018.03.014 |
[16] | G. Papaschinopoulos, B. K. Papadopoulos, On the fuzzy difference equation xn+1=A+B/xn, Soft Computing, 6 (2002), 456–461. https://doi.org/10.1007/s00500-001-0161-7 |
[17] | G. Papaschinopoulos, B. K. Papadopoulos, On the fuzzy difference equation xn+1=A+xn/xn−m, Fuzzy Set. Syst., 129 (2002), 73–81. https://doi.org/10.1016/S0165-0114(01)00198-1 |
[18] |
G. Papaschinopoulos, G. Stefanidou, Boundedness and asymptotic behavior of the solutions of a fuzzy difference equation, Fuzzy Set. Syst., 140 (2003), 523–539. https://doi.org/10.1016/S0165-0114(03)00034-4 doi: 10.1016/S0165-0114(03)00034-4
![]() |
[19] |
D. T. Tollu, I. Yalcinkaya, H. Ahmad, S. Yao, A detailed study on a solvable system related to the linear fractional difference equation, Math. Biosci. Eng., 18 (2021), 5392–5408. https://doi.org/10.3934/mbe.2021273 doi: 10.3934/mbe.2021273
![]() |
[20] |
E. Deeba, A. De Korvin, E. L. Koh, A fuzzy difference equation with an application, J. Differ. Equ. Appl., 2 (1996), 365–374. https://doi.org/10.1080/10236199608808071 doi: 10.1080/10236199608808071
![]() |
[21] | S. Elaydi, An introduction to difference equations, New York: Springer, 1999. https://doi.org/10.1007/978-1-4757-3110-1 |
[22] | K. L. Kocic, G. Ladas, Global behavior of nonlinear difference equations of higher order with applications, Springer Science & Business Media, 1993. |
[23] | B. Bede, Fuzzy sets, In: Mathematics of fuzzy sets and fuzzy logic, Berlin: Springer, 2013. https://doi.org/10.1007/978-3-642-35221-8_1 |
[24] | G. Klir, B. Yuan, Fuzzy sets and fuzzy logic, New Jersey: Prentice Hall, 1995. |
[25] |
C. Wu, B. Zhang, Embedding problem of noncompact fuzzy number space E∼ (I), Fuzzy Set. Syst., 105 (1999), 165–169. https://doi.org/10.1016/s0165-0114(97)00218-2 doi: 10.1016/s0165-0114(97)00218-2
![]() |
[26] |
Q. Zhang, W. Zhang, On a system of two higher-order nonlinear difference equations, Adv. Math. Phy., 2014 (2014), 729273. https://doi.org/10.1155/2014/729273 doi: 10.1155/2014/729273
![]() |
1. | İbrahim Yalçınkaya, Hamdy El-Metwally, Mustafa Bayram, Durhasan Turgut Tollu, On the dynamics of a higher-order fuzzy difference equation with rational terms, 2023, 27, 1432-7643, 10469, 10.1007/s00500-023-08586-y | |
2. | Mehmet Gümüş, İbrahim Yalçinkaya, Durhasan Turgut Tollu, Dynamic analysis of high-order fuzzy difference equation, 2024, 1598-5865, 10.1007/s12190-024-02280-4 | |
3. | Asifa Tassaddiq, Muhammad Tanveer, Muhammad Usman, Dalal Khalid Almutairi, Rabab Alharabi, On the Global Dynamics of a Fourth-Order Riccati-Type Exponential Fuzzy Difference Equation, 2025, 14, 2075-1680, 118, 10.3390/axioms14020118 | |
4. | Osman Topan, Yasin Yazlik, Sevda Atpinar, Dynamics of fuzzy difference equations system with higher-order, 2025, 1598-5865, 10.1007/s12190-025-02440-0 |