《1 引言》
1 引言
灰系统理论是我国著名学者邓聚龙教授1982年创立的一门新兴横断学科[1 ] 。灰数尤其是区间灰数表征及其运算问题, 在灰系统理论的研究与应用过程中起到了重要作用。然而, 目前的区间灰数的表征及其运算存在着较为严重的问题[2 ,3 ,4 ,5 ,6 ,7 ,8 ] 。例如, 在某种情况下, 它会对计算结果灰度产生不正常的放大。
例1 若给定区间灰数a 11 (g ) =1, a 12 (g ) =5, a 21 (g ) =[2 ,3 ] , a 22 (g ) =[0, 1], 其中g 表示灰数 (文献[1 ] 中, 灰数符号用⨂表示——笔者注) , 对这些给定灰数, 用已有的灰数运算规则对其进行计算[1 ,2 ,3 ,4 ,5 ] , 其计算结果为
x ( g ) = [ a 2 2 ( g ) − a 2 1 ( g ) ] / [ ( a 1 1 ( g ) + a 2 2 ( g ) ) − ( a 1 2 ( g ) + a 2 1 ( g ) ) ] = [ 1 / 7 , 3 / 5 ] ( 1 ) max { x ( g ) } = 3 7 ∣ ∣ a 2 1 = 3 , a 2 2 = 0 , min { x ( g ) } = 1 5 ∣ ∣ a 2 1 = 2 , a 2 2 = 1 ( 2 ) x ( g ) = [ a 2 2 ( g ) - a 2 1 ( g ) ] / [ ( a 1 1 ( g ) + a 2 2 ( g ) ) - ( a 1 2 ( g ) + a 2 1 ( g ) ) ] = [ 1 / 7 , 3 / 5 ] ( 1 ) max { x ( g ) } = 3 7 | a 2 1 = 3 , a 2 2 = 0 , min { x ( g ) } = 1 5 | a 2 1 = 2 , a 2 2 = 1 ( 2 )
然而, 用普通的数学计算方法所算得式 (1) [8 ] 的最大和最小值如式 (2) 所示。也就是说, 如果按普通的经典数学计算方法, 式 (1) 的可能值域 (区间灰数) 应为x ′ (g ) =[1/5, 3/7]。
显然, 由式 (1) 和 (2) 的运算结果可知, 1/7<[1/5, 3/7]<3/5, 它表明:采用目前的区间灰数运算规则对区间灰数进行运算, 会造成运算结果与经典数学的运算结果的不一不致。
根据灰数的灰度定义[1 ] , 可计算出x (g ) 和x ′ (g ) 的灰度分别为:g °1 (g ) =1.2308和g °2 (g ) =0.7273。
由此可见, 采用目前的区间灰数运算规则对区间灰数进行运算, 其运算结果的灰度大于经典数学的运算结果, 主要原因在于:采用传统的灰数计算方法可能会造成同一个灰数在相同的计算条件下取了不同的值。例如, 在式 (1) 中, 当分子中取a 22 (g ) =1, a 21 (g ) =2;而分母中取a 22 (g ) =0, a 21 (g ) =3时, 可得式 (1) 中x (g ) 的左端点的值为1/7;这使得x (g ) 的左端点值被不正常地缩小。同样的理由, 传统的灰数计算方法, 在某种情况下也会使其右端点的值被不正常地放大。
为了克服这种灰代数计算方法的蔽端, 笔者仅对这一问题进行一些讨论。
《2 标准区间灰数及其运算》
2 标准区间灰数及其运算
定义1 (标准区间灰数) 若某区间灰数可表示成式 (3) 所示的形式, 且在式 (3) 中, a i 称为G i 的白部, 而
G i = g i + c i γ i , i = 1 , 2 , ⋯ ( 3 ) G i = g i + c i γ i , i = 1 , 2 , ⋯ ( 3 )
c i γ i 称为G i 的灰部, 其中c i 称为灰系数, γ i 称为单位灰数 (或灰数单位) , 则称式 (3) 所表示的灰数形式为标准区间灰数[6 ,7 ,8 ,9 ,10 ] , 或简称为标准灰数。
定义1规定了区间灰数的标准表达形式。根据定义1, 若存在2个标准区间灰数G K =a K +c K γ K 和G L =-a L -c L γ L , 且G K +G L =0, 则称G K 和G L 互为相反灰数。
若存在2个标准区间灰数G K =a K +c K gγ K 和G L =1/G K , 则称G K 和G L 互为灰倒数。
定理1 (灰数的标准表达式) 任一区间灰数G i ∈[a i , b i ], a i ≤b i , i =1, 2, … 均可表示成定义1中的标准灰数的形式[2 ,3 ] , 见式 (3) 。
证明 为不失一般性, 对任一区间灰数G i ∈[a i , b i ], a i ≤b i , i =1, 2, … 进行如下变换:
G i ∈ [ a i , b i ] = [ a i , b i ] + [ a i , a i ] − [ a i , a i ] = [ a i , a i ] + [ 0 , b i − a i ] = a i + ( b i − a i ) ⋅ [ 0 , 1 ] ‚ G i ∈ [ a i , b i ] = [ a i , b i ] + [ a i , a i ] - [ a i , a i ] = [ a i , a i ] + [ 0 , b i - a i ] = a i + ( b i - a i ) ⋅ [ 0 , 1 ] ‚
令c i = (b i -a i ) , γ i ∈[0, 1];则由上式可得
G i = a i + c i γ i ( 4 ) G i = a i + c i γ i ( 4 )
为标准区间灰数的表达形式, 式中a i 称为G i 的白部, 而c i γ i 称为G i 的灰部, 在灰部c i γ i 中c i =1/ (b i -a i ) 称为灰系数, γ i 称为单位灰数 (或灰数单位) 。
定理1构造性地证明, 阐述任一区间灰数都可通过标准变换将其变为式 (4) 所示的标准形式。
定义2 (标准区间灰数的运算规则) 若任意若干个标准区间灰数之间的代数运算F (g ) =g (G 1 , G 2 , …, G n ) 仍采用经典数学的运算规则[2 ,3 ] , 且将这若干个灰数的取数γ i , γ i ∈[0, 1], i =1, 2, … 定为某确定的常数时, 所得到的最小值和最大值:
min { F ( g ) } = min g ( G 1 , G 2 , ⋯ , G n ) | γ i = c i , c i ∈ [ 0 , 1 ] , i = 1 , 2 , ⋯ n ‚ max { F ( g ) } = max g ( G 1 , G 2 , ⋯ , G n ) | γ i = c i , c i ∈ [ 0 , 1 ] , i = 1 , 2 , ⋯ n ‚ min { F ( g ) } = min g ( G 1 , G 2 , ⋯ , G n ) | γ i = c i , c i ∈ [ 0 , 1 ] , i = 1 , 2 , ⋯ n ‚ max { F ( g ) } = max g ( G 1 , G 2 , ⋯ , G n ) | γ i = c i , c i ∈ [ 0 , 1 ] , i = 1 , 2 , ⋯ n ‚
分别作为F (g ) =g (G 1 , G 2 , …, G n ) 的左右端点值的运算过程称为标准区间灰数的运算过程。
《3 第一和第二标准灰数》
3 第一和第二标准灰数
任给某一灰数G i ∈[a i , b i ], (a i ≤b i , i =1, 2, …) , 若对其左端点值进行标准化变换, 即
G ( 1 ) i ∈ [ a i , b i ] = a i − a i + [ a i , b i ] = a i + [ 0 , b i − a i ] = a i + ( b i − a i ) [ 0 , 1 ] = a i + ( b i − a i ) γ ( 1 ) i , ( 0 ≤ γ ( 1 ) i ≤ 1 ) ( 5 ) G ( 2 ) i ∈ [ a i , b i ] = b i − b i + [ a i , b i ] = b i − [ 0 , b i − a i ] = b i − ( b i − a i ) [ 0 , 1 ] = b i − ( b i − a i ) γ ( 2 ) i , ( 0 ≤ γ ( 2 ) i ≤ 1 ) ( 6 ) G i ( 1 ) ∈ [ a i , b i ] = a i - a i + [ a i , b i ] = a i + [ 0 , b i - a i ] = a i + ( b i - a i ) [ 0 , 1 ] = a i + ( b i - a i ) γ i ( 1 ) , ( 0 ≤ γ i ( 1 ) ≤ 1 ) ( 5 ) G i ( 2 ) ∈ [ a i , b i ] = b i - b i + [ a i , b i ] = b i - [ 0 , b i - a i ] = b i - ( b i - a i ) [ 0 , 1 ] = b i - ( b i - a i ) γ i ( 2 ) , ( 0 ≤ γ i ( 2 ) ≤ 1 ) ( 6 )
任一标准灰数, 均可标准化为第一和第二标准化灰数的形式。尽管这两种不同表达形式的标准灰数在形式上是有差异的, 但本质上应该是一致的, 因为他们所表达的是同一个灰数。
定理2 (第一和第二标准灰数的取数之和为1) 任给某一灰数G i =[a i , b i ], (a i ≤b i , i =1, 2, … ) , 若将其表达成第一标准灰数 (G ( 1 ) i i ( 1 ) ) 和第二标准灰数 (G (2) i ) , 其中 γ ( 1 ) i i ( 1 ) , (0≤γ ( 1 ) i i ( 1 ) ≤1) 和γ ( 2 ) i i ( 2 ) , (0≤γ ( 2 ) i i ( 2 ) ≤1) 分别表示这两标准灰数的取数, 则对于同一个灰数而言, 这两种不同表达形式的标准灰数取数之和为1, 即γ (1) i + γ ( 2 ) i i ( 2 ) =1。
证明 为不失一般性, 任给某一灰数G i ∈[a i , b i ], (a i ≤b i , i =1, 2, … ) , 按定义1, 将其表示成第一和第二标准灰数的形式, 如式 (6) 和式 (7) 所示。
由于式 (6) 和式 (7) 所表达的是同一个灰数, 因此, 尽管这两种不同表达形式的标准灰数在形式上是有差异的, 但本质上是应该相等的, 即
G ( 1 ) i = G ( 2 ) i , a i + ( b i − a i ) γ ( 1 ) i = b i − ( b i − a i ) γ ( 2 ) i , γ ( 1 ) i + γ ( 2 ) i = 1 ( 7 ) G i ( 1 ) = G i ( 2 ) , a i + ( b i - a i ) γ i ( 1 ) = b i - ( b i - a i ) γ i ( 2 ) , γ i ( 1 ) + γ i ( 2 ) = 1 ( 7 )
由定理2可知, 同一灰数的不同标准的表达形式不同, 但并没有本质上的差别[2 ,3 ,4 ] , 可以用图1来表示。尽管这两种标准灰数的表达形式没有差别, 但是, 在实际问题的计算过程中, 为了便于对问题的理解和避免取数的混淆, 一般情况下, 在同一问题的计算过程中, 应尽量采用同一种形式的标准灰数, 尤其是对于同一灰数而言, 更应如此。
《图1》
图1 第一和二标准灰数取数γ (1) i 和γ (2) i 之间的关系
Fig.1 Relationship between γ (1) i and γ (2) i
《4 标准灰数的大小判定》
4 标准灰数的大小判定
灰数的大小比较是灰系统理论与应用研究中涉及到的最基本也是最重要的问题之一。首先讨论单个标准灰数的大小判定问题, 然后再研究由单个标准灰数所复合的复合标准灰数的大小判定问题。
任意2个灰数G 1 ∈[a 1 , b 1 ]=a 1 +b 1 γ 1 和G 2 =[a 2 , b 2 ]=a 2 +b 2 γ 2 , 若满足:
若G 1 -G 2 ≥0, 则G 1 ≥G 2 ;若G 1 -G 2 <0, 则G 1 <G 2 。
根据定义4, 可以将2个单个标准灰数G 1 ∈[a , b ], (a ≤b ) 和G 2 ∈[c , d ] (c ≤d ) 之间的大小判定问题归结为以下情况:
1) 当b <c 时, 称[a , b ]<[c , d ], 即灰数[a , b ] (a ≤b ) 小于灰数 [c , d ] (c ≤d ) ;
2) 当b ≤c 时, 称 [a , b ]≤[c , d ], 即灰数 [a , b ] (a ≤b ) 小于等于灰数 [c , d ] (c ≤d ) ; 且仅当b =c 时, 等号成立;
3) 当a =c , b =d , 且两灰数的取数一致 (γ 1 =γ 2 ) 时[1 ] , 即满足
G 1 ∈ [ a 1 , b 1 ] = γ 1 a 1 + ( 1 − γ 1 ) b 1 , γ 1 = [ 0 , 1 ] , G 2 ∈ [ a 2 , b 2 ] = a 2 + b 2 γ 2 , γ 2 = [ 0 , 1 ] , 且 γ 1 = γ 2 ‚ 则 称 [ a , b ] = [ c , d ] , 即 灰 数 [ a , b ] ( a ≤ b ) G 1 ∈ [ a 1 , b 1 ] = γ 1 a 1 + ( 1 - γ 1 ) b 1 , γ 1 = [ 0 , 1 ] , G 2 ∈ [ a 2 , b 2 ] = a 2 + b 2 γ 2 , γ 2 = [ 0 , 1 ] , 且 γ 1 = γ 2 ‚ 则 称 [ a , b ] = [ c , d ] , 即 灰 数 [ a , b ] ( a ≤ b )
4) 当a = c , b = d , 且两灰数的取数不一致 (γ 1 ≠ γ 2 ) 时, 若 (γ 1 < γ 2 ) , 则G 1 > G 2 ;若γ 1 > γ 2 , 则G 1 <G 2 。
以上给出了在一般情况下较易于对其大小进行判别的几种最简单的灰数形式及其判定方法。
定义5 (复合标准灰数) 若某灰数 (G K ) 是由若干个标准灰数 (G ij ) 经过若干次加、减、乘、除等其他形式的数学运算而得到的, 则称该灰数 (G K ) 为复合标准灰数, 简称为复合灰数, 表示为
G K = f ( g i j ) γ i j = [ min γ i j f ( ⋅ ) , max γ i j f ( ⋅ ) ] , ( 0 ≤ γ i j ≤ 1 , i , j = 1 , 2 , ⋯ ) 。 ( 8 ) G Κ = f ( g i j ) γ i j = [ min γ i j f ( ⋅ ) , max γ i j f ( ⋅ ) ] , ( 0 ≤ γ i j ≤ 1 , i , j = 1 , 2 , ⋯ ) 。 ( 8 )
由定义5可知, 由单个标准灰数经过若干次数学运算而得到的灰数为复合标准灰数;由复合标准灰数, 或者复合标准灰数与单个标准灰数经过若干次数学运算而得到的灰数也称为复合标准灰数。复合标准灰数与第一、第二标准灰数也具有相同的表达形式, 具有白部、灰部和灰系数, 只是其灰部的灰系数系由若干单个灰数的灰系数复合而成。
定义6 (复合灰数大小的判定) 任给2个复合灰数:G K =f (g ij ) γ ij , (0≤γ ij ≤1, i , j =1, 2, … ) 和G L =f (g uv ) γ uv (0≤γ uv ≤1, u , v =1, 2, …) , 若无论γ ij 和γ uv 取任何取数, 都能满足G K ≥G L , 即G K |γ ij =[0, 1] ≥G L |γ uvj =[0, 1] , 则称复合灰数G K 不比G L 小; 若无论γ ij 和γ uv 取任何取数, 都能满足G K <G L , 即G K |γ ij =[0, 1] <G L |γ uv =[0, 1] , 则称复合灰数G K 比G L 小; 若无论γ ij 和γ uv 取任何取数, 都能满足G K =G L , 即G K |γ ij =[0, 1] =G L |γ uv =[0, 1] , 则称复合灰数G K 等于G L ; 否则, 则称复合灰数G K 和G L 之间的大小无法判定[1 ,2 ,3 ,4 ] 。
例2 给定2个单个标准灰数, G 1 ∈a + (b -a ) γ 1 , (0≤γ 1 ≤1) ; G 1 =c + (d -c ) γ 2 , (0≤γ 2 ≤1) ;则, 这两个标准灰数通过加法和减法而复合的复合标准灰数分别为:G ( 1 ) 1 2 1 2 ( 1 ) , G ( 2 ) 1 2 1 2 ( 2 ) 。
G ( 1 ) 1 2 = G 1 + G 2 = a + ( b − a ) γ 1 + c + ( d − c ) γ 2 = ( a + c ) + ( ( b − a ) γ 1 + ( d − c ) γ 2 ) ( 9 ) G ( 2 ) 1 2 = G 1 − G 2 = a + ( b − a ) γ 1 − c − ( d − c ) γ 2 = ( a − c ) + ( ( b − a ) γ 1 − ( d − c ) γ 2 ) ( 1 0 ) G 1 2 ( 1 ) = G 1 + G 2 = a + ( b - a ) γ 1 + c + ( d - c ) γ 2 = ( a + c ) + ( ( b - a ) γ 1 + ( d - c ) γ 2 ) ( 9 ) G 1 2 ( 2 ) = G 1 - G 2 = a + ( b - a ) γ 1 - c - ( d - c ) γ 2 = ( a - c ) + ( ( b - a ) γ 1 - ( d - c ) γ 2 ) ( 1 0 )
《5 案例研究》
5 案例研究
根据标准灰数的表征与运算规则, 以基于灰矩阵形式的2×2的零和矩阵博弈为例, 计算其混合策略的解, 如例3所示。
例3 求解灰矩阵博弈G (g ) ={S 1 (g ) , S 2 (g ) , A (g ) }, 其中:
A ( g ) = ( [ 0 ‚ 1 ] 4 [ 2 ‚ 3 ] 2 ) = ( γ 1 1 4 2 + γ 1 2 2 ) ‚ ( γ 1 1 ∈ [ 0 ‚ 1 ] ) ‚ ( γ 1 2 ∈ [ 0 ‚ 1 ] ) A ( g ) = ( [ 0 ‚ 1 ] [ 2 ‚ 3 ] 4 2 ) = ( γ 1 1 2 + γ 1 2 4 2 ) ‚ ( γ 1 1 ∈ [ 0 ‚ 1 ] ) ‚ ( γ 1 2 ∈ [ 0 ‚ 1 ] )
解:运用2×2的零和矩阵博弈混合策略的求解公式[4 ] 与本文提供的标准灰数的表征与计算方法, 并由题意知:a 11 (g ) =γ 11 , a 12 (g ) =2+γ 12 , a 21 (g ) =4, a 22 (g ) =2, 故局中人1和2的最优灰混合策略及其最优灰博弈值分别如式 (11) 至式 (15) 所示。原灰数计算方法与标准灰数计算方法的比较情况, 如表1所示。
x ∗ 1 ( g ) = a 2 2 ( g ) − a 2 1 ( g ) ( a 1 1 ( g ) − a 2 2 ( g ) ) − ( a 1 2 ( g ) − a 2 1 ( g ) ) = { 2 3 , i f , γ 1 1 = 1 , γ 1 2 = 0 2 5 , i f , γ 1 1 = 0 , γ 1 2 = 1 , 即 ∶ x ∗ 1 ( g ) = [ 2 5 , 2 3 ] ( 1 1 ) x ∗ 2 ( g ) = a 1 1 ( g ) − a 1 2 ( g ) ( a 1 1 ( g ) − a 2 2 ( g ) ) − ( a 1 2 ( g ) − a 2 1 ( g ) ) = { 1 3 , i f , γ 1 1 = 1 , γ 1 2 = 0 3 5 , i f , γ 1 1 = 0 , γ 1 2 = 1 , 即 ∶ x ∗ 2 ( g ) = [ 1 3 , 3 5 ] ( 1 2 ) y ∗ 1 ( g ) = a 2 2 ( g ) − a 1 2 ( g ) ( a 1 1 ( g ) − a 2 2 ( g ) ) − ( a 1 2 ( g ) − a 2 1 ( g ) ) = { 0 , i f , γ 1 2 = 0 1 4 , i f , γ 1 1 = 1 , γ 1 2 = 1 ′ , 即 ∶ y ∗ 1 ( g ) = [ 0 , 1 4 ] ( 1 3 ) y ∗ 2 ( g ) = a 1 1 ( g ) − a 2 1 ( g ) ( a 1 1 ( g ) − a 2 2 ( g ) ) − ( a 1 2 ( g ) − a 2 1 ( g ) ) = { 1 , i f , γ 1 2 = 0 4 5 , i f , γ 1 1 = 0 , γ 1 2 = 1 ′ , 即 ∶ y ∗ 2 ( g ) = [ 4 5 , 1 ] ( 1 4 ) v ∗ G ( g ) = a 1 1 ( g ) ⋅ a 2 2 ( g ) − a 1 2 ( g ) ⋅ a 2 1 ( g ) ( a 1 1 ( g ) − a 2 2 ( g ) ) − ( a 1 2 ( g ) − a 2 1 ( g ) ) = { 2 , i f , γ 1 2 = 0 5 2 , i f , γ 1 1 = 1 , γ 1 2 = 1 ′ , 即 ∶ v ∗ G ( g ) = [ 2 , 5 2 ] ( 1 5 ) x 1 * ( g ) = a 2 2 ( g ) - a 2 1 ( g ) ( a 1 1 ( g ) - a 2 2 ( g ) ) - ( a 1 2 ( g ) - a 2 1 ( g ) ) = { 2 3 , i f , γ 1 1 = 1 , γ 1 2 = 0 2 5 , i f , γ 1 1 = 0 , γ 1 2 = 1 , 即 ∶ x 1 * ( g ) = [ 2 5 , 2 3 ] ( 1 1 ) x 2 * ( g ) = a 1 1 ( g ) - a 1 2 ( g ) ( a 1 1 ( g ) - a 2 2 ( g ) ) - ( a 1 2 ( g ) - a 2 1 ( g ) ) = { 1 3 , i f , γ 1 1 = 1 , γ 1 2 = 0 3 5 , i f , γ 1 1 = 0 , γ 1 2 = 1 , 即 ∶ x 2 * ( g ) = [ 1 3 , 3 5 ] ( 1 2 ) y 1 * ( g ) = a 2 2 ( g ) - a 1 2 ( g ) ( a 1 1 ( g ) - a 2 2 ( g ) ) - ( a 1 2 ( g ) - a 2 1 ( g ) ) = { 0 , i f , γ 1 2 = 0 1 4 , i f , γ 1 1 = 1 , γ 1 2 = 1 ′ , 即 ∶ y 1 * ( g ) = [ 0 , 1 4 ] ( 1 3 ) y 2 * ( g ) = a 1 1 ( g ) - a 2 1 ( g ) ( a 1 1 ( g ) - a 2 2 ( g ) ) - ( a 1 2 ( g ) - a 2 1 ( g ) ) = { 1 , i f , γ 1 2 = 0 4 5 , i f , γ 1 1 = 0 , γ 1 2 = 1 ′ , 即 ∶ y 2 * ( g ) = [ 4 5 , 1 ] ( 1 4 ) v G * ( g ) = a 1 1 ( g ) ⋅ a 2 2 ( g ) - a 1 2 ( g ) ⋅ a 2 1 ( g ) ( a 1 1 ( g ) - a 2 2 ( g ) ) - ( a 1 2 ( g ) - a 2 1 ( g ) ) = { 2 , i f , γ 1 2 = 0 5 2 , i f , γ 1 1 = 1 , γ 1 2 = 1 ′ , 即 ∶ v G * ( g ) = [ 2 , 5 2 ] ( 1 5 )
表1 两种不同计算方法的结果比较
Table 1 Compare data of two arithmetic
《表1》
指标 原灰数计算方法标准灰数计算方法 计算结果灰度 计算结果 灰度 x * 1 (g ) [5/2, 2/3] 0.5000 [5/2, 2/3] 0.5000 x * 2 (g ) [1/5, 1] 1.330 [1/3, 3/5] 0.5700 y * 1 (g ) [0, 1/3] 2.0000 [0, 1/4] 2.0000 y * 2 (g ) [3/5, 4/3] 0.4889 [5/4, 1] 0.2222 v * G (g ) [6/5, 4] 1.0769 [2, 5/2] 0.4444
《6 结语》
6 结语
笔者定义了标准区间灰数与第一和第二标准区间灰数的概念, 设计了普通区间灰数与标准区间灰数之间的转换规则, 提供了标准区间灰数之间的比较与运算法则, 从而较好地解决了区间灰数之间的大小比较与运算问题, 并将这一研究成果应用于基于区间数的GM (1, 1) 模型预测问题, 取得了良好的效果。