function [lamata,w] = zhenghufajuzhen_hefa(a)
%求解正互反矩阵的最大特征值和特征向量的和法算法
n = size(a,2);
for i = 1:n
    for j = 1:n
        aa(i,j) = a(i,j)/sum(a(:,j));
    end
end
aaa = [sum(aa(1,:)),sum(aa(2,:)),sum(aa(3,:))];
for i = 1:n
    w(i) = aaa(i)/sum(aaa);
end
aw = a*w';
lamata = sum(aw'./w)/n;
w = w';

 

运行：

a =

    1.0000    2.0000    6.0000
    0.5000    1.0000    4.0000
    0.1667    0.2500    1.0000

>> [lamata,w] = zhenghufanjuzhen_hefa(a)

lamata =

    3.0092

w =

    0.5869
    0.3238
    0.0893