Abstract

This work introduces a computational method for solving the linear two-dimensional fuzzy Fredholm integral equation of the second form (2D-FFIE-2) based on triangular basis functions. We have used the parametric form of fuzzy functions and transformed a 2D-FFIE-2 with three variables in crisp case to a linear Fredholm integral equation of the second kind. First, a method based on the use of two m-sets of orthogonal functions of triangular form is implemented on the integral equation under study to be changed to coupled algebraic equation system. In order to solve these two schemes, a finite iterative algorithm is then applied to evaluate the coefficients that provided the approximate solution of the integral problems. Three examples are given to clarify the efficiency and accuracy of the method. The obtained numerical results are compared with other direct and exact solutions.

1. Introduction

Several methods have been developed to estimate the solution of integral equation systems [13]. Many simple functions are used to approximate the solution of integral equations, such as orthogonal bases dependent on wavelets [4]. In addition, Maleknejad and Mirzaee developed the rationalized Haar functions [5] to approximate the solutions of the Fredholm linear integral equation method. In addition, second-type Fredholm integral equations are solved using direct triangular functions method as seen in [6] and using iterative algorithm-hybrid triangular functions method presented by Ramadan and Ali [7] where this hybrid method treats Fredholm integral equation of one dimension. More recently, Ramadan et al. [8] implemented such hybrid method to tackle system of two linear Fredholm integral equations of one dimension.

Furthermore, Maleknejad et al. [9] suggested by block pulse functions a numerical solution of the integral second-type equation.

It is explained using a series of orthogonal triangular functions, derived from the series of block pulses. Nevertheless, the fuzzy integral equations (FIEs) are required to solve and research a wide number of problems in various applied mathematics subjects, such as connection to physics, spatial, medical, and biology. FIEs therefore require approximate numerical solutions, as they are typically difficult to analytically solve. This thesis introduces a methodology used by the triangular functions (TFs) to solve the fuzzy linear FIE method of the second kind. In various implementation problems, certain parameters are typically represented by a fuzzy number rather than a crisp state, which involves the creation of mathematical models and computational algorithms to handle and solve the general fuzzy integral equations. A general method for solving the fuzzy Fredholm second-type integral equation is proposed in [10]. Recently, numerical methods have been developed to solve linear fuzzy Fredholm integral equation of the second kind in one-dimensional space (1D-FFIE-2) and two-dimensional space (2D-FFIE-2). Also, Fredholm fuzzy integral equations of the second kind are solved using the triangular functions [11], and numerical solution of linear Fredholm fuzzy equation of the second kind by block pulse functions is considered in [12]. Barkhordary et al. and Ramadan et al. [13, 14] presented a numerical technique for solving the fuzzy Fredholm integral equation of second kind. Numerical solution of two-dimensional fuzzy Fredholm integral equations of the second kind is presented via direct method using triangular functions [15]. Nouriani et al. [16] proposed a quadrature iterative method for solving the two-dimensional fuzzy Fredholm integral equations. Ezzati and Ziari [17], Hengamian Asl and Saberi-Nadjafi [18], and Bica and Popescu [19] illustrated a solution of the two-dimensional fuzzy Fredholm integral equations. A modified homotopy perturbation method for solving the two-dimensional fuzzy Fredholm integral equation is detailed in [20]. A two-dimensional nonlinear Volterra–Fredholm fuzzy integral equation is solved by using the Adomian decomposition method [21] and fuzzy bivariate triangular functions [22].

The aim of paper is to generalize the work proposed in [7] and [8] of these basis orthogonal triangular functions on (0, 1) to solve two-dimensional fuzzy Fredholm integral equations.

Section 2 presents some definitions and properties of the orthogonal triangular functions (TFs) (1D-TFs and 2D-TFs). Also, it expands functions by TFs. In Section 3, the definitions and properties of fuzzy function are given while a finite iterative algorithm is presented to solve coupled system of matrix equations in Section 4. The two-dimensional fuzzy integral equation is demonstrated and explained in Section 5 while the suggested method and the proposed iterative algorithm are detailed in Section 6. The illustrative examples and numerical results obtained are presented and discussed in Section 8.

2. Review of Triangular Functions (TFs)

2.1. Triangular Functions (TFs) of One Dimension

Definition 1. Two m-sets of triangular functions (TFs) are defined over the interval [0, T) [5]:where has a positive integer value; ; is the ith left-handed triangular function; and is the ith right-handed triangular function. Assuming , the TFs are defined over and . Based on this definition, it is clear that TFs are disjoint, orthogonal, and complete [5]. Therefore, one may writeThe first terms in the left-hand triangular functions and in the right-hand triangular functions can be written concisely in -vectors format aswhere and are called left-handed triangular function (LHTF) vector and right-handed triangular function (RHTF) vector, respectively. The product of two TF vectors yields the following properties:where 0 is the zero matrix. Also,in which is an identity matrix.

