Applications of the Fuzzy Sumudu Transform for the Solution of First Order Fuzzy Differential Equations

In this paper, we study the classical Sumudu transform in fuzzy environment, referred to as the fuzzy Sumudu transform (FST). We also propose some results on the properties of the FST, such as linearity, preserving, fuzzy derivative, shifting and convolution theorem. In order to show the capability of the FST, we provide a detailed procedure to solve fuzzy differential equations (FDEs). A numerical example is provided to illustrate the usage of the FST.


Introduction
Integral transforms constitute fundamental tools in operational calculus.They are mathematical operators that have been used widely in solving many practical problems in applied mathematics, physics and engineering [1][2][3][4].The precursor of the integral transforms is the Fourier transform, which is used to express functions in a finite interval.Thenceforward, there are a number of works on the theories and applications of integral transforms, some of which are Laplace, Mellin and Hankel transforms [5][6][7].Subsequently, the concept of integral transforms was expanded to remove the necessity of finite intervals.Watugala [8,9] has proposed a new integral transform called the Sumudu transform.The Sumudu transform has been used to solve ordinary differential equations in control engineering problems.
Weerakon [10,11] has extended the Sumudu transform on partial differential equations.The work was then continued by Asiru [12], who studied the convolution theorem of the Sumudu transform, which can be expressed in terms of polynomial and convergent infinite series.Next, Belgacem et al. [13] have emphasized the Laplace-Sumudu duality, which is a vital step in establishing new results on the Sumudu transform.For instance, the duality property has been used to invoke a complex inverse of the Sumudu transform, as a Bromwich contour integral formula [14].Furthermore, in [15], the applications of the Sumudu transform on Bessel functions and equations have been explored.The theories and applications of the Sumudu transform have been examined and explored by many authors (see in [16][17][18][19][20][21]).
In general, the Sumudu transform is considered a popular integral transform for solving differential equations.This is due to its unity property, which eases the process of finding solutions.It is also more powerful compared to other integral transforms, as the function transformed is a similitude of the resulting function.
Many real-world problems are modeled by differential equations.However, we cannot be sure that the models is perfect.For example, the initial value of the models might not be known accurately.The initial value may contain some uncertainty quantities, such as "less than x 0 ", "about x 0 " or "more than x 0 ".If this is the case, the classical differential equations cannot be used to handle this situation.Therefore, it is necessary to study other theories in order to overcome this problem.One of the most popular theories for describing this situation is the fuzzy set theory [22].By incorporating fuzziness into classical mathematics, many authors studied fuzzy derivatives [23][24][25][26][27][28], fuzzy differential equations (FDEs) [29][30][31] and fuzzy fractional differential equations (FFDEs) [32][33][34].
Recently, Allahviranloo and Ahmadi [35] have proposed the fuzzy Laplace transform (FLT) and showed its applications to solve FDEs.The FLT is then used to solve second order linear FDEs and the state-space description of fuzzy linear continuous time systems [36].This work has motivated many authors to expand the theories and applications of FLT in the mathematics and engineering fields [37][38][39][40][41].The work has also motivated a few researchers to study the classical Sumudu transform in the fuzzy setting.The first effort was initiated by Ahmad and Abdul Rahman [42] and further studied by Alam Khan et al. [43].In this paper, we add some new results on the Sumudu transform in the fuzzy setting, especially on the linearity and preserving properties.Some other results may parallel the ones proposed in [43].However, our definition of the Sumudu transform in the fuzzy setting is quite general, and the results are presented in different ways.
This paper is organized as follows.In Section 2, we recall several basic definitions and concepts of fuzzy numbers.In Section 3, we provide a general definition of the fuzzy Sumudu transform (FST) and investigate the duality property between FST and FLT.We also provide some theorems and properties regarding the FST.In Section 4, we construct detailed procedures to solve FDEs.Later, in Section 6, we give conclusions.

