DEMATEL-ISM

方法概述

DEMATEL-ISM 将 DEMATEL(决策实验室分析法)与 ISM(解释结构模型)相结合,充分利用两种方法的优势:DEMATEL 能够量化因素间的直接影响强度,计算中心度与原因度,识别关键因素;ISM 则通过传递闭包和层级划分,揭示系统的层次结构。融合方法弥补了各自不足,既能反映因素的重要性与因果属性,又能构建清晰的层级关系。

该方法的核心思想是:

  • 多位专家对因素间的直接影响程度进行打分,经加权聚合得到直接影响矩阵。
  • 采用 DEMATEL 计算综合影响矩阵,并通过设定阈值 λ 二值化得到邻接矩阵。
  • 基于该邻接矩阵,利用 ISM 的传递闭包算法计算可达矩阵,进而计算可达集、前因集、交集。
  • 通过迭代分层,将因素分配到不同层级,构建骨架矩阵并绘制层次结构图。
  • 同时保留 DEMATEL 的因果散点图(中心度-原因度图),从“重要性-原因性”和“层级结构”两个维度全面分析系统。

该方法广泛应用于供应链风险识别、技术创新因素分析、系统结构建模等复杂决策场景。

计算步骤

1. 构建原始数据矩阵

设有 \(n\) 个因素,\(K\) 位专家。每位专家给出一个 \(n \times n\) 的直接影响矩阵 \(X^{(k)}\),其中元素 \(x_{ij}^{(k)}\) 表示专家 \(k\) 认为因素 \(i\) 对因素 \(j\) 的直接影响强度(通常采用 0~4 分制:0=无影响,1=低影响,2=中影响,3=高影响,4=极高影响),且 \(x_{ii}^{(k)} = 0\)

2. 聚合多个专家评价

设专家 \(k\) 的权重为 \(w_k\)\(\sum_{k=1}^{K} w_k = 1\)),则综合直接影响矩阵 \(A = [a_{ij}]\) 为:

\[ a_{ij} = \sum_{k=1}^{K} w_k \cdot x_{ij}^{(k)} \]

3. DEMATEL 计算综合影响矩阵

  • 规范化直接影响矩阵:选择规范化方法(行和最大值法、列和最大值法、行和列和最大值法),计算规范化因子 \(s\),得到规范影响矩阵 \(B = A / s\)
  • 计算综合影响矩阵\[ T = B + B^2 + B^3 + \cdots = B (I - B)^{-1} \] 其中 \(I\) 为单位矩阵。实际计算中先求 \(I-B\) 的逆矩阵,再左乘 \(B\)

4. DEMATEL 关键指标计算

  • 影响度 \(D_i = \sum_{j=1}^{n} t_{ij}\)
  • 被影响度 \(C_i = \sum_{j=1}^{n} t_{ji}\)
  • 中心度 \(M_i = D_i + C_i\)(反映因素重要性)
  • 原因度 \(R_i = D_i - C_i\)\(R_i>0\) 为原因因素,\(R_i<0\) 为结果因素)
  • 因素属性:以中心度均值为界,将因素分为核心原因、一般原因、核心结果、一般结果四类。

