灰色DEMATEL

方法概述

灰色DEMATEL是一种将灰色系统理论与决策实验室分析法（DEMATEL）相结合的多准则决策方法。在现实决策中，专家评价往往存在信息不完全、认知不确定等问题，传统DEMATEL使用精确数值难以反映这种模糊性。灰色DEMATEL 通过引入区间灰数（上界、下界）表示专家评价，并采用灰数清晰化算法将区间数转换为精确值，再执行标准 DEMATEL 分析，从而更合理地处理评价信息的不确定性。

该方法的核心理念：

  • 每位专家对因素间直接影响程度的判断用区间灰数表示（例如“影响程度在 2~4 之间”），分别记录上界矩阵和下界矩阵。
  • 支持多专家评价，可通过加权平均聚合各专家的上下界矩阵。
  • 采用基于区间数的清晰化公式（考虑每列最大上界），将聚合后的灰数矩阵转换为单一精确直接影响矩阵。
  • 支持两种分析模式：
    • 聚合模式：先对聚合后的灰数矩阵进行清晰化，得到一个精确矩阵，再执行 DEMATEL。
    • 区间模式：分别对聚合后的上界矩阵和下界矩阵执行 DEMATEL，得到两套结果（上界结果和下界结果），反映不确定性的影响范围。
  • 提供六种规范化方法，包括传统的行和最大值、列和最大值、行和列和取最大/最小，以及基于几何平均和算术平均的“弦”方法。
  • 计算综合影响矩阵、影响度、被影响度、中心度、原因度，并绘制中心度-原因度散点图。

该方法适用于专家评价具有明显不确定性、需要考察评价区间对结果影响的问题。

计算步骤

1. 构建灰数直接关系矩阵

设有 \(n\) 个因素，\(K\) 位专家。每位专家 \(k\) 给出一个 \(n \times n\) 的上界矩阵 \(U^{(k)} = [u_{ij}^{(k)}]\) 和一个下界矩阵 \(L^{(k)} = [l_{ij}^{(k)}]\)，其中 \(0 \le l_{ij}^{(k)} \le u_{ij}^{(k)}\)，且 \(u_{ii}^{(k)} = l_{ii}^{(k)} = 0\)。通常采用 0~4 分制或 0~10 分制，专家给出一个区间范围。

2. 聚合多专家灰数矩阵

设专家 \(k\) 的权重为 \(w_k\)（\(\sum_{k=1}^{K} w_k = 1\)），则聚合后的上界矩阵 \(\bar{U}\) 和下界矩阵 \(\bar{L}\) 的元素分别为：

\[ \bar{u}_{ij} = \sum_{k=1}^{K} w_k \cdot u_{ij}^{(k)}, \quad \bar{l}_{ij} = \sum_{k=1}^{K} w_k \cdot l_{ij}^{(k)} \]

3. 灰数清晰化（聚合模式）

对于每个非对角元素 \((\bar{l}_{ij}, \bar{u}_{ij})\)，计算清晰化值 \(s_{ij}\)。记 \(M_j = \max_{i} \bar{u}_{ij}\)（第 \(j\) 列上界的最大值），清晰化公式为：

\[ s_{ij} = \frac{\bar{l}_{ij} \cdot M_j - \bar{l}_{ij}^2 + \bar{u}_{ij}^2}{M_j - \bar{l}_{ij} + \bar{u}_{ij}} \]

若分母为零，则取 \(s_{ij} = (\bar{l}_{ij} + \bar{u}_{ij})/2\)。对角线元素 \(s_{ii}=0\)。最终得到清晰直接影响矩阵 \(S = [s_{ij}]\)。

说明：该公式源自灰色系统理论中的白化权函数，能够将区间数转化为一个“最可能”的清晰值，同时考虑了列上的最大可能影响。

4. 规范化直接影响矩阵