Preliminaries
In this section, we recall some definitions and theorems needed in order to understand the contribution in this paper.The definition of a fuzzy number is as follows.Definition 1. [44] By R, we denote the set of all real numbers.A fuzzy number is a mapping U : R → [0, 1] with the following properties: (1) U is upper semi-continuous, (2) U is fuzzy convex, i.e., U (λx + (1 − λ)y) ≥ min{U (x), U (y)} for all x, y ∈ R, λ ∈ [0, 1], (3) U is normal, i.e., ∃x 0 ∈ R for which U (x 0 ) = 1, (4) supp U = {x ∈ R|U (x) > 0} is the support of the U , and its closure, i.e. cl(supp U ) is compact.
The following definition is the α-level set of fuzzy numbers.Definition 2. [44] Let F(R) be the set of all fuzzy numbers on R. The α-level set of a fuzzy number U ∈ F(R), α ∈ [0, 1], denoted by U α , is defined as: It is clear that the α-level set of a fuzzy number is a closed and bounded interval, i.e., [u α , u α ], where u α and u α denote the lower bound and the upper bound of U , respectively.Since each y ∈ R can be regarded as a fuzzy number y defined by: R can be embedded in F(R).
Remark.[45] Let X be the Cartesian product of universes X = X 1 × . . .× X n and A 1 , . . ., A n be n fuzzy numbers in X 1 , . . ., X n , respectively.A fuzzy function f maps from X to a universe Y , y = f (x 1 , . . ., x n ).Then, the extension principle allows us to define a fuzzy set B in Y by: where: where f −1 is the inverse of f .For n = 1, the extension principle is reduced to: where: Referring to the extension principle, addition on F(R) can be defined by: and scalar multiplication is defined by: where 0 ∈ F(R).
Furthermore, for all α-levels, and: In this paper, the notation U α represents the α-level set of a fuzzy number.We may conclude that the fuzzy number is determined by the endpoints of the intervals U α .This leads to other representation of a fuzzy number, which will be defined by two endpoint functions u α and u α .Friedman et al. [46] and Ma et al. [47] defined the representation as: which satisfy the following requirements: (1) u α is a bounded non-decreasing left continuous function in (0, 1] and right continuous at zero, (2) u α is a bounded non-increasing left continuous function in (0, 1] and right continuous at zero, (3) u α ≤ u α .
A fuzzy number can be represented as a fuzzy membership function.One of the most commonly-used fuzzy membership functions in the literature is the fuzzy triangular membership function.It is defined as follows.
Definition 4. Let U ∈ F(R).U is called a triangular fuzzy number if its membership function has the following form: and its α-level is simply and k > 0, we define addition, subtraction and multiplication by k for U α and V α as: (1) addition, (2) subtraction, (3) scalar multiplication, Definition 6. [48] The distance D(U, V ) between two fuzzy intervals U and V is defined as: where: is the Hausdorff distance between U α and V α .
Thus, we can conclude that D is a metric space and has the following properties: (1) Then, f (x) is improper fuzzy Riemann-integrable on [a, ∞), and the improper fuzzy Riemann-integrable is a fuzzy number.Furthermore, we have: Proposition 1. [51] If each of f (x) and g(x) is a fuzzy-valued function and fuzzy Riemann-integrable on [a, ∞), then f (x) ⊕ g(x) is fuzzy Riemann-integrable on [a, ∞).Moreover, we have: The next definition is Hukuhara's differentiability, also known as H-derivatives.The definition is about H-differences of sets, and it is introduced as follows.Definition 8. [27] Let x, y ∈ F(R).If there exists z ∈ F(R), such that x = y ⊕ z, then z is called the H-difference of x and y, and it is denoted by x − H y. Definition 9. [52,53] Let f : (a, b) → F(R) and x 0 ∈ (a, b).We say that f is strongly generalized differentiable at x 0 , if there exists an element f (x 0 ) ∈ F(R), such that: (1) for all h > 0 sufficiently small, there exist f (x 0 + h) − H f (x 0 ), f (x 0 ) − H f (x 0 − h) and the limits (in the metric D): for all h > 0 sufficiently small, there exist f ) and the limits (in the metric D): (h and −h in the denominators mean 1 h and − 1 h , respectively.) Note.In this paper, we only consider Cases ( 1) and ( 2) in the strongly generalized differentiability proposed by Bede and Gal [52].Chalco-Cano and Roman-Florés [54] stated that Cases (1) and ( 2) are more important, since Cases (3) and ( 4) occur only on a discrete set of points.(1) if f is (1)-differentiable, then f α (x) and f α (x) are differentiable functions and

Fuzzy Sumudu Transform
In order to establish results, some definitions are needed.G(u) and S[f (x)] will be used as the notations for the fuzzy Sumudu transform throughout this paper.Definition 10.Let f : R → F(R) be a continuous fuzzy-valued function.Suppose that f (ux) e −x is improper fuzzy Riemann-integrable on [0, ∞), then ∞ 0 f (ux) e −x dx is called the fuzzy Sumudu transform and is denoted by: where the variable u is used to factor the variable x in the argument of the fuzzy-valued function.
From Theorem 1, we obtain: From the classical Sumudu transform, we have: and: Finally, we have:

Duality Properties of the Fuzzy Laplace and Fuzzy Sumudu Transform
FLT has a close relationship with FST.It is necessary for us to be able to link between the two transforms in order to prove theorems and properties of the FST.The definition for FLT is given as follows.
Definition 11. [35] Let f (x) be a continuous fuzzy-valued function.Suppose that f (x) e −px is improper fuzzy Riemann-integrable on [0, ∞), then ∞ 0 f (x) e −px dx is called the fuzzy Laplace transform and is denoted by: Theorem 3. Let f (x) be a continuous fuzzy-valued function.If F is the fuzzy Laplace transform of f (x) and G is the fuzzy Sumudu transform of f (x), then: Substituting w = ux or x = w u , then we have: In the following corollary, we show that the roles of F and G in Equation ( 3) can be interchanged.
Corollary 1.Let f (x) ∈ F(R), having F and G for the fuzzy Laplace transform and fuzzy Sumudu transform, respectively.Then: Proof.The proof of Equation ( 4) can be obtained by changing u to 1 p in Equation (3).
Equations ( 3) and ( 4) form the fuzzy Laplace-Sumudu duality and serve as a mean of changing between those two transforms when needed.

