TFS-2008-0342
1
Representation of Uncertain Multi-Channel Digital Signal Spaces and Study of Pattern Recognition Based on Metrics and Difference Values on Fuzzy n − cell Number Spaces Guixiang Wang, Peng Shi, Senior Member, IEEE, and Paul Messenger
Abstract—In this paper, we discuss the problem of characterization for uncertain multi-channel digital signal spaces, propose using fuzzy n − cell number space to represent uncertain n − channel digital signal space, and put forward a method of constructing such fuzzy n − cell numbers. We introduce two new metrics and concepts of certain types of difference values on fuzzy n − cell number space, and study their properties. Further, based on the metrics or difference values appropriately defined we put forward an algorithmic version of pattern recognition in an imprecise or uncertain environment, and we also give practical examples to show the application and rationality of the proposed techniques. Index Terms—Uncertain multi-channel digital signals, Fuzzy n − cell numbers, n− dimensional fuzzy vectors, Metrics, Difference values, Pattern recognition
I
I. INTRODUCTION
T is known that in a precise or certain environment, multi-channel digital signals can be represented by elements of multi-dimensional Euclidean space, i.e., crisp multi-dimensional vectors. If however we wish to study multi channel digital signals in an imprecise or uncertain environment, then the signals themselves are imprecise or have no certain bound, and it becomes unwise to use crisp multidimensional vectors to represent them. In this paper, we recommend using fuzzy n − cell numbers to represent imprecise or uncertain multi-channel digital signals, and put forward a method of constructing such fuzzy n − cell numbers.
Manuscript received October 15, 2008. This work is supported in partially by Natural Science Foundations of China (No. 60772006) and Natural Science Foundations of Zhejiang Province, China (No. Y7080044), and Engineering and Physical Sciences Research Council, UK (No. EP/F0219195). G. Wang is with Institute of Operations Research and Cybernetics, Hangzhou Dianzi University, Hangzhou, 310018, China (
[email protected]). P. Shi is with Department of Computing and Mathematical Sciences, University of Glamorgan, Pontypridd, CF37 1DL, UK (
[email protected]). He is also with ILSCM, School of Science and Engineering, Victoria University, Melbourne, 8001, Australia, and School of Mathematics and Statistics, University of South Australia, Adelaide, 5095, Australia. P. Messenger is with Department of Computing and Mathematical Sciences, University of Glamorgan, Pontypridd, CF37 1DL, UK (
[email protected]).
The concept of general fuzzy numbers was introduced by Chang and Zadeh [2] in 1972 with the consideration of the properties of probability functions. Since then both the numbers and the problems in relation to them (see for example [3, 4, 5, 6, 11, 16, 19, 20, 21]) have been widely studied. With the development of theories and applications of fuzzy numbers, this concept becomes more and more important. In [7] Kaleva ever used a special type of n − dimensional fuzzy number, whose sets of cuts are all hyper-rectangles. In 2002 we carefully studied the special type of n − dimensional fuzzy number, and call it fuzzy n − cell number in [14,15]. It has been demonstrated that fuzzy n − cell number is used much more conveniently than general n − dimensional fuzzy numbers in theoretical investigations and some fields of application in [14, 15, 17]. On the other hand, n − dimensional fuzzy vector is also an important concept, which is the Cartesian product of n 1 − dimensional fuzzy numbers. In 1985, Kaufmann and Gupta [8] already studied fuzzy vectors, soon afterwards, Miyakawa and Nakamura et al. [9,10,12] also studied the problems of theories and applications in relating to fuzzy vectors. In 1997, Butnariu [1] studied Methods of solving optimization problems and linear equations in the space of fuzzy vectors. Recently, we [14] showed that fuzzy n − cell numbers and n − dimensional fuzzy vectors can represent each other, and obtained the representations of the joint membership function and the edge membership functions of a fuzzy n − cell number of each other. In a previous paper [15], we defined a metric DL on the fuzzy n − cell number space, and studied its properties. And in paper [14], we again studied this type of metric in regard to two fuzzy n − cell numbers as the form of n − dimensional fuzzy vectors. Although metric DL can be more conveniently used in applications and theoretical investigations, it has some shortcomings. That is, it has a tendency to be rougher, and can not really characterize the degree of difference of two fuzzy n − cell numbers in some applications (see Example 3.1 in Section 3 of this paper). In this paper, in order to discuss the problem of pattern recognition in an imprecise or uncertain environment based on degree of difference, we define two new metrics and some concepts of difference values on fuzzy n − cell number space, which may better characterize the
TFS-2008-0342
2
degree of difference of two fuzzy n − cell numbers in some applications, and study their properties. It is well known that pattern recognition is an important field of research. In this aspect many research achievements have been obtained (for example, see [13]). In this paper, as applications of the metrics and difference values (defined by us), we also study the problem of pattern recognition in an imprecise or uncertain environment, put forward an algorithmic version of pattern recognition based on the metrics or difference values (defined by us) of fuzzy n − cell numbers, and also give examples to show the application and rationality of the method. The organization of the paper is as follows. In Section 2, we give an example to show how to set up fuzzy n − cell numbers to represent imprecise or uncertain multi-channel digital signals. In Section 3, we define two new metrics, and study their properties. In Section 4, we introduce concepts of difference values of two fuzzy n − cell numbers, and examine their properties. In Section 5, an algorithmic version of pattern recognition is given based on the metrics or the difference values defined by us, and examples are also given to show the application and rationality of the method. Finally, in Section 6, we give a brief conclusion of this paper. II. REPRESENTATIONS OF UNCERTAIN MULTI-CHANNEL DIGITAL SIGNALS
A fuzzy set of the Euclidean space R n is a function u : R n → [0,1] . For fuzzy set u , we denote [u ]r = {x ∈ R n : u ( x) ≥ r} for r ∈ [0,1] , and [u ]0 = {x ∈ R n : u ( x) > 0} (the
closure of {x ∈ R n : u ( x) > 0} ). If u is a normal and fuzzy convex fuzzy set of R n , u ( x) is upper semi-continuous, and [u ]0 is compact, then we call u a n − dimensional fuzzy numbers, and denote the n − dimensional fuzzy number space by E n . If u ∈ E , and for each r ∈ [0,1] , [u ] r is a hyper
rectangle, i.e., there exist u i (r ), u i (r ) ∈ R with u i (r ) ≤ u i (r ) , ( i = 1,2, L, n ) such that [u ] r = ∏i =1 [u i (r ), u i (r )] , then we call n
u is a fuzzy n − cell number, and denote the fuzzy n − cell number space by L( E n ) . A n − dimensional fuzzy vector is an ordered class (u1 , u 2 , L , u n ) , where u i ∈ E (i.e., E 1 ), i = 1,2, L, n . In [14], We have shown that fuzzy n − cell numbers and n − dimensional fuzzy vectors can represent each other, and as the representation is unique, L( E n ) and the n − dimensional fuzzy vector space (i.e., the Cartesian product n 64 47 44 8 E × E × L × E ) may be regarded as identical. When exploring and discussing some quantity, properties or laws of movement of phenomena/objects in the physical world, it is essential for us to establish the description space of them. For instance, when the quantity in question is only the one with a single factor, we can take it as a dot in real number field R , that is, the space of quantities corresponding to single factor
can be described by 1 − dimensional Euclidean space R . Similarly, we can describe the quantities with n factors, using n − dimensional Euclidean space R n . However, in the physical world, many phenomena are imprecise or uncertain (such as, have no certain bound). When the quantity discussed by us possesses some imprecise or uncertain attributes, it is unsuitable that we use still R n to represent the space of the quantities (see Remark 2.1). It is our opinion that using the fuzzy n − cell number space discussed in [14, 15] to describe the quantities with some uncertain factors and discuss these quantities in this n − dimensional fuzzy vector space is a more suitable method to reveal the objective laws of things in physical world (see Remark 2.1). In the following example, we demonstrate how we construct a fuzzy n − cell number to represent a quantity that possesses some uncertain attributes based on statistical data. About the algorithmic version of such fuzzy n − cell numbers, we can see the first or second step of the algorithmic version in Section 5. Example 2.1. It is well known that different kinds of terrain or landcover possess the different reflections of the electromagnetic spectrum. Based on this principle, one can set up a method to recognize the category of landcover, a challenging remote sensing classification problem, using spectral and terrain features for vegetation classification in some zone. In remote sensing classification, the colligation of all species covering a zone of 4500 m2 can be boiled down to an element of remote sensing space. We use “Korean Pine accounts for the main part” to denote forest that mainly contains Korean Pines. Because in different “Korean Pine accounts for the main part” areas, there are many different factors such as the difference of the density of Korean Pines, of the species and quantity of other plants, of the physiognomy and so on, the values of reflections of the electromagnetic spectrum are also different. Therefore “Korean Pine accounts for the main part” should not be a certain crisp value but a fuzzy set without certain bound. So, using a fuzzy number to represent the spectral sensitivity level of the “Korean Pine accounts for the main part” is more suitable than using a crisp number. Suppose that we use 4 wave bands: MSS-4, MSS-5, MSS-6, MSS-7. We take 10 samples, and acquire the following data for some zone of “Korean Pine accounts for the main part”: MSS − 4
MSS − 5
MSS − 6
MSS − 7
Sample 1
15 . 01
13 . 30
40 . 50
19 . 37
Sample 2 Sample 3
15 . 60 15 . 82
12 . 56 12 . 79
38 . 81 37 . 70
16 . 35 18 . 16
Sample 4
14 . 90
11 . 70
35 . 50
14 . 75
Sample 5 Sample 6
16 . 10 13 . 80
13 . 80 11 . 94
42 . 10 32 . 10
20 . 75 15 . 54
Sample 7
15 . 90
10 . 98
30 . 87
14 . 29
Sample 8 Sample 9
16 . 82 15 . 50
13 . 67 12 . 58
37 . 64 36 . 10
18 . 62 18 . 02
Sample 10
15 . 38
12 . 48
34 . 08
17 . 45
We can directly work out the following means μ i ( i = 1,2,3,4 ) and standard deviations σ i ( i = 1,2,3,4 ) from the data:
TFS-2008-0342
3
MSS-4 μ i : μ1 = 15.46
MSS-5 μ 2 = 12.58
MSS-6 μ 3 = 36.54
MSS-7 μ 4 = 17.33
σ i : σ 1 = 1.22
σ 2 = 0.88
σ 3 = 3.55
σ 4 = 2.08
From the means and the standard deviations, with ⎧ xi − (μ i − 2σ i ) if xi ∈ [μ i − 2σ i , μ i ] I (0,+∞) ⎪ 2σ i ⎪ ⎪ (μ + 2σ i ) − xi if xi ∈ (μ i , μ i + 2σ i ] I (0,+∞) u i ( xi ) = ⎨ i 2σ i ⎪ ⎪ ⎪0 ⎩
if
xi ∉ [μ i − 2σ i , μ i + 2σ i ] I (0,+∞)
( i = 1,2,3,4 ) we can define 4 triangular model one-dimensional fuzzy numbers u1 , u 2 , u3 and u 4 that respectively correspond with MSS-4, MSS-5, MSS-6 and MSS-7: ⎧ x 1 − 13 . 02 ⎪ ⎪ 17 2. 9. 44 − x1 ⎪ u 1 ( x1 ) = ⎨ 2 . 44 ⎪ ⎪ 0 ⎪⎩ ⎧ x 2 − 10 . 82 ⎪ ⎪ 14 .134. 76− x ⎪ 2 u2 (x2 ) = ⎨ 1 . 76 ⎪ ⎪ 0 ⎩⎪
⎧ x3 ⎪ ⎪ 43 ⎪ u3 ( x3 ) = ⎨ ⎪ ⎪ 0 ⎪⎩
− 29 . 44 7 .1 . 64 − x 3 7 .1
⎧ x 4 − 13 . 17 ⎪ ⎪ 21 .449. 16− x ⎪ 4 u4 (x4 ) = ⎨ 4 . 16 ⎪ ⎪ 0 ⎪⎩
if
x 1 ∈ [13 . 02 , 15 . 46 ]
if
x 1 ∈ (15 . 46 , 17 . 9 ]
if
x 1 ∉ [13 . 02 , 17 . 9 ]
if
x 2 ∈ [10 . 82 , 1 2.58 ]
if
x 2 ∈ (12 . 58 , 1 4.34 ]
if
x 2 ∉ [10 . 82 , 1 4.34 ]
if
x 3 ∈ [ 29 . 44 , 3 6.54 ]
if
x 3 ∈ ( 36 . 54 , 4 3.64 ]
if
x 3 ∉ [ 29 . 44 , 4 3.64 ]
if
x 4 ∈ [13 . 17 , 1 7.33 ]
if
x 4 ∈ (17 . 33 , 2 1.49 ]
if
x 4 ∉ [13 . 17 , 2 1.49 ]
By Theorem 3.1 and 3.2 in [14], we know that u1 , u 2 , u3 and u 4 determine a fuzzy 4 − cell number u = (u1 , u 2, u3 , u 4 ) , and
the membership function of u is u ( x1 , x2 , x3 , x4 ) = min{u1 ( x1 ), u 2 ( x2 ), u 3 ( x3 ), u 4 ( x4 )} ( x1 , x2 , x3 , x4 ) ∈ R 4
Then u can be used to represent “Korean Pine accounts for the main part”. Likewise, from the means and the standard deviations, according to ⎧ (x − μ )2 if xi ∈ (0,+∞) ⎪vi ( xi ) = exp(− i 2i ) i = 1,2,3,4 ⎨ 2σ i ⎪0 if xi ∉ (0,+∞) ⎩ we can also define 4 Gaussian model one-dimensional fuzzy numbers v1 , v2 , v3 and v 4 that respectively correspond with MSS-4, MSS-5, MSS-6 and MSS-7: ⎧ ( x − 15.46) 2 ⎪v1 ( x1 ) = exp(− 1 ) if x1 ∈ (0,+∞) ⎨ 2.98 ⎪⎩ 0 if x1 ∉ (0,+∞)
⎧ ( x − 12.58) 2 ⎪v 2 ( x2 ) = exp(− 2 ) if ⎨ 1.56 ⎪⎩ 0 if ⎧ ( x − 36.54) 2 ⎪v3 ( x3 ) = exp(− 3 ) if ⎨ 25.21 ⎪⎩ 0 if
x2 ∈ (0,+∞) x2 ∉ (0,+∞) x3 ∈ (0,+∞) x3 ∉ (0,+∞)
⎧ ( x − 17.33) 2 ⎪v 4 ( x4 ) = exp(− 4 ) if x 4 ∈ (0,+∞) ⎨ 8.65 ⎪⎩ 0 if x 4 ∉ (0,+∞) and obtain the membership function of the fuzzy 4 − cell number v = (v1 , v2, v3 , v 4 ) determined by v1 , v2 , v3 and v 4 as v ( x 1 , x 2 , x 3 , x 4 ) = min{ v 1 ( x 1 ), v 2 ( x 2 ), v 3 ( x 3 ), v 4 ( x 4 )} , ( x1 , x2 , x3 , x4 ) ∈ R 4 . Then the fuzzy 4 − cell number v can also
be used to represent the “Korean Pine accounts for the main part”. Remark 2.1. Of course, if the quantity to describe is precise and certain, we should use a crisp multi-dimensional vector to represent it. However, if the quantity to describe is imprecise and uncertain, such as “Korean Pine accounts for the main part”, then using a fuzzy n − cell number to represent it is better than using a crisp n − dimensional vector. If we narrowly use a crisp multi-dimensional vector, such as (15.46, 12.58, 36.54 , 17.33) (i.e., the mean vector), to represent “Korean Pine accounts for the main part”, then it can not clearly tell us the relationship of “Korean Pine accounts for the main part” and the zone whose value of reflection of electromagnetic spectrum is (15 . 16 , 12 .80 , 37 .50 , 16 .79 ) since (15.16, 12.80, 37.50, 16.79) ≠ (15.46, 12.58, 36.54, 17.33) . If we use fuzzy n − cell number v = (v1 , v2, v3 , v4 ) to represent it, then we can almost affirm that the zone whose value is (15.16, 12.80, 37.50, 16.79) is part of “Korean Pine accounts for the main part” since v(15.16, 12.80, 37.50, 16.79) = min(0.94, 0.94, 0.93, 0.93) = 0.93 , i.e., the degree of the zone which is “Korean Pine accounts for the main part” is 0.93 . III. METRICS ON FUZZY n − CELL NUMBER SPACE In [3], the authors studied the metric d p ( ⋅ , ⋅ ) (note that in this paper we rewrite d p ( ⋅ , ⋅ ) as D p ( ⋅ , ⋅ ) ) on general n − dimensional fuzzy number space E n , which is defined by 1
D p (u , v) = ( ∫ [d ([u ]r , [v]r )] p dr )1 / p for any u, v ∈ E n , and 0
point out that the metric D p is complete. In [15], we studied the metrics D and DL on L( E n ) , but the two metrics seem to be ‘rough’ in certain applications (see Example 3.1). In the following, other metrics are defined on L( E n ) , which better reveal the difference between two different uncertain quantities (see Example 3.1). Their properties are also discussed such that they may be used appropriately.
TFS-2008-0342
4 Dˆ α , p : L( E n ) × L( E n ) → [0,+∞)
We denote LC ( R n ) = { A : there exist ai ≤ bi , i = 1, 2,L , n such that A = ∏i =1[ai , bi ] } , where n
∏
n
[ai , bi ] is the Cartesian
i =1
and
product [a1 , b1 ] × [a2 , b2 ] × L × [an , bn ] . Theorem 3.1. We define mappings dˆ : LC ( R n ) × LC ( R n ) → [0,+∞)
by
and
and
~ Dα , p : L( E n ) × L( E n ) → [0,+∞) 1 Dˆ α , p (u, v) = ( ∫ [ r ⋅ dˆα ([u ]r , [v]r )] p dr )1 / p
α
0
~ dα : LC(R n ) × LC(R n ) → [0,+∞)
1 ~ ~ Dα , p (u, v) = ( ∫ [r ⋅ dα ([u ]r , [v]r )] p dr )1 / p 0
by
i.e.,
n dˆα ( A, B) = ∑i=1αi ⋅ max{|ai − bi |, | ai − bi |}
Dˆ α , p (u, v)
and
= ( ∫ [ r ∑i =1α i ⋅ max{| ui (r ) − vi (r ) |, | ui (r ) − vi (r ) |}] p dr )1 / p 1
| ai − bi | + | ai − bi | ~ n dα ( A, B) = ∑i=1αi 2
and
for any A = ∏i=1[ai , ai ] ∈ LC(Rn ) and B = ∏i=1[bi , bi ] ∈ LC(Rn ) ,
~ Dα , p (u, v)
n
n
where α = (α1,α2 ,L,αn ) satisfies
∑i=1α i = 1 and α i ≥ 0 , n
i = 1,2, L , n . Then for any A = ∏ i = 1 [ai , ai ], B = ∏ i =1 [bi , bi ], ~ n C = ∏i =1 [c i , c i ] in LC ( R n ) and each k ∈ R , dˆα and d α n
n
satisfies ~ ~ (1) dˆα ( A, B) = dˆα ( B, A) and dα ( A, B) = dα ( B, A) ; ~ (2) dˆα ( A, B) ≥ 0 and d α ( A, B ) ≥ 0 ; ~ (3) dˆα ( A, B) = 0 ⇔ A = B ⇔ dα ( A, B) = 0 ; (4) dˆα ( A, B) ≤ dˆα ( A, C ) + dˆα (C , B) and
n
0
~ ~ d α ( A, B ) ≤ d α ( A, C )
~ + dα (C , B) ;
~ (5) dˆ α ( A + C , B + C ) = dˆ α ( A , B ) and d α ( A + C , B + C )
~ = dα ( A, B ) ; ~ ~ (6) dˆ α (kA, kB) =| k | dˆ α ( A, B) and d α (kA, kB) =| k | d α ( A, B) .
Proof. We only show Proofs (4), (5) and (6) (the other proofs are easy). From
1 1 n ( [ r ∑i =1α i ⋅ (| ui ( r ) − vi (r ) | + | ui ( r ) − vi ( r ) |)] p dr )1 / p 2 ∫0 for any (u, v) ∈ L( E n ) × L(E n ) , where p ≥1, and α = (α1 ,α 2 ,L, =
α n ) satisfies
n
= dˆα ( A, C ) + dˆα (C , B) n dˆα ( A + C , B + C ) = ∑i =1 α i ⋅ max{| ai + ci −bi + ci |, | ai + ci − bi + ci | }
= ∑i =1 α i ⋅ max{| ai + ci −bi − ci |, | ai + ci − bi − ci | }
α i = 1 and α i > 0 , i = 1,2, L , n . Then for
Proof. It is obvious that (1) and (2) of the theorem hold. By the definition of Dˆ α , p , it is obvious that u = v ⇒
n
≤ ∑i =1 α i ⋅ (max{| ai − ci |, | ai − ci |} + max{| ci − bi |, | ci − bi |})
n
i =1
~ any u, v, w ∈ L( E n ) and each k ∈ R , Dˆ α , p and Dα , p satisfy ~ ~ (1) Dˆ α , p (u, v) = Dˆ α , p (v, u ) and Dα , p (u, v) = Dα , p (v, u ) ; ~ (2) Dˆ α , p (u, v) ≥ 0 and Dα , p (u, v) ≥ 0 ; ~ (3) Dˆ α , p (u , v) = 0 ⇔ u = v ⇔ Dα , p (u , v) = 0 ; ~ (4) Dˆ α , p ( u , v ) ≤ Dˆ α , p ( u , w ) + Dˆ α , p ( w , v ) and D α , p ( u , v ) ~ ~ ≤ Dα , p (u , w) + Dα , p ( w, v) ; ~ (5) Dˆ α , p (u + w, v + w) = Dˆ α , p (u , v ) and Dα , p (u + w , v + w ) ~ = Dα , p (u, v) ; ~ (6) Dˆ α , p ( ku , kv ) = | k | Dˆ α , p ( u , v ) and D α , p ( ku , kv ) ~ =| k | Dα , p (u , v) .
n dˆα ( A, B) = ∑i =1 α i ⋅ max{| ai − bi |, | ai − bi |}
≤ ∑i =1 α i ⋅ max{| ai − ci | + | ci − bi |, | ai − ci | + | ci − bi |}
∑
Dˆ α , p (u , v) = 0 . Otherwise, let Dˆ α , p (u , v) = 0 . Then we have
∫ [r ∑ 1
n
0
i =1
α i ⋅ max{| u i ( r ) − v i ( r ) |, | u i ( r ) − v i ( r ) |}] p d r = 0 .
Taking note of αi > 0 , we see that r (∑i =1α i ⋅ max{| ui (r ) − n
n
= dˆα ( A, B) n dˆα (kA, kB) = ∑i =1 α i ⋅ max{| kai − kbi |, | kai − kbi |}
= ∑i =1 α i ⋅ max{| k ai −k bi |, | k ai − k bi |} n
=| k | dˆα ( A, B)
~ we see that (4), (5) and (6) of the theorem hold for dˆα . For dα , we can similarly prove that (4), (5) and (6) of the theorem also hold. Theorem 3.2. We define mappings
vi (r ) |, | ui (r ) − vi (r ) |}) = 0 holds for r almost everywhere on
[0,1] . Further, we have that
∑
n
α i ⋅ max{| ui (r ) − vi (r ) |,
i =1
| ui ( r ) − vi ( r ) |} = 0 holds for r almost everywhere on [0,1] , so
we can see that ui (r ) = vi (r ) and ui (r ) = vi (r ) holds for r almost everywhere on [0,1] for i = 1,2, L , n . Therefore, we obtain that [u ]r = [v]r holds for r almost everywhere on [0,1] , so we know that u = v holds by Lemma 2.1 in [18]. Thus, Dˆ (u , v) = 0 ⇔ u = v holds. Likewise, we can prove α,p
TFS-2008-0342
5
~ Dα , p (u , v) = 0 ⇔ u = v , so (3) of the theorem holds. From
Theorem 3.1, we have 1 n Dˆ α , p (u , v ) = ( ∫ [r ∑i =1 α i ⋅ max{| ui (r ) − vi (r ) |, | ui (r ) − vi (r ) |}] p dr )1 / p 0
1
≤ [ ∫ (r ⋅ [dˆα ([u ]r , [ w]r ) + dˆα ([ w]r , [v]r )] ) p dr ]1 / p
~ 1ˆ 1 ~ Dα, p ≤ Dα, p ≤ Dˆα, p ≤ Dp ≤ D , i.e., Dˆ α , p (u , v) ≤ Dα , p (u , v) 2 2 ≤ Dˆ α , p (u , v) ≤ D p (u , v) ≤ D(u , v) for any u, v ∈ L( E n ) ( D is (1)
discussed in [15]). (2) Dˆ α , p ≤
0
1
1 ≤ ( ∫ [r ⋅ dˆα ([u ]r , [ w]r )] p dr )1 / p + ( ∫ [r ⋅ dˆα ([ w]r , [v]r ) ]p dr )1 / p 0
0
= Dˆ α , p (u , w) + Dˆ α , p ( w, v)
~ ~ ~ The proofs of Dα , p (u , v) ≤ Dα , p (u , w) + Dα , p ( w, v) , (5) and (6)
can be similarly proved.
~ Remark 3.1. From Theorems 3.1 and 3.2, we know dˆα , dα ~ and Dˆ , D are metrics on LC( R n ) and L( E n ) , respectively, α,p
α,p
and satisfy translation invariance, absolute homogeneity. Also, from the factor r of the integrands in the definitions of ~ Dˆ α , p (u, v) and Dα , p (u, v) , we can see that the bigger the degrees of the points are, which belong to the fuzzy n − cell numbers u and v , the greater the effects on the metric of u and v . This is true in reality. Example 3.1. Let u, v, w be the 2 − cell numbers defined by: u = (u1 , u 2 ) , v = (v1 , v2 ) and w = ( w1 , w2 ) , where, ⎧x ⎪ ui ( x) = ⎨ 2 − x ⎪0 ⎩
if if if
x ∈ [ 0 ,1] ⎧x − 2 ⎪ x ∈ (1, 2 ] , v i ( x ) = ⎨ 4 − x ⎪0 x ∉ [0,2 ] ⎩
if if
x ∈ [ 2,3] x ∈ (3, 4]
if
x ∉ [ 2, 4 ]
( i = 1,2 ) if ⎧x ⎪ w1 ( x) = ⎨ 2 − x if ⎪0 if ⎩
⎧ x − 2 if ⎪ x ∈ (1,2] , w2 ( x) = ⎨4 − x if ⎪0 x ∉ [0,2] if ⎩ x ∈ [0,1]
x ∈ [2,3] x ∈ (3,4] x ∉ [2,4]
Then we know that ui (r ) = r , ui (r ) = 2 − r , vi (r ) = 2 + r , vi (r ) = 4 − r ( i = 1,2 ), w1 (r ) = r , w1 (r ) = 2 − r , w2 (r ) = 2 + r
and w2 (r ) = 4 − r for r ∈ [0,1] . From the definitions of DL , we have DL (u , v) = 2 = DL (u , w) , i.e., DL can not tell us the difference of DL (u , v) and DL (u , w) , so we say that DL seems to be ‘rough’ (similar proof for D ). However, as a matter of fact, DL (u , v) and DL (u , w) should have some difference. ~ Taking α = (1 / 2, 1 / 2) , from the definitions of Dα , p , we can ~ obtain Dα , p (u, v) =
2 (1 + p )1 / p
1 ( p + 1)1 / p
≤
~ 1 > = Dα , p (u , w) , this (1 + p )1 / p
accord with fact. If we restrain the metric D p (i.e. d p ( ⋅ , ⋅ ) defined by Diamond in [3], see paragraph 1 of this section) on general n − dimensional fuzzy number space E n into on L( E n ) , then it also becomes a metric on L( E n ) . In the following, we give ~ the relationships of the metrics Dˆ α , p , Dα , p and D p . ~ Theorem 3.3. Metrics Dˆ α , p , Dα , p and D p satisfy
1 1 DL ≤ D , i.e., Dˆ α , p ( u , v ) ( p + 1) 1 / p ( p + 1) 1 / p 1 DL (u, v) ≤ D(u, v) for any u, v ∈ L( E n ) . ( p + 1)1 / p
Proof. For any u, v ∈ L( E n ) and r ∈ [0,1] , by the definitions ~ of dˆα and dα , we have 1 ˆ dα ([u ]r , [v]r ) 2 1 n = ∑ i =1α i ⋅ max{| ui (r ) − vi (r ) |, | ui (r ) − vi (r ) |} 2 =
| ui (r ) − vi (r ) | + | ui (r ) − vi (r ) | + | ui (r ) − vi (r ) | − | ui (r ) − vi (r ) | 1 n α ⋅ ∑ i =1 i 2 2
≤ ∑i =1α i ⋅ n
≤ ∑i =1α i ⋅ n
| ui (r ) − vi ( r ) | + | ui ( r ) − vi (r ) | + | ui ( r ) − vi ( r ) | + | ui (r ) − vi ( r ) | 4 | ui (r ) − vi ( r ) | + | ui ( r ) − vi (r ) |
≤ ∑ i =1 α i ⋅ n
2
~ ( = dα ([u ]r , [v]r ))
2 max{| ui ( r ) − vi ( r ) |, | ui ( r ) − vi (r ) |} 2
= dˆα ([u ]r , [v]r )
From this, we can directly obtain
~ 1 ˆ Dα , p (u , v) ≤ Dα , p (u , v) 2
≤ Dˆ α , p (u , v) .
By Theorem 4.4 in [15], we know d L ≤ d ≤ n d L , where, d L ( A, B) = max1≤i≤n {| ai − bi |, | ai − bi |} for any A = ∏i =1[ai , ai ] n
∈ LC ( R n ) and B = ∏i =1 [bi , bi ] ∈ LC ( R n ) . Therefore, for any n
u, v ∈ L( E n ) , we have Dˆ (u , v ) α,p
= ( ∫ [r ∑i =1α i ⋅ max{| ui (r ) − vi (r ) |, | ui (r ) − vi (r ) |}] p dr )1 / p 1
n
0
≤ ( ∫ [r ∑i =1α i ⋅ max1≤ j ≤ n {| u j (r ) − v j (r ) |, | u j (r ) − v j (r ) |}] p dr )1 / p 1
n
0
= ( ∫ [r ⋅ d L ([u ]r ,[v]r )∑i =1α i ] p dr )1 / p 1
n
0
1
≤ ( ∫ [d ([u ]r ,[v]r )] p dr )1 / p (= D p (u , v)) 0
1
≤ ( ∫ [supr∈[ 0,1] d ([u ]r ,[v]r )] p dr )1 / p 0
1
= ( ∫ [ D(u , v)] p dr )1 / p 0
≤ D(u , v)
Thus, we also obtain Dˆ α , p (u , v) ≤ D p (u, v ) ≤ D(u , v) and the proof of (1) of the theorem is complete. From Dˆ α , p (u , v ) 1
= ( ∫ [ r ∑ i =1α i ⋅ max{| ui ( r ) − vi ( r ) |, | ui ( r ) − vi ( r ) |}] p dr )1 / p 0
n
TFS-2008-0342
6 and call Lα ,a (u, v) and Rα ,a (u, v ) a left difference value and a
1
≤ ( ∫ [r ⋅ d L ([u ]r ,[v]r )] p dr )1 / p 0
right difference value of u and v (with respect to the weight α and parameter a ), respectively. And we denote
1
≤ ( ∫ [r ⋅ sup r∈[ 0,1] d L ([u ]r ,[v]r )] p dr )1 / p 0
1
= ( ∫ [ rDL (u , v)] p dr )1 / p
Δ α ,a (u, v) =
0
1
≤ D (u , v)( ∫ r p dr )1 / p 0
1 [ Lα ,a (u, v) + Rα ,a (u, v)] 2
i.e.,
1 = D (u , v ) ( p + 1)1 / p
and Theorem 4.5 in [15], we can see that (2) of the theorem holds. Remark 3.2. From (1) of Theorem 3.3, we see ~ ~ 1 ˆ Dα , p ≤ Dα , p ≤ Dˆ α , p , i.e., Dˆ α , p and Dα , p are equivalent, so we 2 ~ know that Dˆ α , p and Dα , p induce equivalent topologies on
∑
n
i =1
Δ α ,a (u, v) =
1
α i ∫ 0 r | (ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr
⎛ n α 1 r[u (r ) + v (r ) + u (r ) + v (r )]dr ⎞ ⎜ ∑i =1 i ∫ 0 i ⎟ i i i ⎝ ⎠
a
and call Δ α ,a (u, v) a difference value of u and v (with respect to the weight α
α = (α 1 , α 2 , L , α n ) ∈ R
and parameter a ), where,
with
n
∑
n i =1
α i = 1 and α i ≥ 0
L( E ) by the knowledge of topological space.
( i = 1, 2 , L , n ) , and a ∈ (0,+∞) . Specially, we denote
IV. DIFFERENCE VALUES ON FUZZY n − CELL NUMBER SPACE
1 1 1 Δ a (u , v ) = Δ α ,a (u , v ) as α = ( , , L , ) , and Δ a (u , v) n n n = Δ 1, a (u, v) as u, v ∈ E .
n
In Section 3, we discussed metrics on L( E n ) . But sometimes these have some shortcomings demonstrating the difference of two objects. For example, we consider that the degree of difference of 1 and 2 is bigger than the degree of difference of 1010 and 1010 + 1 though their metrics (Euclidean metric) measure both are 1 . A mapping from the Cartesian product X × X of a set X into R needs to satisfy stronger conditions in order that it can become a metric, and this brings limitations in some applications. The measure used to characterize the differences does not need to satisfy all metric conditions, for example, when we set up a method of pattern recognition basing on the principle of minimal difference (i.e. the principle of maximal likelihood), the measure used to characterize the differences does not need to satisfy all metric conditions. To conveniently set up methods of pattern recognition using fuzzy n − cell numbers, we introduce the concepts of difference values on L( E n ) , and study the properties. Let u ∈ L ( E ) and α = (α 1 , α 2 , L , α n ) ∈ R n
∑
n
i =1
αi = 1
αi ≥ 0
and
( i = 1,2,L, n
).
We
n
satisfy denote
M α (u ) = ∑i =1α i ∫ r[u i ( r ) + u i ( r )]dr , denote M (u ) = M α (u ) n
1
0
1 1 1 as α = ( , , L, ) , and denote M (u ) = M 1 (u ) as u ∈ E . n n n Definition 4.1. Let u, v ∈ L( E n ) with M α (u ), M α (v) ≥ 0 and M α (u ) + M α (v) ≠ 0 . We denote Lα ,a (u, v) =
∑
n
i =1
1
α i ∫ 0 2r | ui (r ) − vi (r ) | dr
Remark 4.1. (1) Generally speaking, we consider that the degree of the difference of two numbers is related not only by the metric of them but also by the sizes of them. As the metrics are the same, the bigger the sizes of the two numbers are, the smaller the degree of their difference is. The denominator [Mα (u) + Mα (v)]a in the definition of Δ α ,a just plays the action (see Example 4.1), and the exponent a in [ M α (u ) + M α (v)]a can be properly chosen accordingly to the case in question. (2) Taking the note of that 1 1 n n
Rα ,a (u, v) =
n
=
∑
n i =1
[ M α (u ) + M α (v)]
a
∑
n i=1
a
⎛ M (ui ) + M (vi ) ⎞ ⎟ =1 ⎟ ⎝ Mα (u) + Mα (v) ⎠
αi ⎜⎜
1
α i ∫ r[| u i ( r ) − vi ( r ) | + | u i ( r ) − vi (r ) |]dr 0
1 ⎛ n ⎞ ⎜ ∑ i = 1 α i ∫ 0 r[u i ( r ) + vi (r ) + u i ( r ) + vi ( r )]dr ⎟ ⎝ ⎠
a
1 ⎛ ⎞ r[u i ( r ) + vi (r ) + u i ( r ) + vi ( r )]dr ⎜ ⎟ ∫ n 0 ⎟ = ∑ αi ⎜ 1 i =1 ⎜⎜ n α ⎟⎟ r [ u ( r ) v ( r ) u ( r ) v ( r )] dr + + + i i i ⎝ ∑ i =1 i ∫ 0 i ⎠
⋅
1
α i ∫ 0 2r | ui (r ) − vi (r ) | dr
1 n
1 1 1 does not necessarily hold as α ≠ ( , , L , ) or a ≠ 1 , from n n n Δ α , a (u , v )
and i =1
i =1
⎛ M (ui ) + M (vi ) ⎞ ⎟ = 1 holds as ⎟ ⎝ M α (u ) + M α (v) ⎠
α i ⎜⎜
α = ( , ,L, ) and a = 1 , and
[ M α (u ) + M α (v)]a
∑
∑
a
n
∫
1 0
r[| u i ( r ) − vi ( r ) | + | u i ( r ) − vi ( r ) |]dr
⎛ 1 ⎞ ⎜ ∫ 0 r[u i ( r ) + vi ( r ) + u i ( r ) + vi ( r )]dr ⎟ ⎝ ⎠ a
=
∑
n i =1
⎛ M (u i ) + M (vi ) ⎞ ⎟ Δ a (u i , vi ) ⎟ ⎝ M α (u ) + M α (v) ⎠
α i ⎜⎜
a
a
TFS-2008-0342
7
we can directly see that Δ1 (u , v) is a convex combination of
~ ~ so we have Dα , p (u, v) = Dα , p (u′, v′) and Δα ,a (u, v) > Δα ,a (u′, v′) .
Δ1 (ui , vi ) , i = 1,2,L, n , but Δ α ,a (u, v) is not necessarily a
convex
combination
of
Δ a (ui , vi )
,
i = 1,2,L, n
as
1 1 1 α ≠ ( , ,L, ) or a ≠ 1 . n n n Example 4.1. Let u, v, u ′, v′ be the 2 − cell numbers defined by: u = (u1 , u 2 ) , v = (v1 , v 2 ) , u ′ = (u1′ , u ′2 ) and
Property 4.1. Let u, v ∈ L( E n ) with M α (u ), M α (v) ≥ 0 and M α (u ) + M α (v) ≠ 0 , α = (α1,α2 ,L,αn ) ∈ Rn with
if if if
xi ∈ [ 0,1] x i ∈ (1, 2 ] x i ∉ [ 0, 2 ]
⎧ − 10 + x i ⎪ v i ( x i ) = ⎨ 12 − x i ⎪ 0 ⎩
if if
x i ∈ [10 ,11] x i ∈ (11,12 ]
if
x i ∉ [10 ,12 ]
if
x i ∈ [100 ,101]
if if
x i ∈ (101,102 ] x i ∉ [100 ,102 ]
⎧ − 100 + x i ⎪ u i′ ( x i ) = ⎨ 102 − x i ⎪ 0 ⎩ ⎧ − 110 + x i ⎪ v i′ ( x i ) = ⎨ 112 − x i ⎪ 0 ⎩
if
x i ∈ [10 ,11]
if if
x i ∈ (11,12 ] x i ∉ [10 ,12 ]
(2) Δ α ,a (u, v) = 0 if and only if u = v ; (3) Δ α ,a (u, v) = Δ α ,a (v, u ) ; (4) Δα ,a (u, w) ≥ Δα ,a (v, w) for any w ∈ L( E n ) and u ≤ v ≤ w ; (5) Δ α ,a (u + w, v + w) ≤ Δ α ,a (u , v) for any w ∈ L( E n ) and w ≥ (0ˆ,0ˆ, L ,0ˆ) ;
(6) Δ α ,a (ku, kv) = k 1−a Δ α ,a (u , v) for any k > 0 . Proof. It is obvious that the conclusions (1) and (3) hold. The proof of conclusion (2) can also be completed by imitating the proof of (3) of Theorem 3.2 by using Lemma 2.1 in [18]. From u ≤ v ≤ w , we know that ui (r ) ≤ vi (r ) ≤ wi (r ) and
u i ( r ) ≤ v i ( r ) ≤ w i ( r ) ( i = 1,2,L, n ), so | ui (r ) − wi (r ) | ≥| v i (r ) − wi (r ) | and | u i (r ) − wi (r ) |≥| vi (r ) − wi (r ) | . Therefore,
we have
∑
vi (r ) = 10 + r , vi (r ) = 12 − r , ui′ (r ) = 100 + r , ui′ (r ) = 102 − r , vi′ (r ) = 110 + r , vi′ (r ) = 112 − r ( r ∈ [0,1] , i = 1,2 ). From the ~ definitions of Dα , p and Δ α ,a , we can obtain ~ Dα , p ( u , v )
1 ⎛ 20 p ⎞ ⎟ = ⎜⎜ 2 ⎝ (1 + p ) ⎟⎠
On the other hand, from ui (r ) ≤ vi (r ) ≤ wi (r ) 1
1/ p
1/ p
∑
≥
∑
n
1
α i ∫ r[| ui ( r ) − wi (r ) | + | ui (r ) − wi (r ) |]dr
n
i =1
0
(M α (u ) + M α ( w) )a
1
α i ∫ r[| vi (r ) − wi (r ) | + | vi (r ) − wi (r ) |]dr
i =1
0
(M α (v) + M α (w) )a
= Δα , a (v, w)
]
1/ p
so conclusion (4) holds. From w ≥ 0ˆ , we know M α ( w) ≥ 0 , so we have Δ α , a (u + w, v + w)
Δα , a (u , v) 1
1
0
0
α1 ∫ r (| r − 10 − r | + | 2 − r − 12 + r |)dr + α 2 ∫ r (| r − 10 − r | + | 2 − r − 12 + r |)dr
⎛ ⎞ ⎜ α1 ∫ 0 r (r + 10 + r + 2 − r + 12 − r ) dr + α 2 ∫ 0 r (r + 10 + r + 2 − r + 12 − r ) dr ⎟ ⎝ ⎠ 10 = a 12 1
1
=
a
=
Δα , a (u′, v′) 1
1
0
Δα , a (u , w) =
1 ⎛ 20 p ⎞ ⎜ ⎟ 2 ⎜⎝ (1 + p ) ⎟⎠
1
α1 ∫ r (| 100 + r − 110 − r | + | 102 − r − 112 + r |)dr + α 2 ∫ r (| 100 + r − 110 − r | + | 102 − r − 112 + r |)dr 0
n
so we can obtain 0 ≤ M α (u) ≤ M α (v) ≤ M α (w) . Thus, we have
]
1/ p
=
and
ui (r) ≤ vi (r) ≤ wi (r) ( i = 1,2,L, n ), we can also see that M α (u ) 0
+ α 2 (| 100 + r − 110 − r | + | 102 − r − 112 + r |))] p dr
=
1
0
= ∑i =1α i ∫ r[ui ( r ) + ui ( r )]dr ≤ ∑i =1α i ∫ r[vi ( r ) + vi ( r )]dr = M α (v) ,
1⎡ 1 [ r (α1 (| 100 + r − 110 − r | + | 102 − r − 112 + r |) 2 ⎢⎣ ∫ 0
=
0
n
~ Dα , p (u ′, v′) =
1
α i ∫ r[| ui (r ) − wi ( r ) | + | ui (r ) − wi (r ) |]dr
n
i =1
≥ ∑i =1α i ∫ r[| vi ( r ) − wi ( r ) | + | vi (r ) − wi ( r ) |]dr
n
+ α 2 (| r − 10 − r | + | 2 − r − 12 + r |))] p dr
αi = 1
(1) Δ α ,a (u, v) ≥ 0 ;
and i = 1,2 . Then we know that ui (r ) = r , ui (r ) = 2 − r ,
1 1 = ⎡ ∫ [ r (α 1 (| r − 10 − r | + | 2 − r − 12 + r |) 2 ⎢⎣ 0
n
i =1
and α i > 0 ( i = 1,2,L, n ), and a ∈ (0,+∞) . Then
v′ = (v1′ , v′2 ) , where, ⎧ xi ⎪ u i ( xi ) = ⎨ 2 − xi ⎪0 ⎩
∑
0
1 1 ⎛ ⎞ ⎜ α1 ∫ 0 r (100 + r + 110 + r + 102 − r + 112 − r )dr + α 2 ∫ 0 r (100 + r + 110 + r + 102 − r + 112 − r )dr ⎟ ⎝ ⎠ 10 = 212 a
a
=
∑
n i =1
1
α i ∫ r[| (u + w)i (r ) − (v + w)i (r ) | + | (u + w)i (r ) − (v + w)i (r ) |]dr 0
1 ⎛ n ⎞ ⎜ ∑ i =1 α i ∫ 0 r[(u + w) i (r ) + (v + w) i (r ) + (u + w) i (r ) + (v + w) i (r )]dr ⎟ ⎝ ⎠
∑
a
1
n i =1
α i ∫ r[| ui (r ) + wi (r ) − vi (r ) − wi (r ) | + | ui (r ) + wi (r ) − vi (r ) − wi (r ) |]dr 0
1 ⎛ n ⎞ ⎜ ∑ i =1 α i ∫ 0 r[ui (r ) + wi (r ) + vi (r ) + wi (r ) + ui (r ) + wi (r ) + vi (r ) + wi (r )]dr ⎟ ⎝ ⎠
∑
n i =1
1
α i ∫ r[| ui (r ) − vi (r ) | + | ui ( r ) − vi ( r ) |]dr 0
(M α (u) + M α (v) + 2M α ( w) )a
a
TFS-2008-0342
≤
∑
n i =1
8 Δα ,2 (u, w)
1
α i ∫ r[| ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr 0
=
= Δα , a (u , v)
i.e., conclusion (5) holds. For any k > 0 , we have Δα , a (ku , kv) =
∑
=
1
α r[| (ku )i (r ) − (kv)i (r ) | + | (ku )i (r ) − (kv)i (r ) |]dr i =1 i ∫ n
0
1 ⎛ ⎞ ⎜ ∑i =1α i ∫ 0 r[(ku )i (r ) + (kv)i (r ) + (ku )i (r ) + (kv)i (r )]dr ⎟ ⎝ ⎠
=
a
=
k ∑i =1α i ∫ r[| ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr 0
1 ⎛ ⎞ k a ⎜ ∑i =1α i ∫ r[ui (r ) + vi (r ) + ui (r ) + vi (r )]dr ⎟ 0 ⎝ ⎠ = k 1− a Δα , a (u , v)
a
Example 4.2). Comparing Δ α ,a (u , w) =
n
i =1
1
α i ∫ 0 r[| ui (r ) − wi (r ) | + | ui (r ) − wi (r ) |]dr
(M
(u ) + M α ( w) )
a
α
1 1 ⎛ ⎞ ⎜α1 ∫ 0 r[0 + 2 + r + 0 + 4 − r]dr + α2 ∫ 0 r[0 + 2 + r + 0 + 4 − r]dr ⎟ ⎝ ⎠
1 3
1
α1 ∫ r(| 0 − 1 − r | + | 0 − 3 + r |)dr + α2 ∫ r(| 0 − 1 − r | + | 0 − 3 + r |)dr 0
0
1 1 ⎞ ⎛ ⎜α1 ∫ 0 r(0 + 1 + r + 0 + 3 − r)dr + α2 ∫ 0 r(0 + 1 + r + 0 + 3 − r)dr ⎟ ⎠ ⎝
1
2
⎧19 ⎪10 + xi ⎪⎪ 1 u i ( xi ) = ⎨ − xi ⎪10 ⎪ 0 ⎩⎪
if
19 9 ,− ] 10 10 9 1 xi ∈ ( − , ] 10 10 19 1 xi ∉ [ − , ] 10 10
xi ∈ [ −
if if
Δ α ,a (u, v) =
1
n
i =1
α i ∫ 0 r[| ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr
(M
α (u ) + M α (v ) )
a
we can see that the untruth of Δ α ,a (u , w) ≥ Δ α ,a (u, v) is caused only by (M α (u ) + M α ( w) )
a
(M α (u ) + M α ( w))a
and
≥ (M α (u ) + M α (v) ) , and when a
(M α (u ) + M α (v))a
are properly
smaller, u ≤ v ≤ w can imply Δ α ,a (u , w) ≥ Δ α ,a (u, v) . So, in general, we may choose a in (0,1] such that Δ α ,a can reasonably characterize the degree of the difference of two fuzzy n − cell numbers. (2) Generally speaking, the difference value Δ α ,a does not satisfy the property of the triangular inequality, i.e., the inequality Δ α , a ( u , v ) ≤ Δ α , a ( u , w ) + Δ α , a ( w , v ) does not necessarily hold for u, v, w ∈ L( E n ) . Example 4.3 can show it. Example 4.2. Let u, v, w be the 2 − cell numbers defined by: u = (0ˆ, 0ˆ) , v = (v , v ) and w = ( w , w ) , where 1
⎧−1 + xi ⎪ vi ( xi ) = ⎨ 3 − xi ⎪0 ⎩
if if if
2
1
2
xi ∈ [1,2] ⎧−2 + xi ⎪ xi ∈ (2,3] , wi ( xi ) = ⎨ 4 − xi ⎪ 0 xi ∉ [1,3] ⎩
if if if
xi ∈ [2,3] xi ∈ (3,4] xi ∉ [ 2,4]
and i = 1, 2 . Then we know that u i ( r ) = 0 , u i ( r ) = 0 , vi (r ) = 1 + r , vi (r ) = 3 − r , wi (r ) = 2 + r , and wi (r ) = 4 − r
( r ∈ [0,1] , i = 1,2 ), so we have u ≤ v ≤ w , but we can see Δ α ,a (u , w) ≤ Δ α ,a (u, v ) , from
2
1 2 Example 4.3. Let u, v, w be the 2 − cell numbers defined by: u = (u , u ) , v = (1ˆ,1ˆ) and w = (2ˆ,2ˆ) , where
and i = 1,2 . Then we know that ui (r) = −
to
∑
2
=
so conclusion (6) holds. Therefore, the proof of the theorem is completed. Remark 4.2. (1) Although the conclusion (4) of Theorem 4.1 holds, u ≤ v ≤ w does not imply Δ α ,a (u , w) ≥ Δ α ,a (u, v) (see
∑
0
0
1
1
n
α1 ∫ r[| 0 − 2 − r | + | 0 − 4 + r |]dr + α2 ∫ r[| 0 − 2 − r | + | 0 − 4 + r |]dr
and Δα ,2 (u, v)
n
n
1
1
(M α (u ) + M α (v) )a
19 1 + r , ui (r) = − r , 10 10
vi (r ) = 1 , vi (r ) = 1 , wi (r ) = 2 , wi (r ) = 2 ( r ∈ [0,1] , i = 1,2 ), 29 1 and Δ α ,1 (v, w) = , 11 3 it implies Δ α ,1 (u , v) > Δ α ,1 (u , w) + Δ α ,1 (v, w) .
so we have Δ α ,1 (u, v) = 19 , Δ α ,1 (u , w) =
Example 4.4. Let u, v, w be the 2 − cell numbers (see Fig.1) defined by: u = (u1 , u 2 ) , v = (2ˆ, 2ˆ) and w = ( w1 , w2 ) , where
⎧ xi ⎪ u i ( x i ) = ⎨2 − x i ⎪0 ⎩ ⎧− 1 + x i ⎪ wi ( xi ) = ⎨ 3 − xi ⎪ 0 ⎩
if if if if if if
xi ∈ [0,1] xi ∈ (1,2] xi ∉ [0,2] xi ∈ [1,2] xi ∈ (2,3] xi ∉ [1,3]
and i = 1,2 . Then we know that ui (r ) = r , ui (r ) = 2 − r , vi (r ) = 2 , vi (r ) = 2 , wi (r ) = 1 + r , wi (r ) = 3 − r ( r ∈ [0,1] , ~ 1 ~ i = 1,2 ), so we can obtain that Dα , p (u , v) = = Dα , p (u, w) and 2 2 Δ α ,a (u, v) = a = Δ α ,a (u , w) . However, it is obvious (see Fig. 3 4.1) that the degree of the difference of u and v is different from the degree of the difference of u and w .
TFS-2008-0342
9
u ( x1 , x2 )
∈ R n with
x2
2
2
x1
v( x1 , x2 )
α i = 1 and α i ≥ 0 ( i =1,2,L, n ), and a ∈ (0,+∞) .
Example 4.5. Let u, v, w be the 2 − cell numbers defined in Example 4.4. Then Λα ,1 (u , v)
1 1
n
i =1
1 1 1 We denote Λ a (u , v) = Λ α ,a (u, v) as α = ( , ,L, ) , and n n n Λ a (u, v) = Λ1,a (u , v) as u, v ∈ E .
1
0
∑
1 ⎛ 2 ⎞ = ⎜ ∑i =1α i ∫ r (| r − 2 | + | 2 − r − 2 |) dr ⎟ 0 ⎝ ⎠
x2
1 2 2 2 2 ⋅ exp⎛⎜ ∑i =1α i (| ∫ tdt − ∫ 1dt |+ | ∫ (2 − t )dt − ∫ 1dt |) ⎞⎟ 0 2 1 2 ⎝ ⎠
3
1
= ∫ 2rdr ⋅ exp(1) 0
1
=e
2 1
0
and
1
2
Λ α ,1 (u, w)
3
x1
w( x1 , x2 )
x2
=
3 1
1
2
3
x1
2 − cell numbers u, v
and
w
in Example 4.4
the difference of the degree of the difference u and v with the ~ degree of the difference of u and w since Dα , p (u, v) ~ = Dα , p (u , w) and Δ α ,a (u, v) = Δ α ,a (u , w) , but we see that the two degrees of the differences indeed have some differences. In order to overcome the defect, we introduce the following concept. Definition 4.2. Let u, v ∈ L( E n ) . We denote Λα , a (u , v ) ⎛ n ⎞ = ⎜ ∑i =1α i ∫ r[| (ui (r ) − vi (r ) | + | ui ( r ) − vi ( r ) |]dr ⎟ 0 ⎝ ⎠ u i (0) v i (0) u i * (1) v i * (1) n ⎛ ⎞ ⋅ exp⎜ a ∑i =1 α i (| ∫ ui (t )dt − ∫ vi (t )dt |+ | ∫ ui (t )dt − ∫ vi (t )dt |) ⎟ u i * (1) v i * (1) u i (0) vi ( 0) ⎝ ⎠ 1
and call Λ α ,a (u , v) a difference value of u and v (with respect to the weight α u (1) =
u i (1) + u i (1)
2
2rdr ⋅ exp(0)
difference value Λ α ,a is more suitable than metrics and
In fact, sometimes, the degree of the difference of two fuzzy numbers is not only related with the metric and the sizes of them, but also related with the degree of fuzzy (we call it fuzzy degree) of them. Example 4.4 shows that for the u, v, w given, ~ the metric Dα , p and the difference value Δ α ,a can not tell us
∗ i
1 0
=1
2
Fig. 1. Fuzzy
∫
so we see Λ α ,1 (u, v) > Λ α ,1 (u , w) . Therefore, in this case, the
1 0
1 ⎛ 2 ⎞ = ⎜ ∑ i = 1 α i ∫ r (| r − 1 − r | + | 2 − r − 3 + r |)dr ⎟ 0 ⎝ ⎠ 1 2 2 3 2 ⎛ ⋅ exp⎜ ∑ i = 1 α i (| ∫ tdt − ∫ (t − 1)dt |+ | ∫ (2 − t )dt − ∫ (3 − t )dt |) ⎞⎟ 0 1 1 2 ⎝ ⎠
∗ i
, v (1) =
and parameter a ), where,
vi (1) + vi (1)
2
, and α = (α1 ,α 2 ,L,α n )
difference value Δ α ,a to characterize the degree of the difference of two fuzzy n − cell numbers. Property 4.2. Let u, v ∈ L( E n ) , α = (α1 , α 2 , L , α n ) ∈ R n with
∑
n
i =1
α i = 1 and α i > 0 ( i = 1,2,L, n ), and a ∈ (0,+∞) .
Then (1) Λ α ,a (u , v) ≥ 0 ; (2) Λ α ,a (u , v) = 0 if and only if u = v ; (3) Λ α ,a (u , v) = Λ α ,a (v, u ) ; (4) Λ α ,a (u + bˆ, v + bˆ) = Λ α ,a (u, v ) for any b ∈ R ; (5) Λ α ,a (ku, kv) =| k | Λ α ,|k |a (u, v) for any k ∈ R . Proof. It is obvious that the conclusions (1) and (3) hold. The proof of conclusion (2) can also be completed by imitating the proof of (3) of Theorem 3.2 by using Lemma 2.1 in [18]. For any b ∈ R and i = 1,2,L, n , we have
(ui + bˆ)(t ) = sup{r ∈ [0,1] : t ∈ [(ui + bˆ)(r ), (ui + bˆ)(r )]} = sup{r ∈ [0,1] : t − b ∈ [ui (r ), ui (r )]} = ui (t − b) hence Λα , a (u + bˆ, v + bˆ) 1 ⎛ n ⎞ = ⎜ ∑i =1α i ∫ r | (ui (r ) + b − vi ( r ) − b | + | ui ( r ) + b − vi ( r ) − b |]dr ⎟ 0 ⎝ ⎠
TFS-2008-0342
10
u i *(1) + b vi *(1) + b n ⎛ ⋅ exp⎜ a ∑i =1α i (| ∫ (ui + bˆ)(t )dt − ∫ (vi + bˆ)(t )dt | ui ( 0 ) + b vi ( 0 ) + b ⎝
kui *(1) kvi *(1) n u i (t / k ) dt − ∫ vi (t / k ) dt | ⋅ exp⎛⎜ a ∑i =1 α i (| ∫ ( 0 ) k u k vi ( 0 ) i ⎝ k ui ( 0 ) k vi ( 0 ) u i (t / k )dt − ∫ vi (t / k )dt |) ⎞⎟ +|∫ * ( 1 ) kui kvi *(1) ⎠
vi (0) +b ⎞ (ui + bˆ)(t )dt − ∫ (vi + bˆ)(t )dt |) ⎟ vi *(1) + b ⎠ 1 n ⎛ ⎞ = ⎜ ∑i =1α i ∫ r | (ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr ⎟ 0 ⎝ ⎠ u i *(1) v i *(1) n ⎛ ⋅ exp⎜ a ∑i =1α i (| ∫ (ui + bˆ)( s + b)ds − ∫ (vi + bˆ)( s + b)ds | ui ( 0 ) vi ( 0 ) ⎝
+|∫
u i (0) +b
u i *(1) + b
1 n ⎛ ⎞ = ⎜ | k | ∑i =1 α i ∫ r[| u i (r ) − vi ( r ) | + | u i ( r ) − vi (r ) |]dr ⎟ 0 ⎝ ⎠ ui *(1) vi *(1) n ⎛ ku i ( s )ds − ∫ kvi ( s )ds | ⋅ exp⎜ a ∑i =1 α i (| ∫ ui ( 0 ) vi ( 0 ) ⎝
⎞ (vi + bˆ)( s + b)ds |) ⎟ ⎠ 1 ⎛ n ⎞ = ⎜ ∑i =1α i ∫ r | (ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr ⎟ 0 ⎠ ⎝ u * ( 1 ) v * ( 1 ) n i i ⎛ ⋅ exp⎜ a ∑i =1α i (| ∫ (ui )( s )ds − ∫ (vi )( s )ds | ui ( 0 ) vi ( 0 ) ⎝ +|∫
+|∫
u i (0)
u i *(1)
u i (0) u i *(1)
(ui + bˆ)( s + b)ds − ∫
(ui )(s ) ds − ∫
vi (0) v i *(1)
v i ( 0)
+|∫
vi *(1)
⎛ n ⎞ = ⎜ ∑i =1α i ∫ r[| ( ku ) i (r ) − ( kv) i ( r ) | + | ( ku) i ( r ) − ( kv) i ( r ) |]dr ⎟ 0 ⎝ ⎠ ( ) * ( 1 ) ( ) * ( 1 ) ku kv n i i ⎛ ( ku ) i (t ) dt − ∫ ( kv) i (t ) dt | ⋅ exp⎜ a ∑i =1α i (| ∫ ( ku ) i ( 0 ) ( kv ) i ( 0 ) ⎝ ( ku ) i *(1)
( ku ) i (t ) dt − ∫
( kv ) i ( 0 ) ( kv ) i *(1)
⎞ ( kv) i (t ) dt |) ⎟ ⎠
1 ⎛ n ⎞ = ⎜ ∑i =1α i ∫ r[| k u i ( r ) − k vi ( r ) | + | k u i ( r ) − k vi ( r ) |]dr ⎟ 0 ⎝ ⎠
kui *(1) kvi *(1) n ⎛ u i (t / k )dt − ∫ vi (t / k ) dt | ⋅ exp⎜ a ∑i =1α i (| ∫ k ui ( 0 ) k vi ( 0 ) ⎝
+|∫
k u i (0) kui *(1)
u i (t / k ) dt − ∫
k vi ( 0) kvi *(1)
⎞ vi (t / k ) dt |) ⎟ ⎠
1 n ⎛ ⎞ = ⎜ k ∑i =1 α i ∫ r[| u i ( r ) − vi ( r ) | + | u i ( r ) − vi (r ) |]dr ⎟ 0 ⎝ ⎠
ui *(1) vi *(1) n ⎛ ⋅ exp⎜ ka∑i =1 α i (| ∫ u i ( s ) ds − ∫ vi ( s ) ds | ui ( 0 ) vi ( 0 ) ⎝ u i (0) vi ( 0) ⎞ +|∫ u i ( s )ds − ∫ vi ( s )ds |) ⎟ ui *(1) vi *(1) ⎠ =| k | ⋅[ Λ α ,|k |a (u, v)]
For any k < 0 , we have Λ α ,a (ku, kv) ⎛ ⎞ = ⎜ ∑i =1 α i ∫ r[| (ku ) i (r ) − (kv) i (r ) | + | (ku ) i (r ) − (kv) i (r ) |]dr ⎟ 0 ⎝ ⎠ n
1
( ku ) i *(1) ( kv ) i *(1) n ⎛ (ku ) i (t )dt − ∫ (kv) i (t )dt | ⋅ exp⎜ a ∑i =1 α i (| ∫ ( ku ) i ( 0 ) ( kv ) i ( 0 ) ⎝
+|∫
( ku ) i ( 0 ) ( ku ) i *(1)
⎞ kvi ( s) ds |) ⎟ ⎠
If k = 0 , it is obvious that Λ α ,a (ku, kv) =| k | Λ α ,|k |a (u, v)
1
+|∫
vi ( 0) vi *(1)
ui *(1) vi *(1) n ⎛ ⋅ exp⎜ | k | ⋅a ∑i =1α i (| ∫ u i ( s) ds − ∫ vi ( s )ds | ui ( 0 ) vi ( 0 ) ⎝ u i (0) vi (0) ⎞ +|∫ u i ( s )ds − ∫ vi ( s )ds |) ⎟ ui *(1) vi *(1) ⎠ =| k | ⋅[Λ α ,|k |a (u , v)]
= Λα , a (u , v)
( ku ) i ( 0 )
ku i ( s )ds − ∫
1 n ⎛ ⎞ = ⎜ | k | ∑i =1 α i ∫ r[| u i (r ) − vi ( r ) | + | u i ( r ) − vi (r ) |]dr ⎟ 0 ⎝ ⎠
⎞ (vi )( s) ds |) ⎟ ⎠
i.e., conclusion (4) holds. For any k > 0 , we have Λ α , a ( ku, kv)
ui (0) ui *(1)
(ku ) i (t )dt − ∫
( kv ) i ( 0 ) ( kv ) i *(1)
⎞ (kv) i (t )dt |) ⎟ ⎠
1 ⎛ n ⎞ = ⎜ ∑i =1 α i ∫ r[| k u i (r ) − k vi (r ) | + | k u i (r ) − k vi (r ) |]dr ⎟ 0 ⎝ ⎠
holds, so conclusion (5) holds. Therefore, the proof of the theorem is completed. At the end of the section, combining the definitions of Δ α ,a and Λ α ,a , we give the following definition of difference value Γα ,a .
Definition 4.3. Let u, v ∈ L( E n ) with M α (u ), M α (v) ≥ 0 and M α (u ) + M α (v) ≠ 0 . We denote Γα , a (u , v) u i * (1) v i * (1) u i (0) v i ( 0) n ⎛ ⎞ exp⎜ a2 ∑ i =1 α i (| ∫ ui (t )dt − ∫ vi (t )dt |+ | ∫ ui (t )dt − ∫ vi (t )dt |) ⎟ ui ( 0) vi ( 0) u i * (1) v i * (1) ⎝ ⎠ = a1 1 ⎛ n ⎞ ⎜ ∑i =1 α i ∫ 0 r[ui (r ) + vi (r ) + ui (r ) + vi (r )]dr ⎟ ⎝ ⎠ 1 ⎛ n ⎞ ⋅ ⎜ ∑i =1 α i ∫ r[| (ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr ⎟ 0 ⎝ ⎠
and call Γα ,a (u , v) a difference value of u and v (with respect to α and a = (a1 , a2 ) ), where, α = (α1,α2 ,L,αn ) ∈ Rn with
∑
n
i =1
α i = 1 and α i ≥ 0 ( i = 1,2,L, n ), and a = (a1 , a2 ) ∈ (0,+∞)
1 1 1 × (0,+∞) . We denote Γ a (u , v) = Γα ,a (u, v) as α = ( , ,L, ) , n n n and Γa (u , v) = Γ1,a (u , v) as u, v ∈ E .
Likewise, we have the following properties about the difference value Γα ,a . Property 4.3.
Let u, v ∈ L( E n ) with M α (u ), M α (v) ≥ 0
and M α (u ) + M α ( v ) ≠ 0 , α = (α 1 , α 2 , L , α n ) ∈ R n
∑
n i =1
with
α i = 1 and α i > 0 ( i = 1,2, L , n ), and a = ( a1 , a 2 )
∈ ( 0 , +∞ ) × ( 0 , +∞ ) . Then
(1) Γα ,a (u, v) ≥ 0 ; (2) Γα ,a (u, v) = 0 if and only if u = v ; (3) Γα ,a (u , v) = Γα ,a (v, u ) ;
TFS-2008-0342
11
(4) Γα ,a (u + bˆ, v + bˆ) ≤ Γα ,a (u, v) for any b ∈ [0,+∞) ; (5) Γα ,a (ku, kv) = k
2− a1
μ
means
Γα ,b (u, v ) for any k ∈ [0,+∞) , where
i j
=
1 mi
∑
mi k =1
c kji
1 m ∑ i (ckji − μ ij ) 2 mi − 1 k =1
and
standard
deviations
b = (a1 , ka2 ) .
σ ij =
Proof. The proofs of the properties can be completed similarly with the proofs of Property 4.1 and 4.2, respectively, so we omit it.
character values of Ci ( i = 1,2, L , l ), respectively. We
V. PATTERN RECOGNITION BASED ON METRICS AND DIFFERENCE VALUES
In Section 3 and 4, we discussed metrics and difference values on L( E n ) . In this section, we put forward an algorithmic version of pattern recognition in an imprecise or uncertain environment based on the metrics and difference values defined by us, and give examples to show the application (see Example 5.1) and rationality (see Example 5.2) of the method. Consider a problem to identify an object (denoted by O ) belonging to some one of l classes (denoted by C1 , C 2 , L, Cl ) in an imprecise or uncertain environment. Let the objects have n characteristics. Since the problem discussed by us take on some imprecise or uncertain attributes, it is unsuitable (see Remark 2.1) that we use a crisp n − dimensional vector (i.e., a standard n − dimensional real number vector) to express the n character values of Ci ( i =1,2,L, l ) or O . Therefore, using the method of statistics, we construct l + 1 fuzzy n − cell numbers to express the n character values of C1 , C 2 , L , Cl and O , respectively, and then put forward an algorithmic version of pattern recognition based on the metrics or the difference values defined by us. Algorithmic version of pattern recognition based on metrics The first step: Depending on the practicality, we first find out one domain of the jth character value of Ci for each i ( i = 1,2, L , l ) and
( j = 1,2, L , n ) of the n
construct triangular model 1 − dimensional fuzzy numbers u ij ( i = 1,2, L , l , j = 1,2, L , n ) as the following: ⎧ x − ( μ ij − 2σ ij ) ⎪ 2σ ij ⎪ i ⎪⎪ ( μ j + 2σ ij ) − x u ij ( x ) = ⎨ 2σ ij ⎪ ⎪ ⎪0 ⎩⎪
if
x ∈ [ μ ij − 2σ ij , μ ij ] I D ij
if
x ∈ ( μ ij , μ ij + 2σ ij ] I D ij
if
x ∉ [ μ ij − 2σ ij , μ ij + 2σ ij ] I D ij
( i = 1,2, L , l , j = 1,2, L , n ) or construct Gaussian model one-dimensional fuzzy numbers v ij ( i = 1,2,L, l , j = 1,2,L, n ) as the following ⎧ ( x − μ ij ) 2 − exp( ) if ⎪ 2 v ij ( x) = ⎨ 2σ ij ⎪0 if ⎩
x ∈ D ij x ∉ D ij
( i = 1,2, L , l , j = 1,2, L , n ) We construct fuzzy n − cell numbers u i = (u1i , u 2i , L , u ni )
( u i ( x1 , x2 ,L, xn ) = min{u1i ( x1 ), u 2i ( x2 ),L, u ni ( xn )} ) and vi = (v1i , v2i ,L, vni )
( v i ( x1 , x2 ,L, xn ) = min{v1i ( x1 ), v2i ( x2 ),L, vni ( xn )} ),
i = 1,2, L , l , and use u i or v i to express the ith class Ci
i characters of the mi samples, denote cmi = (cmi 1 , cmi 2 ,L, cmn )
( i = 1,2, L , l ). The second step: For the object O to be recognized, taking t samples in O , we can gain t classes of data about the n characters of O as the following: o11 o12 L o1n o21 o22 L o2 n O: M M M ot1 ot 2 L otn
( m = 1,2, L , mi ), i.e., we gain the following tables:
We work out the following means (denoted by o1 , o2 ,L , on )
j ( j = 1,2, L , n ), and denote the said domain by D ij .
We arbitrarily take mi samples in Ci ( i = 1,2, L , l ), and gain mi measure values (denoted by c1i , c2i L , cmi i ) of the n
C1 :
1 11 1 21
1 12 1 22
1 1n 1 2n
c
c
L
c
c
c
L
c
c
M
M
M
1 m1 1
1 m1 2
1 m1n
L Cl :
c
l 11 l 21
L c l 12 l 22
C2 :
L
c
c
c
L
c
M
M
M
l ml 1
l ml 2
l ml n
L c
2 1n 2 2n
c
c
L
c
c
c
L
c
M
c
c
c
2 12 2 22
2 m2 1
c
M
M
2 m2 2
2 m2 n
L c
l 1n l 2n
c
c
2 11 2 21
For Ci ( i = 1,2, L , l ), we directly work out the following
and standard deviations (denoted by s1 , s 2 , L , s n ) of the n character values of O : 1 t oi = ∑k =1 oki ( i = 1,2,L, n ) t and 1 t (oki − oi ) 2 ( i = 1,2,L, n ) ∑ t − 1 k =1 We construct triangular model one-dimensional fuzzy numbers wi ( i = 1,2,L, n ) as the following: si =
TFS-2008-0342
12
⎧ x − ( oi − 2 s i ) ⎪ 2si ⎪ ⎪ ( oi + 2 s i ) − x wi ( x ) = ⎨ 2si ⎪ ⎪ ⎪0 ⎩
Taking proper α = (α 1 , α 2 , L , α n ) with
if
x ∈ [oi − 2si , oi ] I (U lj =1 Di j )
if
x ∈ (oi , oi + 2si ] I (U lj =1 Di j )
if
x ∉ [oi − 2si , oi + 2 si ] I (U j =1 Di ) l
⎧ ( x − oi ) ) if ⎪vi′ ( x) = exp(− 2 si2 ⎨ ⎪ 0 if ⎩
x ∈ U lj =1 Di j x ∉ U lj =1 Di j
j
( i = 1,2,L, n )
w′2 ,L, w′n ) ( w′( x1 , x2 ,L, xn ) = min{w1′ ( x1 ), w2′ ( x2 ),L, w′n ( xn )} ),
Γα , a ( w, u j ) w i (0) u ij ( 0 ) w i * (1) u ij * (1) ⎛ ⎞ n exp⎜⎜ a2 ∑ i =1 α i (| ∫ wi (t )dt − ∫ j uij (t ) dt |+ | ∫ wi (t )dt − ∫ j uij (t )dt |) ⎟⎟ w i * (1) u i * (1) wi ( 0) u i (0) ⎝ ⎠ = a1 1 ⎛ n ⎞ j j ⎜ ∑ i =1 α i ∫ 0 r[ wi (r ) + ui (r ) + wi (r ) + ui (r )]dr ⎟ ⎝ ⎠ 1 ⎞ ⎛ n ⋅ ⎜ ∑ i =1 α i ∫ r[| ( wi ( r ) − uij ( r ) | + | wi ( r ) − uij ( r ) |]dr ⎟ 0 ⎠ ⎝
( j = 1,2, L , l ) or w i′ * (1) v ij * (1) w ′ i (0) v ij ( 0 ) ⎛ ⎞ n exp⎜⎜ a 2 ∑ i =1 α i (| ∫ wi′ (t )dt − ∫ j vij (t )dt |+ | ∫ wi′ (t )dt − ∫ vij (t )dt |) ⎟⎟ w i′ ( 0 ) vi (0) w i′ * (1) v i * (1) ⎝ ⎠ = a1 1 ⎛ n ⎞ j j ⎜ ∑ i =1 α i ∫ 0 r[ wi′ (r ) + vi (r ) + wi′ (r ) + vi (r )]dr ⎟ ⎝ ⎠
and use w or w′ to express the object O . The third step:
1 ⎛ n ⎞ ⋅ ⎜ ∑ i =1 α i ∫ r[| ( wi′ (r ) − vij (r ) | + | wi′ (r ) − vij (r ) |]dr ⎟ 0 ⎝ ⎠
∑
n
i =1
α i = 1 and
α i ≥ 0 , i = 1,2,L, n , and p ≥ 1 , we compute the metrics ~ Dα , p ( w, u j )
=
or
1 1 n ( [r ∑i =1α i ⋅ (| wi (r ) − uij (r ) | + | wi (r ) − uij (r ) |)] p dr )1/ p 2 ∫0 ( j = 1,2, L , l ) ~ Dα , p ( w′, v j )
=
α i = 1 and
Γα , a ( w′, v j )
We construct fuzzy n − cell numbers w = (w1 , w2 ,L, wn ) ( w( x1 , x2 ,L, xn ) = min{w1 ( x1 ), w2 ( x2 ),L, wn ( xn )} ) and w′ = (w1′,
Taking proper α = (α1 , α 2 , L , α n ) with
n i =1
α i ≥ 0 , i = 1 , 2 , L , n , and a = ( a1 , a 2 ) ∈ (0,+∞ ) × (0,+∞ ) , we compute the difference values
( i = 1,2,L, n ) or construct Gaussian model one-dimensional fuzzy numbers wi′ ( i = 1,2,L, n ) as the following 2
∑
1 1 n ( [r ∑i =1α i ⋅ (| wi′( r ) − vij (r ) | + | wi′(r ) − vij (r ) |)] p dr )1/ p 2 ∫0 ( j = 1,2, L , l )
The fourth step: We choose u j0 in u1,u2 ,L,ul , or v j0′ in v1, v2 ,L, vl such that ~ ~ ~ ~ Dα , p ( w, u j ) = min{Dα , p ( w, u 1 ), Dα , p ( w, u 2 ), L , Dα , p ( w, u l )} 0
or ~ ~ ~ ~ Dα , p ( w′, v j′ ) = min{Dα , p ( w′, v1 ), Dα , p ( w′, v 2 ), L , Dα , p ( w′, v l )} 0
Then we can consider that object O belongs to the j0 th class C j0 , or belongs to the j0′ th class C j0′ . Remark 5.1. In the third and fourth steps of the above ~ method, we can have the metric Dˆ α , p replace the metric Dα , p , as a result of which, we can also set up a method based on the metric Dˆ α , p . Algorithmic version of pattern recognition based on difference values The first step and the second step: They are respectively same with the first step and the second step of the method of pattern recognition based on the metric, as above. The third step:
( j = 1,2, L , l ) The fourth step: We choose u j0 in u 1 , u 2 ,L, u l , or v j0′ in v1 , v 2 , L , v l such that Γα ,a ( w, u j ) = min{Γα ,a ( w, u 1 ), Γα ,a ( w, u 2 ), L , Γα ,a ( w, u l )} 0
or Γα ,a ( w′, v j′ ) = min{Γα ,a ( w′, v1 ), Γα ,a ( w′, v 2 ), L , Γα ,a ( w′, v l )} 0
Then we can consider that object O belongs to the j0 th class C j0 , or belongs to the j0′ th class C j0′ . Remark 5.2. In the third and fourth steps of the method above, we can have the difference value Δ α ,a or Λ α ,a replace the difference value Γα ,a , as a result of, we can also set up a method based on the difference value Δ α ,a or Λ α ,a . In order to be more obvious, we may use the following diagram to illustrate the methods set up by us. Measure values of the jth character of samples in C1, j=1,2,…,n
Measure values of the jth character of samples in Cl , j=1,2,…,n
Measure values of the ith character of samples in O, i=1,2,…,n
Mean μ 1j , j=1,2,…,n
Mean μ lj , j=1,2,…,n
Mean μi , i=1,2,…,n
Standard deviation σ 1j , j=1,2,…,n
Standard deviation σ lj , j=1,2,…,n
Stan-dard deviation σ i , i=1,2,…,n
1-dimensional fuzzy numbers: u1j , j=1,2,…,n
1-dimensional fuzzy numbers: ulj , j=1,2,…,n
1-dimensional fuzzy numbers: wi , i=1,2,…,n
Fuzzy n-cell number u1
Fuzzy n-cell number ul
Fuzzy n-cell number w
Difference value
Γα , a ( w, u1 )
The minimal
O
Difference value
Γα , a (w, ul 0 )
belongs to the j0th class C j
0
Γα , a ( w, ul )
TFS-2008-0342
13
Example 5.1. Suppose that some terrain consists of five different types of land based cover: C1 : Road; C2 : Farm or
Crop; C3 : Korean Pine accounts for the main part; C4 : Boreal and Broad-leaf Mixture Forest; C5 : Birch Forest. For the five types of land cover ( C1 , C2 , C3 , C4 , C5 ) and by using the four wave bands: MSS-4, MSS-5, MSS-6, MSS-7, we take 10 samples, and acquire the following data:
C5 :
MSS - 4 MSS - 5 MSS - 6 MSS - 7 Sample 1 Sample 2
C1 :
18.62 19.76
20 .71 17.01
58 .20 51 .02
26 .72 24 .32
MSS - 5 13 .45
MSS - 6 45 .76
MSS - 7 23 .64
Sample 2 Sample 3 Sample 4
16 .41 17 .04 16 .58
12 .62 12 .88
44 .71 44 .95
23 .19 22 .54
Sample 5 Sample 6 Sample 7
18 .21 17 .23 15 .89
13 .13 13 .90 13 .73 12 .98
43 .43 46 .38 45 .03 44 .32
22 .98 23 .41 23 .76 23 .15
Sample 8 Sample 9 Sample 10
17 .02 15 .97 17 .67
13 .25 12 .70 13 .34
45 .38 43 .96 46 .09
22 .61 23 .11 23 .62
Sample 3
18.24
19 .46
48 .12
26 .33
18.76
15 .95
56 .35
22 .89
Sample 5 Sample 6
18.96 19.90
18 .78 20 .13
45 .32 50 .82
28 .32 25.05
Sample 7
19.16
18 .58
52 .30
22 .21
and
Sample 8
19 .36
17 .32
55 .02
25.28
Sample 9 Sample 10
19 .36 18.48
17 .98 16 .46
48 .26 46 .63
23 .86 27 .46
1 10 ∑ (ckji − μ ij ) 2 ( i = 1,2,3,4,5 , j = 1,2,3,4 ) 10 − 1 k =1 we can work out the following means and standard deviations: μ11 = 19.06, σ 11 = 0.55 , μ 21 = 18.24, σ 21 = 1.58
By
μ ij =
1 10 i ∑ ckj ( i = 1,2,3,4,5 , j = 1,2,3,4 ) 10 k =1
σ ij =
22.43
μ 31 = 51.20, σ 31 = 4.28 , μ 41 = 25.24, σ 41 = 1.98
45.56
23.89
47.80 43.88
18.05 20.43
μ12 = 21.89, σ 12 = 2.88 , μ 22 = 24.68, σ 22 = 4.82
32.28
53.61
24.95
26.06 17.64
24.15 17.08
42.75 46.73
19.38 21.16
Sample 8
22.23
26.02
49.17
19.04
Sample 9
20.26
21.56
41.37
24.39
μ 34 = 42.41, σ 34 = 2.87 , μ 44 = 21.22, σ 44 = 1.50
Sample 10
25.42
30.52
50.46
22.59
μ15 = 17.01, σ 15 = 0.82 , μ 25 = 13.20, σ 25 = 0.42
Sample 1
24.51
28.13
52.39
Sample 2 Sample 3 Sample 4
19.37
23.41
23.12 18.93
24.82 18.87
Sample 5
21.36
Sample 6 Sample 7
Sample 1 Sample 2 Sample 3 Sample 4 C 3 : Sample 5 Sample 6 Sample 7 Sample 8 Sample 9 Sample 10
MSS - 4 15.01 15.60 15.82 14.90 16.10 13.80 15.90 16 .82 15 .50 15 .38
MSS - 5 MSS - 6 13.30 40.50 12.56 38.81 12.79 37.70 11.70 35.50 13.80 42.10 11.94 32.10 10 .98 30 .87 13 .67 37 .64 12 .58 36 .10 12 .48 34 .08
MSS - 7 19.37 16.35 18.16 14.75 20.75 15.54 14 .29 18 .62 18 .02 17 .45
MSS - 4 MSS - 5 MSS - 6 MSS - 7
C4 :
MSS - 4 18 .12
Sample 4
MSS - 4 MSS - 5 MSS - 6 MSS - 7
C2 :
Sample 1
Sample 1 Sample 2
17 .05 16 .09
13 .53 12 .03
43 .32 38 .65
22 .76 19 .47
Sample 3
15 .44
12 .87
42 .21
21 .20
Sample 4
16 .41
12 .58
40 .87
18 .63
Sample 5 Sample 6
16 .32 17 .21
13 .01 13 .58
46 .56 45 .98
23 .58 21 .33
Sample 7
15 .98
12 .14
41 .09
20 .23
Sample 8
16 .21
12 .68
42 .51
21 .65
Sample 9 Sample 10
15 .09 16 .38
12 .10 13 .31
38 .03 44 .87
21.01 22 .34
μ 32 = 47.37, σ 32 = 4.09 , μ 42 = 21.63, σ 42 = 2.39 μ13 = 15.46, σ 13 = 1.22 , μ 23 = 12.58, σ 23 = 0.88 μ 33 = 36.54, σ 33 = 3.55 , μ 43 = 17.33, σ 43 = 2.08 μ14 = 16.22, σ 14 = 0.64 , μ 24 = 12.78, σ 24 = 0.58
μ 35 = 45.00, σ 35 = 0.94 , μ 45 = 23.20, σ 45 = 0.42 Taking Di j = (0,+∞) ( according to ⎧ x − ( μ ij − 2σ ij ) ⎪ 2σ ij ⎪ ⎪⎪ ( μ ij + 2σ ij ) − x u ij ( x) = ⎨ 2σ ij ⎪ ⎪ ⎪0 ⎩⎪
i = 1,2,3,4 and j = 1,2,3,4,5 ) then
if
x ∈ [ μ ij − 2σ ij , μ ij ] I D ij
if
x ∈ ( μ ij , μ ij + 2σ ij ] I D ij
if
x ∉ [ μ ij − 2σ ij , μ ij + 2σ ij ] I D ij
( i = 1,2,3,4,5 , j = 1,2,3,4 ) we have ⎧ x − 17.96 ⎪ ⎪ 1.10 ⎪ 20.16 − x 1 u1 ( x) = ⎨ ⎪ 1.10 ⎪ 0 ⎩⎪ ⎧ x − 42.64 ⎪ ⎪ 8.56 ⎪ 59.76 − x 1 u3 ( x ) = ⎨ ⎪ 8.56 ⎪ 0 ⎪⎩
if if if if if if
⎧ x − 15.08 ⎪ ⎪ 3.16 ⎪ 21.40 − x 1 x ∈ (19.06, 20.16] , u2 ( x) = ⎨ ⎪ 3.16 ⎪ 0 x ∉ [17.96, 20.16] ⎩⎪
x ∈ [17.96, 19.06]
⎧ x − 21.28 ⎪ ⎪ 3.96 ⎪ 29.20 − x 1 x ∈ (51.20, 59.76] , u4 ( x) = ⎨ ⎪ 3.96 ⎪ 0 x ∉ [42.64, 59.76] ⎪⎩
x ∈ [42.64, 51.20]
if
x ∈ [15.08, 18.24]
if
x ∈ (18.24, 21.40]
if
x ∉ [15.08, 21.40]
if
x ∈ [21.28, 25.24]
if
x ∈ (25.24, 29.20]
if
x ∉ [21.28, 29.20]
TFS-2008-0342 ⎧ x − 16.13 ⎪ ⎪ 5.76 ⎪ 27.65 − x 2 u1 ( x) = ⎨ ⎪ 5.76 ⎪0 ⎩⎪
⎧ x − 39.19 ⎪ ⎪ 8.18 ⎪ 55.55 − x 2 u3 ( x) = ⎨ ⎪ 8.18 ⎪ 0 ⎩⎪ ⎧ x − 13.02 ⎪ ⎪ 2.44 ⎪ 17.90 − x 3 u1 ( x) = ⎨ ⎪ 2.44 ⎪0 ⎪⎩ ⎧ x − 29.44 ⎪ ⎪ 7.10 ⎪ 43.64 − x 3 u3 ( x) = ⎨ ⎪ 7.10 ⎪ 0 ⎪⎩
⎧ x − 14.94 ⎪ ⎪ 1.28 ⎪ 17.50 − x 4 u1 ( x) = ⎨ ⎪ 1.28 ⎪0 ⎪⎩ ⎧ x − 36.67 ⎪ ⎪ 5.74 ⎪ 48.15 − x 4 u3 ( x ) = ⎨ ⎪ 5.74 ⎪ 0 ⎪⎩
⎧ x − 15.37 ⎪ ⎪ 1.64 ⎪ 18.65 − x u15 ( x) = ⎨ ⎪ 1.64 ⎪0 ⎪⎩ ⎧ x − 43.02 ⎪ ⎪ 1.98 ⎪ 46.98 − x 5 u3 ( x) = ⎨ ⎪ 1.98 ⎪ 0 ⎪⎩
if if if
if if if if if if
if if if
if if if
if if if
if if if
if if if
14 ⎧ x − 15.04 ⎪ ⎪ 9.64 ⎪ 34.32 − x 2 x ∈ (21.89, 27.65] , u2 ( x) = ⎨ ⎪ 9.64 ⎪ 0 x ∉ [16.13, 27.65] ⎩⎪
x ∈ [16.13, 21.89]
⎧ x − 16.85 x ∈ [39.19, 47.31] ⎪ ⎪ 4.78 ⎪ 26.41 − x 2 x ∈ (47.31, 55.55] , u4 ( x) = ⎨ ⎪ 4.78 ⎪0 x ∉ [39.19, 55.55] ⎩⎪ ⎧ x − 10.82 ⎪ ⎪ 1.76 ⎪ 14.34 − x 3 x ∈ (15.46, 17.90] , u2 ( x) = ⎨ ⎪ 1.76 ⎪0 x ∉ [13.02, 17.90] ⎪⎩
x ∈ [13.02, 15.46]
⎧ x − 13.17 ⎪ ⎪ 4.16 ⎪ 21.49 − x 3 x ∈ (36.54, 43.64] , u4 ( x) = ⎨ ⎪ 4.16 ⎪ 0 x ∉ [ 29.44, 43.64] ⎪⎩
x ∈ [ 29.44, 36.54]
⎧ x − 11.62 x ∈ [14.94, 16.22] ⎪ ⎪ 1.16 ⎪ 13.94 − x 4 x ∈ (16.22, 17.50] , u2 ( x) = ⎨ ⎪ 1.16 ⎪0 x ∉ [14.94, 17.50] ⎪⎩ ⎧ x − 18.22 ⎪ ⎪ 3.00 ⎪ 24.22 − x 4 , ( ) u x = ⎨ x ∈ (42.41, 48.15] 4 ⎪ 3.00 ⎪ 0 x ∉ [36.67, 48.15] ⎪⎩
x ∈ [36.67, 42.41]
⎧ x − 12.36 ⎪ ⎪ 0.84 ⎪ 14.04 − x x ∈ (17.01, 18.65] , u25 ( x) = ⎨ ⎪ 0.84 ⎪0 x ∉ [15.37, 18.65] ⎪⎩
x ∈ [15.37, 17.01]
⎧ x − 22.36 ⎪ ⎪ 0.84 ⎪ 24.04 − x 5 , = u x ( ) x ∈ (45.00, 46.98] ⎨ 4 ⎪ 0.84 ⎪ 0 x ∉ [43.02, 46.98] ⎪⎩
x ∈ [43.02, 45.00]
if
x ∈ [15.04, 24.68]
u33 (r ) = 7.10r + 29.44, u33 (r ) = 43.64 − 7.10r
if
x ∈ (24.68, 34.32]
if
x ∉ [15.04, 34.32]
u 43 (r ) = 4.16r + 13.17, u 43 (r ) = 21.49 − 4.16r u14 (r ) = 1.28r + 14.94, u14 (r ) = 17.50 − 1.28r
if
x ∈ [16.85, 21.63]
if
x ∈ (21.63, 26.41]
if
x ∉ [16.85, 26.41]
u34 (r ) = 5.74r + 36.67, u 34 (r ) = 48.15 − 5.74r
if
x ∈ [10.82, 12.58]
u 44 (r ) = 3.00r + 18.22, u 44 (r ) = 24.22 − 3.00r
if
x ∈ (12.58, 14.34]
u15 (r ) = 1.64r + 15.37, u15 (r ) = 18.65 − 1.64r
if
x ∉ [10.82, 14.34]
u 25 (r ) = 0.84r + 12.36, u 25 (r ) = 14.04 − 0.84r
if
x ∈ [13.17, 17.33]
u35 (r ) = 1.98r + 43.02, u 35 (r ) = 46.98 − 1.98r
if
x ∈ (17.33, 21.49]
if
x ∉ [13.17, 21.49]
if
x ∈ [11.62, 12.78]
if
x ∈ (12.78, 13.94]
u i ( x1 , x2 , x3 , x4 ) = min{u1i ( x1 ), u 2i ( x2 ), u3i ( x3 ), u 4i ( x4 )}
if
x ∉ [11.62, 13.94]
( x1 , x2 , x3 , x4 ) ∈ R 4 , i = 1,2,3,4,5
if
x ∈ [18.22, 21.22]
if
x ∈ (21.22, 24.22]
if
x ∉ [18.22, 24.22]
if
x ∈ [12.36, 13.20]
if
x ∈ (13.20, 14.04]
if
x ∉ [12.36, 14.04]
if
x ∈ [22.36, 23.20]
if
x ∈ ( 23.20, 24.04]
if
x ∉ [ 22.36, 24.04]
u 24 (r ) = 1.16r + 11.62, u 24 (r ) = 13.94 − 1.16r
u 45 (r ) = 0.84r + 22.36, u 45 (r ) = 24.04 − 0.84r
Thus the fuzzy 4 − cell numbers u i = ( u 1i , u 2i , u 3i , u 4i ) , i = 1, 2 , 3 , 4 , 5 , ie.,
can be used to represent Ci , i = 1,2,3,4,5 , respectively. Using the same four wave bands: MSS-4, MSS-5, MSS-6, MSS-7, we now proceed to examine some zone (i.e., object, denoted by O ) 12 times, stochastically, at various times or positions, or using various viewers, and obtain the following data: Sample 1 Sample 2 Sample 3 Sample 4 O
and for r ∈ [0,1] u11 (r ) = 1.10r + 17.96, u11 (r ) = 20.16 − 1.10r u 12 (r ) = 3.16r + 15.08, u 12 (r ) = 21.40 − 3.16r u31 (r ) = 8.56r + 42.64, u31 (r ) = 59.76 − 8.56r u 14 (r ) = 3.96r + 21.28, u 14 (r ) = 29.20 − 3.96r
MSS − 4 MSS − 5 MSS − 6 MSS − 7 18.31 13.13 45.41 23.52 17.21 12.71 45.67 23.49 17.00 12.57 44.82 22.56 16.32 13.06 43.53 22.77
Sample 5 Sample 6 Sample 7 Sample 8
18.01 16.11 16.21 17.13
13.91 13.59 12.56 13.32
45.31 46.97 43.96 46.01
23.31 23.61 23.14 22.84
Sample 9 Sample 10 Sample 11 Sample 12
16.20 17.51 17.02 18.02
12.34 13.34 13.78 13.25
43.81 45.89 45.21 44.56
23.13 23.52 23.18 23.23
We can work out the following means and standard deviations: o1 = 17.09, s1 = 0.77 , o2 = 13.17, s 2 = 0.50 o3 = 45.10, s3 = 1.01 , o4 = 23.19, s 4 = 0.33
u 22 (r ) = 9.64r + 15.04, u 22 (r ) = 34.32 − 9.64r
So we can obtain the fuzzy 4 − cell number o = (o1 , o2 , o3 , o4 ) , i.e., o( x1 , x2 , x3 , x4 ) = min{o1 ( x1 ), o2 ( x2 ), o3 ( x3 ), o4 ( x 4 )}
u32 (r ) = 8.18r + 39.19, u 32 (r ) = 55.55 − 8.18r
( x1 , x2 , x3 , x4 ) ∈ R 4
u12 (r ) = 5.76r + 16.13, u12 (r ) = 27.65 − 5.76r
u 42 (r ) = 4.78r + 16.85, u 42 (r ) = 26.41 − 4.78r u (r ) = 2.44r + 13.02, u (r ) = 17.90 − 2.44r 3 1
3 1
u 23 (r ) = 1.76r + 10.82, u 23 (r ) = 14.34 − 1.76r
that can be used to represent O , where, ⎧ x − 15.55 ⎪ ⎪ 1.54 ⎪ 18.63 − x o1 ( x) = ⎨ ⎪ 1.54 ⎪0 ⎪⎩
if if if
⎧ x − 12.17 ⎪ 1.00 ⎪ ⎪ 14.17 − x x ∈ (17.09, 18.63] , o2 ( x) = ⎨ ⎪ 1.00 ⎪0 x ∉ [15.55, 18.63] ⎪⎩
x ∈ [15.55, 17.09]
if
x ∈ [12.17, 13.17]
if
x ∈ (13.17, 14.17]
if
x ∉ [12.17, 14.17]
TFS-2008-0342 ⎧ x − 43.08 ⎪ ⎪ 2.02 ⎪ 47.12 − x o3 ( x) = ⎨ ⎪ 2.02 ⎪ 0 ⎩⎪
if if if
15 ⎧ x − 22.53 if ⎪ ⎪ 0.66 ⎪ 23.79 − x x ∈ (45.10, 47.12] , o4 ( x) = ⎨ if ⎪ 0.66 ⎪ x ∉ [ 43.08, 47.12] 0 if ⎩⎪
x ∈ [ 43.08, 45.10]
x ∈ [22.53, 23.19] x ∈ (23.19, 23.79] x ∉ [22.53, 23.79]
and for r ∈ [0,1] . o1 (r ) = 1.54r + 15.55, o1 (r ) = 18.63 − 1.54r o2 (r ) = 1.00r + 12.17, o2 (r ) = 14.17 − 1.00r o3 (r ) = 2.02r + 43.08, o3 (r ) = 47.12 − 2.02r o4 (r ) = 0.66r + 22.53, o4 (r ) = 23.79 − 0.66r 1 1 1 1 1 1 Taking α = ( , , , ) and a = ( , ) , By 5 5 4 4 4 4
although d (C , A) < d (C , B) , the difference of d (C , A) and d (C , B ) is small. Furthermore, from the statistical data, we can see that the ferrous quantities of minerals coming from A are more coincident, but B and C are not. It is almost impossible that some minerals in C come from A , such as, minerals with ferrous quantities 76.02 and 28.56 . So we may judge that C comes from B . If we use fuzzy 1 − cell numbers to represent A , B and C , then we have ⎧ x − 6.75 ⎪ ⎪ 3.34 ⎪ 13.43 − x A( x) = ⎨ ⎪ 3.34 ⎪0 ⎩⎪
Γα , a (u, v)
4
1
0
we can obtain Γα , a (o, u1) = 3.36 , Γα ,a (o, u 2 ) = 10.34 , Γα , a (o, u 3 ) = 4.27 , Γα ,a (o, u 4 ) = 1.11 , Γα ,a (o, u 5 ) = 0.03 . From Γα , a (o, u 5 )
= min{Γα , a (o, u1 ), Γα , a (o, u 2 ), Γα , a (o, u 3 ), Γα , a (o, u 4 ), Γα , a (o, u 5 )} , we
if if
⎧ x + 14.83 ⎪ ⎪ 54.86 ⎪ 94.89 − x x ∈ (10.09, 13.43] , B( x) = ⎨ ⎪ 54.86 ⎪ 0 x ∉ [6.75, 13.43] ⎩⎪ x ∈ [6.75, 10.09]
⎧ x + 13.20 ⎪ ⎪ 37.72 ⎪ 62.24 − x C ( x) = ⎨ ⎪ 37.72 ⎪ 0 ⎪⎩
u i * (1) v i *(1) u i (0) v i (0) 4 ⎛ ⎞ exp⎜ a2 ∑i =1α i (| ∫ ui (t ) dt − ∫ vi (t ) dt |+ | ∫ ui (t ) dt − ∫ vi (t )dt |) ⎟ ui (0) vi ( 0 ) u i * (1) v i *(1) ⎝ ⎠ = a1 ⎛ 4 α 1r[u ( r ) + v ( r ) + u ( r ) + v (r )]dr ⎞ ⎜ ∑i =1 i ∫ 0 i ⎟ i i i ⎝ ⎠
⋅ ∑i =1α i ∫ r[| (ui (r ) − vi (r ) | + | ui (r ) − vi (r ) |]dr
if
if
x ∈ [0, 40.03]
if
x ∈ (40.03, 94.89]
if
x ∉ [0, 94.89]
if
x ∈ [0, 24.52]
if
x ∈ (24.52, 62.24]
if
x ∉ [0, 62.24]
⎧54.86r − 14.83 if r ∈ (0.27,1] A(r ) = 3.34r + 6.75 , B (r ) = ⎨ , if r ∈ [0,0.27] ⎩0 ⎧37.72r − 13.20 if r ∈ (0.35,1] , A(r ) = 13.43 − 3.34r , C (r ) = ⎨ if r ∈ [0,0.35] ⎩0
B (r ) = 94.89 − 54.86r , C(r ) = 62.24 − 37.72r , so Γ( 0.1, 0.1) (C, A)
know that O belongs to C5 , i.e., the zone measured by us is covered by Birch Forest. Remark 5.3. Although, only for the example, using mean vectors ( μ1i , μ 2i , μ 3i , μ 4i ) , i = 1,2,3,4,5 and (o1 , o 2 , o 4 , o 4 ) to
15 .38 14 .01 = 106 .39 > 21 .77 = e1.178 = Γ( 0.1, 0.1) (C , B ) . 4.019 3 .487 Thus we can affirm that C comes from B .
represent respectively C1 , C 2 , C3 , C 4 and O , we perhaps also
VI. CONCLUSION
judge that O belongs to C5 by the usual Euclidean metrics, we still emphasize that using fuzzy n − cell numbers to deal with imprecise or uncertain quantities is better than using crisp n − dimensional vectors. The following example (to simplify and shorten the problem, we only consider 1 − dimensional case) will show this. Example 5.2. The following two classes of ferrous quantities (kilogram per hundred kilogram) of minerals come respectively two different mine areas (denoted by A and B ).
In this paper, we suggest using fuzzy n − cell numbers to represent imprecise or uncertain multichannel digital signals, and put forward a method (see the first or second step of the algorithmic version in Section 5, or see Example 2.1) of constructing such fuzzy n − cell numbers. Although the metrics D and DL have been studied formerly in [14,15], in view of
A : 10.20, 11.76, 8.31, 9.02, 9.63,
8.33, 11.36, 12.30, 12.03, 7.98
B : 52.33, 79.34, 34.51, 62.34, 82.26, 28.36, 17.37, 25.32, 10.11, 8.34
Suppose that one group (denoted by C ) of minerals comes from the one of A and B . The problem to be solved is to identify if C comes from A or B . We take samples, and acquire the following data for C : C :18.20, 20.31, 76.02, 9.36, 28.56, 23.32, 15.36, 20.51, 13.27, 20.32
We can work out: μ A = 10.09 , σ A = 1.67 , μ B = 40.03 ,
σ B = 27.43 , μ C = 24.52 and σ C = 18.86 . If we use crisp means to represent A , B and C , then we have A = 10.09 , B = 40.03 , C = 24.52 and d (C , A) = 14.43 < 15.51 = d (C , B) . If we regard d (C , A) < d (C , B) as evidence, we can draw a conclusion that C comes from A . However, the conclusion does not accord with fact. We should note that
= e 3.27
the roughness of D and DL , we define two new metrics on fuzzy n − cell number space in order that they can better characterize the degree of the difference of two objects in some imprecise or uncertain environment, and we study their properties (Section 3). In some applications, metrics are unsuitable for use in finding the difference of two fuzzy n − cell numbers, so we introduce the concepts of difference values Δ α ,a , Λ α ,a and Γα ,a , study their properties, and show the rationality for their use in characterizing the degree of the difference of two fuzzy n − cell numbers by Remarks and examples (see Section 4). Finally, in Section 5, we put forward an algorithmic version of pattern recognition in an imprecise or uncertain environment based on the metrics and difference values defined by us, and give examples to show the application (see Example 5.1) and rationality (see Example 5.2) of the methods.
TFS-2008-0342
16 ACKNOWLEDGEMENT
The authors would like to thank the Associate Editor and reviewers for their valuable suggestions and comments which have greatly improved the presentation of the paper. This work is partially supported by the Engineering and Physical Sciences Research Council, UK (EP/F029195). REFERENCES [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15]
[16] [17] [18] [19] [20] [21]
D.Butnariu, “Methods of solving optimization problems and linear equations in the space of fuzzy vectors,” Libertas Math, vol. 17 , pp. 1–7, 1997. S.S.L.Chang, and L.A.Zadeh, “On fuzzy mappings and control,” IEEE Trans. Syst. Man Cybernet, vol. 2, pp. 30–34, 1972. P. Diamond, and P. Kloeden, Metric Space of Fuzzy Sets. World Scientific, Singapore, 1994. D. Dubois, and H. Prade, “Operations on fuzzy numbers,” Internat. J. Systems Sci., vol. 9, pp. 613–626, 1978. D. Dubois, and H. Prade, “Towards fuzzy differential calculus—Part I, II, III,” Fuzzy Sets and Systems, vol. 8, pp. 1–17, 105–116, 225–234, 1982. R.Goetschel, W.Voxman, “Elementary calculus,” Fuzzy Sets and Systems, vol. 18, pp. 31–43, 1986. O.Kaleva, “Fuzzy differential equations,” Fuzzy Sets and Systems, vol. 24, pp. 301–317, 1987. A.Kaufmann, M.M.Gupta, Introduction to Fuzzy Arithmetic-Theory and Applications. Van Nostrand, Reinhold, NY, 1985. M.Miyakawa, K.Nakamura, J.Ramik, I.G.Rosenberg, “Joint canonical fuzzy numbers,” Fuzzy Sets and Systems, vol. 53, pp. 39–49, 1993. K.Nakamura, “Canonical fuzzy number of dimension two and fuzzy utility difference for understanding preferential judgements,” Inform. Sci., vol. 50, pp. 1–22, 1990. M.L.Puri, D.A.Ralescu, “Differentials for fuzzy functions,” J.Math. Anal.Appl., vol. 91, pp. 552–558, 1983. J.Ramik, K.Nakamura, “Canonical fuzzy numbers of dimension two,” Fuzzy Sets and Systems, vol. 54, 167–181, 1993. S.Theodoridis, K.Koutroumbas, Pattern Recognition (Third Edition). Elsevier, USA, 2006. G.Wang, Y.Li, C.Wen, “On fuzzy n − cell numbers and n − dimension fuzzy vectors,” Fuzzy Sets and Systems, vol. 158, pp. 71–84, 2007. G.Wang, C.Wu, “Fuzzy n − cell numbers and the differential of fuzzy n − cell number value mappings,” Fuzzy Sets and Systems, vol. 130, pp. 367–381, 2002. G.Wang, C.Wu, “The integral over a directed line segment of fuzzy mapping and its applications,” International Journal of Uncertainty, Fuzzyness and Knowledge-Based Systems, vol. 12, pp. 543–556, 2004. G.Wang, C.Zhao, “The measurability and (K)-integrability about fuzzy n-cell number value mappings,” The Journal of Fuzzy Mathematics, vol. 13, pp. 117–128, 2005. C.Wu, G.Wang, “convergence of sequences of fuzzy numbers and fixed point theorems for inceasing fuzzy mappings and application,” Fuzzy Sets and Systems, vol. 130, pp. 383–390, 2002. T. Wang, Y. Chen and S. Tong, “Fuzzy reasoning models and algorithms on type-2 fuzzy sets”, Int. J. of Innovative Computing Information and Control, vol.4, no.10, pp. 2451-2460, 2008. M. Shimakawa, “A proposal of extension fuzzy reasoning method”, Int. J. of Innovative Computing Information and Control, vol.4, no.10, pp. 2603-2615, 2008. Kimiaki Shinkai, Decision Analysis of Fuzzy Partition Tree Applying Fuzzy Theory, Int. J. of Innovative Computing Information and Control, vol.4, no.10, pp.2581-2594, 2008.
Guixiang Wang received the BSc degree in Mathematics from Hebei Normal University, China in 1982; the MSc degree in Mathematics from Hebei University, China in 1989; and the PhD degree in Mathematics from Harbin Institute of Technology, China in 2002. He worked as an associate professor in Agricultural University of Hebei, China (1992-1999). He joined in Hangzhou Dianzi University, China in 2002 as a professor. He was a post doctoral fellow in Harbin Engineering University, China (2003-2005), and a visiting senior fellow in the University of Glamorgan, UK (2007-2008). Dr Wang’s recent research interests include fuzzy set theory and application, nonlinear systems, signal and information processing, data fusion, pattern recognition.
Peng Shi received the BSc degree in Mathematics from Harbin Institute of Technology, China in 1982; the ME degree in control theory from Harbin Engineering University, China in 1985; the PhD degree in electrical engineering from the University of Newcastle, Australia in 1994; and the PhD degree in mathematics from the University of South Australia in 1998. He was awarded the degree of Doctor of Science (Higher Doctorate) by the University of Glamorgan, UK, in 2006. Dr Shi worked as a lecturer in Heilongjiang University, China (1985-1989), in the University of South Australia (1997-1999) and a senior scientist in Defence Science and Technology Organisation, Department of Defence, Australia (1999-2005). He joined in the University of Glamorgan, UK, in 2004 as a professor. Dr Shi’s research interests include robust control and filtering, fault detection techniques, Markov decision processes, and optimization techniques. He has published a number of papers in these areas. In addition, Dr Shi is a co-author of the two research monographs, Fuzzy Control and Filtering Design for Uncertain Fuzzy Systems (Berlin, Springer, 2006), and Methodologies for Control of Jump Time-Delay Systems (Boston, Kluwer, 2003). Currently, Dr Shi serves as Editor-in-Chief of Int. J. of Innovative Computing, Information and Control, and Associate Editor for a number of other journals, such as IEEE Trans on Automatic Control, IEEE Trans on Systems, Man and Cybernetics, Part B, and IEEE Trans on Fuzzy Systems. Dr Shi is a Fellow of Institute of Mathematics and its Applications (UK), and a Senior Member of IEEE.