Acessibilidade / Reportar erro

Adapted Fuzzy Integral: An Application in the Finite Element Method

ABSTRACT

In this paper we study and define an adapted fuzzy integral, based on the Sugeno integral. Moreover, we present a numerical integration formula which approximates the value of the adapted fuzzy integral. Thus, we prove that the Riemann integral and the adapted fuzzy integral are equivalent for power functions. Next, we apply the formula proposed in the numerical integration, required in the finite element method, to obtain a numerical solution of a boundary value problem for the one-dimensional Poisson equation. Finally, we observed better results of the approximate solution obtained in the example with the use of our formula when compared with the simple trapezoidal rule.

Keywords:
Fuzzy Measure; Sugeno Integral; Finite Element Method; Boundary Value Problem

RESUMO

Nessa proposta estudamos e definimos uma integral fuzzy adaptada, baseada na integral de Sugeno. Ademais, apresentamos uma fórmula de integração numérica que aproxima o valor da integral fuzzy adaptada. Assim, provamos que a integral de Riemann e a integral fuzzy adaptada são equivalentes para funções potência. Logo, aplicamos a fórmula proposta na integração numérica, requerida no método de elementos finitos, para obter uma solução aproximada de um problema de valor de contorno para a equação de Poisson unidimensional. Finalmente, observamos melhores resultados na solução aproximada obtida com o uso da nossa fórmula quando comparada com a regra simples de trapézio.

Palavras-chave:
Medida Fuzzy; Integral de Sugeno; Método de Elementos finitos; Problema de Valor de Contorno

1 INTRODUCTION

The finite element method (FEM) is a general tool to obtain numerical solutions to differential equations. This method is used in various fields of knowledge, such as mathematics, physics and engineering, and in several applications modeled by boundary value problems (BVPs) 1111. D. Sánchez, L.T. Bassani, L.C. Barros & E. Esmi. Ensaios do Método de Elementos Finitos com Integral Fuzzy. In Proceedings of IV CBSF. Recentes Avanços em Sistemas Fuzzy, SBMAC (2016)..

The Sugeno integral was introduced in 1974 for functions whose co-domain is the interval [0,1], that is, it was created to deal with membership functions of fuzzy subsets 1212. M. Sugeno. Theory of fuzzy Integrals and Applications. Ph.D. thesis, Tokyo Institute of Technology (1974)..

However, analytical or numerical integration of functions, that may not represent fuzzy subsets, are required in FEM to solve a BVP. Here, we introduce a numerical approximation for an adapted fuzzy integration based on the Sugeno integral.

This paper is organized as follows. Section 2 introduces some concepts of fuzzy sets theory. In Section 3, we briefly present the FEM to solve a BVP (in a classical way). In Section 4, we define an adapted fuzzy integral based on the Sugeno integral and we propose a numerical integration formula for monotonic and differentiable functions whose range is [0,k], where k ∈ℝ+. In Section 5 we show the approximate resolution of a BVP using our approach and we compare these results with the ones obtained by employing the well-known trapezoidal rule. Finally, in Section 6, we conclude with some final remarks.

2 BASIC CONCEPTS OF FUZZY SETS THEORY