设清晰矩阵 \(S\) 或上界矩阵 \(\bar{U}\) / 下界矩阵 \(\bar{L}\)（取决于模式）。选择规范化方法，计算规范化因子 \(s\)：

方法名称 公式（\(r_i = \sum_j a_{ij}\)，\(c_j = \sum_i a_{ij}\)）
行和最大值 \(s = \max_i r_i\)
列和最大值 \(s = \max_j c_j\)
行和列和取最大 \(s = \max(\max_i r_i, \max_j c_j)\)
行和列和取最小 \(s = \min(\max_i r_i, \max_j c_j)\)
取两最大值的弦（几何平均） \(s = \sqrt{(\max_i r_i) \cdot (\max_j c_j)}\)
弦的最大值（算术平均） \(s = (\max_i r_i + \max_j c_j) / 2\)

规范化矩阵 \(B = A / s\)。

5. 计算综合影响矩阵

\[ T = B + B^2 + B^3 + \cdots = B (I - B)^{-1} \]

其中 \(I\) 为单位矩阵。若 \(I-B\) 不可逆，则使用广义逆（MASS::ginv）。

6. 计算指标

  • 影响度 \(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\) 为结果因素）

以中心度均值 \(\bar{M}\) 为界，将因素分为四类： - \(R_i > 0, M_i > \bar{M}\)：核心原因 - \(R_i > 0, M_i \le \bar{M}\)：一般原因 - \(R_i < 0, M_i > \bar{M}\)：核心结果 - \(R_i < 0, M_i \le \bar{M}\)：一般结果

7. 结果可视化

  • 聚合模式：基于清晰化矩阵的 DEMATEL 结果绘制散点图。
  • 区间模式：可选择基于上界结果或下界结果绘图，观察不确定性范围。

案例分析

案例背景：某企业拟分析影响产品质量的三个因素：员工技能（F1）、设备精度（F2）、原材料质量（F3）。邀请一位专家采用区间打分（0~4 分），上下界矩阵如下：

上界矩阵 \(U\)： \[ U = \begin{bmatrix} 0 & 3 & 2 \\ 2 & 0 & 4 \\ 1 & 3 & 0 \end{bmatrix} \]

下界矩阵 \(L\)： \[ L = \begin{bmatrix} 0 & 2 & 1 \\ 1 & 0 & 3 \\ 0 & 2 & 0 \end{bmatrix} \]

（示例中仅一个专家，因此无需聚合。选择聚合模式，规范化方法采用“行和最大值”。）

计算过程

1. 灰数清晰化

计算每列上界的最大值：

  • 列1：\(\max(0,2,1)=2\)
  • 列2：\(\max(3,0,3)=3\)
  • 列3：\(\max(2,4,0)=4\)

对非对角元素应用清晰化公式：

  • \(s_{12}\)：\(l=2, u=3, M_2=3\)
    \(s = (2×3 - 2^2 + 3^2)/(3 - 2 + 3) = (6-4+9)/4 = 11/4 = 2.75\)
  • \(s_{13}\)：\(l=1, u=2, M_3=4\)
    \(s = (1×4 - 1 + 4)/(4 - 1 + 2) = (4-1+4)/5 = 7/5 = 1.4\)
  • \(s_{21}\)：\(l=1, u=2, M_1=2\)
    \(s = (1×2 - 1 + 4)/(2 - 1 + 2) = (2-1+4)/3 = 5/3 ≈ 1.667\)
  • \(s_{23}\)：\(l=3, u=4, M_3=4\)
    \(s = (3×4 - 9 + 16)/(4 - 3 + 4) = (12-9+16)/5 = 19/5 = 3.8\)
  • \(s_{31}\)：\(l=0, u=1, M_1=2\)
    \(s = (0×2 - 0 + 1)/(2 - 0 + 1) = 1/3 ≈ 0.333\)
  • \(s_{32}\)：\(l=2, u=3, M_2=3\)
    \(s = (2×3 - 4 + 9)/(3 - 2 + 3) = (6-4+9)/4 = 11/4 = 2.75\)

