层次分析法matlab程序.doc

文档介绍：dispC请输入判断矩阵A(n阶门;
A=input('A=');
[n,n] = size(A);
x= ones( n,100);
y= ones( n,100);
m = zeros(1,100);
m(1 )= max(x(:, 1));
y(：」)= x(：,i)；
x(:,2) = A*y(:,1);
m(2)= max(x(: ,2));
y(:,2) = x(:,2)/m(2);
p= ;i=2;k= abs(m(2)-m(1)); while k> p
i=i+1;
x(:,i) = A*y(:,i-1);
m(i)= max(x(: ,i));
y(:,i) = x(:,i)/m(i);
k=abs(m(i)-m(i-1));
end
a=sum(y(:,i));
w=y(:,i)/a;
t=m(i);
disp(w);disp(t);
%以下是一致性检验
CI = (t-n)/(n-1);RI = [0 0
];
CR= CI/RI(n);
if CR<
dispC此矩阵的一致性可以接受门;
disp(,CI = ');disp(CI);
disp(,CR=,);disp(CR);
end
function AH PI nit1 (x,y)
%层次分析的初始化
%默认只有两层x为准则数，y为方案数
% CToT为准则对目标生成的比较阵
%EigOfCri为准则层的特征向量
%EigOfOpt为选项层的特征向量
EigOfCri= zeros(x,1);%准则层的特征向量
EigOfOpt= zeros(y,x);
dim = x;% 维度
RI = [0 0 ];%RI 标准
%生成成对比较阵
for i= 1:dim
CToT(i,:)= input('if输入数据:')；
end
CToT %输出
pause, tempmatrix= zeros(x+ 1);
tempmatrix= AHP1 (dim,CToT);
EigOfCri= tempmatrix(1 :x);
ci1 = tempmatrix(1 + x);
EigOfCri
ci1
pause,
matrix= cell(x);% 元胞数组
ci= zeros(1 ,x);
dim = y;
for k= 1 :x
matrix{ k} = zeros(dim,dim);
%生成成对比较阵
for i= 1:dim
matrix{ k} (i,:) = input('请输入数据:*);
end
%判断该比较阵是不是一致阵
tempmatrix= zeros(y+ 1);
tempmatrix= AHP1 (dim,matrix{ k}):
EigOfOpt(:, k) = tempmatrix(1: y);
ci(k)= tempmatrix(y+ 1);
EigOfOpt(:,k)
c