基于博弈论-TOPSIS 法的船型综合评价方法
2020-12-27张明霞韩丹赵桐鸣
张明霞,韩丹,赵桐鸣
大连理工大学 船舶工程学院,辽宁 大连 116024
船舶的建造投资大、使用周期长、风险高,其设计成功与否,很大程度上取决于方案的设计。因此就要求设计者在船型方案的初步设计阶段,对各可行方案进行准确的船型综合评价。如此,赋权方法是否可以兼顾专家主观意见和指标客观性、评价方法的优劣,将直接影响决策者的判断。因此使用全面综合的赋权方法以及合理的评价模型,能够改善船型综合评价工作的质量,更好地为决策者提供建议。
以往的船舶综合评价工作往往采用一些传统评价方法。李劲松等[1]、汪敏等[2]、姚雷等[3]采用层次分析法对船型方案选型提出决策建议,但此方法依赖专家意见,评价结果过于主观。刘元丰等[4]、高丹等[5]将模糊综合评价法(fuzzy comprehensive evaluation,FCE)分别应用于船舶航行安全评价、船舶溢油事故定级方面,验证了模糊综合评价法在船舶工程领域的适用性,但模糊评价法同样对专家的经验和判断依赖性较强。熊云峰[6]综合考虑了客观因素和主观因素,将得到的复合权重与TOPSIS 结合,建立了复合权重 TOPSIS 船舶性能综合评价方法;张明霞[7]对主客观两类赋权法采用动态组合得到改进复合权重,并将之与TOPSIS 方法结合,建立了改进复合权重TOPSIS评价模型,完成了算例中7 艘修井船的评价,证明了此方法在船型综合评价中的适用性。Pelorus[8]将TOPSIS 与AHP 相结合,确定了一种适用于船舶运营商的成本效益决策工具。张明霞[9]对上述几种评价方法做了对比,发现改进复合权重TOPSIS法区分度最高。
在船型综合评价工作中,评价指标权重的确定是否科学合理对于评价工作结果影响很大。目前常用的船型综合评价方法,例如AHP、FCE 及EWM 等,以上单一方法确定的权重存在片面性。复合权重是将主观赋权法与客观赋权法2 种方法以一定方式进行组合,一定程度上兼顾了决策者主观意愿与评价对象的客观信息。而基于博弈论的组合权重赋权原理,是对采用多种不同赋权方法得到的权重进行博弈集结,来协调不同赋权方法之间的不一致,最终达到一个均衡满意的结果,能够更科学、全面、客观地解决权重的确定问题。Hu[10]基于博弈论对遥控装置的网络改造和验收技术优化进行了分析和探讨,对其应用条件和技术难点进行了分析,在实际应用中取得了良好的效果;陈衍泰[11]运用合作博弈的原理,基于组合评价以及具有相同属性的单一评价方法,提出了应用合作博弈确定组合评价权重系数的方法;何俊[12]将博弈论和灰色关联分析的评估方法引入到雷达抗干扰能力的综合评估上,得到了科学的权重;周建国[13]采用基于博弈论的组合权重确定方法来确定评价指标的权重,并与灰色关联度理论结合,对我国六大区域电力市场的运营效果进行了综合评价。
本文将基于博弈论的组合权重确定方法应用于船型综合评价工作中,通过主观赋权的层次分析法[2]、客观赋权的熵权法[14]和智能赋权的BP-神经网络法[15]分别得到3 个权重向量,对得到的权重进行博弈集结得到最优权重。然后,利用TOPSIS 评价模型对确定指标权重后的船型方案进行综合评价。
1 基于博弈论的组合权重确定方法
综合权重的集结模型可分为博弈集结模型、团队集结模型及群体集结模型3 类[13]。其中博弈
集结模型的基本思想是根据不同赋权方法分别求得的权重,假设这些权重向量是相互独立的,在其中寻找协调一致的关系,可以有效地减少主观赋权法计算权重的主观性,提高客观赋权法的科学性,寻找不同赋权方法的均衡和协调。寻找组合权重和各个权重之间的最小化偏差,最终得到发反映专家意见及指标本身属性的组合权重。由此推导出本文的对策模型即博弈集结模型。
1.1 基本假设
设U={U1,U2,···,Um}为方案集,其中Ui代表第i(0 <i≤m)个可行方案,xij表示第i(0 <i≤m)个方案的第j(0<j≤n)个指标值,可建立初始决策矩阵B={xij|i=1,2,···,m;j=1,2,···,n}。
1.2 归一化处理
在建立基本假设之后,需要对指标进行归一化处理。
方案的评价指标通常分为2 类,一类是效益型指标,一类是消耗型指标。而指标间具有不同的量纲与数量级,没有可比性,所以需要对评价指标进行标准化处理。首先需对评价指标进行无量纲处理,公式如下。
式中:xij表示第i(0 <i≤m)个方案的第j(0 <j≤n)个指标值;r表示x的标准化指标值;r′表示r的标ijijijij准化指标值。
最终得到的标准化决策矩阵R为
1.3 博弈论赋权法
经过基本假设和归一化指标后,设使用L种不同的赋权方法计算得到L个权重向量wk=[wk1,wk2,···,wkn],(k=1,2,···,L),构建一个基本权重向量集:
这些权重向量的任意线性组合均构成一个可能权重集:
式中: αk为权重系数;w为可能的权重向量。
寻找最满意的权重向量w∗,就是使 w与各wk的离差极小化:
根据矩阵的微分性质得到式(3)的优化一阶导数条件为
2 逼近理想解排序法(TOPSIS)
TOPSIS 的原理是借助评价对象与方案正负理想解的相对距离对其进行排序的。正理想解的每一个指标通常取为评价方案中的最好的值;负理想解为待评价方案中各指标最劣的值。TOPSIS法通过检测待评价方案与正负理想解的相对趋近度来判断方案的好坏,显然,越趋近正理想解方案越优。
张明霞[7]建立了改进复合权重TOPSIS 评价模型,将改进复合权重TOPSIS 法用于船型综合评价工作中,完成了7 艘修井船的评价,证明了此方法在船型综合评价中的适用性。又在文献[9]将此方法与层次分析法、模糊综合评判法以及模糊层次分析法对比,发现TOPSIS 法的方案评价得分区分度最高,即越好的方案得分更高,避免出现忽略最优方案的现象。
由此本文选择TOPSIS 法作为本文评价方法,正(负)理想解公式为
式中:R+为正理想解;R-为负理想解。
则各评价方案与正负理想解的距离表示为
那么,贴近度ci的计算公式可表达为
待评价方案的贴近度ci取值范围为(0,1),越接近1 说明该待评价方案距离正理想解越近,方案越优。
在本文中ci值的大小代表船型方案的好坏,可作为方案评价得分。
3 算例验证
本文按照图1 流程对算例进行计算。
图1 算法计算流程
本文算例为13 艘消防船,部分数据通过查阅文献[16]得到。13 艘船分别是澳门消防船(澳消)、珠海消一号(珠消一)、KASUMI(KAS)、香港5 号(香5)、珠江号、MIYAGO(MIY)、KUSUNOKI(KUS)、布引号、横滨号、NAMIHAYA(NAM)、MINOO(MIN)、深消一二号(深一二)、深消三号(深消三)。
综合考虑消防船的经济性、工作性能与绿色度要求,建立造价(P)、年营运成本(S)、灭火能力(F)、航速(Vs)、能效设计指数(energy efficiercy design index,EEDI)(EE)这5 个评价指标,其中P、EEDI是消耗型指标,值越小越好,其他指标为效益型指标,值越大越好。
航速(Vs)为消防船服务航速,利用海军部公式来估算;而消防船的灭火能力(F)受水炮的流量和射程的直接影响;EEDI 反映了消防船作业与航行时二氧化碳的排放情况。
消防船本身的应用场景导致消防船都有比较相似的技术特点:良好的快速性以及持续的灭火能力。消防船需要较高航速来保证在火灾尚未造成严重损失时到达火灾现场;而很多火灾的燃烧时间都会持续几天几夜,因此,持续的灭火能力也是一个消防船所必须的[16]。
为使图表更加清晰直观,对部分船名及指标进行简化。消防船各项指标初始数据如表1 所示。
表1 消防船各项指标初始数据[16]
3.1 层次分析法(AHP)
3.1.1 构造判断矩阵
判断矩阵为层次分析法的信息基础,其元素值能真实反映出评价对象各个指标之间的相对重要程度。为表述每一层中各要素相对其上层某要素的相对重要程度,构造判断矩阵[2]如下:
式中aij为针对准则层Ck而言,指标Ai相对指标Aj重要程度的数值,通常采用1~9 的相对比较标度[2]。
因在本消防船船型综合评价中,指标的重要度为:P(造价)不太重要;S(年营运成本)不太重要;EEDI 值得重视;F(消防能力)很重要;Vs(航速)重要。
以上述指标重要度为例:两指标同等重要,标度为1;某指标比另一指标略微重要,标度为3;某指标比另一指标明显重要,标度为5;某指标比另一指标强烈重要,标度为7;各标度倒数表示反比较意义。
故建立相应的判断矩阵A如表2 所示。
表2 判断矩阵A
为检验判断矩阵中各元素的协调性,需要进行一致性检验,通常一致性检验公式为
式中:n为评价指标的个数;CR为一致性检验指标;RI为平均随机一致性指标,取值与比较因子的数量有关,当n=5 时,RI=1.12。通常认为,当CR<0.1时,认为层次分析法判断矩阵有满意的一致性,否则必须重新调整判断矩阵直至符合一致性标准为止[2]。
经计算,矩阵A的λmax=5.364 5,RI=1.12,CR=0.081,CR<0.1,矩阵A通过一致性检验。
3.1.2 指标权重向量
w=[w1,w2,···,wj,···,wn]
权重向量为,确定各评价因素的权重,计算公式如下:
计算得到层次分析法的一组指标(P,S,EE,A,Vs)权重向量w1为
3.1.3 评价结果
由式(1)~(5)得到标准化决策矩阵:
将R与 权重向量w1代入式(11)得到层次分析法最终评价得分:
如表3 所示。
表3 层次分析法评价得分
3.2 熵权法(EWM)
熵权法的基本思想是根据待评价方案指标本身属性及特点来确定客观权重。信息熵Ej用来度量信息量的大小,建立在原始数据的基础上,客观性比较强[14]。本文采用了以信息熵Ej确定属性权重的方法,具体定义如下。
根据得到的标准化决策矩阵R计 算信息熵Ej:
因此,第j个评价指标标准化的熵权法权重系数值为
由式(12)~(14),计算得到一组熵权法确定的指标(P,S,EE,F,Vs)权重向量w2=[θ1,θ2,···,θn]为
3.3 BP-神经网络法
神经网络可通过学习和训练获取网络的权值和结构。本文以层次分析法得到的评价结果作为目标值,依靠Matlab 所携带的神经网络工具箱nntool 对样本进行训练,得到权值。
当输入的标准化方案指标值训练结束并达到要求的网络精度后,基于调整后得到的最终输入层到隐含层之间的连接权矩阵V,计算各输入层节点到所有隐含层节点之间连接权的绝对值之和,并归一化,得到m个指标的权重[15]:
式中vjl为各输入层节点到所有隐含层节点间的连接权值。
得到BP-神经网络法确定的指标(P,S,EE,F,Vs)权重向量w3为
3.4 博弈论法
分别通过层次分析法、熵权法和BP-神经网络法得到3 组权重向量w1、w2、w3,由式(1)~(5)可得到相应线性方程组的矩阵表示形式:
利用Matlab 求解线性方程组的最优解为α1=1.039 1,α2=0.474 3,α3=-0.529 7。
根据式(6)、(7)可得到博弈论组合权重向量w∗为
3.5 TOPSIS 法船型方案评价
下面分别采用改进复合权重法和基于博弈论的组合权重确定方法确定指标权重,将各评价指标权重与TOPSIS 方法结合,建立改进复合权重TOPSIS 评价模型与博弈论-TOPSIS 评价模型。将上述2 种方法应用到消防船船型的优选中,对13 艘消防船进行船型综合评价。
3.5.1 改进复合权重TOPSIS 法贴近度ci计算
由文献[7]可知:已经由3.1 与3.2 节求出层次分析法权重向量w1与熵权法权重向量w2,于是复合权重 ρ可表示为
动态权重偏好系数 εj为
wi为求得的组合权重值。
再由式(8)~(10),计算各方案的贴近度即评价得分ci如表4 所示。
表4 改进复合权重TOPSIS 法评价得分(ci)
3.5.2 博弈论-TOPSIS 法贴近度ci计算
根据3.4 节得到的博弈论组合权重,由式(8)~(10),计算各方案的贴近度即评价得分ci如表5所示。
表5 博弈论-TOPSIS 法评价得分(ci)
3.5.3 3 种方法结果比较
由于层次分析法、改进复合权重TOPSIS 法、博弈论-TOPSIS 法的5 组指标重要程度一致,因此将表3~5 的评价结果绘图,对3 种不同方法得到的13 艘消防船的评价结果进行对比分析,如图2 所示。
图2 3 种方法的评价得分比较
1)从图2 可以看出,3 种方法得分最高的均为香港5 号(香5)消防船,得分最低为MINOO(MIN)号消防船。即当各个衡量指标的重要度为:P(造价)不太重要、S(年营运成本)不太重要、EEDI 值得重视、F(消防能力)很重要、Vs(航速)重要时,香港5 号(香5)消防船船型最优,MINOO(MIN)号消防船最劣;
2)从图2 可以看出,3 种方法的评价得分曲线趋势一致,其中博弈论-TOPSIS 法与改进复合权重TOPSIS 法这2 种方法的评价得分曲线起伏较大,即区分度较大;层次分析法的方案评价得分曲线趋于平稳,即区分度较小;
3) 13 艘消防船的博弈论-TOPSIS 法与改进复合权重TOPSIS 法相比,评价得分区分度更大,且具有较优方案得分更高,较差方案得分更低的特点。
4 结论
1)针对主、客观赋权法的不足,构建了结合层次分析法-熵权法-BP-神经网络法的博弈论组合权重赋权模型,得到了基于博弈论确定的组合权重。
2)采用博弈论-TOPSIS 法对方案进行评价,结果表明,与层次分析法和改进复合权重TOPSIS 法得到的评价结果趋势一致,证明了博弈论法在船型综合评价中,指标权重确定方面是合理可行的。
3)改进复合权重TOPSIS 法、博弈论-TOPSIS法这2 种方法的方案评价得分区分度较大,层次分析法的方案评价得分区分度较小。其中,结合层次分析法-熵权法-BP-神经网络法的博弈论赋权法结合了多种赋权法得到组合权重,更具合理性。
4)改进复合权重TOPSIS 法评价得分与博弈论-TOPSIS 法评价得分相比,博弈论-TOPSIS 法的方案区分度更大。并且由于博弈论-TOPSIS法与复合权重TOPSIS 法相比,具有较优方案得分更高,较差方案得分更低的特点。故应用博弈论-TOPSIS 法进行船型选优,可以使最优方案更加突出。所以当实际应用中船型数量过多时,应用博弈论-TOPSIS 法进行船型选优时会有更大优势。