MATLAB——DEMATEL代码
A%直接影响矩阵(直接右侧工作区新建一个A,然后把你的打分数据粘进去就行)
[m,n]=size(A)
for i=1:m
a(i)=sum(A(i,:))
end
M=max(a)
b=A/M%归一化
%%
t=b/(eye(m)-b)%综合影响矩阵
for i=1:m
for j=1:m
D(i)=sum(t(i,:))
R(j)=sum(t(:,j))
end
end
f=R+D %中心度
r=R-D %原因度
%%图不太好看
figure
scatter(f,r)
title('f-d')
xlabel('f')
ylabel('d')