5. DEMATEL 到 ISM 的转换

  • 确定阈值 \(\lambda\):可选择“均值+标准差”、“均值”或“自定义”方式。

  • 二值化得到邻接矩阵

    \[ a_{ij}' = \begin{cases} 1, & t_{ij} \ge \lambda \\ 0, & \text{otherwise} \end{cases} \] 同时添加自反性(对角线为 1)。

6. ISM 计算可达矩阵与分层

  • 可达矩阵:采用布尔运算迭代(Warshall 算法)计算传递闭包,或采用“阈值二值化”方法(\(H = T + I\) 后按 \(\lambda\) 二值化)。
  • 计算可达集、前因集、交集
    • 可达集 \(R(i)\):第 \(i\) 行中值为 1 的因素集合。
    • 前因集 \(Q(i)\):第 \(i\) 列中值为 1 的因素集合。
    • 交集 \(C(i) = R(i) \cap Q(i)\)
  • 层级划分:找出满足 \(R(i) = C(i)\) 的因素作为当前层级,移除后重复,直至所有因素分层完毕。
  • 计算骨架矩阵:删除可达矩阵中的传递冗余关系,保留必要的直接关系。

7. 结果可视化

  • DEMATEL 图:以中心度为横坐标、原因度为纵坐标绘制散点图,用颜色区分因素属性。
  • ISM 层次结构图:根据骨架矩阵和层级绘制有向层次图,顶层为驱动因素,底层为结果因素。

案例分析

案例背景:某企业拟分析影响产品质量的五个因素:员工技能(F1)、设备精度(F2)、原材料质量(F3)、工艺规范(F4)、检验流程(F5)。邀请两位专家(权重相等)采用 0~4 分制打分,原始数据(直接矩阵,仅示专家1)如下:

F1 F2 F3 F4 F5
F1 0 3 2 4 1
F2 2 0 3 3 2
F3 1 2 0 1 3
F4 3 2 2 0 2
F5 2 1 4 2 0

计算过程

1. 聚合专家矩阵(等权)得到直接影响矩阵 \(A\)(略)

2. DEMATEL 计算

采用行和最大值规范化,得到规范矩阵 \(B\),再计算综合影响矩阵 \(T\)(数值已简化):

\[ T = \begin{bmatrix} 0.15 & 0.22 & 0.18 & 0.25 & 0.12 \\ 0.20 & 0.12 & 0.21 & 0.23 & 0.14 \\ 0.10 & 0.15 & 0.08 & 0.12 & 0.20 \\ 0.18 & 0.17 & 0.14 & 0.10 & 0.15 \\ 0.12 & 0.11 & 0.22 & 0.16 & 0.09 \end{bmatrix} \]

3. 计算 DEMATEL 指标

因素 影响度 D 被影响度 C 中心度 M 原因度 R 属性
F1 0.92 0.75 1.67 0.17 一般原因
F2 0.90 0.77 1.67 0.13 一般原因
F3 0.65 0.83 1.48 -0.18 一般结果
F4 0.74 0.86 1.60 -0.12 一般结果
F5 0.70 0.70 1.40 0.00 一般结果

平均中心度 \(\bar{M} \approx 1.56\),因此无核心因素。

4. 设定阈值 \(\lambda\)(采用均值+标准差)

均值 \(m_T = 0.15\),标准差 \(sd_T = 0.05\)\(\lambda = 0.20\)。二值化后得邻接矩阵 \(A'\)

\[ A' = \begin{bmatrix} 0 & 1 & 0 & 1 & 0 \\ 0 & 0 & 1 & 1 & 0 \\ 0 & 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \end{bmatrix} \]

5. ISM 可达矩阵与分层

  • 添加自反性,计算传递闭包,得到可达矩阵 \(R\)(略)。

  • 可达集、前因集计算后分层:

    • 第1层:F5(检验流程)
    • 第2层:F3(原材料质量)、F4(工艺规范)
    • 第3层:F2(设备精度)
    • 第4层:F1(员工技能)

6. 结果解读

  • DEMATEL:F1、F2 是原因因素(驱动),F3、F4、F5 是结果因素;各因素中心度相近,重要性差别不大。
  • ISM:F1(员工技能)为最底层驱动因素,F5(检验流程)为顶层结果因素,中间层为传导因素。

常见问题

Q1: 为什么要融合 DEMATEL 和 ISM?

A: DEMATEL 擅长量化因素间的影响强度和因果属性,但无法给出清晰的层次结构;ISM 擅长构建层级关系,但输入依赖二元关系(0/1),无法利用强度信息。融合方法利用 DEMATEL 的综合影响矩阵生成更合理的邻接矩阵,再通过 ISM 揭示系统层次,兼顾定量与结构化分析。

Q2: 阈值 λ 如何选择?

A: 平台支持三种方式:

  • 均值+标准差:常用方法,λ = mean(T) + sd(T),能筛选出显著影响关系。
  • 均值:λ = mean(T),较为宽松。
  • 自定义:用户根据经验设定(0~1)。建议通过敏感性分析选择。

Q3: 可达矩阵生成方法有何区别?

A:

  • 布尔运算迭代(Warshall):对邻接矩阵 \(A' + I\) 进行传递闭包,严格遵循 ISM 理论,适合无环系统。
  • 阈值二值化:直接对 \(T + I\) 按 λ 二值化,不进行传递闭包,适合 DEMATEL 结果本身已包含间接关系的情况。一般推荐布尔运算迭代。

Q4: 如何设置专家权重?

A: 每位专家的权重应反映其权威性或熟悉程度,权重和应为 1。平台提供等权重默认值,用户可手动调整。

Q5: 支持多工作表吗?

A: 支持。Excel 文件中每个工作表代表一位专家的直接影响矩阵,用户可多选工作表,并分别设置权重,系统自动加权聚合。

平台功能

DEMATEL-ISM 融合分析平台提供以下核心功能:

数据输入

  • 支持 Excel(.xlsx, .xls)格式。
  • 每个工作表为一个专家的直接影响矩阵(数值 0~4),第一行和第一列为因素名称,数据区域为方阵。

参数设置

  • DEMATEL 规范化方法:行和最大值、列和最大值、行和列和最大值。
  • 阈值 λ 计算方法:均值+标准差、均值、自定义。
  • 自定义阈值:当选择自定义时输入(0~1)。
  • 可达矩阵生成方法:布尔运算迭代、阈值 λ 二值化。
  • ISM 分层方法:两种等价选项。
  • 选择分析的工作表:可多选,支持多专家。
  • 专家权重设置:为每个专家分配权重(自动归一化)。
  • 对角线置为 0:强制矩阵对角线为 0。

结果展示

  • DEMATEL 因素分析结果:影响度、被影响度、中心度、原因度、权重、排序、属性分类。
  • ISM 因素分级结果:各因素层级、可达集、前因集、交集。
  • 矩阵展示:直接影响矩阵、规范影响矩阵、综合影响矩阵、验证矩阵(DEMATEL);邻接矩阵、可达矩阵、骨架矩阵(ISM)。
  • 转换参数:阈值 λ、可达矩阵生成方法。
  • 分层迭代过程:每一层迭代的详细数据。
  • 可视化:DEMATEL 散点图(中心度-原因度图)和 ISM 层次结构图,均支持丰富的绘图参数自定义。
  • AI 智能分析:基于 DeepSeek API 自动解读结果,提供关键因素识别与层级结构分析(每日限 3 次)。
  • 多格式导出:支持 Excel 和 HTML 报告下载。

使用建议

  1. 准备阶段:确定系统因素列表(建议不超过 15 个),设计问卷,明确打分尺度(0~4 分)。邀请专家独立填写。

  2. 数据收集:将每位专家的矩阵放入 Excel 的不同工作表,按模板格式填写(第一行第一列为因素名称,数据区域为数值方阵)。

  3. 参数设置

    • 选择规范化方法(推荐行和最大值法)。
    • 选择阈值 λ 计算方法(推荐均值+标准差),或根据经验自定义。
    • 选择可达矩阵生成方法(推荐布尔运算迭代)。
    • 若有多位专家,合理设置权重。
  4. 结果解读

    • 首先关注 DEMATEL 结果:识别核心原因因素(应优先管理)和核心结果因素(应作为监测指标)。
    • 其次分析 ISM 层次结构:顶层因素为最终输出或结果,底层因素为根本驱动。
    • 结合两者:原因度高的底层因素是最关键的管理杠杆点。
    • 利用 AI 分析获取更深入的解读。
  5. 迭代优化

    • 若结果不符合预期,可调整阈值 λ 或规范化方法,进行敏感性分析。
    • 对比不同专家权重组合,检验稳定性。
    • 可将融合结果作为 ANP、SD 等方法的输入。

平台界面

官方地址:https://superr.online

DEMATEL-ISM融合工具界面

平台界面包含:数据上传区、参数设置区、多工作表预览、DEMATEL 与 ISM 结果展示、双重可视化图表和 AI 分析模块


参考文献

  1. Gabus A, Fontela E. World problems, an invitation to further thought within the framework of DEMATEL[R]. Battelle Geneva Research Centre, 1972.
  2. Warfield J N. Developing interconnection matrices in structural modeling[J]. IEEE Transactions on Systems, Man, and Cybernetics, 1974, SMC-4(1): 81-87.
  3. 周德群,章玲. 集成 DEMATEL-ISM 的复杂系统因素分析[J]. 系统工程理论与实践,2011, 31(3): 481-488.
  4. 基于 DEMATEL-ISM 的供应链风险因素研究[J]. 中国管理科学,2015, 23(12): 78-85.