可能度矩阵MATLAB代码演示

function value = p_ij(A, B)

    s = (A(2) - A(1)) * (B(2) - B(1));

    if ((A(1) >= B(1)) && (A(2) >= B(2)) && (A(1) <= B(2)))
        s2 = (B(2) - A(1)) * (B(2) - A(1)) / 2;
    end

    if ((A(1) <= B(1)) && (A(2) <= B(2)) && (A(2) >= B(1)))
        s2 = s - (A(2) - B(1)) * (A(2) - B(1)) / 2;
    end

    if (A(1) >= B(2))
        s2 = 0;
    end

    if (A(2) <= B(1))
        s2 = s;
    end

    if ((A(1) >= B(1)) && (A(2) <= B(2)))
        s2 = (B(2) - A(2)) * (A(2) - A(1)) + 1/2 * (A(2) - A(1)) * (A(2) - A(1));
    end

    if ((A(1) <= B(1)) && (A(2) >= B(2)))
        s2 = (B(1) - A(1)) * (B(2) - B(1)) + 1/2 * (B(2) - B(1)) * (B(2) - B(1));
    end

    value = (s - s2) / s;

美团函数调用

wp1 = [1.0413, 1.3711];
wp2 = [0.9969, 1.1858];
wp3 = [0.777, 0.9234];
p_ij(wp1, wp1);
p_ij(wp1, wp2);
p_ij(wp1, wp3);
p_ij(wp2, wp1);
p_ij(wp2, wp2);
p_ij(wp2, wp3);
p_ij(wp3, wp1);
p_ij(wp3, wp2);
p_ij(wp3, wp3);
P = [p_ij(wp1, wp1), p_ij(wp1, wp2), p_ij(wp1, wp3);
   p_ij(wp2, wp1), p_ij(wp2, wp2), p_ij(wp2, wp3);
   p_ij(wp3, wp1), p_ij(wp3, wp2), p_ij(wp3, wp3); ];

美团结果对不上

---

参考文献[15]，参数调用

w1 = [0.2167, 0.2833];
w2 = [0.3083, 0.3167];
w3 = [0.225, 0.2583];
w4 = [0.175, 0.2167];
p_ij(w1, w1);
p_ij(w1, w2);
p_ij(w1, w3);
p_ij(w1, w4);
p_ij(w2, w1);
p_ij(w2, w2);
p_ij(w2, w3);
p_ij(w2, w4);
p_ij(w3, w1);
p_ij(w3, w2);
p_ij(w3, w3);
p_ij(w3, w4);
p_ij(w4, w1);
p_ij(w4, w2);
p_ij(w4, w3);
p_ij(w4, w4);
P = [p_ij(w1, w1), p_ij(w1, w2), p_ij(w1, w3), p_ij(w1, w4);
   p_ij(w2, w1), p_ij(w2, w2), p_ij(w2, w3), p_ij(w2, w4);
   p_ij(w3, w1), p_ij(w3, w2), p_ij(w3, w3), p_ij(w3, w4);
   p_ij(w4, w1), p_ij(w4, w2), p_ij(w4, w3), p_ij(w4, w4); ];

这个结果对的上

---

本文章使用limfx的vscode插件快速发布