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 报告下载。
使用建议
准备阶段:确定系统因素列表(建议不超过 15 个),设计问卷,明确打分尺度(0~4 分)。邀请专家独立填写。
数据收集:将每位专家的矩阵放入 Excel 的不同工作表,按模板格式填写(第一行第一列为因素名称,数据区域为数值方阵)。
参数设置:
- 选择规范化方法(推荐行和最大值法)。
- 选择阈值 λ 计算方法(推荐均值+标准差),或根据经验自定义。
- 选择可达矩阵生成方法(推荐布尔运算迭代)。
- 若有多位专家,合理设置权重。
结果解读:
- 首先关注 DEMATEL 结果:识别核心原因因素(应优先管理)和核心结果因素(应作为监测指标)。
- 其次分析 ISM 层次结构:顶层因素为最终输出或结果,底层因素为根本驱动。
- 结合两者:原因度高的底层因素是最关键的管理杠杆点。
- 利用 AI 分析获取更深入的解读。
迭代优化:
- 若结果不符合预期,可调整阈值 λ 或规范化方法,进行敏感性分析。
- 对比不同专家权重组合,检验稳定性。
- 可将融合结果作为 ANP、SD 等方法的输入。
平台界面

平台界面包含:数据上传区、参数设置区、多工作表预览、DEMATEL 与 ISM 结果展示、双重可视化图表和 AI 分析模块
参考文献:
- Gabus A, Fontela E. World problems, an invitation to further thought within the framework of DEMATEL[R]. Battelle Geneva Research Centre, 1972.
- Warfield J N. Developing interconnection matrices in structural modeling[J]. IEEE Transactions on Systems, Man, and Cybernetics, 1974, SMC-4(1): 81-87.
- 周德群,章玲. 集成 DEMATEL-ISM 的复杂系统因素分析[J]. 系统工程理论与实践,2011, 31(3): 481-488.
- 基于 DEMATEL-ISM 的供应链风险因素研究[J]. 中国管理科学,2015, 23(12): 78-85.