# Fuzzy Number Addition with the Application of Horizontal Membership Functions

Authors: Andrzej Piegat; Marcin Pluciński
Journal: The Scientific World Journal
 (2015)

---

## Abstract

The paper presents addition of fuzzy numbers realised with the application of the multidimensional RDM arithmetic and horizontal membership functions (MFs). Fuzzy arithmetic (FA) is a very difficult task because operations should be performed here on multidimensional information granules. Instead, a lot of FA methods useα-cuts in connection with 1-dimensional classical interval arithmetic that operates not on multidimensional granules but on 1-dimensional intervals. Such approach causes difficulties in calculations and is a reason for arithmetical paradoxes. The multidimensional approach allows for removing drawbacks and weaknesses of FA. It is possible thanks to the application of horizontal membership functions which considerably facilitate calculations because now uncertain values can be inserted directly into equations without using the extension principle. The paper shows how the addition operation can be realised on independent fuzzy numbers and on partly or fully dependent fuzzy numbers with taking into account the order relation and how to solve equations, which can be a difficult task for 1-dimensional FAs.


---

## Body

## 1. Introduction

Fuzzy arithmetic [1–14] is used in uncertainty theory [15–20], grey systems [21, 22], granular computing [9, 23], computing with words [23–26], decision-making [27, 28], and other sciences and engineering branches [4, 26, 29, 30]. Authors of first concepts of the fuzzy arithmetic (shortly FA) based on L-R (left-right) fuzzy numbers were Dubois and Prade [1]. With years, FA has been improved and its new versions have been introduced, for example, the popular α-cuts’ version [4, 6, 31], which in this paper will be called μ-cuts’ FA. In general, all versions of FA can be divided [4] into elementary FA, standard FA, and advanced FA versions. Examples of advanced FA methods can be the generalized vertex method [32], constrained FA [7, 8, 21], algorithmic FA [9], transformation method and extended transformation method [4], and inverse FA [4]. It seems that mostly used FA versions are based on α-cuts [2] and on Moore’s interval arithmetic (shortly IA) [33–35]. In the FA literature, many examples of practical FA applications can be found, for example, [4, 26, 27, 36, 37]. Practical problems require effective FAs which would enable solving uncertain linear and nonlinear equation systems, differential equations, integral calculations, and so forth [38, 39]. Therefore, many scientists investigate these problems and publish achieved results from this area [40–64]. However, almost always they emphasise that the achieved level of FA is not satisfactory and that further investigations are necessary. Therefore, investigations on IA and FA have been continued nonstop. It is testified by new publications in journals, conferences, and new books. Though FA has achieved many application successes, it further on has many weak points. For example, Dymova and Sevastjanov report in their papers [27, 61, 62] the present FA has considerable difficulties in solving even simple equations with one unknown variable. There exist also difficulties in defining neutral and inverse elements for addition and multiplication.The most popular version of FA is the arithmetic based onα-cuts, where an α-cut of a fuzzy set A, denoted as A
α, is defined by (1), where X means domain of the set and x means an element of this domain [4]. Consider(1)
A
α
=
c
u
t
α
A
~
=
x
∈
X
∣
μ
A
~
x
≥
α
.FA based onα-cuts uses principles of the classical Moore’s interval arithmetic [33, 34] for realisation of elementary arithmetic operations such as addition, subtraction, multiplication, and division. Though this arithmetic has many applications, its possibilities are limited, because of some drawbacks, which are commented on also by Moore himself in his books. For example, if X, Y, and Z are intervals, the distributive law,(2)
X
·
Y
+
Z
=
X
·
Y
+
X
·
Z
,does not always hold. Additionally, an additive inverse of an interval is not defined. If X
=
[
x
_
,
x
¯
], then X
+
(
-
X
)
=
[
x
_
-
x
¯
,
x
¯
-
x
_
]
≠
[
0,0
].As Dymowa shows in [27], there are great difficulties in solving even simple interval equations with one unknown. Let us assume that in a system there exists the dependence a
+
x
=
c. This dependence can be presented in few equivalent forms:(3)
a
+
x
=
c
,
a
=
c
-
x
,
x
=
c
-
a
,
a
+
x
-
c
=
0
.If only approximate values ofa and c are known, a
∈
[
1,3
] and c
∈
[
2,5
], and value x is not known, then, on the basis of (3), four interval equations (4)–(7) can be written:(4)
1,3
+
x
_
,
x
¯
=
2,5
,
(5)
1,3
=
2,5
-
x
_
,
x
¯
,
(6)
x
_
,
x
¯
=
2,5
-
1,3
,
(7)
1,3
+
x
_
,
x
¯
-
2,5
=
0
.Solving (4), two equations, 1
+
x
_
=
2 and 3
+
x
¯
=
5, are obtained. They give the solution [
x
_
,
x
¯
]
=
[
1,2
]. Solving (5), two equations, 1
=
2
-
x
¯ and 3
=
5
-
x
_, are obtained. They give the solution [
x
_
,
x
¯
]
=
[
2,1
] being an inverse interval. Solving (6), the solution [
-
1,4
] is obtained and from (7) we get the solution [
4
,
-
1
], which is also an inverse interval. Thus, various solutions, [
1,2
], [
2,1
], [
-
1,4
], and [
4,1
], are obtained depending on extension form of the original dependence. In the case of more complicated mathematical dependencies, the solution number can be considerably higher.Another paradox of FA based onμ-cuts and Moore arithmetic (shortly μ-cuts’ FA) [33–35] is not satisfying the cancellation law for multiplication. For example, equation X
Z
=
Y
Z in the general case does not mean that X
=
Z. It can be testified by an example shown below in which notation [
1,2
,
3
] means the triangle membership function (shortly MF) with support beginning, core position, and support end. Consider(8)
X
Z
=
1,2
,
3
·
-
1,0
,
1
=
-
3,0
,
3
=
Y
Z
=
2,2.5,3
·
-
1,0
,
1
=
-
3,0
,
3
,but(9)
X
=
1,2
,
3
≠
Y
=
2,2.5,3
.Because of this feature ofμ-cuts’ FA, transformations of formulas are not allowable. Why? It will be shown further on. Next important paradox of μ-cuts’ FA is the observation that, during calculation of results of nonlinear formulas, for example, C
=
A
-
A
2, we obtain different, nonunique solutions depending on which form of the formula is used: C
1
=
A
-
A
2,  C
2
=
A
(
1
-
A
), or C
3
=
(
A
-
1
)
+
(
1
-
A
)
(
1
+
A
). For A
=
[
0,1
,
2
], three different solutions are obtained: C
1
=
[
-
4,0
,
2
],  C
2
=
[
-
2,0
,
2
], and C
3
=
[
-
4,0
,
4
]; see Figure 4. Which solution is correct? The above phenomenon means that each transformation of an equation form, in the case of μ-cuts’ FA, can change its solution and that solutions are not unique. Further on, it will be shown that in the case of the multidimensional RDM FA such paradoxes do not occur.The above examples show that the classical IA, which is a basis of theμ-cut version of FA, is not ideal, though it can solve certain problems. Therefore, it can and should be further developed. Further on, a version of FA that is based on horizontal membership functions (MFs) and on μ-cuts will be presented. It also applies the multidimensional RDM arithmetic (M-RDM arithmetic) which has been elaborated by Andrzej Piegat. It has been developed together with coworkers: Marek Landowski, Marcin Pluciński, and Karina Tomaszewska [10, 37, 65–70].


