%層次分析法:plan 最終方案；rule1 準(zhǔn)則層比較矩陣；rule2 方案層比較矩陣(為cell形式) b為各方案對(duì)目標(biāo)權(quán)重，plan為最優(yōu)方案
%erros為rule2中出錯(cuò)的比較矩陣
function [b,plan,erros]=flatanalyse(rule1,rule2)
b=[];
RI=[0.00001 0.00001 0.58 0.9 1.12 1.24 1.32 1.41 1.45 1.49 1.51 1.54 1.56 1.58 1.59 1.605];
%一致性檢驗(yàn)
[x1,a1]=eigs(rule1);%x1為特征向量，a1特征根
c1=imag(a1(1));
b1=abs(x1(:,1));
b1=b1./sum(b1);
[m,n]=size(rule1);
if (a1(1)-n)/(n-1)/RI(m)>0.1| abs(c1)>0.001
    tt=1;
    sprintf('you need to modify your rule1')
end 
nn=size(rule2);
erros=zeros(nn);
tt=0;
for i=1:nn
    a=rule2{i};
    [m2,n2]=size(a);
    [x2,a2]=eigs(a);  
    c2=imag(a2(1));
    b2(:,i)=abs(x2(:,1));
    b2(:,i)=b2(:,i)./sum(b2(:,i));
    if (a2(1)-n2)/(n2-1)/RI(n2)>0.1| abs(c2)>0.001
       tt=1;
       erros(i)=1
       sprintf('you need to modify your rule2')
    end

end
if tt~=1
     for i=1:nn
         b(:,i)=b2(:,i).*b1;
     end
     p=sum(b,2);
     plan=find(p==max(p));
 end