yeqing.m

D-S证据理论数据融合

% 基于权重系数与冲突概率重新分配的证据合成方法,叶 清,系统工程与电子技术。
clear;clc

m1A = 0.98; m1B=0.01; m1C=0.01;
m2A = 0   ; m2B=0.01; m2C=0.99;
m3A = 0.9;  m3B=0;    m3C=0.1;

%构成冲突向量
k12 = 1 - m1A*m2A - m1B*m2B - m1C*m2C;
k13 = 1 - m1A*m3A - m1B*m3B - m1C*m3C;
k23 = 1 - m3A*m2A - m3B*m2B - m3C*m2C;

K1 = [k12,k13];
K2 = [k12,k23];
K3 = [k13,k23];

%冲突向量的归一化处理
K1 = K1/sum(K1);
K2 = K2/sum(K2);
K3 = K3/sum(K3);

%冲突向量的熵
H1 = K1(1)*log(K1(1))+K1(2)*log(K1(2));
H2 = K2(1)*log(K2(1))+K2(2)*log(K2(2));
H3 = K3(1)*log(K3(1))+K3(2)*log(K3(2));

%取熵的倒数
H1 = -1/H1;
H2 = -1/H2;
H3 = -1/H3;

%计算证据的权重
w1 = H1/(H1+H2+H3);
w2 = H2/(H1+H2+H3);
w3 = H3/(H1+H2+H3);
W = [w1,w2,w3];
%折扣率
alpha =  W/max(W);

%调整后的基本概率分配(加权)
m1A = W(1)*m1A; m1B = W(1)*m1B; m1C = W(1)*m1C;
m2A = W(2)*m2A; m2B = W(2)*m2B; m2C = W(2)*m2C;
m3A = W(3)*m3A; m3B = W(3)*m3B; m3C = W(3)*m3C;

% %调整后的冲突
% k = 1-m1A*m2A*m3A-m1B*m2B*m3B-m1C*m2C*m3C;
% PA = m1A*m2A*m3A;
% PB = m1B*m2B*m3B;
% PC = m1C*m2C*m3C;
% QA = (m1A+m2A+m3A)/3;
% QB = (m1B+m2B+m3B)/3;
% QC = (m1C+m2C+m3C)/3;
% 
% mA = PA + k*QA;
% mB = PB + k*QB;
% mC = PC + k*QC;

%------------------------------直接加权
mA = m1A + m2A + m3A;
mB = m1B + m2B + m3B;
mC = m1C + m2C + m3C;
mA = mA/(mA+mB+mC);
mB = mB/(mA+mB+mC);
mC = mC/(mA+mB+mC);