## 2. Horizontal Membership Function

Fuzzy systems use vertical membership functions which were introduced by Zadeh [23]. They have the following form: μ
=
f
(
x
), where x is an independent variable and μ is a dependent one. Figure 1 shows an example of the trapezoidal MF and formula (10) gives its description:(10)
x
∈
a
,
b
:
μ
x
=
x
-
a
b
-
a
,
x
∈
b
,
c
:
μ
x
=
1
,
x
∈
c
,
d
:
μ
x
=
d
-
x
d
-
c
.Figure 1
Trapezium membership function.Vertical MF realises the mappingx
→
μ. The present fuzzy arithmetic is based on just such MFs. The idea of horizontal MFs has been elaborated by Andrzej Piegat. In this paper, an example of a trapezium MF will be presented but horizontal MFs can be used for all types of MFs. A function μ
(
x
) is unambiguous in the direction of the variable μ (Figure 1) and ambiguous in the direction of x. Therefore, it seems impossible to define a membership function in the x-direction. The function from Figure 1 assigns two values of x, x
L
(
μ
) and x
R
(
μ
), for one value of μ. However, let us introduce the RDM variable: α
x
∈
[
0,1
]. This variable has meaning of the relative-distance-measure and allows for determining of any point between two borders x
L
(
μ
) and x
R
(
μ
) of the function (Figure 2). RDM variable α
x takes a value of zero on the left border and a value of 1 on the right border of the function. Between the left border and the right border it takes fractional values. The idea of RDM variables was successfully used in the multidimensional IA [37, 67–70]. The multidimensional IA has shown that full and precise solutions of granular problems have form of multidimensional granules that cannot be explained and understood in terms of 1-dimensional approaches.Figure 2
Visualization of the horizontal approach to fuzzy membership functions.
(a)
(b)Formula (11) defines the left border and the right border of the trapezium MF (Figure 2). Consider(11)
x
L
=
a
+
b
-
a
μ
,
x
R
=
d
-
d
-
c
μ
.RDM variableα
x allows for a gradual transition of points between the left border and the right border. The interval x
(
μ
) in Figure 2 is defined by the following formula:(12)
x
=
x
L
+
x
R
-
x
L
α
x
,
α
x
∈
0,1
,
x
=
a
+
b
-
a
μ
+
d
-
a
-
μ
b
-
a
+
d
-
c
α
x
.It can be noted in formula (12) that x
=
f
(
μ
,
α
x
) is an unambiguous function existing in the 3D space, which can be seen in Figure 3.Figure 3
The horizontal membership functionx
=
(
1
+
2
μ
)
+
(
4
-
3
μ
)
α
x, where α
x
,
μ
∈
[
0,1
], corresponding to the function from Figure 2, in the 3D space as the unambiguous function.Figure 4
Three different membership functions of “results” of the fuzzy formulaC
=
A
-
A
2 written in three equivalent forms C
1
=
A
-
A
2,  C
2
=
A
(
1
-
A
),  C
3
=
(
A
-
1
)
+
(
1
-
A
)
(
1
+
A
) obtained with μ-cuts’ fuzzy arithmetic, where A
=
[
μ
,
2
-
μ
].It should be also noticed that RDM variablesα
i introduce the continuous Cartesian coordinate system in interval and fuzzy arithmetic calculations. In Moore’s arithmetic and in μ-cuts’ FA based on it, only borders of intervals or μ-cuts, for example, [
a
,
b
]
+
[
c
,
d
]
=
[
a
+
b
,
c
+
d
], take part in calculations. Insides of intervals are not taken into account. This fact hinders solving of more complicated problems and results in many paradoxes observed in IA and FA calculations. It also deprives IA and FA of many mathematical properties which the conventional mathematics of crisp numbers has. Few of these properties will be presented further on.Thanks to the fact that RDM variablesα
i introduce in interval calculations local and continuous Cartesian coordinate system, RDM FA possesses almost the same mathematical properties as arithmetic of crisp numbers.


### 2.1. Additive Inverse Element-
X of Element X

Inμ-cuts’ FA, trapezoidal membership function can be defined as a fuzzy interval with left border x
L
(
μ
)
=
a
+
(
b
-
a
)
μ and right border x
R
(
μ
)
=
d
-
(
d
-
c
)
μ (formula (13) and Figure 2). Consider(13)
X
μ
=
x
L
μ
,
x
R
μ
.A negative element-
X is determined by the following formula:(14)
-
X
μ
=
-
x
L
μ
,
x
R
μ
=
-
x
R
μ
,
-
x
L
μ
.The subtraction resultX
(
μ
)
-
X
(
μ
) determined by (15) is a fuzzy interval of nonzero span which means that, in μ-cuts’ FA additive, inverse element does not exist. Consider(15)
X
μ
-
X
μ
=
-
x
R
μ
-
x
L
μ
,
x
R
μ
-
x
L
μ
≠
0,0
.In the case of RDM FA, the elementX
(
μ
,
α
x
) is determined by (16) and its negative element -
X
(
μ
,
α
x
) by (17). Consider(16)
X
μ
,
α
x
=
x
μ
,
α
x
:
x
μ
,
α
x
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
,
α
x
∈
0,1
,
(17)
-
X
μ
,
α
x
=
x
μ
,
α
x
:
x
μ
,
α
x
=
-
x
L
μ
-
α
x
x
R
μ
-
x
L
μ
,
α
x
∈
0,1
.It is easy to check that the subtraction resultX
-
X is exactly equal to zero:(18)
X
μ
,
α
x
-
X
μ
,
α
x
=
0
.Thus, in the case of RDM FA, there exists the inverse additive element-
X
(
μ
,
α
x
).


### 2.2. Distributive Law of Multiplication

In theμ-cuts’ FA, similarly as in Moore’s interval arithmetic, the distributive law holds only in the limited form:(19)
X
μ
Y
μ
+
Z
ν
⊆
X
μ
Y
μ
+
X
μ
Z
μ
.In the multidimensional, continuous RDM FA, the distributive law holds fully. Consider(20)
X
μ
,
α
x
Y
μ
,
α
y
+
Z
μ
,
α
z
=
X
μ
,
α
x
Y
μ
,
α
y
+
X
μ
,
α
x
Z
μ
,
α
z
,
α
x
,
α
y
,
α
z
∈
0,1
.Proof.
For any three intervals expressed in RDM notation, withα
x
,
α
y
,
α
z
∈
[
0,1
],(21)
X
μ
,
α
x
=
x
L
μ
,
α
x
,
x
R
μ
,
α
x
=
x
μ
,
α
x
:
x
μ
,
α
x
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
,
Y
μ
,
α
y
=
y
L
μ
,
α
y
,
y
R
μ
,
α
y
=
y
μ
,
α
y
:
y
μ
,
α
y
=
y
L
μ
+
α
y
y
R
μ
-
y
L
μ
,
Z
μ
,
α
z
=
z
L
μ
,
α
z
,
z
R
μ
,
α
z
=
z
μ
,
α
z
:
z
μ
,
α
z
=
z
L
μ
+
α
z
z
R
μ
-
z
L
μ
,we have(22)
X
μ
,
α
x
Y
μ
,
α
y
+
Z
μ
,
α
z
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
·
y
L
μ
+
α
y
y
R
μ
-
y
L
μ
+
z
L
μ
+
α
z
z
R
μ
-
z
L
μ
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
y
L
μ
+
α
y
y
R
μ
-
y
L
μ
+
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
z
L
μ
+
α
z
z
R
μ
-
z
L
μ
=
X
μ
,
α
x
Y
μ
,
α
y
+
X
μ
,
α
x
Z
μ
,
α
y
.