2.2. Two-Dimensional Triangular Functions and Their Properties [15]

An -set of 2D-TFs on the region is defined bywhere: ; ; ; and and are arbitrary positive integers. Therefore,

Furthermore,where is the th block pulse function defined on and as

Each of the sets and is obviously disjoint:

For , and .

Also, the 2D-TFs are orthogonal, that is,where denotes the Kronecker delta function and

On the other hand, ifthen , the 2D-TF vector, can be defined as

These relations are also satisfied for , similarly. Hence,

Finally by the orthogonality of , we havewhere denotes the Kronecker product defined for two arbitrary matrices P and Q as

The same equations are implied for , , and , by similar computations. Hence, we can carry out double integration of :where is matrix as follows:where and .

2.3. Function Expansion with 1D-TFs and 2D-TFs

The expansion of functions using triangular functions occurs in four situations.(1)The expansion of function over with respect to 1D-TFs is compactly written aswhere we may put and for i= 0, 1, …, m − 1.(2)The expansion of the function defined over Ω by 2D-TFs is as follows:where F is a vector given byand is defined in equation (21). The 2D-TF coefficients in can be computed by sampling the function at grid points and such that and , for various and . So, we havewhere and , . The vector F is called the 2D-TF coefficient vector.(3)The expansion of the function of three variables on with respect to 2D-TFs and 1D-TFs is as follows:where and are 2D-TF vector and 1D-TF vector of dimension and , respectively, and is a 2D-TF coefficient matrix. This matrix can be represented aswhere each block of F is an -matrix that can be computed by sampling the function at grid points such thatLet ; then,(4)The expansion of the function of four variables on with respect to 2D-TFs is as follows:where and are 2D-TF vectors of dimension and respectively, and is a 2D-TF coefficient matrix. This matrix can be represented aswhere each block of is an matrix that can be computed by sampling the function at grid points such thatLet and ; then,

In this paper, we suppose that for convergence.

3. Fuzzy Functions

We now remember through the paper some definitions that are required.

Definition 2. A fuzzy number is a fuzzy set u: R1 ⟶ [0, 1] that conforms to the following condition [23]:(a) is upper semicontinuous(b) outside some interval [c, d](c)There are real numbers a and b, c ≤ a ≤ b ≤ d, for which(i) is increasing in monotonic manner on [c, a](ii) is decreasing in monotonic manner on [b, d](iii) for a ≤ x ≤ b

Definition 3. A fuzzy number u is a pair of functions and , 0 ≤ r ≤ 1, satisfying the following requirement [5]:(a) is bounded monotonic increasing left continuous function(b) is bounded monotonic decreasing left continuous function(c), 0 ≤ r ≤ 1For arbitrary , , and k > 0, we define addition and multiplication by k as

4. Solving Coupled System of Matrix Equations Using Finite Iterative Algorithm [5]

Matrix equations can be solved using various forms of the finite iterative algorithms example [13, 5]. We consider iterative solutions to coupled system similar to the forms of Sylvester matrix equations [5].and second algorithm to solve coupled system of Sylvester matrix equations:

Algorithm 1 (see [5]). A finite iterative algorithm is developed to solve equation (35) as follows:(1)Input A, B, C.(2)Pick arbitrary matrices and .(3)Set(4)If , then stop and and are the final solutions; else, let and go to step 5.(5)Calculate

Algorithm 2 (see [5]). The following finite iterative algorithm is proposed to solve coupled system of Sylvester matrix equation (36):(1)Input matrices: A1; B1; A2; B2; C1; C2.(2)Pick arbitrary matrices and .(3)Set(4)If , then stop and and are the solutions; else set and then go to step 5.(5)Calculate

5. Two-Dimensional Fuzzy Fredholm Integral Equation

Two-dimensional FIE of the second kind is defined as follows [24]:

The linear (2D-FFIE-2) is defined aswhere and are fuzzy real functions on , is an arbitrary kernel function over , and is unknown on .

Throughout this paper, we consider 2D-FFIE-2 with , , and .

Now, introduce parametric form of a 2D-FFIE-2 with respect to Definition 3. Let and , be parametric form of and , respectively. Then parametric form of 2D-FFIE-2 is as follows:for each and . We can see that equations (43) and (44) are system of Fredholm integral equation of the second kind with three variables in crisp case.

6. Proposed Hybrid Iterative Technique

6.1. Converting Linear Two-Dimensional FIEs of Second Kind to Two Crisp Coupled Systems

This section presents an efficient method for soling a 2D-FFIE-2 by using 2D-TFs.

First, consider the following equation:

Now, the problem is to find the TF coefficients of from the known functions and kernel . 2D-TFs are applied for equations

To describe the approach of equation (47), first expand , and by 2D-TFs as follows:where and are defined in equations (3) and (21), respectively, and F are matrix of 2D-TF coefficients of and , respectively, and is -matrix 2D-TF coefficients of .