Definition 2.1 (Fuzzy subset). 33. L.C. Barros, R.C. Bassanezi & W.A. Lodwick. A First Course in Fuzzy Logic, Fuzzy Dynamical Systems, and Biomathematics. Springer, New York (2017).),(44. L.T. Gomes, L.C. Barros & B. Bede. Fuzzy Differential Equations in Various Approaches. SBMAC SpringerBriefs. Springer, New York (2015).), (1313. L.A. Zadeh. Information and control. Fuzzy sets, 8(3) (1965), 338-353.A fuzzy subset A of a (classical set or) universal set 𝕌 is characterized by its membership function

ϕ A : 𝕌 [ 0,1 ] . (2.1)

If the range of the function ϕ A is the set {0;1} in (2.1), then A is said to be a crisp subset.

Definition 2.2 (α-level). Let A be a fuzzy subset of 𝕌 and α ∈ [0;1]. The α-level of the fuzzy subset A is the classical subset of 𝕌 defined by

[ A ] α = { x 𝕌 : ϕ A ( x ) α } , for α [ 0,1 ] .

Unless otherwise stated, we focus on the fuzzy subsets A of ℝ, such that their α-levels are given by [A]α=[aα,aα+],aαaα+, or [A]α= for each α ∈ (0;1]. We use the symbol ℝα to denote the class of these fuzzy subsets.

Definition 2.3 (Fuzzy measure).11. G. Arenas-Diáz & E.R. Ramírez-Lamus. Medidas difusas e integrales difusas. Universitas Scientiarum, 18 (1)(1) (2013), 7-32.),(88. H. Nguyen & E. Walker. A First Course in Fuzzy Logic. CRC Press Taylor & Francis Group (2006).Let 𝒜 be a σ-algebra of a (classical) set Ω. A map µ : 𝒜 → [0; ∞) is called a fuzzy measure when it satisfies:

  1. µ(∅) = 0 and

  2. if A;B𝒜 and AB, then µ(A) ≤ µ (B).

The definition of fuzzy measure proposed by Sugeno in 1212. M. Sugeno. Theory of fuzzy Integrals and Applications. Ph.D. thesis, Tokyo Institute of Technology (1974). adds the boundary condition µ(Ω) = 1 in i), that is, a normalization of the fuzzy measure, which implies that µ is a function from 𝒜 to [0,1] 11. G. Arenas-Diáz & E.R. Ramírez-Lamus. Medidas difusas e integrales difusas. Universitas Scientiarum, 18 (1)(1) (2013), 7-32.. For this work, we denote μ˜ to the normalized fuzzy measures.

Definition 2.4 (Usual Lebesgue measure).Let A ∈ ℝα and α ∈ [0;1]. The usual Lebesgue measure µ of the α-level of A is given by

μ ( [ A ] α ) = a α + a α . (2.2)

The usual Lebesgue measure µ is a fuzzy measure 33. L.C. Barros, R.C. Bassanezi & W.A. Lodwick. A First Course in Fuzzy Logic, Fuzzy Dynamical Systems, and Biomathematics. Springer, New York (2017)..

3 GALERKIN FINITE ELEMENT METHOD

Here, let us consider a linear boundary value problem (BVP) with homogeneous Dirichlet boundary values, in x ∈ [0,1], that describes Poisson’s one-dimensional equation given by 1111. D. Sánchez, L.T. Bassani, L.C. Barros & E. Esmi. Ensaios do Método de Elementos Finitos com Integral Fuzzy. In Proceedings of IV CBSF. Recentes Avanços em Sistemas Fuzzy, SBMAC (2016).:

