Skip to main content

Fuzzy approximate solutions of second-order fuzzy linear boundary value problems

Abstract

In this paper, approximate solutions of second-order linear differential equations with fuzzy boundary conditions, in which coefficient functions maintain the sign, are investigated. The fuzzy linear boundary value problem is converted to a crisp function system of linear equations by the undetermined fuzzy coefficients method. The fuzzy approximate solution of the fuzzy linear differential equation is obtained by solving the crisp linear equations. Some numerical examples are given to illustrate the proposed method.

1 Introduction

Nowadays, fuzzy differential equations (FDEs) is a popular topic studied by many researchers since it is utilized widely for the purpose of modeling problems in science and engineering. Most of the practical problems require the solution of a fuzzy differential equation (FDE) which satisfies fuzzy initial or boundary conditions, therefore a fuzzy initial or boundary problem should be solved. However, many fuzzy initial or boundary value problems could not be solved exactly, sometimes it is even impossible to find their analytical solutions. Thus, considering their approximate solutions is becoming more important.

Prior to discussing fuzzy differential equations and their associated numerical algorithms, it is necessary to present an appropriate brief introduction to derivative of the fuzzy-valued function. The concept of a fuzzy derivative was first introduced by Chang and Zadeh [1], followed up by Dubois and Prade [2] who used the extension principle in their approach. Other fuzzy derivative concepts were proposed by Puri and Ralescu [3] and Goetschel and Vaxman [4] as an extension of the Hukuhara derivative of multivalued functions. Kandel and Byatt [5, 6] applied the concept of fuzzy differential equation to the analysis of fuzzy dynamical problems.

The numerical methods for solving fuzzy differential equation are introduced in [79]. In 2001, Buckley and Feuring [10] presented two analytical methods for solving an n th-order fuzzy linear differential equation with fuzzy initial conditions. Their first method of solution was to fuzzify the crisp solution and then check to see if it satisfies the fuzzy differential equations with fuzzy initial conditions. The second method was the reverse of the first method; in that they firstly solved the fuzzy initial value problem and then checked to see if it defined a fuzzy function. In 2008, Allahviranllo et al. [11] obtained the approximate solution of n th-order linear differential equations with fuzzy initial conditions by using the collocation method. In 2003, O’Regan et al. [12] proved a super-linear result for fuzzy boundary value problems relying on a general Schauder theorem in the metric space. Meanwhile Lakshmikantham et al. [13] investigated the solution of two-point boundary value problems associated with nonlinear fuzzy differential equations by using the extension principle. In 2008, Chen Minghao et al. [14] obtained the conclusion: two-point boundary value problems have the analytic solution only on condition that the new structure and properties to the fuzzy number are given. But for second-order fuzzy linear boundary value problems

y +p(t) y +q(t)y=g(t),t[a,b]
(1.1)

associated with

y(a)= α ˜ ,y(b)= β ˜ ,
(1.2)
y (a)= α ˜ , y (b)= β ˜ ,
(1.3)
y (a) α 0 y(a)= α ˜ , y (b)+ β 0 y(b)= β ˜ ,
(1.4)

it is not the case. Once the coefficient functions p(t), q(t), g(t) are continuous on [a,b] and p(t), q(t) maintain the sign on [a,b], a unique solution must exist.

In this paper, we consider the approximate solution of a class of second-order linear differential Eq. (1.1) under fuzzy boundary value conditions (1.2), (1.3) and (1.4). Based on the undetermined fuzzy coefficients method, we convert a second-order linear differential equation to the crisp system of linear equations. Secondly, we investigate their cases according to indifferent cases of coefficient functions p(t) and q(t) from the original systems and build the corresponding crisp systems of linear equations. Then we derive a fuzzy approximate solution of the fuzzy linear differential equation from solving crisp function systems of linear equations. Finally, some examples are given to illustrate the proposed method. The structure of this paper is organized as follows.

In Section 2, we recall some basic definitions and results about fuzzy numbers as well as fuzzy derivative of the fuzzy-valued function. In Sections 3, 4 and 5, we build crisp function systems of linear equations via analyzing different cases based on the coefficient functions of the fuzzy linear differential equation in detail. The proposed algorithms are illustrated by solving some examples in Section 6 and the conclusion is drawn in Section 7.

2 Preliminaries

2.1 Fuzzy numbers

Definition 2.1 [1]

A fuzzy number is a fuzzy set like u:RI=[0,1] which satisfies:

  1. (1)

    u is upper semi-continuous,

  2. (2)

    u(x)=0 outside some interval [c,d],

  3. (3)

    there are real numbers a, b such that cabd and

(3.1) u(x) is monotonic increasing on [c,a],

(3.2) u(x) is monotonic decreasing on [b,d],

(3.3) u(x)=1, axb.

Let E 1 be the set of all real fuzzy numbers which are normal, upper semi-continuous, convex and compactly supported fuzzy sets.

Definition 2.2 [2]

A fuzzy number u in a parametric form is a pair ( u ̲ , u ¯ ) of functions u ̲ (r), u ¯ (r), 0r1, which satisfies the following requirements:

  1. (1)

    u ̲ (r) is a bounded monotonic increasing left continuous function,

  2. (2)

    u ¯ (r) is a bounded monotonic decreasing left continuous function,

  3. (3)

    u ̲ (r) u ¯ (r), 0r1.