Fundamental Theorems and Properties of the Fuzzy Sumudu Transform
In this section, we provide some theorems and properties associated with FST.Please note that the theorems and properties proposed in this section are extension of the classical Sumudu transform, as studied in [13,14].
Theorem 4. Let f, g : R → F(R) be two continuous fuzzy-valued functions.Suppose that c 1 and c 2 are arbitrary constants, then: Proof.Assume that f (x) = [f α (x), f α (x)] and g(x) = [g α (x), g α (x)].First, we proof for the lower bound of f (x) and g(x).
Secondly, we proof for the upper bound of f (x) and g(x).
Finally, we conclude that: The proof is complete.
In the following theorem, we provide the first preserving theorem.
Theorem 5. Let f : R → F(R) be a continuous fuzzy-valued function and a an arbitrary constant, then: Proof.From Definition 10, The theorem is proven to be true.
Next, we provide the second preserving theorem.Theorem 6.Let f : R → F(R) be a continuous fuzzy-valued function, then: Proof.From the definition of FST, Hence, Multiplying both sides with u, we obtain: For Case 2 from Theorem 2, equivalent to: u . Since: and: Since f is (2)-differentiable, it follows that: The proof is now complete.
In the next theorem, we provide the first shifting theorem.
Theorem 8. Let f : R → F(R) be a continuous fuzzy-valued function and a an arbitrary constant, then: Proof.From Definition 10, By using substitution w = (1 − au)x, we then obtain: The convolution theorem is provided below.
Theorem 9. Let f, g : R → F(R) be two continuous fuzzy-valued functions.Let F (p) and G(p) be fuzzy Laplace transforms, and let M (u) and N (u) be fuzzy Sumudu transforms for f and g, respectively.
Then, the Sumudu transform of the convolution of f and g, Proof.The FLT for (f * g) as in [55] is given by: By the fuzzy Laplace-Sumudu duality relation, and since: the FST for (f * g) is as follows: = uM (u)N (u).

Procedure for Solving Fuzzy Differential Equations
We consider a crisp differential equation given by: where f : [t 0 , T ] × R → R. Suppose that the initial value in Equation ( 5) is not precisely known and modeled with a fuzzy number, we have the following fuzzy initial value problem [56]: where f : [t 0 , T ] × F(R) → F(R) is a continuous fuzzy mapping.By referring to Kaleva [57], we observe that Theorem 2 provides a procedure to solve the Equation ( 6).As a matter of fact, By using FST on Equation (6), we have the following equation.
Case 1: If we consider Y (t) by using a (1)-differentiable function, then from Theorem 2, we get Y (t) = [y α (t), y α (t)].Now, we obtain the following FDE to be solved.
From Theorem 7, for Case 1, To solve Equation ( 9), first we assume that: where L 1 α (u) and U 1 α (u) are solutions of Equation (9).By using the inverse of FST, we then compute y α (t) and y α (t) as follows.
Case 2: If we consider Y (t) by using a (2)-differentiable function, then from Theorem 2, we get . Now, we obtain the following FDE to be solved.
From Theorem 7, for Case 2, To solve Equation ( 11), first we assume that: where L 2 α (u) and U 2 α (u) are solutions of Equation (11).By using the inverse of the FST, we then compute y α (t) and y α (t) as follows.

A Numerical Example
In this section, we provide a numerical example of solving an FDE using the FST.
Example 1.Consider the following initial value problem: By using FST, we have: and: First, we consider the condition where Y (t) is (1)-differentiable.Therefore, from Theorem 2, From Equation (9), we have: u , and: Then, we obtain and: Thus, and:  The results obtained using the FST for both cases proposed in this paper are shown in Figures 1 and 2, respectively.We can see that for Case 1, the result diverges as t increases; while for Case 2, the result implicates that the solution converges as t increases.

Remark.
From the example, we notice that the solutions depend on the differential equation we chose.For Case 1, the solution has the property that the diameter, i.e. diam(supp y(t)) = 2ae t , which is unbounded as t approaches infinity.Comparing to Case 2, the diam(supp y(t)) = 2ae −t → 0 as t approaches infinity, which leads to much more intuitive results.

Conclusions
In this paper, we have studied the classical Sumudu transform in the fuzzy setting.We have also proposed detailed procedures to solve FDEs.In the last part, we have conducted a numerical example of solving first order linear FDE using the FST.