1992 Fuzzy Multiple Attribute Decision Making Methods and Applications Compress
1992 Fuzzy Multiple Attribute Decision Making Methods and Applications Compress
1992 Fuzzy Multiple Attribute Decision Making Methods and Applications Compress
Managing Editors:
Prof. Dr. M. Beckmann
Brown University
Providence, RI 02912, USA
Prof. Dr. W. Krelle
Institut flir Gesellschafts- und Wirtschaftswissenschaften
der Universitat Bonn
Adenauerallee 24-42, W-5300 Bonn, FRG
Shu-Jen Chen Chin-Lai Hwang
In Collaboration with Frank P. Hwang
Springer-Verlag
Berlin Heidelberg New York
London Paris Tokyo
Hong Kong Barcelona
Budapest
Authors
Shu-Jen Chen
HTX International Inc.
115 North 4th Street
Manhattan, KS 66502, USA
Ching-Lai Hwang
Department of Industrial Engineering
Kansas State University
Manhattan, KS 66506, USA
This work is subject to copyright. All rights are reserved, whether the whole or part of
the material is concerned, specifically the rights of translation, reprinting, re-use of
illustrations, recitation, broadcasting, reproduction on microfilms or in any other way,
and storage in data banks. Duplication of this publication or parts thereof is permitted
only under the provisions of the German Copyright Law of September 9, 1965, in its
current version, and permission for use must always be obtained from Springer-Verlag.
Violations are liable for prosecution under the German Copyright Law.
© Springer-Verlag Berlin Heidelberg 1992
problems.
The basic concepts and algorithms from the classical MADM methods
have been used in the development of the fuzzy MADM methods. We give
attributes for each alternative, and (2) the rank ordering of the
methods are concerned with comparing fuzzy numbers and can be applied
example.
Most of the real-world MADM problems contain a mixture of fuzzy
and crisp data and may have a large number of alternatives and up to
is presented.
fuzzy set theory done by Dr. L.A. Zadeh. We are indebted to all the
fine scholars listed in the references who have carried out and
published their research results. Special thanks are due to
I. INTRODUCTION 1
tics, and applicability. We will also present our new fuzzy MADM
approach.
whole nation must face problems of this type, i.e., MADM problems.
Xl X2 Xn
Al xII x 12 X 1n (1.1)
[
A2 X 21 X 22 x 2n
0
Am x m1 x m2 xmn 1
where Ai' i = 1, ... , m are possible course of actions (referred to as
sources:
terms such as good, fair, poor, etc. They are qualitative data
incomplete information.
but the cost is too high, and the OM may wish to get an "approxima-
tion" of that crisp data. When the data is very sensitive (i.e.,
I LI 2 . 2 . 1 Lexicographic Method
2.2 Ordinal I 2.2.2 Elimination by Aspects
2.2.3 Permutation Method
Multiple
Attribute 2. Information 2.3.1 Linear Assignment Method
Decision on Attribute 2.3.2 Simple Additive
Weighting Method (SAW)
Making 2.3 Cardinal
2.3.3 Analytic Hierarchical
Process (AHP)
2.3.4 ELECTRE UI
2.3.5 TOPSIS
ignorance of the phenomenon since one knows only part of the facts.
individual(s).
Literature Surveys
The first attempt at applying fuzzy set theory to multi-attribute
analysis was done by Bellman and Zadeh [B8], who outlined one possible
Zadeh [Z5] who outlined the possibility of using the max-min rule to
decision matrix using a singleton, i.e., a fuzzy set with only one
element in it (see also Zadeh [Z5]). In this approach, the best
methods developed prior to 1979. The studies of Yager and Basson [Y1],
7
Yager [Y4], Jain [J2,J5], Baas and Kwakernaak [B1], and Baldwin and
fuzzy set theory and its application was done by Dubois and Prade
[D26]. They classified the fuzzy MADM into a fuzzy rating phase, in
The most updated summary of the fuzzy MADM study was done by
two-phase process. The first phase derives the fuzzy utilities which
are then compared in the second phase using a fuzzy ranking method.
Applications" [Z30] places more emphasis on fuzzy set theory and its
"Fuzzy Sets, Decision Making, and Expert Systems" [Z31] devotes itself
Classification of Methods
scores with respect to all the attributes for each alternative, and
scores. We will refer to the results of the first and second phase
using the terms "final rating" and "ranking order," respectively.
For a crisp MADM problem, the final ratings are expressed as real
numbers. The ranking order can be easily obtained by comparing these
real numbers. In this case, the main focus of MADM problem solving is
the first phase. In a fuzzy MADM problem, the performance scores of
an alternative with respect to all attributes may be expressed by
fuzzy sets. As a result, the final ratings are expressed by fuzzy
sets. obtaining the ranking order of these fuzzy sets is not a
trivial task. In this case, both phase one and phase two are
important in solving the MADM problem. That is, algorithms for
computing and comparing fuzzy final ratings must be used in sequence
to solve the fuzzy MADM problem.
As indicated by Zimmermann [Z31] fuzzy MADM methods are different
from each other in that they either consider phase one, phase two, or
both phases of MADM. It seems appropriate, therefore, to classify
fuzzy MADM methods into two categories, (1) one that focuses on the
second phase, and (2) another that focuses on either the first phase
or both the first and the second phases together. We shall refer to
the methods in the first category as "fuzzy ranking methods," and
those in the second category as "fuzzy MADM methods." The reason for
such a classification scheme is explained below. Fuzzy ranking
methods are concerned with comparing fuzzy numbers, and can be applied
to virtually all aspects of fuzzy applications (not just fuzzy MADM).
Fuzzy MADM methods are specifically designed to solve MADM problems
which contain fuzzy data. The nature of the general applicability of
fuzzy ranking methods makes it appropriate to separate them from fuzzy
MADM methods. Fuzzy ranking methods and fuzzy MADM methods are
systematically presented in Chapters IV and V, respectively.
10
(1) the comparison medium used, and (2) the technique (fuzzy or
nonfuzzy) needed to develop the comparison medium. A taxonomy of
been made in three stages: stage I, the type of comparison medium used
methods are suitable for solving a problem that has either less than
stage II. The data type allowed by each method can be: (1) all
fuzzy, (2) all fuzzy singleton, (3) all crisp, or (4) a mixture of
fuzzy and crisp. Real world MADM problems contain a mixture of fuzzy
stage III. The basic concepts of fuzzy MADM methods were derived
mainly from classical MADM methods. The classical MADM methods whose
I.Comparison Medium II.Technique Involved Ill.Approaches
Comparison
function
Yager [Y18)
Intuition
linguistic Efstathiou and Tong [E5]
expression Linguistic
approximation Tong and Bonissone (T11)
All crisp I
I Ranking methods I
General MADM Fuzzy arithmetic r Negi [N7)
method
methods.
Stage IV. The technique required to apply each fuzzy MADM
the readers and the researchers if we have omitted any relevant papers.
interested in applications other than MADM (see Table 1.3). They are:
(1) Gaines and Kohout [G2] list 1150 articles and/or books.
(2) Kandel and Yager [K11] list some 1799 entries which overlap
before 1980 which deal with fuzzy set theory and its applications.
(4) Maiers and Sherif [M4] list some 450 entries. They are
(7) Kaufmann and Gupta [K16] give 57 books pertaining to fuzzy set
theory and its applications. They also list the current major sources
where one may look in for important articles: (a) BUSEFUL, (b) Journal
(e) Japanese Working Group in Fuzzy Systems, (f) Chinese Working Group
regarding fuzzy set theory and its applications which can be seen in
example, the number of missile systems the Air Force can choose from
may be less than five; while an elite college may have over thousands
350. For example, one may use price, gas mileage, safety, leg room,
workmanship, and style to evaluate cars; while there may be over 100
a hierachical tree for about 350 attributes may look like this:
Sub Sub-Sub
X111
Xl12
X11
X12
X117
Xl
X17
X2
Attributes--- X3
X4
X5
X6
X7
Note that the number seven appears in the major and subattribute
categories. The number seven is based on Miller's [M16] theory that
seven plus or minus two represents the greatest amount of information
an observer can give us about an object on the basis of an absolute
judgement.
Conflict among attributes. Multiple attributes usually conflict
with each other. For example, in selecting a car, the higher gas
mileage might reduce the comfort rating because of the smaller
passenger space.
Incommensurable units. Each attribute has a different unit of
measurement. In the car selection case, gas mileage is expressed by
miles per gallon (MPG), comfort is expressed by cubic feet if it is
measured by passenger space, cost is indicated by dollars, safety may
be indicated in a nonnumerical way, etc.
Decision Weights. Almost all methods and/or MACM problems
require information regarding the relative importance of each
attribute. The relative importance is usually given by a set of
weights which are normalized to sum to one. In the case of n
attributes, a weight set is
and 1.
four models which may be sold to that country. The Air Force analyst
are: mach, miles, pounds, dollars (in millions), high-low scale, and
Xl X2 X3 X4 Xs X6
A1 2.0 lS00 20000 S.S average very high
[
A2 2.S 2700 18000 6.S low average
1
0 A3 1.8 2000 21000 4.S high high
A4 2.2 1800 20000 S.O average average
There are dozens of MADM methods available. Each method has its
The taxonomy of MAOM methods by Hwang and Yoon [H13] was modified
by Hwang [H15] to Fig. 2.2. Six methods were removed and three new
methods were added. The three methods added were the lexicographic
semiorder method, the weighted product method, and the distance from
target method. The methods removed were the permutation method, the
computations are quite involved and not practical for more than about
four or five alternatives. The AHP method is not a new technique, but
The other three methods removed were LINMAP, the interactive SAW
method, and the MOS with ideal point method. These three all belong
to the third major branch of the original taxonomy in which the type
one reason why these three methods were removed. The problems which
LINMAP and the MOS with ideal point method were designed for involve
reason people buy a particular car; that is, which attributes are most
Dominance
Maximin
No Information Maximax
Conjunctive Method
Standard level (Satisficing Method)
Disjunctive Method
Multiple
Attribute
Decision Lexicographic Method
Elimination By Aspect I\:)
Making Col
Ordinal Lexicographic
Semiorder
Information on
Linear Assignment
Attribute
Method
Simple Additive
Cardinal Weighting Method
ELECTRE
TOPSIS
Weighted Product
Distance from Target
Fig. 2.2 A taxonomy of MADM methods (Hwang [H15]).
24
received from the decision makers, and the salient feature of the
stage, and one of the other methods can be used for the evaluation,
One more way of classifying methods is by the data type which the
lexicographic semiorder, and EBA method apply; of data type rank where
only the dominance, lexicographic, lexicographic semiorder and Linear
the conjunctive method, SAW, ELECTRE and TOPSIS, for example, apply.
Dominance
Screen Conjunctive method
Disjunctive method
Maximin
[MADM Maximax
Lexicographic
Lexicographic
Semiorder
EBA
Evaluate! LAM
Prioritize! SAW
Select ELECTRE
TOPSIS
Weighted Product
Distance from Target
Dominance
Lexicographic
Lexicographic
Semiorder
LAM
Dominance
Maximin
Maximax
Conjunctive
Disjunctive
Lexicographic
Lexicographic
SAW Semiorder
ELECTRE
TOPSIS
Weighted Product
Distance from Target
presentation.
1) Dominance
2) Maximin
3) Maximax
5) Disjunctive method
6) Lexicographic method
11) ELECTRE
12) TOPSIS
Procedure:
Requirement: None
alternatives.
alternatives.
Procedure:
such that
{A;
.
I max min x;)'}, j
i j •
1,2, ••. ,n; i 1,2, ..• ,m.
Procedure:
Procedure:
1, or 2 or ... or n,
must be specified.
alternatives.
attributes is noncompensatory.
W1 »>W 2 »> ... »>Wn where wi is the weight of the ith most important
attribute.
Advantages: Simple, easy to use and understand.
Disadvantages: The tradeoff among attributes is noncompensatory.
References: Hwang and Yoon [H13].
33
Procedure:
value.
and select the alternative(s) with the highest or near highest value
in that attribute.
w1 »>w 2 »> •.. »>w n ' where wi is the weight of the ith most important
attribute.
Logic and Basic Principle: An alternative which has many high ranked
Procedure:
alternative such that the summation of the scores for that assignment
is maximized.
Requirement: None
on each attribute.
data.
an alternative ranked first may have a score of 100, while one ranked
A-
where x ij is the outcome of the ith alternative about the jth attribute
with a numericallY comparable scale, Wj is the importance weight of
the jth attribute.
Requirement: The attributes must be both numerical and comparable.
The decision maker assigns importance weights to attributes.
Advantages: The best known and most widely used method. Simple.
easy to use and understand. The tradeoff among attributes is
compensatory.
Disadvantages: If attributes are complementary (a high score on
one attribute always occurs with a high score on another attribute),
the computed score violates the assumption of separable utility of
each attribute.
References: Hwang and Yoon [H13], MacCrimmon [M2].
37
Procedure:
elaborate.
References: Hwang and Yoon [H13], Benayoun, Roy, and Sussman [Bga],
Roy [R10].
38
to Ideal Solution)
Logic and Basic Principle: The chosen alternative should have the
shortest distance from the ideal solution and the farthest distance
Procedure:
1) Calculate the normalized decision matrix. The normalized
/
2) Calculate the weighted normalized decision matrix. The
n
where wJ' is the weight of the jth attribute and L w. = 1.
j=1 J
3) Determine the ideal and negative-ideal solution .
A
• {(max vijl j E J) , (min vijl j E J') I i =1,2, ... ,m}
.
i
•
{V 1 ' v 2 ' ... , .
vj '
i
... , v
n
} ,
is given as:
S ,-
~ I ~
J!;'l
(V IJ - V-)2
J ' i=l,2, ... ,m.
attributes is compensatory.
Procedure:
such that
A• max
i
attribute values.
attributes is compensatory.
Logic: For some attributes, the best value may be located in the
Procedure:
distance
I i:
j =1
W2
j
(x
i j
- t)2
j'
i=1,2, ... ,m,
where x .. is the outcome of the ith alternative about the jth attribute
~]
preferred here.
attributes is compensatory.
3.1 Introduction
easily assign a person seven feet tall to the "class of tall men".
six-foot tall person to that class, because the term "tall" does not
To cope with this difficulty, Zadeh [Zl] proposed the fuzzy set
between zero and one. Precise membership grades do not convey any
subjectively assessed.
1, iff x e A,
{ (3.2)
0, iff x E A.
When U is a finite set {Xl' ... , Xn }, the fuzzy set on U may also
n
A [ xi/~A(xi)· (3.3)
i=l
44
A
(3.4)
characterized as:
Example 2 (zimmermann [Z26]). Let U = {10, 20, 30, 40, 50, 60,
70, SO, 90, 100}, the possible speed (mph) at which cars can cruise
over a long distance. Then the fuzzy set "comfortable speed for long
where
1 / {1 + [ 1/5 (x - 10)]2}.
45
fuzzy set A.
(3.5)
Since one has the most comfortable feeling at 70 mph, i.e., ~A(70) 1,
grades. These elements are the support of that fuzzy set. That is,
SeA) (3.6)
degree of a. That is, for fuzzy set A its a-cut is defined as:
(3.7)
The a-cut is a more general case of the support of a fuzzy set. When
S (A) •
AO.9 = {70}.
Example (infinite fuzzy set case): Consider again the fuzzy set
real number between [4,16]. And its a-cut at degree of 0.55 is any
real number between [7,12]. That is, the set of real numbers that
is convex if
(3. a)
where xl' x 2 e U, and A e [0,1]. Fig. 3.2 gives a convex fuzzy set
and a nonconvex fuzzy set. All the fuzzy sets in the following
48
chapters are assumed convex. For simplicity, we will use the term
and only if there are one or more x, values such that ~A(x') = 1. All
fuzzy sets in Figs. 3.1 and 3.2 are normal fuzzy sets. This property
the phenomenon that the fuzzy set applies to. Unless otherwise
stated, all the fuzzy sets in the following chapters are assumed
normal.
L ~A(x), x E U. (3.9)
IIAII (3.11)
its cardinality IAI and relative cardinality I IAI I are computed as:
49
5.15
IIAII 10
0.515.
3.2.2.7 The mth Power of a Fuzzy Set: The mth power of fuzzy
(3.12)
x 3 4 5 6 7 8 9 10
o o .2 .4 .6 .8 1.0 o
x 3 4 5 6 7 8 9 10
union can be carried over to fuzzy sets, too. Zadeh [Zl] was the
first to extend intersection and union to fuzzy sets. He used the min
50
by Bellman and Giertz [B9]. They also pointed out that from a logical
solely on its best attribute, even if all of its other attributes are
The max operator is the most suitable operator for that case.
poor attribute value, even when all its other attributes are well
operators allow the OM to express how strong s/he means "and" for a
the OM the flexibility to show how strong s/he means "or" for a
~.3.1 No Compensation H3.3.1.1 The Min Operator
Operators
~.3.2.1Algebraic Product
3.3.2.2 Bounded Product
3.3.2.3 Hamacher's Min
~.3.2 Compensation-Min ~ Operator
Operators 3.3.2.4 Yager's Min Operator
3.3.2.5 Dubois and Prade's Min
Operator
Full Compensation
----f.3.3.1 The max Operator
IFuzzy Operators
V~.3.3
-
Operators
operation.
The following operators are used on fuzzy sets A and B. For
3.3.1.1 The Min Operator: For fuzzy sets A and B, the min
(3.13 )
(3.14)
when taking the intersection of fuzzy sets in the decision space. The
this section.
(3.15)
(3.18 )
where
54
IlA, i f IlS 1,
TW(IlA'IlS)
{ IlS'
0,
if IlA
otherwise.
1, (3.19)
IlAqS (3.20)
q IlAqB
IlAAB (3.21)
55
IlAAB
° min (IlA,Il B )
1 (IlA) (Il B )
Note:
The min operator, the algebraic product, and the bounded product
may be seen as three basic operators for intersection. Each measures
different degree of "and" in the decision space. Sometimes, it is
not so obvious as to which one to use to combine fuzzy sets. This
difficulty may be alleviated by using Hamacher's, Yager's, or Dubois
and Prade's min operators because they have parameters which allows
the DM to define how strong s/he means "and" in the decision space.
By setting parameter values, the DM automatically determines an
intersection operator that may very well represent how s/he wants
to combine fuzzy sets in the decision space. The DM specified
operator may not be any of the basic intersection operators.
x 3 4 5 6 7 8 9 10
IlA 0 0 .2 .4 .6 .8 1.0 0
IlS 0 .5 .7 1.0 .7 .5 0 0
min (IlA, Il S ) 0 0 .2 .4 .6 .5 0 0
IlA( o)S 0 0 0 .4 .3 .3 0 0
Il AqS 0 0 0 .4 0 0 0 0
IlAAS 0 0 .2 .4 .6 .5 0 0
IlA (o)s(7) = max (0, IlA (7) + IlS (7) - 1) = max [0, .3] .3;
IlA1'S (7) IlA (7)IlS (7)/[.5 + .5(IlA (7) + Il S (7) - IlA (7)IlS (7»] .45;
3.3.3.1 The Max Operator: For fuzzy sets A and B, the max
or
~AvB may be seen as the union of fuzzy sets A and B (Zadeh [Zl). It
is interpreted as the "smallest" fuzzy set containing both fuzzy sets
A and B.
full) compensation when taking the union of fuzzy sets in the decision
space. The operators in this category are the algebraic sum, the
bounded sum, Hamacher's max operator, Yager's max operator, and Dubois
section.
~A+B (3.24)
/.LA7'B (3.26)
/.LA + /.L B
(3.27)
q' /.LAq'B
1 min (1, /.LA + /.LB)
CD max (/.LA,/.LB)
59
For i\' = 1, ~Ai\'B becomes the algebraic sum ~A+B. The membership
value ~Ai\'B decreases when i\' approaches 0.
Note:
The max operator, the product sum, and bounded sum can be seen as
the three basic forms of the union operator. Each measures different
degree of "or" in the decision space. It is sometimes difficult to
decide which basic form to use to combine fuzzy sets. The
parameterized operators such as Hamacher's max operator provides the
DM the flexibility to specify how strong s/he means "or" for a given
decision problem. The DM specified union operator may not be any of
the basic union operators.
Assume that 7' = q' = i\' = .5. The results of various union operators
can be summarized as follows.
60
x 3 4 5 6 7 8 9 10
o o .2 .4 .6 .8 1.0 o
o .5 .7 1.0 .7 .5 o o
1.0;
JlA (7) + Jl B (7) - JlA (7)Jl B (7) - [.5 A JlA (7) A Jl B (7)]
max [.5, 1 - JlA (7) , 1 - Jl B (7)]
.76.
the best or the worst attribute value alone. However, when the DM
wants a compromised solution (i.e., a course of action that lies
between the results of "logical and" and "logical or"), the operators
proposed in the previous sections are not applicable. Therefore, some
general operators must be used to encompass the tradeoff information.
This is a new breed of operators which aggregate fuzzy sets in a more
"generalized" format.
~AaB (3.30)
One may use the algebraic product and the algebraic sum for
intersection and union, respectively. Eq.(3.30) then becomes:
(3.31)
Example: Let
Assume that the algebraic product and the algebraic sum are used as the
intersection and union operations, respectively. When we set a = 0.5,
(i.e., the actual operator is located exactly in between "logical and"
x 3 4 5 6 7 8 9 10
~A 0 0 .2 .4 .6 .8 1.0 0
~B 0 .5 .7 1.0 .7 .5 0 0
.61.
simple answer. Yager [Y21], Dubois and Prade [030], and Zimmermann
[Z30,Z31] provide some useful rules that can be used as decision aids
restriction is better.
2, Empirical fit: The operator must be an appropriate model of
to compute.
5. Compensation: Logical "and" and logical "or" are too
fuzzy sets by the product operator will decrease the membership value
each time a new fuzzy set is introduced. Such behavior mayor may not
be desirable.
8. Required scale level of membership function: Different
which requires the lowest scale level is the most preferred one.
3.4.1 The Extension Principle (Dubois and Prade [026], Zadeh [Z6])
defined as:
x 3 4 5 6 7 8 9 10
o o .2 .4 .6 .8 1.0 o
o .5 .7 1.0 .7 .5 o o
i. e. ,
(3.33)
65
M {(y,~(y» I y (3.34)
where
-1
SUp min (~A (x1)'···'~A (x n », if f (y)~O,
{
Y = f(x 1 , .•. x n ) 1 n
~(y)
0, otherwise. (3.35)
Note that Eq.(3.35) is true only when the inverse of f is not zero,
-1 -1
i.e., f (y) ~ o. When f (y) = 0, ~(y) = o. ~(y) is the greatest
{
~A(f-1(y», if f-1(y) ~ 0,
~(y) (3.36)
0, otherwise.
2 3 4 5 6 7
o .4 1 .7 o o
o .1 .8 1 .3 o
6 789 10 11 12 13 14 15 16 17
o o .1 .1 .4 .4 .8 1 .7 .7 .3 o
that 12 = 2x 1 + x 2 are:
4 5
4 2
~A (Xl) .4 1 .7
1
~A (X 2 ) .3 .8 0
2
~A (Xl)" ~A (X 2 ) .3 .8 0
1 2
number operations are from Jain [J2], Mizumoto and Tanaka [M17,M18],
Haeringen, and DeLange [09], Gupta [G14], Kaufmann and Gupta [K15]
P.(x)
1................. .------,.
=
2. Constant on (-CD,a]: ~(X)
° vx e (-CD,a];
",(x)
1.
.9
.8
.7
.6
.5
.4
.3-
.2
.1
0 x
0 10 20 30 40 50 60 70 80 90 100
x 30 40 50 60 70 80 90
~(x) .3 .5 .8 1.0 .S .5 .3
",(x)
1.
.9-
.8
.7
.6
.5
.4-
.3
.2
.1
o ~~~~~~--~~--~-T--~~
o 10 20 30 40 50 60 70 80 90
For example, let M be the fuzzy number "about 60" which may be given
U = {10,20,30,40,SO,60,70,SO,90,100},
~(x)
different ways.
(i) Use of a-cut: Let's define the a level sets for M and N
(3.37)
and
(3.38)
(3.39)
(3.40)
One can see that Eq. (3.40) is an example of the extension principle.
For fuzzy numbers similar to the one in Fig. 3.5, we would use
numbers similar to the one in Fig. 3.6, we would use a-cut to get
70
(3.42)
x,y o 1 2 3 4 5 6
~(x) o .3 .8 1. .5 .1 o
o .6 1. .9 .4 o
z = x+y 1 2 3 4 5 6 7 8 9 10
o .3 .6 .8 1 .9 .5 .4 .1 o
x 0 1 2 3 4 5 6
Y 6 5 4 3 2 1 0
~M(x) 0 .3 .8 1 .5 .1
°
~N(Y) 0 .4 .9 1 .6
°
~M(X)A~N(Y) 0
° .4 .9 .5 .1
°
Thus, ~(+)N(Z=6) = max [0,0,.4,.9,.5,.1,0] = .9.
to 2" and N represent "real numbers close to 8" (see Fig. 3.7), where
72
~(X) 0, x :s 0,
x/2, 0 < X :s 2,
(4-x)/2, 2 < X :s 4,
0, X > 4,
J1 N (y) 0, y :s 3,
[ (y-3)/5, 3 < Y :s 8,
(11-y)/3, 8 < Y :s 11,
0, Y > 11.
JL(X) M N
1.
.9 '.M{+)N
.8
.7
.6
.5
.4
.3-
.2
.1
0 x
0 2 4 6 8 10 12 14 16 18 20
M and N are:
That is, at some a level, the X can be either m1 or m2 , and y can take
0, z :s 3
(z-3)/7, 3 < z :s 10
(-z+15)/5, 10< z :s 15
0, z > 15
(3.43)
(3.44)
(3.45 )
74
of N to M, M(+) (-N).
Because a negative number may appear as a result of subtraction,
z -5 -4 -3 -2 -1 o 1 2 3 4 5
o o .3 .4 .8 .9 1 .6 .5 .1 o
Note that (-5,0) may be dropped from the fuzzy set since, by the
x 0 1 2 3 4 5
Y 1 2 3 4 5 6
~M(Y) 0 .3 .8 1 .5 .1
~N(Y) .6 1 .9 .4 0
~M(x)A~N(Y) 0 .3 .8 .1 0
[5a-11, -7a+1].
0, x ~ -11,
(z+11)/5, -11 < x ~ -6,
(1-z)/7, -6 < x ~ 1,
0, x > 1.
lUx)
1.
-12-10 -8 -6 -4 -2 0 2 4 6
~(x) = 0, V x < 0
and
~N(Y) = 0, V Y < O.
1. At the left, we take into account all pairs (x,y) such that
= max (3.47a)
xy~z
2. At the right, we take into account all pairs (x,y) such that
= max (3.47b)
xy~z
for Which ~(.)N(Z) = 1. This will show what value of z occurs when
we pass from the left to the right of the peak.
(3.48)
77
Eqs.(3.47) and (3.48) are equivalent. This can be easily proved (as
a. When both M and N have the same sign, M(o)N can also be a
M(o)N = N(o)M,
and
(M(o)N) (o)K = M(o) (N(o)K).
1 1
[--m-' -m--]·
2 1
Z = M(o)N:
z 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
o .3 .6 .6 .8.8 1.9.9.9 .5 .5 .5 .4 .4 .4 .4 0
1.
78
x\y 1 2 3 4
1 .3 .3 .3 .3
2 .6 .8
3 .6
4 .S
we have:
x 4 S S
y 4 4 3
.4 .1 .1
Thus, ~M(.)N(14) =max [.4,.1,.1] = .4. Note that the pairs such as
14:Sxy
(~M(3) A ~N(S» are dropped from evaluation because ~N(S) = o.
10a 2 + 6a - z 0, (3.49)
and
6a 2 - 34a + 44 - z 0. (3.50)
and
a = (34 - (100 + 24Z)0.5)/12,
0, z :s 0,
Clearly, ~(')N(Z) is still a fuzzy number even though its left and
P.(X)M N
1.
.9
.8
.7
.6- 'M(.)N
.5
.4
.3
.2
.1 ""
00 5 10 15 20 25 30 35 40 45 50 X
(3.53)
be written as:
or
[~, ___
1 __ ] (Kaufmann and Gupta [K15]). (3.55)
2 n1
(3.56)
and
1 1
[m 1 (-n--)' m2 (-n--) ] .
2 1
altogether.
pretty well in fuzzy cases. It has been noted by Dubois and Prade
[022,026] and Zimmermann [Z30] that given three positive fuzzy numbers
M, Nand K,
Gupta [K15]):
Dubois and Prade [022,026]] pointed out that the fuzzy max is the
dual operation with respect to union, while the fuzzy min is the dual
foregoing statement the fuzzy max and the fuzzy min as:
(3.59)
(3.60 )
respectively, or
(3.61)
(3.62 )
respectively. Graphically, the fuzzy max and the fuzzy min are
p.(x)
M N
1.
JJ.(X)
M N
1.
The properties of the fuzzy max and the fuzzy min are summarized
as follows:
operations.
e. Idempotence:
x,y 1 2 3 4 5 6
~(x) 0 .7 1 .4 .2 0
JlN(y) 0 .3 1 .6 0
z = xvY 1 2 3 4 5 6
Jlmax(Z) 0 .3 1 .7 .2 0
For instance, to obtain Jlmax (Z=3), the (x, y) pairs that satisfy
3 = xvY are:
x 1 2 3 3 3
y 3 3 3 2 1
85
~(x) 0 .7 1 1 1
~N(Y) 1 1 1 .3 0
z = x"y 1 2 3 4 5 6
o .7 1 .4 o
except that the fuzzy min uses z = x"y but the fuzzy max uses z = xVY.
Ma [4a + 1, 8 - 3a],
and
Na [a + 2, 9 - 6a].
''(x)
1.
.9-
.8
.7
.6
-\" N
0
0
M
~(V)N(Z) 0, z ~ 2
(z-2)/l, 2 < z S 2.33
(z-1)/4, 2.33 < z ~ 5
(8-z)/3, 5 < z S 7
(9-z)/6, 7 < z ~ 9
0, z > 9.
L«m-x/a),xsm,a>O,
~(x) { (3.63)
R«x-m/~),x~m,~>O,
where m is the "mean" of fuzzy number M and a, ~ are the left and
right "spreads", respectively. When a = ~ = 0, M is considered a
crisp number m. It is often written as (see Fig.3.l3a)
M (m, a, ~).
If the peak is not unique, the L-R number M has a flat region.
It can be written as (see Fig. 3.l3a):
88
P.lx) . M M'
1....................................... .
P.lx) M M'
1.' ..................................... .
~(x) = 0, X :s 1,
(x - 1) / (m - 1), 1 < x s m, (3.64)
(u - x) / (u - m), m < x:s u,
0, x > u.
89
In Fig. 3.13b, M = (1, m, u) with 1 and u being the lower and upper
M (a, b, c, d)
with the [b,c] interval being the most likely values for M and any
more restricted form than the L-R fuzzy number, in that all "legs"
M (l,m,u) (m,ex,/3)
Tables 3.5 and 3.6 summarize the algebraic operations for triangular
Since the L-R fuzzy number (m,ex,/3) is a special case of the L-R
Inverse of N: N- 1 (n -1 , on -2 , '1n
-2
) (3.66 )
Addition: M (+) N (m+n, a+'1, ~+o) (3.67)
subtraction: M (-) N = (m-n, a+o, ~+'1) (3.68)
Multiplication:
Scalar Multiplication:
M < 0, N < 0: M ( : ) N
m
(ii'
-n~ - m'1 -na - mo)
) (3.76)
2
n2 n
91
a b ca - a7 d~ - b05
M < 0, N > 0: M ( : ) N (c' d' (3.85)
C(C-7)' d(d+o»
b a -b7 - c~ , -a05 - da
M < 0, N < 0: M ( : ) N (c' d' d(d+05) ) (3.86)
C(C-7)
92
(3.99 )
Inverse of N: N- 1 (3.100)
Scalar Multiplication:
Multiplication:
Division:
a1 b1 c1 d1
M > 0, N > 0: M ( : ) N (d' c 2 , b 2 , -
a 2)
(3.108)
2
d1 c1 b1 a1
M < 0, N > 0: M ( : ) N (d' c 2 , b 2 , -
a2
) (3.109)
2
d1 c1 b1 a1
M < 0, N < 0: M ( : ) N (a' b 2 , c 2 , -
d2
) (3.110)
2
94
Using the same figure, we can define the left slope as:
PIx)
-M 1. M
-d -b--c -a
S ( dy 5y 1.
- 1) , Y s -d
-1
Notice that N is no longer a fuzzy number. If we consider only the
dy - 1 ~ 1 5
5y (y - d I d(d+5»'
96
-1
then N can be approximated as
1 1 c'l 7
(a' ~, d(d+c'l)' C(C-7))·
left slopes first, let x and y be two unique real numbers, such that
S«a-x)/a) w S«c-y)/-r)
x a - a S-l(W),
which implies
-1
z = x + y a + c - (a+7) S (w)
where S( (a+c) - z) = w. The same reasoning holds true for the right
a + 7
slopes of M and N, where
S(z - (b+d) w.
(Hc'l
IlM (+)N(Z) ! S(
S(
1,
(a+c)-z)
a+7 '
z ~ a+c, a+7 > 0,
otherwise.
97
(a). We can use the same reasoning to obtain the formula for
subtraction:
Using the same reasoning as for addition, when M > 0 and N > 0,
-1
z = xy ac - (ar + ca) S (w) + ar (S -1 (w» 2
.
One may neglect the second-order equation in S-l(w) when a and rare
The membership function defined on the right hand side of Eq. (3.81)
«a-a) (c-r) ,S(l», and «b+/3) (dH), S(l». Similarly, when M < 0,
98
the identity
M(:) N
where lIN is the inverse of (N defined in part (b». The formula for M
The formulas for M < 0, N > 0 and M < 0, N < 0 can all be derived in
a similar manner.
Remarks
may preserve certain properties very well but some other equations may
preserve other properties very well. The equations are summarized as:
Properties Equations
1
M (+) N M (.) N Eqs. (3.78), (3.81), (3.82), (3.83)
-1 1
M = M
Eq.(3.78)
a b a~+ba
, a~+ba
M > 0
(
b' a' b(b+~) a(a-a»'
M (.) (!)
M { a' b a
h'
-a~-b
a(a-a)
, -a~-b
b(b+~»'
M < 0
(3.112)
3.5 Conclusions
Fuzzy sets are used to model imprecise or vague concepts in the
set-theoretic operators, one can combine various fuzzy sets and derive
mathematical models.
A fuzzy number is a fuzzy set with its elements confined in real
4.1 Introduction
Recall that in Chapter 2 we define the Multiple Attribute
Xl X2 Xn
A1 x 11 x 12 x 1n
[
A2 x 21 x 22 x 2n
D
Am x m1 x m2 x~ 1
~ (w 1 , w2 '···, Wn )
final ratings.
In reality, the alternative performance rating x ij can be crisp,
among various ranking methods. Lee and Li [L3,L7] and McCahone [Mc1]
theory. Similar remarks can also be found in Dubois and Prade [D29].
Li's [L3] classifications. Tseng et al. [T13] also gave some general
ranking results.
three major classes according to the means (or media) each method
uses. The readers should refer to the taxonomy shown in Fig. 4.1.
There are preference relation methods, a fuzzy mean and spread method,
each section.
Basic Notations
Comparison
function
Yager [Y18J
Intuition
linguistic Efstathiou and Tong [ESJ
expression Linguistic
approximation Tong and Bonissone IT11J
L-
Fig. 4.1 A taxonomy of fuzzy ranking methods.
105
that fuzzy final ratings Ui , vi, are readily available for comparison.
The manner in which the fuzzy final ratings were derived shall be
presented in Chapter 5.
Example Set
important issue because it is not clear whether "M > Nil or liN > M."
P.(x)
N M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
p.(x)
N' M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
p.(x)
M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
cannot discriminate between M and N. However, Lee and Li's [L3] and
McCahone's [Mel] methods give N > M, based on the assumption that our
intuition favors fuzzy sets with higher means and smaller spreads.
/L(X)
N M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~--~-+~~~~~--~--~~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
(5) Example 5: When two fuzzy sets share the same "right leg",
/L(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~--~~--~--~~--~--~~ x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
/L(x) M1 M2 M3
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
/L(x) M1 M3
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
ranking methods.
/L(x) U2 U3 U1
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
(9) Example 9: This is a case where all fuzzy sets have the same
spread.
/L(x)
M1 M2 M3
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
(10) Example 10: The fuzzy sets U1 ' U2 , and U3 are from
p,(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~~~~~~--~~~-4--~--~~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
(11) Example 11: In this case, four fuzzy sets are compared.
P,(x)
Ml
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
P,(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~-+--~~~~--4---~-+--~~~~ x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
methods, such as Adamo [A1] and Buckley and Chanas's [B40) methods are
best for screening. Any fuzzy set whose peak is lower than the a
level set by the DM should be eliminated. The remaining ones are then
p,(x)
~ M3 M5
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
(14) Example 14: The fuzzy sets in this example are discrete.
P.(x)
1. M •
.9 N *
.8
.7
.6
.5
.4
.3
I
.2
.1
0 ~-+--~~~-+--~--~-+--~--r-~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
p.(x)
1. M 0
.9 N *
.8
.7
.6
.5 *0
.4
.3
.2
.1
0 x
0 2 3 4 5 6 7 8 9 10
Remark
Many of the methods presented here are not perfect. But they
have shown the process of human efforts to find ways to solve
problems. There are always some good points coming out of each
method. We wish to recognize the effort and process of improvement.
A flawless ranking method may possibly be obtained by combining some
of these good ideas into one algorithm.
114
preferred. Baas and Kwakernaak [B1] first determines the best fuzzy
set and then establishes a preference index which denotes the degree
of preferability of this fuzzy set over other fuzzy sets. Baldwin and
derived for each preferred fuzzy set. The preference score denotes
Algorithm
1, if r i ~ r j , Vj E I,
{ (4.1)
0, otherwise.
(4.2)
This membership function indicates the degree to which the ith fuzzy
set is ranked first. For example, in Fig. 4.2 ~0(1) = 0.5, the number
0.5 characterizes the extent to which M1 is ranked first while ~0(2)
(2) which fuzzy set is better when there is more than one fuzzy
116
IL(X) M M
1 ............................ t ......... 2
. 1Lo(2)
.9
.8
.7
.6 IL. (1)
.5 ···0···- ................. .
.4
.3
.2
.1
o ~~--~-4--~~~~--+_~--4_~x
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.2 Example of ~O(i), i = 1,2.
defined as:
where
m
L (4.5)
j=l
j"i
over all other fuzzy sets. When the number of fuzzy sets to be
compared is small, (less than 10), step 3.1 may be used. Otherwise,
step 3.2 is recommended. Note that either step 3.1 or step 3.2 should
step 3.1
(4.6)
For fuzzy set Mi , there may be more than one ~i value that satisfies
Eq.(4.6). For example, for Ml in Fig.4.3, the ~l values that satisfy
~ (r 1 ) a o are 4.5 and 5.5; while when ~2 = 6.5 and 7.5, we have
1
~M (r 2 ) ao ·
2
To obtain the greatest lower bound and the smallest upper bound
[v. ,v. ], we simply apply Eq.(4.5) on all possible
l.min l.max
P.lx)
1.
.9
.8 ~Q •••....•••••••.••••••
.7
.6
.5
.4
.3
.2
.1
o x
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
•
Fl.g. 4.3 The rA 1 , rA 2 values at a o .
118
1lM. (r j ) d~M
j
(r.) /dr. , Vj .. i,
J J
J
and
~. (~i) d~
i
(r.),
l.
l.
Numerical Example
The described algorithm will be applied to Example 10. The
fuzzy final ratings to be compared are: U1 = (.20,.30,.50), U2 = (.17,
.32, .58), and U 3 = (.25, .40, .70).
P.(x) U, U2 U3
1.
0(3)
.9
.8
.7
0(1)
.6
.5
.4-
.3
.2
.1
0 X
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
i 1 2 3
1 1 1
For example, in Fig. 4.4, the variable r 1 may take values from the
interval [0.20,0.50], Le., r 1 may be .20, .25, .30, .35, .40, .45,
take any real number in [0.17,0.58] and r3 can take any real number in
i 1 2 3
.71 .81 .1
In Fig. 4.4, the dashed line indicates the intersection of ~ (r.), Vj.
j J
That is, the dashed line represents
deriving fuzzy set (v 3 ,Jl p (V 3 )). This fuzzy set may be seen as the
3
preference score of U3 .
since all three fuzzy final ratings are piecewise linear and the
score shall be obtained using step 3.1. We shall set a o = .50 as the
A A
r2 r3 V3
.25
[.245 [
.325
.55
.0775
.3025 (= V3
max
)
.325 .0200
.46 [ .55 .2450
.40
[.245 [
.325
.55
.0025
.2275
.325 -.1050 (= v3 . )
.46 [ .55 .1200
m1n
r-~--~--r-~--4---r-~--~--r--4 x
-.5 -.4 -.3 -.2 -.1 0 .1 .2 .3 .4 .5
and ~U(~l) and ~U(~2) have the same sign and ~U(~3) has the opposite
1 2 3
sign.
~U (r 2 ) 6.6,
2
~U (r 3 ) 6.6.
3
v3 . or v3
m~n max
10 6.6 -3.3
Note
1. McCahone [Mc1] found that when ~o(i), Vi, are piecewise
linear, ~o(i) can be calculated by finding the highest membership
function value associated with the points on the x-axis where a fuzzy
set's membership function values exceeded the others. For example, in
Fig. 4.2, ~o(l) = 1 because the value 1 is the highest membership
function value associated with the points on the x-axis where MI'S
membership values exceed the membership values of M2 . The ~o(i) index
compares fuzzy sets based solely on the locations of their peaks.
This may cause an illogical result (See Note 3).
2. In the numerical example presented earlier, all ~U (x.) are
i l.
piecewise linear. The ~o(i) index, vi, can be calculated by taking
the intersection of the "left leg" of one alternative with the "right
leg" of another. The result is U3 > U2 > U1 • This is the same as
completing Eqs.(4.I) through (4.3).
3. Baldwin and Guild [B3] pointed out that Baas and Kwakernaak's
method may generate counter-intuitive results. For example, in Fig.
4.7a MI is ranked first by this method but our intuition would favor
M2 because M2 's mean value is greater than the mean value of MI.
Baas and Kwakernaak's method also does not have high
124
discriminating ability. Lee and Li [L3] pointed out (see Fig. 4.7b)
that Ml and M2 (which share the same peak), are considered equivalent
suggest that our intuition would favor the fuzzy set which has a
/L(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~-+--~--~-+--~~~~--~--r-~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
/L(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~--~~--~--r--T--~--r--+x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
of the implication:
x => Y (4.7)
where
min
1,
{
0, otherwise.
Numerical Example
obtain
where
126
min
1,
{
0, otherwise.
X1 \X 2 1 2 3
1 .19 o
ll(X - > Y) 2
3
[ .29
.29 .19
o
1
Let Zij e [0,1] be an element of the matrix. It denotes the strict
dominance of Ui over Uj. When Zij = 0, this means that Ui would never
dominate Uj' while Zij = 1 shows that Ui dominates Uj' crisply. Based
Note
where X : "fuzzy set M. and fuzzy set M.", Y: "fuzzy set M; strictly
~ J ~
where IJ. p (x.,x.) varies according to the OM's attitude toward risk.
ij 1 J
When the OM is risk-averse, we have
m
min ~o. (i) , (4.16)
j 1 )
j '" i
m
~O(i) min
j 1
j '" i
~ (x.)] 1\ (4.18 )
j )
z ~. (Xi) ,
~
Z ~.(Xj)'
)
Z ~P .. (Xi,X j ),
~)
129
and then
J1 0 ( i) min (Z).
j
If the membership functions J1M (x.), ~ (x.), and J1 p (x.,x.) are all
i ~ j ) ij ~ )
piecewise linear functions which may be characterized by Fig. 4.8,
0 - x.~
Z
0 - '1
J1 M • (xi»
~
,
Xj - a
Z J1 M (x.»,
f3 - a j ]
Z x. - x. J1 p (x. , x . ) ) ,
~ ) ij ~ )
Z
o - a (4.19 )
1 + (o-'r) + (f3-a) ,
and then
/L(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 Ci
i 1 2 3
6 - a .5 - .17
.24.
1 + (6-7) + (~-a) 1 + .2 + .15
.5 - .25
1 + .2 + .15 .18.
P.(x) U1 U2
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
OJ-....,..-'"-I----;...:.....-........---'I--'-r--r----r-........-...x
a 'Y{3
P.(x) U1 U3
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1 .
0 x
Note
1. The Up (x.,x.) function represents the OM's attitude toward
ij ~ J
risk. Recall that when the OM is risk-neutral, Up (x.,x.) is
ij ~ J
defined by Eq.(4.13). Subsequently, the solution Uo(i), Vi, is
computed based on Eqs.(4.19) and (4.20).
When the OM takes a risk-averse attitude, Up (x.,x.) is defined
ij ~ J
by Eq.(4.12). The solution Uo(i), Vi, is computed using
132
The ranking order for the same set of fuzzy sets may be
very high values even though very low x values are present.
of ~p (v.) (the preference index) is done using the a-cut and by trial
i 1
This does not comply with our intuition. The anti-intuition case can
i 1 2
.10 .31
i 1 2
.33 .28
+..,
d(M,N) .f I ~(x) - Il N(X) I dx (4.23)
-..,
(4.24)
:!.l...._ _ _ _ x
x·
Fig. 4.10 The Hamming distance between fuzzy sets M and N.
Yager [Y11] defines a fuzzy maximum first and then computes the
Hamming distance between each fuzzy set and the fuzzy maximum. The
fuzzy set with the smallest distance from the fuzzy max is considered
the best. Kerre [K17] follows Yager's ranking concept -- except that
he defines a fuzzy max differently from Yager. Nakamura [N2] derives
135
fuzzy minimums from the greatest lower set and the greatest upper set,
respectively. The Hamming distances for a fuzzy min from the greatest
upper set and a fuzzy min from the greatest lower set are computed. A
fuzzy order relation, p, between fuzzy sets M and N is then
constructed. The membership value ~p(M,N) denotes the weighted sum of
M's unique superiority over N for the best and worst possible
situations. Kolodziejczyk [K27] extends Orlovsky's fuzzy preference
relation [02] and the notion of Hamming distance to construct several
preference indices, P1' P 2 , and P 3 • These indices follow very closely
with Kerre's index and Nakamura's preference relation.
Notice that in Fig. 4.10 the Hamming distance is the summation
of the areas not commonly occupied by both M and N. The methods in
this category compare fuzzy sets with the fuzzy max (the fuzzy min)
using the Hamming distance and the fuzzy set that has shortest
(longest) distance to the fuzzy max (the fuzzy min) is considered
better. The ideas of comparison are good. However, since the
comparison is based solely on area measurement and the fuzzy set's
relative location on the x-axis is ignored, the logic of the methods
in this category is not sound.
For example, let us consider fuzzy sets M1 , M2 , and M3 in Fig.
4.11a. Keen observation and common sense indicate M3 > M2 > M1 .
Yager's method, which defines a unique fuzzy max for all comparison
cases, would result in: d(max,M 3 ) = 0.5, d(max,M 2 ) .452, and
d(max,M 1 ) = .436. The ranking order is M1 > M2 > M3 . This ranking
order is against human intuition. The problems with Yager's method
are that crisp numbers will always be ranked the lowest, and a fuzzy
set that overlaps more with Yager's fuzzy max will get a higher rank,
regardless of their relative locations on the x-axis.
Kerre's [K17] method defines a fuzzy max which is problem-
dependent. The fuzzy max in Fig 4.11a is M3 (by Kerre's definition).
136
P.(x) M1 M2 M3
1.
.9
.8
\
Yager's fuzzy max .'
..........
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
p.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
shows that Kerre's method would favor a fuzzy set with smaller area
137
ranked is determined. Then, each fuzzy set is compared with the fuzzy
maximum using the Hamming distance measurement. The fuzzy set(s) that
have the smallest Hamming distance to the fuzzy max are ranked as the
best.
x, (4.25)
fuzzy max for discrete functions, and Fig. 4.12b shows Yager's fuzzy
Once the fuzzy max is defined, we shall use the Hamming distance
to measure the closeness of each fuzzy set to the fuzzy max. For
example, given two fuzzy sets M and N as shown in Fig. 4.13, the
p.(x)
1.
.9
.8 Yager's fuzzy max
.7
.6
.5
.4
.3
:~o ,I
t--+-+--+--4-~-+-~'---f-~~-f. X
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
,0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0
.0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0
P.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
p.(x)
1. M •
.9 N *
.8
.7
.6
.5
.4
.3
.2
.1
0 f x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
x, 0 .1 .2 .3 .4 • 5 .6 .7 .8 .9 1
J-LM(X) 0 0 0 .3 .8 1 .8 .3 0 0 0
J-LN(y) 0 .2 .8 1 .8 .3 0 0 0 0 0
Yager's fuzzy max and M can be seen as designated by the shaded areas.
1
d(max,M) f I J-Lmax(x) - ~(x) I dx
o
.375 .375 3 .5 3
f
o
x dx - f
.3.
(~) dx + f
2
«~)
.375 .2
- x) dx
.587 .70 7
+ f «~) - x) dx + f (x - (~» dx
.50. 2 .58.2
1.0
+ f x dx .433
.7
similarly, the Hamming distance between Yager's fuzzy max and fuzzy
set N in Fig. 4.14b is d(max,N) = .50. The result is M > N.
140
JJ.(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.14a The Hamming distance between Yager's fuzzy max and M.
JJ.(x)
N M yager'\~uzzY m~~
1.
.9 " ,,,,, ,,,, .
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
JJ.(X)
N' M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~~~~--~~~~--~--~-Tx
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Numerical Example
Given three fuzzy final ratings u1 ' u2 ' and U3 as shown in Fig.
4.15, the Hamming distances between Yager's fuzzy max and u 1 ' u 2 ' and
U3 are calculated using Eq. (4.23) as:
P,(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1 -
0 ~~~~~--~~~T-~--T--T--~X
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Note
fuzzy sets with only a single element in them. For example, there are
two fuzzy sets M and N defined as:
M {(.2,.2}} N {(.2,.8}}.
P.(x)
M1 M2 M3
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
We shall use the fuzzy numbers M and N presented in Fig. 4.13 to show
the derivation of ~max(z). For z = 0.6, the possible (x,y) pairs that
satisfy 0.6 = (x y y) are
x .6 .6 .6 .6 .6 .6 .6 .5 .4 .3 .2 .1 0
Y 0 .1 .2 .3 .4 .5 .6 .6 .6 .6 .6 .6 .6
The corresponding ~(x), ~N(Y)' and ~(x) A ~N(Y) are summarized as:
~(X) .8 .8 .8 .8 .8 .8 .8 1 .8 .3 0 0 0
~N(Y) 0 .2 .8 1 .8 .3 0 0 0 0 0 0 0
~(x) A ~N(Y) 0 .2 .8 .8 .8 .3 0 0 0 0 0 0 0
We can follow the same process to get ~max(Z) for other z values.
The resulting ~max(z) is (see Fig. 4.17) :
x,y,z .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0
~max(z) 0 0 .3 .8 1 .8 .3 0 0 0
~(x) 0 0 .3 .8 1 .8 .3 0 0 0
~N(Y) .2 .8 1 .8 .3 0 0 0 0 0
144
p.(z)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~~--~--~-+--4-~--~~z
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.17 The fuzzy max for fuzzy sets M and N in Fig. 4.13.
When ~(x) and ~N(x) are piecewise linear, their fuzzy max can be
shown as illustrated by the dotted line in Fig. 4.18. The fuzzy max
x .4 [.3, .4]
y [.1,.4] .4
z =x V y .4 .4
Consequently, we obtain
~N(Y) [0,1] .5
s up ( (x) 1\ ~ (y» .5 .5
N
That is, sup (~M(.4) A ~N(Y» = 0.5, and sup (~N(.4) A ~(x» 0.5.
P,(x)
N M
1.
.9 The fuzzy max 000
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.18 The fuzzy max for continuous fuzzy sets M and N.
7-1 .45
[.10, .45]
1[.30, .45]
.45 7-1
.55
[.10,.50]
sup (IlM(x) "/olN(Y» .75 .25 sup (~(X) " Il N (y» .75
Thus, Ilmax (Z=.45) = .75 and Ilmax (Z=.55) = .75, respectively. We can
use the same procedure to get /olmax(Z) for other z values. In this
case, the fuzzy max is the fuzzy set M.
The Hamming distance between the fuzzy max and fuzzy set M,
d(max,M), is 0; while the Hamming distance between the fuzzy max
and fuzzy set N, d(max,N), is represented by the shaded area in Fig.
4.10. It is obvious that M > N.
146
P,(x)
N' M The fuzzy max 000
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 "'YI'---.---r X
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.19 The fuzzy max for continuous fuzzy sets M and N'.
When N becomes N', the fuzzy max (see Fig 4.19) is determined
.4 3 .7 8 7 .8 8
d(max,M) J (x-.) dx + J (. -x - . -x) dx + J (. -x) dx
.3. 2 .6. 4 •2 .7 .4
.05
Numerical Example
The fuzzy max for u 1 ' u 2 ' and U3 is the fuzzy final rating u 3 . Since
Note
Referring to Fig. 4.16, the fuzzy max is the fuzzy set M3 . The
Definitions
To understand the derivation of the two fuzzy min and the four
Hamming distances, the following definitions are given.
Fuzzy Min: Given two fuzzy numbers, M and N, in a universe, U,
the fuzzy min is defined as
lL(x)
N M The fuzzy min 000
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.20 The fuzzy min for continuous fuzzy sets Hand N.
determines a ~(y) value. Among the ~(y) values, their maximum gives
~-(Xo). In this case,
[.37,1.0]) 1.0.
iff
149
Consider the point XO in Fig. 4.21 again, there will be some y values
such that y :s xo. In this case Xo 0.37, Y = [0,.37] :s xo. Each y
determines a ~(y) value. Among the ~(y) values, their maximum gives
~+(Xo). In this case,
[0,.37]) .33.
p.(x)
M
1.
.9 GU of M - M+ _"''''''''
.8
.7
GLofM-M - 000
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
The fuzzy sets M- and M+ are presented in Fig. 4.21. The fuzzy
set M- may be seen as the right leg of M plus some ~-(x) = 1 for
x :s x*; while the fuzzy set M+ consists of the left leg of M plus some
distance.
150
mins and four Hamming distances are defined first. Given two fuzzy
Eqs. (4.28) and (4.29). Nakamura defines two fuzzy mins: min(M+,N+)
(see Fig. 4.22a) and min(M-,N-) (see Fig. 4.22b). The four Hamming
Fig. 4.22c):
as:
~p(M, N ) (4.30)
1
2 , if Kw = 0,
151
where
(l-W)[d(M-,min(M-,N-» + d(N-,min(M-,N-»).
,.,,(x)
1. m"1nl"1ITlnnnnl"1nl"'--'
.9
.8 d,
.7
.6
.5
.4
.3
.2
.1
o ~~--~--~~--~--~-+--~--~~x
weighted combination of the M's unique advantages over N for the best
and worst state, to the sum of such weighted combinations of the M's
J.L(x)
N M
1.
.9·
.8
.7
.6
.5
.4
.3
.2-
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
N is determined to be
of "the M and N comparison based on their best states" and "the M and
When ~p(M,N) 2: ~p(N,M), we can conclude that "M 2: N"; otherwise, "M <
resolve this difficulty, Nakamura has proved that, given two normal
K'w (4.33)
1 K' O. (4.34)
2"' w
Eq.(4.34).
Numerical Example
Three fuzzy final ratings in Example 10 are to be compared. They are
Final Ratings
o.
p.(x) U
1.
.9 <12'-
.8
.7
.6
.5
.4
.3
.2
.1
Of--~-'-+-":""""~-.----'\--':""'--.,....---..--.---,.x
54' respectively, as
a b
(area under U1 ) - ( J x - .2 dx + J
X-.17 dX
.1 .15
.2 a
.5 .5 - x
+ J dx)
b .2
.084
1.
156
As a result,
1 - .068 .932, 1, 1.
U1 U2 U3
0.5 .068 o
[ .932
1
.5
1
o
0.5
1
Notice that U3 dominates all the others because
Note
and Nakamura's method is that Yager's and Kerre's methods define fuzzy
maximums and measure the Hamming distance, but Nakamura's method uses
Even though the equations for computing the Hamming distance are
d{M,max{M,N» + d{N,max{M,N»
d{M,N) .
157
p,(x)
1.
.9
.8 M1
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
It defines the Greatest Upper set (GU) and the Greatest Lower set
(GL). Then a fuzzy max is derived from the GU and the GL instead of a
fuzzy min (as in Nakamura's method). Some preference indices are then
defined using the GL, the GU, the fuzzy max, and the Hamming distance.
where d(',') is the Hamming distance defined by Eq.(4.23), and MAN and
(4.38)
,,(x)
N M The fuzzy max 000
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
P 1 (M,N) (4.40)
160
P 1 (N,M)
the property,
1, (4.41)
into account the comparison of M and its fuzzy max at the worst and
d(M+,N+) = 5 1 + 5 2'
+ + +
d(M ,max(M ,N » = 51 + 52.
where M- and M+ are defined as before. 5ince all the Hamming distance
measurements in Eq. (4.43) have been calculated when deriving PI and
P 2 , the P 3 index is calculated as:
P3 (M,N)
Numerical Example
We have three fuzzy final ratings as U1 : (.20, .30, .50), U2
(.17, .32, .58), and U3 = (.25, .40, .70) which may be seen in
Fig. 4.1j. We will rank them using the PI index for its simplicity.
To obtain the P 1 (U 1 ,U 2 ) index, refer to Fig. 4.27. 5ince
d(U1nU 2 ,0) = 53' d(max(u 1 ,U2 ),U1 ) : 52 + 54' d(U1 ,0) = 52 + 53'
and d(U 2 ,0) = 51 + 53 + 54' we get:
.562.
P.(x) U, U2
1.
.9-
.8
.7
.6
.5
.4
.3-
.2
.1 -
0 X
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
.767, .333,
U1 U2 U3
U1 .50 .562 .767
U2
U3 [ .438
.333
.50
.309
.691
.50
1
163
Note
1. The preference index P 1 is the same as Kerre's index. Recall
that P 1 (M,N) (Eq.(4.37» is given as:
where d«MnN),O), d(M,O), and d(N,O) are constant for both P 1 (M,N) and
P1 (N,M). That is the P 1 (M,N) preference index can be determined using
d(max(M,N),M) which resembles Kerre's index. The denominator
guarantees P 1 falls in [0,1].
2. The preference indices P 2 and P 3 resemble Nakamura's index in
the following way. First of all, the property
1, i 1,2,3 (4.44)
holds for both methods. Secondly, the notions of the greatest lower
and upper sets and Hamming distance measurement are used.
The differences between these two methods may be summarized as
follows. .In Nakamura's method, P(M,N) denotes the proposition of "M ~
Nil with IIp representing the degree of truth of such an argument. But
in Kolodziejczyk's method, Pi(M,N), i 1,2,3, denotes the proposition
"M is not better than N," hence, larger Pi(M,N) indicates a stronger
degree of truth of "M :s N." Secondly, Nakamura's index defined by
Eq.(4.30) or Eq.(4.34) explicitly sets IIp = 21 where M N. That
notion is implicitly preserved by introducing Hamming distance d(',O)
such that even when M = N the denominator will not be zero and IIp,
1.
possible states. .
In Kolodziejczyk's method, equal importance of the
proposition "M :S Nil is assumed in the best and worst possible states.
Lastly, the fuzzy min is explicitly used in Nakamura's method while it
is only implicitly used in Kolodziejczyk's method.
3. One or more of the fuzzy preference relations can be used to
derive the ranking order. Kolodziejczyk has given various numerical
examples to show that there is no conflict among the results generated
by the three preference relations. If all three preference indices
always give the same ranking order, the usefulness of P 2 and P 3 (from
a practical viewpoint) is in doubt. Indeed, if one can easily apply
P 1 to derive ranking order, why bother with P 2 and P 3 at all?
165
x 123 456
.2 .5 .8 1 .7 .3
M.6 = p, 4, 5}.
P.(x)
M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
fuzzy sets with the higher a-cut values are considered better.
Similar notions can be seen in Buckley and Chanas [B40], with minor
dominance of one fuzzy set over another. For each pair of fuzzy sets,
one can have either a functional type index or a single value index.
also given.
/L(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .5 .7 .8 .9 1. x
.44 .52
(4.47)
Note
The use of the a-cut in Adamo's approach does not give reliable
ranking order. For example, referring to Fig. 4.29, when a > 0.5,
procedure. Given fuzzy sets as shown in Fig. 4.30 where not all
membership functions assume their peaks at 1, any fuzzy set that has
its peak lower than the a level will not be accepted. In this case,
p,(x)
1.
.9
MS
,---------,
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~~~~~~--~--~~-L4-~~~ x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
the proposition "M > M2 .11 Given fuzzy numbers Ml and M2 and their
\x ex
corresponding ex level sets,
(4.48)
Note
method, and should only be used for screening rather than ranking
purposes.
since a l = .40 < b2 = .56 and a 2 = .44 < bl = .52, we cannot tell
which fuzzy set is better by using this method. When ex = 0.3, we
obtain
Again, using this method, we cannot tell which fuzzy set is better.
169
.----.----.--...1-.......---4--+-.......-----.--..----. J..
-1. -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1. II
Algorithm
D = {((V'~D(v»} where
v* sup v (4.50)
veD
a
v inf v. (4.51)
veD
a
(4.52)
1, if v * v > 0,
Jij(a)
{ -1, i f v *
0, if v *
v
v
< 0,
O.
(4.53)
"M.~ ~ M.".
)
Let
~M.>M.(a). (4.54)
~ )
It follows that
a = ~M.>M.(J). (4.55)
~ )
The membership function in Eq. (4.55) shows that the statement "Mi > Mj
is at J" has a degree of truth of a.
be made.
h(D)
2 J a J(a) da, JO e [-1,1], (4.56)
o
where h(D) is the upper bound of ~D(v). This index is considered the
.--....,.....~.__-+-----.-_+___.-_.____.,.-_.___, J ..
-1. -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1. II
difference determines JOj when both 51 and 52 are on the positive side,
the summation of 51 and 52 gives JOj and when 51 and 52 are on the
negative side, -(5 1 +5 2 ) gives JO. Thus, JO can be a positive or
negative number.
K . *h
2 [ L i*J(~ ) ~*J(h)], (4.57)
K2 i=o K
achieved.
Numerical Example
Given two fuzzy sets M and N (see Fig. 4.33) ,we can compare them
P.(x)
1. M •
.9 N *
.8
.7
.6
.5 *•
.4
.3
.2
.1
0 * 10 X
0 2 3 4 5 6 7 8 9
x,y 3 4 5 6 7 8
/..LM(X) 0 1 .75 .50 .25 0
J.b(v)
1.0
0.5
~-r--~·-4--4---~~--+--4--~~V
-5 -4 -3 -2 -1 0 2 3 4 5
v -3 -2 -1 o 1 2 3 4
o .5 1 .75 .5 .5 .25 o
174
-;--1---:----:----:----:
v* sup v = 4,
veOa=o
and
v inf v = -3.
veoa=o
which indicates the degree of "M > Nil is only 0.14 at the a = ° level.
Similarly, if we set a = 0.25, a .50, a = .75, and a = 1.0, we can
a
° .25 .50 . 75 1.0
Ci
-------------------------1.0
.---.-~.___.____.-+_.J.......,.-.....____._-....___. J ..
-1. -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1. II
Fig. 4.35 a-J ij curve for fuzzy sets M and N in Fig. 4.34.
Note that J(a 1.) = -1 can also be directly obtained from Eq.(4.53)
.
s~nce v* v < o. The curve shape of a versus J ij is approximated in
Fig. 4.35. It shows that "M is strongly dominated (by N) with much
complication" (see Note 3 for details).
step 3. At different a levels, the preference orders vary.
Thus, a unified index, JO, is needed. It is computed using Eq.(4.57)
(given h = 1 and K = 4):
4
~ [E i*J({) - i*J(1)] -0.775.
(4) i=o
Note
1. Although the Jij(a) index shown previously is for two fuzzy
sets, problems with more than two fuzzy sets can also be solved using
Mabuchi's method. If there are U1 , ... , Um fuzzy final ratings, we
can use Jij(a) for each pair of (Ui,U j ) and construct a preference
relation matrix. The ranking order can then be obtained.
Or we can define a fuzzy min as
176
(4.58)
and compare each ~u{ (xi)' Vi, with ~min(v). For each ~U (Xi)' Vi, we
... i
get a Jmi(a). The final ratings with larger Jmi(a) indicate higher
ranking position.
There are two fuzzy sets, Mi and Mj , (as seen in Fig. 4.36). For
M{~
... ~
= [a.,b.] and M.
J
~ ~ a
There are different equations for calculating Jij(a) under differnt
situations.
Ib i -
ajl - la i - b·1
J
Jij(a) (b i - a.)
J
- (a.
~
-
b. )
J
(4.59 )
1, bi > a j ,
Jij(a)
{ 0, b.
~
aj ,
-1, b i < a ..
J
(4.60)
(a i + bi) - (a j + b j )
J .. (a) (b. (4.61)
~J
~
a i ) + (b j aj )
177
JL(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
o ~~--~~--~~~~-+--~~~ x
o aja j b j b j 1.
1
Given a difference set, D, that has a continuous membership
(4.65 )
Note that the results from both Eqs.(4.61) and (4.65) are the
same. This verifies that the comparison ofM. and M. can indeed be
l. )
over the J ij axis, (b) the average position along the J ij axis, and
(c) the height of the curve. Table 4.2 gives the correspondence
Nearly 1 Credible
Over 0.5 Moderate credibility
Under 0.5 Little credibility
Nearly 0 No credibility
180
ex curve 1
functions are used to derive IT (M<N), IT(M=N), and IT(M>N) which denote
the possibility of "M < N," "M = N," and "M > N," respectively. The
three indices are compared, and the index with the highest value
functions. Such functions are then used to define two ranking indices
for fuzzy sets M and N. The first one measures the proposition "M ~
Nil while the second one measures the proposition "M :s N." When the
first index value is greater than zero, the second index value must be
zero, and vice versa. The non-zero index determines the relationship
between M and N.
Definitions
Possibility Measure
from the power set of U, 2U, to the unit interval [0,1]. That is,
184
U
Again, let U = {2,5,8}. 2 = {0, {2}, {5}, {8}, {2,5}, {2,8}, {5,8},
~F(x) as:
A = {2,5} is:
(4.71)
To show Eqs. (4.70) and (4.71) give the same possibility measure,
F = {( 2 , .4), (S,. 7), (8, 1)} and A = {( 2, • S) , (5, . 6)} are used again.
<Xl
.1 .2 .3 .4 .5 .6 .7 .8 .9 1.
.1 1 1 1 1 1 1 1 0 0 0
.2 1 1 1 1 1 1 1 0 0 0
.3 1 1 1 1 1 1 1 0 0 0
.4 1 1 1 1 1 1 1 0 0 0
<X2 .5 1 1 1 1 1 1 1 0 0 0
.6 1 1 1 1 1 1 1 0 0 0
.7 0 0 0 0 0 0 0 0 0 0
.8 0 0 0 0 0 0 0 0 0 0
.9 0 0 0 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0 0 0 0
Given the same F and A, Eqs.(4.70) and (4.71) will result in the same
possibility measure.
Necessity Measure
N: 2 U --+ [0,1],
and
and II (A) .
(4.85)
having M as a fuzzy bound: (i) (M,+co), (ii) (-co,M], (iii) ]M,co), and
189
(iv) (-oo,M[.
(i) The set of numbers possibly greater than or equal to M is
[M,+oo) (see Fig. 4.40a) with membership function
/\
o
P,(x) j
The fuzzy set (-oo,M[
1. "\T"o"'o"'V"'(F~" """""':' -, - - - - - - - - -
involves two fuzzy numbers M and R. The fuzzy sets [R,+~) (the set of
numbers possibly greater than R) and ]R,~) (the set of numbers
necessarily greater than R) are used to derive four indices which
access the proposition "M ~ R." The four indices are IIM([R,+ao»,
IIM(]R,+ao», NM(]R,+ao) and NM([R,+ao» (see Fig.4.41).
JL(X) M R
(iii) By applying Eqs. (4.77) and (4.86), we get the third index,
(iv) By applying Eqs. (4.77) and (4.88) we get the fourth index,
proposed indices. Table 4.3 clearly indicates that the four indices
tions shown in Table 4.3. Thus, the four indices must work in a group
IMI
I
I
R I 1 1 1 1
---I: 1 1 1 1 0
! I
--r--r-i
1 1 1 0 0
1 0 1 0
1 :I--! 1 0 0 0
'---I
I 0 0 0 0
I
shall extend the pairwise indices (as given by Eqs. (4.90), (4.91),
v=max(u 1 ,···,um)
This fuzzy set is very easily obtained from the shapes of M1 , ... , Mm.
That is,
n
max U [( (\ [M. ,+00» (\ M.], (4.95 )
i=l j .. i J 1.
Numerical Example
Dubois and Prade's ranking indices will be applied to Example 10.
First, we shall derive the fuzzy maximum of U1 and U2 using
Eq. (4.94) as:
/L(X) U1 max 23
1. ----------- ----- - -- ----------1-----------------
I'
--.~-
.9 \, ~
~
I.'.'
.8
.7 .'.'.'
,.'
~
.6 I
.5 .'
.4
.3 NO(U)
.2 . __ ...1.. ... __ .................... ~~~~~1
.1
o ~~--~~--~~--~~--~-T~X
o .1 .2 .3 .4 .5 .6 .7 1.
Similarly, we can get the four indices for U2 and the four
indices for u 3 . The results are summarized in Table 4.4. All four
Note
1. Dubois and Prade [D29] gave physical meaning to the four
indices as follows.
leg of R) with the best part of M (the right leg of M). The index
value will be high when the left leg of R is smaller than or equal to
the right leg of M.
index will be large when the right leg of R is smaller than the right
leg of M.
index will be large if the left leg of R is smaller than the left leg
of M.
(iv) NM(]N,oo» compares the right leg of R with the left leg of
2. The equations for deriving PD, PSD, ND, and NSD are very
197
L (m~u), u:!5 m
JIM(u) {
R (u;m), u ~ m
L (n~u), U:!5 n
{
R (u~n), u ~ n
,.,,(x)
M N
1.
a· m- a
~. b - m
1· n - c
o· d - n
The ranking problem has come down to solving the following equations.
u"-m u"-n
1 - R(-{3-) R(-c')-)
and Kwakernaak's [B1] ranking index, while the NSD index is identical
to watson et al.'s [W3] ranking index. Dubois and Prade [029] pointed
out that unless all four index values are used, counter-intuitive
extent that "H ~ R." The IIR[H,m) index measures the extent that "R ~
method.
199
P,(X) M
1. -------
M' R
indices are IT (N<M), IT(N=N), and IT(N)M). They denote the possibilities
of the three events liN < M," liN M," and liN > M," respectively, given
that M is the reference set. The possibility measure with the highest
Notice that if Dubois and Prade's [029] notations are used, ML = [M,m)
The three indices on Nand M are then defined (see Fig. 4.46) as:
/L(x) 1 - to;,. M N 1- ~
1. +----.;;;....,
JI(N"ML_
-mlil>Mr---
ll(N<M)
p.(x) U1 1 - U1R
1. -!-__.,----
,
,,,,,
· , r ; , - - - - - - - - - l ' f f I I f J l l l l l l l f l " " ' l I f l f l l r t . I I ! I I I I I f { l f f l l l l l f I f I l I lJ
.9
.8
.7 '-. f
.6
.5
.4
.3
.2
.1 1- I1L
0 ~~=W~~~~~~~~~~~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Numerical Example
Similarly, we obtain
Note
and
N(N)M) inf max (1 - ~ (u) , 1 - IlN(U». (4.107)
u
Definitions
The Fuzzy Sets "numbers possibly ~ M" and "numbers possibly s M"
A fuzzy number, M, defined by 4-tuple (a,b,c,d) is shown in Fig.
4.48. The membership function is defined (Dubois and Prade [026]) as:
S (b-x I b-a), x s b,
~(X) { S (x-c I d-c), x ~ c, (4.111)
1, otherwise.
204
#L(X)
1. -............... ,.------.
greater than M, f M, and the set of numbers possibly less than M, gM'
0, if x a,
{
:5
1, if x b,
{
:5
gM(x) 0, if x ?! d, (4.113)
/.lM(X) , i f d > x > b.
property "?! M." Other possibility distributions such as f M' and fM"
tive matter.
205
JL(X)
1.
Fig. 4.48a The possible distribution functions: fM' fM" and f M".
JL(X)
1.' """".,","
Fig. 4.48b The possible distribution functions: gM' gM', and gM"·
measure.
0, if x s a,
~(x)
{ 1, if x ~ c, (4.115)
~(x), if c > x > a.
(4.116)
\I x e R, f M , A (x) (4.117)
or
1, if x ~ d,
{ 0, if x s c,
1 - gM(x), if c > x > d.
(4.119)
207
which has the same properties as those of ~(X). Let nACo) be the
A-measure corresponding to PM(X). The associated comparison function
is
(4.120)
Y1 Y2 Y3 zl z2 z3 z4
xl .3 .7 .1 Y1 .9 0 .4 .4
R1 x2 .6 .5 .2 R2 Y2 .4 .2 .3 0
x3 .8 1 .1 Y3 0 1 .7 .6
zl z2 z3 Z4
xl .4 .2 .3 .3
x2 .6 .2 .4 .4
x3 .8 .2 .7 .4
max { . 4 , . 3 , . 2 } .4
Max-min Transitive
all x E X.
pairs.
identified).
where ~ is the set of all real numbers. Two properties are observed
The property (b) suggests that a T (·,·) is not a strict fuzzy order
relation, and it is not adequate to serve as a comparison relation.
max-min transitive (the reader should refer to the original paper for
(4.125)
know "N ~ M" is true. When 0T > 0, we know "N :5 M" is true.
Algorithm
where
0, if x a2 ,
{
:5
1, if x ~ c2, (4.127 )
~N(x), if c 2 > x > a 2 ·
(4.128)
211
where
0, if x s C2 '
{ 1, if x
1 -
~
~N(X),
d2,
if d 2 > x > c 2 •
(4.129)
If i3 T = 0 then cST > 0, and vice versa. When i3 T > 0, we know liN ~ Mil
is true. When cST > 0 we know liN s Mil is true.
The algorithm must be repeated for every pair of fuzzy numbers in
question. It should be noted that the comparison relations are more
"crisp" (i.e., i3 T !!E 1) when A' is the smallest, and cST !!E 1 when A" is
the greatest.
Numerical Example
Delgado et al.'s method will be applied to Example 10 in which
the final fuzzy ratings are: u 1 = (.20,.30,.50), u 2 = (.17,.32,.58),
and u 3 = (.25, .40, .70).
We shall give the detailed computational procedure on the
comparison of u 2 and u 3 to demonstrate the algorithm.
step 1. Set A' = 0 and A" = 0 as the requirement levels for the
accomplishment degree of the property ,,~ u 3 " and "s u 3 ," respectively.
Step 2. Compute comparison functions, fu A'(x), and gu A"(X),
3' 3'
respectively. Given the membership function of u 3 '
212
P.lx)
1.
.9
.8-
.7
.6
.5-
.4
.3
.2
.1
0 ~~~~~~--~~~--~~~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
x - .25 ,
.40 .25 - .25 :s x :s .40,
0, otherwise,
0, ifx :s .25,
hu (x)
3 { 1, ifx
Jl u
'"
(x), if x e (.25,.40),
.40,
3
and
1, if x '" .70,
Pu (x)
3 { 0, ifx :s .40,
1 - Jl u (x), if x e (.40,.70).
3
or
0, if x .70
'"
gu
3'
A" (x)
{ 1, ifx
IlU (x),
;!;
if x
.40
E (.40,.70)
3
.9 j
~
,8 I
I
,7 I
~
,6
,I'
I'
.5
.4 ,/'
.3
,"
~
.2
.1
f
~'
o """"'''''''''''''''''''''' x
o ,1 ,2 ,3 .4 .5 .6 .7 ,8 .9 1.
step 3. The comparison relation for the "", u 3 " case is computed
as:
1 - .82 = .18,
where
214
~u (x) .58 - X
, .32 < X ::s .58,
2 .26
0, otherwise.
Similarly, we can compute the comparison relation for the case of "::s
Fig. 4.51).
.7
.6
.5
.4-
.3
.2
.1
0 ~~~=L~~~-i~~~~~~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
o o o
Note
and
be generalized as:
computation efforts:
found, ~T(M,N) = o.
c. The computed ~T value tends to be close to 1 when A' is set
to its minimum, e.g., A' = -.999. However, the ranking order will not
The method in this category compares the fuzzy numbers using two
criteria: the mean value and the spread of a fuzzy number. The
number with a higher mean value and a lower spread, i.e., higher
Lee and Li [L3] propose the use of generalized mean and standard
deviation based on the probability measures of fuzzy events to rank
J
SCM) x ~M(x) dx
(4.137)
The denominator measures the area under fuzzy number M (see Fig.
CJ'U(M) (4.138)
where SCM) is the support of fuzzy number M (see Section 3.2 Basics of
Fuzzy Sets).
218
P.(x)
M
1.
.9
.8-
.7
.6
.5
.4
.3
.2
.1
0
P.(x) M
1.
.9
.8
.7-
.6
.5
.4
.3
.2
.1
o !--........-.-~-..---i-m-..----.--+n-.....-----.. x
1
'3 (1 + m + n), (4.139)
(]"U(M) (4.140)
S x 2 JlM(X) dx
SCM)
Xp(M) (4.141)
S [~(X)]2 dx
SCM)
S x 2 [JlM(x) ] 2dx
1/2
SCM)
- [x"
(jp(M)
[ S
SCM)
[ ~(X)]2dX
(M)) ,] (4.142)
rewritten as:
i (1 + 2m + n) (4.141a)
(jp(M) (4.142a)
distributions, respectively.
case.
220
Relation of Relation of
X(M i ) and X(M j ) (T (M i ) and (T(M j ) Ranking order
Notice that the spreads, CT, are used only when the mean values are
equal; a smaller spread indicates a better fuzzy number.
Numerical Example
There are three investment projects to be evaluated according to
attributes X1 (cost) , X2 (environmental impact), X3 (estimated profit),
and X4 (maintenance cost/yr). The fuzzy final ratings are given
(see Fig. 4.54) as:
i 1 2 3
#,(X)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
o X
o 2 3 4 5 6 7 8 9 10
since X(U 1 ) > X(U 2 ) > X(U 3 ), the ranking order is U1 > U2 > U3 .
of U1 as:
1
"31 (l+m+n) "3 (S + 6 + 8.4) 6.47,
.S1.
Note
1. This method ranks fuzzy numbers based on two different
criteria, namely, the fuzzy mean and the fuzzy spread of the fuzzy
and (4.138) can sometimes be conflicting with those from Eqs. (4.141)
and (4.142). Lee and Li suggest that when conflict does occur, the
Lee and Li [L3] give the following example. Referring to Fig. 4.SSa,
222
X(M) (4.143)
and
1 b4 ab 4 a4
<reM) { [b-a (4 -3- + 12) + !3 (c 3 _b 3 ) + (4.144)
1 d4 c 3d c4
d-c (12 -3- + 4 )]/[~ (-a-b+c+d)]} -
{(-a 2 -b 2 +c 2 +d 2 -ab+cd)/[3(-a-b+c+d)]}2,
"(x)
1.
"Ix)
1.
-I
p.(x) M,
1.
P.(x) Mj
1. The f uzzy max 000
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Given two fuzzy numbers Mi and Mj (as seen in Fig. 4.56), the
contribution of fuzzy number Mi to the fuzzy max is determined by:
226
(4.145)
where
P.lx) Mj
~
1. The fuzzy min 000
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
where
~min(x)
227
For example, If P(M i ) > P(M j ) and N(M i ) < N(M j ), then Mi > Mj .
In case of unknown ranking order, McCahone proposed two composite
indices to discriminate Mi and Mj . The two composite indices are
defined as:
(4.147)
•
and
(4.148)
The equations and comparison rule sets are good for ranking two
fuzzy numbers. When more than two fuzzy numbers are compared, the
following steps must be taken:
step 1. Find the fuzzy max and the fuzzy min based on all fuzzy
numbers in question (see Fig. 4.58).
step 2. Compute P(M i ) and N(M i ), vi.
step 2.1. Rank the fuzzy numbers according to P(M i ) alone (the
bigger the better).
step 2.2. Rank the fuzzy numbers according to N(M i ) alone (the
smaller the better).
step 3. Compare the two ranking orders found in steps 2.1 and
2.2.
229
#L(X) M1 M3 M2
1. o -,
'. . .,
, ,
.9 ~ ~;. *...
. . ."
"
,
.8 "
'
.7
.6 .".'
.5
.4
.'"
"
.3
.2
..
'
.1
0 ~~~~--~-+--~~~~--~--~~
.'. x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.58 The fuzzy max and the fuzzy min for more than
two fuzzy numbers.
pick those fuzzy numbers that contribute to the conflict and perform
algorithm stops.
Numerical Example
investing a large sum of capital. The tools are commodity, stock, and
1,2,3, as follows:
i 1 2 3
.670 1 .560
230
.9
.8
u
.7
.6
.5
.4
.3
.2
.1
o ~~--~--~-+--~--~~--~~~~x
o 2 3 4 5 6 7 8 9 10
9
S ~max(x) A ~U (x) )dx
3 1
9
.61
S ~U (x) ) dx
3 1
9
S [~min(x) A ~U (x)) dx
3 1
.67.
9
S [~U (x)) dx
3 1
Based on the P(U i ) index alone, U3 = U2 > U1i while on the N(U i )
index alone, we have U3 > U1 > U2 . The two ranking orders are
further comparisons.
.48
and
.52.
Similarly, CP(U 2 ) = .47 and CN(U2 ) = .53. Since CP(U 1 ) > CP(U 2 ) and
have completed the algorithm and have come to the conclusion that
U3 > U1 > U2 ·
Note
ideals, (the fuzzy max and the fuzzy min), as the basis for comparison
that are not overlapping with the fuzzy ideals, as seen in measuring
the Hamming distance. A fuzzy number that overlaps more with the
fuzzy max and less with the fuzzy min is considered better in
McCahone's method.
min and Mn is the fuzzy max, and Ml and Mn do not overlap with other
method.
.9 0
o
.8
.7
.6
.5
.4
.3
.2
.1
oo .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
x
Given fuzzy numbers M1 , ... , Mn' the left and right scores refer
to the intersections of a fuzzy number Mi with the fuzzy min and the
notions (using Chen's [C12] fuzzy max and fuzzy min) where R i , i 1,
methods in this category utilize either the right score or both scores
to derive the total score for each fuzzy number. The fuzzy number
\
IL(X)Chen's fuzzy min Chen's fuzzy max
1, :',
' ,
/
,,
.,
,,
Fig. 4.61 Illustration of the left score and the right score.
where x is the support of Jain's fuzzy max and x max is the biggest
a fuzzy number and Jain's fuzzy max is the crisp score of the fuzzy
number. The fuzzy number with the highest score is the best choice.
234
rank them equally, even though one is obviously better than the other.
except that a fuzzy min is used in addition to a fuzzy max. Note also
that the fuzzy max defined in Chen's method is slightly different from
that in Jain's method.
Chen and Hwang (C13] pointed out that Chen's method considered
answers. Thus, they redefined a fuzzy max and a fuzzy min (which are
procedure for obtaining the left and right scores is the same as that
fuzzy number a crisp score. The fuzzy numbers with higher scores are
section of the nonincreasing part of a fuzzy number with the fuzzy max.
should consider both the fuzzy max associated with various fuzzy
k
i-Lmax(X) (x / X max ] , k > 0, (4.150)
235
n
sup [ V S (M . ) ] (4.151)
i=l 1
and
P.(x)
Jain's fuzzy max
1.
1'0(1)
00:::0(0::)0:::::::::::::: : ::::0 ;
1'0 2 ./
Numerical Example
i 1 2 3
.56 .60 .70
leg of ~U1 (x) and ~max(x) to illustrate this point. Bo.th U1 and
.50 - x • 30:5 :5 . 50 .
~U (x) . 20
X
1
x
~max(x) = 7i
P,(x)
Jain's fuzzy max
1.
.9
.8
0(3)
.7 --~(2r-------
.6
.5 0(1)
.4
.3
.2-
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Note
and Degani [B34], Chen [C12], Zimmermann [Z31], etc.) that this method
and M2 (see Fig. 4.64) which share the same "right leg," Jain's method
contain negative support (i.e., for some x < 0 in {xl gM. (x) > O}),
~
IL(X)
Jain's fuzzy max
1. , .. :
.9
.8
.7
.6
.5
.4
.3
.2
.1 .
0
0 .1 .2 .3 .4 .5 .6 .7
should include the minimum and maximum x points associated with the
x max and xmin are the maximum and minimum numbers in the support set
n
sup [V S(M i )], (4.155 )
i=l
239
n
inf [V SeMi)]' (4.156)
i=l
and
{x I ~,(x) > O}. (4.157)
1.
x - x
max ]k, (4.158)
~min(x)
xmin - x max xm1.'n s X ::s xmax '
where xmax ' xmin' and k are defined as those in Eq.(4.154). Chen's
fuzzy max and fuzzy min are presented in Fig. 4.65.
1, !'-"/ ,;
, I+R(2)
"-------------:------
min
There are two scores associated with each fuzzy number Mi'
namely, the right score ~R(i) and the left score ~L(i) (see Fig.4.65).
They are defined as:
(4.161)
Numerical Example
Chen's method will be applied to Example 10. Assuming k 1,
Chen's fuzzy max and fuzzy min are defined as:
x - x .
m~n x - .17
tlmax(X) .53 ], .17:s x:s .7,
x - .7
tlmin(x) x. - xmax --~--.'5~3~--]' .17 :s x :s .7.
m~n
x - .2 .2 :s x :s .3 (left leg),
{
.1
J.l U (x) .5 - x
1 • 3 < x :s .5 (right leg) .
.2
The intersection of the left leg and J.lmin(X) is given as J.l L (l) = .79
The intersection of the right leg and J.lmax(X) is J.lR (l) = .45. The
summarized as:
i 1 2 3
J.lR{i) .45 .52 .64
Note
1. It is noted by Chen [C12] that the J.lT(i) index alone may not
numbers assume the same total score, Chen's method cannot discriminate
point such that J.lM. (x') = 1. But in this method, Chen removes that
~
242
L,:
Chen's fuzzy max
M1~
JL(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
o ~~ __ ~-L~~__~__~~~~__~~ X
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.66 Fuzzy numbers with equal total scores by Chen's method.
JL(x)
M1
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
restriction by assuming ~M. (x) e [O,h] with 0 < h s 1. The reason for
~
and M2 (as shown in Fig. 4.67), where M1 is a normal fuzzy number and
M2 is not. M1 and M2 have the same meaning, eg., about 0.4. These
evaluating the same character. The person who expresses M1 has more
M2 does.
243
JL(X) Mn
1. -_. -- ;~- ----------- --------------------------------
x X
max
Given fuzzy numbers Mi , i = 1, ... ,n, each has three vertices with
coordinates (ai,O), (mi,h i ), (bi,O) as shown in Fig. 4.68. The
membership functions of Mi , Chen's fuzzy max and fuzzy min are defined
(set k = 1) as:
!
hi{x - ai)/{m i - a i ), ai ::s x ::s mi ,
(x) (4.162)
~i hi{x - bi)/{m i - b i ) , mi ::s x ::s bi ,
0 otherwise,
h (x - x min )
(xmax - x min ) (4.163)
~max{X)
o , otherwise,
244
h (x - xmax)
(xmin - x max ) (4.164)
/.Lmin(x)
o , otherwise.
The left, right, and total scores of fuzzy number Mi , Vi, are then
given as:
h h.
1
(b i - x min )
/.LR (i) (4.165)
hi (Xmax - Xmin ) - h (m. - b i )
1
h h.
1
( x max - ai )
/.LL (i) (4.166)
hi (Xmax - x min ) + h (m i - ai )
(4.167)
/.LT(i), however, are much more complex Since they don't provide
considering the case of Fig. 4.69, both M1 and M2 not only have
the same crisp scores but also the same m = .5 value such that
245
./!
Chen's fuzzy min Chen's fuzzy max
:'" \ \ M, M2
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~--4---~~--~--+---r-~--~ X
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
means.
Another difficulty may arise when applying Chen's method. That
The equality is attributed to the fact that this method considers only
P.(x) M2
M1
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
p.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--,---~~--~--~-+--,-~~~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Chen and Hwang [C13] found that Jain's method may result in
".(x) M1
1 Chen and Hwang's ~Chen and Hwang'
. ""fuzzy min fuzzy max . /
..... .)
L"/~R(2)
<. - •••••••• -.iLR(1i .
.5
~~--~--~~--~--r_~~~~~~ X
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.71 The left and right scores by Chen and Hwang's method.
J.Lmax(x) { x, o :s X :s 1,
0, otherwise,
(4.168)
{
1 - x, O:sx:s 1,
J.Lmin (x) 0, otherwise.
(4.169)
Chen and Hwang's fuzzy max is the same as Yager's fuzzy max defined in
Section 4.3.1. Their fuzzy max and fuzzy min are defined in a manner
such that absolute locations of fuzzy numbers can be automatically
incorporated in the comparison process. Note that their fuzzy max and
fuzzy min are used for all fuzzy comparison cases.
The computational procedure is the same as that of Chen's except
that the definition of J.Lmax(x) and J.Lmin(X) is different. The left
utility score of each fuzzy number Mi is defined (see Fig. 4.71) as:
(4.172)
Numerical Example
We shall apply Chen and Hwang's method to Example 10. The left
and right scores are presented in Fig. 4.72. The total score of
fuzzy number U1 is calculated as follows:
P.(x)
1. ' .
.9
.8 .73
.7
.6 , .................. _.. ,9.4~
.5
.4
.3
.2
.1 .' ••.......
o ~'-'--.--L.,"--'---r----.--+--:...-+---.--....-~ X
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig 4.72 The left and right scores for the fuzzy final ratings.
249
i 1 2 3
J,LR (i) .42 .46 .54
Note
1. We shall compare the results of the three methods. Given
three fuzzy numbers M1 , M2 , and M3 as shown in Fig.4.73a, the ranking
order should be M3 > M2 > M1 , based on human intuition. Table 4.5a
summarizes the crisp scores obtained by the three approaches.
Ml .J75 .400
M2 .415 .440
MJ .786 .777
P.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~-+--T-~--~~~~--+--T--~-TX
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Fig. 4.7Ja comparison among the fuzzy score methods - part (a).
251
~~
1.
.9
~
~
~
.5
A
~
.2
.1
0 x
0 .1 .2 ~ .4 .5 .6 ~ .8 ~ 1.
spreads will be rated equivalent (see Fig. 4.69). In other words, the
Hwang's method.
252
P.(x)
M
1.
Yo ----------------- ---------,
1 1
J g(x) ~. (x) dx / J ~M. (x) dx, (4.173)
o ~ 0 ~
whose value is equal to the area under the membership function ~.'
1
geometric center (see Fig. 4.74). The value of xo may be seen as the
weighted mean value of fuzzy number Mi' The fuzzy numbers with higher
xo values are considered better.
Numerical Example
x - .2 0.2 x
,; ,; 0.3,
{
.1
lJ. u (x) =
.5 - x
1 0.3 < x ,; 0.5.
.2
1 1
f x IJ. U (x) dx / f lJ. u (x) dx
o 1 o 1
.30 2 .50 5
f (x * ~) dx + f (x * ~) dx
.20 .1 .30 .2
0.3 .50 5
f (x-.2) dx + f (~) dx
0.2 .1 .30 .2
.013 + .037
.333 .
. 150
Note
It is not surprising to see that the Xo index alone provides very poor
P.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
00 2 3 4 5 6 7 8 9 10 x
P.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~--~--~~~4-~r-~--~--r-~x
0 2 3 4 5 6 7 8 9 10
and both have the same fuzzy mean but different fuzzy spreads. Since
smaller spread should be ranked higher. If one follows Lee and Li's
1 1
J x ~M(x) dx / J ~M(X) dx, (4.174)
o 0
1 1
J x ~(X) d~(X) / J ~(X) dx. (4.175)
o 0
indices, the fuzzy number ~ is said to be optimal if x ok= max [Xoi ]'
i
256
/L(X) M
1........... --... --_. _......... -
/L(X)
1... _..... __ ._-_._-_._- --- M
There is not always a unique optimal choice. For example, Mk may have
a maximum xoi value but Ml has a maximum Yoi value. Thus the choice
of Mk or Ml is not clearly defined. Murakami et al. suggested that
the OM may make a decision based on his subjective judgment of which
Numerical Example
Murakami et al.'s method is applied to Example 10.
{
.1
Jl u (x)
.5 - x 0.3 < x ::s 0.5.
1
.2
1 1
x o1 ~ x JlU1 (X) dx / ~ JlU1 (x) dx
.30 .50
f (x X-.2) dx + f (x .5-X) dx
.1 .2
.20 .30
0.3 2 .50 5
f (x-.) dx + f (. -x) dx
0.2 .1 .30 .2
.013 + .037
.333 •
. 150
1
Since U1 is a fuzzy triangular number (see Note 1), Y01 = 3·
be summarized as:
i 1 2 3
.33 .37 .45
since all Yoi' i = 1, 2, 3 are the same, the ranking order depends
(4.176)
F(M) f (4.177)
o
where a max = sup ~M(x) and X(Ma) represents the average value of the
x
elements having at least a degree of membership .
/L(x)
1.
p.(x)
1.
.9-
.8
.7
.6
.5
.4
.3
1
.2
.1
0 1
~-+--4---~-+--4---~-T--~--~-T
.3 .5 .6 .7 .8 .9 1.
X
0 .1 .2 .4
.2 .4 1
F(M) I .35 da + I .45 da + ..• + I .60 da
o .2 .8
x - 3, 3 :S X :S 4,
IlM(x)
{ 1
9 - x
-3-
4 < X
6 < X
:S
:S
6,
9.
X(Ma) [a+3+9-3a] / 2 -a + 6
1
F(M) I X(Ma) da 5.5.
o
The shaded area between X(Ma) and the vertical axis in Fig. 4.79
denotes F(M). The higher F(M) index value denotes a better fuzzy
number M.
262
P.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 -......:..,-'-~-'>------.. X
0 2 3 4 5 8 9 10
Numerical Example
F ranking index.
x - .2 x
.2 .3,
{
~ ~
.1
IlU (x)
1
.5 - x
.3 < x ~ .5.
.2
Hence,
1
J (.35 - .05a) da .325.
o
summarized as:
i 1 2 3
.325 .350 .440
263
since the higher F index value indicates a better fuzzy number, the
Note
1 Yager [Y18] points out that when M is a normal, trapezoidal
F(M) (4.178)
x values having a = 1.
2. The counter-intuitive example of the F(M) index is demonstrated
same F index value. Yager [Y18] interprets the three fuzzy numbers as
suggests that since all three fuzzy numbers represent the same fuzzy
This argument is not shared by Lee and Li [L3] who favored the
assumption that human intuition would choose a fuzzy number which has
a higher mean and a smaller spread. In this case, the ranking order
264
/L(x)
1. --- ----------- --------.-- ---
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 2 3 4 5 6 7 8 9 10
Tong and Bonissone [T10] explore the second idea of Freeling, and
various decision statements that take the format of "M is P over all
example, given fuzzy numbers M and N, when the decision maker feels
preference."
transitive if
final ratings.
The questions are asked surrounding the following five factors:
Algorithm
Given fuzzy numbers M1 , ... ,Mm, the major steps of this algorithm
are:
step 1. Define a set of terms that contains the various degrees
of preference the OM requires, such as {strongly, slightly, very much,
no preference, •.• }. This term set is problem- as well as
OM-dependent.
M1 Mi Mm
M1 r 11 r 1i r 1m
L Mi r i1 r ii rim
Mm r m1 rmi r mm
~(xl
1.0 ,...,
NO SL SO VE
-4 -2 0 2 4 6 8 10 12 14
be:
transitive if
r ij l!: (ror) ij
where
Numerical Example
An example from Efstathiou and Tong [ES] is used to illustrate
some of the features of this approach. Three fuzzy final ratings for
alternative A1 , A2 , and A3 are given linguistically U1 'Slightly
below medium,' U2 = 'Poor,' and U3 = 'Good but not very good' (see
Fig. 4.82). Assume the decision maker's preference is obtained as:
U1 U2 U3
L NO SL NO
NO NO NO
SO VE NO
1
where NO, SL, SO, and VE represent 'no preference,' 'slightly,'
'somewhat,' and 'very much,' respectively. Fig. 4.81 shows a possible
way to interpret the above linguistic terms.
similarly, the peaks are well separated between Ul and U3 . U3 has the
best peak position. When we analyze the trade-off between the good
outcome vs. the bad outcome, we see that Ul has a low possibility of
Assume the decision maker has decided that the alternative must at
here since the peaks are well separated. However, the lower peak
value of U3 makes its preference over the others less strong. Thus,
L(U 3 ,U l ) = 'somewhat.'
1
[
NO+NO SL+NO
NO+NO
LoL NO+NO NO+NO NO+NO
SO+NO SO+SL NO+NO
adding "NO" to any term will not alter the meaning of that term.
273
[
NO SL NO
NO NO NO
SO SL+SO NO 1
All the elements in LoL equal their corresponding elements in L except
ranking order is
Somewhat> Slightly>
~(x) SO+SL VE
1.0
-4 -2 0 2 4 6 8 10 12 14
Note
Tong and Bonissone [T10] assume that the final ratings for m
ranking order can be derived from the dominance set. This dominance
set alone may result in several fuzzy numbers being ranked the same.
That is, we may see some indiscrimination case when using the
language.
Dominance Set
{
1, if x s xo,
J..L s U. (x) o (4.182)
1
J..L U . (x), if x > x ,
1
with Xo being the lowest value of x for which J..L u . (x) 1; and liS U."
1
1
275
is the notion for "less than or equal to Ui " formed from Ui • In Fig.
4.84a where the peak of Ui is to the left of Uj' Eq.(4.181) will
result in the intersection of ":s U." and U.; while in Fig. 4. 84b,
~ J
Eq.(4.181) will always give a value of 1.
~
3
~
.5
A
~
~
.1
0
0 2 3 4 5 6 7 8 9 10
U1 U2 U3 U4
~D U1 1 .50 0 0
[
U2 1 1 .50 .80
U3 1 1 1 1
U4 1 1 1 1
1
Eq.(4.181) cannot tell if U3 is better than U4 or the other way
(4.183)
277
is calculated as:
This indicates that the overall degree to which U2 dominates all other
Xk - [ (4.184)
n
~Z (z) max [ min ~U. (xi)]' (4.185)
k (x 1 , . . . ,x n ) i=l 1
alternatives.
One can easily obtain Z3 and Z4 using Eq. (4.185) if ~U (x.), ~i,
i 1
278
are discrete membership functions. When "u. (xi) are piecewise linear
1.
(4.186)
(4.187)
P.(x)
1.' ........ . M
n
L V(U i ) Ui
i=1
iook , k 3,4., i 1,2,3,4, (4.190)
Zk Uk - n
L V(U i )
i=1
iook
(9,9,2,2) - 0*(1,1,.5,2)+.5*(3,6,2,4)+1*(7,10,1,.5)
o + .5 + 1
(9,9,2,2) - (8.5,13,2,2.5)/1.5
(.33,3.33,3.67,3.33),
P.(z)
1
-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6
a term of the linguistic variable "true" which takes its value in the
term set = {true, not true, very true, ... false, not false, very
false, ... , not very true and not very false, ... } (Zadeh [Z5]).
The next step is to find terms in the term sets that approximate
Algorithm
The algorithm can now be summarized by the following steps.
min
j
281
Numerical Example
Three investment tools, commodity (A 1 ), stock (A 2 ), and real
estate (A 3 ) are to be evaluated according to four attributes: risk of
losing capital (Xl)' inflation impact (X 2 ), interest received (X 3 ),
and cash availability (X 4 ). The fuzzy final ratings are given as:
U1 = (1.26, 1.34, .62, .40), U2 = (1.46, 1.46, .86, .80), and U3 =
(2.32, 2.42, .94, .58) (see Fig. 4.88).
Given the fuzzy final ratings, we shall use the stated algorithm
to find the best solution and corresponding linguistic expression.
step 1. Obtain the dominance matrix using
as
U2 U3
[
o .90 .27
1 .51
1 1 1
The element ~D(U1,U2) is obtained as follows. The fuzzy set ~<U (x)
- 1
is obtained using Eq. (4.182) as (see Fig. 4.89):
~<U
- 1
(x)
{ 1, if
~U
1
x s 1.26
(x), if x > 1.26
282
p.(x) U1
1.
LLa
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 2 3
P.(x)
1.0000000000
.9
.8
.7
.6
.5
.4
.3
.2
.1
o ~~~~~--~--~~~~~--~~x
o 2 3
The dashed line in Fig. 4.89 represents the result of [~<u (x) A
- 1
]
V(U 1 ) U1 + V(U 2 ) U2
V(U 1 ) + V(U 2 ) (0.90,1.03,1.60,1.36)
-1 o 2 2.5
A3 = real estate;
P highly preferable;
II = true.
all other alternatives is true." See the Note for the complete
Note
1. Step 4 of this method is actually an important research topic
set is to define some primary terms for P and IT, respectively, and
then use the context-free grammar of Zadeh [Z5] to obtain the term
With the primary terms and context-free grammar, we can get a term set
S A
S A C A
S U
A B
A NB
B T
B HT
B VVT
N Not
H Marginally
H More Or Less
H Sort of
V Very
V Highly
T Indifferent
T Better
T Worse
C Or
C And
286
true, true, very true, ... , not true, not very true, ... } can be
step 4.2. with two term sets, one for P and one for IT, we can
very true."
are defined as
~A(x)
LPiJ f [f(x) ], (centroid point of fuzzy set L i ) ,
LP i l
x
(4.193)
K
LP i l E ~L. (X j ), (4.195)
j=1 ~
K
LP i2 E S (ILL (x.», (4.196)
j=1 i J
K ILL. (X j )
~
LPiJ E (j-1) LP i l (4.197)
j=1
K ILL. (X j )
3
LP i4 E {[ (j-1) - (LPiJ)] [
~
LP i l ]}. (4.198)
j=1
(4.199)
defined by:
(4.200)
with Z using
K
L
j=1
* (4.201)
function is defined as
d * (L *
i , Z) (4.202)
statement as the label of Z. The algorithm stops here and the ranking
5.1 Introduction
Xl X2 X
n
Al xII x 12 X ln
[
D A2 X 2l X 22 x 2n
A
m
x ml x m2 xmn 1
~ (WI' w2 '···, wn )
the final utilities are real numbers, the preferred alternatives are
mixture of fuzzy and crisp data. Most of the real world MADM problems
fuzzy data. Bellman and Zadeh [BS] were the first to relate fuzzy set
theory to decision making problems. In 1977, Baas and Kwakernaak [B1]
work of fuzzy MADM method. During the past 10 years, several fuzzy
MADM methods have been proposed. The only systematic reviews of fuzzy
MADM methods have been conducted by Kickert [K1S] and Zimmermann [Z30,
Z31]. Zimmermann [Z31] among others treated the fuzzy MADM method as
the existing fuzzy MADM methods. There are a total of 1S fuzzy MADM
allowed, (3) the classical MADM method each fuzzy MADM method relates
to, and (4) the technique each method uses. Theoretical backgound as
examples are also given so that each method is easily understood. The
proposed.
Problem Correspond ing Technique
Data Type Approaches
Size MADM Methods Involved
n n
U.=
1 r w). r i )· / r w). (5.1)
j=l j=l
where r .. is the rating of the ith alternative under the jth attribute
1)
A * (5.2)
and
r i ). {(x.1)·,11 r .. »}, lJi,j,
(x 1)
ij
where y). and xi)' take their numbers on the real line ~ and 11 (y.)
Wj )
and 11 (x .. ) take values in [0,1], the utility of alternative Ai'
r ij 1)
n n
r y). xi)' / r y).. (5.3)
j=l j=l
293
n n
sup {[ 1\ ~w. (YJ')] 1\ [ 1\ ~ (x .. )]} (5.4)
v j=l J j=l
r ij ~J
the SAW method, several approaches have been proposed by Baas and
Kwakernaak [B1], Kwakernaak [K32], Dubois and Prade [028], Cheng and
McInnis [C17], and Bonissone [B27]. The first four approaches utilize
Baas and Kwakernaak [B1] were the first to identify the computa-
they proposed the use of the a-cut to obtain the fuzzy utility Ui . An
Algorithm
It is assumed that ~ (y.) and ~ (x .. ) are normalized
Wj ) r ij 1)
There may be more than one Yj value for ~w. (Yj) a o and more than
)
JL(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
O~I-",,--T-.....--+-....--i-l.-.---i-..,..........,.....l......
o 6 7 8 9 10 11 1l
Fig. 5.2 a-level set for Wj and r ij at a o 0.7.
295
J.l.(x)
1.
Fig. 5.3 Case of two fuzzy attributes and two fuzzy weights.
of Yj and x ij ' i,j 1,2. We simply pick the highest u 1" u and
i max '
the lowest u i ' u i , ' and drop all other ui's.
m1n
If the size of the problem increases, such as five attributes and
five weights, there will be 210 = 1024 u 1' values. To identify u
i max
and u, is a bit tedious without the help of a computer. When
1min
problem size increases to 10 attributes and 10 weights, there will be
Ji' (x .. ) dJi r (x 1)
.. ) / dX 1.). (5.6)
r i 1) i
and
(5.7)
Eqs. (5.6) and (5.7). If Jir:(~ij) and [Ji~j(~j)/(~ij-Ui)] have the same
The algorithm loops back to step 1 for another a value. One must
give several a values in order to get an approximated Ji U (u.)
i 1
function will be closer to the real one but will require much more
computational effort.
Numerical Example
Xl X2
o good
fair
fair
good ]
297
P.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~~~~~-+--~~~~~4-~~-+ X
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
~r
12
"
(x 12 ) ~w (Y1)
1
" ~w (Y2)
2
" = 0.75. The values are summarized as:
II II II
XII X12 Y1 ~2 U1
0.25 .8100
1 . 0 - L 0.15 .8239 (= u 1
0.25 .8083 max
. 95-L 0.15 .8227
'85_["65_[
0.25 .7900
1 . 0 - L 0.15 .8109
.55-[ 0.25 .7875
. 95-L 0.15 .8091
0.25 .7300
1 . 0 - L 0.15 .7370
'75_["6~[ 0.25
. 95-L 0.15
0.25
.7292
.7364
.7100
1 . 0 - L 0.15 .7239
.55-[ 0.25 .7083 (= u 1 . )
. 95-L 0.15 .7227 m1n
u
1
= .8100
r 11 (.85)
-5, -5,
~' ~'
r 12(.65)
±oo/ • 04 ±oo,
-5/[-.16] 31.25.
299
In this case ~ I(~ .. ) and [~I (~.)/(~ .. - ui)],Vi,j, do not h~ve the
r i 1) Wj ) 1)
same sign.
either u 1 nor u 1 . .
max m1n
On the other hand,
±00/.0261 too,
and
5/[-.1731] -28.88,
all take the negative sign. (Note that ~' (1) ±oo can be considered
wI
Similarly, we can get the fuzzy utility U2 using the same algorithm.
p.(x) U2 U1
1.......................................... .
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
Ranking of Alternatives
Note
sign. If they do, we can then compute u i using Eq. (5.3) with (X i1 ,
A
.•• , X in ' Y1' ... , y n )· Otherwise, we can proceed with another set of
combinations.
and then picking the maximum and minimum from among these values.
The example presented earlier has two attributes. Together with the
obtained. Even when there are, say, five attributes and five weights
r i ). = {(x ..
1.)
,~
r ij
(x .. »} for alternative Ai' the following steps are
1.)
-
x ij min { X ij E R I ~ (x .. ) ii!: ao }, Ifj , (5.8)
r ij 1.)
*
X ij max { X ij E R I ~r .. (X ij ) ii!: ao }, Ifj, (5.9)
1.)
Yj
- min { Yj E R I ~w.(Yj) ii!: ao } , Ifj , (5.10)
)
302
step 3. - *ij ].
At the a o level, the r ij may be represented by [Xij'X
Put X-ij ' vj, in an order such that
-
m1 s m2 s ... s mn , - (5.12)
-
where m1 min x ij and - The corresponding Yj- values
j
will also be put in an order such that
(5.13)
Similarly, x *
ij are rearranged in an order such that
(5.14)
where m*
1 min x *
ij and *
mn * .
max Xij The corresponding Yj* will
j j
be put in an order such that
(5.15)
lUx)
1.
I'(x)
1.
m1 x 12 , m2 x 13 , m3
-
x 11 ·
*
Similarly, since x 13 ~ *
x 12 ~ * , we set
x 11
m*1 * ,
x 13 m*2 * ,
x 12 m3* * ·
x 11
j *
L zk m + L
n
z- m
k=l k k=j+1 k k
min j * n ],
(5.16)
o ~ j ~ n
L zk + L zk
k=l k=j+1
j n
L zk mk* + L z * mk*
k=l k=j+1 k
ba max j ] . (5.17)
0 0 ~ j ~ n
L
-
zk + L
n
z*
k=l k=j+1 k
304
n n
E YJ" x ij / E YJ", (S.18)
j=l j=l
values. However, u" cannot be guaranteed, even when YJ", Vj, take
~min
Yj as their values. Thus, we need only be concerned with the
Numerical Example
- ll
The interval, [Xll,X * ], is obtained as:
step 3. * as:
Rearrange x ij and x ij
m1 - m2 m1* m2*
zl
- z2 zl* z2*
.15 .95 .25 1.0
zl m1* + z2* m*
2
zl + z2*
Ranking of Alternatives
Algorithm
to be:
(5.19 )
(5.20)
Recall that the actual computation of Eq. (5.1) is tedious (or even
us to focus on two points of each fuzzy number at any given time, thus
a-level sets of wj , - *
[Yj'Yj]' j = 1, ... ,n, we can obtain n a-level sets
- * of the normalized fuzzy weights P j , vj, as:
[Pj'Pj]
1, (5.23)
must hold.
m1 - ~
-
m2 ~ ~
mn (5.24)
m*
1 ~ m2* s ... s m*
n (5.25)
in which m*1 *
min x ij and mn* = max x ij
* • Eqs.(5.24) and (5.25)
j j
facilitate the construction of Eqs.(5.24) and (5.25) in the later
steps.
step 4. The smallest upper and the largest lower bound of Ui ,
u. (
d-1
E Pj*
-
mj ) + [1
d-1
- E p.* -
n
Ep. - -
md +
n
E - -
p. mj
l.min j=1 j=1 J j=d+1 J j=d+1 J
(5.26)
u imax
e-1
(E - e-1
Pj mj* ) + [1 - L p: - E
n n
* me* + E p.* mj*
p.]
j=1 j=1 J j=e+1 J j=e+1 J
(5.27)
d-1 * n
1 - E
j=1
PJ' - E
j=d+1 J
p: Zd e [P~ , P~ ], (5.28)
d-1 * n
r
j=l
PJ' + r p.
j=d+1 J
+ zd 1
which satisfies Eq. (5.23). Dubois and Prade [028] have shown that
e-1 n *
1 - r
j=l
p. -
] j=e+1
r Pj (5.29)
The values assumed by wj ' Vj, must satisfy Eq. (5.23), i.e.,
e-1 n *
r
j=l
PJ' +
j=e+1
r Pj + z e 1.
by the interval [u. ,u. ]. The OM can set several a levels and
~min ~max
repeat the algorithm several times to derive an approximated fuzzy
utility Ui .
Xl X2 X3 X4
[
A1 high mol high v. high fair
D A2 fair fair fair mol good
A3 low v. low mol high poor
1
The weight vector assigned by the DM is
Xl X2 X3 X4
A1 (0,0,.1,.3) (0, .2, .2, .4) (.8,1,1,1) (.3, .5, .5, .7)
D A2 (.3,.5,.5,.7) (.3, .5, .5, .7) (.3, .5, .5, .7) (.6,.8,.8,1)
A3 (.7,.9,1,1) (.8,1,1,1) (.6,.8,.8,1) (0,0, .1, ~3)
I·/,(x)
1.
the ..
normal~zed we~ghts
- * are computed using Eq. (5.21) and
P j , [Pj'Pj]'
(5.22). For example, the a-level set of P 1 is computed as
Y1* .85
.315
Y1* + Y2 + Y3 + Y4 .85 + .75 + .95 + .15
Y1 .75
-
Y1 3 + Y4*
+ Y2* + Y* .75 + .85 + 1. + .25
.263.
Step 3. Given the a-level sets, we can order the four lower
.355
.244
where m-1 = 0, m2 - -
.15, m3 .45, and m4 = .95 as obtained in step 3.
.215
.267
Hence, we conclude that e 2. The values which w., Vj, assume can be
J
summarized as:
.407
P.(x) U1 U2 U3
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Ranking of Alternatives
Note
most. Comparing this number with (2)2n (in Baas and Kwakernaak's
using the Baas and Kwakernaak algorithm. It takes 252 tries using the
Kwakernaak algorithm, while only 10 tries are needed using the Dubois
Algorithm
The following steps are taken for deriving fuzzy utilities.
widths of intervals.
P.(x) M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
P.(x) M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
p.(x) M
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
",(x) W1
Ci ---------------------------------------------------
/3 ---------------------------------------------
'Y ----------------------------------------
l) --------- -------------------------
"'(x) r11
Ci --------------------------------------- r - - - - - ,
~~--~--~~--~--~~--~--~-+ x
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
The first a level to be considered is the largest one among all the Wj
and r ij graphs.
step 3. Given a o ' we can obtain the a-level set for each r ij and
each Wj as:
318
(5.30)
and
(5.31)
- ij
That is, at a o ' r ij can take any value in the interval [Xij'X * ] and Wj
- *
can take any value in [Yj'Yj].
a o level, U *
imaX ' the upper bound of r ij , vj, i.e. , x ij ' must be used.
(5.33)
value will guarantee a larger u i value; on the other hand, since Yj'
m1* :S *
m2:s :S *
mn (5.34)
- as its value.
corresponding wk should take Yk *
Assume mn xiI' the
corresponding wI should take YI* as its value. * such that
For some mp
n
n L w.
j=l J
jE1Wj r ij j .. p
n n
L w. L w. + (w +A) (5.35)
j=l J j=l J P
j .. p
where A is any positive real number, holds, then the upper bound of
.
1.e., * should be selected.
yp' otherwise, y~ is selected.
Given the right combinations of y; and yj, Vj, we can easily
whose upper bounds were used for deriving u imax ' we will use their
levels are exhausted. The resulting fuzzy utilities are also discrete
and have several "steps" in it.
Numerical Example
specifications.
p.(x)
1.
.9
.8
.7
.6
.5 r-- =-
.4
.3
.2
.1
o x
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
P.(X) w2 w1
1. r--
.9
.8
.7
.6
.5 r--- '--
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
1 [.7,.9] [.9,1.]
2 [.5,.7] [.1,.3]
in Eq.(5.33), r 11 and r 12 .
w~ll * and x 12
take x 11 * ' .
respect~vely, as
their values in computing u 1
max
As for the appropriate wj ' the following steps must be taken.
.
S~nce there are two "
attr~butes, ~.e., n = * > x 12
2, and x 11 * ' we know
where m*
1 * and m2*
x 12 * ·
x 11 Based on Cheng and McInnis's search
process, w1 should take Y1* as its value while w2 should take y; as its
value.
.882.
.650.
,,(x)
1. U2 U1
.9
.8
.7
.6
.5
.4
.3
~ ~
- '----
.2
.1
o x
o .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Ranking of Alternatives
The ranking of U1 and U2 can be carried out by using any proper
ranking methodes) presented in Chapter 4. Here by simple observation,
Note
The process of finding u imax and u ioutperforms those in
min
sections 5.2.1, 5.2.2, and 5.2.3 because it requires less
computational efforts. It is doubtful, though, that Cheng and
McInnis's stepwise discrete membership functions can represent fuzzy
information better than the simple piecewise linear membership
functions as many ( Tong and Bonissone [T11], Laarhoven and Prdeycz
[L1], Buckley [B39], Godo and Sierra [G7]) pointed out that fuzzy data
323
is called the L-R type trapezoidal number (a,b,a,~) (see Fig 5.14).
Let fuzzy numbers M = (a,b,a,~) and N = (c,d,7,~), and M > 0 and N >
M(o)N (5.38)
a b a~+da b7+C~ )
M(+)N
(d' c' d(d+~)' C(C-7) (5.39)
Note that only Eq. (5.36) is exact. The remaining equations are only
P.(x)
1.
results are not exact. However, they introduce very little error.
with the help of the approximated algebraic operations, we can
quickly compute the performance of alternative Ai with respect to
attributes, x j ' j=l, ... , n, using
n
E w. r .. (5.40)
j=l J 1.J
xl X2 X3 X4
A1 high mol high v. high fair
Table 5.2 presents the fuzzy set associated with each linguistic term
as:
4
E w)' Xl)'
j=l
(1.26,1.34, .62, .64)
~~
1.
~
B
.7
~
.5
A
~
2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 B .9 1.
~(x)
1.
o
Fig. 5.16 The fuzzy utilities.
where additions and multiplications are conducted using Eqs. (5.36) and
(5.38), respectively.
Ranking of Alternatives
Note
are larger than those which were ~enerated by other fuzzy SAW methods.
Xl X2
o
good
fair
fair
good ]
! Very important, rather unimportant]
Fig. 5.4):
U2 = (.76,.76,.44,.36).
Fig. 5.17 presents U1 and U2 . Comparing Fig. 5.6 with Fig. 5.17, we
find the ordering for A1 and A2 remains the same but the spreads of Ui
in Fig. 5.17 tend to be larger. That is attributed to the fact that
Bonissone's arithmetic formulas are only approximations of the real
formulas.
,.,.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
Ol-.....-L....r--'-...--........-,......L...,........l....,-...-........--.-x
o 2
(5.41)
might equal 3/1, or 5/1, or 7/1, or 9/1. Since the numbers for the
ratio are usually taken from the set {1, 2, ... , 9}, a 15 could be sl/s2
where sl' s2 E {1,2, ... ,9}. Note that if a 15 = 3/1, then a 51 must be
The AHP method uses the pairwise comparison matrices for each
Xl X2 X
n
A1 r 11 r 12 r 1n
D A2
A
m
[ r 21
r m1
r 22
r m2
r 2n
r
mn 1
(5.42)
w (w 1 , w2 ' ... , wn )
330
where r ij and wj , Vi,j, take their numbers on the real line R. The
n n
L wJ' riJ' / L wJ'. (5.43)
j=l j=l
performance scores are also crisp, real numbers. The classical SAW
(k = O,l, ... ,Pij). The term Pij denotes the number of persons who
expressed their comparison ratios.
There are many methods one can use to derive performance scores
of the presence of fuzzy, multiple comparison ratios for the same pair
fuzzy numbers are used to obtain the fuzzy utilities, Ui , Vi' where
P.(x)
1.
framework.
Eigenvector Method
(5.44)
where
1 / a ji , 'v'i,j, (5.45)
a .. (5.46)
1)
(5.47)
nnT
A W Wi /W 2 nw
1[~:l
Wi/Wi Wi/W n
[ W2 /W i
Wn/W i
W2 /W 2
Wn /W 2
W2 /W n
Wn/Wn n
(5.48)
or
(A - n I) W o. (5.49)
Due to the consistency property of Eq.(5.46), the system of
homogeneous linear equations, Eq.(5.49), has only trivial solutions.
In general, the precise values of wi/w j are unknown and must be
estimated. In other words, human judgments cannot be so accurate that
Eq.(5.46) be satisfied completely. We know that in any matrix, small
permutations in the coefficients imply small permutations in the
eigenvalues. If we define A' as the DM's estimate of A and w'
corresponds to A', then
1: equally important
3: weakly more important
5: strongly more important
7: demonstratively more important
9: absolutely more important
location, and reputation. Since the entries of this MADM problem are
The data of this MADM problem can be derived from this hierarchy
decision structure.
The first level consists of a single objective, to have a
in the first level. The third level consists of the three jobs being
performance scores are then combined using the SAW method. The
results are the final ratings of the three jobs. The job with the
highest final rating can best fulfill the ultimate goal (successful
career) .
2 nd level
3 rd level
k
r
j=l
w (x.) Wz(YJ.)' i= 1,2, .•• ,k+1.
Yj ~
(5.52)
k
L bi]'W]~' i 1,2, ..• ,k+l. (5.53)
j=l
Thus, the priority vector of the lowest level with respect to the
element b is given as:
W (5.55)
reputation. His criteria for selecting the jobs and their pairwise
comparison matrix are given in Table 5.3. Due to the vague nature of
jobs with respect to each criteria, rather than the decision matrix.
The eigenvalue of the matrix of Table 5.3 is Amax 6.35 and the
corresponding eigenvector is
Res G B C L Rep
[ .16 .19 .19 .05 .12 .30]T
given as:
Res G B C L Rep
[
B3 A .14 .10 .32 .28 .47 .77
1
B .63 .33 .22 .65 .47 .17
C .24 .57 .46 .07 .07 .05
A
B
[ .40
.34
1
C .26
Res. G. B. C. L. Rep.
1
Research 1 1 1 4 1
2
1
Growth 1 1 2 4 1
2
1 1
Benefits 1 1 5 3
2 2
1 1 1 1 1
Colleagues 4 4 5 1
3 3
1
Location 1 1 3 1 1
3
Reputation 2 2 2 3 1 1
A 1 1 7 A 1 7 9
1
B 1 1 7 B 1 5
7
1 1 1 1
C 1 C 1
7 7 9 5
Note
saaty's AHP method does not directly use fuzzy numbers or
all a l2 a ln
[ [
A a 21 a 22 a 2n
anI a n2 ann
is derived by minimizing
(5.56)
p.(x)
1.
L (5.57)
i<j
where a ijk , k = 1,2, ... ,Pij' are Pij estimates for wi/w j . Note that
Pij can be 0 (if no comparison ratios are expressed), equal to one, or
341
greater than one, in which case there are several DMs who have
expressed their comparison ratios.
If we put Yijk = ln a ijk , zi ln wi' and Zj ln wj , we can
minimize
(5.58)
p
n n n ij
z.1 E Pij - E Pij Zj E E Yijk' vi, (5.59)
j=l j=l j=l k=l
j .. i j"i j .. i
(5.60)
(5.61)
342
Algorithm
The algorithm is shown in the following steps.
step 1. Consult with the DMs and obtain n+1 fuzzy reciprocal
matrices that take the following form (see Table 5.5).
D
(1,1,1)
(1,1,1)
a a
--~~~~!-------------------------~~~~~---
.
----------------------------------------
(1,1,1)
n n P ij
l·(EP .. ) - E E [lnliJ·k],Vi. (5.62)
~ j=l ~J j=l k=l
j .. i j .. i
343
P ..
n n n 1J
m.( I: P .. ) I: P .. m. I: I: [In mijk ] , vi. (5.63)
1 j=l 1J j=l 1J J j=l k=l
j"i j"i j"i
P ij
n n n
u.( I: P .. ) - I: P . . 1. I: I: [In U ijk ], vi. (5.64)
1 j=l 1J j=l 1J J j=l k=l
j"i j"i j"i
As In(lijk) and In(u ijk ) are lower and upper values of In(a ijk )
-In (a jik ), the following must hold true.
0, Vi,j,k.
Thus Eqs. (5.62) and (5.64) are linear dependent. The same holds for
is given as
(5.65)
(5.66)
n -1
n -1
where Al ( I: exp(u i » , A2 ( I: exp(m i » ,
i=l i=l
n -1
A3 ( I: exp(li» .
i=l
344
n
r
j=l
wJ' riJ' (5.67)
Al A2 A3
2 2
(1, 1, 1) (3' 1, 2~) (3' 1, 2~)
2 2 1 2
(3' 1, 2~) (5' 2' 3)
----2-----3--------------------2--1--2----
(3' 1, 2) (1, 1, 1) (5' 2' 3)
2 3
(3' 1, 2)
----2-----3------3-----5------------------
(3' 1, 2) (2' 2, 2) (1, 1, 1)
(~, 2, ~)
(1, 1, 1)
(1, 1, 1)
(1, 1, 1)
Al A2 A3
(1, 1, 1)
5
(2' 3, 2)
2
5
(2' 3 2)
2
5
(2' 3, 2)
2
3
(2' 2, 2)
2
3' 2)
2 1 2
(7' (1, 1, 1) (3' 1, ~)
2 2
2 1
(7' 3' 5~)
2 1
(5' 2' 3~)
2 1 2
(7' 3' 5) (1, 1, 1)
346
A1 A2 A3
3
A1 (1, 1, 1) (2' 2, ~)
2
2 1 ~)
(5' 2' 3
(1, 1, 1)
(2 1 ~) (1, 1, 1)
5' 2' 3
(3 2 ~)
2' '2
Xl X2 X3 X4
2 1 ~) 2 2 1
(1, 1, 1) (3' 2' 2 (3' 1, ~) (7' 3' ~)
2 5
2 1 2 1
Xl (5' 2' ~)
3 (7' 3' ~)
5
3 ~) 2 1
(2' 2. 2 (5' 2' ~)
3
------------------------------------------------
(~, 1, ~) (1, 1, 1) (5 3 Z) (~, 1, ~)
2' '2
(~, 2, ~) (~, 3, ~) (~, 1, ~)
(2 1 ~)
(~, 2, ~)
5' 2' 3
------------------------------------------------
5
(2' 3, Z)
2
2
(3' 1, 2~)
3
(2' 2, ~)
2
(1, 1, 1)
X4
5
(2' 3, Z)
2
2
(3' 1, 2~)
3 2 1
(2' 2, ~) (5' 2' ~)
2 3
347
3 3 3 Pij
11 E P1j E P 1j U j E E In (1 1jk ) (5.68)
j=2 j=2 j=2 k=1
3 3 3 Pij
12 E P2j - E P 2j U j E E In(1 2jk ) (5.69)
j=1 j=1 j=1 k=1
j"2 j"2 j"2
2 2 2 Pij
13 E P3 j - E P 3j U j E E In(1 3jk ) (5.70)
j=1 j=1 j=1 k=1
3 3 3 Pij
In1 E P1 j - E P 1j In j E E In(In 1jk ) (5.71)
j=2 j=2 j=2 k=1
3 3 3 Pij
In2 E P2 j - E P 2j In j E E In(m 2jk ) (5.72)
j=l j=1 j=l k=l
j"2 j"2 j"2
2 2 2 Pij
m3 E P3j - E P3jmj E E In(m 3jk ) (5.73)
j=1 j=1 j=1 k=l
3 3 3 Pij
U1 E P1 j - E P 2j l j L E In(u 1jk ) (5.74)
j=2 j=2 j=2 k=1
3 3 3 Pij
U2 E P2 j - E P 2j l j E E In(u 2jk ) (5.75)
j=1 j=1 j=1 k=1
j"2 j"2 jOO2
2 2 2 Pij
u3 E P3 j - E P3jlj E E In(u 3jk ) (5.76)
j=1 j=l j=1 k=1
31 3 - 2U 1 - u2 .405
4m 1 - 2m 2 - 2m3 - .693
3m2 - 2m 1 - m3 - .693
3m 3 - 2m 1 - m2 1.386
4U 1 - 212 - 21 3 .811
3U 2 - 211 - 13 .405
3U 3 - 211 - 12 2.238
i Ii mi ui
1 0 .087 .4812
2 -.0063 0 .3228
( .1959,.2890,.4310 )
349
3 -1
( 1: exp(u.» .1959,
i=l 1
Xl X2 X3 X4
A1 [(.196, .289, .431) (.405, .546, .714) ( .54, .579, .603) (.162, .25, .394)]
A2 (.195,.265,.368) (.162,.182,.204) (.163,.217,.292) (.313, .5,.763)
A3 (.344,.446,.561)(.277,.273, .34)(.158,.205,.267)(.209, .25,.305)
w [(.149,.194,.256),(.235,.319,.431),(.112,.140,.180),
(.263,.347,.451)]
3
1: w). r 1)' (.227, .398, .705)
j=l
350
p.(x)
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~~~----~~~~--~----~--Tx
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Note
Laarhoven and Pedrycz's approach suffers the following drawbacks.
1. There is not always a solution to Eqs.(5.62), (5.6J), and
(5.64). Indeed, as In(lijk) and In(u ijk ) are lower and upper values
of In(a ijk ) = -In (a jik ), the following must hold:
Thus Eqs.(5.62) and (5.64) are linear dependent. The same holds true
351
for Eq. (5.63). There may not be any solution to Eqs.(5.62) through
(5.64) .
z3 = (.5634,.5200,.7439)
the OMs. For example, in Fig. 5.22b (5,6,8,9) represents the fuzzy
[
~ a .. ] lin
j=l 1.)
",(x)
1." ---------------~---"'"
",(x)
about 5 to 1
1.
.9
.8
.7
.6
.5-
.4
.3
.2
.1
00 2 3 4 5 6 7 8 9 10 X
presented.
Let M1 = (a 1 , b 1 , c 1 , d 1 ) and M2
(1) Addition
Q (5.83)
(2) Multiplication
Q (5.84)
where
follows. For any unique x on the horizontal axis, ~Q(x) can be:
x ~Q(x)
:s a 0
d
'"
b:sx:sc
0
1
a:sx:sb ex E [0,1]
c:sx:sd ex E [0,1]
[a 2 ,b 2 ] where
1,2, (5.85)
355
the product x
x (5.86)
(5.87)
(d 1 +d 2 ) [R11+R21,R12+R22])} (5.88)
x ~Q' (x)
~ (a 1 +a 2 ) 0
z:: (d 1 +d 2 ) 0
(b 1 +b 2 ) :S x ~ (c 1 +c 2 ) 1
(a 1 +a 2 ) ~ x ~ (b 1 +b 2 ) a E [0,1]
(c 1 +c 2 ) ~ x ~ (c 1 +c 2 ) a E [0,1]
the addition x xl + x 2 will take the form
x (5.89)
x (5.90)
Algorithm
OM is assumed for the following steps. The case of multiple OMs shall
be explained in the Note section.
(5.91)
is given as:
(5.92)
The following will detail the derivation of fuzzy weight wi. Let
n l/n
[ II «b" - a 1'J') 0: + a 1'J')] , 0: e [0,1], (5.93)
j=l 1J
n
g,(o:) = [ II «c .. - d 1,),) 0: + b, ,)]l/n, 0: e [0,1], (5.94)
1 j=l 1) 1)
n
[ II a .. ] l/n (5.95)
j=l 1)
and
m
a = L ai · (5.96)
i=l
b, c, di
1 1
c:,~, a ), vi, (5.97)
x gw, (x)
1
~ (ai/d) ) 0
;,: (di/a) 0
[bi/c, ci/b] 1
x (5.98)
x (5.99)
where
m
f(a) E fiCa) (5.100)
i=l
m
g(a) E gi(a) (5.101)
i=l
n
E wJ" r iJ", vi. (5.102)
j=l
step 1. The expert estimates the relative weight ratios for each
pair of alternatives under every attribute as well as the relative
weight ratios for the attributes. This results in four reciprocal
comparison matrices:
Xl A1 A2 A3
1 1 1 1 1 1
A1 1 (4' 3' 3' 1)
2 (2' 2' 2' 1)
2
A2 (2, 3, 3, 4) 1 (1, 1, 2, 2)
1 1
A3 (2, 2, 2, 2) (2' 2' 1, 1) 1
X2 A1 A2 A3
A1 1 ( 6, 6, 6, 7) (2, 2, 4, 4)
1 1 1 1) 1 1
1 (2' 1, 1)
A2 (7' 6' 6' 6 2'
1 1 1 1)
A3 (4' 4' 2' 2 (1, 1, 2, 2) 1
X3 A1 A2 A3
A1 1 (1, 2, 2, 3) (7, 8, 8, 8)
1 1 1
A2 (3' 2' 2' 1) 1 (3, 3, 4, 4)
1 1 1 1 1 1
A3 (8' 8' 8' 1)
7 (4' 4' 3' 1)
3
1
Xl X2 X3
1 1 1 1 1 1
Xl 1 (7' 6' 6' 1)
5 (3' 2' 2' 1)
X2 (5, 6, 6, 7) 1 (3, 3, 3, 3)
1 1 1
X3 (1, 2, 2, 3) (3' 3' 3' 1)
3
1
360
step 2. For the first reciprocal matrix, the geometric mean is:
3
( II a .)1/3 )1/3 (1* 1* 1)1/3 .5.
a1 (all * a 12 * a 13 4 2
j=l 1J
and
3
a2 ( II a .)1/3 1. 2599
j=l 2J
3
a3 ( II a .)1/3 1. 0000.
j=l 3J
Hence,
3
a = E a. 2.7599
i=l 1
summarized as:
i 1 2 3
a1 b1 c1 d1
r 11 «(i' c
, "'1) , -a) (.1285, .1517, .1839, .2283),
a2 b2 c2 d2
r 21 «(i' c
, "'1)' -a) (.3239, .3976, .6072, .7247),
a3 b3 c3 d3
r 31 «(i, -C' "'1) , -a) (.2571, .2757, .4210, .4565) •
361
Xl X2
Al[
(.1285, .1517, .1839, .2283) (.4991, .5162, .8594, .9107 )
X3
(.3957, .6106, .6495, .8936) A1
W3
(.1596, .2222, .2222, .2839)]
{.0107[.00064,.00549],.0168,.0204,.0379[.00244,-.0199]}
Rl (d 1 - c 1 ) (d 2 - c 2 )
R2 - [d 2 (d 1 - c 1 ) + d 1 (d 2 - c 2 )]·
summarized as:
1 {.0107[.00064,.00549],.0168,.0204,.0379[.00244,-.01990]}
2 {.2834[.00169,.05907],.3441,.5730,.7134[.00598,-.14637]}
3 {.0632[.01345,.05907],.1357,.1443,.2537[.01506,-.12444]}
U1 = {.3573[.0158,.1236],.4966,.7377,1.0050[.0235,-.2907]}.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 ~~~~--~~--~--~~--~--r-~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Note
The proposed approach can be easily extended to incorporate
making case.
Jl J2 J3 J4 J5
[
Xl (7,8,8,9) (5,5,6,6) (6,7,7,8) (8,9,9,9) (5,5,5,5']
X2 (4,4,5,6) (7,7,7,7) (0,0,0,1) (2,2,3,3) (5,6,6,7)
Jl J2 J3 J4 J5
Al [ (5,5,5,5, (2,2,2,2) (7,8,9,9) (4,4,5,5) (0,0,0,"]
A2 (5,5,5,5) (7,8,9,9) (5,6,6,7) (9,9,9,9) (7,8,8,9)
Al
A2
r Jl
4 '4'S'S'
(5,5,5,5)
J2
(5,5,6,6)
(5,6,6,7)
J3
(5,5,6,7)
(5,5,5,7)
J4
(4,4,6,6)
(9,9,9,9)
J5
(4,5,6,7']
(5,6,6,7)
as:
i r il r i2
determined using
n Wj
[,r (--yo- ) riJ,]/n (5.103)
J=l
(5.104)
by Buckley that values like 0.7, 0.8, 0.9 might be appropriate for u.
e ij e ji = 1.
We shall calculate e 23 and e 32 to demonstrate the computational
JL(x)
1.
~~------~--~----------~r-~~x
2.968 3.080 3.534 3.828
summarized as:
i\j 1 2 3
1 .112
e ..
1.J
2
3
[ 1
1 .402
Let's assume that ~ 0.9. The rank of the fuzzy utilities can
be determined based on
Remark
good as its fuzzy method peers), we shall modify the numerical example
xl A1 A2 A3 x2 A1 A2 A3
A1 1 1/3 1/2 Al 1 6 3
A2 3 1 2 A2 1/6 1 1/2
A3 2 1/2 1 A3 1/3 2 1
X3 Al A2 A3 w xl x2 x3
Al 1 2 8 Xl 1 1/6 1/2
A2 1/2 1 3 x2 6 1 3
A3 1/8 1/3 1 x3 2 1/3 1
Xl X2 X3
0 Al
[.53 96 .6667 .0964 ~
~ [.UU]
A2
A3
.1634
.2969
.1111
.2222
.7006
.2029
] .6667
.2222
.5153
u o x W [ .2479
.2262 1
369
approach be represented by
1, i 1,2,3.
only part of the fuzzy utility, may be as good as the fuzzy utility in
less computations.
fuzzy number)
xO: ij
k ..
1J
, x > 0
{x + d ij ) (3 1J
..
1
/.la .• (x) (5.106)
1J 0 , X :s 0
370
a }, Va. (5.107)
28 (see Fig. 5.25a). Consequently, ITQ (X=7) = .20 and ITQ (X=28) 0.0.
~PIQ(a=o.o) = sup [.20, 0.0] = .20 based on Eq. (5.107) (see Fig.
5.25b). By setting different a levels, we can obtain ~PIQ(a) as
denoted by the solid line and the point (1,1) on Fig. 5.25b. Note
that the dashed line from point (0,0) to (1,1) on Fig. 5.25b indicates
",(x)
1.
.9
.8
.7
.6
.5
~
.4
. .'
~
~
.3 ~\
~\~
.2 ·-------I--
I-
.1 .' -
l:
o 0~-5~-1...0------~--- x
15 20 25 30
/L. «~)
Pia (11)
1. ---- ---- ------- --- ------------ --- --- -------------------;:,
.8 .... _______________________________________ ::;,6: _____
.9 _- ----- i.
.7 ;"i'
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1. a
The two indices are (1) the possibility of matching IT(P;Q), and (2)
which is also a crisp number in [0,1]. When both P and Q are fuzzy,
x
(5.113)
and
N(PjQ) inf max (X,/.lPIQ(X» S max (~,1 - /.lPIQ(~» = ~.
x
(5.114 )
,.,.(x)
1.
.9
.8 U(~!~L. __ _
.7
.6
.5
.4
.3
min (S .11S)
j=l, ... ,n
min (S.116)
j=l, ... ,n
o 0 0 0 0
where A = (Xl'···' Xn ), Ai = (X i1 ,···, x in ), and Xj and Xij are
defined on the same domain U. The vector AO is the cutoff vector
specified by the DM, while Ai' i = 1, ••• , m, is the vector that
contains the performance scores of the ith alternative under all
attributes.
Eqs.(S.llS) and (S.116) suggest that the matching is done
attribute by attribute. These matching results are to be aggregated
using the min operator to preserve the respective semantics of
possibility and necessity of the indices. Eqs.(S.llS) and (S.116)
376
unequal weights are used, one of the following formulas can be used.
max Wj 1, (5.119)
j=l, ... ,n
what extent we are certain that the fuzzy set of importance is includ-
0
Sj II(Xj,X ij ), j 1, ... ,n, (5.120)
or
Sj
0
N(Xj,X ij ) , j 1, ... , n. (5.121)
0 0
min max ( 1-Wj,N(XjiXij» N(A iAi). (5.123)
j
n
L wJ' I-Lp, (x), "</ X e U. (5.126)
j=l J
Note that Eq. (5.126) does not have the intersection or union
event ~x defined on the crisp set { 1, ... , n}. ~x is the fuzzy set
(5.127)
If Then
However, when the stated rule is not followed, we can't say that Ai >
Aj nor Aj > Ai. Also note that the necessity index, N, is more
important than the possibility index, II, because when the N index is
positive we can be certain that the alternative (more or less) matches
the requirements set by the OM.
p.<x)
1.
.7
.6
.5
.4
.3
.2
.1
0
0 2 4 6 8 10 12 14 16 18 20 X
.3
.2 f
l
expensiv
.1 §~
°0~~~2~~3~~4~~5--~6~~7~~8--~9--+10X
~ : f
.6 ~ :
~ .
.5
~ ~
.4
.3 ~
.2 ~
~
.1 ~
o ~~--~~~~~~~--~~--~~x
~
o 2 4 6 8 10 12 14 16 18 20
between 140-160
/L(x)not very fast ~rather fast
1. \ .. '" 'I""""" fast
.9 : l"
. \ ~~ /
8 ! ~ : ~
• : ~ ". ~$
.7 : \ :, ~
.6 ! ': ~
" • $
.5 ,: ~
.4 ." ~\
.3
,2 \ ~
.1 \
~
Ol--...........L...,i-i-..,....I.....I-...;....-+--I---.--....--..... x
o 100 120140 160180200220240260
Fig. 5.2Sd The fuzzy set of used car maximum speed (mph).
382
o 0 0 0
«xl' 0.8), (X 2 ' 0.5), (x 3 , 1), (x 4 , 0.2»
where the numbers 0.8, 0.5, 1, and 0.2 are weights associated with
each attribute, and the xj, j =1,2,3,4, are summarized as:
o
The necessity of matching between Xl and x 11 is (see Fig. 5.29):
P.(x)
1.
.9
.8
.7
.6
.4 (
(
.3 ~
.2 ~~
.1 ~
o O~~-2~~4--~6--~8---1~O--1~2--1~4--~16--~18--~20X
i 1 2 3 4 5 6 7
Il(A O ; Ai) .5 .3 o o o .2 .8
N(Ao ; A.) o o o o o o .5
Note
This method utilizes the min operator to combine the
single-attribute index values for the conjunctive case (Eq.(S.117»
and utilizes the 'max' operator to combine the single-attribute index
values for the disjunctive case (Eq.(S.118». In both cases, only a
small portion of the available information is used in deriving the
final conclusions. That is, the matching of one alternative over
another with a cutoff vector is determined solely on one attribute.
Thus, even if an alternative matches fairly well on all but one
attribute, another alternative that matches about average on all
attributes would be chosen over it (for the conjunctive case).
Similar remarks apply to the disjunctive case.
The applicability of this method is rather limited. In general,
it would only be reasonable to use this method if the OM is assumed to
have a pessimistic nature in the decision process (for Eq. (5.117»,
while Eq. (5.118) should only be used when the OM is assumed to have an
optimistic stance toward the decision situation.
385
represented in the "IF .•• THEN ... " decision rule format.
Algorithm
decision environment.
and analysts takes place; and in Step 4, where the decision rules are
Special
Price Downtime Features UtilitX
intolerable
::s $400 low poor poor
sa! $600 moderate medium
$800 high good good
'"
very good very good
ideal
IF THEN
Price Downtime special Features utility
:s $400 moderate good good
moderate, and the special features are good, then the utility is
Xl X2 X3
D A1 $520 moderate poor
[
A2 $455 moderate-high very good
A3 $430 low good
A4 $875 moderate good 1
where the term "moderate-high" may be interpreted as a fuzzy set
{ ( . 5/moderate), (. 5/high) } .
can simply match the performance data with the decision rules in
Table 5.14 and determine the fuzzy utility of each alternative. For
example, because one of the rules in Table 5.14 is matched, the fuzzy
IF THEN
Price Downtime Special Features Utility
~ $800 moderate good very poor
388
needed:
389
or
moderate, poor] with the decision rules in Table 5.14, we obtain fuzzy
U1 {(.4/medium), (.6/medium-poor)}.
respectively. Note that the membership values .725 and 275 are
{(l/very good)}
good] is
1 {(.4/medium),(.6/medium-poor)}
3 {(.85/good),(.15/medium-good)}
4 {(l/very poor)}
391
Note that these fuzzy utilities are fuzzy sets which contain
linguistic terms and the real numbers e [0,1] associated with each
linguistic term may be treated as the membership value. The member-
ship values of a fuzzy set (e.g., fuzzy utility) need not sum to one.
step 7. Rank the fuzzy utilities. Let's assume the meaning of
the atomic linguistic terms {good, medium good, medium, medium poor,
poor} can be shown in Fig. 5.30 or expressed in discrete membership
functions as:
x 0 .1 .2 .25 .3 .4 .5 .6 .7 .75 .8 .9 1
good 0 .16 .3 .7 1
med.goOd 0 .3 .7 1 .7 .3 0
medium 0 .5 1 .5 0
med.poor 0 .3 .7 1 .7 .3 0
poor 1 .7 .3 .16 0
,.,.(x)
med.poor medium med.good
1.
. 9-
.8
.7
.6-
.5
.4-
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
2
Ilvery good(X) {Ilgood (x)} •
.7 .75 .8 .9 1
very
Xgoodl 0 .0256 .09 .49 1
x o .1 .2 .25 .3 .4 .5 .6 .7
(I) 0.4/medium o .4 .4 .4 o
(II) 0.6/med. poor o .3 .6 .6 .6 .3 o
Using the max operator to combine (I) and (II), we can obtain U1 as:
o .1 .2 .25 .3 .4 .5 .6 .7
o .3 .6 .6 .6 .4 .4 .4 o
1
x 0 .1 .2 .25 .3 .4 .5 .6 .7 .75 .8 .9 1
/.lU (x) .275 .275 .275 .16 0 .275 .275 .275 0 .16 .3 .5 .5
2
/.lU (x) 0 .15 .15 .16 .3 .7 .85
3
/.lU (x) 1 .49 .09 .0256 0
4
JL(X)
1.
.9
.8 ~,4
.7 "
.6
.5
.4
.3
.2
.1
0 ~-T--~~~-T--~--r-~~~--r-~x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Note
(5.133)
394
{(l/very good)}.
JiA (x..
(n) ) min (Ji A (x .. » (5.134)
i 1) j=l, ... ,n i 1)
the help of the OM's heuristics. For a small problem like the
= 1,048,576 decision rules, and even with the help of the OM's
the OMs and analysts is time consuming. When the OMs are not
Negi [N7] believes that the MADM problem containing fuzzy data
can be solved using classical MADM algorithm in which fuzzy data are
applied to any type of MADM problem involving fuzzy data. Negi left
should be used.
We shall use one of the well known classical MADM methods, TOPSIS
matrix, D, of m x n dimension:
Xl X. X
) n
D A1 x 11 x 1j x 1n
A. Xu x .. X.
1 1) 1n
Am x x x
m1 mj mn
Algorithm
P.(x) xij Wj
1.
.9
.8
.7
.6-
.5
.4
.3
.2
.1
0 x
{ r ij -
Xj / x ij , vj, X. is a cost attribute
J
(5.135)
Xl Xj X
n
D' A1 r 11 r 1j r 1n
Ai r i1 r ij rin
A r m1 r rmn
m mj
d ij ) and Xj* = (a *j , b *
j , c*
j , d*
j ), we have
397
a .. b ij c ij d ..
r ij Xij (+) Xj* (21. , -¥) (5.136a)
d*
j ' c *j b*
j ' aj
- -
- a.
(~, c
bj , c. d:
r ij Xj (+) x ij
~J
ij b~j' ....l.-)
a ij (5.136b)
(5.137)
fuzzy operations:
a .. b ij c ij d ..
v ij r ij ( . ) Wj (-¥ aj , ~j' T.,-¥ .s.) (5.138)
dj c *j b*
j J a. J
J
-
a. - d:
V ij r ij (.) Wj (....l.- a j ,
d ij
bj
c ij ~j' b ij
l T.,.:.:.L
J a
.s.)
J
(5.139)
ij
v
398
A * *
[VI' ... , Vn* ) ,
A - [V~, ... , Vn ) ,
(5.140)
For each column j, we find a v ij which has the greatest mean as the Vj*
and a v ij which has the lowest mean as the v j .
n
Si* r * i 1, ..., m• (5.141)
j=1 °ij'
and
n
si r °ij' i 1, ... , m• (5.142)
j=1
399
* v .. - *
Dij 1.) Vj
Di j v .. - Vj
1.)
The computation is straight forward. For fuzzy data, Negi defined the
difference between two fuzzy numbers g (x) and g *(x) (based on Zadeh
v ij Vj
[Zl]) as:
(5.145)
p.(x)
1.
Xl X2 X3 X4 X5 X6
(Mach) (NM) (lbs) ($M)
D A1 1500 20000 5.5 avg. v. high
[ 2.0
A2 2.5 2700 18000 6.5 low avg.
1.8 2000 21000 4.5 high high
1
A3
A4 2.2 1800 20000 5.0 avg. avg.
The decision matrix and the weight set are then rewritten as:
401
Xl X2 X3 X4 X5
D A1 2.0 1500 20000 5.5 (.3, .5, .5, .7)
[
A2 2.5 2700 18000 6.5 ( .1, .2, .2, .3)
A3 1.8 2000 21000 4.5 (.7, .8, .8, .9)
A4 2.2 1800 20000 5.0 (.3, .5, .5, .7)
X6
(.9, .95, .95, 1. 0) A1
(.3, .5, .5, .7) A2
( .7, .8, .8, .9) A3
( .3, .5, .5, .7) 1 A4
X2 , X3 , and X4 using
r ij -
Xj I x ij , for cost attribute X4 •
Xl X2 X3 X4 X5
R A1 .80 .55 .95 .82 ( . 3, .5, .5, .7)
[
A2 1.00 1. 00 .86 .69 ( .1, .2, .2, .3)
A3 .72 .74 1.00 1.00 (.7, .8, .8, .9)
A4 .88 .67 .95 .90 ( . 3, .5, .5, .7)
X6
( .9, .95, .95, 1. 0) A1
( .3, .5, .5, .7) A2
( .7, .8, .8, .9) A3
( .3, .5, .5, .7) 1 A4
402
the weights are fuzzy data, we shall apply fuzzy multiplication on Xl'
.675,.75), we have
By the same token, we can obtain all the vij's and present them as:
Xl X2
v (.48,.54,.54,.60) (.22, .275, .275, .33)
(.6,.675,.675.75) (.4,.5,.5,.6)
(.432, .486, .486, .54) (.296,.37, .37.444)
(.528, .594, .594, .66) (.268, .335, .335, .402)
X3 X4
(.38, .475, .475, .57) (.328, .41, .41, .492)
(.336, .43, .43, .516) (.276, .345, .345, .414)
(.4,.5,.5,.6) (.4, .5,.5,.6)
(.38, .475, .475, .57) (.36,.45,.45,.54)
X5 X6
(.225, .412, .412, .63) (.81,.903,.903,1) A1
(.075, .165, .165, .27) (.27, .475, .475, .70) A2
( .428, .66, .66, .81) (.63, .76, .76, .90)
1
A3
(.225, .412, .412, .63) (.27, .475, .475,.70) A4
summarized as:
403
Xl X2 X3 X4 X5 X6
M A1 .54 .275 .475 .41 .422
.904 :]
[
A2 .675 * .500 * .43 .345 .170 .481
A3 .486 .37 .50 * .50 * .633 * .763
A4 .594 .335 .475 .450 .422 .482
Once the v ij ' vi,j, are obtained, Vj* and Vj- are easily identified
as:
A* (V 1* , V2* ' V3* ' V4* ' V5* ' V6* ) (V 21 , V22 ' V33 ' V34 ' V35 ' V16 )
A - - -
(V~, V2 ' v 3 ' V4 ' V5 ' V~)
- - (V 31 , V12 ' V23 ' V24 ' V25 ' V26 )
Since all the data are fuzzy, the difference measures between v ij
Xl X2 X3 X4 X5 X6
D* A1 1 1 .128 .495 .551 0
[ ]
A 0 0 .376 .917 1 1
A3 0 .747 0 0 0 .571
A4 .575 .988 .128 .263 .551 1
404
Xl X2 X3 X4 X5 X6
-
D A1 .474 0 .249 .430 .846 1
[
A 1 1 0 0 0 0
1
A3 0 .736 .376 .917 1 .803
A4 .900 .492 .249 .660 .846 0
D~l 1 - {sup
x
Jl v
11
(x) 1\ Jl
v1 - (x) ] } 1 - L~l .474
6
*
Sl *
L D1j 1 + 1 + .128 + .495 + .551 + 0 3.174,
j=l
6
S~ L D~j .474 + 0 + .249 + .430 + .846 + 1 2.999.
j=l
i 1 2 3 4
i 1 2 3 4
Note
1. The algorithm of TOPSIS is not the only classical MADM method
available. There are other approaches such as Simple Additive Weight-
ing method, Linear Assignment Method, ELECTRE, and so on. Note that
the original algorithms must be modified slightly so as to incorporate
fuzzy mathematics or to simplify the computation. For example, in
classical TOPSIS normalization of the decision matrix is done using
m 2 1/2
x .. /
~J
[r
i=1
x .. ]
~J
(5.146)
surely better than Al (Roy [R9], Hwang and Yoon [H13]). Through the
successive assessments of the outranking relations of the other
ELECTRE method (see Roy [R9,R10], Hwang and Yoon [H13]). Since then
Takeda [T2], Siskos et al. [S20], Brans et al. [B35], and Martel et al.
outranking relations.
multilevel graph. Brans et al. [B35] assist the analyst with six
the decision data. A confidence index and a doubt index are computed
first which in turn are used to calculate fuzzy outranking relations.
We shall present Roy's approach for its originality, siskos et
aI's approach and Brans et aI's methods for their simplicity. We will
only present the basic concepts and the algorithm of Takeda's approach
because it is too cumbersome to use. Martel et al. approach is more
difficult to use than Takeda's approach and will not be presented
here. Interested readers should refer to the original paper [Mll].
Definitions
attribute, respectively.
1. reflexive: Sd(k,k) = 1, Yk
2. nonsymmetrical: Sd(k,l) does not necessarily equal Sd(l,k).
d d
Ak is preferred to Al iff S (k,l) > S (l,k);
Roy [R11] proposed the use of the degree of concordance and the
proposition "Ak over AI" with respect to the jth attribute, is defined
i, ~ r
1 ,r
kj + t J lj (5.147)
r (r + tl?
lj - kj J ti,
-..::::...<'---:.---....:..:;~----"''---, rkJ' + J $ r l J'
t~ - tl?
J J
dj(k,l) (5.148)
r + t P,
t~ - tl? kj J
J J
n
C(k,l) L wJ' CJ,(k,l), Vk,l, (5.149)
j=l
where wj , vj, are the weights assigned by the OM. The concordance and
1 n
O(k,l) n L f(dJ,(k,l),C(k,l» (5.150)
j=l
where
elk,1)
1"1---------,.
Algorithm
Roy's approach will be presented as a series of successive steps
as follows.
step 1. Initialization. The decision maker must set the
413
i
thresholds: (1) indifference threshold, t j , (2) preference threshold,
order.
Xl X2 x3
A1 5.2 5.6 2.6
[
A2 7.1 4.5 1.8
3.8 8.0 6.5 (5.153)
1
A3
A4 6.4 4.2 9.0
414
Al A2 A3 A4
A1 1 .067 1 .53
C1 A2
A3
A4
[ .40
1
1
1
0
.867
1
1
1
1
0
1 1
For example, we can obtain C1 (1,2) using Eq.(5.l47) as:
7.1 - 7.2
C1 (1,2) i p
.067,
.5 - 2
A1 A2 A3 A4
A1 1 1 0 1
[
A2 .60 1 0 1
C2
1 1 1 1
1
A3
A4 .40 1 0 1
A1 A2 A3 A4
A1 1 1 0 0
[
A2 .80 1 0 0
C3
1 1 1 0
1
A3
A4 1 1 1 1
A1 A2 A3
0 0 0
.~o 1
A1
[
A2 0 0 0
d1
A3 0 .433 0
A4 0 0 0
d 1 (3,2)
7.1 - (3.8+2)
.433,
5 - 2
because r 31 + ti ~ r 21 ~ r 31 + t~ is found.
A1 A2 A3 A4
A1 0 0 .133 0
[
A2 0 0 .500 0
d2
1
A3 0 0 0 0
A4 0 0 .600 0
A1 A2 A3 A4
A1 0 0 .633 1
[
d3 A2 0 0 .90 1
A3 0 0 0 .167
A4 0 0 0 0 1
step 4.1 Construct the aggregated concordance relation, C, using
Eq. (4.149) as:
A1 A2 A3 A4
A1 1 .767 .25 .533
[
A2 .77 1 .25 .65
C
1
A3 .85 .750 1 .40
A4 .76 .967 .60 1
C(1,2) .767.
A1 A2 A3 A4
A1 1 1 .830 .667
[
A2 1 1 .600 .667
D
A3 1 1 1 1
A4 1 1 1 1 1
417
For example, we can obtain 0(1,2) using Eqs.(5.150) and (5.151) as:
f(d 2 (1,2) ,C(1,2» = f(d 3 (1,2),C(1,2) = 1 because C(1,2) > d 2 (1,2) and
A1 A2 A3 A4
A1 1 .767 .208 .356
[
sd A2 .77 1 .15 .434
.85 .750 1 .40
1
A3
A4 .76 .967 .60 1
relation sd as:
A1 A2 A3
A1 1 0 0
[
sd A2 0 1 0
1 0
A3 0 1
A4 0 0 0 1
where none of the alternatives outrank the others.
418
when the Ir kj - rIjl differences are very small and can be considered
nonsignificant.
1, i f r Ij - r kj :s 0,
Cj(k,l}
{ 0, i f r Ij - r kj tj , (5.154)
~
r Ij - r kj
[1 - ]L, L > 0, otherwise
t,
J
p,(x)
C.(k,1)
1.1----+-'---..
attribute.
relations into a single one that takes into account the weight of each
n
C(k,l) L WJ. C J' (k,l), (5.155)
j=l
Al.
as:
I
1, i f r Ij - r kj ~
rj
dj(k,l) 0, i f r Ij - r kj ~
tj (5.156)
r Ij - r kj
]L, L >0, otherwise
v. - t.
J J
p.(x)
1.
incomparability.
If Then
C(k,l) 1, Vk,l Cj(k,l) 1
and
dj(k,l) 0
C(k,l) (5.157)
{IT*[l - d *(k,l)]} 1-C(k,1) ,
j j
for { jl C(k,l) < dj*(k,l)}.
424
The second part of this definition allows us to take into account only
the most significant discordances in the calculation of Sd (k,l). In
the case of dj*(k,l) =1 for at least one attribute j*, the attribute
Xj * plays the role of "dictator" to the outranking of Ak over AI.
{ 0, otherwise. (5.158)
since ~D(k,l) is the fuzzy set of alternatives AI' VI, that are
dominated by ~, it is natural to define the nondominance relation
as the complement of ~D(k,l). THat is a nondominance relation can be
defined as: (see Section 3.2.2.1 Complement of a Fuzzy Set)
1 - ~D(k,l). (5.159)
425
Vk, that are not dominated by AI. The intersection of all ~ND(I,k)
for all alternatives will result in the fuzzy set of alternatives that
1 - max ~D(I,k)
Algorithm
outranking relation.
stee 8. * (A )
Obtain a set of nondominated alternatives using, Ilk .
stee 9. Identify the best alternative, A* , with
Il * (A * ) max Il * (~).
k=l, ... ,m
Xl X2 X3
step 1. Initialization.
veto threshold vector ~, and the weight set ware given by the OM as:
A1 A2 A3 A4
A1 1 .05 1 .40
A2 1 1 1 1
A3 .30 0 1 0
A4 1 .65 1 1
(7.1 - 5.2)
C1 (1,2) = 1 - 2
0.05.
A1 A2 A3 A4
A1 1 1 0 1
C2 A2 .45 1 0 1
A3 1 1 1 1
A4 .30 .85 0 1
A1 A2 A3 A4
A1 1 1 0 0
C3 A2 .60 1 1 0
A3 1 1 1 0
A4 1 1 1 1
428
A1 A2 A3 A4
C(l,2) .76.
A1 A2 A3 A4
A1 0 0 0 0
d1 A2 0 0 0 0
A3 0 .43 0 .20
A4 0 0 0 0
d 1 (3,2) 0.43.
A1 A2 A3 A4
0 0 .13 0
0 0 .50 0
0 0 0 0
0 0 .60 0
A1 A2 A3 A4
0 0 .63 1
0 0 .90 1
0 0 0 .17
0 0 0 0
A1 A2 A3 A4
A1 1 .76 .12 0
sd A2 .64 1 .02 0
greater than both d 1 (1,3) = 0 and d 2 (1,3) = .13 but is smaller than
C(1,3) .12.
1-C(1,3)
A1 A2 A3 A4
A1 0 .12 0 0
/.1 0 A2 0 0 0 0
A3 .71 .73 0 0
For example, since C(4,3) .60 > C(3,4) 0.40, the value of /.10(4,3)
is
• 60 - .40 .20 •
A1 A2 A3
A1 1 .88 1
/.I NO A2 1 1 1 1
A3 .29 .27 1 1
min [1,.29,.28]
.28
summarized as:
431
k 1 2 3 4
Since the alternative that has the highest j.1* value is considered
Note
For example,
3
E wJorlJo (.25) (.73)+(.40) (.70)+(.35) (.29) .564.
j=1
rO1J°
the best value among x ij . The final ranking order is A3 > A4 > Al >
A2 •
clear, however, that the classical SAW method is much easier to apply,
Generalized criterion
Let Xj be an benefit attribute, i.e, bigger better. When two
alternatives, Ak and AI' are compared with respect to Xj ' the result
of the comparison must be expressed in terms of preferences.
Therefore, a preference function P(k,l) which gives the intensity of
the preference of ~ over Al is defined. The preference structures
pertaining to P(k,l) and their meaning are summarized below.
There are four possible preference relations between alternatives
Ak and AI· They are:
433
P(k,l) , if d ~ 0,
R(d)
{ P(l ,k) , if d ~ o.
............. _;;,;;_._---
H(d)
PREFERENCE 0
AtOVER"
~------~~~------~d
0, d = 0,
B(d) { 1, Idl > 0,
(see Fig. 5.39) (5.162)
0, Idl :s u,
B(d) { 1, otherwise,
(see Fig. 5.40) (5.163)
I v , Idl :s v,
B(d) { Idl1, otherwise,
(see Fig. 5.41) (5.164)
435
0, Idl :S u,
R(d) { 0.5, u :S Idl :S v, (see Fig. 5.42) (5.165)
1, otherwl.se,
0, Idl :s u,
R(d) { (Idl - u) I(v - u),u :S Idl :s v, (5.166)
1, otherwise,
( see Fig. 5 . 43 )
d 2
R(d) 1 - exp{- -----}. (see Fig. 5.44) (5.167)
2 0- 2
H(d)
~-----------±-------------d
H(d)
~------_~q----~I~--~q------~d
H(d)
~----=-----~----~----~d
H(d)
~------~--~------------~d
-p -q q p
H(d)
~---7---_~q--~--~q~--~p----d
H(d)
~~------~~~----~--~d
a
Fig. 5.44 General criterion Type 6: Gaussian criterion.
Each of the shapes can be very easily defined because only one or two
parameters are to be fixed. The possible parameters are:
1. u j is an indifference threshold--the largest d j value below
which the DM considers there is no difference.
2. v. is a strict preference threshold--the lowest value of d.
J J
above which the DM considers there is strict preference.
3. u j is the standard deviation of a normal distribution.
These formulas do not represent an exhaustive list. other shapes of
the H(d) function could be considered, too.
If criterion Xg is to be minimized, then the formula for
calculating the difference between ~ and Al must be defined as:
(5.168)
Algorithm
step 1. Initialization.
are determined.
'Vi (5.169a)
'Vi (5.169b)
A1 A2 Am
I
sd A1 Sd(1,2) sd(1,m) (5.171)
A2 sd (2,1) sd(2,m)
Am sd(m,1) sd(m,2)
439
F(k) (5.174)
If Then
F+(k) = F+(l) 1+
~ Al
If Then
(5.175)
(5.176)
(5.178)
Note that the center of [Zk'Y k ] is the mean of F(k) and the
all the intervals be less than the mean distance between two
If Then
~ is indifferent to Al
matrix, o.
442
Step 1. Initialization.
After serious discussion with the DM, the analyst concludes that
in Table 5.15.
Xl II u 10
X2 III v 30
X3 V u 0.5; v 4.5
X4 IV u 1; v 5
X5 I
X6 VI (j = 5
d j or d g Xl X2 X3 X4 X5 X6
(A 1 , A2 ) -15 32 -4 4.3 -7 4
(A 1 ' A3 ) 3 30 -2 1.8 -4 -2
(A 1 , A4) -40 10 4 2.1 -1 -5
(A 1 , A5) -28 18 0 -3.4 -5 -3
(A 2 , A1 ) 15 -32 4 -4.3 7 -4
(A 2 , A3 ) 18 -2 2 -2.5 3 -6
(A 2 , A4) -25 -22 8 -2.2 6 -9
(A 2 , A5) -13 -14 4 -7.7 2 -7
(A 3 , A1 ) -3 -30 2 -1.8 4 2
(A 3 , A2 ) -18 2 -2 2.5 -3 6
. (A 3 ,A4) -43 -20 6 0.3 3 -3
(A 3 , A5) -31 -12 2 -5.2 -1 -1
(A 4 , A1 ) 40 -10 -4 -2.1 1 -5
(A 4 , A2 ) 25 22 -8 2.2 -6 9
(A 4 , A3 ) 43 20 -6 -0.3 -3 3
(A 4 , A5) 12 8 -4 -5.5 -4 2
(A 5 , A1 ) 28 -18 0 3.4 5 3
(A 5 , A2 ) 13 14 -4 7.7 -2 7
(A 5 , A3 ) 31 12 -2 5.2 1 1
(A 5 , A4) -12 -8 4 5.5 4 -2
65 - 80 -15.
90 - 58 32.
444
the H(d} functions defined in step 1. From Table 5.15, we know that
Eq. (5.163) is used for H(d 1 }; Eq.(5.164} is used for H(d 2 }; Eq. (5.166)
is used for H(d 3 }; Eq. (5.165) is used for H(d 4 ); Eq. (5.162) is used
1
Sd(1,2}
n {[1:H(d.}] + [1: H(d g } ]}
j J g
1 .296
"6 [H(d 4 ) + H(d 2 } + H(d 6 }]
d6
- 1 - exp (- 42
H(d 6 ) 1 - exp( .274
20"2 2 ( . 5) 2
A1 A2 A3 A4 A5
Sd A1 .296 .250 .285 .100
A2 .479 .396 .333 .313
A3 .242 .180 .333 .063
A4 .333 .506 .305 .224
A5 .444 .515 .570 .479
445
AI' the degree to which Al outranks others and the degree of dominance
F (1) -.567.
similarly, we can calculate all the F+, F-, and F for all the
F+ (1) -
F (1) F (1)
has the highest degree of domination over other alternatives but has
the lowest degree of dominance by other alternatives. It is also
observed that A2 and A4 are incomparable. A2 is very cheap, generates
a very small amount of electricity, and has a low degree of safety,
while A4 generates a large amount of electricity but costs more to
construct and has a high degree of safety. The ranking order is a
partial order because one cannot tell whether A2 or A4 is better. It
is given as:
.0IS4.
Assume that a = O.IS, the lower and upper bounds of the interval
Al A2 A3 A4 AS
..-----~ A2 Al
AS ---+
A4 A3
Note
can set a = O.OS level and get new intervals for Al and A3 as [-.119,
448
S~ S~ S~
.5173 .8371 .6928 .6853 .8784
Algorithm
For a MAOM problem, the OM must initially specify a set of
weights pertaining to the attributes. The discordance relation is
constructed first. The set of weights is modified iteration by
.iteration until the maximum (minimum) degree of credibility ~(~,AI)
n
{Wlw1~· .. ~w.~ ... ~w ~O, L wk 1}.
J n k=l
= max
reD kl
This equation indicates that if d kl > ~ then Ak can never outrank Ai.
l/r, if lsjsr,
{
0, if j>r.
452
,
if dkl=O,
,
if dkl>O,
max
l:sr:sn
step 5. Ask the OM: "Is Ak at least as good as AI? respond yes
or no." If yes, constraints
n
1: w.
j=l ]
n
1: w. 1 and 1: Wt:s a-c (c > 0 and small)
t=l ] eC kl
453
fuzzy multilevel graph in which the numbers associated with arc (k,l)
Note
Takeda's method shares the same disadvantage found in the other
outranking approaches.
Takeda's procedure.
454
Xl Xj x
n
A1 IJ. 1 (X 1 ) IJ. 1 (X j ) IJ. 1 (X n )
the values in the decision matrix are all given as degrees of "how one
normalization when the decision data are fuzzy. The decision data
seen in Bellman and Zadeh [BS]. Although its original intention was
associates with each alternative the gain or loss resulting from the
fuzzy set {(x'~G(X)) I xeU} where U is the universe of the fuzzy set G.
~G(x) is the membership function of the fuzzy goal and takes its
values in [0,1]. For example, the fuzzy goal expressed in words such
4 -1
(1 + (X - 15» E [0,1].
The x value that makes the highest ~G(x) value is the preferred one.
denominator for the various fuzzy goals and fuzzy constraints and
m -1
(1 + a(x - 6» E [0,1],
then at x=2 and x=10 we have approximately ~C(X) = 0.71; while at x=l
and x=ll, ~C(x) = 0.50; and at x=o and x=12, ~C(x) is about 0.32.
The above definitions of goals and constraints in a fuzzy
from the intersection of the goals and constraints. That is, given
D (5.187)
The selection of the most appropriate x value for fuzzy set D is then
given as:
For example, given the universe U {1, 2, .•• , 10} and the fuzzy
goals G1 and G2 , and fuzzy constraints C1 and C2 :
x 1 2 3 4 5 6 7 8 9 10
~G (x) 0 .1 .4 .8 1.0 .7 .4 .2 0 0
1
~G (x) .1 .6 1.0 .9 .8 .6 .5 .3 0 0
2
~C (x) .3 .6 .9 1.0 .8 .7 .5 .3 .2 .1
1
~C (x) .2 .4 .6 .7 .9 1.0 .8 .6 .4 .2
2
458
1 2 3 4 5 6 7 8 9 10
° .1 .4 .7 .8 .6 .4 .2
° °
By Eq.(5.189), we obtain
which indicates that when x = 5, the best decision has been reached.
We can conclude that the value 5 is the best selection among all
possible candidates, i. e., {1, 2, ... , 10}.
Note that none of the ~D(X), Vx, have full membership, i.e., all ~D(X)
are less than one. This indicates that a certain degree of conflict
exists among goals and constraints and none of the x values fully
Xl X2 X3 X4
[
.7 .3 .3 .5
1
.5 .8 .3 .1
.4 .6 .8 .2
where ~A (X 2 ) = .30 indicates that candidate Al satisfies the
1
experience requirement at a degree of 0.30; while ~A (X 3 ) = 0.8 shows
3
that candidate A3 'S communication ability is subjectively rated at .8.
The intersection of all attribues' values for each candidate is:
x
.3 .1 .2
Note
1. We include Bellman and Zadeh's approach in our study because
it is the pioneering work on fuzzy decision analysis. In their
approach, although the data are expressed by crisp, real numbers,
their intrinsic characteristics are fuzzy. The decision variable x in
this algorithm is treated as an alternative set. The values that x
may take are the alternatives to be evaluated. The decision problem
modeled by Bellman and Zadeh can be treated as a classical MADM
problem in that all decision data are crisp and the decision is to
select one out of several possible alternatives.
2. It has long been recognized that the applicability of the
maximin method in MADM is very limited (Hwang and Yoon [H13). The
maximin method utilizes only a small portion of the available
information in making a choice. It may happen that if an alternative
is clearly superior in all but one attribute which is below average,
another alternative which is only average on all attributes will be
chosen over it. We do not recommend this approach in any circumstance
except when the OM takes a defensive strategy, i.e., a pessimistic
attitude toward risk.
460
values, but since only single values are provided, they can also be
MADM problem.
saaty's method [Sl). When the fuzzy decision takes place, the
performance data under all attributes for each alternative are raised
is preferred.
Algorithm
Xl x. Xn
J
A1 ~1(X1) ~l(Xj) ~l(xn)
A ~m (Xl) ~m(Xj) ~m (X n )
m
Xl Xj Xn
(W 1 ) (W j ) (Wn )
Al ~11 ~1j ~ln
{wI) (W j ) (wn )
Ai ~i1 ~ij ~in (5.191)
{wI) (W j ) (wn )
Am ~m1 ~mj ~~
Wj
max [min ~ij ]. (5.192)
i j
462
Xl X2 X3 X4
[
.7 .3 .3 .5
1
.5 .8 .3 .1
.4 .6 .8 .2
Xl X2 X3 X4
Xl 1 3 7 9
1
X2 3" 1 6 7
1 1
X3 1 3
"1 6
1 1 1
X4 1
9 "1 3"
~I = (.773,.400,.106,.053)
~ (2.32,1.20,.32,.16)
Xl X2 X3 X4
[
A1 .44 .24 .60 .90
.20 .76 .68 .69
1
A2
A .12 .54 .93 .77
( ~12 ) 1.20 24
..
The first candidate, A1 , has the highest ~min value and thus will be
selected.
Note
that those which have higher membership values are reduced much less
than those which have smaller membership values. The effect of raising
the requirement of the fuzzy set, i.e., attributes that are less
x .1 .3 .5 .7 1.0
.1 .4 .6 .8 1.0
464
.1 .3 .5 .7 1.0
~:2(X) I .01 .16 .36 .64 1.0
which may be treated as a fuzzy set "very close to 1". On the other
hand, F may also be raised to the power of 1/2 to represent another
fuzzy set "sort of close to 1" as:
x .1 .3 .5 .7 1.0
~ 5 (x) .32 .63 .76 .89 1.0
F·
approach is that the MACH problem may contain fuzzy and crisp data.
Furthermore, fuzzy data may be expressed in linguistic terms or in
fuzzy numbers.
The proposed approach is composed of two major phases. The first
phase converts fuzzy data into crisp scores. The fuzzy data can be
linguistic terms, fuzzy sets, or fuzzy numbers. If the fuzzy data
are linguistic terms, they are transformed into fuzzy numbers first.
Then all the fuzzy numbers (or fuzzy sets) are assigned crisp scores.
The result of the first phase is a decision matrix which contains only
crisp data. In the second phase, classical MACH methods (see Chapter
2) can be utilized to determine the ranking order of alternatives. In
general, mathematical computations are reduced to a minimum. A SO-by-
SO decision matrix can be solved as easily as a S-by-S one. The
easy-to-use and easy-to-understand characteristics of this new
approach make it valuable to management and system analysts.
",(x)
medium high
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
",(x)
low medium high
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
",(x) very
1. low low medium
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
x
.1 .2 .3 .4 .5 .6 .7 .8 .9 1.
medium high to
P.(x~Ow to medium
high v.high
1. ,V.lOW low low medium
.9 \
\ $~
~\ ~
.8 \
\
I~
I ,
\
.7 v. ~l lv. ~
.6 lo~\ lhi~h
.5
\ t ~
~
~
.4 ~
~
.3 ~
~
.2 ~
.1 ~
l
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
P.(x)
low 10
v.low
medium
low
1
medium
medium
high
high 10
v.high
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
Scale 1 2 3 4 5 6 7 8
No. of
terms used two three five five six seven nine eleven
none yes
v. low yes yes yes yes yes
low-v. low yes yes
low yes yes yes yes yes yes yes
Example
The concept behind our system is to match the linguistic terms
under each attribute with one of the conversion scales, assuming that
the OM is not available for consulting (a non interactive situation).
Once a scale is found which matches all the linguistic terms used, it
will be employed to convert linguistic terms to their corresponding
fuzzy numbers. If more than one scale is found to match all the
linguistic terms used, the scale with the least number of terms is
used for conversions. The following example will illustrate our
concept and procedure.
Four commercial nonlinear programming software packages are
evaluated with respect to their ease of use, cost, editing facilities,
problem storage and retrieval, and supporting documents. The decision
matrix (modified from Golden and Wasil [GG]) is given as:
473
Jain's [J2,J5] and Chen's [C12] fuzzy ranking approaches. The crisp
X,
{ (5.193)
0, otherwise
1 - x, 0 :s X :s 1
{ (5.194)
0, otherwise
(5.195)
(5.196)
Given the left and right scores of M, we can compute the total score
of Musing:
(5.197)
P.(x)
Ml M2
1.
.9
.8
.5
.4
.3
.2
.1 .' .
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
{
x, o " x " 1
/.lmax(X)
0, otherwise
and
{ o " x "
I - x, 1
/.lmin(X)
0, otherwise
x - .4
, 0.4 " x < 0.6
.2
1
/.lM (x)
1
.8 - x , 0.6 x 0.8
.2 " "
x - .6
.2
, 0.6
" x < 0.8
1
/.lM (x)
2 I - X
, 0.8 x 1.0
.2 " "
[~R(M1) + 1 - ~L(M1)]
~T(M1) .584.
2
one scale is identified, the one with the fewest linguistic terms will
with fuzzy elements into one with real numbers. It can be solved
using any classical MADM method (see Chapter 2). We have chosen to
use the TOPSIS method for the following numerical examples because of
(Hwang and Yoon [H13]) will give cardinal order of the alternatives.
477
surrounded by mountains, and covers about 30,000 sq. km. The main
where v.good = very good, exl. excellent, v.easy very easy, and
diff. = difficult.
Assume the weights for all attributes are equal. The problem is
step 1. We have identified that attributes X3 ' Xs ' X6 ' X8 ' X9 '
X10 ' and X11 contain linguistic terms. These will be converted to
fuzzy numbers.
match the linguistic terms with Scale 8. For X8 ' Scale 3 is used
again, while Scale 5 is used for the 10th attribute. The results are
shown in Figs. S.47, 5.48, 5.49, and 5.S0.
x+
1
x+
2 X3
• X4 Xs
• •
X6 X7 •
X8
A1 99.6 4 .909 .7 .909 .667 90 .909
A2 85.7 19 .717 .5 .717 .954 80 .909
A3 101.1 SO .283 .01 .500 .500 80 .717
A4 9S.1 50 .909 .1 .283 .954 60 .500
AS 101.8 SO .SOO .01 .283 .333 70 .500
X9* X10 • X*
ll
.909 .917 .909 A1
.717 .7S0 .717 A2
.283 .416 .500 A3
.717 .2S0 .283 A4
.500 .416 .500 AS
where attributes with a "+" sign are cost attributes (the smaller the
,.,.(x) very
poor fair good good
1. ,
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
,.,.(x)
poor fair good exel
1.
,
.9
.8
.7
.6
.5
.4
.3
.2
.1
0
0 .1 .2 .3 .4 .5 .6 .7 .8 .9
p,(x) very
fair good good
1.
.9
.8
.7
.6
.5
.4
.3
.2
.1
0 x
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.
x+
1
x+
2 X3* X4 X5* X6*
A1 0.4599 0.0451 0.5753 0.8082 0.6870 0.4116
A2 0.3958 0.2141 0.4538 0.5773 0.5419 0.5887
A3 0.4669 0.5634 0.1791 0.0115 0.3779 0.3085
A4 0.4392 0.5634 0.5753 0.1155 0.2139 0.5887
A5 0.4701 0.5634 0.3164 0.0115 0.2139 0.2055
X7 x8* X9* *
X10 *
X11
0.5249 0.5566 0.6150 0.6812 0.6559 A1
0.4666 0.5566 0.4851 0.5572 0.5174 A2
0.4666 0.4391 0.1915 0.3090 0.3608 A3
1
0.3499 0.3062 0.4851 0.1857 0.2042 A4
0.4082 0.3062 0.3383 0.3090 0.3608 A5
formula
482
11 2
{ ~ ( ) }1/2, ~ 1,2,3,4,5,
LVi]' - v]' ...
j=1
i 1 2 345
0.1884 0.4154 1.2890 1.2386 1.3134
11 _ 2
s, { ~
LVi]'
(
- v]'
) } 1/2,
~-
j=1
i 1 2 345
s, 1. 4314 1. 0845 0.3294 0.6339 0.2890
~-
1,2,3,4,5,
We can obtain
i 1 2 345
0.8837 0.7230 0.2035 0.3385 0.1804
following pitfalls.
1. Size of Problem.
attributes. And the alternatives evaluated are less than 10. While
the algorithms of these methods may be valid, real world MADM problems
seldom come in such small dimensions. As a result, existing methods
cannot effectively (from labor and cost aspects) solve real world
problems.
The first example is the selection of the most promising
admission problem summarized in Table 5.23. Notice that there are only
up to several thousand.
484
The majority of the fuzzy MADM methods assume the decision data
such data should not be forced into a fuzzy format because it will
Attributes
Attribute A6 A7 AS
Xl X2 X3
A1 .5 .5 .2
.7 .4 .01
[
A2
A3 .3 .8 .6
A4 .6 .4 .9
1
In this problem, Xl = candidate should be young, X2 = should be
experienced, X3 = be able to communicate. The datum may be explained
as: "Candidate Ai satisfies attributes Xj at degree of a e [0,1]".
Theoretically, the data is fuzzy. But from an operational aspect, the
decision matrix is no more than a classical MADM model which can be
easily solved by an appropriate classical MACM method.
VI. CONCLUDING REMARKS
Xl X2
A1 x 11 x 12
x 21 x 22 (6.1)
[
A2
D
~ xm1 xm2 x~ 1
where Ai represents an alternative, Xj denotes an attribute, and
x ij is the performance rating of alternative Ai with respect to
attribute Xj • x ij is also referred to as the decision data.
In most of the real-world problems, some of the decision data x ij
can be precisely assessed while others cannot. We use real (crisp)
numbers to represent data which can be precisely measured. For those
data which cannot be precisely assessed, we use Zadeh's [Zl] fuzzy
sets (numbers) to denote them. The use of fuzzy set theory allows us
to incorporate unquantifiable information, incomplete information,
nonobtainable information, and partially ignorant facts into the
decision model.
decision model. When decision data are precisely known, they should
they are transformed into fuzzy numbers first. Then by using fuzzy
decision matrix are assigned crisp scores. The result of the first
The existing works on this topic are few and seem quite arbitrary. We
4~
The system contains eight scales. There are generic verbal terms
(ranging from 2 to 11) in the system where Scale 1 contains only two
After all the linguistic terms have been converted to fuzzy sets
(numbers), a decision matrix which contains either fuzzy numbers, real
to numbers E [0,1) which are called the fuzzy scores. This is done
First, a total score is determined for each fuzzy number. The fuzzy
scoring system, the total scores are not compared, rather, they
problems and the majority of them are capable of handling large size
MADM problems. Our new approach has successfully extended that
ability to the fuzzy problem domain.
2. Our approach allows MADM problems to take data in the forms of
Creative confrontation
and stimulating ideas
Expert judgement!
group participation Systematic structuring
~
Simulation
Implementation and
controlling
The voting process is carried out by all the voters who may be in a
the general run of electors. On the other hand, the primary concern
effectiveness.
Let us use an example to illustrate the voting and counting
following manner:
Condorcet effect.
How do we solve this Condorect effect problem? It is a vital
social choice theory. The question is: What kind of decisions are
ways: (1) the social choice function, and (2) the social welfare
functions.
judgment/group participation.
1.1 Nonranked Voting System Condorcet Principle 3.1 Arrow's Conditions for Social Welfare Function
1.1.1 One Member Elected From Two Candidates 2.1 Condorcet's Function 3.2 Arrow's Possibility Theorem for Two Alternatives
1.1.2 One member Elected From Multi-Candidates .
(a) The First Past the Post System 2.2 Borda's Function 3.3 Arrow's General Possibility Theorem
(b) Majority Representation System
Repeated Ballots 2.3 Copeland's Function 3.4 Black's Single-Peaked Preferences
The Second Ballods
2.4 Nanson's Function 3.5 Bowman and Colantoni's Approach
1.1.3 Election of Two or More Members
1.1.3.1 The Single Non-Transferable Vote 2.5 Dodgson's Function 3.6 Goodman and Markowitz's Approach
1.1.3.2 Multiple Vote
1.1.3.3 Limited Vote 2.6 Kemeny's Function 3.7 Cardinal Social Welfare Function
1.1.3.4 Cumulative Vote
1.1.3.5 List Systems 2.7 Cook and Seiford's Function 3.7.1 Value Function for Certainty Case
(a) Highest Average $
(b) Greatest Remainder 2.8 Fishburn's Function 3.7.2 Utility Function for Uncertainty Case
1.1.3.6 Approval Vote
2.9 Eigenvector Function 3.7.2.1 Additive Group Utility Function
1.2 Preferential Voting System
2.10 Bernardo's Assignment Approach 3.7.2.2 Multiplicative Group Utility Function
1.2.1 Simple Majority Decision Rule
1.2.1.1 Two-Alternative Case 2.11 Cook and Seiford's Ordinal 3.7.3.1 Gymnastics Competitions
1.2.1.2 More Than Two Alternatives Case Intersection Method
(a) Paradox of Voting 3.7.3.2 Extended Contribution Rule (ECR) Method
(b) The Condorcet Effect
Fig. 6.3 Approaches in Social Choice Theory [Hwang and Lin, H14].
Phases Extraction Exploration Selection Execution
/",/,,~t~t
Idea Issue Problem Problem Planning
Stimulation Clarification Structuring Solving Solution
Activities
t :> <: t ::::>c< t>< t ~Ietentlng
Creative Polling of Systematic Simulation and
Confrontation Experts/participant Structuring A Controlling
three basic abstract forms and approaches of game theory are shown in
Fig. 6.5.
Core concept
Characteristic
Function form of
game
-E Shapley value
Parametrical approach
linguistic term, fuzzy ranking methods, fuzzy MADM methods, and MADM
manner. This system has not been tested vigorously, nor theoretically
following factors: (1) the shape, (2) the height, (3) relative
location on x-axis and y-axis, (4) the spread, (5) the centroid point,
may be attributed to the fact that they use only one or two factors to
ranking methods are not perfect, they have shown the process of human
efforts to find ways to solve problems. There are always some good
MADM problem. The many factors that characterize fuzzy numbers may be
solve this type of MADM problem is not difficult. But to identify the
is obvious.
503
the number m can be 5, 10, 50, 100, 1000, or over thousands, and the
in Fig. 2.1a). The decision data x ij can take fuzzy or crisp data as
simply apply the method that they are familiar with (or willing to
use). This practice will often result in an ad hoc decision.
which states that we must let the problems which we are trying to
solve determine the methods we need to use rather than having the
5~
have been conducted during the past decade. For example, Hwang and
entire armory of both the classical and fuzzy MADM methods. The next
Inference rules must be built to match the problem types and the
decision methods.
VII. BIBLIOGRAPHY
A1. Adamo, J.M., Fuzzy decision trees, Fuzzy Sets and Systems,
Vol.4, No.3, 1980, pp. 207-220.
A4. Alley, H.A., C.P. Bacinello, and K.W. Hipel, Fuzzy set
approaches to planning in the grand river basin,
Advances in water Resources, Vol. 2, 1979, pp. 3-12.
B16. Blin, J.M., and A.B. Whinston, Fuzzy sets and social
choice, Journal of Cybernetics, Vol. 3, 1974, pp. 28-36.
B17. Blin, J.M., K.S. FU, A.B. Whinston. and K.B. Moberg,
Pattern recognition in micro-economics, Journal of
Cybernetics, Vol. 3, 1974, pp. 17-27.
B22. Bolanos, M., M.T. Lamata, and S. Moral, Decision making problems
under a Shaffer's evidence, First International Fuzzy System
Association (IFSA) Congress, July 1985.
B30. Borisov, A.N., Y.Y. Luns., V.A. Popov, Fuzzy decision analysis
in the tasks of electric network development, Proc. IFAC
symposium on Fuzzy Information, Knowledge Representat~on and
Decision Analys~s, Marse~lle, Pergamon Press, New York, 1983,
pp. 289 294.
C11. Chang, W.K. A Study on the Ranking of Fuzzy Alternatives and its
Apelication to Dec1sion Making, Ph.D. dissertation, Tamkang
un1versity, Ta~an, 1982.
04. Delgado, M., J.L. Verdegay, and M.A. Villa, A procedure for
ranking fuzzy numbers using fuzzy relations, Fuzzy Sets and
Systems, Vol. 26, No.1, 1988.
07. Deutsch, S.J. and C.J. Malborg, A fuzzy set approach to data
set evaluation for decision support, IEEE Trans On Systems,
Man, and Cybernetics, Vol. SMC-15, No.6, 1985, pp. 777-783.
026. Dubois, D., and H. Prade, Fuzzy Sets and Systems: Theory
and Applications, Academic Press, New York, 1980.
028. Dubois, D., and H. Prade, The use of fuzzy numbers in decision
analysis, In: Fuzzy Information and Decision Processes, M.M.
Gupta and E. Sanchez (eds.), North Holland, 1982, pp.309-321.
513
E8. Enta, Y., Fuzzy decision theory, In: Fuzzy Set and Possibily
Theory - Recent Developments, R.R. Yager (ed.), Pergamon
Press, New York, 1982, pp. 439-449.
E13. Esogbue, A.O., and R.C. Elder, Fuzzy sets and the modeling
of physician decision process, part II: Fuzzy diagnosis
decision models, Fuzzy Sets and Systems, Vol, 3, No.1,
1980, pp. 1-9.
E15. Esogbue, A.O., and Z.M. Ahipo, Fuzzy sets and water
resources planning, In: Fuzzy Set and Possibility Theory
Recent Development, R.R. Yager(ed.), Pergamon Press, New
York, 1982, pp. 450-465.
515
G9. Gupta, M.M., and E.H. Mamdani, Second IFAC round table on fuzzy
automata and decision processes, Automatica, Vol. 12, 1976,
pp.291-296.
G10. Gupta, M.M., G.N. Saridies, and B.R. Gaines (eds.), Fuzzy
Automata and Decision Processes, North-Holland, Amsterdam, 1977.
G11. Gupta, M.M., R.K. Ragade, and R.R. Yager (eds.), Advances in
Fuzzy Set Theory and Applications, North-Holland, New York,
1979.
G16. Gusev, L.A., and I.M. Smirnova, Fuzzy sets - Theory and
applications (A survey), Automation and Remote Control, No.5,
May 1973, pp. 66-85.
H14. Hwang, C.L., and M.J. Lin, Group Decision Making Under
Multiple Criteria, springer-Verlag, New York, 1987.
K19a. Kim, K., and K.S. Park, Ranking fuzzy numbers with index of
optimism, Fuzzy Sets and Systems, Vol. 35, 1990, pp. 143-150.
K23. Klir, G.J., and T.A. Folger, Fuzzy Sets, Uncertainly, and
Information, Prentice Hall, Englewood cliffs, New Jersey, 1988.
K31. Kuzmin, V.B., and S.I. Travkin, Fuzzy choice, In: optimization
Models Using Fuzzy Sets and Possibility Theory, J. Kacprzyk
(eds.), D. Reidel Publishing Co., Boston, 1987, pp. 99-121.
L3. Lee, E.S., and R.L. Li, Comparison of fuzzy numbers based on
the probability measure of fuzzy events, Computer and
Mathematics with Applications, Vol. 15, 1988, pp. 887-896.
521
L7. Li, R.J., and E.S. Lee, Ranking fuzzy numbers - A comparison,
Proceedings of the North American Fuzzy Information Processing
society Workshop, J.L. Chameau and James T.P. Yao (eds.),
NAFIPS 1987, Purdue university 1987, pp. 169-204.
M5. Mamdani, E.H., and B.R. Gaines (eds.), Fuzzy Reasoning and
its Applications, Academic Press, London, 1981.
M16. Miller, G.A., The magic number seven, plus or minus seven,
Psychological Review, Vol.63, 1965, pp 81-97.
M20. Mizumoto, M., Fuzzy sets and their operations, Part II,
Information and Control, Vol. 50, 1981, pp. 160-174.
N5. Nauta Lemke, H.R. van, T.G. Dijkman, H. van Haeringen, and M.
Pleeging, A characteristic optimism factor in fuzzy decision-
making, Proceedings of IFAC Symposium on Fuzzy Information,
Knowledge Representation and Decislon Analysls, Pergamon Press,
New York, 1983, pp. 283-288.
04. Osgood, C.E., G.J. Suci, and P.H. Tannenbaum, The Measurement
of Meaning (9th ed.), University of Illinois Press, Urbana,
1975.
R13. Roy, B., Selektieren, sortieren und orden mit hilfe von
Pravalenzrelationen, zeitschr. F. Betriebsw. Forschung,
Vol.32, 1980, pp. 465-496.
S17. Simon, H.A., Two heads are better than one: The collaboration
between AI and OR, Interfaces, Vol. 17, No.4, 1987, pp. 8-15.
S19. siskos, J., and Ph. Hubert, A survey and a new comparative
approach, European Journal of operation Research, Vol. 13,
1983, pp. 278-299.
S20. siskos, J.L., J. Lochard, and J. Lombard, A multicriteria
decision making methodology under fuzziness: Application to the
evaluation of radiological protection in nuclear power plants,
In: TIMS/studies in the Management sciences, Vol. 20,
H.J. Zimmermann (eds.), Elsevier Science Publishers,
North-Holland, 1984, pp. 261-283.
T4. Tanaka, H.T., T. Okuda, and K. Asai, Decision making and its
goal in a fuzzy environment, Proceedings of US-Japan Seminar On
Fuzzy Sets and their Applicat10ns, Berkeley, CA, July 1974.
V1. Van Gigch, J.P., and L.L. Pipino, From absolute to probable
and fuzzy in decision-making, Kybernetes, Vol. 9, No.1, 1980.
W1a. Wang, P.P., and S.K. Chang, Fuzzy Sets - Theory and Applications
to Policy Analysis and Information Systems, Plenum Press,
New York, 1980.
W2. Watada, J., and H. Tanaka, The perspective of possibilis-
tic models in decision making, Proceedings of the 7th
International Conference on Multiple criteria Decision
Making, Kyoto, Japan, 1986, pp. 328-337.
W3. Watson, S.R., J.J. Weiss, and M.L. Donnell, Fuzzy decision
analysis, IEEE Trans. On Systems, Man, and cybernetics,
SMC-9, 1979, pp. 1-9.
Y6. Yager, R.R., Ranking fuzzy subsets over the unit interval,
Proc. 1978 CDC, 1978, pp.1435-1437.
Y18. Yager, R.R., A procedure for ordering fuzzy subsets of the unit
interval, Information Sciences, Vol. 24, 1981, pp. 143-161.
Y31a. Ying, M., The alternative measures of fuzzy sets, Fuzzy Sets
and Systems, Vol. 37, 1990, pp. 105-110.
Y32. Yoon, K., System Selection by Multiple Attribute Decision
Making, Ph.D. Thesis, Dept. of Industrial Engineering,
Kansas State University, 1980.
Y33. Yoon, K., A reconciliation among discrete compromise solutions,
Journal of Operational Research Society, Vol. 38, 1987,
pp.277-286.
Vol. 211: P. van den Heuve1, The Stability of a Macroeconomic Vol. 233: G. Wagenhals, The World Copper Market. XI, 190
System with Quantity Constraints. VII, 169 pages. 1983. pages. 1984.
Vol. 212: R. Sato and T. Nono, Invariance Principles and the Vol. 234: B.C. Eaves, A Course in Triangulations for Solving
Structure of Technology. V, 94 pages. 1983. Equations with Deformations. III, 302 pages. 1984.
Vol. 213: Aspiration Levels in Bargaining and Economic Vol. 235: Stochastic Models in Reliability Theory Proceedings,
Decision Making. Proceedings, 1982. Edited by R. Tietz. VIII, 1984. Edited by S. Osaki and Y. Hatoyama. VII, 212 pages.
406 pages. 1983. 1984.
Vol. 214: M. Faber, H. Niemes und G. Stephan, Entropie, Um- Vol. 236: G. Gandolfo, P.C. Padoan, A DiseqUilibrium Model
we1tschutz und Rohstoffverbrauch. IX, 181 Seiten. 1983. of Real and Financial Accumulation in an Open Economy. VI,
Vol. 215: Semi-Infinite Programming and Applications. Pro- 172 pages. 1984.
ceedings, 1981. Edited by A. V. Fiacco and K.O. Kortanek. XI, Vol. 237: Misspecification Analysis. Proceedings, 1983. Edited
322 pages. 1983. by T.K. Dijkstra. V, 129 pages. 1984.
Vol. 216: H.H. Miiller, Fiscal Policies in a General Equilibrium Vol. 238: W. Domschke, A. Drexl, Location and Layout
Model with Persistent Unemployment. VI, 92 pages. 1983. Planning. IV, 134 pages. 1985.
Vol. 217: Ch. Grootaert, The Relation Between Final Demand Vol. 239: Microeconomic Models of Housing Markets. Edited
and Income Distribution. XIV, 105 pages. 1983. by K. Stahl. VII, 197 pages. 1985.
Vol 218: P. van Loon, A Dynamic Theory of the Firm: Vol. 240: Contributions to Operations Research. Proceedings,
Production, Finance and Investment. VII, 191 pages. 1983. 1984. Edited by K. Neumann and D. Pallaschke. V, 190 pages.
Vol. 219: E. van Damme, Refinements of the Nash Equilibrium 1985.
Concept. VI. 151 pages. 1983. Vol. 241: U. Willmann, Das Konzept rationaler Preiser-
Vol. 220: M. Aoki, Notes on Economic Time Series Analysis: wartungen. XI, 310 Seiten. 1985.
System Theoretic Perspectives. IX, 249 pages. 1983. Vol. 242: Decision Making with Multiple Objectives.
Vol. 221: S. Nakamura, An Inter-Industry Trans10g Model of Proceedings, 1984. Edited by Y.Y. Haimes and V. Chankong.
Prices and Technical Change for the West German Economy. XI, 571 pages. 1985.
XIV, 290 pages. 1984. Vol. 243: Integer Programming and Related Areas. A Classified
Vol. 222: P. Meier, Energy Systems Analysis for Developing Bibliography 1981-1984. Edited by R. von Randow. XX, 386
Countries. VI, 344 pages. 1984. pages. 1985.
Vol. 223: W. Trockel, Market Demand. VIII, 205 pages. 1984. Vol. 244: Advances in Equilibrium Theory. Proceedings, 1984.
Edited by C.D. Aliprantis, O. Burkinshaw and N.J. Rothman.
Vol. 224: M. Kiy, Ein disaggregiertes Prognosesystem fur die 11,235 pages. 1985.
Bundesrepublik Deutschland. XVIII, 276 Seiten. 1984.
Vol. 245: J.E.M. Wilhelm, Arbitrage Theory. VII, 114 pages.
Vol. 225: T.R. von Ungern-Sternberg, Zur Analyse von Miirk- 1985.
ten mit unvollstiindiger Nachfragerinformaton. IX, 125 Seiten.
Vol. 246: P.W. Oller, Dynamic Feature Space Modelling,
1984.
Filtering and Self-Tuning Control of Stochastic Systems. XIV,
Vol. 226: Selected Topics in Operations Research and Math- 177 pages. 1985.
ematical Economics. Proceedings, 1963. Edited by G. Hammer
and D. Pallaschke IX, 478 pages. 1984. Vol. 247: Optimization and Discrete Choice in Urban Systems.
Proceedings, 1983. Edited by B.G. Hutchinson, P. Nijkamp and
Vol. 227: Risk and Capital. Proceedings, 1983. Edited by G. M. Bally VI, 371 pages. 1985.
Bamberg and K. Spremann VII, 306 pages. 1984.
Vol. 248: Pural Rationality and Interactive Decision Processes.
Vol. 228: Nonlinear Models ofFluclUating Growth. Proceedings, Proceedings, 1984. Edited by M. Grauer, M. Thompson and A.P.
1983. Edited by R.M. Goodwin, M. Kriiger and A. Vercelli. Wierzbicki. VI, 354 pages. 1985.
XVII, 277 pages. 1984.
Vol. 249: Spatial Price Equilibrium: Advances in Theory,
Vol. 229: Interactive Decision Analysis. Proceedings, 1983. Ed- Computation and Application. Proceedings, 1984. Edited by P.T.
ited by M. Grauer and A.P. Wierzbicki. VIII, 269 pages. 1984. Harker. VII, 277 pages. 1985.
Vol. 230: Macro-Economic Planning with Conflicting Goals. Vol. 250: M. Roubens, Ph. Vincke, Preference Modelling. VIII,
Proceedings, 1982. Edited by M. Despontin, P. Nijkamp and J. 94 pages. 1985.
Spronk. VI, 297 pages. 1984.
Vol. 251: Input-Output Modeling. Proceedings, 1984. Edited
Vol. 231: G.F. Newell, The M/M/8 Service System with Ranked by A. Smyshlyaev. VI, 261 pages. 1985.
Servers in Heavy Traffic. XI, 126 pages. 1984.
Vol. 252: A. Birolini, On the Use of Stochastic Processes in
Vol. 232: L. Bauwens, Bayesian Full Information Analysis of Modeling Reliability Problems. VI, 105 pages. 1985.
Simultaneous Equation Models Using Integration by Monte
Carlo. VI, 114 pages. 1984.
Vol. 253: C. Withagen, Economic Theory and International Vol. 281: Ch.-L. Hwang, M.-J. Lin, Group Decision Making
Trade in Natural Exhaustible Resources. VI, 172 pages. 1985. under Multiple Criteria. XI, 400 pages. 1987.
Vol. 254: S. Miiller, Arbitrage Pricing of Contingent Claims. Vol. 282: K. Schittkowski, More Test Examples for Nonlinear
VIII, 151 pages. 1985. Programming Codes. V, 261 pages. 1987.
Vol. 255: Nondifferentiable Optimization: Motivations and Ap- Vol. 283: G. Gabisch, H.-W. Lorenz, Business Cycle Theory.
plications. Proceedings, 1984. Edited by V.F. Demyanov and VII, 229 pages. 1987.
D. Pallaschke. VI, 350 pages. 1985. Vol. 284: H. Liitkepohl. Forecasting Aggregated Vector ARMA
Vol. 256: Convexity and Duality in Optimization. Proceedings, Processes. X, 323 pages. 1987.
1984. Edited by J. Ponstein. V, 142 pages. 1985. Vol. 285: Toward Interactive and Intelligent Decision Support
Vol. 257: Dynamics of Macrosystems. Proceedings, 1984. Edited Systems. Volume 1. Proceedings, 1986. Edited by Y. Sawaragi,
by J.-P. Aubin, D. Saari and K. Sigmund. VI, 280 pages. 1985. K. Inoue and H. Nakayama. XII, 445 pages. 1987.
Vol. 258: H. Funke, Eine allgemeine Theorie der Polypol- und Vol. 286: Toward Interactive and Intelligent Decision Support
Oligopolpreisbildung. III, 237 pages. 1985. Systems. Volume 2. Proceedings, 1986. Edited by Y. Sawaragi,
Vol. 259: Infinite Programming. Proceedings, 1984. Edited by K. Inoue and H. Nakayama. XII, 450 pages. 1987.
E.J. Anderson and A.B. Philpott. XIV, 244 pages. 1985. Vol. 287: Dynamical Systems. Proceedings, 1985. Edited by
Vol. 260: H.-J. Kruse, Degeneracy Graphs and the A.B. Kurzhanski and K. Sigmund. VI, 215 pages. 1987.
Neighbourhood Problem. VIII, 128 pages. 1986. Vol. 288: G.D. Rudebusch, The Estimation of Macroeconomic
Vol. 261: Th.R. Gulledge, Jr., N.K. Womer, The Economics of Disequilibrium Models with Regime Classification Informati-
Made-to-Order Production. VI, 134 pages. 1986. on. VII,128 pages. 1987.
Vol. 262: H.U. Buhl, A Neo-Classical Theory of Distribution Vol. 289: B.R. Meijboom, Planning in Decentralized Firms. X,
and Wealth. V, 146 pages. 1986. 168 pages. 1987.
Vol. 263: M. Schafer, Resource Extraction and Market Vol. 290: D.A. Carlson, A. Haurie, Infinite Horizon Optimal
Struucture. XI, 154 pages. 1986. Control. XI, 254 pages. 1987.
Vol. 264: Models of Economic Dynamics. Proceedings, 1983. Vol. 291: N. Takahashi, Design of Adaptive Organizations. VI,
Edited by H.F. Sonnenschein. VII, 212 pages. 1986. 140 pages. 1987.
Vol. 265: Dynamic Games and Applications in Economics. Vol. 292: I. Tchijov, L. Tomaszewicz (Eds.), Input-Output
Edited by T. Basar. IX, 288 pages. 1986. Modeling. Proceedings, 1985. VI, 195 pages. 1987.
Vol. 266: Multi-Stage Production Planning and Inventory Vol. 293: D. Batten, J. Casti, B. Johansson (Eds.), Economic
Control. Edited by S. Axsater, Ch. Schneeweiss and E. Silver. Evolution and Structural Adjustment. Proceedings, 1985. VI,
V, 264 pages.1986. 382 pages.
Vol. 267: R. Bemelmans, The Capacity Aspect of Inventories. Vol. 294: J. Jahn, W. Knabs (Eds.), Recent Advances and
IX, 165 pages. 1986. Historical Development of Vector Optimization. VII, 405 pages.
1987.
Vol. 268: V. Firchau, Information Evaluation in Capital Markets.
VII, 103 pages. 1986. Vol. 295. H. Meister, The Purification Problem for Constrained
Games with Incomplete Information. X, 127 pages. 1987.
Vol. 269: A. Borglin, H. Keiding, Optimality in Infinite Horizon
Economies. VI, 180 pages. 1986. Vol. 296: A. Borsch-Supan, Econometric Analysis of Discrete
Choice. VIII, 211 pages. 1987.
Vol. 270: Technological Change, Employment and Spatial Dy-
namics. Proceedings, 1985. Edited by P. Nijkamp. VII, 466 Vol. 297: V. Fedorov, H. Lauter (Eds.), Model-Oriented Data
pages. 1986. Analysis. Proceedings, 1987. VI, 239 pages. 1988.
Vol. 271: C. Hildreth, The Cowles Commission in Chicago, Vol. 298: S.H. Chew, Q. Zheng, Integral Global Optimization.
1939-1955. V, 176 pages. 1986. VII, 179 pages. 1988.
Vol. 272: G. Clemenz, Credit Markets with Asymmetric Infor- Vol. 299: K. Marti, Descent Directions and Efficient Solutions
mation. VIII,212 pages. 1986. in Discretely Distributed Stochastic Programs. XIV, 178 pages.
1988.
Vol. 273: Large-Scale Modelling and Interactive Decision
Analysis. Proceedings, 1985. Edited by G. Fandel, M. Grauer, Vol. 300: U. Derigs, Programming in Networks and Graphs.
A. Kurzhanski and A.P. Wierzbicki. VII, 363 pages. 1986. XI, 315 pages. 1988.
Vol. 274: W.K. Klein Haneveld, Duality in Stochastic Linear Vol. 301: J. Kacprzyk, M. Roubens (Eds.), Non-Conventional
and Dynamic Programming. VII, 295 pages. 1986. Preference Relations in Decision Making. VII, 155 pages. 1988.
Vol. 275: Competition, Instability, and Nonlinear Cycles. Pro- Vol. 302: H.A. Eiselt, G. Pederzoli (Eds.), Advances in
ceedings, 1985. Edited by W. Semmler. XII, 340 pages. 1986. Optimization and Control. Proceedings, 1986. VIII, 372 pages.
1988.
Vol. 276: M.R. Baye, D.A. Black, Consumer Behavior, Cost of
Living Measures, and the Income Tax. VII, 119 pages. 1986. Vol. 303: F.X. Diebold, Empirical Modeling of Exchange Rate
Dynamics. VII, 143 pages. 1988.
Vol. 277: Studies in Austrian Capital Theory, Investment and
Time. Edited by M. Faber. VI, 317 pages. 1986. Vol. 304: A. Kurzhanski, K. Neumann, D. Pallaschke (Eds.),
Optimization, Parallel Processing and Applications.
Vol. 278: W.E. Diewert, The Measurement of the Economic Proceedings, 1987. VI, 292 pages. 1988.
Benefits of Infrastructure Services. V, 202 pages. 1986.
Vol. 305: G.-J.C.Th. van Schijndel, Dynamic Firm and Investor
Vol. 279: H.-J. Biittler, G. Frei and B. Schips, Estimation of Behaviour under Progressive Personal Taxation. X, 215
Disequilibrium Modes. VI, 114 pages. 1986. pages. 1988.
Vol. 280: H.T. Lau, Combinatorial Heuristic Algorithms with Vol. 306: Ch. Klein, A Static Microeconomic Model of Pure
FORTRAN. VII, 126 pages. 1986. Competition. VIII, 139 pages. 1988.
Vol. 307: T.K. Dijkstra (Ed.), On Model Uncertainty and its Vol. 333: N. Dellaert, Production to Order. VII, 158 pages. 1989.
Statistical Implications. VII, 138 pages. 1988. Vol. 334: H.-W. Lorenz, Nonlinear Dynamical Economics and
Vol. 308: J.R. Daduna, A. Wren (Eds.), Computer-Aided Tran- Chaotic Motion. XI, 248 pages. 1989.
sit Scheduling. VIII, 339 pages. 1988.
Vol. 335: A.G. Lockett, G. Islei (Eds.), Improving Decision
Vol. 309: G. Ricci, K. Velupillai (Eds.). Growth Cycles and Making in Organisations. Proceedings. IX, 606 pages. 1989.
Multisectoral Economics: the Goodwin Tradition. III, 126 pages.
Vol. 336: T. Puu, Nonlinear Economic Dynamics. VII, 119
1988. pages. 1989.
Vol. 310: J. Kacprzyk, M. Fedrizzi (Eds.), Combining Fuzzy Vol. 337: A. Lewandowski, I. Stanchev (Eds.), Methodology
Imprecision with Probabilistic Uncertainty in Decision Making. and Software for Interactive Decision Support. VIII, 309 pages.
IX, 399 pages. 1988. 1989.
Vol. 311: R. Fare, Fundamentals of Production Theory. IX, 163
Vol. 338: J.K. Ho, R.P. Sundarraj, DECOMP: an Implementation
pages. 1988. of Dantzig-Wolfe Decomposition for Linear Programming. VI,
Vol. 312: J. Krishnakumar, Estimation of Simultaneous Equation 206 pages.
Models with Error Components Structure. X, 357 pages. 1988.
Vol. 339: J. Terceiro Lomba, Estimation of Dynamic
Vol. 313: W. Jammernegg. Sequential Binary Investment Deci- Econometric Models with Errors in Variables. VIII, 116 pages.
sions. VI. 156 pages. 1988. 1990.
Vol. 314: R. Tietz. W. Albers. R. Selten (Eds.). Bounded Ratio- Vol. 340: T. Vasko, R. Ayres, L. Fontvieille (Eds.), Life Cycles
nal Behavior in Experimental Games and Markets. VI. 368 and Long Waves. XIV, 293 pages. 1990.
pages. 1988. Vol. 341: G.R. Uhlich, Descriptive Theories of Bargaining. IX,
Vol. 315: I. Orishimo, G.J.D. Hewings, P. Nijkamp (Eds), In- 165 pages. 1990.
formation Technology: Social and Spatial Perspectives. Vol. 342: K. Okuguchi, F. Szidarovszky, The Theory of
Proceedings 1986. VI. 268 pages. 1988. Oligopoly with Multi-Product Firms. V, 167 pages. 1990.
Vol. 316: R.L. Basmann, D.J. Slottje, K. Hayes. J.D. Johnson, Vol. 343: C. Chiarella, The Elements of a Nonlinear Theory of
D.J. Molina, The Generalized Fechner-Thurstone Direct Utility Economic Dynamics. IX, 149 pages. 1990.
Function and Some of its Uses. VIII, 159 pages. 1988.
Vol. 344: K. Neumann, Stochastic Project Networks. XI, 237
Vol. 317: L. Bianco, A. La Bella (Eds.), Freight Transport pages. 1990.
Planning and Logistics. Proceedings, 1987. X, 568 pages. 1988.
Vol. 345: A. Cambini, E. Castagnoli, L. Martein, P Mazzoleni,
Vol. 318: T. Doup, Simplicial Algorithms on the Simplotope. S. Schaible (Eds.), Generalized Convexity and Fractional
VIII. 262 pages. 1988. Programming with Economic Applications. Proceedings, 1988.
Vol. 319: D.T. Luc. Theory of Vector Optimization. VIII, 173 VII, 361 pages. 1990.
pages. 1989. Vol. 346: R. von Randow (Ed.), Integer Programming and
Vol. 320: D. van der Wijst, Financial Structure in Small Related Areas. A Classified Bibliography 1984-1987. XIII, 514
Business. VII, 181 pages. 1989. pages. 1990.
Vol. 321: M. Di Matteo, R.M. Goodwin. A. Vercelli (Eds.). Vol. 347: D. Rios Insua, Sensitivity Analysis in Multi-objective
Technological and Social Factors in Long Term Fluctuations. Decision Making. XI, 193 pages. 1990.
Proceedings. IX, 442 pages. 1989. Vol. 348: H. Stormer, Binary Functions and their Applications.
Vol. 322: T. Kollintzas (Ed.), The Rational Expectations VIII, 151 pages. 1990.
Equilibrium Inventory Model. XI, 269 pages. 1989. Vol. 349: G.A. Pfann, Dynamic Modelling of Stochastic Demand
Vol. 323: M.B.M. de Koster. Capacity Oriented Analysis and for Manufacturing Employment. VI, 158 pages. 1990.
Design of Production Systems. XII, 245 pages. 1989. Vol. 350: W.-B. Zhang, Economic Dynamics. X, 232 pages.
Vol. 324: I.M. Bomze, B.M. Potscher. Game Theoretical Foun- 1990.
dations of Evolutionary Stability. VI, 145 pages. 1989. Vol. 351: A. Lewandowski, V. Volkovich (Eds.), Multiobjective
Vol. 325: P. Ferri, E. Greenberg, The Labor Market and Business Problems of Mathematical Programming. Proceedings, 1988.
Cycle Theories. X, 183 pages. 1989. VII, 315 pages. 1991.
Vol. 326: Ch. Sauer, Alternative Theories of Output, Unem- Vol. 352: O. van Hilten, Optimal Firm Behaviour in the Context
ployment, and Inflation in Germany: 1960-1985. XIII. 206 of Technological Progress and a Business Cycle. XII, 229 pages.
pages. 1989. 1991.
Vol. 327: M. Tawada, Production Structure and International Vol. 353: G. Riccil (Ed.), Declslon Processes In Economics.
Trade. V. 132 pages. 1989. Proceedings. 1989. III. 209 pages 1991.
Vol. 328: W. GUth. B. Kalkofen, Unique Solutions for Strategic Vol. 354: M. Ivaldi, A Structural Analysis of Expectation For-
Games. VII, 200 pages. 1989. mation. XII, 230 pages. 1991.
Vol. 329: G. Tillmann, Equity, Incentives. and Taxation. VI, Vol. 355: M. Salomon. Deterministic Lotsizlng Models for Pro-
132 pages. 1989. duction Planning. VII, 158 pages. 1991.
Vol. 330: P.M. Kort, Optimal Dynamic Investment Policies of Vol. 356: P. Korhonen, A. Lewandowski, J . Wallenius (Eds.),
a Value Maximizing Firm. VII, 185 pages. 1989. Multiple Crltena Decision Supporl. Proceedings, 1989. XII, 393
Vol. 331: A. Lewandowski, A.P. Wierzbicki (Eds.), Aspiration pages. 1991.
Based Decision Support Systems. X, 400 pages. 1989. Vol. 358: P. Knottnerus, Linear Models with Correlaled Distur-
Vol. 332: T.R. Gulledge, Jr., L.A. Litteral (Eds.), Cost Analysis bances. VIII, 196 pages. 1991.
Applications of Economics and Operations Research. Vol. 359: E. de Jong, Exchange Rate Determination and Opti-
Proceedings. VII, 422 pages. 1989. mal Economic Policy Under Various Exchange Rate Regimes.
VII, 270 pages. 1991.
Vol. 360: P. Stalder, Regime Translations, Spillovers and Buffer
Stocks. VI, 193 pages. 1991.
Vol. 361: C. F. Daganzo, Logistics Systems Analysis. X, 321
pages. 1991.
Vol. 362: F. Gehreis, Essays In Macroeconomics of an Open
Economy. VII, 183 pages. 1991.
Vol. 363: C. Puppe, Distorted Probabilities and Choice under
Risk. VIII, 100 pages. 1991
Vol. 364: B. Horvath, Are Policy Variables Exogenous? XII,
162 pages. 1991.
Vol. 365: G. A Heuer, U. Leopold-Wildburger. Balanced
Silverman Games on General Discrete Sets. V, 140 pages. 1991.
Vol. 366: J. Gruber (Ed.), Econometric Decision Models. Pro-
ceedings, 1989. VIII, 636 pages. 1991.
Vol. 367: M. Grauer, D. B. Pressmar (Eds.), Parallel Computing
and Mathematical Optimization. Proceedings. V, 208 pages.
1991.
Vol. 368: M. Fedrizzi, J. Kacprzyk, M. Roubens (Eds.),
Interactive Fuzzy Optimization. VII, 216 pages. 1991.
Vol. 369: R. Koblo, The Visible Hand. VIII, 131 pages. 1991.
Vol. 370: M. J. Beckmann, M. N. Gopalan, R. Subramanian
(Eds.), Stochastic Processes and their Applications. Proceedings,
1990. XLI, 292 pages. 1991.
Vol. 371: A. Schmutzler, Flexibility and Adjustment to Infor-
mation in Sequential Decision Problems. VIII, 198 pages. 1991.
Vol. 372: J. Esteban, The Social Viability of Money. X, 202
pages. 1991.
Vol. 373: A. Billot, Economic Theory of Fuzzy Equilibra. XIII,
164 pages. 1992.
Vol. 374: G. Pflug, U. Dieter (Eds.), Simulation and Optimi-
zation. Proceedings, 1990. X, 162 pages. 1992.
Vol. 375: S.-I. Chen, Ch.-L. Hwang, Fuzzy Multiple Attribute
Decision Making. XII, 536 pages. 1992.