{ d 2 u d x 2 = f ( x ) ,0 < x < 1, u ( 0 ) = 0, u ( 1 ) = 0. (3.1)

Using the traditional way to find an approximate solution for this BVP by FEM, we multiply the differential equation by test functions vV 0, where V 0 is a Hilbert space that considers the boundary values and is defined by 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...

V 0 = { v C 1 [ 0,1 ] : 0 1 ( v ( x ) 2 + v ' ( x ) 2 ) d x < , v ( 0 ) = v ( 1 ) = 0 } ,

considering the weak derivative for v(x).

In this context, we can find uV 0 such that

0 1 u " ( x ) v ( x ) d x = 0 1 f ( x ) v ( x ) d x , v V 0 ,

and using the equivalent of Riemann’s integration by parts we obtain that

0 1 u " ( x ) v ( x ) d x = 0 1 u ' ( x ) v ' ( x ) d x [ u ' ( x ) v ( x ) ] 0 1 = 0 1 u ' ( x ) v ' ( x ) d x .

Thus, we obtain the variational formulation 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
for problem (3.1), which consists of determining uV 0 that satisfies

0 1 u ' ( x ) v ' ( x ) d x = 0 1 f ( x ) v ( x ) d x , v V 0 . (3.2)

Solving the BVP (3.1) is equivalent to solving a problem in a variational formulation, according to equation (3.2) 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
.

Since uV 0 (and V 0 is an infinite dimensional Hilbert space), we want to build an approximate solution u h on a finite dimensional subspace Vh0 of V h66. J. Li & Y. Chen. Computational Partial Differential Equations Using MATLAB. CRC Press Taylor & Francis Group (2008).. To this end, we consider a partition of the interval [0,1]:

τ h : 0 = x 0 < x 1 < < x n < x n + 1 = 1,

where h = x i -x i-1 , i =1,2,..., n+1, In this work, we use the finite dimensional vector space given by:

Vh0:= {vh: vh is a piecewise linear and continuous function on τ h ; v h (0) = v h (1) = 0}.

Moreover, we consider the hat-functions (see Figure 1) as the piecewise linear basis functions {φi}i=1n given by

φ i ( x ) = { x x i 1 h , x i 1 x x i , x i + 1 x h , x i x x i + 1 , 0, otherwise . (3.3)

Figure 1:
Piecewise linear basis functions of Vh0.

Using this finite element formulation (or this discrete variational formulation) 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
), (66. J. Li & Y. Chen. Computational Partial Differential Equations Using MATLAB. CRC Press Taylor & Francis Group (2008)., approximating the solution of (3.1) consists in obtaining uhVh0, such that

0 1 u h ' ( x ) v ' ( x ) d x = 0 1 f ( x ) v ( x ) d x , v V h 0 .

Using the basis functions of Vh0(φjVh0), u h is given by a linear combination of functions φj22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
),(66. J. Li & Y. Chen. Computational Partial Differential Equations Using MATLAB. CRC Press Taylor & Francis Group (2008)., with coefficients ξj such that

u h ( x ) j = 1 n ξ j φ j ( x ) and u h ' = j = 1 n ξ j φ j ' ( x ) . (3.4)

Since u h (x) in (3.4) is an approximation to u(x), we have

j = 1 n ξ j ( 0 1 φ j ' ( x ) v ' ( x ) d x ) = 0 1 f ( x ) v ( x ) d x , v V h 0 .

Taking v(x) = φi (x), for each i =1,2,..,n, we can find ξj , for j =1,2,...,n, by solving the following system of linear equations

j = 1 n ( 0 1 φ j 1 ( x ) φ i ' ( x ) d x ) ξ j = 0 1 f ( x ) φ i ( x ) d x , i = 1,2,..., n .

Thus, we have a discrete problem represented by a linear system of equations 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
),(66. J. Li & Y. Chen. Computational Partial Differential Equations Using MATLAB. CRC Press Taylor & Francis Group (2008)., that can be given in the matrix form as

A ξ = b , (3.5)

where A=[aij]i,j=1nn×n, with

a i j = 0 1 φ i ' ( x ) φ j ' ( x ) d x , (3.6)

ξ=[ξj]j=1nn and b=[bi]i=1nn, with

b i = 0 1 f ( x ) φ i ( x ) d x . (3.7)

If we assume the basis functions (3.3), then {φi'}i=1n are defined, for i =1,2,...,n, as follows:

φ i ' ( x ) = { 1 h , x i 1 x x i , 1 h , x i x x i + 1 , 0, otherwise .

Thus, the coefficients of matrix A, described above in (3.6), are given by

a i i = 0 1 φ i ' ( x ) φ i ' ( x ) d x = x i 1 x i ( 1 h 1 h ) d x + x i x i + 1 ( 1 h 1 h ) d x = 2 h , (3.8)

for i =1,2,...,n, and

a i , i 1 = 0 1 φ i ' ( x ) φ i 1 ' ( x ) d x = x i 1 x i ( 1 h 1 h ) d x = 1 h , (3.9)

for i =2,3,...,n.

Since A is symmetric, we have that a i-1 ,i = a i,i-1 and by the definition of the basis we conclude that a ij = 0 for |i - j| >1.

Thus, we may obtain vector b in (3.7), by using the simple trapezoidal rule to approximate a numerical integration:

b i = 0 1 f ( x ) φ i ( x ) d x = x i 1 x i f ( x ) φ i ( x ) d x + x i x i + 1 f ( x ) φ i ( x ) d x = x i 1 x i f ( x ) x x i 1 h d x + x i x i + 1 f ( x ) x i + 1 x h d x x i x i 1 2 f ( x i ) + x i + 1 x i 2 f ( x i ) = h f ( x i ) . (3.10)

Thus, if f(x i ) = f i for i =1,2,...,n, then the final configuration of system (3.5) is

A = 1 h [ 2 1 0 0 1 2 1 0 1 2 0 1 0 0 1 2 ] , ξ = [ ξ 1 ξ 2 ξ n 1 ξ n ] and b = h [ f 1 f 2 f n 1 f n ] .

Consequently, system (3.5) has an unique solution, because A is a sparse, tridiagonal, symmetric and positive definite square matrix, and therefore it has an inverse matrix.

Finally, coefficients ξ j are obtained and the approximate solution u h , which depends strongly on the the function f, is determined. The function u h corresponds to the numerical solution of a BVP (3.1) by the FEM 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
.

4 ADAPTATION OF THE SUGENO INTEGRAL

Definition 4.1 (Sugeno Integral). 99. W. Pedrycz & F. Gomide. Fuzzy systems engineering toward human-centric Computing. IEEE Press. John Wiley & Sons, New Jersey, EUA (2007).), (1212. M. Sugeno. Theory of fuzzy Integrals and Applications. Ph.D. thesis, Tokyo Institute of Technology (1974).Let f: Ω → [0;1] be a function andμ˜a normalized fuzzy measure on Ω. The Sugeno integral of f with respect toμ˜is given by

Ω f d μ ˜ = sup α [ 0,1 ] [ α H ˜ ( α ) ] , (4.1)

whereis the minimum operator andH˜(α)=μ˜({xΩ:f(x)α})is called the level function of f33. L.C. Barros, R.C. Bassanezi & W.A. Lodwick. A First Course in Fuzzy Logic, Fuzzy Dynamical Systems, and Biomathematics. Springer, New York (2017)..

Theorem 1.33. L.C. Barros, R.C. Bassanezi & W.A. Lodwick. A First Course in Fuzzy Logic, Fuzzy Dynamical Systems, and Biomathematics. Springer, New York (2017).), (1212. M. Sugeno. Theory of fuzzy Integrals and Applications. Ph.D. thesis, Tokyo Institute of Technology (1974).Let f: Ω → [0,1] be a function (typically a membership function) andμ˜a normalized fuzzy measure on Ω: If

H ˜ ( α ) = μ ˜ ( { x Ω : f ( x ) α } )

has a fixed pointα¯, then

Ω f d μ ˜ = α ¯ = H ˜ ( α ¯ ) .

In 11. G. Arenas-Diáz & E.R. Ramírez-Lamus. Medidas difusas e integrales difusas. Universitas Scientiarum, 18 (1)(1) (2013), 7-32., a theoretical and applied study in fuzzy measures and integrals is presented. It also states that the Sugeno integral can be interpreted geometrically as the side of the greatest inscribed square between the integrated function and the x-axis (see Figure 2b).

Figure 2:
Sugeno integral: (2a) The fixed point of H˜(α) and (2b) Geometric interpretation of (4.1).

For example, we can consider the function f: Ω →[0,1] defined by f(x) = - 4x 2 + 4x. Thus, f is a typical membership function of a fuzzy subset F of ℝ whose α-levels are given by

[ F ] α = { x : 4 x 2 + 4 x α } = [ 1 1 α 2 , 1 + 1 α 2 ] .

If µ is the usual Lebesgue measure (2.2) on Ω = ℝ, then the level function H˜(α) is

H ˜ ( α ) = μ ( [ F ] α ) = 1 + 1 α 2 1 1 α 2 = 1 α .

Therefore, according to (4.1), the Sugeno integral can be written as:

Ω f d μ ˜ = sup α [ 0,1 ] [ α 1 α ] .

Since H˜(α) is a decreasing function, we have

Ω f d μ ˜ = α ¯ = H ˜ ( α ¯ ) = 1 + 5 2 ,

since α¯=1α¯α¯=1+520.61803.

The function f, the level function H˜(α) and the result obtained by the Sugeno integral are shown in Figure 2a and 2b.

As mentioned in 77. T. Murofushi & M. Sugeno. Fuzzy measures and fuzzy integrals. Fuzzy Measures and Integrals: Theory and Applications, (2000), 3-41., the Sugeno integral was defined only for functions whose range is contained in [0,1] and for normalized fuzzy measures. As we can note in the previous example, the Sugeno integral is considered a good approximation to the Riemann integral (because 01fdx=2/30.66667). However, in 88. H. Nguyen & E. Walker. A First Course in Fuzzy Logic. CRC Press Taylor & Francis Group (2006)., it was mentioned that the use of the Sugeno integral for functions whose range is contained in [0,k], k > 0. But, when k ≠ 1 or the fuzzy measure µ is not normalized, the Sugeno integral may not be a good approximation to the Riemman integral (see examples in 1010. H. Román-Flores, A. Flores-Franulic & Y. Chalco-Cano. The fuzzy integral for monotone functions. Applied Mathematics and Computation, 185(1) (2007), 492-498.).

We established an adapted fuzzy integral based on the Sugeno integral and with respect to a finite fuzzy measure. This integral approximates the Riemann integral for any functions whose range is in [0,k], with k ∈ ℝ+.

Definition 4.2 (Adapted fuzzy integral). Let f : Ω [0;k] be a measurable function such that k = supx∈ Ω f (x). The adapted fuzzy integral of f on Ω, with respect to the finite fuzzy measure µ : 𝒜 [0; ∞), that is, µ) < ∞, is given by

^ Ω f d μ = { k μ ( Ω ) Ω f ^ d μ ˜ , i f k > 0 0 , i f k = 0 (4.2)

wheref^(x)=f(x)kfor all x ∈ Ω, andμ˜(B)=μ(B)μ(Ω)for all B ∈ 𝒜.

Remark 1.If µ is a normalized fuzzy measure and f a function whose range is [0;1], we have that µ(Ω) = 1 and k = 1. In these conditions, the adapted fuzzy integral (4.2) coincides with the (classical) Sugeno integral (4.1).

Next, we establish another representation of (4.2) and we present a numerical integration formula to calculate it.

First, from (4.2) and (4.1), we have that

Ω f ^ d μ ˜ = sup α [ 0,1 ] [ α H ˜ ^ ( α ) ] , (4.3)

where

H ˜ ^ ( α ) = μ ˜ { x Ω : f ^ ( x ) α } = μ ˜ { x Ω : f ( x ) k α } = μ { x Ω : f ( x ) k α } μ ( Ω ) = H ( k α ) μ ( Ω ) (4.4)

and H(β) = µ ({x∈Ω: f(x) ≥ β }).

Thus, according to Theorem 1 and Equation (4.4), we have that the value of (4.3) is given by

Ω f ^ d μ ˜ = α ¯ = H ˜ ^ ( α ¯ ) = H ( k α ¯ ) μ ( Ω ) .

Taking β¯=kα¯, we have that

α ¯ = β ¯ k = H ( β ¯ ) μ ( Ω ) . (4.5)

Therefore, another representation of (4.2), for k > 0, is given as

^ Ω f d μ = k μ ( Ω ) α ¯ = μ ( Ω ) β ¯ = k H ( β ¯ ) .

Thus, from Equation (4.5), the calculus of the adapted fuzzy integral boils down to solving the equation β=kμ(Ω)H(β).

In what follows we provide a numerical integration formula to approximate the Riemann integral, over [a, b] ⊆ ℝ+ of differentiable and increasing functions.

Let f: Ω + be a differentiable and increasing function. The Riemann integral of f over Ω = [a,b] is approximated by

S Ω ( f ) = d f ( a ) + d c 2 c + d f ' ( b ) , (4.6)

where c = f(b) - f(a) and d = b - a = µ (Ω).

Note that

a b f ( x ) d x = a b f ( a ) d x + a b ( f ( x ) f ( a ) ) d x ( b a ) f ( a ) + ^ Ω g d μ = d f ( a ) + c μ ( Ω ) Ω g ^ d μ ˜ . d f ( a ) + d c 2 c + d f ' ( b ) = S Ω ( f ) .

The value ^Ωgdμ=dβ¯dc2c+df'(b) is obtained as follows. From (4.5), we have

β ¯ = c H ( β ¯ ) μ ( Ω ) = c d μ ( [ g 1 ( β ¯ ) , b ] ) = c d ( b g 1 ( β ¯ ) )

From g(g1(β¯))=g(bdβ¯c) and using first-order Taylor series about b, we obtain β¯g(b)dβ¯cg'(b) and, therefore,

β ¯ g ( b ) 1 + d g ' ( b ) c = f ( b ) f ( a ) c + d f ' ( b ) c = c 2 c + d f ' ( b ) .

The next theorem establishes that formula (4.6) coincides with the Riemann integral of certain power functions.

Theorem 2. Let f : Ω + be a power function given by f (x) = c 1 x n +c 2 with c 1 , c 2 , n ∈ ℝ+ , and let µ be an usual Lebesgue measure in Ω = [0;b]. We have

0 b f ( x ) d x = S Ω ( f ) .

Proof. On the one hand, we have that

0 b f ( x ) d x = 0 b ( c 1 x n + c 2 ) d x = ( c 1 x n + 1 n + 1 + c 2 x ) | 0 b = c 1 b n + 1 n + 1 + c 2 b .

On the other hand, from (4.6) we get

S Ω ( f ) = b f ( 0 ) + b ( c 1 b n ) 2 c 1 b n + b n c 1 b n 1 = b c 2 + b ( c 1 b n ) 2 c 1 b n ( 1 + n ) = b c 2 + c 1 b n + 1 1 + n .

A consequence of Theorem 2 for Ω = [a, b] ⊆ ℝ+ is that

a b ( c 1 x n + c 2 ) d x = 0 b ( c 1 x n + c 2 ) d x 0 a ( c 1 x n + c 2 ) d x = S [ 0, b ] ( c 1 x n + c 2 ) S [ 0, a ] ( c 1 x n + c 2 ) = b f ( 0 ) + b ( f ( b ) f ( 0 ) ) 2 ( f ( b ) f ( 0 ) ) + b f ' ( b ) a f ( 0 ) a ( f ( a ) f ( 0 ) ) 2 ( f ( a ) f ( 0 ) ) + a f ' ( a ) .

5 NUMERICAL APPROACH AND SIMULATIONS

In the following example, we use the numerical integration formula (4.6) to calculate all the required elements for the matrix system (3.5) and, thus, to approximate the solution of the BVP (3.1) via FEM using an adapted fuzzy integral (4.2) based on the Sugeno integral.

First, in the elements a i,i of the matrix A, with Ωi = [x i-1 ,x i ], and using equidistant nodes d = h and c = 0, we have

a i i = x i 1 x i 1 h 2 d x + x i x i + 1 1 h 2 d x = S Ω i ( 1 h 2 ) + S Ω i + 1 ( 1 h 2 ) = 1 h + 1 h = 2 h . (5.1)

Thus, we observe that the entries a ii of the matrix A given in (3.8) coincide with those obtained by our proposal in (5.1).

Second, using the adapted fuzzy integral to determine the entries a i,i-1 = a i-1,i of matrix A, we have d = h, c = 0, and

a i , i 1 = x i 1 x i 1 h 2 d x = S Ω i ( 1 h 2 ) = 1 h . (5.2)

One can observe that the coefficients a i,i-1 = a i-1,i obtained via the classical integration in (3.9) and the our formula (5.2) are equal. By definition of the basis, we have that a ij = 0 for |i - j| >1 since SΩi(0)=0.

Finally, the load vector b is obtained using the numerical integration formula (4.6) for increasing functions f(x) and g(x) = xf(x), for I = 1,2,...,n:

b i = x i 1 x i f ( x ) ( x x i 1 h ) d x + x i x i + 1 f ( x ) ( x i + 1 x h ) d x = 1 h [ x i 1 x i x f ( x ) d x ( x i 1 ) x i 1 x i f ( x ) d x + ( x i + 1 ) x i x i + 1 f ( x ) d x x i x i + 1 x f ( x ) d x ] 1 h [ S Ω i ( g ) ( x i 1 ) S Ω i ( f ) + ( x i + 1 ) S Ω i + 1 ( f ) S Ω i + 1 ( g ) ] .

In order to illustrate the good results of formula (4.6), we solve the BVP in (3.1) taking the load function f(x)=15x2x. This yields the exact solution which satisfies the boundary values given by

u ( x ) = ( 60 63 ) ( x x 9 / 2 ) . (5.3)

In Figure 3a and Figure 3b, we can see the analytical solution (5.3) as well as the numerical solutions obtained by numerical integration based on the trapezoidal rule (3.10) and by applying the formula (4.6). In this case, for i =1,2,...,n, we consider x i = ih to divide the domain [0,1] for both cases.

Figure 3:
Graphs of the analytic solution for the BVP in (3.1) with f(x)=15x2x, and the numerical solutions obtained by FEM using the trapezoidal rule and S Ω(f).

Table 1 shows the errors E(x i ) = U aprox(x i ) - u(x i ), for different values of n, between the numerical solutions U aprox(x i ) and the exact solution u(x i ). Here, we consider the infinity norm and the 2-norm 55. R.J. Leveque. Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems. SIAM (2007).:

E = max 1 i n | E ( x i ) | , and E 2 = i = 1 n | E ( x i ) | 2 .

Table 1:
Errors between the numerical solutions and the exact solution of BVP in (3.1), with f(x)=15x2x

6 FINAL CONSIDERATIONS

In this paper, we introduced the notion of an adapted fuzzy integral based on the Sugeno integral. In addition, we propose a numerical integration formula (4.6) for differentiable, increasing, and non-negative functions whose range is contained in [0,k]. We verified that, under certain conditions, our proposal yields error equal to zero for power functions, according to Theorem 2. In initial simulations, the numerical solution of a BVP via FEM using our approach obtained better results than a well-known numerical integration method, namely trapezoidal rule (see Figure 3 and Table 1).

ACKNOWLEDGMENT

This work was partially supported by CONICYT of Chile, CNPq under grant n. 306546/2017-5 and n. 142380/2016-4, and Fapesp under grant n. 2016/26040-7.

REFERENCES

  • 1
    G. Arenas-Diáz & E.R. Ramírez-Lamus. Medidas difusas e integrales difusas. Universitas Scientiarum, 18 (1)(1) (2013), 7-32.
  • 2
    M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf
    » http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf
  • 3
    L.C. Barros, R.C. Bassanezi & W.A. Lodwick. A First Course in Fuzzy Logic, Fuzzy Dynamical Systems, and Biomathematics. Springer, New York (2017).
  • 4
    L.T. Gomes, L.C. Barros & B. Bede. Fuzzy Differential Equations in Various Approaches. SBMAC SpringerBriefs. Springer, New York (2015).
  • 5
    R.J. Leveque. Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems. SIAM (2007).
  • 6
    J. Li & Y. Chen. Computational Partial Differential Equations Using MATLAB. CRC Press Taylor & Francis Group (2008).
  • 7
    T. Murofushi & M. Sugeno. Fuzzy measures and fuzzy integrals. Fuzzy Measures and Integrals: Theory and Applications, (2000), 3-41.
  • 8
    H. Nguyen & E. Walker. A First Course in Fuzzy Logic. CRC Press Taylor & Francis Group (2006).
  • 9
    W. Pedrycz & F. Gomide. Fuzzy systems engineering toward human-centric Computing. IEEE Press. John Wiley & Sons, New Jersey, EUA (2007).
  • 10
    H. Román-Flores, A. Flores-Franulic & Y. Chalco-Cano. The fuzzy integral for monotone functions. Applied Mathematics and Computation, 185(1) (2007), 492-498.
  • 11
    D. Sánchez, L.T. Bassani, L.C. Barros & E. Esmi. Ensaios do Método de Elementos Finitos com Integral Fuzzy. In Proceedings of IV CBSF. Recentes Avanços em Sistemas Fuzzy, SBMAC (2016).
  • 12
    M. Sugeno. Theory of fuzzy Integrals and Applications. Ph.D. thesis, Tokyo Institute of Technology (1974).
  • 13
    L.A. Zadeh. Information and control. Fuzzy sets, 8(3) (1965), 338-353.

APPENDIX

Table 2 presents some numerical results using the numerical integration S Ω( f ) (formula (4.6) based on the adapted fuzzy integral) for differentiable and increasing functions in the interval Ω = [0;1]. The results are compared with others numerical integration rules such as midpoint I MP (f), trapezoidal I Tra (f) and Simpson I Sim (f). The reference value was obtained using the Gauss Quadrature rule I GQ ( f ) ≈ abf(x)dx, which is an optimal numerical approximation which requires a resolution of a linear system of equations for each integration 22. M. Asadzadeh. An Introduction to the Finite Element Method (FEM) for Differential Equations (2012). URL http://www.math.chalmers.se/~mohammad/teaching/PDEbok/draft_ FEM_version4.pdf.
http://www.math.chalmers.se/~mohammad/te...
),(55. R.J. Leveque. Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems. SIAM (2007)..

Table 2:
Absolute errors obtained using several rules of numerical integration, and by (4.6).

Publication Dates

  • Publication in this collection
    Jan-Apr 2018

History

  • Received
    28 Mar 2017
  • Accepted
    22 Jan 2018
Sociedade Brasileira de Matemática Aplicada e Computacional Rua Maestro João Seppe, nº. 900, 16º. andar - Sala 163 , 13561-120 São Carlos - SP, Tel. / Fax: (55 16) 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br