MATLAB——DEMATEL代碼
A%直接影響矩陣(直接右側(cè)工作區(qū)新建一個(gè)A,然后把你的打分?jǐn)?shù)據(jù)粘進(jìn)去就行)
[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')
標(biāo)簽: