ahpcalc.m

阿瑟玩儿玩儿温热温热我认为而我热热特特让他人味儿我说的撒发生地方的是发生地方反思

function [finalrank,rank,rc]=ahpcalc(order,altvcriteria,critvcriteria)
%AHPCALC    Perform a complete AHP calculation
%    [FINALRANK,RANK,RC]=AHPCALC(ORDER,ALTVCRITERIA,CRITVCRITERIA)
%
%    For N alternatives and M criteria for ranking:
%
%     ORDER is a matrix with M rows and N columns.
%       The elements of the I_th row are the integers 1 to N,
%       indexing the alternatives in descending order of preference
%       with respect to criterion I.
%
%     ALTVCRITERIA is a matrix with M rows and N*(N-1)/2 columns.
%       The I_th row contains the upper triangular entries of the AHP
%       reciprocal matrix ranking the N alternatives with respect to
%       criterion I, in the descending order of preference indicated
%       in the I_th row of ORDER. 
%
%     CRITVCRITERIA may be a vector of length M*(M-1)/2 giving the
%       upper triangular entries of the AHP reciprocal matrix
%       comparing the criteria, or it may be a normalized vector of
%       length M giving the AHP vector weighting the criteria.
%
%    FINALRANK is a normalized vector of length N giving the overall
%     ranking of the alternatives.
%
%    The rows of RANK give the ranking of the alternatives with
%     respect to the individual attributes.
%    
%    RC is the vector giving the normalized weights of the
%     criteria. Note that if CRITVCRITERIA is specified as a
%     normalized vector, then RC=CRITVCRITERIA
%
% Copyright (C) 2001 Michael J. Scott
ordersize = size(order);
avcsize = size(altvcriteria);
cvcsize = size(critvcriteria);

if ordersize(1) ~= avcsize(1)
 disp 'ORDER and ALTVCRITERIA must have same number of rows.';
 disp 'Exiting.';
 return;
end
if (ordersize(2)*(ordersize(2)-1)/2) ~= avcsize(2)
 disp 'Alternative comparisons of wrong order.';
 disp 'Exiting.';
 return;
end
% If the criteria comparison is already an appropriate normalized
% vector, then use it.
if (cvcsize(2) == ordersize(1))
 if (sum(critvcriteria) <2)
  rc = critvcriteria';
 end
else   %otherwise calculate the normalized vector
 if (ordersize(1)*(ordersize(1)-1)/2) ~= cvcsize(2)
  disp 'Attribute comparisons of wrong order.';
  disp 'Exiting.';
  return;
 end

 critmat = ahpmakemat(critvcriteria);
 [rc,cc] = ahpvector(critmat);
 if cc > 0.1
  s=sprintf('   C.R.=%.3f from criteria-criteria comparisons',cc);
  disp([s]);
 end
end

for i=1:ordersize(1)
 amat = ahpmakemat(altvcriteria(i,:));
 [ra,ca] = ahpvector(amat);
 if ca > 0.1
  s= sprintf('  C.R.=%.3f from comparing alternatives on criterion %d',ca,i);
  disp([s]);
 end
 rank(i,:)=ahpreorder(order(i,:),ra);
end
finalrank = rank'*rc;