清晰矩阵 \(S\)： \[ S = \begin{bmatrix} 0 & 2.75 & 1.4 \\ 1.667 & 0 & 3.8 \\ 0.333 & 2.75 & 0 \end{bmatrix} \]

2. 规范化（行和最大值法）

计算行和：\([4.15, 5.467, 3.083]\)，最大值 \(s=5.467\)，规范化矩阵 \(B = S / 5.467\)：

\[ B = \begin{bmatrix} 0 & 0.503 & 0.256 \\ 0.305 & 0 & 0.695 \\ 0.061 & 0.503 & 0 \end{bmatrix} \]

3. 综合影响矩阵 \(T = B(I-B)^{-1}\)

经计算（过程略），得 \(T\)：

\[ T = \begin{bmatrix} 0.210 & 0.604 & 0.423 \\ 0.388 & 0.275 & 0.810 \\ 0.127 & 0.557 & 0.217 \end{bmatrix} \]

4. 计算指标

  • 影响度 \(D\)：\([1.237, 1.473, 0.901]\)
  • 被影响度 \(C\)：\([0.725, 1.436, 1.450]\)
  • 中心度 \(M\)：\([1.962, 2.909, 2.351]\)
  • 原因度 \(R\)：\([0.512, 0.037, -0.549]\)
  • 平均中心度 \(\bar{M} = 2.407\)

因素属性：

  • F1：\(R>0, M<2.407\) → 一般原因
  • F2：\(R>0, M>2.407\) → 核心原因
  • F3：\(R<0, M<2.407\) → 一般结果

结论：设备精度（F2）为核心原因因素，应优先改善；员工技能（F1）为一般原因；原材料质量（F3）为结果因素，受前两者影响。

常见问题

Q1: 什么是灰色 DEMATEL？与传统 DEMATEL 有何区别？

A: 传统 DEMATEL 要求专家给出精确的直接影响评分，而灰色 DEMATEL 允许专家用区间灰数（上界、下界）表达判断的不确定性。它通过灰数清晰化算法将区间转换为精确值，再执行 DEMATEL，更符合实际决策中的信息不完全情况。

Q2: 如何理解灰数清晰化公式？

A: 该公式基于灰色系统理论的白化权函数，将区间数 \([l, u]\) 转化为一个“最可能”的清晰值，同时考虑该列的最大上界 \(M_j\)。其数学形式保证了清晰值介于 \(l\) 和 \(u\) 之间，且当 \(l=u\) 时退化为该数值。公式推导源自灰色关联分析中的“最大熵”原理。

Q3: 两种分析模式（聚合模式 vs 区间模式）如何选择？

A:

  • 聚合模式：先清晰化后分析，得到一个确定的结果。适用于需要唯一决策结论的场景。
  • 区间模式：分别对上界和下界矩阵进行 DEMATEL，得到两个极端结果，反映不确定性对结论的影响范围。适用于需要评估风险或敏感性的场景。

Q4: 六种规范化方法有何区别？

A: 前三种是 DEMATEL 中常用的（行和最大、列和最大、行和列和最大）。后三种是扩展方法：

  • 行和列和取最小：用两个最大值的较小者，可能得到较小的规范化因子，使矩阵元素相对较大。
  • 几何平均（弦） 和 算术平均（弦的最大值）：取两个最大值的平均值，是一种折中策略。几何平均对极端值更稳健。

Q5: 如何处理多专家数据？

A: 平台要求每位专家的数据包含成对的上下界工作表（例如“专家1_Upper”、“专家1_Lower”）。用户可勾选多个专家对，并设置权重，系统会自动加权平均聚合上下界矩阵。

Q6: 矩阵对角线是否必须为 0？

A: 是。DEMATEL 要求因素不自影响，因此对角线应始终为 0。平台提供“对角线置为 0”选项，可自动修正。

