层次分析法matlab代码实现

disp('请输入准则层判断矩阵A(n阶)');
A=input('A=');
[n,n]=size(A); % 得到矩阵的阶数
[V,D]=eig(A); % 求得特征向量和特征值
% 求出最大特征值和它所对应的特征向量
tempNum=D(1,1); % 特征值的初值
pos=1; % 标记的第一个数
for h=1:n
    if D(h,h)>tempNum
        tempNum=D(h,h);
        pos=h; % 最大特征值的位置
    end
end
w=abs(V(:,pos)); % 找最大特征值对应的特征向量
w=w/sum(w); % 归一化处理
t=D(pos,pos); % t指的是最大的特征值
disp('准则层特征向量w=');disp(w);disp('准则层最大特征根t=');disp(t);
% 一致性检验
CI=(t-n)/(n-1);RI=[0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59 1.60 1.61 1.615 1.62 1.63];
CR=CI/RI(n);
if CR<0.10
    disp('此矩阵的一致性可以接受!');
    disp('CI=');disp(CI);
    disp('CR=');disp(CR);
else disp('此矩阵的一致性验证失败，请重新进行评分!');
end

以下是对每一句MATLAB代码的详细解释：

matlabCopy code

disp('请输入准则层判断矩阵A(n阶)'); A = input('A=');

这段代码首先通过 disp 函数输出提示信息，要求用户输入准则层判断矩阵 A 的阶数 n。然后使用 input 函数接受用户输入的矩阵 A。

matlabCopy code

[n, n] = size(A); % 得到矩阵的阶数

获取矩阵 A 的阶数 n，因为输入的是正方形矩阵，所以两个维度的大小是相等的。

matlabCopy code

[V, D] = eig(A); % 求得特征向量和特征值

使用 eig 函数计算矩阵 A 的特征向量矩阵 V 和特征值矩阵 D。

matlabCopy code

% 求出最大特征值和它所对应的特征向量 tempNum = D(1, 1); % 特征值的初值 pos = 1; % 标记的第一个数 for h = 1:n if D(h, h) > tempNum tempNum = D(h, h); pos = h; % 最大特征值的位置 end end

在这部分代码中，通过循环找到特征值矩阵 D 中最大的特征值及其位置。

matlabCopy code

w = abs(V(:, pos)); % 找最大特征值对应的特征向量 w = w / sum(w); % 归一化处理

找到最大特征值对应的特征向量 w 并进行归一化处理，使其元素之和等于1。

matlabCopy code

t = D(pos, pos); % t指的是最大的特征值 disp('准则层特征向量w='); disp(w); disp('准则层最大特征根t='); disp(t);

输出最大特征值对应的特征向量 w 和最大特征值 t。

matlabCopy code

% 一致性检验 CI = (t - n) / (n - 1); RI = [0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59 1.60 1.61 1.615 1.62 1.63]; CR = CI / RI(n);

进行一致性检验，计算一致性指标 CI 和一致性比率 CR。 RI 是随机一致性指标表，其值根据矩阵的阶数 n 取定。

matlabCopy code

if CR < 0.10 disp('此矩阵的一致性可以接受!'); disp('CI='); disp(CI); disp('CR='); disp(CR); else disp('此矩阵的一致性验证失败，请重新进行评分!'); end

根据一致性比率 CR 的值进行判断，如果小于0.10，则输出一致性可以接受的信息，否则输出一致性验证失败的信息。