Analytic hierarchy process (matlab)

  This article draws on mathematical modeling Teacher Qingfeng's courseware and ideas , You can click the view link to view the video explanation of teacher Qingfeng ： Qingfeng mathematical modeling ：https://www.bilibili.com/video/BV1DW411s7wi

Just look at the following blog , Here is a supplement to this blog .

  Learning notes of Qingfeng mathematical modeling —— Analytic hierarchy process （AHP）

What if the consistency inspection fails ？ Here's the picture ：

  Program

clc,clear
disp(' Please enter the judgment matrix A')
A=input('A=');
[n,n] = size(A);
% % % % % % % % % % % % % Method 1：  Calculate the weight by arithmetic average % % % % % % % % % % % % %
Sum_A = sum(A);
SUM_A = repmat(Sum_A,n,1);
Stand_A = A ./ SUM_A;

disp(' The result of calculating the weight by arithmetic average method is ：');
disp(sum(Stand_A,2)./n)
% % % % % % % % % % % % % Method 2：  Calculate the weight by geometric average method % % % % % % % % % % % % %
Prduct_A = prod(A,2);
Prduct_n_A = Prduct_A .^ (1/n);
disp(' The result of weight calculation by geometric average method is ：');
disp(Prduct_n_A ./ sum(Prduct_n_A))
% % % % % % % % % % % % % Method 3：  Calculate the weight by eigenvalue method % % % % % % % % % % % % %
[V,D] = eig(A);
Max_eig = max(max(D));
[r,c]=find(D == Max_eig , 1);
disp(' The result of calculating the weight by eigenvalue method is ：');
disp( V(:,c) ./ sum(V(:,c)) )
% % % % % % % % % % % % % Here is how to calculate the consistency ratio CR Link % % % % % % % % % % % % %
CI = (Max_eig - n) / (n-1);
RI=[0 0.0001 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];  % Watch out! , there RI Most support  n = 15
%  here n=2 when , It must be a uniform matrix , therefore CI = 0, In order to avoid denominator being 0, Change the second element here to be very close to 0 Positive number of 
CR=CI/RI(n);
disp(' Consistency indicators CI=');disp(CI);
disp(' Consistency ratio CR=');disp(CR);
if CR<0.10
    disp(' because CR<0.10, So the judgment matrix A The consistency of is acceptable !');
else
    disp(' Be careful ：CR >= 0.10, Therefore, the judgment matrix A Need to change !');
end

The variables inside n Is the number of indicators , In this case n by 5 .

In practice , Just input the judgment matrix into the program . In this case , As follows 6 A judgment matrix ：