### 2.3. Different Forms of Nonlinear Formulas

Let us consider results ofμ-cuts’ FA and RDM FA for the nonlinear formula C
=
A
-
A
2, where A is the fuzzy interval A
=
[
0,1
,
2
] (Figure 4). This formula can take at least 3 forms: C
1
=
A
-
A
2, C
2
=
A
(
1
-
A
), and C
3
=
(
A
-
1
)
+
(
1
-
A
)
(
1
+
A
). Fuzzy interval A has the form A
=
[
a
L
(
μ
)
,
a
R
(
μ
)
]
=
[
μ
,
2
-
μ
]. Calculating value of C for particular forms, C
1, C
2, and C
3, with the application of μ-cuts’ FA, three different results are obtained:(23)
C
1
=
-
4
+
3
μ
-
μ
2
,
2
-
μ
-
μ
2
,
with  support  span
-
4,2
,
C
2
=
-
2
+
3
μ
-
μ
2
,
2
-
3
μ
+
μ
2
,
with  support  span
-
2,2
,
C
3
=
-
4
+
5
μ
-
μ
2
,
4
-
5
μ
+
μ
2
,
with  support  span
-
4,4
.Each of the results,C
1, C
2, and C
3, is different. For μ
=
0, results mean spans (widths) of MF supports. They are equal to [
-
4,2
], [
-
2,2
], and [
-
4,4
]. The differences are considerable. Thus, the formula C
=
A
-
A
2 cannot be calculated uniquely with μ-cuts’ FA. If we use the multidimensional RDM FA, then interval A is expressed by (24) and its square A
2 by (25), where a is generic variable of A. One should notice that functions expressing A and A
2 are fully coupled (correlated); that is, values of a and α
a have to be identical both in A and in A
2. Consider(24)
A
:
a
=
a
L
+
α
a
a
R
-
a
L
=
μ
+
2
α
a
1
-
μ
,
(25)
A
2
:
a
2
=
μ
+
2
α
a
1
-
μ
2
,
α
a
∈
0,1
.FormulaC
1
=
A
-
A
2 can be determined by(26)
C
1
=
A
-
A
2
:
a
-
a
2
=
μ
+
2
α
a
1
-
μ
-
μ
+
2
α
a
1
-
μ
2
,
α
a
∈
0,1
.FormulaC
2
=
A
(
1
-
A
) can be determined by(27)
C
2
=
A
1
-
A
:
a
1
-
a
=
μ
+
2
α
a
1
-
μ
·
1
-
μ
+
2
α
a
1
-
μ
=
μ
+
2
α
a
1
-
μ
-
μ
+
2
α
a
1
-
μ
2
=
a
-
a
2
,
α
a
∈
0,1
.As can be seen, result for the second formC
2 is the same as for C
1. Result for the third form C
3 can be determined by(28)
C
3
=
A
-
1
+
1
-
A
1
+
A
:
a
-
1
+
a
-
1
1
+
a
=
μ
+
2
α
a
1
-
μ
-
1
+
1
-
μ
+
2
α
a
1
-
μ
·
1
-
μ
+
2
α
a
1
-
μ
=
μ
+
2
α
a
1
-
μ
-
1
+
1
-
μ
+
2
α
a
1
-
μ
2
=
μ
+
2
α
a
1
-
μ
-
μ
+
2
α
a
1
-
μ
2
=
a
-
a
2
,
α
a
∈
0,1
.The resultC
3 is identical to C
1 and C
2. As the analysed example shows the multidimensional RDM FA provides unique solutions for nonlinear formulas independently of their mathematical form, it cannot be said about μ-cuts’ FA. To better realise this fact, particular solutions will be visualised in Figure 4.One can easily see that all three solutions,C
1, C
2, and C
3, are incorrect. According to the solution C
1, it is possible that the difference c
=
a
-
a
2 for a
∈
[
0,2
] could be equal to, for example, -
4 or +
2. However, these values cannot be obtained with any value of a
∈
[
0,2
]. According to the “result” C
2, possible values of c
=
a
-
a
2 could be equal to +
1 or +
2. But, in fact, these values are impossible. According to the “result” C
3, possible values of c are equal to -
4 or +
4. However, it is not true.The full and precise 3D solution, in which formula (obtained with use of the RDM FA) is given by(29)
C
:
c
=
μ
-
μ
2
+
2
α
a
1
-
μ
-
4
α
a
1
-
μ
-
4
α
a
2
1
-
μ
2
,is shown in Figure 5.Figure 5
Full and complete 3D result of formulaC
=
A
-
A
2, where A
=
[
μ
,
2
-
μ
], obtained with use of the multidimensional RDM FA (29), with marked upper (right) and lower (left) constraint of 2D membership function.If we are interested in the simplified, traditional MF of the full 3D result in the horizontal version, it can be easily determined analytically on the basis of formula (30) for particular membership levels (μ-cuts). Consider(30)
C
=
A
-
A
2
=
c
L
μ
,
c
R
μ
=
min
α
a
⁡
c
μ
,
α
a
,
max
α
a
⁡
c
μ
,
α
a
,
μ
,
α
a
∈
0,1
,where c
(
μ
,
α
a
) is determined by formula (29). Extremes of c
(
μ
,
α
a
) can be detected by usual, analytical function examination. However, one should remember that function extremes can lie not only on domain borders of variable α
a (α
a
=
0, α
a
=
1) but also inside of the domain (fractional values of α
a). Then, the extremes can be detected by identification of critical points of the derivative d
c
(
μ
,
α
a
)
/
d
α
a
=
0.In the analysed problem, examination of the functionc
(
μ
,
α
a
) has shown that its left border c
L
(
μ
) and right border c
R
(
μ
) have mathematical form expressed by the following formula:(31)
C
=
A
-
A
2
=
c
L
μ
,
c
R
μ
,
c
L
=
c
μ
,
α
a
=
1
=
-
2
+
3
μ
-
μ
2
,
c
R
=
0.25
for
μ
∈
0,0.5
c
μ
,
α
a
=
0
=
μ
-
μ
2
for
μ
∈
0.5,1
.Membership function representing the span (width) of the full resultc
(
μ
,
α
a
) of the fuzzy formula C
=
A
-
A
2 is shown in Figure 6.Figure 6
Membership function representing the span of the full 3D result of the fuzzy formulaC
=
A
-
A
2.Comparison of MFs of results obtained with use of theμ-cuts’ FA (Figure 4) and with use of the multidimensional RDM FA (Figure 6) shows that “results” provided by the popular μ-cuts’ FA are excessively broad and false, because they suggest that certain values of the result c
=
a
-
a
2 are possible (e.g., -
4, +
2, +
4, and others), whereas they in fact are impossible. It can be easily checked analytically or numerically for various allowed values of a
∈
[
0,2
].Similarly, as the fuzzy formulaC
=
A
-
A
2 has been solved with use of the multidimensional RDM FA, other linear and nonlinear formulas or equations can be solved also. Their solutions usually will be less uncertain than solutions obtained with not-multidimensional methods.


