1. Introduction
The fuzzy set theory was pioneered by Lotfi Zadeh [
1]. Fuzzy sets are widely used in a variety of mathematical modeling, and fuzzy derivative was introduced in [
2,
3,
4,
5,
6,
7,
8,
9] and is an extension of ordinary derivatives. This provides the foundation of fuzzy differential equations, and uncertain systems are modeled using fuzzy differential equations.
Fuzzy differential equations are solved using many methods, such as differential or fuzzy inclusion [
10,
11,
12], the extension principle and characterization theorem [
13,
14,
15], and artificial neural networks [
16]. Numerical methods such as the differential transformation method [
17,
18], extension of numerical solution methods for ordinary differential equations (ODEs) [
19,
20,
21,
22,
23,
24], interactive derivatives [
25], multistep methods [
26], and Runge–Kutta (RK) methods [
27] are also applied to solving these equations. The metric properties of fuzzy functions are studied in [
28,
29]. The existence of solutions for second-order fuzzy differential equations is proved in [
21,
30]. As a result, higher-order fuzzy differential equations can be solved by using generalized fuzzy derivatives. A system of differential equations and second-order differential equations are solved in [
6,
31].
A first-order differential equation with an initial value is called an initial value problem, and this sort of equation has a unique solution. If the initial values are fuzzified, then the problem is called a fuzzy initial value problem (FIVP). The differential equation with fuzzy initial values has two solutions for every fuzzification parameter. These solutions depend upon the definition of fuzzy differentiability, and the definition changes as the sign of the product of function and its first derivative change [
32].
Fuzzy differential equations are interpreted as a family of differential inclusions. Fuzzy interval arithmetic is very simple to compute because it is applicable to the endpoints of an interval, and this leads to an estimate of the confidence interval. These fuzzy solutions are constructed from the solution of the crisp problem generated by the fuzzy initial value problem. Fuzzy solutions are developed using the parametric representation of fuzzy functions. In this representation, fuzzy sets are represented by lower and upper approximations: the lower approximation is increasing and the upper approximation is a decreasing function. The function’s property of increasing and decreasing is utilized with the definition of switching points. The fuzzy differentiability of type-(i) switches to that of type-(ii) at the point where the function changes its increasing nature to decreasing or vice versa. This property is also used for the first-order derivative function, and similar logic is implemented for developing the solution of second-order differential equations.
The differential transformation method (DTM) [
33] is a way of finding the numerical solution of linear and nonlinear differential equations [
34]. The DTM is also used for solving fuzzy differential equations (FDEs) with FIVPs [
10,
35]. In this paper, a solution of a second-order fuzzy differential equation is proposed. The second-order differential equation is converted to an equal system of first-order differential equations, and initial values are fuzzified. The switching points are verified for the solution and first-order derivative functions. The DTM is used to solve the system of equations, and the solution leads to the approximation of a fuzzy solution of the problem.
Two-point initial value problems are considered in [
14,
36,
37,
38], and a fuzzy boundary value problem is solved in [
39]. In these approaches, the definition of differentiability in the given interval is chosen initially without switching. In this paper, the FIVP is approached by using FDTM, and the idea is extended to second-order fuzzy differential equations with the adaptability of min and max operators to handle the switching points. The solution set produces the upper and lower approximations for the bounded solution of the second-order differential equation. The first-order FIVP represents a single solution with the application of two definitions of differentiability, and for second-order FIVP, the four definitions of differentiability are combined with the application of min, max operators, where the definition of differentiability depends upon the interval of the validity of the solution and the switching point within that interval. This leads to a unique solution of the differential equation. The rest of the paper is organized as follows: First, we start with preliminaries in
Section 2. The switching points and fuzzy differentiability for second-order derivatives are explained in
Section 3. In
Section 4, we describe the DTM, and its extension is derived in
Section 4.1. In
Section 5, numerical examples are discussed, and finally,
Section 7 consists of the conclusion with a summary of our contribution, as well as a discussion of future directions of investigation.
2. Preliminaries
Numerous techniques and terms are used to define a fuzzy number and fuzzy membership functions [
40]. In this paper, the following definitions are used.
Definition 1 (Fuzzy Set). An ordered pair of a set and a membership function is called a fuzzy set, and a membership function is defined as . The fuzzy set has the following properties:
- 1.
A set of all the points where the membership function has a nonzero value is called the support of the fuzzy set.
- 2.
The set of points where is called the core of a fuzzy set η.
- 3.
A fuzzy set η is normal if ∃ at least one such that .
- 4.
A fuzzy set η is convex if for all and .
- 5.
A fuzzy set is upper semi-continuous in R
- 6.
Compactness of a fuzzy set is defined as , where stands for a closure operator.
The fuzzy real numbers, numbers with a triangular sendograph, are represented by such that . These fuzzy number are represented as three points, , with ; then, the -level fuzzy set is represented by and . The cuts for all are defined as . These are called the level sets, and the zero-level compact set is represented by . The level sets are nested sets, and the cuts of any fuzzy number are always closed and bounded intervals.
Definition 2. Fuzzy numbers are represented as an ordered pair of functions , where , and and are called the lower and upper boundaries of . The is bounded, left continuous, and decreasing, and is bounded, right continuous, and an increasing function on , and for all . A crisp number is represented by an equivalent fuzzy notation with an equal approximation of lower and upper . Let h and g be two fuzzy numbers defined in γ-level notation, where and ; then, the algebra of fuzzy numbers is defined as
- Addition:
Two fuzzy numbers h and g are added as follows: and , and its γ-level addition will be
- Subtraction:
The subtraction of fuzzy triangular numbers is computed as follows: and . The subtraction of fuzzy numbers is defined as addition, and the parameterized subtraction is defined as
- Scalar Multiplication:
The scalar multiplication of a scalar a with a fuzzy number is computed as follows: when and when , and the scalar multiplication will become
- Multiplication:
Let be two fuzzy numbers; then, multiplication of fuzzy numbers is computed as follows: , where and
- Division:
Division between two fuzzy numbers is defined conditionally. If , then is defined as and
Definition 3 (H-Difference)
. The set of real numbers is defined by fuzzy membership function E, and are three fuzzy numbers in E such that . Then, the H-difference of two fuzzy real numbers is denoted as . Its γ-level representation will become Definition 4 (g-Difference)
. The set of real numbers is defined by fuzzy membership function E, and are three fuzzy numbers in E such that . Then, the g-difference of two fuzzy real numbers is defined as , and its γ cuts are Definition 5 (Hausdorff Distance)
. Let be a function. Then, the Hausdorff distance is defined aswhere a complete metric space is represented by D. The g-difference of intervals always exists. This definition is equivalent to the usual definitions for metric spaces of fuzzy numbers. Definition 6 (Generalized Hukuhara Difference)
. The generalized Hukuhara difference between two fuzzy numbers is defined as a and b as follows:where and with the conditions that is a nondecreasing function and is a nonincreasing function such that . Definition 7 (Generalized Hukuhara Differentiability [
41])
. Let f be a fuzzy function and , and let t be any point in the domain. Then, f is a strongly generalized differentiable function f at t and will be defined by the following.for allIf h is sufficiently small and , then
Differentiable type-(i): There exist and ; then,or Differentiable type-(ii): If and exist and
for allIf h is sufficiently small and , then
Differentiable type-(i): There exist and ; then, Differentiable type-(ii): There exist and ; then,
Reasoning as above, there exists the Hukuhara difference . Now, if we suppose , then we see that we cannot use the above kind of reasoning to prove that the H-differences and and the derivative exist.
Definition 8 (Second-Order Generalized Hukuhara Differentiability [
41])
. Let be a strongly differentiable function in the generalized sense on the domain, and let there exist such that, for all , we have- (I)
and and - (II)
and and
If f has differentiability of type-(i) and exists, then there are two possibilities for the differentiability of function . Similarly, if f has differentiability of type-(ii), then has two possibilities. Thus has four definitions of differentiability.
Definition 9 (Generalized Differentiability ([
42]))
. Let and k be such that ; then, the g_derivative of the function at is defined asIf exists, then f is called differentiable in the generalized sense. Theorem 1. Let be such thatwhere and are differentiable real-valued functions with respect to x. Then, is g_ differentiable, and mathematically,The interval is divided into subintervals. The functional value at the end point of the previous interval becomes the initial value of the next subinterval and at each boundary point of the subinterval. The conditions are verified, and if there is a switching point, then the definition of fuzzy derivative changes from type-(i) to type-(ii) and vice versa. The use of min and max operators handles the switching points, which are defined in the next subsection. (Proof of Theorem 1 follows from [42].) Proposition 1. ([43]) Let be Hukuhara differentiable 2 times andsuch that and . If are differentiable, then we can write in a similar fashion:In fact, the second-order derivative of a real-valued function represents the four definitions of fuzzy derivatives of fuzzy functions. The definition of the derivative changes at the switching point. Theorem 2 (Theorem 3.6 in [
41])
. Let be a function and denote for each . Then,If F is (I)-differentiable, then f and g are differentiable functions and
If f is (II)-differentiable, then f and g are differentiable functions and .
Theorem 3 (Theorem 3.9 in [
41])
. Let be an interval-valued function, where f is an increasing function and g is a decreasing function, and let the system and be two fuzzy functions. Then,- 1.
Let and be two differentiable functions. If F and are (I)-differentiable, then .
- 2.
Let and be two differentiable functions. If F is (I)-differentiable and is (II)-differentiable, then .
- 3.
Let and be two differentiable functions. If F is (II)-differentiable and is (I)-differentiable, then .
- 4.
Let and be two differentiable functions. If F and are (II)-differentiable, then .
Theorem 4 (Theorem 3.1 in [
44])
. Let and assume that is continuous. A mapping is a solution to the initial value problem and if and only if x and are continuous and satisfy one of the following conditions:, where and are (i)-differentials,
, where and are (ii)-differentials,
, where is (i)-differentiable and is a (ii)-differential, or
, where is (ii)-differentiable and is a (i)-differential.
Remark 1. (i) This theorem shows that four systems of fuzzy differentiability represent only two systems, and two pairs are equivalent. According to the conditions of the triangular fuzzy function, f is an increasing function and g is a decreasing function. If F follows the first definition, then . This shows that is an increasing function and is a decreasing function. In the other case, , which shows that is now decreasing and is an increasing function. (ii) Theorem 4 is for H-differentiability, and its mode-generalized form is differentiability, so it is also differentiable.
3. Fuzzy Differentiability and Switching Points
Definition 10 (Switching Point). The points in an interval where fuzzy differentiability of type-(i) changes to type-(ii) and vice versa are called switching points.
For the first-order initial value differential equations, the above-cited switching points apply. This logic can be derived from any one of the two functions
f and
g, and the definition of fuzzy differentiability is associated with any of the above functions. In
Table 1,
f is increasing. The definition changes from type-(i) to type-(ii) if
f changes its increasing nature to decreasing and vice versa.
The possibility for
f is described in the
Table 1. If both the function and its derivative are increasing functions (similarly for
g decreasing) then, type-(i) is applicable, if
changes its nature and is a decreasing function then differentiability of type-(ii) is applicable. This behavior of function is verified over the interval. We extend the same logic to the second derivative, and we get the representations listed in
Table 2.
The above table shows that the fuzzy function follows the same differentiability if there is no switch in increasing or decreasing nature of function or there are two switches. The pair with one switching point also have similar behavior. We transform the second-order differential equation into equal system of first-order differential equations.
This second-order differential equation is decomposed into the system of first-order differential equation as
where transformation
and
is used to describe the equivalent system of first-order differential equations.
The solution system described in [
36] will become
This system is used to initiate the fuzzy solution of the differential equations, and switching points are determined by using the fact that
f has differentiability type-(i) when
; otherwise, it has type-(ii). This shows that there is only one possible solution at one point, and the other solution is possible at another point. In between these two points, the definition of differentiability changes from type-(i) to type-(ii) and vice versa. Similarly, if
, then
has type-(i) differentiability [
41]. The point where
changes to
is called the switching point. This concept is similarly applied to second derivatives.
Definition 11. (Fuzzy Power Series ([
45]))
. Let u and be two numbers of the fuzzy number space , and their difference is a nonnegative fuzzy number. Let be a sequence of nonnegative fuzzy numbers. Now, define with and and . Then, the power series of fuzzy numbers with the coefficients is given by , which can be expressed in terms of γ-level sets as follows: Definition 12. The differential transform of the -order derivative of a function is defined asThen, , the inverse transform of the above relation, is defined asEquation (9) represents the Taylor series expansion of function at . 5. Numerical Examples
In this section, we consider the numerical examples of first and second order.
Example 1. Consider the following FIVP: Solution:
The analytical solution of this problem is . Before fuzzifying the differential equation, we first analyze the switching points of the solution function as mentioned in Section 3. For this purpose, the graph of the solution function and its first and second derivatives in Figure 1. The graphs show that and for . This shows that f has differentiability of type-(i), has differentiability of type-(ii), and the second-order differential equation has the solution of the system. This solution is switched for , where and , and the equation has the solution of type within the interval. This shows that the solution switches from type to .
Now, the second-order differential equation is converted into a system of first-order differential equations. Differential Equation (13) is converted into a system of first-order differential equations, represented in Equation (14). The initial conditions also change to . Replacing , then above equation will become Every real number can be represented as the triangular fuzzy number. The initial values are fuzzified, and this problem becomes an FIVP. These values are represented as the interval-valued fuzzy setSimilarly, the second value is fuzzified as The system of ODEs is converted into the equivalent system of ordinary differential equations with fuzzy initial value conditions. This process splits the system of ODEs into two equivalent systems of ODEs, and Equation (15) will become Taking the differential transformation of the above system of equations, we get the recurrence relation for this system of equations between the two grid points and , represented aswith initial conditions, we get and . The whole interval is divided into the subintervals, and switching points are checked at each partition point. The differential transformation coefficients are determined from the above set of equations, and the lower and upper solutions are developed for each subdomain separately using the differential transformation method of order 20, namely, . The final value of the preceding subdomain is used as the initial value of the next subdomain. For example, at , y(t) in the first subinterval will be such that . The results are compared in Table 4, in which the computed solutions for the grid point 0 are shown for the changing fuzzification parameter γ. The point is the switching point, where . For the next input, the values are interchanged, i.e., the upper solution becomes the lower solution and the lower solution changes to the upper solution by applying the min and max operators.
Example 2. Consider the following FIVP:over the interval . Solution:
The analytical solution of the problem is .
In the first step, the domain of the analytic solution function is analyzed, and the Figure 2 of the solution function with its first and second derivatives is drawn. This graph shows that for all . As marked in Section 3, if f has a solution according to type-(i) derivatives, and if for all as a result, then also has a derivative of type-(i). Thus, the problem has the solution in the system. The solution of a second-order differential equation is determined by converting Differential Equation (18) into a system of first-order ODEs according to Equation (19).The initial conditions are also transferred according to the used transformation, . Now, replacing ,The differential transformation of this ODE system at grid point is computed as Equation (21). To solve the second-order initial value problem, the system of initial value problem is converted to the fuzzy initial value problem. For this, the initial conditions are fuzzified as and , and γ-level transformation of these initial conditions will be and . Similarly, the recurrence relation for the upper solution will be written as stated in Equation (21). The computed solution for multiple grid points is shown in Table 5 for . The results show that there is no switching point, and the definition of the derivative is chosen according to the initial conditions.
Example 3. Consider the following FIVP: Solution:
is the analytical solution of Differential Equation (22). The Figure 3 shows that the function and its first and second derivatives are always negative. This results in for all . As a result, f has a solution according to the type-(i) derivative, and for all ; as a result, has a solution according to the type-(i) derivative, and the problem has the solution in the system. For a numerical solution of a higher-order differential equation with fuzzy initial values, first convert Equation (22) in the system of first-order differential equations as Equation (23).The initial conditions become . The system of initial value problem is converted to the fuzzy initial value problem, and these conditions become and . The differential equation for the system between the two grid points and can be represented as Equation (24), the recurrence relation for the differential transform method:Now, for the fuzzy initial value solutions, the lower system of the recurrence relation can be written as Equation (25). The transformation of initial conditions yields and . The series solution is computed for the upper and lower approximation in each subdomain separately. The initial values at the point are the terminal values of the previous series, which is A similar method is followed for other points, and the results are represented in Table 6. These results show that the solution has no switching point, and the system has differentiability in the domain.