Fuzzy Logic An Introductory Course For Eng - Enric Trillas
Fuzzy Logic An Introductory Course For Eng - Enric Trillas
Fuzzy Logic An Introductory Course For Eng - Enric Trillas
Enric Trillas
Luka Eciolaza
Fuzzy Logic
An Introductory Course for Engineering
Students
Studies in Fuzziness and Soft Computing
Volume 320
Series editor
Janusz Kacprzyk, Polish Academy of Sciences, Warsaw, Poland
e-mail: [email protected]
About this Series
Fuzzy Logic
An Introductory Course for Engineering
Students
123
Enric Trillas Luka Eciolaza
European Centre for Soft Computing European Centre for Soft Computing
Mieres, Asturias Mieres, Asturias
Spain Spain
This book was thought as a non-conventional first course textbook in Fuzzy Logic
for engineers ending with an introduction to one of the most fruitful topics arisen
from it, Fuzzy Control. It is from the teaching’s strategy of the authors, summarized
by “Nothing can substitute the own homework of the student” from which it comes
its non-conventional character, partially manifested by the ‘continuous’ form of
presenting the considered topics by joining theoretical explanations and examples,
and not always following the typically mathematical style of ‘theorem-corolaries’.
Behind this strategy is the opinion that, at the university level, students and
professors ought to learn jointly, students do not wait to receive everything from the
professor’s lectures, but should read more than a single recommended textbook.
Consequently, this book is neither a manual with recipes to be uncritically applied,
nor it is directed to those that can be only interested in mathematical subtleties. The
reader should be aware that fuzzy logic is the study and computational management
of imprecision and non-random uncertainty, both with the highest accuracy and
precision possible at each case, that fuzzy logic is not fuzzy in itself.
Each university course requires a particular teaching tactic that not only depends
on the number of lecturing hours, but on the aim of the course and on the audience’s
characteristics. In particular, additional tutorials supplied by the professor are
essential for a good learning process. Tutorials in which other forms of considering
the course’s topics and more sophisticated problems can be proposed. This is at the
own hands of the professor.
The book just presents some basic mathematical models for fuzzy logic but
without the intention to just subordinate it to mathematics. Fuzzy logic is neither a
part of mathematics, nor even of logic, like Physics is not so. Notwithstanding,
what is paramount is the importance and usefulness of mathematical models in
experimental sciences and technology, as well as in computer science and computer
technology and, in particular, in Soft Computing, where fuzzy logic plays a pivotal
role. But the suitability of such models only can come from the success of its testing
against some reality, for instance, in true applications; applications play in the
techno-scientific world an analogous role to that of experimentation in natural
sciences. For instance, if the branch called ‘Fuzzy Control’ served not as a direct
vii
viii Preface
justification of fuzzy logic, the success fuzzy logic has in control applications can be
seen as a kind of experimentation to show its usefulness in the study of dynamical
systems linguistically described by systems of imprecise rules. Fuzzy logic is much
more than what is in this introductory textbook; its applications spread along many
domains of science and technology.
The reader should be early acquainted with the fact that, differently from clas-
sical bi-valuate logic, almost all in fuzzy logic are context-dependent and purpose-
driven. For instance, when representing a system of imprecise linguistic rules in
fuzzy terms, all the predicates, connectives, and conditionals in the rules should be
specified accordingly with its contextual meaning and the type of inference, for-
wards or backwards, to be done. In fuzzy logic not only everything is a matter of
degree, but its practice requires the art of designing systems by means of the
available theoretic armamentarium. The student should be conscious that a mistake
in the design process can conduct to solve a different problem than the targeted one.
It should never be forgotten that it is most important in research to not stop
questioning (Albert Einstein). Posing good questions (Isaac Rabi) whose answers
can result in fertile ones (Karl Menger) is what reveals the relevance of a researcher.
Contents
ix
x Contents
Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
Chapter 1
On the Roots of Fuzzy Sets
As fuzzy sets were introduced by Zadeh in 1965, they were born closely linked with
imprecise predicates, that is, with names of non-precisely defined classes of objects.
Even more, most of the applications of Zadeh’s ideas are made with properties
the objects do verify in some degree between the two classical extremes 0 and 1,
respectively. Because of that, it is not at all odd to introduce fuzzy sets from some
considerations on how predicates are used in language. We will follow Wittgenstein’s
statement “The meaning of a word is its use in language”.
Usually, isolated words mean nothing. To mean something, words are to be used
in a given context and in a known way. Words do serve to describe perceptions, to
translate reasoning, and to show the reasons for judgements.
For example, what it is meant by the predicate unleaty? It is impossible to answer
this question, since nobody has never heard something like “this is unleaty”, “such
is leaty”, etc. Neither unleaty, nor leaty, are English words, nobody has used them
and they don’t appear in English dictionaries. Meaning is inherited by predicates P
only after being used, in some ground, by means of elemental statements ‘x is P’.
By now, unleaty has no meaning.
Words are introduced in a language by using them in concrete ways. For example,
the Spanish word ‘madre’ comes from the latin ‘mater’ that, at its turn, came from one
in older indo-european languages, but always used to name someone’s mother. Later
on, the word could take, by analogy, new meanings as it is, for example, ‘mother
country’, ‘goodmother’, ‘mother in law’, ‘mother of all wars’, ‘mother Nature’,
‘mother-of-pearl’, etc. The predicate leaty neither appears in English, nor in Spanish,
French, German, …, because it was never before used to name a property of the
elements in some class, like it is with tall, young, middle-aged, with people, high
with buildings or mountains, or heavy with metals, for example. It is only through
its use that predicates do acquire ‘meaning’.
© Springer International Publishing Switzerland 2015 1
E. Trillas and L. Eciolaza, Fuzzy Logic, Studies in Fuzziness
and Soft Computing 320, DOI 10.1007/978-3-319-14203-6_1
2 1 On the Roots of Fuzzy Sets
x P y ⇔ x is less P than y,
and suppose P is a preorder (enjoys the reflexive and transitive properties). That is,
• x P x, for all x in X
• If x P y, and y P z, then x P z.
The preordered set (X, P ) reflects the organization P induces in X , and the
preorder P is the primary use of P in X . The relation −1P is defined by ‘x P
−1
y ⇔ y P x’. Note that the relation P is, usually, empirically perceived; it could
also be called the perceptive meaning of P in X .
We will say that ‘x is equally P than y’ whenever x P y and y P x, and write
x = P y, with (= P ) = ( P ∩ −1 P ). Since, obviously,
1.1 A Genesis of Fuzzy Sets 3
• x = P x, for all x in X
• x =P y ⇔ y =P x
• x = P y, and y = P z imply x = P z,
the relation = P is an equivalence in X , and gives the quotient set X/ = P , of the
classes of equally-P elements. The predicate P is semi-rigid in X if X/ = P consists
in a finite number of classes. Of course, all predicates on a finite X are semi-rigid.
Example 1.1.1 Let it be X = {x1 , . . . , x5 }, and P a predicate inducing the preorder
given by the matrix with entries
1, if xi P x j
entr y (i, j) =
0, otherwise,
that is, ⎛ ⎞
1 0 1 0 1
⎜0 1 0 1 0⎟
⎜ ⎟
[ P ] = ⎜
⎜1 0 1 0 1⎟⎟
⎝0 1 0 1 0⎠
1 0 1 0 1
The quotient set X/ = P has the two classes {x1 , x3 , x5 } and {x2 , x4 }.
If, in the same X , it is Q with primary use defined by
xi Q x j ⇔ i j,
that is, = P is the equality. Hence, X/ = P = {{x}; x ∈ [0, 10]}, and P is not semi-
rigid.
1.1.1 L-Degree
then μ P (x) μ P (y)- in the order of the poset. The idea behind this definition is
that μ P (x) ∈ L evaluates up to which extent x is P, up to which extent x verifies the
property named by P. It can be written,
x μ P y ⇔ μ P (x) μ P (y),
L-sets are ‘complex fuzzy sets’. Notice that (C, ) is isomorphic with the set of the
sub-intervals [a, b] ⊆ [0, 1], once partially ordered by
When it is only known that the degree up to which “x is P” belongs to some sub-
interval [a, b], it can be taken μ P (x) = a + ib. In science and technology, complex
quantities are not at all rare.
Remark 1.1.3 Provided L = [0, 1], the relation μ P is linear or total, since given
x, y in X , it is either μ P (x) μ P (y), or μ P (x) ≥ μ P (y); that is, either x μ P y,
or y μ P x.
Remark 1.1.4 The relation P is not always linear, there can exist elements x, y in X
such that it is neither x P y, nor y P x, that is, elements which are not comparable
under P . Of course, in the cases in which P is not linear and L = [0, 1], it cannot
be P = μ P , the degree cannot perfectly reflect the primary use of P, and the
relation μ P enlarges the primary use P with all the links between elements in X
contained in P − μ P , provided this difference-set is non empty.
1.1 A Genesis of Fuzzy Sets 5
Once the use ( P , μ P ) of P in X is given for the poset (L , ), it will be said that
the new object P defined by
∼
• x ∈r P , if and only if r = μ P (x), for r ∈ L,
∼
• P = Q , if and only if μ P = μ Q
∼ ∼
is the L-set labeled P. The set L X = {μ; μ : X → L}, is usually and abusively
called the set of all L-sets in X , since it contains all the possible degrees in L.
Notice that a predicate P could give many L-sets, in dependance on which poset
(L , ), and which function μ P ∈ L X , are chosen.
From now on, it will be supposed that (l, ) has a minimum element α (α
r, ∀r ∈ L), and a maximum element ω(r ω, ∀r ∈ L). With L 0 = {α, ω}, (L 0 , )
is a poset isomorphic to ({0, 1}, ).
Since a classical subset A of X is characterized by its membership function
1, if x∈A
μ A (x) =
0, if x∈
/ A,
the subset A can be viewed as an L-set, whose function only takes the values α or ω,
that is,
ω, if x ∈ A
μ A (x) =
α, if x ∈ / A.
Hence, the classical subsets of X are nothing else than the functions in L 0X , a set
included in L X . Classical subsets are limiting or particular case of L-sets. They are
degenerate L-sets.
Notice that if the classical subset A is labeled by a crisp predicate P, this predicate
is semi-rigid since X/ = P has, at most, two classes. Crisp predicates are rigid.
Very often μ P is perceptually designed, that is, designed from what is perceived
or known on the concrete use of P in X . It can be thought that such design is purely
subjective, in the sense of being made just by what the designer believes on the use
of P. But this would not be the case, except if the designer proceeds in a non rational
way. The designer should try to be as most sure as possible on the correction of his/
her/its perceptions.
6 1 On the Roots of Fuzzy Sets
In the second case, μ P has at most two values, and, for obvious reasons, we will
only consider the situation where this values vx , v y are different. When,
• X/ = P = {[x]}, and either μ P (x) = α for all x in X , or μ P (x) = ω for all x in
X , it results P = ∅ in the first case, and P = X in the second. In both cases, P
∼ ∼
is a rigid or binary predicate in X .
1.1 A Genesis of Fuzzy Sets 7
{x}
{y}
1 2 3 4 X
that is, P is the classical or crisp subset [x] of X , and P is a rigid or binary predicate
∼
in X .
Let us show some examples.
Example 1.1.7 With P = around five, and P as shown in the example [2 in 1.2],
the function
⎧
⎨ 0, if x ∈ [0, 4] ∪ [6, 10]
μ P (x) = x − 4, if x ∈ [4, 5]
⎩
6−x if x ∈ [5, 6],
whose graphics is
8 1 On the Roots of Fuzzy Sets
Example 1.1.8 Consider P = big in X = [0, 10]. Let us show several possible
degrees for it, after agreeing that ‘x P y if and only if x y, in the linear order of
[0, 10]’.
A [0, 1]-degree μ P is any function X = [0, 10] → [0, 1], such that
that is, any non-decreasing function (of which there are many). We can also agree
that μ P (0) = 0, and μ P (10) = 1. With this, it is clear that all degrees for big will
show some family resemblance.
Once fixed (L , ) = ([0, 1], ), and P = , the different uses of big only
depend on which function μ P is chosen to reflect the meaning of the predicate in
[0, 10]. Of course, it is P = μ P if and only if function μ P is strictly non-decreasing,
as it is the case either with μ P (x) = x/10, or with μ P (x) = x 2 /100. Provided μ P
is not strictly non-decreasing as, for example, with
since 6 P 4 (strictly), but μ P (6) = μ P (4) = 0.5, μ P does not perfectly reflect the
primary use of big in [0, 10]. In the same way, the crisp degree
1, if x > 8
μ P (x) =
0, otherwise,
does not perfectly reflect the primary use of big when translated into ‘after eight’.
Another possible model for μ P is
⎧
⎨ 0, if x ∈ [0, 2]
μ P (x) = 6 ,
x−2
if x ∈ [2, 8]
⎩
1 if x ∈ [8, 10],
with graphic.
1.1 A Genesis of Fuzzy Sets 9
All these models are linear, with the exception os μ P (x) = x 2 /100, that is
quadratic. Another quadratic models are given by μ(1)
P = ( 6 ) , and μ P =
x−2 2 (2)
(1)
1 − μ P (10 − x) = 1 − ( 8−x
6 ) , with graphics
2
with graphics.
10 1 On the Roots of Fuzzy Sets
Example 1.1.9 The predicate old, once numerically characterized by Ch old = Age,
can be translated into the interval [0, 120], in years, by μold (x) = μbig (Age(x)),
and a linear model for it can be
⎧
⎪
⎨ 0, if 0 Age(x) 45
μold (x) = μbig (Age(x)) = 30 ,
x−45
if 45 Age(x) 75
⎪
⎩
1 if 75 Age(x) 120.
with graphic.
In cases like this one, it is important to notice that the function μold depends on
the values 45 and 75, as well as on the form of the curve in the sub-interval [45, 75].
It can be supposed, for example, that the above function is supplied by a person
in the range of the fifties, but that one in the seventies would design μold as the
curve
Remark 1.1.10 There are different models for the uses of the same predicate P in X ,
and such uses are reflected in the corresponding models μ P in [0, 1] X . It is because of
this that it is actually important the process of designing the membership functions.
Example 1.1.11 Analogously to the case of big, the predicate A5 = around five in
[0, 10], can have non-linear but quadratic models, as the one given by
⎧
⎨ 0, if x ∈ [0, 4] ∪ [6, 10]
μ A5 (x) = (x − 4)2 , if x ∈ [4, 5]
⎩
(6 − x)2 if x ∈ [5, 6],
whose graphics is
Remark 1.1.12 Since each P in a set X can have different degrees μ P , at each
particular case the meaning of P should be well captured to not represent it by a
mistaken function that will translate a different use of the predicate.
Remark 1.1.13 Given μ P , and the L−set P , the degree is also called the membership
∼
function of the L-set. At this respect,
• x ∈α P , is classically written x ∈
/ P
∼ ∼
• x ∈ω P , is classically written x ∈ P .
∼ ∼
12 1 On the Roots of Fuzzy Sets
1.2.1 Antonyms
a P = −1
P
and verifies,
• μa(a P) = μ(a P) ◦ A = (μ P ◦ A) ◦ A = μ P ◦ (A ◦ A) = μ P ◦ id X = μ P .
Then, for all symmetry A in X , we have an opposite for P. For example, in X =
[0, 10] with μbig (x) = 10
x
, and A(x) = 10 − x, it is μbig (10 − x) = 10−x
10 = 1 − 10 ,
x
If big is represented by
⎧
⎨ 0, if x ∈ [0, 4]
(x−4)
μbig (x) = 4 , if x ∈ [4, 8]
⎩
1, if x ∈ [8, 10],
graphically,
in [0, 1].
The last two examples show a serious trouble. There are no points x ∈ X such
that it is simultaneously μ P (x) = 0 and μa P (x) = 0. The pairs of opposites (P, a P)
for which there is a region in X such that both μ P and μa P take the value 0, called
neutral region, are called regular opposites. Hence, the two above pairs (big, small)
are not regular. Nevertheless, if big is represented by
0, if 0x 7
μbig (x) = (x−7)
3 , if 7 x 10,
14 1 On the Roots of Fuzzy Sets
0, if 3 x 10
μsmall (x) = μbig (10 − x) = (3−x)
3 , if 0 x 3,
with graphics
that shows the neutral region (3, 7). This pair is regular.
1.2.2 Negations
Let it P be a predicate in X , and P = notP its negate. The only we can say about
the relation between P and P is that it is P ⊂ −1
P , since
x P y ⇒ y P x ⇒ μ P (y) μ P (x)
⇒ N (μ P (x)) N (μ P (y)) ⇔ μ P (x) μ P (y).
then
1−a
Nλ (a) = , with λ > −1, for all a ∈ [0, 1].
1 + λa
it results:
• If λ ∈ (−1, 0], then Nλ N0
• If λ ∈ (0, +∞], then N0 < Nλ ,
graphically
μ
Nλ , λ ∈ (−1, 0)
N0
N λ , λ ∈ (0, + ∞)
with graphics.
Example 1.2.2 With μsmall as defined in the second case of 2.2, and with N1 it is
⎧
⎪
⎨ 1, if 6 x 10
1 − μsmall (x)
μnot small (x) = N1 (μsmall (x)) = = 10−x
x−2
, if 2x 6
1 + μsmall (x) ⎪⎩
0 if 0 x 2,
Example 1.2.3 With μbig as defined in the third case of 2.2, and with N0 , it is
1, if 0x 7
μnot big (x) =
3 , 7 x 10,
10−x
if
1.2 Opposite, Negate, and Middle 17
with graphics
With pairs of antonyms (P, a P), it should always be taken into account that con-
ditional statements like “If the bottle is empty, then it is not full”, conduct to the
inequality μa P μnot P , with P = f ull, showing that not P could be taken as the
biggest antonym of P. It is not often the case in which a P = not P, practically it
only happens when aP is such that there is not any linguistic term a P in the language.
When a P and not P are not coincidental, it is said that a P is a strict antonym of P.
When modeling μa P = μ P ◦ A, and μnot P = N ◦ μ P , with a symmetry A of X
and a strong negation N in [0, 1], it results the condition of coherence:
μP ◦ A N ◦ μP ,
results ⎧
⎨ 0, if x ∈ [6, 10]
μsmall (x) = μbig (10 − x) = 4 ,
6−x
if ∈ [2, 6]
⎩
1, if x ∈ [0, 2],
⎧
⎨ 1, if x ∈ [0, 4]
μnot big (x) = 1 − μbig (x) = 4 ,
8−x
if ∈ [4, 8]
⎩
0, if x ∈ [8, 10],
18 1 On the Roots of Fuzzy Sets
whose graphics show that the pair (big, small) is coherent, since μsmall μnot big .
A(x) x
μsmall (x) = μbig (A(x)) = μnot big (x) = 1 −
10 10
hence, A(x) 10 − x is the condition A must satisfy. For example,
• If A1 (x) = 10 − x, it results μsmall (x) = 1 − 10x
= μnot big (x), a non-regular
case.
• If A2 (x) = 10 · 10−x
10+x , for which A2 (x) 10 − x, it results μsmall (x) = μbig (10 ·
10−x
10+x ) = 10−x
10+x , with graphics
10 − x x 10 − x
μsmall (x) = μbig (10 − x) = Nλ (μbig (x)) = Nλ ( ) =
10 10 10 + λx
1.2 Opposite, Negate, and Middle 19
Given a regular pair of antonyms (P, a P), the middle or medium term of them, is
the predicate MP = not P and Not aP, as it was said in Sect. 1.2.1. Although we
don’t yet studied the diverse models for the conjunction ‘and’, by the moment let us
take the model given by the operation minimum (= min) in [0, 1]. Then,
μ small
μ big
6 8
The triplets (P, M P, a P) play an important role in the applications of fuzzy sets.
20 1 On the Roots of Fuzzy Sets
1.3 AND/OR
Let P, Q be two predicates in X . Consider the new predicates ‘P and Q’, and ‘P or
Q’, used by means of:
• ‘x is P and Q’ ⇔ ‘x is P’ and ‘x is Q’
• ‘x is P or Q’ ⇔ Not (Not ‘x is P’ and Not ‘x is Q’) ⇔ Not (‘x is P ’ and ’ x is
Q ’) ⇔ x is ‘(P and Q )’,
with respective primary uses P and Q ⊂ P ∩ Q , and P or Q . Take L-degrees
μP , μQ .
1.3.1 AND
In the case (L , ) is a lattice (see Sect. 1.6) with the minimum operation · = min,
and ∗ ·, it implies
μ P and Q (x) μ P (x) · μ Q (x).
1.3.2 OR
μ P or Q = μP and Q .
Remark 1.3.1 The lattice operation ·(+) is not necessarily the only operation ∗(⊕),
that can exist. In the case (L , ) is not a lattice for ·(+), there can also exist other
operations ∗ and ⊕. For example, L = [0, 1] is not a lattice with ∗ = pr od, but
(a ∗ b) ⊕ (a ∗ b) a ⊕ b,
Remark 1.3.4 It should be noticed that what has been presented is sufficient but not
necessary for the representation of μ P and Q and μ P or Q b means of μ P and μ Q .
μ M P (x) = μ P and (a P) (x) = μ P (x) ∗ μ(a P) (x) = N (μ P (x)) ∗ N (μ P (A(x))),
for all x in X .
‘x is (P is τ )’ := x is P is τ ,
μ P is τ = μτ ◦ μ P
0, if 0 x 0.5
μτ (x) =
1, if 0.5 x 1,
it results
0, if 5 x < 10
μτ ◦ μ P (x) =
1, if 0 x < 5,
Linguistic modifiers or linguistic hedges, m, are adverbs acting on P just in the con-
catenated form m P. For example, with m = very and P = tall, it is m P = very tall.
A characteristic that linguistically distinguishes imprecise predicates from pre-
cise ones, is that in the first case and once P and m are given, mP is immediately
understood. If P is precise (for example, P = even in the set of natural numbers),
mP needs of a new definition to be understood (what it means very even?). Modifiers
only modify, but do not change abruptly imprecise predicates.
If P in X is with the use ( P , μ P ), and m in μ P (x) ⊂ L is with ⊂ m and
μm , provided m P ⊂ P , it can be taken the degree
μm P = μm ◦ μ P ,
(x is P, y is Q) ∈ R(P, Q),
for all (x, y) ∈ X × Y , and it again remains to be tested that μ Q|P is actually a
L-degree for Q|P. For example,
• If Q|P = Q × P , and J is non-decreasing in both variables, it is
(x1 , y1 ) Q|P (x2 , y2 ) ⇔ x1 P x2 , y1 Q y2 ⇒ μ P (x1 ) μ P (x2 ), and
μ Q (y1 ) μ Q (y2 ) ⇒
or
μ Q|P (x1 , y1 ) μ Q|P (x2 , y2 ).
• If Q|P = −1 −1
P × Q , and J is decreasing in both variables, it also follows the
same conclusion,
• If Q|P = −1
P × Q , and J is decreasing in its first variable, and non-decreasing
in the second, it also follows the same conclusion.
Etc.
1.4 Qualified, Modified, and Constrained Predicates 25
ω, if (x, y) ∈ R
μ R (x, y) =
α, otherwise.
The meaning of words is not fixed for all people and all context. For example,
in a dinner with three commensals the deliciousness of the dessert plates could
easily result in three different orderings of such plates. Since language is a social
phenomenon, also meaning is such, and it is possible to talk on the meaning of
predicates for a group of people in, of course, a given context.
For a group of people G = { p1 , . . . , pm }, a predicate P on X can show m primary
meanings P,i , 1 i m. Since
m
( P,i ) = P,G
i=1
26 1 On the Roots of Fuzzy Sets
Notice that provided all P,i are preorders, P,G is also a preorder.
If m L − degr ees μ(i)P are known for each primary meaning P,i , since
• x = P,G y ⇔ x = P,1 y & . . . & x = P,m y,
• x P,G y ⇔ x P,1 y & . . . & x P,m y,
for each function : L m → L, non-decreasing in each place i between 1 and m
(for example, if a b then (a, x2 , . . . , xm ) (b, x2 , . . . , xm )), or aggregation
function, it results
• x P,G y ⇒ (μ(1) (m) (1) (m)
P (x), . . . , μ P (x)) (μ P (y), . . . , μ P (y)),
that allows to take
(1) (m)
μGP (X ) = (μ P (x), . . . , μ P (x)), for all x ∈ X,
1.4.5 Synonims
If P and Q are id X -primary synonyms, it is said that they are exact or perfect
synonyms when μ P = μ Q , and it results ( P , μ P ) = ( Q , μ Q ).
1.4 Qualified, Modified, and Constrained Predicates 27
it follows that
μ Q = μ P ◦ u −1
y1 μ Q y2 ⇔ u −1 (y1 ) μ P u −1 (y2 ),
or
x1 μ P x2 ⇔ u(x1 ) μ Q u(x2 ),
For example, with the before mentioned predicates short and small, it is
y1 μ Q y2 ⇔ y1 /10 μ P y2 /10.
Remark 1.4.5 The definition of primary meaning is just a formal one trying to
approach an important aspect of the meaning of linguistic predicates, when act-
ing on a given universe of discourse. The same can be said about the definition of
u-primary synonyms with which it does not hold, in general, that a pair of linguistic
synonyms are necessarily u−primary synonyms. Anyway, what can be said is that
Q is a migration of P to the universe Y .
Remark 1.4.6 In some way, the current meaning of a predicate, the form in which it
is used today in the plane language, partially inherits its past meanings.
28 1 On the Roots of Fuzzy Sets
Linguistic variables are basic tools in most application of fuzzy sets in the technol-
ogy’s field. They do mainly appear when linguistically describing the behavior of the
physical variables of a system. A linguistic variable explicits a concept by (linguis-
tically) granulating some elemental components of it, by showing the perceptually
distinguishable shades that are relevant for the corresponding application.
A linguistic variable LV is formed after considering
1. Its principal predicate, P
2. One of the opposites of P, a P
3. Some linguistic modifiers m 1 , . . . , m n ,
and by adding:
4. Its negate (not P), or the middle-predicate (MP), or P and Q, or not m 1 P, or P
and m 2 a P, . . .
Then LV, is called the linguistic variable generated by P, and reflects the linguistic
granulation perceived for the concept. For example,
• LV = Age, is Age = { young, old, middle-aged, not old, not very young, …}
• LV = Truth, is Truth = { true, false, very false, not very true, …}
• LV = Temperature, is Temperature = { cold, hot, warm, not cold, not very hot, …}
• LV = Size, is Size = {large, small, medium, very large, …}
• LV = Height, is
1.5 Linguistic Variables 29
– For buildings, Height = {low, high, medium, very high, not very low, …}
– For people, Height = {tall, short, medium, very tall, more or less short, …}
• LV = Speed, is Speed = {fast, slow, very slow, more or less fast, not fast, …}
Although the number of terms in a Linguistic Variable can be large, usually it is
comprised between 5 and 9 (7 ± 2) since, in a lot of cases, less than 5 shades is poor
and more than 9 is excessive. Anyway in a good number of applications there only
appear the three terms P, aP, and MP.
Usually, in the applications, the variables range in the set of real numbers and,
because of this, the predicate P acts in some interval of the real line. For example,
the linguistic variable ‘Temperature’ in the interval between
−10, and 50 degrees Celsius, is often represented by only the three terms μcold , μhot ,
μwar m , with warm = not cold and not hot.
Analogously, ‘Height’ for people, can be represented in [0, 2] meters by the lin-
guistic variable with the four terms μtall , μshor t , μver y shor t , μmor e or less tall in the
following figure,
It is sometimes useful in the applications that, once ordered in some sequence, the
fuzzy sets in a linguistic variable L V = {μ0 , μ1 , . . . , μn } do form what is called a
fuzzy partition (or a unit’s partition), that is, verifying
n
μ j (x) = 1, ∀x ∈ X.
j=0
30 1 On the Roots of Fuzzy Sets
⎧
⎪ 0, if 0 x j − 1,
⎪
⎪
⎨x +1− j if j − 1 x j,
μ j (x) = for 1 j 3
⎪
⎪ j +1−x if j x j + 1,
⎪
⎩
1, if j + 1 x 4,
0, if 0 x 3,
μ4 (x) =
x −3 if 3 x 4,
Graphically,
Obviously,
4
• If 0 x 1, j=0 μ j (x) = μ0 (x) + μ1 (x) = 1 − x + x = 1
• If 1 j 3, j x j + 1, 4j=0 μ j (x) = μ j (x) + μ j+1 (x) = ( j + 1 − x) +
(x + 1 − j − 1) = 1.
• If 3 x 4, 4j=0 μ j (x) = μ3 (x) + μ4 (x) = 4 − x + x − 3 = 1
Hence {μ0 , μ1 , . . . , μn } is a fuzzy partition of [0, 4]. Notice that each μ j can be
labeled by the predicate around j= A j , that is μ j = μ A j .
1.5 Linguistic Variables 31
0, if 0 x 6,
μ2 (x) =
4 , 6 x 10.
x−6
if
Graphically
Since, obviously, μ0 (x) + μ1 (x) + μ2 (x) = 1 for all x ∈ [0, 10], {μ0 , μ1 , μ2 } is
a fuzzy partition of the interval [0, 10].
Example 1.5.3 By proceeding in the same way that in last example, it is easy to
prove that it is
{μcold , μhot , μwar m },
Notice that, with the strong negation N0 = 1−id, if {μ0 , μ1 , μ2 } is a fuzzy partition
of X , it is
μ1 (x) + μ2 (x) = 1 − μ0 (x), or μ1 + μ2 = μ0
μ0 (x) + μ1 (x) = 1 − μ2 (x), or μ0 + μ1 = μ2
μ0 (x) + μ2 (x) = 1 − μ1 (x), or μ0 + μ2 = μ1
The complement (with N0 ) of each element in a fuzzy partition, is just the addition
of the other elements in the partition.
Provided (L , ) only verifies the reflexive and the transitive properties, is called
a preordered set.
2. An algebraic structure (L , ·, +) is a lattice, provided the binary operations · and
+ verify,
• Both operations are commutative and associative, that is, a · b = b · a, a + b =
b + a, a · (b · c) = (a · b) · c, a + (b + c) = (a + b) + c, for all a, b, c in L.
• Both operations are idempotent, that is, a · a = a, a + a = a, for all a in L.
Without proof, let’s state what follows:
1.6.1 Examples
2.1 Introduction
From now on it will be only considered the case in which (L , ) = ([0, 1], ), that
is, of Zadeh’s fuzzy sets, with predicates P in X known through a degree μ P : X →
[0, 1], and without knowing, necessarily, its primary use P . The set of all fuzzy sets
in X , [0, 1] X , will be also denoted by F(X ). In this case, the preorder μ P is linear,
or total, since for all x, y in X it is either μ P (x) μ P (y), or μ P (y) μ P (x), that
is, it is either x μ P y or y μ P x for all x, y in X . Hence, μ P rarely will perfectly
reflect the primary use of P in X , since P is usually not linear.
In the case in which X is finite, X = {x1 , . . . , xn }, the fuzzy sets μ ∈ [0, 1] X ,
will be represented by
with the convention that if some term μ(x j )/x j does not appear, is that it is μ(x j ) = 0.
For example, with X = {1, 2, 3, 4}, the expression
refers to the fuzzy set in X given by μ(x1 ) = 0.5, μ(x2 ) = 0.7, μ(x3 ) = 0, μ(x4 ) = 1.
Analogously, the fuzzy set μ = N0 ◦ μ (N0 = 1 − id) is
If A, B are crisp subsets in X and Y , respectively, that is, A ∈ P(X ) and B ∈ P(X ),
its cartesian product A × B = {(a, b); a ∈ A, b ∈ B} ⊂ X × Y , is with the
© Springer International Publishing Switzerland 2015 35
E. Trillas and L. Eciolaza, Fuzzy Logic, Studies in Fuzziness
and Soft Computing 320, DOI 10.1007/978-3-319-14203-6_2
36 2 Algebras of Fuzzy Sets
it is
μ × σ = 0.9/(x1 , y1 ) + 0.7/(x1 , y2 ) + 0.8/(x2 , y1 ) + 0.7/(x2 , y2 ),
x y
(μbig × μsmall )(x, y) = min( , 1 − ),
5 7
the representation of the cartesian product as a surface contained in the cube [0, 5] ×
[0, 7] × [0, 1].
Of course, if μ = μ A ∈ {0, 1} X , σ = μ B ∈ {0, 1} X , it is not only μ × σ ∈ {0, 1} X
but μ × σ = μ A×B .
f : P(X ) → P(Y ), A → f (A).
2.1 Introduction 37
and f is known as the ‘extension’ of f to the fuzzy parts, and the definition as the
Zadeh’s Extension Principle.
For example, if f : [0, 10] → [0, 1], is given by f (x) = 1 − 10 x
, the
fuzzy set μ(x) = 10 in [0, 1][0,10] extends to the fuzzy set in [0, 1],
x
f (μ)(y) =
sup{μ(x); f (x) = y} = sup{ 10x
; 1 − 10x
= y} = 1 − y, for all y ∈ [0, 1].
If X = {1, 2, 3, 4}, Y = {a, b, c}, the mapping f : X → Y such that
extends the fuzzy set μ = 1/1 + 0.4/2 + 1/3 + 0.7/4 in F(X ), to the fuzzy set
f (μ)
in F(Y ) with values,
f (μ)(a) = max{μ(x); x ∈ f −1 (a)} = max{μ(1), μ(2)} = max(1, 0.4) = 1
f (μ)(b) = max{μ(3), μ(4)} = 1
f (μ)(c) = 0, since f −1 (c) = ∅.
Hence,
f (μ) = 1/a + 1/b,
Like with the cartesian product and with the extension principle, all operations with
fuzzy sets must reproduce, when the data are crisp, the corresponding result obtained
in the crisp theory. This is the principle of preservation of the classical case, that is
forced by the will, and the necessity, of including all ‘the classical’ as a particular
case of the algebras of fuzzy sets.
To illustrate this preservation’s principle, let us show a negative example. With
X = [0, 1], and all μ ∈ [0, 1][0,1] , the function
verifies:
• μ∗0 (x) = 1 − μ0 (1 − x) = 1: μ∗0 = μ1
• μ∗1 (x) = 1 − μ1 (1 − x) = 1 − 1 = 0: μ∗1 = μ0
• μ σ ⇒ 1 − σ(1 − x) 1 − μ(1 − x) ⇒ σ ∗ μ∗
• μ∗∗ (x) = 1 − μ∗ (1 − x) = 1 − [1 − μ(x)] = μ(x): μ∗∗ = μ.
Hence, it could seem that the function μ → μ∗ can be taken as a “strong negation”
for the fuzzy sets in [0, 1], but it is not the case. Notice that if μ ∈ {0, 1}[0,1] , then
it should be also μ∗ ∈ {0, 1}[0,1] , that is, if μ represents a classical subset of [0, 1],
also μ∗ should represent not only a classical subset but precisely the complement of
μ. But with A = [0, 21 ] ⊂ X ,
1, 0 x 0.5,
μ A (x) =
0, 0.5 < x 1,
follows,
1, 0.5 < x 1, 0, 0.5 < x 1,
μ∗A (x) = 1 − μ A (1 − x) = 1 − =
0, 0 x 0.5, 1, 0 x 0.5
that represents the subset [0, 0.5], but not Ac = (0.5, 1]. The unary operation ∗
violates the preservation principle, and hence it cannot be taken into account to
negate fuzzy sets.
2.1.4 Resolution
Let us denote by μr the constant fuzzy sets in [0, 1] X , μr (x) = r , for r ∈ [0, 1]
and all x ∈ X . Notice that in {0, 1} X there are only the “constants” μ0 and μ1 , that
correspond to the sets ∅ and X , respectively.
Given μ ∈ [0, 1] X , let us denote by μ(r ) the fuzzy (crisp) set
(r ) 1, if r μ(x),
μ (x) =
0, otherwise,
for all r ∈ [0, 1], and by [μ(r ) ] the corresponding classical subset {x ∈ X ; r μ(x)}.
These sets are called the r-cuts of μ and it is always [μ(0) ] = X .
For example, in the following figures are shown, respectively, the constant fuzzy
set μ0.5 , and the 0.3-cut of two different fuzzy sets.
2.1 Introduction 39
40 2 Algebras of Fuzzy Sets
Example 2.1.3 With X = {1, 2, 3}, take μ : X × X → [0, 1], given by μ(i, j) =
min(i, j)
3 . This fuzzy set in X × X can be represented either by the matrix
⎛1 1 1 ⎞
3 3 3
⎜ ⎟
⎝ 1
3
2
3
2
3 ⎠,
1 2
3 3 1
or by the graph
1/3 2/3
1/3
1 2
1/3
2/3
1/3 2/3
1/3
2 1
Since the matrices of μ(1) , μ( 3 ) , and μ( 3 ) are respectively
⎛ ⎞ ⎛ ⎞ ⎛ ⎞
000 000 111
⎝0 0 0⎠,⎝0 1 1⎠,⎝1 1 1⎠,
001 011 111
it results⎛ ⎛ ⎛ ⎞⎞ ⎛ ⎛ ⎞⎞ ⎛ ⎛ ⎞⎞⎞
000 000 111
⎜ ⎜ ⎜ ⎟⎟ ⎜ ⎜ ⎟⎟ ⎜ ⎜ ⎟⎟⎟
max ⎝min ⎝1, ⎝ 0 0 0 ⎠⎠ , min ⎝ 23 , ⎝ 0 1 1 ⎠⎠ , min ⎝ 13 , ⎝ 1 1 1 ⎠⎠⎠
001 011 111
⎛⎛ ⎞ ⎛ ⎞ ⎛ 1 1 1 ⎞⎞ ⎛ 1 1 1 ⎞
000 0 0 0 3 3 3 3 3 3
⎜⎜ ⎟ ⎜ ⎟ ⎜ ⎟⎟ ⎜ ⎟
= max ⎝⎝ 0 0 0 ⎠ , ⎝ 0 23 23 ⎠ , ⎝ 13 13 13 ⎠⎠ = ⎝ 31 23 23 ⎠ ,
001 0 23 23 1 1 1
3 3 3
1 2
3 3 1
accordingly with the theorem of resolution.
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 41
2.2.1 Introduction
Functions
μ ∈ F(X ) = [0, 1] X ,
• Notice that F(X ) will be taken as ‘ordered’ (partially) by means of the binary
pointwise relation
μ A μ B ⇔ A ⊂ B; μ A = μ B ⇔ A = B,
and
x ∈ A ⇔ μ A (x) = 1; x∈
/ A ⇔ μ A (x) = 0.
• The height of μ ∈ F(X ) is H (μ) = Sup μ(x) = sup μ. In the last example, it is
x∈X
H (μ) = 1. In the finite example
n
|μ| = μ(xi )
i=1
Remark 2.2.1 The pointwise definition of fuzzy sets inclusion implies that, for exam-
ple, the fuzzy sets
in X , do not verify μ σ although it is σ(x2 ) < μ(x2 ), σ(x3 ) = μ(x3 ), σ(x4 ) <
μ(x4 ), but σ(x1 ) > μ(x1 ), with σ(x1 ) − μ(x1 ) = 0.00001. Pointwise ‘inclusion’ is
strongly affected by very small variations of the membership values. Actually, it is
not a flexible, or fuzzy, concept, but a crisp one.
Because of this, it could be preferable to take the inclusion of fuzzy sets as an
gradable concept r (r ∈ [0, 1]), and a used definition of which is
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 43
| min(μ, σ)|
μ r σ ⇔ r,
|μ|
n
with | min(μ, σ)| = i=1 min(μ(xi ), σ(xi )).
In last example, it is | min(μ, σ)| = 0.7 + 0.7 + 1 + 0.6 = 3, |μ| = 0.7 + 0.8 +
1 + 0.7 = 3.2, and r = 3/3.2 = 0.9375 ≈ 0.94. That is, μ σ : μ ‘is almost
0.9375
included in’ σ.
Since |σ| = 0.70001 + 0.7 + 1 + 0.6 = 3.00001, it is | min(μ,σ)|
|σ| = 0.9999, or
σ 0.9999 μ. That is, σ is more included in μ, than μ is included in σ!
μ σ ⇒ μ σ.
1
7
7
7
min(μ A (i), μ B (i)) = 3, μ B (i) = 4, μ A (i) = 4.
i=1 i=1 i=1
hence
μ A 3 μ B , or A ⊂ 3 B
4 4
μ B 3 μ A , or B ⊂ 3 A
4 4
Once F(X ) = ([0, 1] X ; ; =) is taken, a general algebra of fuzzy sets comes from
endowing F(X ) with three operations:
1. : [0, 1] X → [0, 1] X ,
2. · : [0, 1] X × [0, 1] X → [0, 1] X ,
3. + : [0, 1] X × [0, 1] X → [0, 1] X ,
respectively called the complement μ of μ, the intersection μ · σ of ‘μ and σ’, and
the union μ + σ of ‘μ or σ’. Then ([0, 1] X ; ; =; ·; +; ), is called an algebra of
fuzzy sets, provided the following laws do hold:
44 2 Algebras of Fuzzy Sets
Remark 2.2.3 It is not difficult to prove that no general algebra of fuzzy sets is a
Boolean algebra. The proof comes from the fact that to be a Boolean algebra would
imply μ · μ = μ0 and μ + μ = μ1 for all μ ∈ [0, 1] X , and consists in finding some
μ for which these equalities are not satisfied.
(μ + σ) = μ · σ , (μ · σ) = μ + σ ,
or μ · σ min(μ, σ). Hence, the operation min is the greatest possible intersection
of fuzzy sets. Analogously, μ0 μ, μ0 σ ⇒ μ0 + σ = σ μ + σ, μ + μ0 = μ
μ + σ, and max(μ, σ) μ + σ: The operation max is the smallest possible union of
fuzzy sets.
Proof The operations min, and max are obviously idempotent. Let us show that if ·
is idempotent, it must be · = min. By Theorem 2.2.6, it is always μ · σ min(μ, σ),
and the idempotency of · implies
min(μ, σ) μ · σ,
(μ · σ) = μ + σ , (μ + σ) = μ · σ
result equivalent. Then, it is enough to prove max(μ, σ) = (min(μ , σ )) , for all
μ, σ in [0, 1] X . Since
min(μ , σ ) μ , min(μ , σ ) σ ,
it follows
μ = μ (min(μ , σ )) , σ = σ (min(μ , σ )) ,
and
max(μ, σ) (min(μ , σ )) (2.1)
Theorem 2.2.10 (Kleene’s Law) In all general algebra of fuzzy sets it holds the law
μ · μ σ + σ ,
Remark 2.2.11 Concerning duality, Theorem 2.2.9 only states that the algebra given
by the triplets (min, max, ), with involutive, are dual algebras. But they are not the
only dual algebras. For example, with · = product,
μ + σ = (μ · σ ) ,
Remark 2.2.12 It is easy to prove that, for each algebra of fuzzy sets ([0, 1] X , ·, +, ),
the operation
μ + σ = (μ · σ ) ,
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 47
gives the new algebra ([0, 1] X , ·, + , ). If the complement is involutive (μ = μ),
then (μ + σ) = μ · σ .
Analogously, with the operation μ · σ = (μ + σ ) , one has the new algebra
([0, 1] X , · , + , ) and, if is involutive, (μ · σ) = μ + σ .
A ⊆ Ac ⇒ A ∩ A ⊆ A ∩ Ac ⇒ A ⊆ ∅ ⇒ A = ∅.
Perhaps, this is the reason of the difficulty children do have on accepting that ∅ is a
set!
Within an algebra of fuzzy sets there are many self-contradictory fuzzy sets. For
example, with N = 1 − id it is
1
μ μ ⇔ μ(x) 1 − μ(x) ⇔ μ(x) , ∀x ∈ X,
2
hence: μ is self-contradictory if and only if μ μ 1 . Analogously, with the strong
2
negation N (x) = 1−x
1+x , it is
1 − μ(x)
μ μ ⇔ μ(x) ⇔ μ(x)2 + 2μ(x) − 1 0
1 + μ(x)
√
⇔ μ(x) 2 − 1, ∀x ∈ X,
Proof It is,
• μ max(μ, μ ) μ + μ ⇒ (μ + μ ) μ ⇒ (μ ) ((μ + μ ) )
• μ max(μ, μ ) μ + μ ⇒ (μ + μ ) (μ )
then (μ + μ ) ((μ + μ ) ) .
μ · μ = μ0 , μ + μ = μ1 ,
that is, to find for which intersections · and which unions +, these equations do hold.
Of course, they do not hold in all cases, for example, with N = 1 − id,
• If · = min, it is not always min(μ(x), 1 − μ(x)) = 0,
• If + = max, it is not always max(μ(x), 1 − μ(x)) = 1,
• If · = W, (W (a, b) = max(0, a + b − 1)) it is W (a, 1 − a) = max(0, a + 1 −
a − 1)) = 0, and W (μ(x), 1 − μ(x))) = 0 for all x in X
• If + = W ∗ , (W ∗ (a, b) = min(1, a + b)), it is W ∗ (μ(x), 1 − μ(x))) = 1 for all x
in X .
That is, there are algebras of fuzzy sets where this forms of non-contradiction
or excluded-middle hold, and algebras where this principles do not jointly hold.
In the algebras with the triplet (min, max, 1 − id) do not hold both principles, in
the algebras with (W, max, 1 − id) it holds the principle of non-contradiction but
not that of excluded-middle, in the algebras with (min, W ∗ , 1 − id) it holds the
excluded-middle but not the principle of non-contradiction, and in the algebras with
(W, W ∗ , 1 − id) both principles hold.
Remark 2.2.16 Results in Theorems 2.2.13 and 2.2.14 challenge the usual statement
that in fuzzy sets the basic principles of Non-contradiction and Excluded-middle fail.
A statement that could conduct to believe that fuzzy set algebras are not properly
grounded in a solid ground.
The fact is, notwithstanding, that these two principles were established before
the current ways of considering the problems of logic and, of course, before the
nomenclature of set theory. In set theory (or Boolean algebras), A ∩ Ac = ∅ and
A ∩ Ac ⊂ (A ∩ Ac )c are equivalent formulas since, as it was said, it is
B = ∅ ⇔ B ⊂ Bc,
an equivalence only verified in the setting of ortholattices (of which Boolean algebras
are a particular case), but that does not hold on weaker algebraic structures like it
is, for example, the case of the above defined algebras of fuzzy sets. Let us call
‘restricted’ the principles stated by μ · μ = μ0 and μ + μ = μ1 .
50 2 Algebras of Fuzzy Sets
(μ ∗ σ)(x) =
∗(μ(x), σ(x)),
for all μ, σ in [0, 1] X and all x in X . Of course, by this formula, a numerical operation
∗ allows to define an operation ∗ for fuzzy sets.
( f (μ))(x) =
f (μ(x)),
μ (x) = 1 − μ(x),
There are three functions F : [0, 1]×[0, 1] → [0, 1], G : [0, 1]×[0, 1] → [0, 1],
and N : [0, 1] → [0, 1] with which
It is easy to show that this operation verifies the laws b, d and e in Sect. 2.1.1. Hence,
it is a union for fuzzy sets that, in addition, is commutative. It is not idempotent,
since if μ ∈ [0, 1] X − {0, 1} X , it is (μ + μ)(x) = H (μ) = μ(x). It does not exist a
function G : [0, 1] × [0, 1] → [0, 1] such that
for all x ∈ X and all μ, σ[0, 1] X . Indeed, let us suppose that such a G does exist,
and take μ = μ0.5 . Then
• With σ = μ0 , it is (μ + σ)(x) = max( 21 , 0) = 21 . Hence G( 21 , 0) = 21 .
• With σ(x) = x, is (μ + σ)(x) = max(H (μ), H (σ)) = max( 21 , 1) = 1, and
(μ + σ)(0) = 1 = G( 21 , 0), that is absurd.
To have a not-decomposable ‘intersection’, it is enough to define, with μ = 1 − μ,
the dual operation,
An standard algebra of fuzzy sets is a decomposable algebra of fuzzy sets such that:
1. μ · σ = σ · μ, for all μ, σ in [0, 1] X (· is commutative)
2. μ + σ = σ + μ, for all μ, σ in [0, 1] X (+ is commutative)
3. μ · (σ · λ) = (μ · σ) · λ, for all μ, σ, λ in [0, 1] X (· is associative)
4. μ + (σ + λ) = (μ + σ) + λ, for all μ, σ, λ in [0, 1] X (+ is associative)
5. μ = μ, for all μ in [0, 1] X ( is involutive).
Hence, writing
μ · σ = T ◦ (μ × σ), μ + σ = S ◦ (μ × σ), μ = N ◦ μ,
functions T, S : [0, 1] × [0, 1] → [0, 1] and N : [0, 1] → [0, 1], in addition to the
corresponding general properties stated before, must verify
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 53
• T is commutative, S is commutative
• T is associative, S is associative
• N is involutive,
that is:
• T (a, b) = T (b, a), S(a, b) = S(b, a), for all a, b in [0, 1]
• T (a, T (b, c)) = T (T (a, b), c), S(a, S(b, c)) = S(S(a, b), c), for all a, b, c in
[0, 1]
• N (N (a)) = a, for all a in [0, 1], or N ◦ N = id, or N = N −1 .
Functions T and S are called t-norms and t-conorms, respectively. Functions N
are strong negations. Hence, ([0, 1], T, ) is an ordered semigroup with neutral 1,
and absorbent 0, and ([0, 1], S, ) is also an ordered semigroup but with neutral 0
and absorbent 1. Since N (1) = 0, it seems that this two kind of ordered semigroups
should show some character of duality. This duality goes in the way:
• If T is a t-norm, TN = N ◦ S ◦ (N × N ) is a t-conorm
• If S is a t-conorm, S N = N ◦ S ◦ (N × N ) is a t-norm
that are easy to prove. Of course, from Sect. 2.1.4,
• If T is a t-norm, T min, and min is a t-norm
• If S is a t-conorm, max S, and max is a t-conorm
Hence, for all t-norm T and all t-conorm S:
T min max S,
is a t-conorm such that S Z ∗ for all t-conorm S. Hence, for all t-norm T and all
t-conorm S,
Z T min max S Z ∗ .
1 Notice that T S mean T (a, b) S(a, b), for all (a, b) ∈ [0, 1] × [0, 1].
54 2 Algebras of Fuzzy Sets
W T pr od min,
∗
it follows max T pr ∗
od W , and
∗ ∗
Z W T pr od min max T pr od W .
it follows that t-norm W has zero-divisors, for example, W (0.5, 0.4) = 0; t-norms
min and T pr od do not have zero-divisors:
• T pr od (a, b) = 0 ⇔ a = 0, or b = 0
• min(a, b) = 0 ⇔ a = 0, or b = 0
Proposition 2.2.22 The only idempotent t-norm, i.e. T (a, a) = a, for all a ∈ [0, 1],
is T = min.
Proposition 2.2.23 The only idempotent t-conorm, i.e. S(a, a) = a, for all a ∈
[0, 1], is S = max.
Proof If S is idempotent, max(a, b) = S(max(a, b), max(a, b)) ≥ S(a, b), since
max(a, b) ≥ a, and max(a, b) ≥ b. Hence, max(a, b) ≥ S(a, b) ≥ max(a, b)
implies S = max.
Remark 2.2.25 Analogous considerations can be made for t-conorms. There are
∗
discontinuous t-conorms like Z ∗ , and continuous ones like T pr ∗
od and W . The only
for which all elements in [0, 1] are idempotent is S = max. Since S(0, 0) = 0
and S(1, 1) = 1, 0 and 1 are always idempotent, and there are t-conorms that only
have these two idempotents (like T pr∗ ∗
od and W ), as well as those that have some
idempotents different from 0,1. There are t-conorms without one-divisors, like max
∗ , and t-conorms with one-divisors like W ∗ , for example, W ∗ (0.5, 0.5) =
and T pr od
min(1, 1) = 1.
Remark 2.2.26 There is not a characterization theorem for all t-norms (t-conorms),
but it is a characterization of the continuous t-norms (t-conorms) that will be pre-
sented by means of the following, and easy to prove, results:
• If ϕ : [0, 1] → [0, 1] verifies, (1) If x y, then ϕ(x) ϕ(y), (2) ϕ is bijective,
(3) ϕ(0) = 0, ϕ(1) = 1 (ϕ is an order-automorphism of the ordered interval
([0, 1], )), and T is a t-norm, then Tϕ = ϕ−1 ◦ T ◦ (ϕ × ϕ) is also a t-norm.
Given T , the set {Tϕ ; ϕ an order-automorphism} is called the family of T .
• T is a continuous t-norm if and only if all t-norms Tϕ are continuous.
• If S is a t-conorm, then Sϕ = ϕ−1 ◦ S ◦ (ϕ × ϕ) is also a t-conorm, and S
is continuous if and only if all t-conorms Sϕ are continuous, the set {Sϕ ; ϕ an
order-automorphism} is called the family of S.
In particular,
• The family of T = min, is reduced to the only t-norm min, since ϕ−1 (min(ϕ(a),
ϕ(b)) = min(ϕ−1 (ϕ(a)), ϕ−1 (ϕ(b))) = min(a, b)
• The family of T pr od contains all continuous t-norms of the form prodϕ (a, b) =
ϕ−1 (ϕ(a).(ϕ(b)).
• The family of W contains all t-norms of the form Wϕ (a, b) = ϕ−1 (W (ϕ(a),
ϕ(b)) = ϕ−1 (max(0, ϕ(a) + ϕ(b) − 1)), and all of them are continuous t-norms.
Notice that no t-norm in the family { pr odϕ } has zero-divisors, since pr odϕ (a, b) =
0 ⇔ ϕ(a) · ϕ(b) = 0 ⇔ ϕ(a) = 0 or ϕ(b) = 0 ⇔ a = 0, or b = 0. Instead
all t-norms Wϕ have zero-divisors, since Wϕ (a, b) = 0 ⇔ max(0, ϕ(a) + ϕ(b) −
1) ⇔ ϕ(a) + ϕ(b) 1. Of course, neither t-norms pr odϕ , nor Wϕ , have more
idempotents than 0 and 1:
• a = Wϕ (a, a) = ϕ−1 (max(0, 2ϕ(a) − 1)) ⇔ ϕ(a) = max((0, 2ϕ(a) − 1) ⇔
ϕ(a) = 0 or ϕ(a) = 1 or a = 0 or a = 1.
• a = pr odϕ (a, a) = ϕ−1 (ϕ(a).(ϕ(a)) ⇔ ϕ(a) = ϕ(a).(ϕ(a) ⇔ ϕ(a) = 0 or
ϕ(a) = 1 or a = 0 or a = 1.
Analogously,
56 2 Algebras of Fuzzy Sets
pr odϕ∗ (a, a) = ϕ−1 ( pr od ∗ (ϕ(a), ϕ(b))) = ϕ−1 (ϕ(a) + ϕ(b) − ϕ(a) · ϕ(b))
but
√
Pr odϕ∗ (a, b) = ϕ−1 (ϕ(a) + ϕ(b) − ϕ(a) · ϕ(b)) = a + br − a r · br .
r r
As it was said before, an strong negation is a function N : [0, 1] → [0, 1] such that
• N (0) = 1
• If a b, then N (b) N (a)
• N (N (a)) = a, for all a ∈ [0, 1], or N 2 = id.
Notice that N 2 = id is equivalent to N = N −1 , that shows N is a continuous
function: It is N (1) = N (N (0)) = 0, and if a < b it should be N (b) < N (a) since
N (b) = N (a) would imply N (N (b)) = N (N (a)), or a = b. Hence, N is strictly
decreasing.
Since N is continuous, the equation N (x) = x has solutions, but there is only
one. Suppose N (x1 ) = x1 and N (x2 ) = x2 . Either x1 x2 , or x2 < x1 . In the
first case, it follows N (x2 ) N (x1 ), or x2 x1 , and x1 = x2 . In the second case,
N (x1 ) < N (x2 ), or x1 < x2 , that is absurd. Then, each strong negation has a single
fixed point N (n) = n, in the open interval (0, 1), since N (0) = 1, N (1) = 0 show
that 0 and 1 are not fixed points.
Remark 2.2.28 In the classical case (a Boolean algebra L, or a power set P(X )), the
transformation
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 57
the equation μ = μ , N (μ(x)) = μ(x) for all x in X , has the only solution μ(x) = n
for all x in X , that is μ = μn , with n = N (n) the fix point of N .
In the fuzzy case, that mapping F shows a kind of symmetry that is not in the
crisp case.
An order-automorphism of the ordered unit interval ([0, 1], ) ϕ : [0, 1] →
[0, 1], verifies by definition,
• ϕ(0) = 0, ϕ(1) = 1
• If a < b, then ϕ(a) < ϕ(b).
Hence, ϕ is continuous, and its inverse function ϕ−1 verifies,
• ϕ−1 (0) = 0, ϕ−1 (1) = 1
• If a < b, then ϕ−1 (a) < ϕ−1 (b).
Let us denote by Nϕ the function Nϕ : [0, 1] → [0, 1] defined by
Proof Let it be n = N (n) ∈ (0, 1) the fixed point of N , and consider an strictly
non-decreasing function h : [0, n] → [0, 21 ] such that h(0) = 0 and h(n) = 21 . With
h define the function ϕ : [0, 1] → [0, 1] by
h(x), if x ∈ [0, n]
ϕ(x) =
1 − h(N (x)), if x ∈ (n, 1].
58 2 Algebras of Fuzzy Sets
This function ϕ is, obviously, continuous, strictly increasing,2 and verifies ϕ(0) =
h(0) = 0, ϕ(1) = 1 − h(N (1)) = 1 − h(0) = 1. Then
• If x ∈ [0, n), or N (x) ∈ (n, 1], ϕ(N (x)) = 1 − h(x) = 1 − ϕ(x), and N (x) =
ϕ−1 (1 − ϕ(x)).
• If x = n, N (n) = n = h −1 ( 21 ) = ϕ−1 ( 21 ), or N (n) = ϕ−1 (1 − ϕ(x)
• If x ∈ (n, 1], or N (x) ∈ [0, n), ϕ(N (x)) + ϕ(x) = h(N (x)) + 1 − h(N (x)) = 1
In conclusion, N (x) = ϕ−1 (1 − ϕ(x)), for all x in [0, 1], or N = Nϕ .
Notice that the proof of last theorem shows clearly that the order-automorphism
ϕ such that N = Nϕ is not unique. Notice also that with ϕ = id[0,1] it follows
N (x) = 1 − x, the fundamental strong negation, with which it results N = Nϕ =
ϕ−1 ◦ (1 − id[0,1] ) ◦ ϕ = ϕ−1 ◦ N ◦ ϕ, that is, all strong negations belong to the
family of N0 (x) = 1 − x. Nevertheless, in all cases it is n = ϕ−1 ( 21 ) the fixed point
of Nϕ . √
If ϕ(x) = x 2 , it results Nϕ (x) = 1 − x 2 , called the circular negation. If ϕ(x) =
−1 −1 −1 1−x
1+x , or ϕ (x) = 2−x , it follows Nϕ (x) = ϕ (1 − 1+x ) = ϕ ( 1+x ) = 1+3x ,
2x x 2x 1−x
negations Nλ = 1+λx
1−x
(λ > −1) that only depends on one single parameter.
Remark 2.2.31 The only linear strong negation N is N = N0 , since from N (a) =
αa + β, with N (0) = 1 = β and N (1) = 0 = α + 1, follows α = −1 and
N (a) = 1 − a.
Remark 2.2.32 The only “rational” strong negations N of the form N (x) = ax+b cx+d ,
a, b different of 0, are those Nλ (λ > −1) in the Sugeno’s family. It follows from:
• N (0) = 1 = db , or d = b
c+d , or a = −b
• N (1) = 0 = a+b
that gives
−bx + b b(1 − x) 1−x
N (x) = = = .
cx + b cx + b 1 + bc x
2 For (x < y) is evident that ϕ(x) < ϕ(y) if either x, y ∈ [0, n], or x, y ∈ (n, 1]. If x ∈ [0, n],
y ∈ (n, 1] and x < y, since h(x) + h(N (x)) < 1, it is h(x) < 1 − h(N (x)), or ϕ(x) < ϕ(y).
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 59
As it was said, the only t-norm that is idempotent for all a in [0, 1], is T = min, and
the t-norms in { pr od} ∪ {W } only have the idempotents 0 and 1. As it was also said,
there are t-norms with several (but not all) idempotent elements. For example, the
function
1 1
+ W (3x − 1, 3y − 1), if (x, y) ∈ [ 13 , 23 ]2
T (x, y) = 3 3
min(x, y), otherwise ,
0.8
min
0.6
2/3
0.4 W
0.2 1/3
0
1
1/3 2/3
0.5
0.8 1
0 0.4 0.6
0 0.2
the only idempotent elements are b1 , b2 , b3 , b4 , etc., as well as 0 and 1, that is, the
points giving the partition of [0, 1].
Remark 2.2.38 Which t-norms are strictly non-decreasing in the sense that if 0 <
a < b < 1, then T (a, c) < T (b, c) for all c ∈ [0, 1]?
• If T = min, the answer is negative. For example, 0.3 < 0.5, but min(0.2, 0.3) =
min(0.2, 0.5) = 0.2
• If T = Wϕ , the answer is also negative. For example,0.3 < 0.5, but W (0.2, 0.3) =
W (0.2, 0.5) = 0
• If T = prodϕ , the answer is positive, since: a < b ⇒ ϕ(a) < ϕ(b) ⇒ ϕ(a) ·
ϕ(c) < ϕ(b) · ϕ(c) ⇒ ϕ−1 (ϕ(a) · ϕ(c))) < ϕ−1 (ϕ(b) · ϕ(c))), or pr odϕ (a, c) <
pr odϕ (b, c), because ϕ(c) ∈ (0, 1].
• If T is an ordinal-sum, it can’t be strictly non-decreasing because of the values it
takes with min.
Analogously, the only t-conorms that are strictly non-decreasing are those in
{ pr odϕ∗ }.
As it was said, in all standard algebras ([0, 1] X , T, S, N ) of fuzzy sets the triplet
(T, S, N ) share the following common properties:
1. T and S are commutative and associative
2. 1 is neutral for T , and 0 is neutral for S
3. 0 is absorbent for T , and 1 is absorbent for S
4. For all T and S, it is T min < max S
5. Each T (S) is non decreasing in the two variables
6. N is involutive, strictly decreasing and such that N (0) = 1,
a list of properties that gives some basic laws for fuzzy sets in the standard algebras,
like
• μ · σ = σ · μ, μ + σ = σ + μ,
• μ + (σ + λ) = (μ + σ) + λ = (σ + μ) + λ = λ + (σ + μ)
• μ · μ1 = μ, μ + μ0 = μ, μ + μ1 = μ1 , μ · μ0 = μ0
• If μ σ, then μ · λ σ · λ, and λ + μ σ + λ
• etc.
62 2 Algebras of Fuzzy Sets
Anyway, a lot of laws typical of classical sets are not always valid in all standard
algebras of fuzzy sets. For example, (P(X ), ∩, ∪, c ) is a Boolean algebra and no
one ([0, 1] X , T, S, N ) is a Boolean algebra. In particular, (P(X ), ∩, ∪) is a lattice
and the only standard algebra that is a lattice is that with T = min and S = max.
Let us study in which standard algebras some laws of crisp sets do hold.
Proof With b = c = 1, is T (a, S(1, 1)) = S(T (a, 1), T (a, 1)) or a = S(a, a). That
is S = max.
Provided S = max, the equation is T (a, max(b, c)) = max(T (a, b), T (a, c)). If
either b c or c b, it is immediate to check its validity for all t-norm T.
Hence,
• In all standard algebras with (T, max), it holds μ · (σ + λ) = μ · σ + μ · λ,
• In all standard algebras with (min, S), it holds μ + (σ · λ) = (μ + σ) · (μ + λ)
• The two distributive laws (2.3) and (2.4) do jointly hold if and only if T = min
and S = max.
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 63
μ · μ = μ0 ?
Proof With the fixed point n ∈ (0, 1) of N , it follows T (n, n) = 0, that is, T has zero-
divisors. Hence, T = Wϕ , and Wϕ (a, N (a)) = ϕ−1 (max(0, ϕ(a)+ϕ(N (a))−1)) =
0, or max(0, ϕ(a) + ϕ(N (a)) − 1) = 0, or ϕ(a) + ϕ(N (a)) − 1 0, that implies
ϕ(N (a)) 1 − ϕ(a), or N (a) ϕ−1 (1 − ϕ(a)) = Nϕ (a), for all a ∈ [0, 1]. Hence
N Nϕ . The reciprocal is a simple calculation.
64 2 Algebras of Fuzzy Sets
Proof With N (n) = n ∈ (0, 1), it follows S(n, n) = 1. That is S = Wψ∗ , and 1 =
Wψ∗ (a, n(a)) = ψ −1 (min(1, ψ(a) + ψ(N (a)))), or 1 = min(1, ψ(a) + ψ(N (a))).
Hence, 1 ψ(a)+ψ(N (a)), or Nψ (a) = ψ −1 (1−ψ(a) N (a)). That is, Nψ N .
The reciprocal is a simple calculation.
Proof If T = min, since a S(a, b), it follows min(a, S(a, b)) = a. With b = 0,
the equation gives T (a, a) = a, and = min.
Proof Since T (a, b) a, it follows max(a, T (a, b)) = a. With b = 1, the equation
gives S(a, a) = a, and = max.
Hence,
• The law μ · (μ + σ) = μ, holds for all S and T = min
• The law μ + (μ · σ) = μ, holds for all T and S = max
• The two laws hold jointly if and only if T = min and S = max.
With classical sets it always holds the law of von Neumann, or law of the perfect
repartition,
A = (A ∩ B) ∪ (A ∩ B c ),
A = (A ∪ B) ∩ (A ∪ B c ),
μ = μ · σ + μ · σ , μ = (μ + σ) · (μ + σ ),
66 2 Algebras of Fuzzy Sets
Lemma 2.2.46 The equation a = S(T (a, b), T (a, N (b))) holds if and only if T =
pr odϕ , S = Wϕ∗ , N = Nϕ .
It can be also proven that a = T (S(a, b), S(a, N (b))) if and only if T = Wϕ ,
S = pr odϕ∗ , N = Nϕ . Notice only that a = 0 gives T (b, N (b)) = 0, or T = Wϕ
and N = Nϕ .
Notice that the verification of von Neumann’s law require, in the case of fuzzy
sets, non-dual theories, like those given by the triplets ( pr odϕ , Wϕ∗ , Nϕ ), and
(Wϕ , pr odϕ∗ , Nϕ ).
The results in last section can be summarized in the following Table 2.1.
Hence, the algebras with the triplets (min, max, N ) are the ones that preserve more
structural Boolean properties. Indeed, these algebras preserve all the basic Boolean
laws except those of non-contradiction and excluded-middle. They are distributive
pseudo-complemented lattices, that is, De Morgan algebras that, in addition and like
all algebras of fuzzy sets, verify the law of Kleene,
for all a, b in [0, 1]. The algebras given by the triplets (min, max, N ) are De Morgan-
Kleene algebras.
It can be considered, in addition to the structural Boolean laws, the cases that can be
derived from them, for example,
(A ∩ B c )c = B ∪ (Ac ∩ B c ),
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 67
(μ · σ ) = σ + μ · σ
or
N (T (a, N (b))) = S(b, T (N (a), N (b))),
μ · μ + μ · μ = μ0 , []
there are no triplets (T, S, N ) for which it can hold S(T (a, a), T (a, N (a))) = 0 for
all a in [0, 1].
In the same vein, there are some laws that have solutions when different t-norms,
t-conorms and strong negations are considered. For example, (μ + μ) · (μ · μ ) = μ0 ,
that comes from (A ∪ A) ∩ (A ∩ Ac ) = ∅, translated in the form T1 (S(a, a),
T2 (a, N (a))) = 0, has infinite solutions like, for example, with an strong nega-
tion N , such that N N0 , T1 = min, T2 = W and any t-conorm S, since
min(S(a, a), T2 (a, N (a)) = T2 (a, N (a)) = W (a, N (a)) = max(0, a + N (a) −
1) = 0, because of T2 (a, N (a)) a S(a, a), and N (a) 1 − a, or
a + N (a) − 1 0.
Another case is given by the classical (derived) laws
A ∩ (Ac ∪ B) = A ∩ B, A ∪ (Ac ∩ B) = A ∪ B,
μ · (μ + σ) = μ · σ, μ + (μ · σ) = μ + σ,
T1 (a, S(N (a), b)) = T2 (a, b), S1 (a, T (N (a), b)) = S2 (a, b),
2.2.9 Examples
Which one of the following linguistic labels: cold, hot, warm, more or less cold,
more or less warm, is the more adequate for the temperatures of 20, 21 and 22 ◦ C?
Solution. With C = cold, it is μC (20) = 10 5
= 0.5 and μmor e or less C (20) =
√
0.5 = 0.71. To obtain H = hot, we can compute μ H as the opposite of μC :
⎧
⎨ 1, if 45 x 50
μ H (x) = μC (50 + 10 − x) = μC (60 − x) = 10 ,
x−35
if 35 x 45
⎩
0, if 10 x 35
graphically
μw = μcold · μhot
Then:
• x = 20, gives μc (20) = 5
10 = 0.5, μ H (20) = 0, μW (20) = 5
10 = 0.5,
• x = 21, gives μc (21) = 4
10 = 0.4, μ H (21) = 0, μW (21) = 6
10 = 0.6,
• x = 22, gives μc (22) = 3
10 = 0.3, μ H (22) = 0, μW (22) = 7
10 = 0.7,
70 2 Algebras of Fuzzy Sets
and
√
• μmor e or less cold (20) = √ μc (20) = 0.71,
μmor e or less war m (20) = μW (20) = 0.71
• μmor e or less ;cold (21) = 0.63, μmor e or less ;war m (21) = 0.77,
• μmor e or less ;cold (22) = 0.55, μmor e or less ;war m (22) = 0.84.
Hence
• The more adequate linguistic label for x = 20, cannot be decided but it could be
either ‘not cold’, or ‘not warm’. Since, it is not hot at all, we can take ‘not cold’.
• For x = 21, is ‘mol warm’ (mol = more or less)
• For x = 22, is ‘mol warm’
37 Age( p) 41,
and neither Age( p) 32, nor 43 Age( p). What can be said on the degree up to
which it could be Age( p) = 35, and Age( p) = 42?
Solution. What is unknown is the variation of Age( p) between 32 and 37, as well
as between 41 and 43. Since Age varies continuously, we can suppose there are two
functions
f
g Age ( p )
32 37 41 43
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 71
3 1
μ Age( p) (35) = , μ Age( p) (42) = .
5 2
As it will be seen later on, 0.6 is the possibility that Age( p) = 35, and 0.5 that
of Age( p) = 42. Hence, it seems a little bit more possible that it be ‘Age( p) = 42’
than ‘ Age( p) = 35’.
Example 2.2.49 Knowing that Height(John) = 175 cm, and Height(Peter) = 180 cm,
consider the two statements:
p = It is false that John is not very tall or is more or less short
q = It is false that Peter is not very tall or is more or less short.
which is more true?
Solution. Both statements can be written by
Is false that x is P,
with P = ‘(not very tall) or (more or less short)’.
Hence
μ P (x) = S(μver y tall (x), μnot shor t (x)) = S(N (μtall (x)2 ), μtall (A(x)) ),
with, perhaps, μtall (x) = 0.025x −3.75, x ∈ [150, 190], if we need to have numbers.
Hence, with A(x) = 210 − x, it is A(175) = 210 − 175 = 35, and μtall (35) = 0,
as well as A(180) = 220 − 180 = 30, and μtall (30) = 0, because of that
μ P (175) = S(N (μtall (175)), 0) = N (μtall (175)2 )
μ P (180) = S(N (μtall (180)), 0) = N (μtall (180)2 ).
Since μtall is strictly non-decreasing between 150 and 190, it is μtall (175) <
μtall (180), and N (μtall (180)2 ) < N (μtall (175)2 ). Finally,
N (N (μtall (175)2 )) < N (N (μtall (180)2 )), or μtall (175)2 < μtall (180)2 ,
μ + μ · σ = μ, μ · (μ + σ) = μ,
as well as that the negation is linear, Determine the triplet (T, S, N ) and, with X =
[0, 10], and ⎧
⎨ 0, if 0 x 5
x
μ(x) = , σ(x) = x−5 , if 5 x 8
10 ⎩ x−6
6
4 , if 8 x 10,
compute μ · σ, μ + σ, and μ + σ.
Solution. The first law of absorption μ + μ · σ = μ, implies S = max for
any T . The second law of absorption μ · (μ + σ) = μ, implies T = min for
any S. Hence (T, S) = (min, max), and the only linear N is N = N0 . Hence,
(T, S, N ) = (min, max, 1 − id). With the graphics of μ and σ in the figure,
μ′ μ
7 80:9 :
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 73
Nλ ,λ ∈(−1,0)
N0
Nλ ,λ ∈(0, + ∞)
Solution. It is ⎧
⎪ 1, if x ∈ [0, 3] ∪ [0, 4]
⎪
⎪
⎨ 0, if x ∈ [4, 6]
μ (x) = 4−x
⎪
⎪ x , if x ∈ [3, 4]
⎪
⎩ x−6
x , if x ∈ [6, 7],
Hence,
74 2 Algebras of Fuzzy Sets
Example 2.2.52 The negation is linear, and the standard algebra must verify the
law μ = μ · σ + μ · σ . Determine the triplet (T, S, N ), and with μ(x) = 10
x
(in
X = [0, 10]) and ⎧
⎪
⎨ 1, if x ∈ [0, 5]
σ(x) = 2 ,
7−x
if x ∈ [5, 7]
⎪
⎩
0, if x ∈ [7, 10]
⎧ ∗ x ⎫
⎨ W ( 10 , 1 − 10 ) = 1
x
⎪ ⎪
⎬
(σ · μ + σ · μ )(x) = W ∗ ( x(7−x)
20 , 2 (1 −
7−x
10 )
x
= 7−x = σ(x).
⎪
⎩ ∗
2 ⎪
⎭
W (0, 0) = 0
Finally, since,
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 75
⎧ ⎧
⎨1 ⎪
⎨0
(μ · σ ) (x) = 1 − (μ · σ )(x) = 1 − x(x−5) , and (μ · σ )(x) = x ) x−5 ,
(1 − 10
⎩ x
20 ⎪
⎩
2
1 − 10 1− x 10
it results
⎧ ∗
⎨ W (1, 0) = 1
⎪
(σ + (μ · σ ))(x) = W ∗ ( 7−x
2 , 2 (1 − 10 ) = 1 −
x−5 x x(x−5)
= (μ · σ ) (x) .
⎪
⎩ ∗
20
W (0, 1 − 10
x
) = 1 − 10
x
Example 2.2.53 Predicate F = high fever refers to the interval [37, 42] in a clini-
cal thermometer, in which the values {37, 37.5, 38, . . . , 41.5, 42} are significative.
Asking an expert one obtains the following fuzzy set
where it is clear that 0/37 + 0/37.5 + 0/38 is avoided since this values of the body’s
temperature are not significative for F. With all that, give the membership function
of P = very high fever, Q = more or less high fever, R = low fever, S = not high
fever.
√
Solution. With the usual definition μver y F = μ2F , μmol F = + μ F , μlow F =
μ F (37 + 42 − x) = μ F (79 − x), μnot F = 1 − μ F , it results:
• μ P = 0.09/38.5 + 0.25/39 + 0.49/39.5 + 0.64/40 + 0.81/40.5 + 1/41 +
1/41.5 + 1/42.
• μ Q = 0.55/38.5 + 0.7/39 + 0.84/39.5 + 0.89/40 + 0.95/40.5 + 1/41 +
1/41.5 + 1/42.
• μ R = 1/37 + 1/37.5 + 1/38 + 0.9/38.5 + 0.8/39.5 + 0.7/39.5 +
0.5/40.5 + 0.3/40.5.
• μ S = 1/37 + 1/37.5 + 1/38 + 0.7/38.5 + 0.5/39.5 + 0.3/39.5 +
0.2/40.5 + 0.1/40.5.
Notice the incoherence produced by μ S = μnot F μlow F = μ R . An incoherence
showing that it cannot be taken the representation μnot F = 1 − μ F , but some
μnot F = N ◦ μ F with N ≥ N0 .
For example, if N (x) = 1−0.9x
1−x
, it is
μ S = 1/37 + 1/37.5 + 1/38 + +0.96/38.5 + 0.91/39.5 + 0.81/39.5 +
0.71/40 + 0.53/40.5 + 1/41 + 1/41.5 + 1/42, showing μlow F μnot F .
Look that
μ F & low F = 0.3/38.5 + 0.5/39.5 + 0.7/39.5 + 0.5/40 + 0.3/40.5
provided μ F & low F = min(μ F , μlow F ).
76 2 Algebras of Fuzzy Sets
with convenient continuous t-norm T and t-conorm S. This formulas are the descrip-
tion of p and q in fuzzy terms.
For example, if a = 20, b = 50, μ P = 50−x 30 , if 20 x 50, and 40 − a =
30, 40 + b = 50, with μ A40 piece-wise linear, T = min, S = max, with
⎧
⎪
⎨ 0, if 0 x 50
μa P (x) = μ P (100 − x) = 30 ,
x−50
if 50 x 80
⎪
⎩
1, if 80 x 100,
the graphics is
The slashed function describes p, and the continuous one describes q. Of course,
Degree( p) min(μ P (x), μ A40 (x)) μ P (35) = 0.5.
Remark 2.2.55 To select T and S, the following points could be taken into account,
• It could be perfectly the case that ‘John is young and not young’ with a positive
degree. Hence, the laws of Non-contradiction can be avoided, and T ∈ / {W }.
2.2 The Concept of an ‘Algebra of Fuzzy Sets’ 77
• Since it is reasonable to accept that ‘John is old’ and ‘John is old’ does coincide
with ‘John is old’, and ‘John is young’ or ‘John is young’ does coincide with
‘John is young’ , we can decide to take either T = min, S = max, or T and S as
ordinal-sums.
• Provided idempotency is avoided i.e., ‘John is young’ and ‘John is young’ does
coincide with ‘John is very young’, instead of T = min, we can take T = prod,
that is more interactively than min. In this case, because it does not seem that
duality should be avoided, we could take √
S = Prod∗ , and
√
then
∗
μver y young or mol old (45) =Prod ( 36 , 6 ) = 6 (1 + 36 ) = 0.5636,
1 6 1 35 6
√
that is greater than the value 6
6
= 0.408 obtained with = max.
Hence,
• If N = N0 , the set is [0, 5], and the threshold (of selfcontradiction) of big is 5.
• If N = N1 , the set is [0, 4.142], and the threshold is 4.142
• If N = N2 , the set is [0, 3.66], and the threshold is 3.66
Notice that changing big by not big, the thresholds do remain but the sets are, respec-
tively, [5, 10], [4.142, 10], and [3.66, 10].
The kind of problems this section will deal with are like the following. An exam is
corrected by three referees R1 , R2 , R3 , each
3one with a different weight of strongness
W (Ri ) ∈ [0, 1], 1 i 3, such that i=1 W (Ri ) = 1. Each referee assigns a
numerical qualification pi ∈ [0, 10] to the exam delivered by a given student. How
these qualification can be “aggregated” to obtain final qualification for the student’s
exam? A recognized usual way of doing it is by the weighted mean:
1 p1 p2 p3
Q= · W (R1 ) + · W (R2 ) + · W (R3 ),
10 10 10 10
pi
with 10 ∈ [0, 10]. For example, if W = (0.5, 0.3, 0.2) and P = (7, 6, 5), it follows
78 2 Algebras of Fuzzy Sets
1
Q = 0.7 × 0.5 + 0.6 × 0.3 + 0.2 × 0.5 = 0.63
10
that implies Q = 6.3. Provided the three referees have the same weight, it is W =
( 13 , 13 , 13 ), and then, Q = p1 · 13 + p2 · 13 + p3 · 13 = p1 + p32 + p3 = 7+6+5
3 = 6, is
just the arithmetic mean of the three qualifications.
Another way of obtaining the final qualification, this time by ignoring the referee’s
character, is by the geometric mean
√ √
Q= p1 · p2 · p3 = 7 × 6 × 5 = 5.94,
3 3
showing
√ that in a problem with p1 = p2 = 10, p3 = 0, it results Q =
3
10 × 10 × 0 = 0, when the arithmetic mean is 20
3 = 6.67.
such that
1. A is continuous in all variables
2. A(0, . . . , 0) = 0, and A(1, . . . , 1) = 1
3. If x1 y1 , . . . , xn yn , then A(x1 , . . . , xn ) A(y1 , . . . , yn ).
Sometimes it is said that A is an n-dimensional aggregation function. Continuous
t-norms and continuous t-conorms are 2-dimensional aggregation functions.
Of the many types of aggregation functions, a particular and important type are
the quasi-linear means,
n
M(x1 , . . . , xn ) = f −1 pi . f (xi )
i=1
n
with ( p1 , . . . , pn ) in [0, 1], verifying pi = 1, and f : [0, 1] → R, continuous,
i=1
one-to-one, and monotonic. Function f is called the generator of M.
Notice that if f is the identity f (x) = x, we get the weighted means:
n
M(x1 , . . . , xn ) = pi · xi ,
i=1
that with p1 = 1
n (1 i n) is the arithmetic mean
2.3 On Aggregating Imprecise Information 79
1
n
M(x1 , . . . , xn ) = xi ,
n
i=1
Because they are associative, continuous t-norms and continuous t-conorms can be
extended to n-dimensional aggregation functions. For example, with n = 3,
2.3.4 Examples
The pointwise aggregation of classical sets is not, in general, a classical set, but a
fuzzy one. For example, the arithmetic mean verifies
1
M(0, 0) = 0, M(0, 1) = M(1, 0) = , M(1, 1) = 1
2
and, if A, B are crisp subsets, M(A, B) is not a crisp subset if given by M(μ A , μ B )
(x) = M(μ A (x), μ B (x)). On the contrary, with the geometric mean G, it is
Solution.
M(μ, σ) = 0.75/1 + 0.35/2 + 0.5/3 + 0.85/4 + 0.4/5
G(μ, σ) = 0.735/1 + 0/2 + 0.5/3 + 0.837/4 + 0/5
O(μ, σ) = (0.4 × 0.6 + 0.6 × 0.9)/1 + (0.4 × 0 + 0.6 × 0.7)/2 + (0.4 ×
0.5 + 0.6 × 0.5)/3 + (0.4 × 0.7 + 0.6 × 1)/4 + (0.4 × 0 + 0.6 × 0.8)/5 =
0.72/1 + 0.42/2 + 0.5/3 + 0.88/4 + 0.48/5.
Notice that G(μ, σ) M(μ, σ), but that neither G(μ, σ) and O(μ, σ), nor M(μ, σ)
and O(μ, σ), are order-comparable.
82 2 Algebras of Fuzzy Sets
Example 2.3.2 A linguistic variable has the fuzzy values H = high, S = short and
M = medium, with H represented by
⎧
⎨ 0, if 0 x 5
μ H (x) = 1, if 7 x 10
⎩ x−5
2 , if 5 x 7.
The solution for the first supposition appears in the following sequence of figures
(with · = min, = 1 − id).
graphically,
0.7
M
0.3
3 5 7
Chapter 3
Reasoning and Fuzzy Logic
The question is, in fact, a philosophical one whose discussion does not correspond
to this text, and that received a lot of comments and discussions by philosophers.
Instead of such question, there is the more particular,what is a logic?, that can be
answered not philosophicaly but in terms of the mathematical definition of what is a
consequence’s operator. A definition that corresponds to an abstraction of the term
“deduction”.
Remark 3.1.2 Instead of a lattice, let us take the set [0, 1] X endowed with a fuzzy
intersection
μ∧ = μ1 · · · μn = T ◦ (μ1 × · · · × μn )
(T a continuous t-norm), the partial order μ σ ⇔ μ(x) σ (x), for all x ∈ X , and
the empty set μ0 = μ∅ . Take the set P0 ([0, 1] X ) that consists of the finite subsets
P = {μ1 , . . . , μn } ⊂ [0, 1] X such that μ∧ = μ0 . The definition
Cons(P) = {σ ∈ [0, 1] X ; μ∧ σ },
Example 3.1.3 Let it be the set X = {1, 2, 3, 4, 5} endowed with the preorder
1 2 3 4 5
It is, for example C({1, 5}) = {1, 5, 2, 3, 4}, C({4, 5}) = {4, 5}, C({2, 3, 4}) =
{2, 3, 4, 5}, C({1, 2, 4, 5}) = {1, 2, 4, 5, 3} = X , C({5}) = {5}, and C({3}) =
{3, 4, 5}.
Hence, it is possible to identify a logic in a set with a preordering of it.
3.1.2 Conjecturing
P σ ⇔ σ ∈
/ C(P) ⇔ N ◦ σ ∈
/ C(P),
P σ ⇔ ∀μ ∈ P : μ σ .
In what follows, we will only take into account the operator Cons(P) = {σ ; μ∧
σ }, provided μ∧ = μ0 . Consequently, the set of the conjectures that, through
Cons, is associated to any finite set P = {μ1 , . . . , μn } of premises such that
μ∧ = μ1 · · · μn = T ◦ (μ1 × . . . × μn ) = μ0 , for some continuous t-norm T , is
where
• H yp(P) = {σ ∈ Con j (P); μ0 < σ < μ∧ }, is the set of hypotheses of P.
• Sp(P) = {σ ∈ Con j (P); μ∧ N C σ }, is the set of speculations of P, where NC
means non-comparable under ,
and verifying
3.1 What Does It Mean “Logic”? 87
Notice that
• Consequences follow (in the partial order ) from μ∧ : all the premises explain
the consequences.
• Hypotheses explain all the premises, since σ < μ∧ μi , all the premises follow
from each hypothesis.
• Speculations, are the conjectures for which it is neither μ∧ σ , nor σ < μ∧ .
With all that, processes,
• P σ , with σ ∈ Con j (P), is a guessing, or conjectural reasoning
• P σ , with σ ∈ Cons(P), is a deduction, or deductive reasoning
• P σ , with σ ∈ H yp(P), is an abduction, or abductive reasoning
• P σ , with σ ∈ Sp(P), is an speculation, or speculative reasoning
Concerning this types of conclusive reasonings, it should be pointed out what follows.
1. The set Con j (P) is not always in A = P0 (X ). Hence, it can’t be taken as a set
of premises, it has not sense to consider Con j (Con j (P)), or Cons(Con j (P)).
2. If P ⊂ Q, it is Con j (Q) ⊂ Con j (P), because if σ ∈ Con j (Q), from I n f Q
I n f P, if I n f Q σ , it is I n f P σ . Hence, if σ ∈ Con j (Q) it is σ ∈
Con j (P), and Con j is anti-monotonic.
3. The sets H yp(P) and Sp(P) are not always in P0 (X ). Hence, they can’t be taken
as sets of premises.
4. If P ⊂ Q, it is H yp(Q) ⊂ H yp(P), since σ ∈ H yp(Q), or μ0 < σ < I n f Q
implies μ0 < σ < I n f P. Hence, the operator H yp is anti-monotonic.
5. Concerning the operator Sp, if P ⊂ Q it can be Sp(P) Sp(Q), and Sp(Q)
Sp(P). Hence, Sp is a non-monotonic operator.
If x is P, then y is Q,
where x ∈ X , y ∈ Y , P is a predicate (precise or imprecise) in X , and Q is a predicate
(precise or imprecise) in Y . For example, with X = [0, 1], Y = [0, 10],
If x is small, then y is big,
or,
If x is small and y is big, then z is not small,
with x, y ∈ [0, 1] and z ∈ [0, 10].
What it means to represent a conditional statement like “If x is P, then y is Q”?
It means to translate it in fuzzy terms. For example, ‘x is P’ and ‘y is Q’ will be
translated by μ P (x) and μ Q (y) with adequate fuzzy sets μ P , μ Q ∈ [0, 1] X , adequate
in the sense that they capture the use of P on X and Q on Y .
But how to represent the full statement “If x is P, then y is Q”:= μ P (x) →
μ Q (y)? It is always done, in fuzzy logic, by means of a function J : [0, 1]×[0, 1] →
[0, 1], such that
for all x ∈ X , and y ∈ Y . But, which function J can be taken? It depends on the
‘meaning’ of the conditional statement, and this requires to look at what happens in
general.
Remark 3.2.1 Imprecise conditionals are very useful in ordinary life, for instance,
• If the turn is not so far, and the car’s speed is not high, then press softly the break,
• If the the food was well-cooked and of quality, and the service was good, the the
tip should be higher than 15 %.
Both in common life and in technology, a lot of imprecise conditionals (rules) are
considered. The need for its representation will be obvious in the next section.
Let us consider, in the first place, the case in which the statements are crisp and
belong to a Boolean algebra (B, ·, +, ; 0, 1), where ·, +, stand, respectively, for
the intersection (and), the union (or), and the complement (negation, not), 0 is the
minimum element, and 1 is the maximum. As it is well known, this Boolean algebra
is naturally ordered by means of the partial order
a b ⇔ a · b = a ⇔ a + b = b ⇔ b = a + a · b.
because of:
1. a · z b ⇒ a + a · z = a + z a + b, and since z a + z, follows z a + b
2. z a + b ⇒ a · z a · (a + b) = a · b b.
Then, a · (a → b) b ⇔ a → b a + b.
A conditional function is a mapping →: B × B → B, such that a · (a → b) b
for all a, b ∈ B. Hence, in a Boolean algebra, the biggest conditional is a + b, the
so-called material conditional, and any smaller function is also a conditional. For
example, from
a · b b a + b, a a + b
a · (b + b ) + a · b, a + a · b, b + a · b
since a · (b + b ) + a · b = a + a · b = (a + a ) · (a + b) = a + b, and
b + a · b = (b + a ) · (b + b ) = a + b.
Notice that from z 1 a + b, z 2 a + b, follows z 1 · z 2 (a + b) · (a + b) =
a + b, z 1 + z 2 (a + b) + (a + b) = a + b, hence, the union and the intersection
Remark 3.2.2 If the Boolean algebra B is complete, that is, for any A ⊂ B, A = ∅,
it exists Sup A ∈ B, then
Sup{z ∈ B; a · z b} = Sup{z ∈ B; z a + b} = a + b.
coincidental with a + b. The Sasaki and the Dishkant hooks are used as models for
the conditional statements in the reasoning in Quantum Logic.
Remark 3.2.5 The scheme of Modus Ponens
If a, then b
a
b,
corresponds to forwards reasoning, that is, goes from the antecedent a to the con-
sequent b thanks to the conditional a → b, through a · (a → b) b. Backwards
reasoning goes from the consequent to the antecedent (also thanks to a → b), ad it
is modeled by the Modus Tollens scheme.
If a, then b
not b
not a,
Let us return to the case of fuzzy logic, that is, to a conditional linguistic expression,
or rule, like ‘If x is P, then y is Q’, represented in fuzzy terms by
3.2 Reasoning with Conditionals: Representation 91
for all x ∈ X , y ∈ Y . The problem is which function J : [0, 1] × [0, 1] → [0, 1],
is to be taken at each case if, of course, it gives a conditional. That is, if from the
premises {‘x is P’, ‘If x is P, then y is Q’} follows ‘y is Q’ as a logical consequence.
Formally speaking, it should exist a continous t-norm T0 such that
for all x ∈ X , y ∈ Y . This condition, that should hold for any μ P (x) ∈ [0, 1], and
any μ Q (y) ∈ [0, 1], conducts to the inequality
for all a, b in [0, 1], J should verify to represent conditional statements. It is called
the Modus Ponens Inequality, since it allows the scheme of reasoning
If x is P, then y is Q
x is P
y is Q
T0 (a, J (a, b)) b ⇔ J (a, b) JT0 (a, b) = sup{z ∈ [0, 1]; T0 (z, a) b}.
Hence, for each T0 , the greatest T0 -conditional is the function JT0 , since, it verifies
the MP-inequality T0 (a, JT0 (a, b)) = min(a, b) b.
Remark 3.2.6 For reasons that will be latter on presented, T -conditionals JT are
called R-implications (R shorting residuated). They come directly from the Boolean
equation a + b = sup{z; a · z b}.
that coincides with the values of max(1 − μ(x), σ (y)). That is, all R-implications
do coincide with the Boolean material conditional μ + σ in the case that μ and σ
92 3 Reasoning and Fuzzy Logic
are crisp sets. R-implications generalize the material conditional. Of course, this will
happen with any J such that
Example 3.2.8 Let us see how is JT , when T is, respectively, the continuous t-norm
min, pr odϕ , Wϕ .
1, if a b
• T = min, Jmin (a, b) = sup{z ∈ [0, 1]; min(z, a) b} = (G ödel
b, if a > b
implication).
T = pr odϕ , JT (a, b) = sup{z ∈ [0, 1]; ϕ(a).ϕ(z) ϕ(b)} =
•
1, if a b
ϕ(b) (Goguen implication).
ϕ −1 ( ϕ(a) ), if a > b
• T = Wϕ , JT (a, b) = sup{z ∈ [0, 1]; ϕ −1 (W (ϕ(a), ϕ(z)) b} = ϕ −1 (min(1, 1−
ϕ(a) + ϕ(b))) (Łukasiewicz implication).
Since each JT is a T -conditional, G ödel’s is a min-conditional, Goguen’s are prodϕ -
conditionals, and Łukasiewicz’s are Wϕ -conditionals. Notice that the S-implications
of the form
are exactly the Łukasiewicz’s R-implications: the only R-implications that are
S-implications are the Łukasiewicz’s ones. If, for instance, it were
3.2 Reasoning with Conditionals: Representation 93
1, if a b
Jmin (a, b) = S(N (a), b) =
b, if a > b
it will result
1, if N (a) b
S(a, b) =
b, if N (a) > b
a function that is not a t-conorm, since S(a, 0) = 0 = a, if a > 0. Hence, Jmin is not
an S-implication. An analogous reasoning shows that J pr odϕ are not S-implications.
with which J2 (N (b), N (a)) = S(N (a), T (a, b)) = J1 (a, b) or, equivalently,
J2 (a, b) = J1 (N (b), N (a)): D-operators are the contrasymmetricals of Q-operators.
Hence, it can be repeated all that has been said for J1 . For example,
If S = max, T = min, N = N0 , it is J2 (a, b) = J1 (1 − a, 1 − b) =
max(b, min(1 − b, 1 − a)) = max(b, 1 − max(a, b)), that verifies
94 3 Reasoning and Fuzzy Logic
functions with the inconvenience of the property J (a, b) = J (b, a), but verifying,
that is, all of them are conditionals for any t-norm T0 and in particular, for the greatest
of them. They are always taken as min-conditionals. For example,
• If T = min, J (a, b) = min(a, b), is called the Mamdani conditional
• If T = prod ϕ , J (a, b) = ϕ −1 (ϕ(a) · ϕ(b)), are called Larsen conditionals
• T = Wϕ is never used, since it can be J (a, b) = 0 with a > 0 and b > 0.
For example, with ϕ(x) = x(1+x)
(an order automorphism), it is ϕ −1 (x) =
√ 2
8x+1−1
2 , and
J (a, b) = T (a r , bs ),
Example 3.2.13 Once given a conditional statement ‘If x is P, then y is Q’, and
represented ‘x is P’ by μ P (x), and ‘y is Q’ by μ Q (y), it remains to be understood
what it is meant by the ‘statement’ μ P (x) → μ Q (y). It could be, or not to be,
μ P (x) → μ Q (y) = (μ P → μ Q )(x, y), with μ P → μ Q a fuzzy set in X × Y ,
identified with some expression involving the connectives and (·), or (+), not ( ).
In the affirmative case, it is said that μ P → μ Q is expressible in material form, for
example, μ P → μ Q = μP +μ Q , or μ P → μ Q = μP +μ P ·μ Q , etc. These material
3.2 Reasoning with Conditionals: Representation 95
forms are called protoforms; for instance, the protoform of the Kleene-Diennes fuzzy
conditional is a + b, and that of Mamdani fuzzy conditional is a · b.
If μ P → μ Q does not correspond with a protoform, one can try to represent it
by means of an R-implication, that is, by Jmin or by some J pr odϕ , since all the JWϕ
do correspond to a protoform a + b, or μP + μ Q , with + represented by Wϕ∗ and
by Nϕ .
In addition, there is a problem that should be taken into account when representing
μ P → μ Q . The problem is the following. Suppose that we know μ P → μ Q
should be represented by a function J that is a min-conditional, but that we are
not able to decide a protoform and we take Jmin . Since J Jmin , we will reach
the biggest possible output. This should be known. Analogously, if J should be a
prod-conditional, from J J pr od , follows the same comment.
Example 3.2.14 Let’s stop for a while in the above mentioned concept of implication
function, a concept that comes directly from the properties shown by the Boolean
conditional a → b = a + b, whose truth value is usually represented by
that is, for instance 0.1 < 0.4 but J D (0.3, 0.1) > J D (0.3, 0.4).
96 3 Reasoning and Fuzzy Logic
{μ, μ → σ } σ,
or backwards,
{σ , μ → σ } μ .
The first type of inference corresponds to search for the solutions of μ·(μ → σ ) σ ,
that is,for J and T0 such that,
Hence, given J , we need to know T0 such that T0 (a, J (a, b)) b, for forward
inference, and given J and N , we need to know T1 such that T1 (N (b), J (a, b))
N (a), for backwards inference.
Notice that the two t-norms in (∗ ) and (∗∗ ) are not necessarily coincidental. For
example, given J (a, b) = max(1 − a, b), with N = N0 , can we do backwards
inference? To answer this question we just need to know if there is a continuous
t-norm T1 such that T1 (1 − b, max(0, max(1 − a, b)) 1 − a, with a = 1 it results
T1 (1 − b, b) = 0, and T1 = W . Then, since,
W (1 − b, max(1 − a, b)) = max(0,
1 − b + max(1 − a, b)− 1)
1 − a − b, if a + b 1
= max(0, max(1 − a − b, 0)) = 1 − a, because
0, if a + b > 1
b 0 implies −b 0, and 1 − a − b 1 − a. Finally, the answer is: Yes, with
T1 = W . It can also be done backwards inference in the following cases,
1 − b, if a b
1. With Jmin , since W (1−b, Jmin (a, b)) = 1−a.
W (1 − b, b) = 0, if a > b
3.2 Reasoning with Conditionals: Representation 97
1 − b, if a b
2. With J pr od , since W (1 − b, J pr od (a, b)) = 1 − a.
a , if a > b
b 1−b
3. With JW , since W (1 − b, min(1, 1 − a + b)) = min(1 − b, 1 − a) 1 − a.
4. With J (a, b) = 1 − a + ab, since W (1 − b, min(1, 1 − a + ab)) =
(1 − b)(1 − a) 1 − a.
5. With J (a, b) = max(1−a, min(a, b)), since W (1−b, max(1−a, min(a, b)))
W (1 − b, max(1 − a, b)) 1 − a.
6. With J (a, b) = pr od ∗ (1 − a, a · b) = 1 − a + a 2 b, since W (1 − b, pr od ∗ (1 − a,
a · b)) W (1 − b, pr od ∗ (1 − a, b)) = W (1 − b, 1 − a + a · b) 1 − a.
Remark 3.2.17 The name R-implication, or residuated implication, comes from the
idea of ‘residuum’ that clearly appear in the case of J pr od when
b
If a > b, then J pr od (a, b) = .
a
Remark 3.2.18 In the same vein under which it was proven that R-implications J
with T = Wϕ are not S-implications, it is easy to show that they are not expressible
in material protoform, that is, by an expression with logical connectives. Take the
perhaps more general material protoform μ · (σ + σ ) + μ · σ . Is it possible that
or, S1 (N1 (a), 0) = 0. This means that S1 is not a t-conorm. Hence, the decision of
representing an R-implication can’t be taken from a material protoform interpretation
of it.
98 3 Reasoning and Fuzzy Logic
{μ, μ → σ } σ
is classically called Modus Ponendo Ponens (from the Latin, mode of starting the
truth (of σ ) by placing the truth (of μ)), or, for short Modus Ponens. That given by
the scheme {σ , μ → σ } μ is classically called Modus Tollendo Tollens (from
the Latin, mode of stating the falsity (of μ) by placing the falsity (of σ )), or, for
short Modus Tollens. They correspond to what we called forwards and backwards
reasoning. But there are again other modes of reasoning that can be considered, for
example,
• Modus Tollendo Ponens, given by the scheme {μ , μ+σ } σ and also called Mode
of Disjunctive Reasoning, and classically proven by μ ·(μ+σ ) = μ·μ +μ ·σ =
μ · σ μ (in a Boolean algebra).
• Modus Ponendo Tollens, given by the scheme {μ, (μ·σ ) } σ , classically proven
by μ · (μ · σ ) = μ · (μ + σ ) = μ · σ σ (in a Boolean algebra).
• Constructive Dilemma, given by the scheme {μ + λ, μ → σ, λ → η} σ + η,
classically proven by (μ + λ) · (μ → σ ) · (λ → η) = (μ + λ) · (μ + σ ) · (λ + η) =
μσ λ +ηλμ +ημσ +ηλσ σ +η (in a Boolean algebra where a → b = a +b).
• Destructive Dilemma, given by the scheme {μ + σ , λ → μ, η → σ, } λ + σ ,
classically proven by (μ + σ )(λ + μ)(η + σ ) = λ (μ · η + μ σ ) + η (σ · λ +
μ · σ ) λ + η (in a Boolean algebra where a → b = a + b).
What in the fuzzy case? For example, in the case of the Disjunctive Mode we
need to find all the possibilities for μ · (μ + σ ) σ , that is, to solve the functional
equation
Remark 3.3.1 It should be pointed out that the Modus Ponendo Tollens (MPT) can
be reduced, in the case of duality, to the disjunctive mode by means of the change
μ = α , σ = β , in which case since (μ · σ ) = μ + σ it follows μ · (μ · σ ) =
μ · (μ + σ ) = α · (α + β) β = σ . Hence, it holds with the triplet (Wϕ , Nϕ , Wϕ∗ ).
A central topic fuzzy logic deals with are non-rigid, dynamic systems involving
‘variables’ x1 , . . . , xn , y taking values in, respectively, universes X 1 , . . . , X n , Y ,
and constrained by imprecise rules ri of the type
P* Q*
is the greatest function verifying the GMP (3.1). This formula is known as the Com-
positional Rule of Fuzzy Inference (CRI, for short), and was introduced by Lotfi A.
Zadeh as the output fuzzy logic considers in the systems that are described by rules. It
is not to be forgotten that T0 is the continuous t-norm that makes J a T0 -conditional.
Sometimes, the input is just numerical, crisp in the form ‘x = x0 ’, that is, ‘x is
P ∗ ’ is ‘x is x0 ’, or ‘x ∈ {x0 }’ and then μ P ∗ = μ{x0 } , with
1, if x = x0
μ{x0 } (x) =
0, if x = x0 .
In this case,
Example 3.4.1 Take X = [0, 1], Y = [0, 10], and the rule ‘If x is small, then y is
big’, with the observation that ‘x is big’ and J (a, b) = max(1 − a, b). With the
y
membership functions μ P (x) = 1 − x, μ Q (y) = 10 , μ P ∗ (x) = x, it results
y y
μ Q ∗ (y) = sup W (x, max(x, )) = W (1, max(1, )) = 1,
x∈[0,1] 10 10
Example 3.4.2 With the same rule of last example and the input x0 = 0.5, it is
y
μ Q ∗ (y) = max(0.5, max(0.5, )), ∀y ∈ [0, 1]
10
graphically,
Q*
Example 3.4.3 With the rule ‘If x is small, then y = 8’, and the input x = 0.5, is
1, if y = 8
μ Q ∗ (y) = max(0.5, μ{8} (y)) = ,
0.5, if y = 8
graphically,
102 3 Reasoning and Fuzzy Logic
Remark 3.4.4 In general, it is not easy to assign a name to the functional output μ Q ∗ ,
that is, to express Q ∗ linguistically. In Example 3.4.2, it could be said Q ∗ = big
after 5 and constantly equal to 0.5 before 5. In Example 3.4.3, it could be said Q ∗ =
almost always 0.5.
Let’s consider the particular case where both universes X and Y are finite sets. If
X = {x1 , . . . , xn }, Y = {y1 , . . . , ym },
= max T0 (ri∗ , ai j ), 1 i n,
1i n
s ∗j = max T0 (ri∗ , ai j ), 1 j m,
1i n
in which the elements of the classical product of matrices (rows by columns) ri∗ ·
1i n
ai j are substituted by max T0 (ri∗ , ai j ).
1i n
This composition is called the max-T0 product of matrices, instead of the classical
sum-prod composition. Hence,
[μ Q ∗ ] = (s1∗ , . . . , sm∗ ) = [μ P ∗ ] ⊗ [J ],
a11 = J (0.7, 0.9) = min(1, 1 − 0.7 + 0.9) = 1; a12 = J (0.7, 0.6) = 0.9; a13 =
J (0.7, 0) = 0.3; a14 = J (0.7, 0.8) = 1
a21 = J (0.8, 0.9) = 1; a22 = J (0.8, 0.6) = 0.8; a23 = J (0.8, 0) = 0.2; a24 =
J (0.8, 0.8) = 1
a31 = J (1, 0.9) = 0.9; a32 = J (1, 0.6) = 0.6; a33 = J (1, 0) = 0; a34 =
J (1, 0.8) = 0.8.
Hence, ⎛ ⎞
1 0.9 0.3 1
[μ Q ∗ ] = (0.6 0.7 1) ⊗ ⎝ 1 0.8 0.2 1 ⎠ = (0.9 1 0.9 0.8),
0.9 0.6 0 0.8
since: (max(W (0.6, 1), W (0.7, 1), W (1, 0.9)), max(W (0.6, 0.9), W (0.7, 0.8),
W (1, 0.6)), max(W (0.6, 0.3), W (0.7, 0.2), W (1, 0)), max(W (0.6, 1), W (0.7, 1),
W (1, 0.8))) = (max(0.6, 0.7, 0.9), max(1, 1, 0.6), max(0.9, 0.9, 0),
max(0.6, 0.7, 0.8)) = (0.9 1 0.9 0.8). That is
Actually, there are no systems described by a single rule. What to do when a system
is described by, at least, two rules? With, for example
• r1: If x is P1 , then y is Q 1
• r2: If x is P2 , then y is Q 2 ,
an input μ P ∗ gives
104 3 Reasoning and Fuzzy Logic
and analogously for more than two rules. For example, in the case of p rules
r1 , . . . , r p , the result will be
μ Q ∗ = max(μ Q ∗1 , . . . , μ Q ∗p ),
where μ Q i∗ (1 i p) is the output obtained with the rule ri and the input μ P ∗ .
Example 3.4.6 With X = [0, 1], Y = [0, 10], consider the rules
• If x is big, then y = 2
• If x is small, then y = 8
• If x is ar ound 0.5, then y = 6,
and the input x0 = 0.4. Which is the final output of this system if the rules are
represented by J (a, b) = a · b?
0.4, if y = 2
• μ Q ∗1 (y) = J (μ B (0.4), μ{2} (y)) = 0.4.μ{2} (y) =
0, if y = 2
0.6, if y = 8
• μ Q ∗2 (y) = J (μs (0.4), μ{8} (y)) = 0.6.μ{2} (y) =
0, if y = 8
• μ Q ∗3 (y) = J (μ A0.5 (0.4), μ{6} (y)) = μ A0.5 (0.4).μ{6} (y)
μ A0.5 (0.4), if y = 6
=
0, if y = 6,
with μ B (x) = x, and μ S (x) = 1 − x. Taking as μ A0.5 the triangular function
and since the left side equation is y = x−0.35
0.15 , it is μ A0.5 (0.4) = 0. 3. Hence,
μ Q ∗3 (y) = 0. 3, if y = 6
0, if y = 6,
3.4 Inference with Fuzzy Rules 105
Finally
⎧
⎪
⎪ 0.4, if y = 2
⎪
⎨
μ Q ∗ (y) = max(μ Q ∗1 (y), μ Q ∗2 (y), μ Q ∗3 (y)) = 0. 3, if y = 6
⎪
⎪ if y = 8
⎪ 0.6,
⎩
0, otherwise
Q*
Hence, the numerical output that corresponds to the input x0 = 0.4, is y0 = 5.7.
Remark 3.4.7 Notice that once a system of rules linguistically describing the
behavior of a system is given, and where the consequents of the rules are numerical,
at each numerical input x0 in X does correspond a numerical output y0 in Y . In that
way, a function CRI: X → Y is defined. As it will be later on commented, were
the system’s behaviour previously known by a continuous function f : X → Y , the
function CRI approaches, under some additional conditions, the function f .
Remark 3.4.8 Look how important is to properly select the T-conditionals repre-
senting the rules.
Given the rule ‘If x is small, then y is big’, with X = Y = [0, 1], and μ S (x) =
1 − x, μ B (y) = y, J (a, b) = max(1 − a, b),
it follows J (μ S (x), μ B (y)) = max(x, y), that could be interpreted as ‘x is big or y
is big’.
With J (a, b) = min(1, 1 − a + b), it follows J (μ S (x), μ B (y)) = min(1, x + y) =
W ∗ (x, y), also interpretable as ‘x is big or y is big’. But with J (a, b) = min(a, b),
is J (μ S (x), μ B (y)) = min(1 − x, y), interpreted as ‘x is small and y is big’.
106 3 Reasoning and Fuzzy Logic
3.4.3 Examples
Example 3.4.9 Rule ‘If x is big, then y = 0.8’, with x, y in [0, 1], and the observation
x ∈ [0.4, 0.6]. Hence:
x, if y = 0.8
J (μ B (x), μ{0.8} (y)) = x · μ{0.8} (y) =
0, if y = 0.8
0.6, if y = 0.8
Then, μ Q ∗ (y) = Sup min(μ[0.4,0.6] (x), xμ{0.8} (y)) =
x∈[0,1] 0, if y = 0.8,
since,
Example 3.4.10 Rule: ‘If x is big, then y is small’, with the same observation as that
in the last example and with μ B (x) = x, μ S (y) = 1 − y, and J (a, b) = min(a, b),
follows:
μ Q ∗ (y) = Sup min(μ[0.4,0.6] (x), min(x, 1 − y)) = Sup min(min(μ[0.4,0.6] (x),
x∈[0,1] x∈[0,1]
x), 1 − y) = Sup min(x, 1 − y) = min(0.6, 1 − y).
x∈[0,1]
Example 3.4.11 X = {1, 2, 3}, Y = {6, 7}. Rule: ‘If x is around 2, then y = 6’, and
μ P ∗ (x) = 0.6/1 + 0.9/2 + 0.7/3, with ⎛ J (a, ⎞
b) = ab (Larsen). It results
0.5 0
(μ Q ∗ (6), μ Q ∗ (7)) = (0.6 0.9 0.7) ⊗ ⎝ 1 0 ⎠ = (0.9 0.7), that is μ Q ∗ = 0.9/6 +
0.5 0
0/7.
µ Q*
µ Q1* µ Q2*
Finally, ⎧
⎨ 0.4(1 − y), if 0 y 0.463
μ Q ∗ (y) = max(μ Q ∗1 (y), μ Q ∗2 (y)) = y 2 , if 0.463 y 0.6
⎩
0.36, if 0.6 y 1.
Example 3.4.13 Let’s find the function CRI:X → Y , in the case with X = [0, 1],
Y = [0, 1], and
• r1: If x is small, then y = 9
• r2: If x is big, then y = 2,
1 − x, if y = 9
it follows μ Q ∗1 (y) = (1 − x)μ{9} (y) =
0, if y = 9,
x, if y = 2
μ Q ∗2 (y) = xμ{2} (y) =
0, if y = 2,
and
⎧
⎨ x, if y = 2
μ Q ∗ (y) = max(μ Q ∗1 (y), μ Q ∗2 (y)) = 1 − x, if y = 9 that gives,
⎩
0, otherwise
2x + 9(1 − x)
CRI(x) = = 9 − 7x,
x +1−x
108 3 Reasoning and Fuzzy Logic
as the “theoretical” (linear) behavior of the system (x, y). For each x0 ∈ X , the value
CRI(x0 ) ∈ Y is the defuzzified value that corresponds to x0 .
This theorem (whose proof is here avoided) is simply an existential one, since
there is no general method for obtaing neither a fuzzy representation of the system
of rules, not the defuzzifiction method. It simply shows that it is possible to find a
CRI approaching enough well f for all points in [a, b].
3.5 Deffuzification
How to defuzzify non discrete outputs μ Q ∗ ? Let us proceed with two examples
without computational difficulties.
1st Example. Rules,
• r1: If x is big, then y is small
• r2: If x is small, then y is big
with X = [0, 1], and Y = [0, 10]. Take,
y y
μ B (x) = x, μ S (y) = 1 − , μ S (x) = 1 − x, μ B (y) = ,
10 10
and J (a, b) = min(a, b) -Mamdani-. Notice that, with the observation x 0 = 0.5,
y y
μ Q ∗1 (y) = min(0.5, 1 − ), μ Q ∗2 (y) = min(1 − 0.5, ).
10 10
3.5 Deffuzification 109
y y y
Then, μ Q ∗ (y) = max(min(0.5, 1 − 10 ), min(0.5, 10 )) = max(0.5, min(1 − 10 ,
y y y
10 )) = 0.5, since min(1 − ,
10 10 ) 0.5.
The area below μ Q ∗ (y) = 0.5, is A = 0.5 × 10 = 5 square units. Hence, a way to
defuzzify μ Q ∗ consists of searching the center of area, that is, a point y0 ∈ [0, 10]
such that
y0 y0
A
0.5dy = = 2.5, or dy = 5.
2
0 0
y
Hence, [y]00 = y0 = 5. The defuzzified value corresponding to x0 = 0.5, is y0 = 5.
The method, when the conditional is Mamdani, is graphically reflected as follows.
2nd Example. Identical to the first example, but with the input x0 = 0.3. It is
y y
μ Q ∗1 (y) = min(0.3, 1 − ), μ Q ∗2 (y) = min(1 − 0.3, ).
10 10
Hence,
⎧
⎨ 0.3, if 0 y 3
y
μ Q ∗ (y) = 10 , if 3 y 7
⎩
0.7, if 7 y 10
3 y0
y
0.3dy + dy = A/2 = 2.5.
10
0 3
y0 y0
Then 3 × 0.3 + 1
10 3 ydy = 2.5, or 3 ydy = 10(2.5 − 0.9) = 16. Thus,
y 2 y0
[ ] = 16 ⇒ y02 − 9 = 32 ⇒ y02 = 41 ⇒ y0 = 6.4.
2 3
The defuzzified value that corresponds to x0 = 0.3, is y0 = 6.4.
Defuzzifying with the centre of the area we obtained an output for all values x0
in [0, 1]. Let’s see it by means of the function CRI with defuzzification made by the
centre of the area.
1. The graphics, for any input x0 1/2, is
(1 − x0 − x0 )(10(1 − x0 ) − 10x0 )
A = 10(1 − x0 )x0 + (10 − 10(1 − x0 ))(1 − x0 ) +
2
3.5 Deffuzification 111
then,
and
A = 20x0 (1 − x0 ) + 5 − 20x0 + 20x02 = 5.
Hence,
y0
y
10x02 + dy = 2.5 ⇒ y0 = 50 − 100x02 ,
10
10x0
and
C R I (x) = 50 − 100x 2 , if x 1/2.
√
For example, C R I (0.5) = 0.5, and C R I (0.3) = 41 = 6.4, as it was shown. It
is also C R I (0.1) = 7.
√
Last formula, C R I (x) = 20 − 100x 2 , gives real values provided 20−100x 2
0, or x 2 1/5. Since, it is x 1/5, that implies x 2 1/4 1/2, and it follows
that the formula is useful for all x ∈ [0, 1] such that x 1/2.
2. For any input x0 1/2, the graphic is
y0 y0
y
10(1 − x0 ) +2
dy = 2.5, or 100(1 − x0 )2 + ydy = 2.5,
10
10(1−x0 ) 10(1−x0 )
112 3 Reasoning and Fuzzy Logic
giving y02 = 50 − 100(1 − x0 ), or y0 = 50 − 100(1 − x0 )2 , that gives real
values provided
√ 50 − 100(1 − x0 )2 0, equivalent to 1/2 (1 − x0 )2 , or to
x0 1 − 1/2. Since, 1/2 x, it is (1 − x)2 1/4 1/2. Then,
C R I (x) = 50 − 100(1 − x)2 , if 1/2 x.
For example,√ √ √
CRI(0.7) = 41 = 6.4, CRI(0.8) = 46 = 6.78, CRI(0.9) = 49 = 7.
3. Finally, with defuzzification by the centre of area, is:
50(1 − 2x ), if x 1/2
2
CRI(x) =
50(1 − 2(1 − x) ), if x 1/2.
2
√ √
Notice that√CRI(0) = 50, CRI(0.5) = 50(1 − 20.52 ) = 25 = 5, and
CRI(1) = 50.
The graphic of CRI is
As it was said before, the output is a logical consequence of the premises given by
the input and the rule. Notwithstanding, the situation is different if, taking the rule as
defining the system, only the input is considered as a premise. But, before to consider
this question, let us consider what happens when there is more than a single rule.
1. If μ Q ∗1 , μ Q ∗2 ∈ Cons({μ P ∗ }), from μ P ∗ max(μ Q ∗1 , μ Q ∗2 ) = μ Q ∗ , follows
μ Q ∗ ∈ Cons({μ P ∗ }).
2. If μ Q ∗1 or μ Q ∗2 is a conjecture of {μ P ∗ }, then μ Q ∗ ∈ Con j ({μ P ∗ }). The proof
follows in this way, provided it is, for instance, μ Q ∗1 ∈ Con j ({μ P ∗ }),
• μ Q ∗ = max(μ Q ∗1 , μ Q ∗2 ) ⇒ μ Q ∗1 μ Q ∗ , μ Q ∗2 μ Q ∗ ⇒ μQ ∗ μQ ∗ , μQ ∗
1
μQ ∗ ⇒ μQ ∗ min(μQ ∗ , μQ ∗ ).
2 1 2
• If μ P ∗ μQ ∗ , then μ P ∗ μQ ∗ , that is absurd. Hence, it is μ P ∗ μQ ∗ , and
1
μ Q ∗ ∈ Con j ({μ P ∗ }).
3.6 Rules and Conjectures 113
In conclusion
• If all the partial outputs μ Q ∗1 , μ Q ∗2 , ..., μ Q ∗n , are consequences of the input μ P ∗ ,
also the final output μ Q ∗ is a consequence of μ P ∗ .
• If at least one of the partial outputs μ Q i∗ (1 i n) is just a conjecture of the
input μ P ∗ , also the final output μ Q ∗ is a conjecture of μ P ∗ .
Nevertheless, it is not usual that μ Q ∗ results to be a consequence of the single input
μ P ∗ . Let’s us introduce a necessary and sufficient condition for it in the particular
case in which there is only one rule represented by J (a, b) = ab (Larsen).
Hence, in this particular case, the necessary and sufficient condition for being μ Q ∗ ∈
Cons({μ{x0 } }) is that μ P ∗ (x0 ) = μ Q ∗ (x0 ) = 1. Nevertheless, what happens in
most of the cases is that μ Q ∗ ∈ Con j ({μ{x0 } }), with μ Q ∗ ∈ Sp({μ{x0 } }), or μ Q ∗ ∈
H yp({μ{x0 } }).
Let’s show an example in which the output is a speculation of the input and other in
which the output is a hypothesis.
Example. Rule, “If x is small, then y is big”, with X = Y = [0, 10] and J (a, b) = ab
y
(Larsen), with μ S (x) = 1 − 10x
, μ B (y) = 10 , and x0 = 5. It is
5 y y
μ Q ∗ (y) = J (1 − , )= ,
10 10 20
and from the graphic
114 3 Reasoning and Fuzzy Logic
{5}
it is clear that μ Q ∗ is not comparable with μ{5} (μ Q ∗ NC μ{5} ), and that μ{5}
μQ ∗ = 1 − μ Q ∗ . Hence,
Example. Rule, “If x is big, then y is very big”, with X = Y = [0, 10] and the
observation “x is constantly equal to 0.8” for all x ∈ [0, 1]. Taking
follows
Graphically,
√
since y 2 = 0.8 means y = 0.8. Hence, it is μ0 = μ Q ∗ μ P ∗ , and μ P ∗ = μQ ∗ =
1 − μ Q ∗ , that imply μ Q ∗ ∈ Con j ({μ P ∗ }), and namely μ Q ∗ ∈ H yp({μ P ∗ }).
Notice that this second example contains the observation that the input is a con-
stant.
Last Remark
For any continuous t-norm T , the function μ : Y → [0, 1], defined by
does verify
3.7 Two Final Examples 115
For example, with the rule ‘If x is small, then y is big’ (X = Y = [0, 1]),
and J (a, b) = max(1 − a, b) that is a W -conditional, taking μ S (x) = 1 − x and
μ B (y) = y, follows:
• With T = W , μ(y) = Sup W ((1 − x), max(x, y)) = Sup (0, y − x) = y =
x∈[0,1] x∈[0,1]
μ B (y).
• With T = pr od, μ(y) = Sup (1 − x) max(x, y) = Sup max((1 − x)x, (1 −
x∈[0,1] x∈[0,1]
x)y) = y = max(1/4, y/2), not coincidental with μ B .
• with T = min, μ(y) = Sup min((1 − x), max(x, y)) = 1, or μ = μ1 , also not
x∈[0,1]
coincidental with μ B .
μ R (xi , y j ) = ri j , 1 i n, 1 j m, or,
μ R = r11 /(x1 , y1 ) + · · · + rnm /(xn , ym ),
In the finite case there is again another representation of a fuzzy relation by means
of a directed graph.
For example, if X 1 = {x1 , x2 } and X 2 = {y1 , y2 , y3 }, the fuzzy
0.5 0.7 1
relation [R] = , corresponds to the directed graph
0.8 0 0.8
4.2 How to Compose Fuzzy Relations? 119
0.5 y1
0.8
x1
0.7
y2
x2 1
0.8
y3
Given two fuzzy relations μ : X × Y → [0, 1], and σ : Y × Z → [0, 1], how can
we obtain a relation λ : X × Z → [0, 1] through μ and σ? To solve this problem,
there is the Sup − T product of fuzzy relations, given by
giving the matrix [λ] = (ti j ) as the Max-T product, or composition, of the matrices
[rik ] and [sk j ], that was introduced before by: (ti j ) = (rik ) ⊗T (sk j ).
120 4 Fuzzy Relations
• X = { p1 , . . . , p4 }, a set of patients
• Y = {s1 , s2 , s3 }, a set of symptoms
• Z = {d1 , . . . , d5 }, a set of deceases,
and the fuzzy relation σ
⎛ ⎞
0.7 0 0 0.3 0.6
[σ] = ⎝ 0.5 0.5 0.8 0.4 0 ⎠
0 0.7 0.2 0.9 0
showing the medical knowledge of how strongly each symptom is associated with a
decease. Suppose also that, by examining the patients, the doctors conclude the matrix
⎛ ⎞
0 0.3 0.4
⎜ 0.2 0.5 0.3 ⎟
[μ] = ⎜
⎝ 0.8 0
⎟
0 ⎠
0.7 0.7 0.9
that describes numerically how strongly the symptoms are manifested in the patients.
Then,
[λ] = [μ] ⊗min [σ]
t43 = max(min(0.7, 0), min(0.7, 0.8), min(0.9, 0.2)) = max(0, 0.7, 0.2) = 0.7.
The matrix [λ] results from a mixing between knowledge and observation.
1 From, Fuzzy Set Theory, by G.J. Klir, U.H. St. Clair, B. Yuan, Prentice/Hall, N.J., 1997.
4.2 How to Compose Fuzzy Relations? 121
with the matrices [σ]t , [μ]t , defined by σ t (x, y) = σ(y, x), μt (x, y) = μ(y, x),
giving
([μ]t )t = [μ].
The most relevant properties of a fuzzy relation μ : X × X → [0, 1] are the following,
1. Reflexive property, μ(x, x) = 1, for all x ∈ X .
2. Symmetric property, μ(x, y) = μ(y, x), for all x, y ∈ X , implies x = y.
3. Antisymmetric property, μ(x, y) > 0, μ(y, x) > 0 implies x = y.
4. T-transitive property T (μ(x, y), μ(y, z)) μ(x, z), for all x, y, z ∈ X , and
some continuous t-norm T .
In the finite case, for what concerns properties reflexive and symmetric, the matrix
[μ] = (ti j ) shows the respective properties,
1 . It is tii = 1, for all 1 i n, that is, the main diagonal of [μ] is constituted
by n numbers equal to 1.
2 . It is ti j = t ji , for all 1 i, j n, that is, the elements of [μ] are placed
symmetrically with respect to the main diagonal.
For example, the matrix
1 0.7
is reflexive, but not symmetric,
0.6 1
For what concerns the antisymmetric property, ti j > 0 and t ji > 0, implies i = j.
For example, the matrix ⎛ ⎞
1 0 0.7 0
⎜ 0.6 1 0 0.7 ⎟
⎜ ⎟
⎝ 0 0.5 1 0.8 ⎠
0.7 0 0 1
is antisymmetric.
or T (ti j , t jk ) tik , it is Max T (ti j , t jk ) tik . That is, [μ] ⊗T [μ] [μ].
1 j n
b. If [μ] ⊗T [μ] [μ], or Max T (ti j , t jk ) tik , follows T (ti j , t jk ) tik , for all
1 j n
1 i, j n. That is, μ is T -transitive.
implies
tik = Max T (ti j , t jk ), or [μ] = [μ] ⊗T [μ].
1 j n
Remark 4.3.1 The definitions given in this section contain the case of the corre-
sponding classical crisp definitions,
• If R ⊂ X × X is a classical reflexive relation in X , its membership function μ R
reflects (x, x) ∈ R for all x ∈ X , by μ R (x, x) = 1.
• If R ⊂ X × X is a classical symmetric relation in X,
• If R ⊂ X × X is antisymmetric,
for all t-norms T . Notice that if μ R (x, y) = 0 or μ R (y, z) = 0, then, for example,
Remark 4.3.3 If the fuzzy relation μ R is T-transitive, and T1 is a t-norm such that
T1 T , then μ R is also T1 -transitive, since
T1 (μ R (x, y), μ R (y, z)) T (μ R (x, y), μ R (y, z)) μ R (x, z),
for all x, y, z.
All constant fuzzy sets μk (x) = k for all x ∈ X , and k ∈ [0, 1], are T-states of
any fuzzy relation μ : X × X → [0, 1]: T (μk (x), μ(x, y)) μk (x) = k = μk (y).
For example, μ0 = μ∅ and μ1 = μ X , are always T -states. Hence, the set T (μ) of
all T -states σ of μ is never empty. From now on, in general we will only refer to non
constant T -states σ.
1
μx2 (x1 )μ(x1 , y) = μ(x1 , y),
8
showing,
• 8 μ(x 1 , x 1 )
1
= 8 = μx2 (x 1 )
1
• 8 μ(x 1 , x 2 )
1
= 1
82
< 18 = μx2 (x2 )
• 8 μ(x 1 , x 3 )
1
= 2
82
< 38 = μx2 (x3 ),
etc. That is, the three fuzzy sets μx1 , μx2 , μx3 are pr od-states of μ.
Remark 4.4.1 When the fuzzy relation μ represents a conditional statement Q/P
(a fuzzy rule ‘If x is P, then y is Q’), the T -states of μ are among the fuzzy sets
verifying the Modus Ponens with respect to the continuous t-norm T .
Given a fuzzy relation μ : X × X → [0, 1], the α-cuts of μ are the classical (crisp)
relations μ(α) defined by,
since,
Hence, only if μ is min-transitive, it is sure that all α-cuts are classical preorders,
and μ results decomposed in the family of preorders {μα ; α ∈ (0, 1)}.
is obviously reflexive but not symmetric, and verifies [μ] ⊗min [μ] = [μ]. Hence, μ
is a min-preorder. Its different α-cuts are
⎛ ⎞ ⎛ ⎞ ⎛ ⎞
1 0 1 0 1 0 1 0 1 1 1 1
⎜0 1 0 0⎟ ⎜0 1 0 0⎟ ⎜0 1 0 0⎟
[μ(1) ] = ⎜ ⎟ ⎜ ⎟ ⎜
⎝ 1 0 1 0 ⎠ , [μ(0.8) ] = ⎝ 1 0 1 0 ⎠ , [μ(0.6) ] = ⎝ 1 1 1 1 ⎠ ,
⎟
0 0 0 1 0 1 0 1 0 1 0 1
⎛ ⎞ ⎛ ⎞
1 1 1 1 1 1 1 1
⎜0 1 0 0⎟ ⎜1 1 1 1⎟
[μ(0.4) ] = ⎜ ⎟ ⎜
⎝ 1 1 1 1 ⎠ , [μ(0.3) ] = ⎝ 1 1 1 1 ⎠ ,
⎟
1 1 1 1 1 1 1 1
that give the classical (α) preorders, that follows:
• (1) : 1 (1) 1, 2 (1) 2, 3 (1) 3, 4 (1) 4, 1 (1) 3, 3 (1) 1.
• (0.8) : 1 (0.8) 1, 2 (0.8) 2, 3 (0.8) 3, 4 (0.8) 4, 1 (0.8) 3, 3 (0.8)
1, 4 (0.8) 2.
• (0.6) : 1 (0.6) 1, . . . , 4 (0.6) 4, 1 (0.6) 3, 3 (0.6) 1, 4 (0.6) 2, 1 (0.6)
2, 1 (0.6) 4, 3 (0.6) 2, 3 (0.6) 4.
• (0.4) : 1 (0.4) 1, . . . , 4 (0.4) 4, 1 (0.4) 3, . . . , 4 (0.4) 1, 4 (0.4) 3.
• (0.3) : 1 (0.3) 1, . . . , 4 (0.3) 3, 2 (0.3) 1, . . . , 4 (0.3) 2,
and, since, 0.3 0.4 0.6 0.8 1, verify (1) ⊂(0.8) ⊂(0.6) ⊂(0.4) ⊂(0.3) .
Example 4.5.2 The fuzzy relation μ : X × X → [0, 1], with X = {1, 2, 3, 4, 5, 6},
given by
⎛ ⎞
1 0.2 1 0.6 0.2 0.6
⎜ 0.2 1 0.2 0.2 0.8 0.2 ⎟
⎜ ⎟
⎜ 1 0.2 1 0.6 0.2 0.6 ⎟
[μ] = ⎜
⎜ 0.6 0.2 0.6 1 0.2 0.8 ⎟
⎟
⎜ ⎟
⎝ 0.2 0.8 0.2 0.2 1 0.8 ⎠
0.6 0.2 0.6 0.8 0.8 1
4.5 Fuzzy relations and α-cuts 127
is reflexive, symmetrical and min-transitive, since [μ] ⊗min [μ] = [μ]. Hence, all the
α-cuts are classical equivalence relations, each one defining a partition of X . The
different values of α are 0.2, 0.6, 0.8 and 1 (levels of crispness of μ), and it is easy
to see that the corresponding partitions π0.2 , π0.6 , π0.8 , and π1 , can be located as the
partition tree:
1 1 3 4 6 2 5
0.8 1 3 4 6 2 5
0.6 1 3 4 6 2 5
0.2 1 3 4 6 2 5
that only connects the pairs (1, 1), (2, 2), . . . (6, 6), and is called the disconnected
partial order.
128 4 Fuzzy Relations
5 6
⎛ ⎞
1 1 0 0 0 0
⎜0 1 0 0 0 0⎟
⎜ ⎟
⎜0 0 1 0 0 0⎟
2. [μ(0.8) ] = ⎜
⎜0
⎟ 2 3 4
⎜ 0 0 1 0 0⎟⎟
⎝0 0 0 0 1 0⎠
0 0 0 0 0 1
1
5 6
⎛ ⎞
1 1 0 1 1 0
⎜0 1 0 0 1 0⎟
⎜ ⎟
⎜0 0 1 0 0 0⎟
3. [μ(0.6) ] = ⎜
⎜0
⎟ 2 3 4
⎜ 0 0 1 1 0⎟⎟
⎝0 0 0 0 1 0⎠
0 0 0 0 0 1
1
5 6
⎛ ⎞
1 1 0 1 1 0
⎜0 1 0 0 1 0⎟
⎜ ⎟
⎜0 0 1 0 1 0⎟
4. [μ(0.5) ] = ⎜
⎜0
⎟ 2 3 4
⎜ 0 0 1 1 0⎟⎟
⎝0 0 0 0 1 0⎠
0 0 0 0 0 1
1
5 6
⎛ ⎞
1 1 0 1 1 1
⎜0 1 0 0 1 0⎟
⎜ ⎟
⎜0 0 1 0 1 0⎟
5. [μ(0.4) ] = ⎜
⎜0
⎟ 2 3 4
⎜ 0 0 1 1 0⎟⎟
⎝0 0 0 0 1 0⎠
0 0 0 0 0 1
1
4.5 Fuzzy relations and α-cuts 129
All this crisp partial orderings come from the fuzzy ordering given by [μ], and
which valued directed graph is
5 6
0.6
0.5 0.4
2 3 4
0.2
0.8 0.6
where some arrows are avoided because of the min-transitivity of μ. For instance,
the degree between 1 and 5, is
shortly written ET = min(JT , JTδ ), with JTδ (a, b) = JT (b, a). It is obvious that
as well as that ET (a, a) = 1, and ET (a, b) = ET (b, a). What it is not so obvious is
that relations JT are T-Transitive:
© Springer International Publishing Switzerland 2015 131
E. Trillas and L. Eciolaza, Fuzzy Logic, Studies in Fuzziness
and Soft Computing 320, DOI 10.1007/978-3-319-14203-6_5
132 5 T-Preorders and T-Indistinguishabilities
To avoid some difficult technicalities, we will exemplify this general result in the
particular case JW (a, b) = min(1, 1 − a + b):
T (ET (a, b), ET (b, c)) = T (min(JT (a, b), JT (b, a)), min(JT (b, c), JT (c, b)))
≤ T (JT (a, b), J(b, c)) = ET (a, c),
T (R(a, b), R(b, c)) = T (Inf Ri (a, b), Inf Ri (b, c))
i∈I i∈I
≤ Inf T (Ri (a, b), Ri (b, c))
i∈I
≤ Inf Ri (a, c) = R(a, c),
i∈I
and
E(a, b) = Inf Ei (a, b) = Inf Ei (b, a) = E(b, a).
i∈I i∈I
5.1 Which Is the Aim of This Section? 133
Finally,
T (E(a, b), E(b, c)) = T (Inf Ei (a, b), Inf Ei (b, c))
i∈I i∈I
≤ Inf T (Ei (a, b), Ei (b, c))
i∈I
≤ Inf Ei (a, c) = E(a, c)
i∈I
If σ ∈ [0, 1]X , define JTσ (x, y) = JT (σ(x), σ(y)). Obviously JTσ is a T-Preorder in
X × X, and if F ⊂ [0, 1]X it is also
Inf JTσ
σ∈F
a T-Preorder.
Given a fuzzy relation μ : X × X → [0, 1], consider the set T (μ) of its T-states,
that is the fuzzy sets σ : X → [0, 1], such that
thus,
σ ∈ T (μ) ⇔ μ ≤ JTσ
Hence,
σ ∈ T (μ) ⇔ μ ≤ Inf JTσ ,
σ∈T (μ)
For example,
• T = min, all min-preorders are of the form,
σ 1, if σ(x) ≤ σ(y)
μ(x, y) = Inf Jmin = Inf
σ∈T (μ) σ∈T (μ) σ(y), if σ(x) > σ(y)
Remark 5.2.1 As an immediate consequence of all that has been said, for any fam-
ily of fuzzy sets F ⊂ [0, 1]X , the T-preorder Inf σ∈F JTα has the elements in F as
T-states. If μ(x, y) = Inf σ∈F JTσ (x, y), it follows μ ≤ JTσ , ∀σ ∈ F, or, equivalently
σ ∈ T (μ), ∀σ ∈ F, or F ⊂ T (μ). For example, with X = [0, 1] and the two functions
σ1 (x) = x, σ2 (x) = 1 − x, it is
a T-preorder. With T = W ,
is a W-preorder.
ETσ (x, y) = min(JTσ (x, y), JTσ (y, x)) = min(JT (σ(x), σ(y)), JT (σ(y), σ(x))).
1, if σ(x) ≤ σ(y) 1, if σ(y) ≤ σ(x)
• T = min, μ(x, y) = Inf min ,
σ∈T (μ) σ(y), if σ(x) > σ(y) σ(x), if σ(y) > σ(x)
1, if σ(x) = σ(y)
= Inf
σ∈T (μ) min(σ(x), σ(y)), otherwise
1, if σ(x) = σ(y)
= Inf min(σ(x), σ(y)), otherwise
σ∈T (μ)
1, if σ(x) ≤ σ(y) 1, if σ(y) ≤ σ(x)
• T = prod, μ(x · y) = Inf min σ(y) , σ(x)
σ∈T (μ) σ(x) if σ(x) > σ(y)
, σ(y) , if σ(y) > σ(x)
= Inf
1,
if σ(x) = σ(y)
σ(y) σ(x)
σ∈T (μ) min σ(x) , σ(y) , otherwise
⎧
⎨ 1,
if σ(x) = σ(y)
= σ(y) σ(x)
⎩ Inf min σ(x) , σ(y) , otherwise
σ∈T (μ)
Remark 5.3.1 Like in the case of T-Preorders, for any family of functions F ⊂
[0, 1]X , the T-indistinguishability Inf σ∈F ETσ has the elements in F as T-states. Notice
that with μ(x, y) = Inf σ∈F ETσ (x, y) = Inf σ∈F min(JTσ (x, y), JTσ (y, x)), it follows
μ(x, y) ≤ JTσ (x, y), for all x, y in X, that is equivalent to σ ∈ T (μ). Hence,
∀σ ∈ F ⇒ σ ∈ T (μ), or F ⊂ T (μ).
a T-indistinguishability.
• With T = W , is
Since, f.e.,
σ1 1 2 σ1 3 1
JW (1, 2) = min 1, 1 − + = 1, JW (3, 1) = min 1, 1 − + = 1/2,
4 4 4 4
σ1 1 σ1 3
JW (4, 1) = min 1, 1 − 1 + = 1/4, JW (3, 4) = min 1, 1 − + 1 = 1, etc.
4 4
5.3 The Characterization of T-Indistinguishabilities 137
It is also ⎛ ⎞
1 3/4 1/2 1/4
σ2 ⎜ 1 1 3/4 1/2 ⎟
⎜
[JW ] = ⎝ ⎟
1/2 1 1 3/4 ⎠
1 1 1 1
since, f.e.,
σ2 3 1 σ2 1 3
JW (1, 2) = min 1, 1 − + = 3/4, JW (3, 1) = min 1, 1 − + = 1/2,
4 2 4 4
σ2 3 σ2 1
JW (4, 1) = min 1, 1 − 0 + = 1, JW (3, 4) = min 1, 1 − + 0 = 3/4, etc.
4 4
Hence,
⎛ ⎞
1 3/4 1/2 1/4
⎜ 3/4 1 3/4 1/2 ⎟
σ1 σ2 ⎜ ⎟
[EW ] = min([JW ], [JW ]) = ⎜ ⎟
⎝ 1/2 3/4 1 3/4 ⎠
1/4 1/2 3/4 1
1/4 3/4
1/2
1 1/2 3
3/4
3/4
2
138 5 T-Preorders and T-Indistinguishabilities
ET (μ, μ) = 1, or Sup T ◦ (μ × μ) = 1,
T (ET (μ, σ), ET (σ, λ)) = T (Sup T (μ, σ), Sup T (σ, λ))
= Sup T (T (μ, σ), T (σ, λ))
≤ Sup T (μ, λ) = ET (μ, λ),
it follows
μ · σ = 0.6|1 + T (0.4, 0.5)|2 + T (0.8, 0.8)|3 + 0.7|4,
and
• Emin (μ, σ) = 0.8, since min(0.8, 0.8) = 0.8, min(0.4, 0.5) = 0.4
• Eprod (μ, σ) = 0.64, since prod(0.8, 0.8) = 0.64, prod(0.4, 0.5) = 0.2
• EW (μ, σ) = 0.6, since W (0.8, 0.8) = 0.6, W (0.4, 0.5) = 0
Remark 5.3.4 Provided T = min or T = prodϕ , if ET (μ, σ) > 0, and ET (σ, α) > 0,
it is
0 < T (ET (μ, σ), ET (σ, α)) ≤ ET (μ, α),
what follows is that, from the equivalence given by the existence of r > 0 such that
E(x, y) > r, E(y, z) > r, it is
Wϕ (r, r) = ϕ−1 (max(0, 2ϕ(r) − 1)) ≤ Wϕ (E(x, y), E(y, z)) ≤ E(x, z).
Hence, to have 0 < E(x, z), it is necessary that 0 < Wϕ (r, r), that is,
6.1 Introduction
This extension includes the crisp subsets A ⊂ R, since μA ∈ {0, 1}R ⊂ [0, 1]R .
For example, with A1 = {1, 2, 3}, and A2 = {1, 3, 5}, and only taking into account
the numbers in N ⊂ R, it is
Analogously, it results
• [a, b] [c, d] = [a − d, b − c]
• [a, b] ⊗ [c, d] = [min(ad, ac,bd, bc),
max(ad,ac, bd,bc)]
• If 0 ∈ ÷ [c, d] = min c , d , c , d , max c , d , c , d
/ [c, d], [a, b] a a b b a a b b
.
For example,
• [7, 8] ⊕ [−1, 9] = [6, 17]
• [7, 8] [−1, 9] = [−2, 9]
• [3, 4]
[2, 2] = [6, 8]
• ÷ [1, 2] = [2, 10]
[4, 10]
• 2 ⊕ [7, 8] = [2, 2] ⊕ [7, 8] = [9, 10]
• 2
[7, 8] = [2, 2]
[7, 8] = [min(14, 16, 14, 16), max(14, 16)] = [14, 16]
7 8 7 8 7
÷ 2 = [7, 8]
• [7, 8] ÷ [2, 2] = min 2 , 2 , max 2 , 2 = 2, 4
In short, and accordingly with the ‘preservation of the classical case’, through
the extension principle both the numerical and the interval arithmetics are preserved.
What follows is a yet larger arithmetic with fuzzy concepts.
For all μ ∈ [0, 1]R , and all r ∈ [0, 1], it can be computed that:
• (r ⊕ μ)(t) = (μr ⊕ μ)(t) = Sup min(μr (x), μ(y))
t=x+y
1, μ(y), if x = r μ(y), x = r
= Sup min = Sup
t=x+y 0, μ(y), if x = r t=x+y 0, x = r
= μ(t − r).
Hence ((−1) ⊕ μ)(t) = μ(t − 1), (μ1 ⊕ μ)(t) = μ(t − 1), (μ 0 ⊕ μ)(t) = μ(t).
• (r
μ)(t) = (μr
μ)(t) = Sup min(μr (x), μ(t)) = μ rt , if r = 0. Hence,
t=x·y
1
r
μ (t) = μ(rt), (1
μ)(t) = μ(t), (μ1
μ)(t) = μ(t).
• (μ0
μ)(t) = Sup min(μ0 (x), μ(t)) = 0 = μ0 (t).
t=x·y
1
μ t , if t = 0
• μ1 (t) = Sup min(1, μ(y)) = Supμ(y) =
t= x t= x 0, if t = 0
y y
compute 2
μ and 2 ⊕ μ.
It is (2
μ)(t) = μ(t/2), and (2 ⊕ μ)(t) = μ(t − 2). Hence,
6.1 Introduction 143
It should be pointed out that, since the cartesian product μ × σ of μ ∈ [0, 1]X and
σ ∈ [0, 1]Y , is defined by
the operations , with ∗ ∈ {+, −, ·, :}, can also be represented in the form
Example 6.1.2 If X = {1, 2, 3}, and μ = 0.7|1 + 0.9|2 + 1|3, σ = 0.8|1 + 0.9|3,
since X × X = {(1, 1), (1, 2), . . . , (3, 2), (3, 3)}, it results
μ ⊕ σ(2) = 0.7, μ ⊕ σ(3) = max(0, 0.8) = 0.8, μ ⊕ σ(4) = max(0.7, 0.8, 0) = 0.8,
Hence,
μ ⊕ σ = 0.7|2 + 0.8|3 + 0.8|4 + 0.9|5 + 0.9|6.
In the usual scientific computation √is rather unusual to consider exact numbers, as in
the case of taking π ≈ 3.1416, or 2 ≈√ 1.4142. In some cases, an interval containing
the number is considered, for example, 2 ∈ [1.412, 1.413]. Fuzzy numbers are just
a “fuzzification” of this last idea. For example, the fuzzy number “approximately 2”
can be taken either as the fuzzy set,
1
0 1 1.9 2 2.1
That kind of fuzzy sets are known as fuzzy numbers. Their definition is the fol-
lowing.
Definition 6.2.1 A fuzzy set μ ∈ [0, 1]R is a fuzzy number, provided there is a
closed interval [a, b] ∈ R, such that
1. μ(x) = 1, for all x ∈ [a, b]
2. It exists a function L : (−∞, a) → [0, 1], that is continuous and non-decreasing,
such that μ(x) = L(x), for all x ∈ (−∞, a).
3. It exists a function R : (b, +∞) → [0, 1], that is continuous and decreasing, such
that μ(x) = R(x), for all x ∈ (b, +∞).
6.2 Fuzzy Numbers 145
It can be proven that, if ∗ ∈ {+, −, ×, :}, and if μ, σ are fuzzy numbers, then μ σ
is also a fuzzy number. For example, with the fuzzy number,
⎧
⎨ x − 2, if x ∈ (2, 3)
μ3 (x) = 4 − x, if x ∈ (3, 4)
⎩
0, otherwise,
since for t < 4, or x+y < 4, and for t > 8, or x+y > 8, it should be (μ3 ⊕μ3 )(t) = 0;
for t = 6, or x + y = 6, it should be (μ3 ⊕ μ3 )(6) = 1, and:
• For t ∈ [4, 6], L is the segment joining (4, 0) and (6, 1), that is
x y 1
x−4
0 = 4 0 1 = −x + 2y + 4, or y =
6 1 1 2
146 6 Fuzzy Arithmetic
• For t ∈ [6, 8], R is the segment joining (6, 1) and (8, 0), that is
x y 1
8−x
0 = 6 1 1 = x + 2y − 8, or y =
8 0 1 2
Graphically
Notice that the result is also a fuzzy number μ6 , such that the amplitude of it
([8, 4], 8 − 4 = 4) is twice than that of the initial μ3 ([2, 4], 4 − 2 = 2). When
operating with fuzzy numbers, the amplitude grows. If there is uncertainty about the
value 3, the uncertainty about 3 + 3 is twice, and such uncertainty will grow each
time we ‘add’ more numbers.
Let’s show a systematic way of operating with triangular fuzzy numbers, that is,
those represented by ⎧
⎨ L(x), if x ∈ (a, n)
μn (x) = R(x), if x ∈ (n, b)
⎩
0, otherwise,
Namely, take
and
t = x + y = (α + β)z + n + m − α − β,
giving
t − (n + m) + (α + β)
z= ,
α+β
that is
t − (n + m) + (α + β)
(μn ⊕ μm )(t) = , if t ∈ (m + n − (α + β), m + n).
α+β
t − (n + m) − (α + β)
(μn ⊕ μm )(t) = , if t ∈ (n + m, m + n − (α + β)).
α+β
it results
⎧ t+2 ⎫
⎨ 2 ,
⎪ if x ∈ (−2, 0) ⎪
⎬
(μ3 μ3 )(t) = 2−t
2 , if x ∈ (0, 2)
⎪
⎩ ⎪
⎭
0, otherwise.
= μ0 (t), with the interval [a, b] = [−2, +2].
Finally, ⎧√ ⎫
⎨ t −√2, if t ∈ [4, 9] ⎬
(μ3
μ3 )(t) = 4 − t, if t ∈ [9, 16] = μ9 (t).
⎩ ⎭
0, otherwise.
Graphically
Remark 6.2.2 It should be pointed out that, in the case of the product of triangular
fuzzy numbers, the result is not a triangular (linear) fuzzy number, since functions L
and R are not linear. In addition the interval [a, b] is not symmetrical. Look that for
μ9 is [4, 16] whose mid point is not 9 but 10.
μm 1x , if x=0
or μn ÷ μm = μn
μ , with μ
1 1
(x) = , and
m m 0, if x = 0
m = 1.
1 1
μm
÷ μ3 , with μ3 as in Sect. 6.2.2. It will be:
Let’s compute the example μ3
⎧ 1
⎪ 4 − 1y , if y ∈ 4, 3
1
⎨ 1 1
1
μ3 (y) = 1y − 2, if y ∈ 3, 2
⎪
⎩
0, otherwise.
Hence, for μ3
÷ μ3 = μ3
1
μ3 , it will be (μ3
÷ μ3 )(1) = 1, and
α+2
• L, α = 4 − 1
y =x−2⇒t =x·y = 4−α ⇒α= 4t−2
t+1
• R, α = 4 − x = 1
y −2⇒t =x·y = 4−α
α+2 ⇒α= 4−2t
t+1 .
That is ⎧ 4t−2 ⎫
⎨ t+1 ,
⎪ if t ∈ [ 21 , 1] ⎪
⎬
(μ3
÷ μ3 )(t) =
t+1 ,
4−2t
if t ∈ [1, 2] = μ1 (t)
⎪
⎩ ⎪
⎭
0, otherwise
Graphically,
Remark 6.2.3 Obviously, like in the product’s case, the result is not a linear-
triangular fuzzy number, and the basic interval [a, b] is not symmetrical respect
to the point mn .
6.2.3 Note
As it is easy to see, if μ, σ ∈ [0, 1]R , and ∗ ∈ {+, −, ·, :}, it follows the relation
between the corresponding α-cuts:
6.2 Fuzzy Numbers 151
(μ σ)α = μα σα ,
It is for this equality that the before hand computations were made. For example,
with
⎧ x+1 ⎧ x−1
⎨ 2 ,
⎪ if t ∈ (−1, 1) ⎨ 2 ,
⎪ if t ∈ (1, 3]
(μ1 )(x) = 3−x
2 , if t ∈ (1, 3) , (μ3 )(x) = 5−x
2 , if t ∈ (3, 5)
⎪
⎩ ⎪
⎩
0, otherwise 0, otherwise
with which
[μ1 ⊕ μ3 ]α = [4α, 8 − 4α], for α ∈ (0, 1].
Hence,
⎧x
⎨ 4,
⎪ if t ∈ (0, 4]
(μ1 ⊕ μ3 )(t) = 4 ,
8−x
if t ∈ (4, 8] = μ4 (t)
⎪
⎩
0, otherwise
As it is well known, (R, min, max) is a distributive lattice that come from the totally
ordered set (R, ≤). The order ≤ is definable from the lattice operations min, max by
(μ
∧ σ)(t) = Sup (μ × σ)(x, y), and (μ
∨ σ)(t) = Sup (μ × σ)(x, y).
t=min(x,y) t=max(x,y)
μ ≤∗ σ ⇔ μ
∧σ = μ ⇔ μ
∨ σ = σ.
6.3.1 Example
With X = {1, 2, 3}, and μ = 0.8|1 + 0.7|2 + 1|3, σ = 0.9|1 + 1|2 + 0.6|3, compute
∧ σ, μ
μ ∨ σ.
Since, t = min(x, y) and t = max(x, y) belong to {1, 2, 3}, it results:
∧ σ)(t) =
• (μ Max [min(μ(x), σ(y))], and:
t=min(x,y)
– (μ
∧ σ)(1) = max(min(μ(1), σ(1)), min(μ(1), σ(2)), min(μ(2), σ(1)),
min(μ(1), σ(3)), min(μ(3), σ(1)) = max(min(0.8, 0.9), min(0.8, 1),
min(0.7, 0.9), min(0.8, 0.6), min(1, 0.9)) = max(0.8, 0.8, 0.7, 0.6, 0.9) = 0.9
∧ σ)(2) = max(min(μ(2), σ(3)), min(μ(3), σ(2)))
– (μ
= max(min(0.7, 0.6), min(1, 1)) = max(0.6, 1) = 1
∧ σ)(3) = min(μ(3), σ(3)) = min(1, 0.6) = 0.6
– (μ
that is μ
∧ σ = 0.9|1 + 1|2 + 0.6|3.
This fuzzy set is different from μ · σ = 0.8|1 + 0.7|2 + 0.6|3, with the t-norm min.
∨ σ)(t) = Max [min(μ(x), σ(y))], and:
• (μ
t=max(x,y)
– (μ
∨ σ)(1) = min(μ(1), σ(1)) = min(0.8, 0.9) = 0.8
∨ σ)(2) = max(min(μ(1), σ(2)), min(μ(2), σ(1)))
– (μ
= max(min(0.8, 1), min(0.7, 0.9)) = max(0.8, 0.7) = 0.8
6.3 A Note on the Lattice of Fuzzy Numbers 153
– (μ
∨ σ)(3) = max(min(μ(1), σ(3)), min(μ(3), σ(1)), min(μ(2), σ(3)),
min(μ(3), σ(2)), min(μ(3), σ(3))) = max(min(0.8, 0.6), min(1, 0.9),
min(0.7, 0.6), min(1, 1), min(1, 0.6)) = 1,
that is μ
∨ σ = 0.8|1 + 0.8|2 + 1|3.
This fuzzy set is different from μ + σ = 0.9|1 + 1|2 + 1|3, with + the t-conorm
max.
(μ
∧ σ)(t) = max[Sup min(μ(x), σ(t)), Sup min(μ(t), σ(y))],
t≤x t≤y
it results:
1. If t ≤ a1 , (μ
∧ σ)(t) = 0
2. If a2 ≤ t, (μ ∧ σ)(t) = 0
3. If a1 ≤ t ≤ b1 , (μ
∧ σ)(t) = μ(t)
4. If b1 ≤ t ≤ c, (μ ∧ σ)(t) = σ(t)
5. If c ≤ t ≤ a2 , (μ
∧ σ)(t) = μ(t).
∧ σ = μ, and μ ≤∗ σ, although it is not μ ≤ σ.
Hence, μ
154 6 Fuzzy Arithmetic
In classical logic, only two quantifiers are considered. The universal quantifier ∀ (for
all), and ∃ (it exists, or for some), the existential quantifier, with the addition of ∃!
(exists only one). For example, given a sequence of real numbers (an ), it is said that
the real number a is its limit, when
A question which answer is, obviously, 35 × 25 = 875, that is, at least the 8.75 %
of the employees are young.
Another example is
Between 15 and 25 employees are married
Between 5 and 10 married employees are young
What can be said on the employees that are young?
in which what matters is the length of the two intervals [15, 25] (l = 10), and
[5, 10] (l = 5), that give at least, the interval [5, 10 + (10 − 5)] = [5, 15] of
employees that are young.
In natural language imprecise quantifiers like ‘about five’, ‘about half’, ‘most’,
etc., appear and can be represented by means of fuzzy numbers. For example,
6.4 A Note on Fuzzy Quantifiers 155
and the scores of LF of the (supposedly) five students {1, 2, 3, 4, 5} in the class,
given by
μLF = 0|1 + 0|2 + 0.75|3 + 1|4 + 0.5|5,
since the equation of the line joining the points (2, 0) and (3, 1), in the figure of μA3 ,
is y = x − 2.
There are about 6 employees in the company that are young and whose
computer skills are high.
with H1 (x) = μP1 (H1 (x)), H2 (x) = μP2 (H2 (x)), for all x ∈ X, and that correspond
with the rewritting:
There are about 6 employees that are young and with high computer skills,
n
of the given statement. Finally, with Z = |H1 ∩ H2 | = i=1 min(μP1 (F1 (xi )),
μP2 (F2 (xi ))), and Q(Z) = Q(|H1 ∩ H2 |), results the more compressed form
Z is Q
6.4 A Note on Fuzzy Quantifiers 157
corresponding to:
The number of employees that are young and with high computer skills, is
about 6,
Z is Q,
|H|
with Z = |X| , that gives the truth value t = Q(Z).
Z1 is Q1
Z2 is Q2
A rule often used to obtain Q is the following. If its exists f : Rn → R such that
Z = f (Z1 , . . . , Zn ), take Q = f (Q1 , . . . , Qn ). In the example, it is Z1 ·Z2 = |X2 | = Z,
hence Q = Q1 · Q2 , and the inference is
As a final example,
Most of the workers are young
About half of the young workers are women
Most × About half workers are young women.
Chapter 7
Fuzzy Measures
7.1 Introduction
Fuzzy sets not only appear by representing imprecise predicates, but also partial or
incomplete information. This is the case, for example, of a function X taking values
in [0, 10], of which it is only known that
5 X 7, 8 X , and X 2.
measures, as it is, for example, the number of apples in a basket. But, how can the
concept of measure be formalized?
Example 7.2.1 Take F = [0, 1] X , (L , ) = ([0, 1], ) the unit interval with the
partial linear order of the real line, and the qualitative relation “μ is less fuzzy than
σ”, translated by the so-called sharpened order S (=) defined by
μ(x) σ(y), if σ(x) 1/2
μ S σ ⇔
σ(x) μ(x), if σ(x) > 1/2,
that is a reflexive, transitive and antisymmetric, crisp relation. The fuzzy set μ0.5 is
the highest one, and all crisp sets μ ∈ {0, 1} X are minimals in (F, S ). Any mapping,
m : [0, 1] X → [0, 1] such that
• If μ is crisp, then m(μ) = 0
• m(μ0.5 ) = 1
• If μ S σ, then m(μ) m(σ),
is a ([0, 1], )-measure since m(μ0 ) = 0 because of μ0 ∈ {0, 1} X . These measures
are called measures of fuzziness, or fuzzy entropies.
If X = {x1 , . . . , xn } is finite, the following mappings are examples of fuzzy
entropies:
1. m(μ) = 1 − 2 max |μ(xi ) − 21 |
1i n
2. m(μ) = 2 max μ(xi ) · (1 − μ(xi ))
1i n
n
3. m(μ) = i=1 σ(μ(xi )), with σ(x) = x ln x − (1 − x) ln(1 − x) (logarithmic
entropy).
1 n 1, if μ(x) > 0.5
4. m(μ) = 2n 1=1 |μ(xi ) − μCμ (xi )|, with μCμ (x) = , the
0, if μ(x) 0.5
closest crisp set to μ (linear index of fuzziness)
n
5. m(μ) = 2n 1
i=1 (μ(x i ) − μCμ (x i )) (quadratic index of fuzziness).
2
7.2 The Concept of a Measure 161
Remark 7.2.2 1. For some specific problems, measures of fuzziness are selected
verifying the additional property of symmetry:
• For some negation N , it is m(μ) = m(N ◦ μ) = m(μ ).
For example, with N = 1 − id, measures 1, 2, 3, 4 do verify this property of
symmetry m(μ) = m(1 − μ).
Example 7.2.3 Take F ∈ [0, 1] X , with the partial pointwise order ‘μ σ ⇔ μ(x)
σ(x), for all x ∈ X ’ (and such that μ0 , μ1 ∈ F). A mapping m : [0, 1] X → [0, 1] is
a fuzzy measure provided m verifies:
1. m(μ0 ) = 0
2. m(μ1 ) = 1
3. If μ σ, then m(μ) m(σ).
μ · σ μ, μ · σ σ, μ μ + σ, σ μ + σ.
m(μ · σ) m(μ), m(μ · σ) m(σ), m(μ) m(μ + σ), m(σ) m(μ + σ),
and
m(μ · σ) min(m(μ), m(σ)) max(m(μ), m(σ)) m(μ + σ).
for all A, B ∈ P(X ). Notice that fuzzy measures can be applied to both [0, 1] X and
P(X ).
162 7 Fuzzy Measures
Given a triplet (X, F, m), where m is a fuzzy measure, if for some negation (N ) and
some union +(S), is
1. When μ σ , then m(μ + σ) m(μ) + m(σ), m is sub-additive
2. When μ σ , then m(μ + σ) m(μ) + m(σ), m is super-additive,
and when m is both sub-additive and super-additive, that is
This classification (once completed with those measures that are neither sub-
additive, nor super-additive), in the case in which μ, σ ∈ {0, 1} X , particularizes
to:
• If A ∩ B = ∅, and m(A ∪ B) m(A) + m(B), m is sub-additive
• If A ∩ B = ∅, and m(A ∪ B) m(A) + m(B), m is super-additive
• If A ∩ B = ∅, and m(A ∪ B) = m(A) + m(B), m is additive.
|A|
Example 7.3.1 The measure m(A) = n , in a finite set X = {x1 , . . . , xn }, is addi-
tive.
7.4 λ-Measures
A ∩ B = ∅ : m λ (A ∪ B) m λ (A) + m λ (B),
that is, if λ ∈ (−1, 0), all the corresponding λ-measures are sub-additive. As it
is easy to prove, if λ ∈ (0, +∞), m λ is super-additive.
• As it is well known, if X is a finite set X = {x1 , . . . , xn }, all probabilities m 0 :
P(X ) → [0, 1] are defined by choosing n numbers m 0 ({xi }) ∈ [0, 1], 1 i n,
n
verifying i=1 m 0 ({xi }) = 1, because 1 = m 0 ({x1 , . . . , xn }) = m 0 (x1 ) + . . . +
m 0 (xn ). Something analogous happens with λ-measures when X = {x1 , . . . , xn }.
For example, if X = {x1 , x2 , x3 }, it follows
1 = m λ (X ) = m λ ({x1 , x2 , x3 })
= m λ ({x1 , x2 } ∪ {x3 }) = m λ ({x1 , x2 }) + m λ (x3 ) + λm λ ({x1 , x2 })m λ (x3 )
3
= m λ (xi ) + λ m λ (xi )m λ (x j ) + λ2 m λ (x1 )m λ (x2 )m λ (x3 ).
i=1 1=i< j=3
1 = m λ (x2 )+m λ (x3 )+λ[m λ (x2 )m λ (x3 )] = m λ (x2 )+m λ (x3 )+λm λ (x2 )m λ (x3 )
1−m λ (x2 )
that is: m λ (x3 ) = 1+λm λ (x2 ) . With m λ (x 2 ) = 0.7, results m λ (x 3 ) = 1+0.3λ . With
0.3
forcing that some of the values π(xi ) should equal 1. For example, if X = {x1 , x2 , x3 },
the three values π(x1 ) = 0, π(x2 ) = 0.5, π(x3 ) = 1, define a particular mea-
sure of possibility on P(X ). It is, for example, π({x1 , x2 }) = max(0, 0.5) = 0.5,
π({x1 , x3 }) = max(0, 1) = 1, etc.
Theorem 7.5.2 For each μ ∈ [0, 1] X such that Sup μ = 1, the mapping πμ : F →
[0, 1] defined by
πμ (A) = Sup min(μ(x), μ A (x)), A ∈ F,
x∈X
is a possibility measure.
Hence, all non self-contradictory fuzzy sets, those μ such that Sup μ = 1, can be
viewed as possibility distributions.
Theorem 7.5.4 Given a possibility measure π : F → [0, 1], the function Nπ (A) =
1 − π(AC ), for all A ∈ F, is a necessity measure.
Theorem 7.5.5 Given a necessity measure N : F → [0, 1], the function π N (A) =
1 − N (AC ), for all A ∈ F is a possibility measure.
Proof
1 = max(π(A), π(AC )), and 0 = min(N (A), N (AC )),
Remark 7.5.8 Although the proof will not be presented, let’s show the following
important notice. In the case X is finite, for any possibility measure π it exists a (non
unique!) fuzzy set μ ∈ [0, 1] X with Sup μ = 1 such that π = πμ . In the finite case,
all possibility measures come from possibility distributions.
7.6 Examples
Example 7.6.1 On the age of a person p it is only available the incomplete informa-
tion given by
1. 37 Age( p) 41
2. It is neither Age( p) 32, nor Age( p) 43.
What can be said about the possibility and the necessity of “Age( p) 42”,
“Age( p) 40”, and “Age( p) 33”?
Solution
The available incomplete information can be represented by the following possi-
bility distribution μ:
Hence,
• πμ (Age( p) 42) = πμ ([42, 100]) = Sup min(μ(x), μ[42,100] (x)) =
x∈[0,100]
Sup μ(x) = μ(42) : ∈ (0, 1). Hence Nπμ (Age( p) 42) = 0.
x∈[42,100]
The value μ(42) can be computed as follows. The segment between (41, 1) and
(43, 0), verifies
168 7 Fuzzy Measures
x y 1
43 − x
0 = 41 1 1 = x + 2y − 43 ⇒ y = ,
43 0 1 2
Example 7.6.2 It is only known that a function F : X → [0, 1] doesn’t take any
value below 1/4 but takes values above 3/4. What can be said on the possibility and
necessity of the imprecise statements “F is small”, and “F is not small”?
Solution
The available but incomplete information can be translated by the possibility
distribution μ
7.6 Examples 169
Remark 7.6.3 Notice that this example is, like the following, not with questions
related to precise or crisp sets, but to impreciseness (fuzzy sets). Although Possibility
Theory is introduced with crisp sets, it is also applicable to fuzzy sets within the theory
given by the triplet (min, max, 1 − id).
Find the possibility and the necessity of the statement “John is around 35 years
old”.
Solution. The graphics of μY and μ P , with P = Ar ound 35, are
Hence, πμ (John is around 35 years old) = Sup x∈[0,100] min(μY (x), μ P (x)) =
2
3 , since it does correspond with the intersection of μY and μ P , that is, of the straight
10 , y = 5 . It results x = 3 and y = 3 . Since
lines respectively given by y = 40−x x−3 100 2
πμ (John is around 35 years old) < 1, it results Nπμ (John is around 35 yearsold)
= 0.
170 7 Fuzzy Measures
provided the pair (N , π) is dual. A probability p is said consistent with the dual pair
(N , π) if
In this hypothesis, if N (A) > 0, it is also p(A) > 0, and π(A) > 0, and if
π(A) = 0 it is N (A) = p(A) = 0. That is,
• If something is just a little bit probable, it is possible, but not necessarily necessary.
Example 7.7.1 Let is X = {1, 2, 3}, and μ = 0.7|1 + 1|2 + 0.5|3. Then, with
(A)
πμ (A) = Sup min(μ(i), μ A (i)), and Nμ = 1 − πμ (A), we get:
i∈X
1. πμ (1) = μ(1) = 0.7, πμ (2) = μ(2) = 1, πμ (3) = μ(3) = 0.5.
πμ ({1, 2}) = max(μ(1), μ(2)) = 1, πμ ({1, 3}) = max(0.7, 0.5) = 0.7,
πμ ({2, 3}) = 1, πμ {X } = 1.
2. Nμ (1) = 1 − πμ ({2, 3}) = 1 − 1 = 0, Nμ (2) = 1 − πμ ({1, 3}) = 1 − 0.7 = 0.3,
Nμ (3) = 1 − πμ ({1, 2}) = 0
Nμ ({1, 2}) = 1 − πμ (3) = 0.5, Nμ ({1, 3}) = 1 − πμ (2) = 0, Nμ ({2, 3}) =
1 − πμ (1) = 0.3, Nμ (X ) = 1.
Hence, the consistent probabilities are given by triplets p(1), p(2), p(3) in [0, 1]
such that p(1) + p(2) + p(3) = 1, and verifying:
For instance, with p(1) = 0.4, p(2) = 0.5, p(3) = 0.1, we have a consistent
probability, as well as with p(1) = 0.5, p(2) = 0.3, p(3) = 0.2. But the probability
given by p(1) = 0.6, p(2) = 0.2, p(3) = 0.2 is not consistent, because of p(2) <
0.3. With it there is an element whose probability is smaller than its necessity. In the
same vein, the probability given by the triplet p(1) = 0.1, p(2) = 0.3, p(3) = 0.6, is
also non-consistent because one of the probabilities is greater than the corresponding
possibility.
• If A, B ∈ F ⇒ A ∩ B, A ∪ B, AC , B ⊂ ∈ F.
• ∅ ∈ F, X ∈ F.
Proof Items (1) and (2) just follow from the fact that p is a 0-measure. p(X ) =
p(∅⊂ ) = 1 − p(∅) = 1. Finally, since
A ∪ B = (A ∩ B) ∪ (B − A) ∪ (A − B),
with (A ∩ B) ∩ (B − A) = ∅, (A ∩ B) ∩ (A − B) = ∅, and (B − A) ∩ (A − B) = ∅,
follows
p(A ∪ B) = p(A ∩ B) + p(B − A) + p(A − B).
Of special importance for the applications are the probabilities defined in the real
line R, with F = B the so-called Borel’s algebra, given by all the unions, comple-
ments and intersections of the open, closed, semi-open, and semi-closed intervals of
R. Then, if A ∈ B, the probability of A is defined by the Lebesgue-Stieltjes integral
p(A) = dP = μ A (x)d x = E(μ A ),
A R
Example 7.8.2 Which is the probability of the fuzzy set μ (fuzzy event) given by
Solution
10 · p(μ) = [0,10] μd x = [3,7] μd x = [3,4] μ1 d x + [4,6] 1d x + [6,7] μ2 d x =
2 + [3,4] (x − 3)d x + [6,7] (7 − x)d x = 2 + 21 + 21 = 3. Then p(μ) = 10
3
= 0.3.
7.8 Probability of Fuzzy Sets 173
Solution
5
6
10· p(μ) = [0,10] μd x = [4,5] μ1 d x + [5,6] μ2 d x = 4 (4−x)d x + 5 (6−x)d x =
− 21 + 23 = 1. Then p(μ) = 0.1.
σ μ
0 1/4 3/4 1
1/2
1 7/16
Then p(σ · μ) = 1/4 (4 − x)d x + 1/2 xd x = 7
16 , and p(σ|μ) = 1/2 = 78 .
Chapter 8
An Introduction to Fuzzy Control
8.1 Introduction
Probably one of the most successful developments of fuzzy reasoning, from the
industrial point of view, is the design of fuzzy control systems, also called linguistic
control systems, or simpler, the applications of fuzzy controllers.
A fuzzy control system is based on a set of fuzzy “if-then” rules of behavior
that consider the kind of stimuli from the environment, that the system will receive,
meanwhile at a given time, the values of these stimuli represent the facts, that the
rules have to consider to offer proper actions.
As it has been proved from its origins, fuzzy control should be useful in situations
where (a) There is no acceptable mathematical model for the plant, (b) There are
experienced human operators who can satisfactorily control the plant and provide
qualitative control rules in terms of vague and fuzzy sentences, and (c) In applica-
tions where there is a large uncertainty or unknown variation in plant parameters and
structures. In this cases, fuzzy control can be considered as a model-free approach
and it does not require a mathematical model of the objective plant. It is referred to
as a knowledge-based control approach, and it makes an effective use of all available
information related to the system, from sensors which provide numerical measure-
ments of key variables to human experts who provide linguistic descriptions about
the system and control instructions.
However, currently fuzzy control approach is mainly devoted to model-based
methods. On the one hand, there are cases in which operators cannot precisely tell
their action in a particular situation, or simply, operator’s control may not be always
optimal with respect to some performance objective. In this sense, identification for
obtaining fuzzy models from process data is very important.
On the other hand, although the uncertainty in dynamic systems can be handled
using appropriate control actions (e.g., their sensitivity to external disturbances and
parameter changes can be reduced stabilizing an unstable system, or the systems
dynamic behavior could be modified speeding up a slow system), the application
of control action over a system can potentially destabilize stable plants. Thus,
Data Rule
Base Base
IN OUT
Inference
Fuzzification Defuzzification
Engine
The data base contains all information needed to specify a particular configuration
of the fuzzy controller. It will have information related to the number, shape and dis-
tribution of the fuzzy sets specifying the meaning of the linguistic terms of each
linguistic variable involved in the application. The data base also contains infor-
mation related to the kind of crisp-fuzzy and fuzzy-crisp conversions, and the type
of operations for numerical calculations of the conjunction or disjunction premises.
Similarly with respect to the availability of operations for the implication, that real-
izes the “if-then” connective of the rules, as well as for the operation that computes
the aggregation of the conclusions of several rules that might be simultaneously
activated by the prevailing conditions of the environment.
The rule base contains the set of rules that will govern the behavior of the controller.
Through the rule base the task of modeling and control design is decomposed into a
group of local tasks, which tend to be easy to handle.
8.1 Introduction 177
The processing unit cares for the compatibility of data (input and output interfaces)
and for the execution of the rules (inference engine), which is mostly done through
pointwise numerical calculations of implications.
As it has been mentioned in previous chapters, the basic rules of reasoning used
in classical logic are the modus ponens for forward reasoning and modus tollens for
backwards reasoning. The symbolic expression of the modus ponens,
A B
A
B,
means that if the rule “if A then B” is given and the event A is observed, then the
event B should also be observed.
These processes are referred to in the literature as inference. As it has been
presented in this book, in the case of fuzzy logic, a generalization of modus ponens
is used, based on fuzzy sets. Given a universal set X , a fuzzy set A on X is defined by
its membership function μ A : X → [0, 1] and for all x ∈ X , μ A (x) gives the degree
of membership of x to A, or the degree with which x fulfills the concept represented
by A. Fuzzy sets have a semantic role: they represent the way in which a statement
“x is A” is used in a given context.
The generalized modus ponens used in fuzzy control may be given in its simplest
expression as
A B
A
B ,
where A, A∗ , B and B ∗ are fuzzy sets, A and A∗ are defined on a same universe, but
they are not necessarily equal. Similarly for B and B ∗ . The meaning in this case is
the following: given a rule “if A then B” and observing an event A∗ which is similar
to A, an event B ∗ is expected, which should also be similar to B.
To allow more specified situations to define a rule base in a multivariate system,
in the expression A → B, A may stand for a set of conditions that have to be fulfilled
at the same time. The formal representation is a conjunction of fuzzy premises.
For the computation of conjunctions, operations belonging to the class of trian-
gular norms, or simply t-norms, are used. If T : [0, 1]2 → [0, 1] is a t-norm, then it
is non-decreasing, associative, commutative and has 1 as identity. The other needed
operation is the “then” in “if A then B”. In the case of fuzzy logic this operation is
not only an extension of the classical implication and its characterization has been
thoroughly studied.
It becomes apparent that the meaning of the predicates freezing, cold and medium
is different if the system is intended for Ottawa, Madrid or Dakar, since they would
not be used in exactly the same way. Similarly, the meaning of a large, average or low
heating demand (in a KW scale) is different for a system meant to heat an office room,
a conference hall for 200 people or a 25 stores office building. The relative ordering
of the terms large, average and low will certainly be the same, and their shapes will
probably be the same, in all three mentioned cases. Considering one instance of the
problem as illustrated in the figure where it is assumed that the external temperature
is 6 ◦ C.
freezing cold medium low av. high
R1
R2
R3
-15 6 15
It is fairly obvious that if the external temperature is 6 ◦ C, the first rule does
not apply, since 6 ◦ is not considered to be freezing, but “more cold than medium”.
The degree of satisfaction of the premises or conditions stated in the “if” part of
rules 2 and 3 will affect the strength of the corresponding conclusions. This will be
specified by the “then” operation. Rules 2 and 3 are activated to a certain degree and
give proportional suggestions for action. These have to be combined into one single
action by means of an aggregation operation. From the many aggregation operations
that may be used for this purpose, the pointwise maximum is usually the first choice.
The aggregated result is also shown in the figure.
Since it was assumed that finally the water temperature should be proportional to
the heating demand, the fuzzy set representing the aggregation of activated heating
demands has to be converted into a real value through the deffuzzification. This
process will imply an information loss, since it is analogous to representing a signal
with only one coefficient of its Fourier power spectrum. However, experimental
results have shown that “approximating” a fuzzy set by the abscise of its gravity
center or the abscise of its center of area leads to an adequate control performance.
8.1 Introduction 179
8.1.1 Note
This chapter is divided in two main parts that on the one hand pretends to put fuzzy
control in the context of the book, and on the other hand introduces the present
situation in the area of fuzzy control.
The first part makes a revision of the fundamentals of approximated reasoning
explained throughout the book. A condensed view of how to apply the presented
theory from a fuzzy control perspective is given.
The second part outlines a glimpse of the state of the art of fuzzy control, where
fuzzy plant models and the procedures to obtain them are introduced. In the area
of fuzzy control there has been a shift from its original motivation of interpretable
fuzzy systems (systems that emulate human control strategy and are easy to use and
understand) towards a much more rigorous analysis where mainstream (nonlinear)
control criteria are considered. The reader must be aware that in order to fully under-
stand this part, he/she should probably refer to more specialized books in the topic
as this section only pretends to present a very general view of the problem.
for all x ∈ X , y ∈ Y . Fuzzy logic works mainly within the positive supposition,
and several families of such functions J exist. Such numerical functions are called
conditional functions, and their diverse types are derived from the linguistic meaning
of the conditional phrase (the rule) “If x is P, then y is Q”, that is, from its use in
the universe X × Y at each particular problem.
Even if it is some repetition of what was earlier presented, let’s summarize the
steps that are necessary.
180 8 An Introduction to Fuzzy Control
The question is the following: If we observe that the input variables are in the
“states” “x1 is P1∗ ”, “x2 is P2∗ ”, . . . , “xn is Pn∗ ”, respectively, what can be inferred
for variable y? That is, supposing that it lies in the “state” “y is Q ∗ ”, what is Q ∗ ?
Without loss of generality let us consider only the case n = m = 2:
R1 : If x1 is P11 and x2 is P12 , then y is Q1
R2 : If x1 is P21 and x2 is P22 , then y is Q2
x1 is P1 and x2 is P2
y is Q , what is Q ?
a → b ⇔ J (a, b)
for all x ∈ X , y ∈ Y .
Hence, for each type of implication function J we need to know which T1 allows
the verification of the Modus Ponens inequality:
For example, with an S-implication, since T1 (a, S(N (a), b)) b implies
(with b = 0) T1 (a, N (a)) = 0, it should be T1 = Wϕ for some automorphism
ϕ of ([0, 1], ).
With R-implications JT , since
also implies T1 (a, N (a)) = 0 for all a ∈ [0, 1], it should be also T1 = Wϕ .
Concerning ML-implications, since
because both T1 Min and T Min, the Modus Ponens inequality is verified for
all t-norms T1 and, hence, for T1 = Min (the biggest t-norm).
If T1 verifies T1 (a, J (a, b)) b, because of the well-known result that for
left-continuous t-norms T1 , T1 (a, t) b is equivalent to t JT1 (a, b), it results
that the inequality is equivalent to J (a, b) JT1 (a, b). Hence, among the func-
tions J verifying the Modus Ponens inequality with a continuous t-norm T1 , the
R-implication JT1 is the biggest one and, consequently, T1 (a, JT1 (a, b)) is closer to
b than T1 (a, J (a, b)). In particular, it is
1 if a b
J M (a, b) = Min(a, b) JMin (a, b) =
b if a > b
and
1 if a b
JL (a, b) = a · b JMin (a, b) JProd (a, b) = b
a if a > b
(since a · b b).
Third Step: Zadeh’s Compositional Rule of Inference CRI
Once the rule “If x is P, then y is Q” is represented by J (μ P (x), μ Q (y)), and a
continuous t-norm T1 such that J JT1 is known, the inference:
It should be pointed out that Zadeh’s CRI is not a “result” but a meta-rule. It is a
“directive” allowing to reach a solution to our problem, and it should be noticed
that when P ∗ = P it is not in general Q ∗ = Q. For example, in the case of ML-
implications it is:
In that case,
for all y ∈ Y .
For example, let J be a ML-implication, then μ Q ∗ (y) = T (μ P (x0 ), μ Q (y)).
If X = [0, 10], Y = [0, 1], P = close to 4, Q = big, with uses as shown in the
following figure and moreover x0 = 3.5, with J (a, b) = Min(a, b), then μ Q ∗ (y) =
Min(μ P (3.5), μ Q (y)) = Min(0.5, μ Q (y)), with μ P (x) = x − 3 between 3 and 4.
1 1
μP μQ
0.5
0
3 4 5 10 0 0.5 1
184 8 An Introduction to Fuzzy Control
0.5
μ Q∗
0 0.5 1
when x0 = 3.5.
Fifth Step: Numerical Consequent
This is the case in which μ Q is y = y0 or y ∈ {y0 }. That is “y is Q” corresponds to
“y is y0 ”, and
1 if y = y0
μ Q (y) = μ y0 =
0 if y = y0
In this case:
J (μ P (x), 1) if y = y0
J (μ P (x), μ y0 (y)) =
J (μ P (x), 0) if y = y0 ,
and the output μ Q ∗ depends on the values of J, namely, on J (a, 1) and J (a, 0).
Notice that if:
• J is an S-implication, J (a, 1) = 1; J (a, 0) = N (a).
• J is an Q-implication, J (a, 1) = S(N (a), a); J (a, 0) = N (a).
• J is an R-implication, J (a, 1) = 1; J (a, 0) = Sup{z ∈ [0, 1]; T (z, a) = 0}.
• J is an ML-implication, J (a, 1) = a; J (a, 0) = 0.
When J is an ML-implication, with J (a, 1) = a and J (a, 0) = 0:
Notice that if the input is also numerical, i.e. x = x0 then the output is:
Sup μ P ∗ (x) if y = y0 and x = x0
μ Q∗ = x∈X
0 if y = y0 and x = x0 .
8.2 Revising Conditional and Implications in Fuzzy Control 185
and let us find μ Q ∗ using both methods by supposing “close-to 4” as in the former
example, μbig (y) = y, μsmall (x) = 1 − 10
x
, and μsmall (y) = 1 − y.
∗ ∗
The outputs Q 1 , Q 2 are:
IF μ close-to 4 THEN μ big
1 1
μ Q∗1
R1: 0.5
0
3 4 5 10 0 0.5 1
0.65 μ Q∗2
R2:
0 1
3.5 10 0 0.35
186 8 An Introduction to Fuzzy Control
0.65 μ Q∗
0.5
0 0.35 0.5 1
μ big μ small
1 1 1
0.65 0.65
Seventh Step: A More Complex Example with Numerical Inputs and Conse-
quences
Let us consider the case:
With Larsen’s method and translating the and in the antecedents also by T = Prod.
• Rule R1 is represented by JL (T (μ P11 (x1 ), μ P12 (x2 )), μ y1 (y))
μ P11 (x1 ) · μ P12 (x2 ), y = y1
= μ P11 (x1 ) · μ P12 (x2 ) · μ y1 (y) =
0, y = y1
8.2 Revising Conditional and Implications in Fuzzy Control 187
Consequently:
⎧
⎨ μ P11 (x1∗ ) · μ P12 (x2∗ ), y = y1
μ Q∗ = Max(μ Q 1 (y), μ Q 2 (y)) μ P21 (x1∗ ) · μ P22 (x2∗ ),
∗ ∗ y = y2
⎩
0, otherwise
the most popular method of defuzzification is that consisting in taking the weighted
mean
α1 · y1 + α2 · y2 + · · · + αn · yn
μ Q∗ =
α1 + α1 + · · · + αn
The case of m rules with numerical consequents and numerical inputs, with
Larsen’s implication function and with defuzzification by the weighted mean, is the
188 8 An Introduction to Fuzzy Control
Example 8.2.1 In the example shown in sixth step where the output μ Q ∗ is obtained
through Mamdani, the area below μ Q ∗ is easily computed by 0.35 × 0.15 +
0.15 × 0.15
2 + 0.5 × 1 = 0.564. Hence, the center of area is a point y0 ∈ (0, 1)
such that, the areas to the left and to the right of y0 are equal, i.e.:
y0
0.564
= 0.282 = μ Q ∗ (y)dy
2
0
0.35 y0
= 0.65dy + (1 − y)dy
0 0.35
y0
= 0.228 + y0 − 0.35 − ydy
0.35
as the line joining the points (0.35, 0.65) and (0.5, 0.5) is z = 1 − y. Hence:
y0
y0 − ydy = 0.282 − 0.228 + 0.35 = 0.404
0.35
y0
y2 y02 0.352
y0 − = y0 − − = 0.404
2 0.35 2 2
μ P11 (x1∗ )μ P12 (x2∗ 2)y1 + μ P21 (x1∗ )μ P22 (x2∗ )y2
y0 =
μ P11 (x1∗ )μ P12 (x2∗ ) + μ P21 (x1∗ )μ P22 (x2∗ )
Provided that X 1 = X 2 = [0, 1], Y = [0, 10], μ P11 (x1 ) = x1 , μ P12 (x2 ) = 1 − x2 ,
y1 = 6, μ P21 (x1 ) = 1 − x1 , μ P22 (x2 ) = x2 , y2 = 4, x1∗ = 0.3 and x2∗ = 0.7, the
calculation will be:
0.3 × (1 − 0.7) × 6 + (1 − 0.3) × 0.7 × 4 0.54 + 1.96 2.5
y0 = = = = 4.31
0.3 × (1 − 0.7) + (1 − 0.3) × 0.7 0.09 + 0.49 0.58
8.2 Revising Conditional and Implications in Fuzzy Control 189
As rule R1 is represented by J1 = μ P11 (x1 ) · μ P12 (x2 ) · μq1 (y) and rule R2 by
J2 = μ P21 (x1 ) · μ P22 (x2 ) · μq2 (y), it follows:
μ Q ∗1 (y) = Sup Min(μx1∗ x2∗ (x1 , x2 ), μ P11 (x1 ) · μ P12 (x2 ) · μq1 (y))
x∈X, y∈Y
μ P11 (x1∗ )μ P12 (x2∗ ), if q1
= μ P11 (x1∗ )μ P12 (x2∗ )μq1 (y) =
0, otherwise
μ P21 (x1∗ )μ P22 (x2∗ ), if q2
μ Q ∗2 (y) = μ P21 (x1∗ )μ P22 (x2∗ )μq2 (y) =
0, otherwise
Hence,
⎧
⎨ μ P11 (x1∗ )μ P12 (x2∗ ), if q1
μ Q ∗ (y) = Max(μ Q ∗1 (y), μ Q ∗2 (y)) = μ P21 (x1∗ )μ P22 (x2∗ ), if q2
⎩
0, otherwise
190 8 An Introduction to Fuzzy Control
And finally:
μ P11 (x1∗ )μ P12 (x2∗ )q1∗ + μ P21 (x1∗ )μ P22 (x2∗ )q2∗
y0 =
μ P11 (x1∗ )μ P12 (x2∗ ) + μ P21 (x1∗ )μ P22 (x2∗ )
If the dynamical system is linear and time invariant, the differential and algebraic
equations may be written in a matrix form, such as,
ẋ(t) = Ax(t) + Bu(t)
y(t) = C x(t)
Example 8.3.1 Taken a mechanical system of a mass with a spring and a damper as
an example,
y(t)
k
u
m
where an external force is the only input u(t) to the system, and the position y(t) is
the output. The physical equation of this single input single output (SISO) system is
8.3 Control of Nonlinear Systems 191
m ÿ + b ẏ + ky = u
which represents a second order linear system. In this case the state variables are
defined as
x1 (t) = y(t)
x2 (t) = ẏ(t)
y(t) = x1 (t).
ẋ1 0 1 x1 0
= + 1 u,
ẋ2 − mk − mb x2 m
x1
y= 10 .
x2
Both state and output equations form the matrix state space representation of the
system
ẋ = Ax + Bu
y = Cx
with
0 1 0
A= , B= , C= 10 .
− mk − mb 1
m
The so-called Takagi-Sugeno (TS) fuzzy model is particularly useful for the control
of nonlinear systems. The TS fuzzy model is described by fuzzy “if-then” rules which
represent local linear input-output relations of a nonlinear system in a state-space
192 8 An Introduction to Fuzzy Control
form. The overall fuzzy model of the system is achieved by fuzzy “blending” of the
linear systems. A TS fuzzy model representing a continuous dynamical system has
the following form,
Rule i:
IF z 1 (t) is Mi1 and · · · and z p (t) is Mi p
ẋ(t) = Ai x(t) + Bi u(t)
THEN
y(t) = Ci x(t)
for i = 1, 2, . . . , r .
Here, Mi j is the fuzzy set and r is the number of model rules; x(t) ∈ R n is the state
vector, u(t) ∈ R m is the input vector, y(t) ∈ R q is the output vector, Ai ∈ R n×n ,
Bi ∈ R n×m , and Ci ∈ R q×n ; z(t) = {z 1 (t), . . . , z p (t)} are known premise variables
that may be functions of the state variables and/or external disturbances.
Each linear consequent equation represented by Ai x(t) + Bi u(t) is called a “sub-
system”.
ẋ(t) = ri=1 h i (z(t)){Ai x(t) + Bi u(t)}
y(t) = ri=1 h i (z(t)) Ci x(t)
where,
wi (z(t))
h i (z(t)) = r ,
i=1 wi (z(t))
p
wi (z(t)) = Mi j (z(t))
j=1
for all t. The term Mi j (z(t)) is the grade of membership of z j (t) in Mi j . We have
r
i=1 h i (z(t)) = 1, and h i (z(t)) 0, for i = 1, 2, . . . , r and all t.
Thanks to the normalized membership functions, the linear dynamic TS model
is in fact a convex combination of local linear models. This property facilitates the
stability analysis of the fuzzy system. Fuzzy systems are universal function approx-
imators and hence can be used to model a wide class of processes.
Example 8.3.2 Consider the following nonlinear dynamic system
⎧
⎨ ẋ1 = −x1 + x1 x2
ẋ2 = x1 − 3x2
⎩
y = x1
with x1 , x2 ∈ [−1, 1]. This system can be exactly represented, using the sector
nonlinearity approach. In this model, the scheduling variable z 1 is chosen as x2 , the
fuzzy sets are M11 = ‘around − 1’, M12 = ‘around 1’, and the corresponding
membership functions are h 11 = (1 − z 1 )/2 and h 12 = (1 + z 1 )/2, respectively. The
following TS fuzzy system with linear consequents is the exact representation of the
original model,
8.3 Control of Nonlinear Systems 193
R1 :
IF z 1 is around − 1
⎧
⎨ −2 0
ẋ = x
THEN 1 −3
⎩
y = x1
R2 :
IF z 1 is around 1
⎧
⎨ 00
ẋ = x
THEN 1 −3
⎩
y = x1
It can be easily seen that this is an exact representationof the nonlinear system
2
in the compact set x1 , x2 ∈ [−1, 1], when expressions i=1 h i (z(t)) Ai x(t) and
2
h
i=1 i (z(t)) C i x(t) are developed,
1 − x2 −2 0 x1 1 + x2 0 0 x1 −x1 + x1 x2
+ =
2 1 −3 x2 2 1 −3 x2 x1 − 3x2
1 − x2 1 + x2
x1 + x1 = x1 = y.
2 2
The main idea of the controller design based on TS fuzzy systems is to derive
each control rule so as to compensate each subsystem or rule of the fuzzy system.
So, for every subsystem
ẋ = Ai x + Bi u
u = −Fi x
ẋ = (Ai − Bi Fi )x.
The design procedure is conceptually simple and natural. The stability analysis and
control design problems can be reduced to linear matrix inequality (LMI) problems.
194 8 An Introduction to Fuzzy Control
with P = P T > 0. If this Lyapunov function is considered, its derivative along the
trajectories of the TS fuzzy subsystem of the type ẋ = Ai x is
and so:
Theorem 8.3.3 The equilibrium of a TS fuzzy system is globally asymptotically
stable if there exists a common positive definite matrix P such that
for i = 1, 2, . . . , r .
That is, a common P has to exist for all subsystems. This theorem reduces to the
Lyapunov stability theorem for continuous linear systems when r = 1. If there exists
a P > 0 such that V (x(t)) proves the stability of the system, it is also said to be
quadratically stable as V (x(t)) is called a quadratic Lyapunov function.
The theorem presents a sufficient condition for the quadratic stability of the TS
system, and this is a common P problem that can be solved efficiently via convex
optimization techniques for LMIs. For systems and control, the LMI optimization is
particularly useful due to the fact that a wide variety of system and control problems
can be recast as LMI problems. Apart from a few special cases these problems do
not have analytical solutions. However, through the LMI framework they can be
efficiently solved numerically in all cases. Therefore, recasting a control problem as
an LMI problem is equivalent to finding a “solution” to the original problem.
8.3 Control of Nonlinear Systems 195
for i = 1, 2, . . . , r .
The fuzzy control rules have a linear controller (state feedback laws in this case)
in the consequent parts. The overall fuzzy controller is represented by
r
u(t) = − h i (z(t)) Fi x(t).
i=1
The fuzzy controller design consists on determining the local feedback gains Fi in
the consequent parts. With PDC we have a simple and natural procedure to handle
nonlinear control systems. Other nonlinear control techniques require special and
rather involved knowledge.
The overall closed-loop system using the PDC controller method is,
r
r
ẋ(t) = h i (z(t))h j (z(t)){Ai − Bi F j }x(t)}.
i=1 j=1
For this case, the following sufficient condition for stability can be obtained:
Theorem 8.3.4 The equilibrium of a fuzzy control system is globally asymptotically
stable if there exists a common positive definite matrix P such that
where x1 (t) is the angle of the pendulum, x2 (t) is the angular velocity; g is the
gravity, m is the mass of the pendulum, M is the mass of the cart, l is the length of
the pendulum, u(t) is the force applied to the cart and a = 1/(m + M).
A fuzzy model which approximates the dynamics of the nonlinear plant for the
range x1 ∈ (−π/2, π/2) can be realized by the following two rules,
where,
0 1 0
A1 = g B1 = g
4l/3−aml 0 − 4l/3−aml
0 1 0
A2 = 2g B2 = aβ
π(4l/3−amlβ 2 )
0 − 4l/3−aml
R2 R1
x1
-90 0 90
Through the PDC control, the stability of the system is guaranteed for example
with
F1 = [−120.67 − 22.67] ,
F2 = [−2551.6 − 764]
where it is obtained
3.62 0.62
P= .
0.62 0.28
8.3 Control of Nonlinear Systems 197
The extension of the pendulum control to the full circle x1 ∈ [−π, π] work space,
would be achieved similarly, only by adding two more rules to the TS fuzzy system.
Recently, Piecewise Bilinear (PB) model is being used for control purpose. The
PB model is a fully parametric model to represent Linear/Nonlinear systems. The
obtained model is built on piecewise rectangular regions, and each region is defined
by four vertices partitioning the state space. As the conventional nonlinear system
control based on TS fuzzy model represents a connection of linear state-space models
by sector nonlinearity, the PB model represents a convex combination of the vertices
defining piecewise regions.
In this approach, bilinear functions are used to regionally approximate any given
function. A bilinear function is a nonlinear function of the form y = a + bx1 +
cx2 + d x1 x2 , where any four points in the three dimensional space are spanned with
a bi-affine plane.
PB model has a good general approximation capability and it has a continuous
crossing over the piecewise regions. Its interpretability, simplicity and visibility facil-
itates the realization of controllers in industrial applications. A local error does not
trigger a global error and its interpolation nature generates robust outputs. A draw-
back of the PB model is that the stability analysis based on Lyapunov is difficult as
bilinear matrix inequalities (BMI) must be solved.
If a general case of an affine two-dimensional nonlinear control system is consid-
ered, ⎧
⎨ ẋ1 = f 1 (x1 , x2 )
ẋ2 = f 2 (x1 , x2 ) + g(x1 , x2 ) · u
⎩
y = h(x1 , x2 )
where i ∈ (1, . . . , n 1 ) and j ∈ (1, . . . , n 2 ) are integers, and d1 (i) < d1 (i + 1),
d2 ( j) < d2 ( j + 1). The PB models are formed by matrices of size (n 1 × n 2 ), where
n 1 and n 2 represent the number of partitions of dimension x1 and x2 respectively.
Each value in the matrix is referred to as a vertex in the PB model. The operational
region of the system is divided into (n 1 − 1 × n 2 − 1) piecewise regions that are
analyzed independently.
198 8 An Introduction to Fuzzy Control
IF x is W στ , THEN ẋ is f (σ, τ )
which in a two-dimensional case, x ∈ R2 is a state vector, W στ = (w1σ (x1 ), w2τ (x2 ))T
is a membership function vector, f (σ, τ ) = ( f 1 (σ, τ ), f 2 (σ, τ ))T ∈ R is a singleton
consequent vector, and σ, τ ∈ Z are integers (1 σ n 1 , 1 τ n 2 ) defined by,
where ⎧ σ
⎪
⎪ w (x1 ) = 1 − α,
⎪ 1σ+1
⎨ w1 (x1 ) = α,
⎪ τ
⎪ w2 (x2 ) = 1 − β,
⎪
⎩ τ +1
w2 (x2 ) = β,
8.3 Control of Nonlinear Systems 199
and
x1 − d1 (σ)
α=
d1 (σ + 1) − d1 (σ)
x2 − d2 (τ )
β=
d2 (τ + 1) − d2 (τ )
j
in which case w1i , w2 ∈ [0, 1].
In every region of the PB models, i.e.: f 1 (x1 , x2 ), the values are computed through
bilinear interpolation of the corresponding four vertexes as shown in the figure.
Note that the approximation is made by only using the values of a nonlinear
function at the vertexes of Ri j ’s.
if a controller σ+1 τ +1
j
u = C(x1 , x2 ) = w1i (x1 )w2 (x2 )C(i, j)
i=σ j=τ
Example 8.3.6 Considering the so called “Van der Pol” oscillator, which is a non-
linear plant expressed as,
x˙1 = x2
x˙2 = −x1 + (1 − x12 )x2 + u
The following figures show both the original surface of the Van der Pol system
and it’s piecewise bilinear representation for the vertex positions d1 and d2 .
The goal is to design a controller C so that the closed-loop system behaves like
the linear plant represented as,
x˙1 = x2
x˙2 = ax1 + bx2 + u
where a = −1 and b = −3. The PB representation of f 2d p for the linear plant is the
following (Table 8.2).
With the controller calculated as u = f 2d p − f 2 p , the following PB model of the
controller C is obtained (Table 8.3).
Next figures show both the controller C and the resulting closed-loop plant
surfaces respectively.
Bibliography
1. B. Bede, Mathematics of Fuzzy Sets and Fuzzy Logic (Springer, London, 2013)
2. R. Belohlavek, Fuzzy Relational Systems: Foundations and Principles (Kluwer Academic
Publishers, New York, 2002)
3. G. Böhme, Fuzzy-Logik: Einführung in die algebraischen und logischen Grundlagen (Springer,
London, 1993)
4. B. Bouchon-Meunier, La logique floue (PUF, Paris, 2003)
5. S.P. Boyd et al., Linear Matrix Inequalities in System and Control Theory (SIAM, Philadelphia,
1994)
6. B. Demant, Fuzzy-theorie oder die faszination des vagen (Vieweg, Brauschweig, Wiesbaden,
1993)
7. D. Dubois et al., Fuzzy Sets: History and Basic Notions (Springer, London, 2000)
8. M. Grabisch et al., Fuzzy Measures and Integrals: Theory and Applications (Springer, New
York, 2000)
9. G. Klir, B. Yuan, Fuzzy Sets and Fuzzy Logic, vol 4 (Prentice Hall, New Jersey, 1995)
10. G.J. Klir, U. St Clair, B. Yuan, Fuzzy Set Theory: Foundations and Applications (Prentice-Hall,
Upper Saddle River, 1997)
11. Z. Lendek et al., Stability Analysis and Nonlinear Observer Design Using Takagi-Sugeno Fuzzy
Models (Springer, London, 2010)
12. B. Liu, Theory and Practice of Uncertain Programming. (Springer, London, 2009)
13. R. Lowen, Fuzzy set theory (Basic Concepts, Techniques and Bibliography, Dordrecht, 1996)
14. R. López de Mántaras, Approximate Reasoning Models (Prentice Hall PTR, Upper Saddle
River, 1990)
15. D. McNeill, Fuzzy Logic: The Revolutionary Computer Technology that is Changing Our World
(Simon and Schuster, New York, 1994)
16. G. Metcalfe, N. Olivetti, D.M. Gabbay, Proof Theory for Fuzzy Logics, vol 36 (Springer,
London, 2008)
17. H.T. Nguyen, E.A. Walker, A First Course in Fuzzy Logic (CRC Press, Boca Raton, 2005)
18. A. Syropoulos, Theory of Fuzzy Computation (Springer, New York, 2014)
19. K. Tanaka, An Introduction to Fuzzy Logic for Practical Applications (Springer, New York,
1997)
20. K. Tanaka, H.O. Wang, Fuzzy Control Systems Design and Analysis: A Linear Matrix Inequality
Approach (Wiley, Danvers, 2004)
21. E. Trillas, Conjuntos Borrosos (Vicens-Vives, Barcelona, 1980)
22. E. Trillas, C. Alsina, J.-M. Terricabras, Introducción a la lógica borrosa. (Ariel Barcelona,
Spain, 1995)
© Springer International Publishing Switzerland 2015 203
E. Trillas and L. Eciolaza, Fuzzy Logic, Studies in Fuzziness
and Soft Computing 320, DOI 10.1007/978-3-319-14203-6
204 Bibliography
23. E. Turunen, E. Turunen, Mathematics Behind Fuzzy Logic (Physica-Verlag, Heidelberg, 1999)
24. X. Wang, D. Ruan, E.E. Kerre, Mathematics of Fuzziness: Basic Issues (Springer, Heidelberg,
2009)
25. J. Yen, R. Langari, Fuzzy Logic: Intelligence, Control, and Information (Prentice-Hall, Upper
Saddle River, 1998)
26. L.A. Zadeh, Fuzzy Sets and Applications: Selected Papers (Wiley, New York, 1987)
27. L.A. Zadeh, Computing with Words: Principal Concepts and Ideas (Springer, Berlin, 2012)
28. L.A. Zadeh, B. Yuan, G.J. Klir, Fuzzy sets, fuzzy logic, and fuzzy systems: selected papers by
Lotfi A (World Scientific, Zadeh, 1996)
29. H.J. Zimmermann, Fuzzy Set Theory and Its Applications, Second Revised Edition. (Springer,
New York, 1992)