## 2.1. Additive Inverse Element-
X of Element X

Inμ-cuts’ FA, trapezoidal membership function can be defined as a fuzzy interval with left border x
L
(
μ
)
=
a
+
(
b
-
a
)
μ and right border x
R
(
μ
)
=
d
-
(
d
-
c
)
μ (formula (13) and Figure 2). Consider(13)
X
μ
=
x
L
μ
,
x
R
μ
.A negative element-
X is determined by the following formula:(14)
-
X
μ
=
-
x
L
μ
,
x
R
μ
=
-
x
R
μ
,
-
x
L
μ
.The subtraction resultX
(
μ
)
-
X
(
μ
) determined by (15) is a fuzzy interval of nonzero span which means that, in μ-cuts’ FA additive, inverse element does not exist. Consider(15)
X
μ
-
X
μ
=
-
x
R
μ
-
x
L
μ
,
x
R
μ
-
x
L
μ
≠
0,0
.In the case of RDM FA, the elementX
(
μ
,
α
x
) is determined by (16) and its negative element -
X
(
μ
,
α
x
) by (17). Consider(16)
X
μ
,
α
x
=
x
μ
,
α
x
:
x
μ
,
α
x
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
,
α
x
∈
0,1
,
(17)
-
X
μ
,
α
x
=
x
μ
,
α
x
:
x
μ
,
α
x
=
-
x
L
μ
-
α
x
x
R
μ
-
x
L
μ
,
α
x
∈
0,1
.It is easy to check that the subtraction resultX
-
X is exactly equal to zero:(18)
X
μ
,
α
x
-
X
μ
,
α
x
=
0
.Thus, in the case of RDM FA, there exists the inverse additive element-
X
(
μ
,
α
x
).


## 2.2. Distributive Law of Multiplication

In theμ-cuts’ FA, similarly as in Moore’s interval arithmetic, the distributive law holds only in the limited form:(19)
X
μ
Y
μ
+
Z
ν
⊆
X
μ
Y
μ
+
X
μ
Z
μ
.In the multidimensional, continuous RDM FA, the distributive law holds fully. Consider(20)
X
μ
,
α
x
Y
μ
,
α
y
+
Z
μ
,
α
z
=
X
μ
,
α
x
Y
μ
,
α
y
+
X
μ
,
α
x
Z
μ
,
α
z
,
α
x
,
α
y
,
α
z
∈
0,1
.Proof.
For any three intervals expressed in RDM notation, withα
x
,
α
y
,
α
z
∈
[
0,1
],(21)
X
μ
,
α
x
=
x
L
μ
,
α
x
,
x
R
μ
,
α
x
=
x
μ
,
α
x
:
x
μ
,
α
x
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
,
Y
μ
,
α
y
=
y
L
μ
,
α
y
,
y
R
μ
,
α
y
=
y
μ
,
α
y
:
y
μ
,
α
y
=
y
L
μ
+
α
y
y
R
μ
-
y
L
μ
,
Z
μ
,
α
z
=
z
L
μ
,
α
z
,
z
R
μ
,
α
z
=
z
μ
,
α
z
:
z
μ
,
α
z
=
z
L
μ
+
α
z
z
R
μ
-
z
L
μ
,we have(22)
X
μ
,
α
x
Y
μ
,
α
y
+
Z
μ
,
α
z
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
·
y
L
μ
+
α
y
y
R
μ
-
y
L
μ
+
z
L
μ
+
α
z
z
R
μ
-
z
L
μ
=
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
y
L
μ
+
α
y
y
R
μ
-
y
L
μ
+
x
L
μ
+
α
x
x
R
μ
-
x
L
μ
z
L
μ
+
α
z
z
R
μ
-
z
L
μ
=
X
μ
,
α
x
Y
μ
,
α
y
+
X
μ
,
α
x
Z
μ
,
α
y
.


## 2.3. Different Forms of Nonlinear Formulas

Let us consider results ofμ-cuts’ FA and RDM FA for the nonlinear formula C
=
A
-
A
2, where A is the fuzzy interval A
=
[
0,1
,
2
] (Figure 4). This formula can take at least 3 forms: C
1
=
A
-
A
2, C
2
=
A
(
1
-
A
), and C
3
=
(
A
-
1
)
+
(
1
-
A
)
(
1
+
A
). Fuzzy interval A has the form A
=
[
a
L
(
μ
)
,
a
R
(
μ
)
]
=
[
μ
,
2
-
μ
]. Calculating value of C for particular forms, C
1, C
2, and C
3, with the application of μ-cuts’ FA, three different results are obtained:(23)
C
1
=
-
4
+
3
μ
-
μ
2
,
2
-
μ
-
μ
2
,
with  support  span
-
4,2
,
C
2
=
-
2
+
3
μ
-
μ
2
,
2
-
3
μ
+
μ
2
,
with  support  span
-
2,2
,
C
3
=
-
4
+
5
μ
-
μ
2
,
4
-
5
μ
+
μ
2
,
with  support  span
-
4,4
.Each of the results,C
1, C
2, and C
3, is different. For μ
=
0, results mean spans (widths) of MF supports. They are equal to [
-
4,2
], [
-
2,2
], and [
-
4,4
]. The differences are considerable. Thus, the formula C
=
A
-
A
2 cannot be calculated uniquely with μ-cuts’ FA. If we use the multidimensional RDM FA, then interval A is expressed by (24) and its square A
2 by (25), where a is generic variable of A. One should notice that functions expressing A and A
2 are fully coupled (correlated); that is, values of a and α
a have to be identical both in A and in A
2. Consider(24)
A
:
a
=
a
L
+
α
a
a
R
-
a
L
=
μ
+
2
α
a
1
-
μ
,
(25)
A
2
:
a
2
=
μ
+
2
α
a
1
-
μ
2
,
α
a
∈
0,1
.FormulaC
1
=
A
-
A
2 can be determined by(26)
C
1
=
A
-
A
2
:
a
-
a
2
=
μ
+
2
α
a
1
-
μ
-
μ
+
2
α
a
1
-
μ
2
,
α
a
∈
0,1
.FormulaC
2
=
A
(
1
-
A
) can be determined by(27)
C
2
=
A
1
-
A
:
a
1
-
a
=
μ
+
2
α
a
1
-
μ
·
1
-
μ
+
2
α
a
1
-
μ
=
μ
+
2
α
a
1
-
μ
-
μ
+
2
α
a
1
-
μ
2
=
a
-
a
2
,
α
a
∈
0,1
.As can be seen, result for the second formC
2 is the same as for C
1. Result for the third form C
3 can be determined by(28)
C
3
=
A
-
1
+
1
-
A
1
+
A
:
a
-
1
+
a
-
1
1
+
a
=
μ
+
2
α
a
1
-
μ
-
1
+
1
-
μ
+
2
α
a
1
-
μ
·
1
-
μ
+
2
α
a
1
-
μ
=
μ
+
2
α
a
1
-
μ
-
1
+
1
-
μ
+
2
α
a
1
-
μ
2
=
μ
+
2
α
a
1
-
μ
-
μ
+
2
α
a
1
-
μ
2
=
a
-
a
2
,
α
a
∈
0,1
.The resultC
3 is identical to C
1 and C
2. As the analysed example shows the multidimensional RDM FA provides unique solutions for nonlinear formulas independently of their mathematical form, it cannot be said about μ-cuts’ FA. To better realise this fact, particular solutions will be visualised in Figure 4.One can easily see that all three solutions,C
1, C
2, and C
3, are incorrect. According to the solution C
1, it is possible that the difference c
=
a
-
a
2 for a
∈
[
0,2
] could be equal to, for example, -
4 or +
2. However, these values cannot be obtained with any value of a
∈
[
0,2
]. According to the “result” C
2, possible values of c
=
a
-
a
2 could be equal to +
1 or +
2. But, in fact, these values are impossible. According to the “result” C
3, possible values of c are equal to -
4 or +
4. However, it is not true.The full and precise 3D solution, in which formula (obtained with use of the RDM FA) is given by(29)
C
:
c
=
μ
-
μ
2
+
2
α
a
1
-
μ
-
4
α
a
1
-
μ
-
4
α
a
2
1
-
μ
2
,is shown in Figure 5.Figure 5
Full and complete 3D result of formulaC
=
A
-
A
2, where A
=
[
μ
,
2
-
μ
], obtained with use of the multidimensional RDM FA (29), with marked upper (right) and lower (left) constraint of 2D membership function.If we are interested in the simplified, traditional MF of the full 3D result in the horizontal version, it can be easily determined analytically on the basis of formula (30) for particular membership levels (μ-cuts). Consider(30)
C
=
A
-
A
2
=
c
L
μ
,
c
R
μ
=
min
α
a
⁡
c
μ
,
α
a
,
max
α
a
⁡
c
μ
,
α
a
,
μ
,
α
a
∈
0,1
,where c
(
μ
,
α
a
) is determined by formula (29). Extremes of c
(
μ
,
α
a
) can be detected by usual, analytical function examination. However, one should remember that function extremes can lie not only on domain borders of variable α
a (α
a
=
0, α
a
=
1) but also inside of the domain (fractional values of α
a). Then, the extremes can be detected by identification of critical points of the derivative d
c
(
μ
,
α
a
)
/
d
α
a
=
0.In the analysed problem, examination of the functionc
(
μ
,
α
a
) has shown that its left border c
L
(
μ
) and right border c
R
(
μ
) have mathematical form expressed by the following formula:(31)
C
=
A
-
A
2
=
c
L
μ
,
c
R
μ
,
c
L
=
c
μ
,
α
a
=
1
=
-
2
+
3
μ
-
μ
2
,
c
R
=
0.25
for
μ
∈
0,0.5
c
μ
,
α
a
=
0
=
μ
-
μ
2
for
μ
∈
0.5,1
.Membership function representing the span (width) of the full resultc
(
μ
,
α
a
) of the fuzzy formula C
=
A
-
A
2 is shown in Figure 6.Figure 6
Membership function representing the span of the full 3D result of the fuzzy formulaC
=
A
-
A
2.Comparison of MFs of results obtained with use of theμ-cuts’ FA (Figure 4) and with use of the multidimensional RDM FA (Figure 6) shows that “results” provided by the popular μ-cuts’ FA are excessively broad and false, because they suggest that certain values of the result c
=
a
-
a
2 are possible (e.g., -
4, +
2, +
4, and others), whereas they in fact are impossible. It can be easily checked analytically or numerically for various allowed values of a
∈
[
0,2
].Similarly, as the fuzzy formulaC
=
A
-
A
2 has been solved with use of the multidimensional RDM FA, other linear and nonlinear formulas or equations can be solved also. Their solutions usually will be less uncertain than solutions obtained with not-multidimensional methods.


