1. 引言
多准则决策理论(MCDM)在经济和管理科学领域具有广泛应用。当决策问题涉及的指标数量和方案数量较多时,建立合理的层次结构(目标层、准则层、方案层)对于问题的分析至关重要。Saaty [1]于上世纪70年代提出了一类定性定量相结合的层次分析法(AHP),该方法通过建立判断矩阵将多个准则的比较转化为两两比较,再通过算术平均值、几何平均值、最大特征根等方法计算得到各个准则的权重。由于其思路明确且操作简单,目前已在经济、金融、教育、管理、能源等诸多领域具有广泛应用(见[2]-[6])。但随着人们关注的问题越来越复杂,其涉及的准则数量也日益增多,从而使得准则之间的相互比较耗时耗力。与人力成本增大相比,更为严重的问题则是:随着判断矩阵的阶数逐渐增大,其一致性很难通过检验,从而不得不反复修正最初的判断矩阵,并且依旧不能确保新修正后的判断矩阵能够通过一致性检验。
针对AHP比较次数多且难以通过一致性检验的问题,Rezaei [7]在2015年提出了一类名为最优最劣法(Best-Worst Method,以下简记为BWM)的多准则决策方法,其基本思想是:首先,在所有准则
中确定一个最优(重要)准则
和一个最劣(不重要)准则
;然后,将
与
依次比较并用1~9之间的整数赋值,将准则间的重要性比较值构成行向量
;接着,再将所有准则
与
依次比较,得到比较向量
;最后,通过求解以下非线性优化模型(1.1)得到准则
的权重
:
(1.1)
注意到模型(1.1)通常存在多解,Rezaei在2016年又提出了一种基于区间优先度的区间分析法,详见[8]。该方法针对准则
和
的权重区间
和
,按下面公式定义了一种区间优先度:
(1.2)
并由此构建一个优先度矩阵:
(1.3)
令
(1.4)
则优先度矩阵(1.3)可转化为下面的0-1矩阵(1.5):
(1.5)
然后,将0-1矩阵
中的第i行所有元素相加,用其和来表示第i个准则
的得分
,即:
(1.6)
并根据该得分给所有准则排序。仔细观察上述方法可知,要得到0-1矩阵(1.5),必须计算优先度矩阵(1.3)。从公式(1.2)进一步可知,要计算区间优先度,必须求出每个准则的权重的上界和下界。然而,在实际应用中,求解非线性优化模型(1.1)通常是比较困难的,区间优先度(1.2)式的计算也是比较复杂的。因此,自然提出以下问题:
(1) 能否不计算准则权重的区间优先度(1.2)和优先度矩阵(1.3),更便捷地得到0-1矩阵(1.5)?
(2) 是否可以“绕过”求解非线性优化模型(1.1),直接根据比较向量
和
估计出准则的权重?
针对上述问题,本文将利用“区间中心比较”和“尺度变换思想”给出一个肯定的回答。余文安排如下:第二节回顾非线性BWM(1.1)的一种变形及权重上下界的求解方法;第三节给出基于“权重区间中心比较”的0-1矩阵(1.5)的计算方法,并给出几种权重区间中心点的估计方法;第四节给出本文所提方法在电网关键物资选取中的实际应用;第五节总结本文的主要结果并探讨未来值得研究的一些新问题。
2. 预备知识
本节首先回顾非线性BWM(1.1)的区间分析法。为进一步简化非线性优化模型(1.1)中的目标函数,Rezaei在[8]中将该模型转化为如下等价的优化问题:
(2.1)
求解(2.1),可获得目标函数
的最小值
。然而,当
达到此最小值
时,其最优权重的取值一般并不唯一。正如文[8]所示,该模型“是否存在多解”与准则之间相互比较的“一致性”密切相关。下面给出“一致性”的定义以及模型(2.1)存在唯一解的条件,相关证明过程详见[8]。
定义1. 设有n个准则
,
表示最优准则
与第j个准则
的重要性比较值,
表示第j个准则
与最劣准则
的重要性比较值。如果对于任意
,都有
,则称该比较具有完全一致性。
命题1. 设有n个准则
。如果准则之间的比较具有完全一致性(即:对任意j,都有
),则模型(2.1)具有唯一解。
如果模型(2.1)存在多解,则可通过构建以下两个优化模型(以
为例)求得其最优权重的下界和上界:
(2.2)
和
(2.3)
令模型(2.2)和(2.3)的最优值分别记为
和
,则准则
的权重区间即为
。接下来,利用(1.2) 式即可计算出任意两个准则
和
的权重区间的优先度。最后,利用(1.3)~(1.5)式即可得到准则
的重要性排序。
3. 区间优先度的简化计算与一致性尺度变换
由上一节(1.2)~(1.5)式可知,要得到0-1矩阵(1.5),必须先根据(1.2)式计算区间优先度,再将其与0.5比较大小。事实上,该计算过程可以得到极大的简化,下面定理表明,上述过程可归结为:直接比较权重区间中心点的大小关系。除非特别说明,以下总假设有n个准则
。
定理1. 设准则
和
的权重区间分别为
和
,且两个权重区间的中心点分别为
,
。如果
,则
。
证明. 下面根据权重区间
和
的端点的相对位置关系分两种情况展开讨论。
情况1:
。
此时,
,因此,
(3.1)
又因为
,所以
,进而有
(3.2)
联合(3.1)和(3.2)可得
(3.3)
因为
,所以
(3.4)
根据区间优先度的计算公式(1.2)式,(3.4)式等价于
(3.5)
情况2:
。
此时,
,由此可知
(3.6)
又因为
,此时仍有
(3.7)
联立(3.6)式和(3.7)式可得
(3.8)
另一方面,因为
,所以有
(3.9)
综合(3.8)式和(3.9)式可得
(3.10)
证毕。
同理可证明下面的结论。
命题2. 设准则
和
的权重区间分别为
和
,且两个权重区间的中心点分别为
和
。如果
,则
。
为说明本文定理1和命题2相较于已有方法的优势,下面选取与[8]中相同例子说明定理1和命题2的实用性和有效性。
例1. 某消费者准备购买一辆汽车,其考虑的属性为:质量(
),价格(
),舒适度(
),安全性(
),款式(
)。该消费者认为最优准则为价格(
),最劣准则为款式(
)。最优准则与各准则的比较向量(BO)为:
;各准则与最劣准则的比较向量(OW)为:
。求解模型(2.1)~(2.3)可知,准则
到
的权重区间分别为:
进一步计算可得,准则
到
的权重区间的中心点分别为:
(3.11)
下面利用定理1和命题2,给准则
到
排序。首先,由(3.11)可知:
所以,由定理1和命题2可得:
,
,
,
。再由(1.4)式可知:
。同理可求得其它
。整理可得如下0-1矩阵
(3.12)
计算各行之和可知,准则
到
的排序为:
。
注1. 在例1中,利用本文定理1和命题2得出的结论与[8]中的结果完全一致,其区别在于:本文“绕过”了区间优先度(1.2)式的复杂计算过程,更便捷的获得了0-1矩阵(3.12)。
然而,需要强调的是,定理1和命题2虽然规避了区间优先度(1.2)式的计算过程,但仍需通过求解优化模型(2.1)~(2.3)来获取每个准则的权重上界和下界,这无疑是限制了该方法的推广应用。那么,能否进一步“绕过”求解优化模型这一复杂过程,而较为准确地估计出每个准则的权重区间的中心点呢?
为此,需要分析模型(2.1)产生多解的原因。回顾本文命题1可知,模型(2.1)之所以产生多解是因为最初在进行准则比较时产生了不一致现象,即:
。反之,如果准则之间的比较是完全一致的,那么,就可以利用比较向量
直接计算出各准则的权重
(3.13)
基于此,现在思考:是否可以通过某种变换,将所有不一致的比较转化为“形式上一致”的比较,从而利用(3.13)式的方法直接算出各个准则的权重。简言之,就是当出现
时,将
和
分别转化为
和
,使得
(3.14)
进而利用下式计算各准则的权重
(3.15)
那么,如何设计变换,将
和
分别转化为
和
使得(3.14)式成立呢?下面对此展开探讨,首先,给出以下定义(相关记号的含义与前文相同)。
定义2. 设对于准则
,有
,则称“
”为准则
的“一致性尺度”。
定义3. 设对于准则
,有
,
,
为实数集。如果存在映射
和
使得
,
且
,则称
和
为准则
的“一致性尺度变换”。
例2. 假设
,
,因为
,所以准则
的一致性尺度为
;由于
,所以与准则
有关的比较并不具有一致性。为了将该比较“打造”成的具有一致性,需要将一致性尺度“9”缩小为“8”,与此同时,
和
也应随之分别变换为
和
,并使其满足
。此时的一致性尺度变换可选为:
,
,其中
。
(1) 一致性尺度变换I
受例2的启发,考虑让
和
乘以同一个比例系数r,使得
,从该式求解出“r”可得
(3.16)
从而
(3.17)
且
(3.18)
下面以例1为例说明上述一致性尺度变换的有效性。在例1中,
,显然,舒适度(
)与安全性(
)两个属性的成对比较存在不一致,因此,将这两个各个准则的相关比较
,
与
,
按(3.17)式和(3.18)式分别进行调整,修正后的结果为:
(3.19)
再利用(3.15)式可得
(3.20)
现在回顾例1的(3.11)式,利用与[1]中的区间分析法和完全相同的数据,得到的各准则的权重区间的中心向量为
(3.21)
对比
与
可知,二者的总偏差为0.0513,最大偏差为0.0198,平均偏差为0.01026。从定量角度看,两种方法的结果平均相差约1%,最大偏差不超过2%;从定性角度看,两种方法得到的结论都是:
。综上,用“一致性尺度变换I”得到的排序对最终结果没有影响。
(2) 一致性尺度变换II
在“一致性变换I”中,
和
是按相同的比例r进行缩放的,下面考虑
和
分别按比例
和
进行缩放(
),此时,只需寻找
和
使得
(3.22)
由此解得
(3.23)
(3) 一致性尺度变换III
在(3.22)式中,
表示缩放的比例系数
和
与
和
的大小成比例。另一种思路与一致性尺度变换II相反,对于成对比较数值较大的给予较小幅度的调整;对于成对比较数值较小的给予较大幅度的调整。据此,对(3.22)修正如下
(3.24)
求解(3.24)可得
(3.25)
(4) 一致性尺度变换IV
如果对
和
的“绝对大小”进行相等幅度的调整,则可得到以下一致性尺度变换
(3.26)
联立(3.26)式中两式可得
(3.27)
求解(3.27)式对应的方程可得
。
(5) 一致性尺度变换V
如果对
与
差的“相对大小”进行相等幅度的调整,则可得到以下一致性尺度变换
(3.28)
联立(3.28)式中的两式可得
(3.29)
求解(3.29)式对应的方程可得
。
(6) 一致性尺度变换VI
受“几何平均算子”的计算方式的启发,下面考虑将
和
按指数r进行缩放,具体算法如下:
(3.30)
求解(3.30)式可得
(3.31)
以上给出了六种一致性尺度变换,哪一种相对更优呢?下面结合例1,将各一致性尺度变换的结果与通过模型(2.1)得到的权重区间的中心点进行对比分析。将
和
依次代入(3.22)~(3.31)式,其结果如表1所示。
Table 1. Deviation analysis of different consistency scale transformations (1)
表1. 不同一致性尺度变换的偏差分析(1)
方法 |
权重值 |
总偏差 |
平均偏差 |
一致性尺度变换I |
(0.2205, 0.4411, 0.1559, 0.1273, 0.0551) |
0.0514 |
0.0103 |
一致性尺度变换Ⅱ |
(0.2264, 0.4528, 0.1132, 0.1509, 0.0566) |
0.0943 |
0.0189 |
一致性尺度变换Ⅲ |
(0.2197, 0.4394, 0.1554, 0.1306, 0.0549) |
0.0482 |
0.0096 |
一致性尺度变换IV |
(0.2198, 0.4395, 0.1554, 0.1303, 0.0549) |
0.0484 |
0.0097 |
一致性尺度变换V |
(0.2214, 0.4428, 0.1565, 0.1239, 0.0553) |
0.0548 |
0.011 |
一致性尺度变换VI |
(0.2211, 0.4429, 0.1566, 0.1238, 0.0554) |
0.0549 |
0.011 |
区间分析法 |
(0.2024, 0.4609, 0.1536, 0.1345, 0.0521) |
0 |
0 |
由表1可知,上述六种一致性尺度变换的平均偏离差均在1%左右,其中方法三与方法四稍优,但差异极小。因此,用上述一致性尺度变换的方法估计权重区间的中心点,不仅具有较好的准确性,还能有效地规避求解非线性优化模型(2.1)。
此外,本文还发现,上述一致性尺度变换的方法除了可以用于近似估计权重区间的中心,在精度要求不高时,亦可以直接用来估计准则权重。下面以[7]中的一个例子说明。
例3. 设有三个准则:
。最优准则为
,最劣准则为
,
与各准则比较的结果(BO)以及各准则与
比较的结果(OW)分别为:
(3.32)
显然,本例存在比较“不一致”问题,但由于仅包含三个准则,因此模型(2.1)存在唯一解,将(3.32)式中的数值分别代入(3.22)~(3.31)式,其结果如表2所示。
Table 2. Deviation analysis of different consistency scale transformations (2)
表2. 不同一致性尺度变换的偏差分析(2)
方法 |
权重值 |
总偏差 |
平均偏差 |
一致性尺度变换I |
(0.0742, 0.3320, 0.5938) |
0.0135 |
0.0027 |
一致性尺度变换Ⅱ |
(0.0769, 0.3077, 0.6154) |
0.0620 |
0.0124 |
一致性尺度变换Ⅲ |
(0.0730, 0.3426, 0.5843) |
0.0112 |
0.0022 |
一致性尺度变换IV |
(0.0730, 0.3431, 0.5839) |
0.012 |
0.0024 |
一致性尺度变换V |
(0.0753, 0.3219, 0.6028) |
0.0336 |
0.0067 |
一致性尺度变换VI |
(0.0753, 0.3222, 0.6025) |
0.0336 |
0.0066 |
区间分析法 |
(0.0714, 0.3387, 0.5899) |
0 |
0 |
可以看到“一致性尺度变换III和IV”仍然具有相对较好的效果,这里利用“一致性尺度变换III”求得的
,由此算得“一致性尺度变换III”的一致性比率为
,与最优一致率0.058相比,相差甚小并且同样小于0.1,这也意味着使用“一致性尺度变换III”来直接估计权重具有较好地可行性。
注2. 当精度要求不高时,通过“一致性尺度变换”估计权重可大幅度减少计算量。需要注意的是,在含有多个指标层的决策问题中,一致性尺度变换法产生的偏差会逐层累积。设第一层指标为
,
下的第二层指标为
,以下各层指标记号以此类推。假设利用一致性尺度变换计算指标
的权重为
,偏差为
;计算指标
的权重为
,偏差为
。那么,
的权重的累计偏差为
(3.33)
因为
是
和
高阶无穷小,因此,
由此可见偏差的叠加效应。以此类推,第三层与第四层的最大偏差可类似分析。当面对的指标层层数过多时,可先在上层用原始优化模型计算,仅在最后一层或几层使用一致性尺度变换的估算方法,以减少最终偏差。
需要注意的是,估算方法与传统BWM求解过程的本质有所不同:估算方法将不一致比较转换为一致比较,并在此基础上计算权重;而BWM求解过程则是调整所有选项权重,使权重间的比值
尽可能接近决策者给出的原比较值
。
4. 一致性尺度变换法在电网关键物资选取中的实际应用
本节主要关注国家电网企业物资采购过程中现存的一个实际问题:关键物资选取。近年来,随着国家电网有限公司逐渐实行物资集中采购的管理模式,很多物资的实际采购价格与中标价格差距出现较大偏差,其主要原因在于:外部经济环境变化无常,各种电网物资的原材料价格在一定的范围内常有较大波动,而电网关键设备的申报、融资、招标需要经历较长的时间,招标时间相对于申报时间有较长的时间滞后,致使申报价格与招标执行价格经常出现较大的偏差,产生了不必要的融资高额利息,提高了设备的采购成本,且对电网的维护成本预算具有较大影响。反之,如果在采购过程中,物资价格被压得过低,则无法保证供应商在物资交付时的质量,轻则降低电力设施使用寿命,重则增大安全事故隐患。因此,针对电网物资的“价格预测”研究日益成为近些年的研究热点,详见[9]-[16]。
Table 3. Characteristic data of some supplies
表3. 部分物资特征数据
物资唯一码 |
需求数量 |
概算价格(元) |
中标单价(元) |
概中差 |
标准物料 |
产权属性 |
项目类型 |
工程性质 |
建设性质 |
资金属性 |
2305E03146 |
120 |
7500.00 |
6975.00 |
525.00 |
2 |
国网 |
技改项目 |
其他 |
改造 |
资本性 |
2308A03747 |
1800.00 |
124.74 |
125.03 |
0.29 |
1 |
国网 |
基建项目 |
线路 |
新建 |
资本性 |
2308C04279 |
1.00 |
113800.00 |
124300.00 |
10500.00 |
2 |
省网 |
基建项目 |
变电 |
新建 |
资本性 |
2308C04318 |
1.00 |
129150.00 |
146900.00 |
17750.00 |
2 |
国网 |
基建项目 |
变电 |
改造 |
资本性 |
2308C04453 |
1.00 |
183723.80 |
242950.00 |
59226.20 |
2 |
国网 |
基建项目 |
变电 |
新建 |
资本性 |
2308H04292 |
150.00 |
128 |
125.03 |
2.97 |
1 |
省网 |
基建项目 |
变电 |
增容 |
资本性 |
2308M04004 |
1.00 |
20000.00 |
14650.00 |
5350.00 |
1 |
国网 |
基建项目 |
变电 |
扩建 |
资本性 |
2311A03990 |
420.00 |
98.51 |
125.03 |
26.52 |
1 |
国网 |
大修项目 |
其他 |
修理 |
成本性 |
2311A03991 |
1330.00 |
440 |
409.69 |
30.31 |
1 |
国网 |
大修项目 |
其他 |
修理 |
成本性 |
2311A03992 |
700.00 |
440 |
409.69 |
30.31 |
1 |
国网 |
大修项目 |
其他 |
修理 |
成本性 |
2311A03995 |
15.00 |
7000.00 |
6969.84 |
30.16 |
2 |
国网 |
大修项目 |
其他 |
修理 |
成本性 |
2311A04026 |
15.00 |
6000.00 |
5969.79 |
30.21 |
2 |
国网 |
大修项目 |
其他 |
修理 |
成本性 |
2311A04114 |
2.00 |
5000.00 |
4660.00 |
340.00 |
2 |
国网 |
大修项目 |
其他 |
修理 |
成本性 |
2311B01226 |
0.45 |
657660.00 |
621582.02 |
36077.98 |
1 |
国网 |
基建项目 |
变电 |
新建 |
资本性 |
2311B01227 |
0.35 |
111870.00 |
152754.37 |
40884.37 |
1 |
国网 |
基建项目 |
变电 |
新建 |
资本性 |
2311B01228 |
0.57 |
15820.00 |
22369.84 |
6549.84 |
1 |
国网 |
基建项目 |
变电 |
新建 |
资本性 |
2311B02836 |
2200.00 |
6 |
5.9 |
0.10 |
2 |
国网 |
基建项目 |
其他 |
新建 |
资本性 |
2311C07528 |
1.00 |
385000.00 |
352560.00 |
32440.00 |
2 |
国网 |
技改项目 |
其他 |
改造 |
资本性 |
2311D01835 |
1.00 |
37197.65 |
36825.89 |
371.76 |
1 |
国网 |
技改项目 |
其他 |
改造 |
资本性 |
然而,目前国家电网采购部门每年需要采购的物资种类有上千种,针对每种物资建立价格预测模型并不现实。事实上,针对其中的关键物资进行价格预测更具现实意义。那么,如何在众多物资中筛选出“关键物资”呢?下面以江苏省电网公司2024年第一季度和第二季度的采购数据为例说明本文提出的“一致性尺度变换法”在关键物资选取中的应用。这两个季度该企业共采购物资4000余种,限于篇幅,表3选取其中的20种物资进行分析,同理可分析更多物资。注意到“关键物资选取”本质上是一个多属性评价问题,接下来,依次采取以下步骤:
第一步:确定关键物资的评价指标
经过与电网企业采购部门研讨,此处选择以下10个评价指标:C1:需求数量;C2:概算单价;C3:中标单价;C4:概算价格与中标单价差值;C5:标准物料;C6:产权属性;C7:项目类型;C8:工程性质;C9:建设性质;C10:资金属性。其中,概算价格与中标单价差值简称为“概中差”。
第二步:确定最优指标和最劣指标,并进行成对比较
通过业内调研,这里确定的最优指标为“概算价格与中标价格差值(
)”,最劣指标为“产权属性(
)”。最优指标与其他指标的比较向量和各指标与最劣指标的比较向量分别为
(4.1)
第三步:确定各指标的权重
由比较向量
和
可知,该比较存在明显的不一致性,例如:
,
,而
。因此,模型(2.1)可能存在多解。下面以本文提出的“一致性尺度变换III”为例,用其直接估计出各准则的权重,可得
第四步:数据标准化
由于准则
至
均是语言性评价,因此,首先将其进行量化处理,对于“产权属性(
)”,设置“国网”为2分,“省网”为1分;对于“项目类型(
)”,设置“技改项目”为1分,“基建项目”为2分,“大修项目”为3分;对于“工程性质(
)”,设置“变电”为3分,“线路”为2分,“其他”为1分;对于“建设性质(
)”,设置“改造”为1分,“新建”为2分,“修理”为3分,“增容”为4分,“扩建”为5分;对于“资金属性(
)”,设置“成本性”为1分,“资本性”为2分。然后利用“线性比例变化”法对表3中数据进行标准化,得到
最后,利用准则权重对标准化后的矩阵进行线性加权求和,可得20中物资的得分分别为:
,
;
,
;
,
;
,
;
,
;
,
;
,
;
,
;
,
;
,
。
若选取得分排名前三的物资作为关键物资,则关键物资为:2308C04453,2311B01227,2311C01835。
注3. 本算例采取的权重估计方法为“一致性尺度变换III”,事实上,也可采用其他一致性尺度变换法。另外,由于本算例需要计算各物资的得分,因此,需要借助相关指标的权重值。如果仅需了解各指标的排序,则可在估计出权重区间中心值之后,再利用0-1矩阵(1.5)和(1.6)式即可获得最终排序。
5. 结语
本文针对非线性BWM的多解问题,改进了传统的区间分析法,通过比较权重区间的中心值规避了区间优先度的复杂计算过程。此外,基于“尺度变换”的思想,给出了几类权重中心的估计方法,从而绕过了求解非线性BWM的复杂过程,而且偏差仅在1%左右。由于本文所给出的方法不要求决策者掌握过多的最优化理论,因此,具有广阔的推广价值。最后,将一致性尺度变换法成功应用于电网企业的关键物资选取问题,这也充分体现了本文所建立的方法的有效性和便捷性。同时,需要指出的是:除本文所建立的几类一致性尺度变换之外,是否还有其它更好的一致性尺度变换,以及一致性尺度变换在其他领域的应用等,都值得后期深入探索。
基金项目
江苏省青蓝工程,国家社科基金项目(项目号:24BH178)。