《1 引言》
1 引言
近年来, 国内外许多学者致力于多目标决策 (MADM, multiple attribute decision making) 问题的研究, 如文献[1 ,2 ] 就较全面地介绍了MADM的理论和方法, 目前MADM问题的研究主要集中在解决具有多个指标、有限个方案的静态决策问题, 关于动态决策和评估问题的研究近几年相继有一些文献[3 ,4 ,5 ] , 但在已有的文献中, 大多是将较成熟的静态多指标决策理论和方法移植到动态决策问题上。由于在工程技术系统、社会、经济、管理中存在众多多时段的综合决策与评价问题, 这类问题实际上是在决策空间和目标空间的基础上, 又增加了一个时间空间, 是具有时间、指标、方案的三维决策排序问题。在对这一类问题的决策时, 由于指标集中, 有些指标具有不同的量纲, 对它们难以进行直接比较, 因而需要对原始决策矩阵进行初始化处理。由于初始化处理方法较多, 使用不同的初始化处理方法可能会得出不同的评价结果。而现有的对决策矩阵处理都是采用在[0, 1]区间上的线性变换法, 这种方法存在只奖不罚的不足[4 ] 。由于在实际中有时某事物的积极因素与消极因素之间可能相互抵消, 这种现象在实际生产和社会生活中是大量存在的, 因此有必要在对决策矩阵进行变换时, 把规范化决策矩阵中的数值从[0, 1]扩充到[-1, 1]上。笔者利用Vague思想[6 ,7 ,8 ] 和集对分析理论的思想[9 ,10 ] , 提出了一种易于计算且实用的[-1, 1]线性变换算子, 该变换把规范化决策矩阵的数据从[0, 1]区间上拓展到[-1, 1]上, 利用此方法得到规范化决策矩阵, 利用规范化决策矩阵寻求各时段的正、负理想方案, 并在此基础上建立了一种新的基于动态多指标灰色关联分析决策模型。
《2 原理与方法》
2 原理与方法
设动态多指标决策问题有n 个被评估对象或拟定的决策方案组成决策方案集S ;m 个评价指标或属性组成指标集A , A = {A 1 , A 2 , …, A m };方案S i 在t 时段 (阶段) (t = 1, 2, …, T ) 对指标A j 的属性值为x ij (t ) (i =1, 2, …, n , j = 1, 2, …, m ) , 则t 时段方案集S 对指标集A 的决策矩阵为
X ( t ) = ⎡ ⎣ ⎢ ⎢ ⎢ ⎢ ⎢ x 1 1 ( t ) x 2 1 ( t ) ⋮ x n 1 ( t ) x 1 2 ( t ) x 2 2 ( t ) ⋮ x n 2 ( t ) ⋯ ⋯ ⋯ x 1 m ( t ) x 2 m ( t ) ⋮ x n m ( t ) ⎤ ⎦ ⎥ ⎥ ⎥ ⎥ ⎥ 。 X ( t ) = [ x 1 1 ( t ) x 1 2 ( t ) ⋯ x 1 m ( t ) x 2 1 ( t ) x 2 2 ( t ) ⋯ x 2 m ( t ) ⋮ ⋮ ⋮ x n 1 ( t ) x n 2 ( t ) ⋯ x n m ( t ) ] 。
指标属性集A ={A 1 , A 2 , …, A m }一般情况下可分为3种类型, 即效益型、成本型和区间型指标 (固定型可视为区间型的特例) 。效益型指标是其值越大越好;成本型指标是其值越小越好;区间型指标是其值落在某一特定区间[A , B ]为最好。
令, Ω = ∪ i = 1 3 Ω i Ω = ∪ i = 1 3 Ω i , 且Ω i ∩Ω j =Ø, i , j = 1, 2, 3, i ≠j , 其中Ω 1 为效益型指标集, Ω 2 为成本型指标集, Ω 3 为区间型指标集。
《2.1 决策矩阵规范化处理方法》
2.1 决策矩阵规范化处理方法
由于指标集中的指标具有不同的量纲, 在决策时, 它们难以进行直接比较, 因而需要对原始决策矩阵进行初始化处理。初始化处理方法较多[11 ] , 如均值化变换、始点零象化变换, 初值化变换、百分比变换、归一化变换、级差最大值化变换、区间值化变换等等, 使用不同的初始化变换处理方法可能会导致评价结果的不同, 而现有的对决策矩阵处理都是采用的[0, 1]区间线性变换法, 这种方法存在只奖不罚的不足。因此, 笔者利用Vague思想和集对分析的理论思想, 把[0, 1]区间上的线性变换拓展到[-1, 1]上的线性变换, 提出了一种易于计算且实用的[-1, 1]线性变换算子, 而且也弥补了文献[5 ] 中提到的不足, 其基本思想是对于评价对象的指标值如果优于平均水平时, 赋予0 ~ 1的正值, 如果劣于平均水平时, 赋予0 ~-1的负值, 具体做法为
首先令Z j ( t ) = 1 n ∑ i = 1 n x i j ( t ) ‚ j = 1 , 2 , ⋯ , m , t = 1 , 2 , ⋯ , T Ζ j ( t ) = 1 n ∑ i = 1 n x i j ( t ) ‚ j = 1 , 2 , ⋯ , m , t = 1 , 2 , ⋯ , Τ 。
r i j ( t ) = [ x i j ( t ) − z j ( t ) ] / [ max ( max j { x i j ( t ) } − z j ( t ) , z j ( t ) − min j { x i j ( t ) } ) ] ; ( 1 ) r i j ( t ) = [ x i j ( t ) - z j ( t ) ] / [ max ( max j { x i j ( t ) } - z j ( t ) , z j ( t ) - min j { x i j ( t ) } ) ] ; ( 1 )
r i j ( t ) = [ z j ( t ) − x i j ( t ) ] / [ max ( max j { x i j ( t ) } − z j ( t ) , z j ( t ) − min j { x i j ( t ) } ) ] ; ( 2 ) r i j ( t ) = [ z j ( t ) - x i j ( t ) ] / [ max ( max j { x i j ( t ) } - z j ( t ) , z j ( t ) - min j { x i j ( t ) } ) ] ; ( 2 )
若A j 为[A , B ]区间型指标 (包括固定型指标, 此时A =B ) , 则
r i j ( t ) = ⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ 1 − 2 ( A − x i j ( t ) ) A − min j ( x i j ( t ) } x i j ( t ) < A , 1 − 2 ( x i j ( t ) − B ) max j { x i j ( t ) } − B x i j ( t ) > B ‚ 1 x i j ( t ) ∈ [ A , B ] 。 ( 3 ) r i j ( t ) = { 1 - 2 ( A - x i j ( t ) ) A - min j ( x i j ( t ) } x i j ( t ) < A , 1 - 2 ( x i j ( t ) - B ) max j { x i j ( t ) } - B x i j ( t ) > B ‚ 1 x i j ( t ) ∈ [ A , B ] 。 ( 3 )
通过此线性变换算子对决策矩阵X (t ) = (x ij ) 进行规范化变换, 则规范化决策矩阵R (t ) = (γ ij ) , R (t ) = (γ ij ) 中的元素都是无量纲的, 并且所有的元素均符合奖优罚劣的标准及Vague思想, 而对任意的γ ij (t ) ∈[-1, 1], (i = 1, 2, …, n , j =1, 2, …, m , t = 1, 2, …, T ) , 由于规范化决策矩阵中的元素无量纲, 因此它们之间可以进行直接比较;并且由规范化决策矩阵容易求得理想方案和负理想方案。理想方案中各元素都是正值, 负理想方案中各元素都是负值, 这种方法所得到的理想方案和负理想方案与人们在实际中的思维印象是一致的。
《2.2 灰色关联度分析》
2.2 灰色关联度分析
定义 设r + j (t ) =max{r ij (t ) 1≤i ≤n }
r − j ( t ) = min { r i j ( t ) 1 ≤ i ≤ n } ( j = 1 , 2 , ⋯ , m ) , r j - ( t ) = min { r i j ( t ) 1 ≤ i ≤ n } ( j = 1 , 2 , ⋯ , m ) ,
S + ( t ) = { r + 1 ( t ) , r + 2 ( t ) , ⋯ , r + m ( t ) } ( 4 ) S + ( t ) = { r 1 + ( t ) , r 2 + ( t ) , ⋯ , r m + ( t ) } ( 4 )
S − ( t ) = { r − 1 ( t ) , r − 2 ( t ) , ⋯ , r − m ( t ) } ( 5 ) S - ( t ) = { r 1 - ( t ) , r 2 - ( t ) , ⋯ , r m - ( t ) } ( 5 )
分别为t 时段 (t =1, 2, …, T ) 决策的理想方案和负理想方案。
显然, 理想方案与负理想方案都不存在, 若理想方案存在, 就直接选择理想方案作为选择对象, 无需挑选。若负理想方案存在, 由于它是所有方案中是最坏的, 因而可以直接从方案集中删除。
设通过Delphi调查法或AHP法确定了每个指标的权重向量为
w = ( w 1 , w 2 , ⋯ , w m ) , w j > 0 ( j = 1 , 2 , ⋯ , m ) , ∑ j = 1 m w j = 1 。 w = ( w 1 , w 2 , ⋯ , w m ) , w j > 0 ( j = 1 , 2 , ⋯ , m ) , ∑ j = 1 m w j = 1 。
由灰色关联分析方法可知[12 ] , 在t 时段第i 方案与理想方案S + (t ) 关于指标A j (t ) 的关联系数为
ξ + i j ( t ) = [ min i min j | r i j ( t ) − r + j ( t ) | + ρ max i max j | r i j ( t ) − r + j ( t ) | ] / [ | r i j ( t ) − r + j ( t ) | + ρ max i max j | r i j ( t ) − r + j ( t ) | ] , ( i = 1 , 2 , ⋯ , n , j = 1 , 2 , ⋯ , m ) ( 6 ) ξ i j + ( t ) = [ min i min j | r i j ( t ) - r j + ( t ) | + ρ max i max j | r i j ( t ) - r j + ( t ) | ] / [ | r i j ( t ) - r j + ( t ) | + ρ max i max j | r i j ( t ) - r j + ( t ) | ] , ( i = 1 , 2 , ⋯ , n , j = 1 , 2 , ⋯ , m ) ( 6 )
在t 时段第i 方案与负理想方案S - (t ) 关于指标A j (t ) 的关联系数为
ξ − i j ( t ) = [ min i min j | r i j ( t ) − r − j ( t ) | + ρ max i max j | r i j ( t ) − r − j ( t ) | ] / [ | r i j ( t ) − r − j ( t ) | + ρ max i max j | r i j ( t ) − r − j ( t ) | ] ‚ ( i = 1 , 2 , ⋯ , n , j = 1 , 2 , ⋯ , m ) ( 7 ) ξ i j - ( t ) = [ min i min j | r i j ( t ) - r j - ( t ) | + ρ max i max j | r i j ( t ) - r j - ( t ) | ] / [ | r i j ( t ) - r j - ( t ) | + ρ max i max j | r i j ( t ) - r j - ( t ) | ] ‚ ( i = 1 , 2 , ⋯ , n , j = 1 , 2 , ⋯ , m ) ( 7 )
式中, ρ 为分辨系数, ρ ∈[0, 1], 一般取ρ = 0.5。
在t 时段, 第i 方案S i 与理想方案S + (t ) 和负理想方案S - (t ) 的关联度分别为
γ + i ( t ) = ∑ j = 1 m w j ξ + i j ( t ) , ( i = 1 , 2 , ⋯ , n , t = 1 , 2 , ⋯ , T ) ( 8 ) γ − i ( t ) = ∑ j = 1 m w j ξ − i j ( t ) , ( i = 1 , 2 , ⋯ , n , t = 1 , 2 , ⋯ , T ) ( 9 ) γ i + ( t ) = ∑ j = 1 m w j ξ i j + ( t ) , ( i = 1 , 2 , ⋯ , n , t = 1 , 2 , ⋯ , Τ ) ( 8 ) γ i - ( t ) = ∑ j = 1 m w j ξ i j - ( t ) , ( i = 1 , 2 , ⋯ , n , t = 1 , 2 , ⋯ , Τ ) ( 9 )
《2.3 建立动态多指标决策优化模型》
2.3 建立动态多指标决策优化模型
γ + i (t ) 越大, 决策方案S i 与理想方案S + (t ) 越接近, 方案则越佳;γ - i (t ) 的意义恰好相反。因此, 决策方案应该距理想方案最近, 而同时又距负理想方案最远。假设决策方案S i 以优属度u i 从属于理想方案S + (t ) , 那么决策方案S i 以1-u i 从属于负理想S - (t ) , 为确定优属度u i , 建立如下优化模型:
min F ( u ) = min ∑ t = 1 T { [ ( 1 − u i ) v ( t ) γ + i ( t ) ] 2 + [ u i v ( t ) γ − i ( t ) ] 2 } , ( 1 0 ) min F ( u ) = min ∑ t = 1 Τ { [ ( 1 - u i ) v ( t ) γ i + ( t ) ] 2 + [ u i v ( t ) γ i - ( t ) ] 2 } , ( 1 0 )
其中u = (u 1 , u 2 , …, u n ) 为系统的最优解向量;v 为决策方案各时段的权重向量, 其中v = (v (1) , v (2) , …, v (T ) ) , v (t ) >0 (t =1, 2, …, T ) , ∑ t = 1 T v ( t ) = 1 ∑ t = 1 Τ v ( t ) = 1 。
∂ F ∂ u i = 2 ∑ t = 1 T ( 1 − u i ) ( v ( t ) γ + i ( t ) ) ( − v ( t ) γ + i ( t ) ) + 2 ∑ t = 1 T u i ( v ( t ) γ − i ( t ) ) v ( t ) γ − i ( t ) = 0 , 故 u i = ( ∑ t = 1 T ( v ( t ) γ + i ( t ) ) 2 ) / ( ∑ t = 1 T ( v ( t ) γ + i ( t ) ) 2 + ∑ t = 1 T ( v ( t ) γ − i ( t ) ) 2 ) ( i = 1 , 2 , ⋯ , n ) 。 ( 1 1 ) ∂ F ∂ u i = 2 ∑ t = 1 Τ ( 1 - u i ) ( v ( t ) γ i + ( t ) ) ( - v ( t ) γ i + ( t ) ) + 2 ∑ t = 1 Τ u i ( v ( t ) γ i - ( t ) ) v ( t ) γ i - ( t ) = 0 , 故 u i = ( ∑ t = 1 Τ ( v ( t ) γ i + ( t ) ) 2 ) / ( ∑ t = 1 Τ ( v ( t ) γ i + ( t ) ) 2 + ∑ t = 1 Τ ( v ( t ) γ i - ( t ) ) 2 ) ( i = 1 , 2 , ⋯ , n ) 。 ( 1 1 )
由于u i 表示决策方案S i 优属于理想方案S + (t ) 的程度, u i 越大, 决策方案S i 越优, 反之, u i 越小, 决策方案S i 越差。
《3 应用实例》
3 应用实例
投资银行欲对某4家企业A 1 , A 2 , A 3 , A 4 进行投资, 选取5项指标进行决策, 分别为投资净产值率, 投资销售率, 投资成本率, 投资利税率, 环境污染程度。投资银行对这4家企业A 1 , A 2 , A 3 , A 4 进行了考察, 具体数据见文献[3 ] , 文献[3 ] 根据Delphi调查法, 按照专家们的意见确定了各指标的权重向量为w = (0.241, 0.225, 0.239, 0.217, 0.078) , 3个时段的权重向量为v = (0.329, 0.326, 0.345) 。
R ( 1 ) = ⎡ ⎣ ⎢ ⎢ ⎢ ⎢ − 0 . 2 3 8 4 0 . 8 6 0 9 − 1 0 . 3 7 7 4 − 0 . 4 3 3 7 0 . 8 4 3 9 − 1 0 . 5 8 9 8 0 . 0 5 7 1 − 0 . 6 8 5 7 − 0 . 3 7 1 4 1 0 . 5 4 3 8 − 0 . 0 8 7 7 0 . 5 4 3 8 − 1 − 0 . 9 3 1 − 0 . 6 5 5 1 1 0 . 5 8 6 2 ⎤ ⎦ ⎥ ⎥ ⎥ ⎥ ‚ R ( 2 ) = ⎡ ⎣ ⎢ ⎢ ⎢ ⎢ 1 − 0 . 1 1 1 1 − 0 . 3 3 3 3 − 0 . 5 5 5 6 1 − 0 . 5 5 1 0 − 0 . 2 2 4 4 − 0 . 2 2 4 4 − 0 . 4 5 4 5 − 0 . 2 5 − 0 . 2 9 5 4 1 0 . 6 4 1 7 1 − 0 . 9 7 0 1 − 0 . 6 7 1 6 − 1 0 . 5 7 5 7 0 . 8 1 8 1 − 0 . 3 9 3 9 ⎤ ⎦ ⎥ ⎥ ⎥ ⎥ ‚ R ( 3 ) = ⎡ ⎣ ⎢ ⎢ ⎢ ⎢ − 0 . 9 8 2 3 1 0 . 1 1 5 0 − 0 . 1 3 2 7 − 0 . 2 6 3 1 1 − 0 . 8 4 2 1 0 . 1 0 5 2 − 0 . 4 4 5 3 − 0 . 9 1 5 9 0 . 3 6 1 3 1 0 . 6 2 5 0 − 1 0 . 3 7 5 − 0 . 1 5 7 8 0 . 4 7 3 6 0 . 6 8 4 2 − 1 ⎤ ⎦ ⎥ ⎥ ⎥ ⎥ 。 R ( 1 ) = [ - 0 . 2 3 8 4 - 0 . 4 3 3 7 0 . 0 5 7 1 0 . 5 4 3 8 - 0 . 9 3 1 0 . 8 6 0 9 0 . 8 4 3 9 - 0 . 6 8 5 7 - 0 . 0 8 7 7 - 0 . 6 5 5 1 - 1 - 1 - 0 . 3 7 1 4 0 . 5 4 3 8 1 0 . 3 7 7 4 0 . 5 8 9 8 1 - 1 0 . 5 8 6 2 ] ‚ R ( 2 ) = [ 1 1 - 0 . 4 5 4 5 0 . 6 4 1 7 - 1 - 0 . 1 1 1 1 - 0 . 5 5 1 0 - 0 . 2 5 1 0 . 5 7 5 7 - 0 . 3 3 3 3 - 0 . 2 2 4 4 - 0 . 2 9 5 4 - 0 . 9 7 0 1 0 . 8 1 8 1 - 0 . 5 5 5 6 - 0 . 2 2 4 4 1 - 0 . 6 7 1 6 - 0 . 3 9 3 9 ] ‚ R ( 3 ) = [ - 0 . 9 8 2 3 - 0 . 2 6 3 1 - 0 . 4 4 5 3 0 . 6 2 5 - 0 . 1 5 7 8 1 1 - 0 . 9 1 5 9 0 0 . 4 7 3 6 0 . 1 1 5 0 - 0 . 8 4 2 1 0 . 3 6 1 3 - 1 0 . 6 8 4 2 - 0 . 1 3 2 7 0 . 1 0 5 2 1 0 . 3 7 5 - 1 ] 。
在第一时段, 理想方案S + (1) 和负理想方案S - (1) 分别为
S + ( 1 ) = { 0 . 8 6 0 9 , 0 . 8 4 3 9 , 1 , 0 . 5 4 3 8 , 1 } ‚ S − ( 1 ) = { − 1 , − 1 , − 0 . 6 8 5 7 , − 1 , − 0 . 9 3 1 } 。 S + ( 1 ) = { 0 . 8 6 0 9 , 0 . 8 4 3 9 , 1 , 0 . 5 4 3 8 , 1 } ‚ S - ( 1 ) = { - 1 , - 1 , - 0 . 6 8 5 7 , - 1 , - 0 . 9 3 1 } 。
因而在第一时段, 所有方案S i 与理想方案S + (1) 和负理想方案S - (1) 的灰色关联度分别为
γ + 1 ( 1 ) = 0 . 5 7 3 4 , γ + 2 ( 1 ) = 0 . 7 1 2 9 , γ + 3 ( 1 ) = 0 . 5 5 3 4 , γ + 4 ( 1 ) = 0 . 7 1 5 8 ; γ − 1 ( 1 ) = 0 . 5 7 3 1 , γ − 2 ( 1 ) = 0 . 5 7 0 9 , γ − 3 ( 1 ) = 0 . 7 5 5 8 , γ 4 ( 1 ) = 0 . 5 1 8 7 。 γ 1 + ( 1 ) = 0 . 5 7 3 4 , γ 2 + ( 1 ) = 0 . 7 1 2 9 , γ 3 + ( 1 ) = 0 . 5 5 3 4 , γ 4 + ( 1 ) = 0 . 7 1 5 8 ; γ 1 - ( 1 ) = 0 . 5 7 3 1 , γ 2 - ( 1 ) = 0 . 5 7 0 9 , γ 3 - ( 1 ) = 0 . 7 5 5 8 , γ 4 ( 1 ) = 0 . 5 1 8 7 。
同理, 可计算在第二、第三时段所有方案S i 与理想方案和负理想方案的灰色关联度分别为
γ + 1 ( 2 ) = 0 . 7 3 9 6 , γ + 2 ( 2 ) = 0 . 5 7 0 8 , γ + 3 ( 2 ) = 0 . 4 3 8 6 , γ + 4 ( 2 ) = 0 . 5 3 3 6 , γ − 1 ( 2 ) = 0 . 7 0 4 9 , γ − 2 ( 2 ) = 0 . 7 0 0 6 , γ − 3 ( 2 ) = 0 . 8 2 3 4 , γ − 4 ( 2 ) = 0 . 7 5 1 0 ; γ + 1 ( 3 ) = 0 . 5 3 5 6 , γ + 2 ( 3 ) = 0 . 7 4 4 9 , γ + 3 ( 3 ) = 0 . 5 1 1 5 , γ + 4 ( 3 ) = 0 . 6 7 1 9 , γ − 1 ( 3 ) = 0 . 6 6 9 5 , γ − 2 ( 3 ) = 0 . 5 3 7 4 , γ − 3 ( ( 3 ) = 0 . 6 8 9 7 , γ − 4 ( 3 ) = 0 . 4 9 5 2 。 γ 1 + ( 2 ) = 0 . 7 3 9 6 , γ 2 + ( 2 ) = 0 . 5 7 0 8 , γ 3 + ( 2 ) = 0 . 4 3 8 6 , γ 4 + ( 2 ) = 0 . 5 3 3 6 , γ 1 - ( 2 ) = 0 . 7 0 4 9 , γ 2 - ( 2 ) = 0 . 7 0 0 6 , γ 3 - ( 2 ) = 0 . 8 2 3 4 , γ 4 - ( 2 ) = 0 . 7 5 1 0 ; γ 1 + ( 3 ) = 0 . 5 3 5 6 , γ 2 + ( 3 ) = 0 . 7 4 4 9 , γ 3 + ( 3 ) = 0 . 5 1 1 5 , γ 4 + ( 3 ) = 0 . 6 7 1 9 , γ 1 - ( 3 ) = 0 . 6 6 9 5 , γ 2 - ( 3 ) = 0 . 5 3 7 4 , γ 3 - ( ( 3 ) = 0 . 6 8 9 7 , γ 4 - ( 3 ) = 0 . 4 9 5 2 。
由式 (11) 可得各方案关于3个时段的综合排序最优决策向量为:
u = ( 0 . 4 7 4 3 , 0 . 5 6 1 0 , 0 . 3 0 7 9 , 0 . 5 4 2 2 ) 。 u = ( 0 . 4 7 4 3 , 0 . 5 6 1 0 , 0 . 3 0 7 9 , 0 . 5 4 2 2 ) 。
由优化模型及排序原则可知, 这4家企业的排序为A 2 , A 4 , A 1 , A 3 。结果与文献[3 ] 的相同, 但分辨程度更高, 它清晰地反映了与平均水平的差距, 物理概念也更加清晰。
《4 结论》
4 结论
笔者通过研究动态多指标系统的决策特点, 利用Vague思想和集对分析的理论思想, 提出了一种易于计算且实用的[-1, 1]线性变换算子, 利用此方法寻求各时段的正、负理想方案, 建立了一种新的基于动态多指标灰色关联分析决策模型, 在模型中充分考虑了各指标在系统中的成长特性, 将此特性用于灰色关联分析, 决策结果较为客观可靠, 且易于在计算机上实现。此方法为解决动态多指标决策问题提供了一条新的途经。