## 3. Addition of Two Independent Triangle Fuzzy Numbers

Triangle numbers (Figure7) are used in practice very frequently. Two numbers “about 1.0” and “about 1.1” can represent weights of two stones, S1 and S2, that were weighed on the scales. The scales have shown x
1
=
1.0 kg for stone 1 and x
2
=
1.1 kg for stone 2. Because the scales have the maximal error ±
0.1 kg, stone weights are uncertain and can be described by fuzzy numbers: X
1 is [0.9, 1.0, 1.1] and X
2 is [1.0, 1.1, 1.2].Figure 7
Membership functions of uncertain weightsX
1 and X
2 of two stones S1 and S2.
(a)
(b)Horizontal MF of the stone S1 is given by (32) and horizontal MF of the stone S2 is given by (33). Consider(32)
x
1
=
a
1
+
μ
b
1
-
a
1
+
α
x
1
d
1
-
a
1
1
-
μ
,
x
1
=
0.9
+
0.1
μ
+
0.2
α
x
1
1
-
μ
,
α
x
1
∈
0,1
,
(33)
x
2
=
a
2
+
μ
b
2
-
a
2
+
α
x
2
d
2
-
a
2
1
-
μ
,
x
2
=
1.0
+
0.1
μ
+
0.2
α
x
2
1
-
μ
,
α
x
2
∈
0,1
.If we want to calculate the sumy of weights of both stones, then the calculation with horizontal MFs can be made similarly as in a classical arithmetic, without use of the extension principle of Zadeh [37]. It is necessary only when vertical MFs are used. According to this principle, the membership function of the addition result is determined by (34), where A
1 and A
2 are fuzzy numbers [4, 10]. Consider(34)
μ
A
1
+
A
2
y
=
⋁
y
=
x
1
+
x
2
μ
A
1
x
1
∧
μ
A
2
x
2
,
∀
x
1
,
x
2
,
y
∈
R
.If horizontal MFs are used, then the sumy
=
x
1
+
x
2 is created directly by adding x
1 and x
2 determined by (32) and (33). Consider(35)
x
1
∈
a
1
,
d
1
,
x
2
∈
a
2
,
d
2
,
y
=
x
1
+
x
2
=
a
1
+
a
2
+
μ
b
1
+
b
2
-
a
1
-
a
2
+
α
x
1
d
1
-
a
1
+
α
x
2
d
2
-
a
2
1
-
μ
,
y
=
1.9
+
0.2
μ
+
0.2
α
x
1
+
α
x
2
1
-
μ
,
α
x
1
,
α
x
2
,
μ
∈
0,1
.Formula (35) shows that the sum y
=
x
1
+
x
2 is not 1-dimensional. It is a function defined in the 4D space because y
=
f
(
μ
,
α
x
1
,
α
x
2
). Thus, it cannot be visualised but it can be shown as the projection onto the 3D space: α
x
1
×
α
x
2
×
Y (Figure 8). The fourth dimension μ is shown in a form of μ-cuts. For μ
=
0, sum (35) takes a form: y
=
1.9
+
0.2
(
α
x
1
+
α
x
2
). For μ
=
1, the sum has value y
=
2.1 independently of values of α
x
1 and α
x
2. Surfaces of both μ-cuts can be seen in Figure 8.Figure 8
Multidimensional granule of the addition result of two fuzzy numbers “about 1.0” and “about 1.1” being the set of quadruples{
μ
,
α
x
1
,
α
x
2
,
y
} determined by formula (35) and projected onto the 3D space α
x
1
×
α
x
2
×
Y.The multidimensional result granule (Figure8) contains an infinitive number of μ-cuts corresponding to particular fractional values of μ
∈
[
0,1
]. Figure 8 shows only two border cuts for μ
=
0 and μ
=
1. Apart from the 3D projection shown in Figure 8, other projections of the full 4D granule are also possible. Figure 9 shows the projection onto the 3D space X
1
×
X
2
×
μ. In this figure, addition results y
=
x
1
+
x
2 are shown in a form of contour lines y
=
x
1
+
x
2
=
c
o
n
s
t.Figure 9
Projection of the full 4D addition result granule onto the 3D spaceX
1
×
X
2
×
μ with slant contour lines corresponding to particular addition results y
=
x
1
+
x
2
=
c
o
n
s
t.Each of the contour linesy
=
c
o
n
s
t is the set of infinite number of tuples {
x
1
,
x
2
} satisfying the condition x
1
+
x
2
=
y. As can be seen in Figure 9, sum y
=
2.3 can occur only for x
1
=
1.1 and x
2
=
1.2. Instead, sum y
=
2.1 has considerably greater possibility of occurrence because the number of tuples x
1
,
x
2
∣
x
1
+
x
2
=
y
=
2.1 is infinitely large. The value y
=
2.1 occurs, for example, for x
1
=
1.0 and x
2
=
1.1 and for x
1
=
1.01 and x
2
=
1.09. A cardinality of particular sets Y
y [69] can be calculated with the following formula:(36)
Abs  Card
y
=
Card
y
=
∫
C
min
⁡
μ
x
1
,
μ
x
2
d
s
,where C is a contour line defined as(37)
C
=
x
1
,
x
2
∣
x
1
+
x
2
=
y
,
x
1
∈
a
1
,
d
1
,
x
2
∈
a
2
,
d
2
.The cardinality of particular setsY
y is equal to an area of y-cuts of the granule (Figure 9). It is easy to calculate. For example, the cardinality of the set Y
2.1 equals(38)
Card
y
=
2.1
=
1
2
d
2
-
a
2
2
+
d
1
-
a
1
2
≈
0.14142
.It is the greatest of all cardinalities (Max Abs Card(y)). It can be shown that the cardinality of particular sets is a square function of the sum y
=
x
1
+
x
2 ((39) and (40)).Fory
<
a
1
+
d
2,(39)
Abs  Card
y
=
Card
y
=
y
-
a
1
+
a
2
a
2
+
d
1
-
a
1
+
a
2
2
·
Max  Abs  Card
y
.Fory
≥
a
1
+
d
2,(40)
Abs  Card
y
=
Card
y
=
d
1
+
d
2
-
y
d
1
+
d
2
-
d
1
+
a
2
2
·
Max  Abs  Card
y
.The absolute cardinality of the sumy
=
x
1
+
x
2 can take various numeric values. Therefore, it should be normalised to interval [
0,1
] for convenience ((41) and (42)). The normalised cardinality N
o
r
m
C
a
r
d
(
y
) can also be called the relative cardinality Rel Card(
y
).Fory
<
a
1
+
d
2,(41)
Norm  Card
y
=
Rel  Card
y
=
Abs  Card
y
Max  Abs  Card
y
=
y
-
a
1
+
a
2
d
1
-
a
1
2
.Fory
≥
a
1
+
d
2,(42)
Norm  Card
y
=
Rel  Card
y
=
Abs  Card
y
Max  Abs  Card
y
=
d
1
+
d
2
-
y
d
2
-
a
2
2
.It should be noticed that formulas (41) and (42) concern fuzzy numbers with an identical support. In the considered example, d
1
-
a
1
=
d
2
-
a
2
=
0.2, so they take the following form:(43)
F
o
r
y
<
a
1
+
d
2
:  Norm  Card
y
=
y
-
1.9
0.2
2
,
F
o
r
y
≥
a
1
+
d
2
:  Norm  Card
y
=
2.3
-
y
0.2
2
.The cardinality distribution of particular values of the sumy
=
x
1
+
x
2 for the considered example is shown in Figure 10.Figure 10
Normalised cardinality distributionN
o
r
m
C
a
r
d
(
y
) of the occurrence of particular values of the sum y
=
x
1
+
x
2 being the 2D representation of the 4D addition result of two independent fuzzy numbers “about 1.0” and “about 1.1” shown in Figure 7.It should be emphasised once more that the distribution determined by formula (43) and shown in Figure 10 is not the addition result of two fuzzy numbers but only 2D information about the cardinality of particular sum values y
=
x
1
+
x
2 derived from the full 4D addition result given by (35). The representation provided by Zadeh’s extension principle (34) can also be obtained on the basis of the exact solution (35). The full 4D solution is given below. Consider(44)
y
=
1.9
+
0.2
μ
+
0.2
α
x
1
+
α
x
2
1
-
μ
α
x
1
,
α
x
2
,
μ
∈
0,1
.The minimal value of the sum (the left border)y
L for various levels of μ
∈
[
0,1
] is obtained for α
x
1
=
α
x
2
=
0:(45)
y
L
=
1.9
+
0.2
μ
.The maximal sum value (the right border)y
R is obtained for α
x
1
=
α
x
2
=
1:(46)
y
R
=
2.3
-
0.2
μ
.The 2D representation of the addition result of two fuzzy numbers “about 1.0” and “about 1.1” obtained with formulae (45) and (46) is shown in Figure 11.Figure 11
2D representation of the full 4D addition result obtained with formulas (45) and (46) corresponding to the “result” obtained with Zadeh’s extension principle (34).The 2D representation obtained according to Zadeh’s extension principle (34) informs only about the maximal spread of the full solution granule on particular μ-levels (Figure 9). Geometrically, the distribution from Figure 11 corresponds to the cross section going from the corner y
=
a
1
+
a
2
=
1.9 to the corner y
=
d
1
+
d
2
=
2.3 through the peak of the membership function in Figure 9.