平台功能

灰色 DEMATEL 分析平台提供以下核心功能：

数据输入

  • 支持 Excel（.xlsx, .xls）格式。
  • 每个专家需提供两个工作表：上界矩阵（后缀 _Upper）和下界矩阵（后缀 _Lower），例如“Exp1_Upper”、“Exp1_Lower”。
  • 每个工作表的第一行和第一列为因素名称，数据区域为方阵。
  • 自动识别成对的工作表，支持多专家。

参数设置

  • 规范化方法：六种可选（行和最大值、列和最大值、行和列和取最大、行和列和取最小、几何平均、算术平均）。
  • 分析模式：聚合模式（先清晰化后分析）、区间模式（分别分析上下界）。
  • 选择分析的因素对：从检测到的专家对中多选。
  • 专家权重：为每个选中的专家分配权重（自动归一化）。
  • 对角线置为 0：强制矩阵对角线为 0。
  • 小数位数：控制输出精度（默认 4 位）。

结果展示

  • 聚合模式：
    • 因素分析结果（影响度、被影响度、中心度、原因度、权重、排序、属性）
    • 直接影响矩阵（清晰化后）
    • 规范影响矩阵
    • 综合影响矩阵
    • 统计分析（规范化因子、平均中心度、矩阵计算方法）
  • 区间模式：
    • 分别展示上界结果和下界结果，每套结果包含上述所有表格。
  • 可视化：
    • 中心度-原因度散点图（可切换为影响度-被影响度图）
    • 区间模式下可选择基于上界或下界绘图
    • 丰富的绘图参数（点状/块状、颜色、字体大小、图例等）
  • AI 智能分析：基于 DeepSeek API 自动解读结果，提供决策建议（每日限 3 次）。
  • 多格式导出：支持 Excel 和 HTML 报告下载。

使用建议

  1. 准备阶段：确定系统因素列表（建议不超过 15 个）。设计问卷，要求专家给出区间打分（例如“该因素对另一因素的影响程度在 2~4 之间”）。明确评分尺度（如 0~4 分或 0~10 分）。

  2. 数据收集：将每位专家的上界矩阵和下界矩阵分别放入 Excel 的不同工作表，命名规则为“专家名_Upper”和“专家名_Lower”。确保矩阵为方阵，对角线为 0，首行首列为因素名称。

  3. 参数设置：

    • 选择规范化方法（推荐行和最大值或几何平均）。
    • 选择分析模式：若需要唯一结论选聚合模式；若需考察不确定性范围选区间模式。
    • 若有多位专家，合理设置权重（如根据专家资历或熟悉程度）。
    • 调整绘图参数以获得清晰图形。

  4. 结果解读：

    • 聚合模式：直接查看核心原因和核心结果因素，作为管理重点。
    • 区间模式：比较上界和下界结果，观察因素属性的稳定性。若某因素在上界和下界中属性一致，则结论稳健；若差异大，则需谨慎处理或补充信息。
    • 利用 AI 分析获取专业解读。

  5. 迭代优化：

    • 若结果不符合预期，可尝试不同规范化方法或调整专家权重。
    • 对区间模式，可观察中心度和原因度的变化范围，评估不确定性对排序的影响。
    • 可将灰色 DEMATEL 结果与经典 DEMATEL（使用精确打分）对比，检验灰数处理的必要性。

平台界面

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

灰色DEMATEL工具界面

平台界面包含：数据上传区、参数设置区、多专家对预览、分析结果展示和AI分析模块

参考文献：

  1. 邓聚龙. 灰色系统理论教程[M]. 华中理工大学出版社，1990.
  2. 刘思峰，谢乃明. 灰色系统理论及其应用[M]. 科学出版社，2013.
  3. 灰色 DEMATEL 方法在供应链风险识别中的应用[J]. 系统工程理论与实践，2015, 35(9): 2287-2294.