層次分析法matlab源程序6頁

1、層次分析法matlab源程序disp('請輸入判斷矩陣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(:,1)=x(:,1);x(:,2)=A*y(:,1);m(2)=max(x(:,2);y(:,2)=x(:,2)/m(2);p=0.0001;i=2;k=abs(m(2)-m(1);while  k>p  i=i+1;  x(:,i)=A*y(:,i-1);  m(i)=ma

2、x(x(:,i);  y(:,i)=x(:,i)/m(i);  k=abs(m(i)-m(i-1);enda=sum(y(:,i);w=y(:,i)/a;t=m(i);disp(w);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;CR=CI/RI(n);if CR<0.10    d

3、isp('此矩陣的一致性可以接受!');    disp('CI=');disp(CI);    disp('CR=');disp(CR);endfunction AHPInit1(x,y)%層次分析的初始化%默認(rèn)只有兩層 x為準(zhǔn)則數(shù)，y為方案數(shù)%CToT為準(zhǔn)則對目標(biāo)生成的比較陣%EigOfCri為準(zhǔn)則層的特征向量%EigOfOpt為選項層的特征向量EigOfCri=zeros(x,1);%準(zhǔn)則層的特征向量EigOfOpt=zeros(y,x);dim=x;%維度RI=0 0 0.58

4、0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51;%RI標(biāo)準(zhǔn)%生成成對比較陣for i=1:dim        CToT(i,:)=input('請輸入數(shù)據(jù):');endCToT %輸出pause,tempmatrix=zeros(x+1);tempmatrix=AHP1(dim,CToT);EigOfCri=tempmatrix(1:x);ci1=tempmatrix(1+x);EigOfCrici1pause,matrix=cell(x);%元胞數(shù)組ci=zeros(1,x

5、);dim=y;for k=1:x    matrixk=zeros(dim,dim);%生成成對比較陣for i=1:dim        matrixk(i,:)=input('請輸入數(shù)據(jù):');end%判斷該比較陣是不是一致陣tempmatrix=zeros(y+1);tempmatrix=AHP1(dim,matrixk);EigOfOpt(:,k)=tempmatrix(1:y);ci(k)=tempmatrix(y+1);EigOfOpt(:,k)ci(k)pause

6、,end%下面進(jìn)行組合一致性檢查RI=0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51;CR=ci1/RI(x)+ci*EigOfCri/RI(y);CRif CR>0.1    disp('組合一致性不通過，請重新評分')    returnend%下面根據(jù)比較陣的結(jié)果進(jìn)行組合result=EigOfOpt*EigOfCri;resultfunction f=AHP1(dim,CmpMatrix)RI=0 0 0.58 0.90 1.12 1.24 1.32 1.

7、41 1.45 1.49 1.51;%判斷該比較陣是不是一致陣%判斷該比較陣是不是一致陣V,D=eig(CmpMatrix);%求得特征向量和特征值%求出最大特征值和它所對應(yīng)的特征向量tempNum=D(1,1);pos=1;for h=1:dim    if D(h,h)>tempNum        tempNum=D(h,h);        pos=h;    endend    eigVector=V(:,pos);maxeig=D(pos,pos);maxeigdimCI=(maxeig-dim)/(dim-1);CR=CI/RI(dim);if CR>0.1    disp('準(zhǔn)則對目標(biāo)影響度評分生成的矩陣不是一致陣，請重新評分')    returnendCI%歸一化sum=0;for h=1:dim &#