## 4. Addition of Two Fuzzy Numbers with Taking into Account the Order Relation

Let us analyse an example of a fuzzy number addition similar to that in Section3, but with small difference, with the additional knowledge that not only MFs of added numbers but also their order relation is known.Thus, there are two stones S1 and S2 which had been weighed on spring scales with maximal error equal to 0.1 kg. Therefore, weightsx
1 and x
2 of stones are uncertain and can be expressed in the form of fuzzy numbers: X
1 is [
0.9,1.0,1.1
] and X
2 is [
1.0,1.1,1.2
]. MFs of uncertain weights are shown in Figure 7. After weighting the stones on the spring scales, their weights were compared on balance scales. The scales showed that the weight x
2, though uncertain, is greater than x
1. Therefore, the weights order x
2
>
x
1 is known and can be taken into account in the weights’ adding. The knowledge of the order relation changes (constraints) the domain of possible weight tuples {
x
1
,
x
2
} in comparison with the example from Section 3. Now less tuples will be feasible (Figure 12).Figure 12
Domain of feasible tuples of weight values{
x
1
,
x
2
} which can take part in adding of fuzzy numbers x
1 and x
2.In the considered example, uncertain weightsx
1 and x
2 are expressed by the same horizontal MFs as given by formulas (32) and (33). However, the addition result y is determined by more constraints (formula (47)) as in the case without the weights’ order relation. Consider(47)
x
1
∈
a
1
,
d
1
,
x
2
∈
a
2
,
d
2
,
x
1
<
x
2
,
y
=
x
1
+
x
2
=
a
1
+
a
2
+
μ
b
1
+
b
2
-
a
1
-
a
2
+
α
x
1
d
1
-
a
1
+
α
x
2
d
2
-
a
2
1
-
μ
,
y
=
1.9
+
0.2
μ
+
0.2
α
x
1
+
α
x
2
1
-
μ
,
α
x
1
,
α
x
2
,
μ
∈
0,1
.Equations (47) determine the full and precise addition result y
=
f
(
α
x
1
,
α
x
2
,
μ
) that has a form of the 4D information granule. This granule cannot be visualised. However, we can visualise its 3D projection onto the space X
1
×
X
2
×
Y that is shown in Figure 13.Figure 13
3D projection of the 4D granule of the addition result of two fuzzy numbers “about 1.0” and “about 1.1” onto the spaceX
1
×
X
2
×
Y.A comparison of addition result granules obtained without (Figure8) and with (Figure 13) taking into account the order relation x
2
>
x
1 shows that it considerably changes the granules shape. Figure 14 shows other projection of the full 4D result granule onto the space X
1
×
X
2
×
μ. Values of the sum y
=
x
1
+
x
2 are shown in this figure in a form of contour lines of constant y values.Figure 14
Projection of the full 4D result granule of the addition of two fuzzy numbers “about 1.0” and “about 1.1” with taking into account the order relationx
2
>
x
1 onto the 3D space X
1
×
X
2
×
μ.The granule of the addition result from Figure14 can be compared to the result granule of addition without the order relation x
2
>
x
1 shown in Figure 9. The comparison shows that the relation x
2
>
x
1 considerably decreases the domain of possible solutions. At the same time, we can observe decrease of cardinalities of particular y values for y
∈
[
2.0,2.2
]. Now, formulas for cardinalities take new forms: (50)–(56). Consider(48)
Max  Abs  Card
=
1
2
·
7
8
·
d
2
-
a
2
2
+
d
1
-
a
1
2
.Fory
∈
[
a
1
+
a
2
,
a
1
+
b
2
]
=
[
1.9,2.0
],(49)
Abs  Card
y
=
8
7
y
-
a
1
+
a
2
d
1
-
a
1
2
·
Max  Abs  Card
,
(50)
Norm  Card
y
=
8
7
y
-
a
1
+
a
2
d
1
-
a
1
2
=
8
7
y
-
1.9
0.2
2
.Fory
∈
[
a
1
+
b
2
,
a
1
+
d
2
]
=
[
2.0,2.1
],(51)
Abs  Card
y
=
8
7
y
-
a
1
+
a
2
d
1
-
a
1
2
-
1
8
y
-
a
1
+
b
2
d
1
-
b
1
2
·
Max  Abs  Card
,
(52)
Norm  Card
y
=
8
7
y
-
a
1
+
a
2
d
1
-
a
1
2
-
1
8
y
-
a
1
+
b
2
d
1
-
b
1
2
=
8
7
y
-
1.9
0.2
2
-
1
8
y
-
2.0
0.1
2
.Fory
∈
[
a
1
+
d
2
,
b
1
+
d
2
]
=
[
2.1,2.2
],(53)
Abs  Card
y
=
8
7
d
1
+
d
2
-
y
d
2
-
a
2
2
-
1
8
b
1
+
d
2
-
y
b
1
-
a
1
2
·
Max  Abs  Card
,
(54)
Norm  Card
y
=
8
7
d
1
+
d
2
-
y
d
2
-
a
2
2
-
1
8
b
1
+
d
2
-
y
b
1
-
a
1
2
=
8
7
2.3
-
y
0.2
2
-
1
8
2.2
-
y
0.1
2
.Fory
∈
[
b
1
+
d
2
,
d
1
+
d
2
]
=
[
2.2,2.3
],(55)
Abs  Card
y
=
8
7
d
1
+
d
2
-
y
d
1
-
a
1
2
·
Max  Abs  Card
,
(56)
Norm  Card
y
=
8
7
d
1
+
d
2
-
y
d
1
-
a
1
2
=
8
7
2.3
-
y
0.2
2
.Figures15(a) and 15(b) show the comparison of 2D representations of full 4D addition results in a form of cardinality distributions N
o
r
m
C
a
r
d
(
y
) determined by (50)–(56).Figure 15
2D representation of the full 4D result in a form of the cardinality distribution for the case of the addition of two independent fuzzy numbers (a) and for the case of the addition of two fuzzy numbers constrained by the order relationx
2
>
x
1 (b).
(a)
(b)It should be noticed that the 2D representation obtained with use of Zadeh’s extension principle (34) is still identical and has the form shown in Figure 16 independently whether the addition is realised with or without taking into account the order relation x
2
>
x
1.Figure 16
2D representation of the full 4D addition result obtained with use of the extension principle (34) both for the case of two independent numbers and dependent numbers constrained by the order relation x
2
>
x
1.The extension principle (34) does not “perceive” the domain loss (compare Figures 8 and 13) which is its important drawback. It provides even less information about the full 4D addition result than the 2D representation determined with formulas (50)–(56) and shown in Figure 15.


