前言
DEMATEL-ISM联用的两种方式
第一、基于综合影响矩阵T求截距的方式。
第二、基于中心度与原因度绝对值求偏序的方式。
DEMATEL、ISM能联用根本原因是两者都是基于有向图网络。
DEMATEL-ISM联用流程图,及基础理论说明
DEMATEL理论基础
DEMATEL是针对复杂有向网络的一种研究方法。其核心有两部分:
影响的传递性
有向网络的度中心性
综合影响矩阵T求解的原理——影响的传递性
有向网络如上图:其中A对C无直接影响。很显然A对C的影响是由A对B施加影响,B对C施加影响。A对C的影响是通过B为中介获得的。
$ T=\mathcal{N}+\mathcal{N}^2+\mathcal{N}^3+\cdots+\mathcal{N}^k $
其中:$ \mathcal{N}\times\mathcal{N} $的意义为增加的间接影响。即矩阵与矩阵相乘得到的值为矩阵要素之间增加的间接影响。
$ \mathcal{N}^\infty $为零阵。意思是影响不停的传递下去,其值为零。
$ \mathcal{N}^\infty $为零阵。前提要求是主对角线的值全部为0
综合影响矩阵$ T= \sum\limits_{k=1}^\infty {\mathcal{N}^k} $的意义为最初的直接影响加上所有的间接影响。
综合影响矩阵$ T=\mathcal{N}(I-\mathcal{N})^{-1} $
有向图的度中心性原理——Degree Centrality
对于有向图中任意一结点发出有向边称为出度(out degree)或者出边,接受有向边称为入度(in degree)或者入边。
DEMATEL中最终的四个参数:影响度D、被影响度C、中心度M、原因度R
影响度D:即综合影响矩阵的行和,对于有向图网络中即所有出边值之和,特别要指出的结点对结点自身有出度。
被影响度C:即综合影响矩阵的列和,对于有向图网络中即所有入边值之和,特别要指出的结点对结点自身有入度。
中心度M:影响度D与被影响度C之和。
原因度R:即出度(影响度)减去入度(被影响度)。原因度也可以叫成结果度,由定义可知,当其值为负数的时候,要素体现为结果,当其值为正数的时候,体现为原因的因素。
基于截距方式获得关系矩阵(布尔矩阵,邻接矩阵)
截距方式获得关系矩阵是最常用的方法。其特点是被处理的矩阵同关系矩阵都是方阵
基于直接影响矩阵$\mathcal{O}$可以通过截距的方式得到关系矩阵。其形式如下:
$\mathcal{O}=[0]_{n \times n}$取截距后(其中截距的值为直接影响矩阵中矩阵值的取值范围)得到布尔邻接矩阵$\mathcal{A}=[a]_{n \times n}$
基于规范化矩阵$\mathcal{N}$可以通过截距的方式得到关系矩阵。其形式如下:
$\mathcal{N}=[n]_{n \times n}$取截距后(其中截距的值为[0,1]之间)得到布尔邻接矩阵$\mathcal{A}=[a]_{n \times n}$
基于综合影响矩阵$\mathcal{T}$可以通过截距的方式得到关系矩阵。其形式如下:
$\mathcal{T}=[t]_{n \times n}$取截距后(其中截距的值为[0,1]之间)得到布尔邻接矩阵$\mathcal{A}=[a]_{n \times n}$
对于任意方阵$\mathcal{Z}=[z]_{n \times n}$取截距$\lambda$后,得到关系矩阵$\mathcal{A}=[a]_{n \times n}$
则有:
$ a_{ij} = \left\{ \begin{array}{ll} \mathcal{0} & \textrm{if $ z_{ij}< \lambda $}\\ \mathcal{1} & \textrm{if $ z_{ij}≥\lambda$} \end{array} \right. $
DEMATEL与ISM基于截距方式的联用,主要是针对综合影响矩阵T取截距
T与直接影响矩阵O与规范化矩阵N最大的区别在于其已经计算完了影响的传递性。如果只是针对直接影响矩阵或者规范化矩阵N取截距则该方法不能称为联用的方法,其实质是模糊解释结构模型的一般方法。
模糊解释结构模型方法的地址:模糊解释结构模型在线计算(收费版本)
综合影响矩阵T取截距得到关系矩阵A的特点
由于直接影响矩阵O与规范化矩阵N的主对角线都为0,而综合影响矩阵T主对角线大都不全部为0
因此T取截距得到的关系矩阵A主对角线未必全部为0
基偏序方式获得关系矩阵(布尔矩阵,邻接矩阵)
偏序相关概念
$$\require{cancel} \require{AMScd} \begin{CD} D=\left[ d_{ij} \right]_{n \times m}@>偏序规则>>A=\left[a_{ij} \right]_{n \times n} \\\end{CD} $$
其中 $D=\left[ d_{ij} \right]_{n \times m}$ 为决策评价矩阵。$n$行$m$列。$n$代表评价对象(要素、方案、样本);$m$代表维度(准则、属性、目标)。
其中 $A=\left[ a_{ij} \right]_{n \times n}$ 为关系矩阵。是一个布尔方阵。$n$代表评价对象(要素、方案、样本)。
对于决策矩阵$D$中 $n$个要素的任何一列都具有严格的可比性。
偏序规则
对于含有m列的评价矩阵D,其中的任意一列即指标维度,具有同属性,可比较的前提。维度的这种优劣的比较至少有着两种属性。
数值越大越优,数值越小越差,称之为正向指标。记作p1、p2……pm。 数值越小越好,数值越大越差,称之为负向指标。记作q1、q2……qm。
对于决策矩阵$D$中的任意两行$x,y$
负向指标有 $d_{(x,p1)} \geqslant d_{(y,p1)} 且d_{(x,p2)} \geqslant d_{(y,p2)} 且 {\cdots}且d_{(x,pm)} \geqslant d_{(y,pm)}$ 同时有
正向指标有 $d_{(x,q1)} \leqslant d_{(y,q1)} 且d_{(x,q2)} \leqslant d_{(y,q2)} 且 {\cdots}且d_{(x,qm)} \leqslant d_{(y,qm)}$
符合上述规则,要素$x$与要素$y$的偏序关系记作:$x ≺ y$
$x \prec y$的意义为$y要素$优于(好于,牛逼于,帅于,猛于)$x要素$ 。
上述规则成为偏序规则。对于决策矩阵通过偏序规则可以得到关系矩阵 $A$
$$a_{xy}= \begin{cases} 1, x \prec y \\ 0, 其它 \end{cases} $$
取偏序简单例子
$$ 示例一: \begin{CD} D=\begin{array}{c|c|c|c|c|c|c}{M_{3 \times1}} &-(负向指标) \\ \hline A1 &1 \\ \hline A2 &2 \\ \hline A3 &3 \\ \hline \end{array} @>取偏序>> A=\begin{array}{c|c|c|c|c|c|c}{M_{3 \times 3}} & A1 &A2 &A3\\ \hline A1 &- & & \\ \hline A2 &1 &- &\\ \hline A3 &1 & 1 & -\\ \hline \end{array} \end{CD} $$把只有1列的决策矩阵$D$中的负向指标想象成排名,A1为第1名。关系矩阵$A$中 A2->A1即A2行A1列对应的单元格意思为A1比A2牛逼,即$A2 \prec A1$
$$ 示例二: \begin{CD} D=\begin{array}{c|c|c|c|c|c|c}{M_{17 \times2}} & 正向指标 & 正向指标 \\ \hline A1 &1.9223 &0.59336 \\ \hline A2 &2.86838 &0.16965\\ \hline A3 &1.38284 &0.22882\\ \hline \end{array} @>取偏序>> A=\begin{array}{c|c|c|c|c|c|c}{M_{3 \times 3}} & A1 &A2 &A3\\ \hline A1 &- & & \\ \hline A2 & &- &\\ \hline A3 &1 & & -\\ \hline \end{array} \end{CD} $$在三组数据中只有A3的两个属性值都小于于A1。关系矩阵$A$中 A3->A1即A3行A1列对应的单元格意思为A1比A3牛逼,即$A3 \prec A1$
五种取偏序的方法
取偏序之前的首要是要先判断每列的属性,即是正向指标(数值越大越优,数值越大越牛逼)还是负向指标(数值越小越优,数值越小越牛逼)
基于DEMATEL的度中心性以及影响的传递性原理。
通常是采取中心度M与原因度的绝对值|R|来求偏序。
其中原因度的绝对值|R|来自于原因度的属性,其值为正的是整体体现为原因,其值为负是整体体现为结果。
使用偏序法前时,其核心步骤是判断每列的属性(即正向指标,还是负向指标)。
DEMATEL-AISM的计算地址:DEMATEL-AISM(基于偏序)计算(收费版本)
上面的计算实例都有讲清楚每列的属性。
为了区分截距法与偏序法,偏序法中DEMATEL-AISM联用又称为DEMATTEL-AHDT
ISM、HDT、DSM三个名称:解释结构模型的三种名称:解释结构模型(ISM),哈斯图技术(HDT),依赖结构模型(DSM)
AHDT的英文全称Adversarial Hasse Diagram Technology 中文对抗哈斯图技术
取偏序的几个重要性质
一般来说,如DEMATEL-AISM联用方法中:
基于截距的方法,最后要素与要素之间的关系是因果关系
基于偏序的方法,最后要素与要素之间的关系一定是比较关系(优劣比较,好坏比较,重要程度比较)
只有一列取偏序得到的必定为全序,即层级拓扑表现为一条棍子(直链),对应的系统称之为完全刚性系统,得到的层级拓扑图不存在活动要素。
上面性质很好理解,在进行综合评价的时候,最终都会形成一列,并由该列进行比较排序。如:通常的考试会以总分来进行优劣排序,并给出结论。
高考选拔排序:高考选拔与AISM
上述定义的偏序规则得到的关系矩阵$A$即为可达矩阵$R$
即有$A=B=R$
上述定义的偏序规则得到的关系矩阵$A$中通常不包含回路。
即有$R=R'$
即有$S'=S$
上述定义的偏序规则得到的关系矩阵$A$如果存在回路,则对应的决策矩阵$D$存在若干行完全相等,完全相等的行构成一条回路。
在决策矩阵$D$中无重复行时,对抗拓扑层级图的最简方式如下:
其中$HS=S'=S$
$HS$称为哈斯矩阵
$S$称为一般性骨架矩阵(回路以菊花链方式表达)
$S'$称为骨架矩阵
运用$S+I$进行层级划分的实例请参看论文:基于对抗解释结构模型的军事训练方法可推广性评价模型
DEMATEL求权重的两种方法
D-ANP求权重的原理
AHP与ANP方法求权重的原理
层次分析法的英语术语为 Analytic Hierarchy Process,可以简称 AHP,是一种运用广泛的定性和定量结合在一起的分析决策研究方法。该方法是Saaty 于 20世纪 70 年代初,在为美国国防部研究一个课题时,提出的一种研究方法。该方法把要素组织成一个树形的层次结构,通过构建不同的层次的成对比较决策矩阵,计算得到各个要素的权值。AHP 的原理是把研究的对象的要素层次化,它首先 确定研究对象的总目标,根据总目标把研究对象分解为不同由若干要素组成的子系统,这些子系统是按照层次的组织方式形成一个层次结构,最终通过子系统里要素之间的成对比较矩阵得出所有要素的权重。
网络分析法的英文术语为 Analytic Network Process,简称 ANP,是美国 Saaty教授在 1996 年提出的一种决策方法,该方法的提出是基于层次分析法,是一种适应非独立递阶层次结构的方法。ANP相对于AHP而言,用网络结构代替了层次结构,同时会将要素间的相关性考虑进去,用非线性结构代替线性层次结构,还加入了反馈机制,并考虑到低层要素对于高层要素的支配作用。此后ANP作为一种多准则决策方法,逐渐成为学术界研究的焦点问题,并被广泛用于各领域,解决了诸多复杂决策问题。
ANP 和 AHP 有相似之处,因为两者的理论基础都是基于成对比较决策矩阵而来。两者应用的领域与解决的问题大多是无法精确计算,需要加入一定的主观判断的复杂的问题中。两种方法在实际操作中也是相近的。ANP 是 AHP 的延伸,而 AHP 是 ANP 的特例,ANP 与 AHP 是包含关系,AHP 其本质就是 ANP 里的最简单的一种情况。但 AHP 中没有考虑同级要素之间的相互影响、相互依赖、 相互支配的关系,而 ANP 是在充分考虑 AHP 的准则的基础上,进一步考虑准则层里各要素间的相互关系
为了解决ANP中复杂与繁琐的建模与计算问题,Saaty所在的美国国家科学基金会,基于 ANP 的计算原理开发了 Super Decisions(简称 SD)软件,随着该软件的不断完善,使得对复杂系统的评价决策更加方便快捷。SD 软件的出现,同 ISM系列相关计算软件的出现一样,大大的降低了研究中计算的工作量,使得研究者能够把主要精力放在建模本身,而不是在计算的细节上。
SD软件也可以进行AHP的运算。
决策矩阵是决策活动中用来分析和求解决策问题的依据,是各决策要素间定性或定量等相互关系的集中表现矩阵。在风险型决策中,决策矩阵又称为“损益表”、“风险型矩阵”或“决策表”,其基本构成元素主要包括决策变量、状态变量、概率和损益值。在多属性或多目标决策中,针对不同的决策问题,决策矩阵通常称为“决策矩阵”、“属性矩阵”、“属性表”和“判断矩阵”。对于判断矩阵,根据其元素的特性,可 进一步细分为“正互反型判断矩阵”、“互补型判断矩阵”、“模糊判断矩阵”和“语义判断矩阵”等
正互反判断矩阵
由Thurstone提出的两两成对比较判断矩阵(通常称为正互反判断矩阵)被广泛用来处理多准则决策(多目标决策和多属性决策)中的主客观判断信息。两两比较决策矩阵是层次分析法 AHP 和网络分析法 ANP 的核心。也是AHP与ANP的基础假设。脱离了正互反判断矩阵,就无从谈起AHP后面的一致性校验。
两两成对比较判断矩阵(正互反判断矩阵)定义描述如下:
$$\begin{array} {c|c|c} {M_{n \times n}} &A_1 & A_2 & \cdots &A_n \\ \hline A_1 & - & a_{12} & \cdots & a_{1n} \\ \hline A_2 & a_{21} & - & \cdots & a_{2n} \\ \hline \vdots & \vdots & \vdots & - & \vdots \\ \hline A_n & a_{n1} & a_{n2} & \cdots & -\\ \hline \end{array}$$
该决策矩阵可简写为 $A=[a_{ij}]_{n \times n}$
决策矩阵 A的元素具有非负性,即 矩阵内任意值不为负数。$a_{ij} >0$
互反性指的是 $a_{ij}= \frac {1}{a_{ji}}$
DEMATEL中的直接影响矩阵$O$
AHP/ANP的基本假设是正互反判断矩阵,对应的是无向图(特殊的有向图)
为了对比把直接影响矩阵定义为$A$其形式描述如下:
$$\begin{array} {c|c|c} {M_{n \times n}} &A_1 & A_2 & \cdots &A_n \\ \hline A_1 & 0 & a_{12} & \cdots & a_{1n} \\ \hline A_2 & a_{21} &0 & \cdots & a_{2n} \\ \hline \vdots & \vdots & \vdots & 0 & \vdots \\ \hline A_n & a_{n1} & a_{n2} & \cdots & 0\\ \hline \end{array}$$
从两者的定义上可以看出,对于一个决策单元来说。
AHP需要进行$\frac {n(n-1)}{2}$次比较判断。主对角线一般取1,表示自己同自己比较,重要程度是一致的。
DEMATEL需要进行${n(n-1)}$次比较判断。主对角线取0。
AHP、ANP、DEMATEL原始矩阵对比
设:三者都分为三大类,每类分别有三个要素。
AHP的结构如下:
$$\begin{array} {c|ccc|ccc|ccc} {M_{9 \times 9}} &A1 & A2 & A3 &B1 &B2 &B3 &C1 &C2 &C3\\ \hline A1 &1 & \frac{1}{2} & 3 & & & & & & \\ A2 &2 & 1 & \frac{1}{4} & & & & & & \\ A3 &\frac{1}{3} & 4 & 1 & & & & & & \\ \hline B1 & & & &1 & \frac{1}{3} & 3 & & & \\ B2 & & & &3 & 1 & \frac{1}{7} & & & \\ B3 & & & &\frac{1}{3} & 7 & 1 & & & \\ \hline C1 & & & & & & &1 & \frac{1}{3} & 3 \\ C2 & & & & & & &3 & 1 & \frac{1}{4} \\ C3 & & & & & & &\frac{1}{3} & 4 & 1 \\ \hline \end{array}$$
$$\begin{array} {c|c|c} {M_{3 \times 3}} &A & B & C \\ \hline A &1 & \frac{1}{2} & 1 \\ \hline B &2 & 1 & 1 \\ \hline C &1 & 1 & 1 \\ \hline \end{array}$$
ANP的结构如下:
$$\begin{array} {c|ccc|ccc|ccc} {M_{9 \times 9}} &A1 & A2 & A3 &B1 &B2 &B3 &C1 &C2 &C3\\ \hline A1 &- & \frac{1}{2} & 3 & 1 & 1 &1 &1 & 2 &2 \\ A2 &2 & - & \frac{1}{4} & & & & & & \\ A3 &\frac{1}{3} & 4 & - & & & & & & \\ \hline B1 & 1& & &- & \frac{1}{3} & 3 & & & \\ B2 & 1 & & &3 & - & \frac{1}{7} & & & \\ B3 & 1& & &\frac{1}{3} & 7 & - & & 1 & \\ \hline C1 & 1 & & & & & &- & \frac{1}{3} & 3 \\ C2 &\frac{1}{2} & & & & & &3 & - & \frac{1}{4} \\ C3 & \frac{1}{2} & & & & 1 & &\frac{1}{3} & 4 & - \\ \hline \end{array}$$
注:ANP中的原始矩阵,后来取消了决策单元内必须两两比较的部分,即A1与A2可以存在不能比较的部分。
DEMATEL的结构如下:
$$\begin{array} {c|ccc|ccc|ccc} {M_{9 \times 9}} &A1 & A2 & A3 &B1 &B2 &B3 &C1 &C2 &C3\\ \hline A1 &0 & 3 & 3 & 1 & 1 &1 &1 & 2 &2 \\ A2 &3 & 0 & 3 & 1& 4 & 2 &4 & 3 &3 \\ A3 &3 & 4 & 0 & 1 & 0 & 0 & 0& 0 & 0\\ \hline B1 & 1& 2 & 1 &0 & 0 & 3 & 0& 4 &2 \\ B2 & 2 &1 & 1&3 & 0& 4 & 1& 3& 4\\ B3 & 1& 1&1 &1& 3 & 0 &3 & 1 &4 \\ \hline C1 & 1 & 4 & 1& 1& 1 & 1&0 & 1 & 3 \\ C2 &1 & 3 &1& 0& 1 &1 &3&0 & 4 \\ C3 & 1 & 1& 1& 1& 1 & 2&3 & 4 & 0 \\ \hline \end{array}$$
注:DEMATEL中的原始矩阵,用得最多的是0-4分制,多人打分直接把分数相加即可。直接影响矩阵中,主对角线的必为0。
从AHP与ANP求权重的过程,看D-ANP求权重的过程
一个决策单元的AHP求权重的例子(不包含一致性校验):
判断矩阵假设如下:$$\begin{array} {c|c|c} {M_{3 \times 3}} &A & B & C \\ \hline A &1 & \frac{1}{2} & 1 \\ \hline B &2 & 1 & 1 \\ \hline C &1 & 1 & 1 \\ \hline \end{array}$$
每列单独归一化(求出比重):
$$\begin{array} {c|c|c} {M_{4 \times 3}} &A & B & C \\ \hline A &1 & \frac{1}{2} & 1 \\ \hline B &2 & 1 & 1 \\ \hline C &1 & 1 & 1 \\ \hline 列和 &4 & 2.5 & 3 \\ \hline \end{array}$$
$$\begin{array} {c|c|c} {M_{3 \times 3}} &A & B & C \\ \hline A & 0.25 & 0.2 & 0.33333 \\ \hline B & 0.5 & 0.4 & 0.33333 \\ \hline C & 0.5 & 0.4 & 0.33333 \\ \hline \end{array}$$
按行求和:
$$\begin{array} {c|c|c} {M_{3 \times 3}} &行和 \\ \hline A & 0.78333333 \\ \hline B & 1.23333333 \\ \hline C & 2.01666667 \\ \hline \end{array}$$
其比重就是权重
$$\begin{array} {c|c|c} {M_{3 \times 3}} & 权重 \\ \hline A & 0.194215 \\ \hline B & 0.305785 \\ \hline C & 0.5 \\ \hline \end{array}$$
由上述求权重的过程注意观察两个问题:
主对角线为什么取1而非0。
第一步是先列归一化。而不是行归一化,再求列和,再计算权重。
ANP求权重的例子及流程:
基于FISM-ANP-灰色聚类的软件项目开发风险评价研究(页码第53页开始)流程图如下:
ANP求权重与D-ANP求权重流程的比较:
从流程图来看,其中的综合影响矩阵$T$等同于ANP中的未加权超矩阵。
D-ANP方法与ANP方法最大的不同在于,充分的考虑到了影响的传递性。
目前很多发表的关于D-ANP方法的论文,相当一部分是错的,它们是基于直接影响矩阵$O$来进行计算的。
加权超矩阵的计算及特殊情况处理
D-ANP方法中极限超矩阵是通过综合影响矩阵$T$每一列单独归一化(求百分比)得到。
当$T$的任意一列值全部为0的时候,对应加权矩阵的列全部为0
加权矩阵等同于AHP方法中单个决策矩阵求权重的第一步,即列归一化。具体见上面AHP中的算例。
通过极限超矩阵求权重的方法,也可以运用到AHP当中。
DEMATEL及联用的各种组合方式
DEMATEL | 常规DEMATEL计算 可选择归一化方法 | 免费下载流程图计算需付费 |
---|---|---|
比较型常规DEMATEL计算 选两个归一化方法 | 免费下载流程图计算需付费 | |
比较型常规DEMATEL计算 选一个归一化方法 | 免费下载流程图计算需付费 | |
比较型GREY-DEMATEL计算 | 免费下载流程图计算需付费 | |
三者比较型GREY-DEMATEL计算 | 免费下载流程图计算需付费 | |
截距型 | 基于综合影响矩阵$T$的联用。 | 截距型的运用 |
截距型DEMATEL-ISM联用 | 免费下载流程图计算需付费 | |
截距型DEMATEL-ISM联用其截距特点为不确定值,即矩阵中的值等于截距值时候可以取0或者1两种情况 | ||
DEMATEL-SecISM 区段截取,即双截距截取。 | ||
DEMATEL-TAISM(综合对抗解释结构模型) 带有统计特征,总能找到一个合适的结果。 | 免费下载流程图计算需付费 | |
DEMATEL-AISM 特征拐点层级拓扑图 | 免费下载流程图计算需付费 | |
DEMATEL-AISM 特征拐点对比 | 免费下载流程图计算需付费 | |
偏序型运用 | 基于{中心度M,原因度R}其中原因度取绝对值 | 具体见上 |
DEMATEL-AISM-Weight含权重计算 | 免费下载流程图计算需付费 | |
dematel-AHDT五种哈斯图 | 免费下载流程图计算需付费 | |
G-D-AHDT方法新 | 免费下载流程图计算需付费 | |
DEMATEL-ANP联用 | 本质还是基于T的运用 | 原理见上 |
DEMATEL-ANPD-ANP方法求权重 | 免费下载流程图计算需付费 |