A crisp number x is simply represented by ( u ̲ (r), u ¯ (r))=(x,x), 0r1. By appropriate definitions, the fuzzy number space {( u ̲ (r), u ¯ (r))} becomes a convex cone E 1 which could be embedded isomorphically and isometrically into a Banach space [15, 16].

Definition 2.3 [2]

Let x=( x ̲ (r), x ¯ (r)), y=( y ̲ (r), y ¯ (r)) E 1 , 0r1 and arbitrary kR. Then

  1. (1)

    x=y iff x ̲ (r)= y ̲ (r) and x ¯ (r)= y ¯ (r),

  2. (2)

    x+y=( x ̲ (r)+ y ̲ (r), x ¯ (r)+ y ¯ (r)),

  3. (3)

    xy=( x ̲ (r) y ¯ (r), x ¯ (r) y ̲ (r)),

  4. (4)

    kx= { ( k x ̲ ( r ) , k x ¯ ( r ) ) , k 0 , ( k x ¯ ( r ) , k x ̲ ( r ) ) , k < 0 .

Definition 2.4 [17]

For arbitrary u=( u ̲ , u ¯ ),v=( v ̲ , v ¯ ) E 1 , the quantity

D(u,v)= ( 0 1 ( u ̲ ( r ) v ̲ ( r ) ) 2 d r + 0 1 ( u ¯ ( r ) v ¯ ( r ) ) 2 d r ) 1 2

is the distance between fuzzy numbers u and v.

2.2 Second-order fuzzy boundary value problem

Definition 2.5 [3]

Let x,y E 1 . If there exists z E 1 such that x=y+z, then z is called the Hukuhara-difference of fuzzy numbers x and y, and it is denoted by z=xy.

In this paper the sign stands always for Hukuhara-difference, and let us remark that xyx+(1)y.

Definition 2.6 [18]

Let f:[a,b] E 1 and t 0 [a,b]. We say that f is Hukuhara differential at t 0 , if there exists an element f ( t 0 ) E 1 such that for all h>0 sufficiently small, f( t 0 +h)f( t 0 ), f( t 0 )f( t 0 h) and the limits (in the metric D)

lim h 0 f ( t 0 + h ) f ( t 0 ) h = lim h 0 f ( t 0 ) f ( t 0 h ) h = f ( t 0 ).

Lemma 2.1 [19]

If g:[a,b]R is differential on [a,b] such that g , g are nonnegative and monotonic increasing on [a,b], then c E 1 , f(x)=cg(x) is differential on [a,b] and

f (x)=c g (x), f (x)=c g (x),x[a,b].

Definition 2.7 The second-order differential equation

y =f ( t , y , y ) ,t[a,b]
(2.1)

with the boundary value conditions

y(a)= α ˜ ,y(b)= β ˜ , α ˜ , β ˜ E 1 ,
(2.2)
y (a)= α ˜ , y (b)= β ˜ , α ˜ , β ˜ E 1 ,
(2.3)
y (a) α 0 y(a)= α ˜ , y (b)+ β 0 y(b)= β ˜ , α 0 , β 0 0, α 0 + β 0 >0,
(2.4)

are called the second-order fuzzy boundary value problems (FBVPs). The differential Eq. (2.1) along with fuzzy boundary value conditions (2.2), (2.3) and (2.4) are said to be second-order fuzzy differential equation No. 1, No. 2 and No. 3 boundary value problems, respectively.

In particular, when f(t,y, y ) is a linear function with respect to y and y , Eq. (2.1) is reduced to Eq. (1.1)

y +p(t) y +q(t)y=g(t),t[a,b],

and it is a linear differential equation. In this paper, we discuss the approximate solution of the second-order fuzzy linear differential function boundary value problem. For simplicity, we only discuss the second-order fuzzy linear differential function with fuzzy boundary value conditions (2.3) and (2.4).

3 Method for solving No. 2 FBVPs

3.1 The undetermined fuzzy coefficients method

The undetermined fuzzy coefficients method is to seek an approximate solution as

y ˜ N (t)= k = 0 N θ ˜ k ϕ k (t),
(3.1)

where ϕ k (t), k=0,1,,N, are positive basic functions whose all differentiations are positive. We compute the fuzzy coefficients in (3.1) by setting the error to zero as follows:

E=D ( y ˜ + p ( t ) y ˜ + q ( t ) y ˜ , g ˜ ( t ) ) +D ( y ˜ ( a ) , α ˜ ) +D ( y ˜ ( b ) , β ˜ ) .
(3.2)

Then we substitute (3.1) in (3.2) and represent them in parametric forms, then

{ y ( t , r ) + p ( t ) y ( t , r ) + q ( t ) y ( t , r ) ̲ = g ( t , r ) ̲ , y ̲ ( a , r ) = α ̲ ( r ) , y ̲ ( b , r ) = β ̲ ( r ) , y ( t , r ) + p ( t ) y ( t , r ) + q ( t ) y ( t , r ) ¯ = g ( t , r ) ¯ , y ¯ ( a , r ) = α ¯ ( r ) , y ¯ ( b , r ) = β ¯ ( r ) .
(3.3)

Lemma 3.1 [11]

Let basic functions ϕ k (t), k=0,1,,N, and all of their differentiations be positive, without loss of generality. Then ( y ̲ N ) ( i ) (t)= y N ( i ) ̲ (t) and ( y ¯ N ( i ) )¯(t)= y N ( i ) ¯ (t), i=0,1,2.

In order to solve Eq. (1.1) with condition (1.3), we need to investigate the system of Eq. (3.3). In this section we consider two cases.

3.2 Case 1 p(t)q(t)0

Suppose that coefficients p(t), q(t) are nonnegative. From (3.3), we have

y ̲ (t,r)+p(t) y ̲ (t,r)+q(t) y ̲ (t,r)= g ̲ (t,r),
(3.4)
y ¯ (t,r)+p(t) y ¯ (t,r)+q(t) y ¯ (t,r)= g ¯ (t,r),
(3.5)
y ̲ (a,r)= α ̲ (r), y ̲ (b,r)= β ̲ (r), y ¯ (a,r)= α ¯ (r), y ¯ (b,r)= β ¯ (r).
(3.6)

And when coefficients p(t), q(t) are negative, we have

y ̲ (t,r)p(t) y ¯ (t,r)q(t) y ¯ (t,r)= g ̲ (t,r),
(3.7)
y ¯ (t,r)p(t) y ̲ (t,r)q(t) y ̲ (t,r)= g ¯ (t,r).
(3.8)

If (3.1) is substituted in (3.4) and (3.5), respectively, then

k = 0 N θ ̲ k (r) ϕ k (t)+p(t) k = 0 N θ ̲ k (r) ϕ k (t)+q(t) k = 0 N θ ̲ k (r) ϕ k (t)= g ̲ (t,r),

and

k = 0 N θ ¯ k (r) ϕ k (t)+p(t) k = 0 N θ ¯ k (r) ϕ k (t)+q(t) k = 0 N θ ¯ k (r) ϕ k (t)= g ¯ (t,r),

also

k = 0 N θ ̲ k ( r ) ϕ k ( a ) = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) ϕ k ( b ) = β ̲ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( a ) = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( b ) = β ¯ ( r ) .

By setting

{ γ k = ϕ k ( t ) + p ( t ) ϕ k ( t ) + q ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = ϕ k ( b ) , k=0,1,,N,

the following system is obtained:

{ k = 0 N θ ̲ k ( r ) γ k = g ̲ ( t , r ) , k = 0 N θ ̲ k ( r ) σ a k = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) σ b k = β ̲ ( r ) k = 0 N θ ¯ k ( r ) γ k = g ¯ ( t , r ) , k = 0 N θ ¯ k ( r ) σ a k = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) σ b k = β ¯ ( r ) .
(3.9)

Equation (3.9) is a system of linear equations S(t)X(r)=Y(r) such that

S= ( S 1 S 2 S 2 S 1 ) ,

where

S 1 = ( γ 0 γ 1 γ N σ a 0 σ a 1 σ a N σ b 0 σ b 1 σ b N ) , S 2 = ( 0 0 0 0 0 0 0 0 0 ) .

And

X = ( θ ̲ 0 , θ ̲ 1 , , θ ̲ N , θ ¯ 0 , θ ¯ 1 , , θ ¯ N ) , Y = ( g ̲ ( t , r ) , α ̲ ( r ) , β ̲ ( r ) , g ¯ ( t , r ) , α ¯ ( r ) , β ¯ ( r ) ) .

The variables θ ̲ 0 , θ ̲ 1 ,, θ ̲ N , θ ¯ 0 , θ ¯ 1 ,, θ ¯ N are obtained by solving (3.9) by setting t=s, s[a,b]. These variables yield the fuzzy approximate solution ( y ̲ (t,r), y ¯ (t,r)).

In the same way, when coefficients p(t), q(t) are negative, we build the corresponding system of linear equations as follows:

( γ 0 γ 1 γ N δ 0 δ 1 δ N σ a 0 σ a 1 σ a N 0 0 0 σ b 0 σ b 1 σ b N 0 0 0 δ 0 δ 1 δ N γ 0 γ 1 γ N 0 0 0 σ a 0 σ a 1 σ a N 0 0 0 σ b 0 σ b 1 σ b N ) ( θ ̲ 0 θ ̲ 1 θ ̲ N θ ¯ 0 θ ¯ 1 θ ¯ N ) = ( g ̲ ( t , r ) α ̲ ( r ) β ̲ ( r ) g ¯ ( t , r ) α ¯ ( r ) β ¯ ( r ) ) ,
(3.10)

where

{ γ k = ϕ k ( t ) , δ k = p ( t ) ϕ k ( t ) q ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = ϕ k ( b ) , k=0,1,,N.

3.3 Case 2 p(t)q(t)<0

Suppose that coefficient p(t) is nonnegative and q(t) is negative. From (3.3), we have

y ̲ (t,r)+p(t) y ̲ (t,r)q(t) y ¯ (t,r)= g ̲ (t,r),
(3.11)
y ¯ (t,r)+p(t) y ¯ (t,r)q(t) y ̲ (t,r)= g ¯ (t,r).
(3.12)

When coefficient p(t) is negative and q(t) is nonnegative, we have

y ̲ (t,r)p(t) y ¯ (t,r)+q(t) y ̲ (t,r)= g ̲ (t,r),
(3.13)
y ¯ (t,r)p(t) y ̲ (t,r)+q(t) y ¯ (t,r)= g ¯ (t,r).
(3.14)

If (3.1) is substituted in (3.11) and (3.12), respectively, then

k = 0 N θ ̲ k (r) ϕ k (t)+p(t) k = 0 N θ ̲ k (r) ϕ k (t)q(t) k = 0 N θ ¯ k (r) ϕ k (t)= g ̲ (t,r),

and

k = 0 N θ ¯ k (r) ϕ k (t)+p(t) k = 0 N θ ¯ k (r) ϕ k (t)q(t) k = 0 N θ ̲ k (r) ϕ k (t)= g ¯ (t,r),

also

k = 0 N θ ̲ k ( r ) ϕ k ( a ) = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) ϕ k ( b ) = β ̲ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( a ) = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( b ) = β ¯ ( r ) .

By setting

{ ζ k = ϕ k ( t ) + p ( t ) ϕ k ( t ) , ξ k = q ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = ϕ k ( b ) , k=0,1,,N,

the following system is obtained:

{ k = 0 N θ ̲ k ( r ) ζ k + k = 0 N θ ¯ k ( r ) ξ k = g ̲ ( t , r ) , k = 0 N θ ̲ k ( r ) σ a k = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) σ b k = β ̲ ( r ) k = 0 N θ ¯ k ( r ) ζ k + k = 0 N θ ̲ k ( r ) ξ k = g ¯ ( t , r ) , k = 0 N θ ¯ k ( r ) σ a k = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) σ b k = β ¯ ( r ) .
(3.15)

Equation (3.15) is a system of linear equations S(t)X(r)=Y(r) such that

S= ( S 1 S 2 S 2 S 1 ) ,

where

S 1 = ( ζ 0 ζ 1 ζ N σ a 0 σ a 1 σ a N σ b 0 σ b 1 σ b N ) , S 2 = ( ξ 0 ξ 1 ξ N 0 0 0 0 0 0 ) .

And

X = ( θ ̲ 0 , θ ̲ 1 , , θ ̲ N , θ ¯ 0 , θ ¯ 1 , , θ ¯ N ) , Y = ( g ̲ ( t , r ) , α ̲ ( r ) , β ̲ ( r ) , g ¯ ( t , r ) , α ¯ ( r ) , β ¯ ( r ) ) .

The variables θ ̲ 0 , θ ̲ 1 ,, θ ̲ N , θ ¯ 0 , θ ¯ 1 ,, θ ¯ N are obtained by solving (3.15) by setting t=s, s[a,b]. These variables yield the fuzzy approximate solution ( y ̲ (t,r), y ¯ (t,r)).

In the same way, when coefficient p(t) is negative and q(t) is nonnegative, we set up the corresponding system of linear equations as follows:

( ζ 0 ζ 1 ζ N ξ 0 ξ 1 ξ N σ a 0 σ a 1 σ a N 0 0 0 σ b 0 σ b 1 σ b N 0 0 0 ξ 0 ξ 1 ξ N ζ 0 ζ 1 ζ N 0 0 0 σ a 0 σ a 1 σ a N 0 0 0 σ b 0 σ b 1 σ b N ) ( θ ̲ 0 θ ̲ 1 θ ̲ N θ ¯ 0 θ ¯ 1 α ¯ N ) = ( g ̲ ( t , r ) α ̲ ( r ) β ̲ ( r ) g ¯ ( t , r ) α ¯ ( r ) β ¯ ( r ) ) ,
(3.16)

where

{ ζ k = ϕ k ( t ) + q ( t ) ϕ k ( t ) , ξ k = p ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = ϕ k ( b ) , k=0,1,,N.

4 Method for solving No. 3 FBVPs

4.1 The undetermined fuzzy coefficients method

Similarly, we compute the fuzzy coefficients in (3.1) by setting the error to zero as follows:

E=D ( y ˜ + p ( t ) y ˜ + q ( t ) y ˜ , g ˜ ( t ) ) +D ( y ˜ ( a ) α 0 y ˜ ( a ) , α ˜ ) +D ( y ˜ ( b ) + β 0 y ˜ ( b ) , β ˜ ) .
(4.1)

Then we substitute (3.1) in (4.1) and represent them in parametric forms:

{ y ( t , r ) + p ( t ) y ( t , r ) + q ( t ) y ( t , r ) ̲ = g ( t , r ) ̲ , y ̲ ( a , r ) + α 0 y ¯ ( a , r ) = α ̲ ( r ) , y ̲ ( b , r ) + β 0 y ̲ ( b , r ) = β ̲ ( r ) , y ( t , r ) + p ( t ) y ( t , r ) + q ( t ) y ( t , r ) ¯ = g ( t , r ) ¯ , y ¯ ( a , r ) + α 0 y ̲ ( a , r ) = α ¯ ( r ) , y ¯ ( b , r ) + β 0 y ¯ ( b , r ) = β ¯ ( r ) .
(4.2)

In order to solve Eq. (1.1) with Eq. (1.4), we need to investigate the system of Eq. (4.3). In this section we also consider two cases.

4.2 Case 1 p(t)q(t)0

Suppose that coefficients p(t), q(t) are nonnegative. From (4.2), we have

y ̲ (t,r)+p(t) y ̲ (t,r)+q(t) y ̲ (t,r)= g ̲ (t,r),
(4.3)
y ¯ (t,r)+p(t) y ¯ (t,r)+q(t) y ¯ (t,r)= g ¯ (t,r),
(4.4)
y ̲ (a,r)+ α 0 y ¯ (a,r)= α ̲ (r), y ̲ (b,r)+ β 0 y ̲ (b,r)= β ̲ (r),
(4.5)
y ¯ (a,r)+ α 0 y ̲ (a,r)= α ¯ (r), y ¯ (b,r)+ β 0 y ¯ (b,r)= β ¯ (r).
(4.6)

And when coefficients p(t), q(t) are negative, we have

y ̲ (t,r)p(t) y ¯ (t,r)q(t) y ¯ (t,r)= g ̲ (t,r),
(4.7)
y ¯ (t,r)p(t) y ̲ (t,r)q(t) y ̲ (t,r)= g ¯ (t,r).
(4.8)

If (3.1) is substituted in (4.3) and (4.4), respectively, then

k = 0 N θ ̲ k (r) ϕ k (t)+p(t) k = 0 N θ ̲ k (r) ϕ k (t)+q(t) k = 0 N θ ̲ k (r) ϕ k (t)= g ̲ (t,r),

and

k = 0 N θ ¯ k (r) ϕ k (t)+p(t) k = 0 N θ ¯ k (r) ϕ k (t)+q(t) k = 0 N θ ¯ k (r) ϕ k (t)= g ¯ (t,r),

also

k = 0 N θ ̲ k ( r ) ϕ k ( a ) + α 0 k = 0 N θ ¯ k ( r ) ϕ k ( a ) = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) ϕ k ( b ) + β 0 k = 0 N θ ̲ k ( r ) ϕ k ( b ) = β ̲ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( a ) + α 0 k = 0 N θ ̲ k ( r ) ϕ k ( a ) = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( b ) + β 0 k = 0 N θ ¯ k ( r ) ϕ k ( b ) = β ¯ ( r ) .

By setting

{ γ k = ϕ k ( t ) + p ( t ) ϕ k ( t ) + q ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = α 0 ϕ k ( a ) , δ ( k ) = ϕ k ( b ) + β 0 ϕ k ( b ) , k=0,1,,N,

the following system is obtained:

{ k = 0 N θ ̲ k ( r ) γ k = g ̲ ( t , r ) , k = 0 N θ ̲ k ( r ) σ a k + k = 0 N θ ¯ k ( r ) σ b k = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) δ k = β ̲ ( r ) k = 0 N θ ¯ k ( r ) γ k = g ¯ ( t , r ) , k = 0 N θ ¯ k ( r ) σ a k + k = 0 N θ ̲ k ( r ) σ b k = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) δ k = β ¯ ( r ) .
(4.9)

Equation (4.9) is a system of linear equations S(t)X(r)=Y(r) such that

S= ( S 1 S 2 S 2 S 1 ) ,

where

S 1 = ( γ 0 γ 1 γ N σ a 0 σ a 1 σ a N δ b 0 δ b 1 δ b N ) , S 2 = ( 0 0 0 σ b 0 σ b 1 σ b N 0 0 0 ) .

And

X = ( θ ̲ 0 , θ ̲ 1 , , θ ̲ N , θ ¯ 0 , θ ¯ 1 , , θ ¯ N ) , Y = ( g ̲ ( t , r ) , α ̲ ( r ) , β ̲ ( r ) , g ¯ ( t , r ) , α ¯ ( r ) , β ¯ ( r ) ) .

The variables θ ̲ 0 , θ ̲ 1 ,, θ ̲ N , θ ¯ 0 , θ ¯ 1 ,, θ ¯ N are obtained by solving (4.9) by setting t=s, s[a,b]. These variables yield the fuzzy approximate solution ( y ̲ (t,r), y ¯ (t,r)).

Similarly, when coefficients p(t), q(t) are negative, we build the corresponding system of linear equations as follows:

( γ 0 γ 1 γ N η 0 η 1 η N σ a 0 σ a 1 σ a N σ b 0 σ b 1 σ b N δ 0 δ 1 δ N 0 0 0 η 0 η 1 η N γ 0 γ 1 γ N σ b 0 σ b 1 σ b N σ a 0 σ a 1 σ a N 0 0 0 δ 0 δ 1 δ N ) ( θ ̲ 0 θ ̲ 1 θ ̲ N θ ¯ 0 θ ¯ 1 θ ¯ N ) = ( g ̲ ( t , r ) α ̲ ( r ) β ̲ ( r ) g ¯ ( t , r ) α ¯ ( r ) β ¯ ( r ) ) ,
(4.10)

where

{ γ k = ϕ k ( t ) , η k = p ( t ) ϕ k ( t ) q ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = α 0 ϕ k ( a ) , δ ( k ) = ϕ k ( b ) + β 0 ϕ k ( b ) , k=0,1,,N,

4.3 Case 2 p(t)q(t)<0

Suppose that coefficient p(t) is nonnegative and q(t) is negative. From (4.3), we have

y ̲ (t,r)+p(t) y ̲ (t,r)q(t) y ¯ (t,r)= g ̲ (t,r),
(4.11)
y ¯ (t,r)+p(t) y ¯ (t,r)q(t) y ̲ (t,r)= g ¯ (t,r).
(4.12)

When coefficient p(t) is negative and q(t) is nonnegative, we have

y ̲ (t,r)p(t) y ¯ (t,r)+q(t) y ̲ (t,r)= g ̲ (t,r),
(4.13)
y ¯ (t,r)p(t) y ̲ (t,r)+q(t) y ¯ (t,r)= g ¯ (t,r).
(4.14)

If (3.1) is substituted in (4.11) and (4.12), respectively, then

k = 0 N θ ̲ k (r) ϕ k (t)p(t) k = 0 N θ ̲ k (r) ϕ k (t)+q(t) k = 0 N θ ¯ k (r) ϕ k (t)= g ̲ (t,r),

and

k = 0 N θ ¯ k (r) ϕ k (t)p(t) k = 0 N θ ¯ k (r) ϕ k (t)+q(t) k = 0 N θ ̲ k (r) ϕ k (t)= g ¯ (t,r),

also

k = 0 N θ ̲ k ( r ) ϕ k ( a ) + α 0 k = 0 N θ ¯ k ( r ) ϕ k ( a ) = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) ϕ k ( b ) + β 0 k = 0 N θ ̲ k ( r ) ϕ k ( b ) = β ̲ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( a ) + α 0 k = 0 N θ ̲ k ( r ) ϕ k ( a ) = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) ϕ k ( b ) + β 0 k = 0 N θ ¯ k ( r ) ϕ k ( b ) = β ¯ ( r ) .

By setting

{ γ k = ϕ k ( t ) + p ( t ) ϕ k ( t ) , ξ k = q ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = α 0 ϕ k ( a ) , δ ( k ) = ϕ k ( b ) + β 0 ϕ k ( b ) , k=0,1,,N,

the following system is obtained:

{ k = 0 N θ ̲ k ( r ) γ k + k = 0 N θ ¯ k ( r ) ξ k = g ̲ ( t , r ) , k = 0 N θ ̲ k ( r ) σ a k + k = 0 N θ ¯ k ( r ) σ b k = α ̲ ( r ) , k = 0 N θ ̲ k ( r ) δ k = β ̲ ( r ) , k = 0 N θ ¯ k ( r ) γ k + k = 0 N θ ̲ k ( r ) ξ k = g ¯ ( t , r ) , k = 0 N θ ¯ k ( r ) σ a k + k = 0 N θ ̲ k ( r ) σ b k = α ¯ ( r ) , k = 0 N θ ¯ k ( r ) δ k = β ¯ ( r ) .
(4.15)

Equation (4.15) is a system of linear equations S(t)X(t)=Y(r) such that

S= ( S 1 S 2 S 2 S 1 ) ,

where

S 1 = ( γ 0 γ 1 γ N σ a 0 σ a 1 σ a N δ 0 δ 1 δ N ) , S 2 = ( ξ 0 ξ 1 ξ N σ b 0 σ b 1 σ b N 0 0 0 ) .

And

X = ( θ ̲ 0 , θ ̲ 1 , , θ ̲ N , θ ¯ 0 , θ ¯ 1 , , θ ¯ N ) , Y = ( g ̲ ( t , r ) , α ̲ ( r ) , β ̲ ( r ) , g ¯ ( t , r ) , α ¯ ( r ) , β ¯ ( r ) ) .

The variables θ ̲ 0 , θ ̲ 1 ,, θ ̲ N , θ ¯ 0 , θ ¯ 1 ,, θ ¯ N are obtained by solving (4.14) by setting t=s, s[a,b]. These variables yield the fuzzy approximate solution ( y ̲ (t,r), y ¯ (t,r)).

Similarly, when coefficient p(t) is negative and q(t) is nonnegative, we extend the corresponding system of linear equations as follows:

( γ 0 γ 1 γ N ξ 0 ξ 1 ξ N σ a 0 σ a 1 σ a N σ b 0 σ b 1 σ b N δ 0 δ 1 δ N 0 0 0 ξ 0 ξ 1 ξ N γ 0 γ 1 γ N σ b 0 σ b 1 σ b N σ a 0 σ a 1 σ a N 0 0 0 δ 0 δ 1 δ N ) ( θ ̲ 0 θ ̲ 1 θ ̲ N θ ¯ 0 θ ¯ 1 θ ¯ N ) = ( g ̲ ( t , r ) α ̲ ( r ) β ̲ ( r ) g ¯ ( t , r ) α ¯ ( r ) β ¯ ( r ) ) ,
(4.16)

where

{ γ k = ϕ k ( t ) + q ( t ) ϕ k ( t ) , ξ k = p ( t ) ϕ k ( t ) , σ a k = ϕ k ( a ) , σ b k = α 0 ϕ k ( a ) , δ ( k ) = ϕ k ( b ) + β 0 ϕ k ( b ) , k=0,1,,N.

Likewise, for Eq. (1.1) with fuzzy boundary conditions (1.2), the following results are obvious.

5 Approximate solutions of second-order FLBVPs

The above model linear Eqs. (3.9), (3.10), (3.15), (3.16), (4.9), (4.10), (4.15) and (4.16) are 6×2(N+1) function systems of linear equations and they have the same form as follows:

S(t)X(r)=Y(r).
(5.1)

In the process of solving Eq. (5.1) by setting t=q, q[a,b], no matter whether it is consistent or inconsistent, we obtain the minimal norm least squares solution [20] by using the generalized inverse of the coefficient matrix S, i.e.,

X(r)= S (q)Y(r).
(5.2)

Thus we get

θ ̲ 0 , θ ̲ 1 ,, θ ̲ N , θ ¯ 0 , θ ¯ 1 ,, θ ¯ N .

Therefore, we obtain the fuzzy approximate solution of the original fuzzy linear differential equation as follows:

y ̲ ( t , r ) = θ ̲ 0 ( r ) ϕ 0 ( t ) + θ ̲ 1 ( r ) ϕ 1 ( t ) + + θ ̲ N ( r ) ϕ N ( t ) , y ¯ ( t , r ) = θ ¯ 0 ( r ) ϕ 0 ( t ) + θ ¯ 1 ( r ) ϕ 1 ( t ) + + θ ¯ N ( r ) ϕ N ( t ) .
(5.3)

6 Numerical examples

Example 6.1 Consider the following second-order fuzzy linear differential equation:

{ y + y = t , t [ 0 , π 2 ] y ˜ ( 0 ) = ( 0.1 + 0.1 r , 0.1 0.1 r ) , y ˜ ( π 2 ) = ( π 2 + 0.1 r , 1 + π 2 0.1 r ) .

The exact solution is as follows:

Y ̲ ( t , r ) = ( 0.1 + 0.1 r ) cos t + ( 0.1 r ) sin t t , Y ¯ ( t , r ) = ( 0.1 0.1 r ) cos t + ( 1 + π 0.1 r ) sin t t .

If ϕ k (t)= t k , k=0,1,2,3, then

y ̲ (t,r)= α ̲ 0 + α ̲ 1 t+ α ̲ 2 t 2 + α ̲ 3 t 3 , y ¯ (t,r)= α ¯ 0 + α ¯ 1 t+ α ¯ 2 t 2 + α ¯ 3 t 3 .

From (3.11), we build the following system:

( 1 t 2 + t 2 6 t + t 3 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 t 2 + t 2 6 t + t 3 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 ) ( θ ̲ 0 ( r ) θ ̲ 1 ( r ) θ ̲ 2 ( r ) θ ̲ 3 ( r ) θ ¯ 0 ( r ) θ ¯ 1 ( r ) θ ¯ 2 ( r ) θ ¯ 3 ( r ) ) = ( t 0.1 + 0.1 r π 2 + 0.1 r t 0.1 0.1 r 1 + π 2 0.1 r ) .

By setting t= 1 3 , the parameters α ̲ 0 , α ̲ 1 , α ̲ 2 , α ̲ 3 , α ¯ 0 , α ¯ 1 , α ¯ 2 , α ¯ 3 are obtained, and by putting them into (5.3), we have

y ̲ ( t , r ) = ( 0.100 + 0.100 r ) + ( 1.5708 + 0.100 r ) t y ( t , r ) = + ( 0.0674 0.0229 r ) t 2 + ( 0.0936 0.03165 r ) t 3 , y ¯ ( t , r ) = ( 0.100 0.100 r ) + ( 4.142 0.100 r ) t y ¯ ( t , r ) = + ( 0.1128 + 0.0227 r ) t 2 + ( 0.1568 + 0.03161 r ) t 3 .

Tables 1, 2, 3 and 4 show the comparisons between the exact solution and the approximate solution at t=0.01 and t=0.001 for some r[0,1]; all data are computed by Matlab7.x.

Table 1 Comparisons between the exact solution and the approximate solution
Table 2 Comparisons between the exact solution and the approximate solution
Table 3 Comparisons between the exact solution and the approximate solution
Table 4 Comparisons between the exact solution and the approximate solution

Example 6.2 Consider the following second-order fuzzy linear differential equation:

{ y 4 y + 4 y = 4 t 4 , t 0 y ˜ ( 0 ) = ( 3 + 2 r , 6 5 r ) , y ˜ ( 1 ) = ( ( 1 + 5 r ) e 2 + 1 , ( 7 5 r ) e 2 + 1 ) .

The exact solution of the equation is

Y ̲ ( t , r ) = ( 1 + r ) e 2 t + ( 1 + r ) t e 2 t + t , Y ¯ ( t , r ) = ( 2 r ) e 2 t + ( 1 r ) t e 2 t + t .

The extended linear equations S(t)X(r)=Y(r) is as follows:

( 4 4 t 2 + 4 t 2 6 t + 4 t 3 0 4 8 t 12 t 2 0 1 0 0 0 0 0 0 0 1 2 0 0 0 0 0 0 4 8 t 12 t 2 4 4 t 2 + 4 t 2 6 t + 4 t 3 0 0 0 0 0 1 0 0 0 0 0 0 0 1 2 0 ) ( α ̲ 0 ( r ) α ̲ 1 ( r ) α ̲ 2 ( r ) α ̲ 3 ( r ) α ¯ 0 ( r ) α ¯ 1 ( r ) α ¯ 2 ( r ) α ¯ 3 ( r ) ) = ( 4 t 4 3 + 2 r ( 1 + 5 r ) e 2 + 1 4 t 4 6 5 r ( 7 5 r ) e 2 + 1 ) .

By setting t= 1 6 , the parameters α ̲ 0 (r), α ̲ 1 (r), α ̲ 2 (r), α ̲ 3 (r), α ¯ 0 (r), α ¯ 1 (r), α ¯ 2 (r), α ¯ 3 (r) are obtained. Tables 5 and 6 show comparisons between the exact solution and the approximate solution at t=0.001 for some r[0,1]; all data were calculated by Matlab7.x.

Table 5 Comparisons between the exact solution and the approximate solution
Table 6 Comparisons between the exact solution and the approximate solution

Form Tables 1, 2, 3, 4, 5 and 6, we know that the approximate solutions obtained from the proposed method are best close to the exact solutions of original linear deferential equations with fuzzy boundary value conditions.

7 Conclusion

In this paper the approximate method similar to the undetermined coefficients method, based on a positive basis for solving second-order fuzzy linear boundary value problems, was discussed. Three classes of boundary conditions and the general case were considered. According to the sign of coefficient functions of the fuzzy linear differential equation, the corresponding function systems of linear equations were set up. Following each other, fuzzy approximate solutions were obtained by solving a crisp function extended system of linear equations. Numerical examples show that our methods are practical and efficient.

References

  1. Zadeh LA: The concept of a linguistic variable and its application to approximate reasoning. Inf. Sci. 1975, 8: 199-249. 10.1016/0020-0255(75)90036-5

    Article  MathSciNet  MATH  Google Scholar 

  2. Dubois D, Prade H: Operations on fuzzy numbers. Int. J. Syst. Sci. 1978, 9: 613-626. 10.1080/00207727808941724

    Article  MathSciNet  MATH  Google Scholar 

  3. Puri ML, Ralescu DA: Differentials for fuzzy functions. J. Math. Anal. Appl. 1983, 91: 552-558. 10.1016/0022-247X(83)90169-5

    Article  MathSciNet  MATH  Google Scholar 

  4. Goetschel R, Voxman W: Elementary calculus. Fuzzy Sets Syst. 1986, 18: 31-43. 10.1016/0165-0114(86)90026-6

    Article  MathSciNet  MATH  Google Scholar 

  5. Kandel A: Fuzzy dynamic systems and the nature of their solutions. In Fuzzy Sets Theory and Application to Policy Analysis and Information Systems. Edited by: Wang PP, Chang SK. Plenum Press, New York; 1980:93-1221.

    Google Scholar 

  6. Kandel A, Byatt WJ: Fuzzy differential equations. Proceedings of the International Conference on Cybernetics and Society 1987, 1213-1216. Tokyo

    Google Scholar 

  7. Abbasbandy S, Allahviranloo T: Numerical solutions of fuzzy differential equations by Taylor method. Comput. Methods Appl. Math. 2002, 2: 113-124.

    Article  MathSciNet  Google Scholar 

  8. Abbasbandy S, Allahviranloo T, Lopez-Pouso O, Otadi M: Numerical methods for fuzzy differential inclusions. Comput. Math. Appl. 2004, 48: 1633-1641. 10.1016/j.camwa.2004.03.009

    Article  MathSciNet  MATH  Google Scholar 

  9. Allahviranloo T, Ahmady N, Ahmady E: Numerical solutions of fuzzy differential equations by predictor-corrector method. Inf. Sci. 2007, 177: 1633-1647. 10.1016/j.ins.2006.09.015

    Article  MathSciNet  MATH  Google Scholar 

  10. Buckley JJ, Feuring T: Fuzzy initial value problem for Nth-order fuzzy differential equations. Fuzzy Sets Syst. 2001, 121: 247-255. 10.1016/S0165-0114(00)00028-2

    Article  MathSciNet  MATH  Google Scholar 

  11. Allahviranloo T: Nth-order fuzzy linear differential equations. Inf. Sci. 2008, 178: 1309-1324. 10.1016/j.ins.2007.10.013

    Article  MathSciNet  MATH  Google Scholar 

  12. O’Regan D, Lakshkantham V, Nieto JJ: Initial and boundary value problems for fuzzy differential equations. Nonlinear Anal. 2003, 54: 405-415. 10.1016/S0362-546X(03)00097-X

    Article  MathSciNet  MATH  Google Scholar 

  13. Lakshkantham V, Murty KN, Turner J: Two points boundary value problems associated with nonlinear fuzzy differential equations. Math. Inequal. Appl. 2003, 4: 527-533.

    Google Scholar 

  14. Chen M, Wu C, Xue X: On fuzzy boundary value problems. Inf. Sci. 2008, 178: 1877-1892. 10.1016/j.ins.2007.11.017

    Article  MathSciNet  MATH  Google Scholar 

  15. Wu CX, Ma M: Embedding problem of fuzzy number space: Part I. Fuzzy Sets Syst. 1991, 44: 33-38. 10.1016/0165-0114(91)90030-T

    Article  MathSciNet  MATH  Google Scholar 

  16. Wu CX, Ma M: Embedding problem of fuzzy number space: Part III. Fuzzy Sets Syst. 1992, 46: 281-286. 10.1016/0165-0114(92)90142-Q

    Article  MathSciNet  MATH  Google Scholar 

  17. He O, Yi W: On fuzzy differential equations. Fuzzy Sets Syst. 1989, 24: 321-325.

    Article  MathSciNet  MATH  Google Scholar 

  18. Bede B: Note on “Numerical solutions of fuzzy differential equations by predictor-corrector method”. Inf. Sci. 2008, 178: 1917-1922. 10.1016/j.ins.2007.11.016

    Article  MathSciNet  MATH  Google Scholar 

  19. Bede B, Gal S: Generalizations of the differentiability of fuzzy-number-valued functions with applications to fuzzy differential equation. Fuzzy Sets Syst. 2005, 151: 581-599. 10.1016/j.fss.2004.08.001

    Article  MathSciNet  MATH  Google Scholar 

  20. Ben-Israel A, Greville TNE: Generalized Inverses: Theory and Applications. 2nd edition. Springer, New York; 2003.

    MATH  Google Scholar 

Download references

Acknowledgements

The work is supported by the Natural Scientific Funds of PR China (71061013, 21175108) and the Youth Research Ability Project of Northwest Normal University (NWNU-LKQN-1120).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Xiaobin Guo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally and significantly in writing this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Guo, X., Shang, D. & Lu, X. Fuzzy approximate solutions of second-order fuzzy linear boundary value problems. Bound Value Probl 2013, 212 (2013). https://doi.org/10.1186/1687-2770-2013-212

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-212

Keywords