## 5. Addition of Two Fully Dependent Fuzzy Numbers

Correlations between variables have no significance in the case of an addition of two crisp numbers. However, in the case of intervals, fuzzy intervals, and fuzzy numbers, mutual dependencies are of great importance. Let us consider another example.Small John did not know exactly how much money he had in his money box after a month of saving. He evaluated that the sum was about $10, at least $7, and no more than $13. This evaluation can be expressed by the fuzzy number [46, 49, 51]. His father promised to double the saved sum after 1 month. He opened the box, checked the sum, and added exactly the same sum to the box. However, father did not inform John of the sum he added. How much money does small John have in the box now?Let us denote byx
1 the uncertain sum of money which small John had in the money box at the beginning: x
1 is “about 10” [46, 49, 51]. The horizontal MF of the set “about 10” can be determined on the base of formula (32):(57)
x
1
=
7
+
3
μ
+
6
α
x
1
1
-
μ
,
μ
,
α
x
1
∈
0,1
.Though the precise value of the sumx
1 is not known for us and for John, his father knew it and added the same sum x
2 to John’s money box. Because x
1
=
x
2, then from (58) we can conclude that α
x
1
=
α
x
2. Consider(58)
x
1
=
7
+
3
μ
+
6
α
x
1
1
-
μ
=
x
2
=
7
+
3
μ
+
6
α
x
2
1
-
μ
,
α
x
1
,
α
x
2
,
μ
∈
0,1
.It means that both RDM variables and also variablesx
1 and x
2 are mutually fully coupled (correlated). The resulting sum y
=
x
1
+
x
2 can be expressed by the following formula:(59)
y
=
x
1
+
x
2
=
14
+
6
μ
+
6
1
-
μ
α
x
1
+
α
x
2
,
α
x
1
=
α
x
2
α
x
1
,
α
x
2
,
μ
∈
0,1
,and it is visualised in Figure 17.Figure 17
Visualization of the addition result of two fully dependent fuzzy numbers: “about 10” + “about 10” for which the conditionα
x
1
=
α
x
2 has to be satisfied.As Figure17 shows, the addition result “about 20” is in this case the flat and not 3D information granule as in the case of the addition of fully independent fuzzy numbers; see Figure 9. As before, the flat 2D representation in the form of the horizontal MF y
=
f
(
μ
) can be determined on the basis of exact solution (59). This function can be characterized with left and right border functions y
L
(
μ
) and y
R
(
μ
). Both functions can be determined from the full result:(60)
y
=
14
+
6
μ
+
12
α
x
1
1
-
μ
.The left, minimal function border is obtained for minimal values of RDM variablesα
x
1
=
α
x
2
=
0. Consider(61)
y
L
μ
=
14
+
6
μ
.The right border of the 2D representation, the border of maximaly values, is obtained for maximal values of RDM variables α
x
1
=
α
x
2
=
1. Consider(62)
y
R
μ
=
26
-
6
μ
.2D representation of the full 4D addition result (59) is shown in Figure 18.Figure 18
2D representation of the full 4D result (59) y
=
f
(
μ
,
α
x
1
,
α
x
2
) of the addition of two fully dependent fuzzy numbers “about 10” + “about 10” satisfying the condition α
x
1
=
α
x
2.