To obtain the solution of equation (47) from equations (49), (51), and (53), we have

Using equation (22), we haveand thenwhere and are -matrix and is - matrix, so is -matrix, where is unknown.

Then, we havewith

Setwhich lead to the following two crisp linear systems:and

Similarly, we expand and by 2D-TFs, and by substituting them into equation (42), two coupled crisp linear systems, similar to (61) and (65), are obtained. It is clear that all matrices in the two coupled crisp linear systems (61) and (65) are square matrices of dimensions . Thus, we need to obtain the coefficient matrices , and in order to get the approximate numerical solution of the form:

6.2. Proposed Iterative Algorithm for Solving Coupled Systems (61) and (65)

An iterative algorithm is proposed here to solve the two coupled systems (61) and (65) as a generalization of Algorithm 2.

Algorithm 3. Algorithm 2 is modified and generalized to work out for systems (61) and (65) as follows.
First, for coupled system (61):(1)Input .(2)Choose arbitrary matrices .(3)For , set(4)If , then stop and are the solutions; else, let and go to step 5.(5)ComputeFor coupled system (65), the algorithm is repeated with replacing by and by where the block matrix is computed.
The approximate crisp numerical solution for equation (51) of the form is then obtained.
In a similar manner, the crisp numerical solution for equation (53) of the form can be obtained by carrying out the above-proposed algorithm for the other coupled crisp systems similar to (61) and (65). Finally, the solution for linear 2D-FFIE-2 by 2D-TFs is then given as

7. Convergence Analysis of the Proposed Method

In this section, we obtain error estimate for the numerical method proposed in previous section.

Theorem 1. The solution of the two-dimensional Fredholm fuzzy integral equations given by equation (1) by using 2D-TFs converges to exact solution if

Proof. Assume and represent the exact and approximate solutions of equation (1), respectively. Therefore,By using equation (1), we can writewhereAlso, we have , so as and since is bounded.
Thus,so the proof of the theorem is completed.

Remark 1. In our theoretical investigation for the proposed method, we take .

8. Numerical Results and Discussion

This section demonstrates the effectiveness and the accuracy of our proposed hybrid method, 2D-TFs and an iterative algorithm, on some examples. The solution of each example is obtained for different values of , and and is compared with the exact solution, the direct method, and the presented method when the tolerance criteria residual is >e−4 and >e−8.

Example 1. The 2D-FFIE-2 given in Mirzaee et al. [15] is considered asIn this case, the exact solution is given byThe number of iterations for solving the two coupled matrix equations to obtain the coefficient matrices taken by our proposed iterative algorithm is when the tolerance criteria residual is >e−4 which indicates that the hybrid proposed method is quite efficient and has good accuracy as seen from Tables 1 and 2.

Remark 2. The numerical results for the approximate solution using the direct method in Tables 1 and 2 are taken from Table 3 in [15], while the numerical results using the direct and the proposed iterative methods in Tables 3 and 4 are obtained using our own program written using MATLAB R2018b. Also, the number of iterations for solving the two coupled matrix equations to obtain the coefficient matrices taken by our proposed iterative algorithm is when the tolerance criteria residual is >e−4 which indicates that the hybrid proposed method is quite efficient and has good accuracy as seen from Tables 3 and 4.

Example 2. Consider the following made up linear two-dimensional fuzzy Fredholm integral equations (2D-FFIE-2):In this case, the exact solution is given by

Remark 3. The number of iterations for solving the two coupled matrix equations to obtain the coefficient matrices taken by our proposed iterative algorithm is when the tolerance criteria residual is >e−4 which indicates that the hybrid proposed method is quite efficient and has good accuracy as seen from Tables 5 and 6.

Example 3. Consider the following 2D-FFIE-2:In this case, the exact solution is given by

Remark 4. In Table 7, the number of iterations for this example by our proposed iterative algorithm is when the tolerance criteria residual is >e−4 which indicates that the hybrid proposed method is quite efficient. Moreover, we can see that our method has good accuracy which can be further improved by increasing the residual.

9. Conclusion

Fuzzy control applications and a large proportion of applied mathematical topics require the solution of the fuzzy integral equations. The paper introduced the 2D-TFs method for approximating the solution of linear 2D-FFIE-2, which is a hybrid of triangular functions and an iterative algorithm. The method is simple, efficient, and accurate and is based on converting the original equation into two crisp systems (2D-FFIE-2). The efficiency and simplicity of the proposed method are demonstrated via numerical examples with known exact solutions, and the results are given. Furthermore, the exceptional value of the proposed method is low cost of the equation setting with no need for any projection method or integration.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

The authors declare that the study was realized in collaboration with the same responsibility. All authors read and approved the manuscript.

Acknowledgments

The authors did not receive any funds or grants from any public or private institutions.