## 6. Solving Additive Fuzzy Equation:A
+
X
=
C

As it was mentioned in Section1, solving equations with one unknown x makes great difficulties even for the interval arithmetic which is considerably simpler than the fuzzy arithmetic. It will be shown in this section that solving fuzzy equations with use of horizontal MFs is possible and not difficult. Let us consider the last example.Corn harvested from field 1 was weighed on the scales that have maximal error equal to 1 ton. The scales indicated 5 tons. Thus, the true weighta of the corn can be expressed in a form of the fuzzy number A
=
[
4,5
,
6
]. Corn harvested from field 2 was not weighed on the field. Therefore, its true weight x was not known even approximately. Both crops were brought to a warehouse. There, they were mixed and weighed together on 3 scales that also have the maximal error equal to 1 ton. The scales indicated 15 tons. Thus, the sum c of crops can be expressed in a form of the fuzzy number C
=
[
14,15,16
]. How large is the weight x of the crop from field 2? The knowledge of a and x is important, because farmers 1 and 2 should be paid fairly for their delivery.This problem cannot be solved uniquely withμ-cut method based on the traditional interval arithmetic. However, it can be solved with an application of horizontal MFs (63) which can be constructed on the basis of the general formula (12). Consider(63)
A
=
4,5
,
6
:
a
=
4
+
μ
+
2
α
a
1
-
μ
,
C
=
14,15,16
:
c
=
14
+
μ
+
2
α
c
1
-
μ
,
μ
,
α
a
,
α
c
∈
0,1
.Valuec, though known only approximately, is the result of the addition of also only approximately known values a and x. They are connected by a dependence a
+
x
=
c. Thus, the unknown value x of crop from field 2 can be calculated from the following formula: (64)
x
=
c
-
a
=
14
+
μ
+
2
α
c
1
-
μ
-
4
+
μ
+
2
α
a
1
-
μ
=
10
+
2
α
c
-
α
a
1
-
μ
μ
,
α
a
,
α
c
∈
0,1
.It can be easily noticed that formula (64) is four-dimensional: x
=
f
(
μ
,
α
a
,
α
c
). It cannot be expressed exactly in the 2D space, which is suggested by the extension principle and various versions of the fuzzy arithmetic. However, it can be visualised in a form of the 3D projection on subspaces, for example, on the subspace A
×
X
×
μ (variable α
a is equivalent to variable a because of the transformation a
=
4
+
2
α
a) (Figure 19).Figure 19
3D projection of the full 4D solution of (64) from the space A
×
C
×
X
×
μ onto the space A
×
X
×
μ.Figure19 shows the support of the membership function of solution (64) being the μ-cut of this function on the level μ
=
0. It can be noticed that this support does not have a shape of a rectangle but of a parallelogram, which is shown additionally in Figure 20.Figure 20
Support of the 3D membership function of the solution from Figure19. TP1 = (
5,9.5
) and TP2 = (
4.5,9
) are examples of test points that can be used for the correctness testing of obtained solutions.Figure20 shows that the solution of the fuzzy equation A
+
X
=
C is not 1-dimensional one and therefore it is not possible to express it in a form of the interval [
x
_
,
x
¯
] which is usually done by these versions of fuzzy arithmetics that are based on the classical IA. Moore’s IA provides us with two possible solutions: (65) and (66). The first result (65) is obtained if the equation A
+
X
=
C is solved for the cut μ
=
0. Consider(65)
A
+
X
=
C
,
4,6
+
x
_
,
x
¯
=
14,16
,
4
+
x
_
=
14
,
x
_
=
10
,
6
+
x
¯
=
16
,
x
¯
=
10
,
x
_
,
x
¯
=
10,10
.Solution (65) represents only one of many possible solution subsets: it is the subset containing tuples a
,
x
∣
x
=
10. However, apart from this subset, there is also an infinitely large set of solution tuples {
a
,
x
} which satisfy the interval equation A
+
X
=
C. Example of such tuple is the test point TP1 (
a
,
x
)
=
(
5.0,9.5
) in Figure 20, because, for a
=
5
∈
[
4,6
], the sum a
+
x
=
c
=
5
+
9.5
=
14.5
∈
[
14,16
].The second possible solution of the equationA
+
X
=
C for μ
=
0 gives formula (66). Consider(66)
A
+
X
=
C
,
4,6
+
x
_
,
x
¯
=
14,16
,
x
_
,
x
¯
=
14,16
-
4,6
=
8,12
.One can easily convince himself that this result is not the precise solution of the considered equation. It is only information about the greatest span[
8,12
] of the precise solution (64), which has a form of a parallelogram (Figure 20). As can be seen in Figure 20, for the interval [
x
_
,
x
¯
]
=
[
8,12
], we can create an infinite number of tuples {
a
,
x
} which are not solutions of the considered equation. For example, the tuple (
4.5,9
) corresponding to the test point TP2 has coordinates a
=
4.5
∈
[
4,6
] and x
=
9
∈
[
8,12
]
=
[
x
_
,
x
¯
]. Above-mentioned examples show that only multidimensional approach to the fuzzy and interval arithmetic can give correct problem solutions. It should be noticed that the maximal solution span [
8,12
] provided by Moore’s arithmetic can be easily obtained from the multidimensional solution (64) of equation A
+
X
=
C. If we are interested in the solution span on the cut level μ
=
0, then this value of μ can be inserted into (64) being 4D solution and, as a result, we obtain(67)
x
=
c
-
a
=
10
+
2
α
c
-
α
a
,
α
a
,
α
c
∈
0,1
.A simple analysis of this equation gives the conclusion that the minimal valuex
_ occurs for α
a
=
1 and α
c
=
0. This value equals x
_
=
8. The maximal value x
¯ occurs for α
a
=
0 and α
c
=
1. It is equal to 12. Hence, the maximal solution span [
x
_
,
x
¯
]
=
[
8,12
]. Figure 21 shows 2D representation of the full 4D solution of the equation A
+
X
=
C in a form of the cardinality distribution of particular x values.Figure 21
2D distribution of the full 4D solution of the fuzzy equationA
+
X
=
C in a form of the cardinality distribution of particular x values, x
∈
[
8,12
].


## 7. Conclusions

The paper presents the new concept of horizontal MFs which were used for the addition as the exemplary operation on fuzzy numbers. The application of such functions in combination with the RDM arithmetic enables multidimensional approach to operations on fuzzy numbers and thereby it allows for removing some drawbacks and weaknesses of the fuzzy arithmetic. The use of horizontal MFs considerably facilitates calculations because now uncertain values can be inserted directly into equations without using the extension principle.Additionally, the RDM arithmetic enables taking into account correlations occurring between arguments of mathematical operations which is not possible using one-dimensional approach offered by classic FA methods.



---

*Source: 367214-2015-06-23.xml*