層次分析法程序代碼

1、層次分析法1)建立層次結(jié)構(gòu)模型:決策目標(biāo)準(zhǔn)則C1準(zhǔn)則C2準(zhǔn)則C3準(zhǔn)則C4P1P2P3P4P5P6(2)構(gòu)造判斷矩陣判斷矩陣A=(aj )應(yīng)為正互反矩陣，而且 3ij的判斷如下(19尺度法)：標(biāo)度aaiJ含義135792, 4, 6, 8Ci與Ci的影響相同Ci比Cj的影響稍強(qiáng)Ci比Cj的影響強(qiáng)Ci比Cj的影響明顯的強(qiáng)e比e的影響絕對的強(qiáng)c與c的影響之比在上述兩個(gè)相鄰的等 級之間1 12，9Ci與Cj的影響之比為上面aiJ的復(fù)反數(shù)(3 )單層排序及一致性檢驗(yàn)1、單層排序求解判斷矩陣 A的最大特征值max，再由最大特征值求出對應(yīng)的特征向量- A 二，max，',并將標(biāo)準(zhǔn)化，即為同一層相對于

2、上一層某一因素的權(quán)重，根據(jù)此 權(quán)重的大小，便可確定該層因素的排序。2、一致性檢驗(yàn)取一致性指標(biāo)CI二-mx n，（ n為A的階數(shù)）n 1取隨機(jī)性指標(biāo) RI如下:n123456789RI000.580.91.121.241.321.411.45令CR二里，若CR ： 0.1，則認(rèn)為A具有一致性。RI否則，需要對 A進(jìn)行調(diào)整，直到具有滿意的一致性為止。（4）層次總排序及一致性檢驗(yàn)假定準(zhǔn)則層Ci,C2,ll（,Cn排序完成，其權(quán)重分別為ai,a2,IH,an，方案層P包含m個(gè)方 案：R,PJII,Pm。其相對于上一層的 Cj（j =1,2,|樸，n ）對方案層P中的m個(gè)方案進(jìn)行單 層排序，其排序權(quán)重記

3、為4）山2）,|,% j=1,2,|, n ,則方案層P中第i個(gè)方案Pi的總n排序權(quán)重為Zj majbj，見卜表：層次CC1C2CnP層總排序權(quán)重層次pa1a2anRbnb12b1nnZ ajdjj#P2b22b2nn為 ajb2jj#Pnbm1bm2bmnn壬 ajbmjj#從而確定P層的排序。例：純文本文件txt3.txt中的數(shù)據(jù)格式如下:1 1 1 4 1 1/21 1 2 4 1 1/21 1/2 1 5 3 1/21/4 1/4 1/5 1 1/3 1/31 1 1/3 3 1 12 2 2 3 3 11 1/4 1/24 1 32 1/3 11 1/4 1/54 1 1/25 2

4、11 3 1/31/3 1 1/73 7 11 1/3 53 1 71/5 1/7 11 1 71 1 71/7 1/7 11 7 91/7 1 11/9 1 1matlab 程序：>> fid=fopen('txt3.txt','r');n1=6;n2=3;a=;for i=1:n1tmp=str2num(fgetl(fid);a=a;tmp; % 讀準(zhǔn)則層判斷矩陣endfor i=1:n1str1=char('b',int2str(i),'=;'); str2=char('b',int2str(i)

5、,'=b',int2str(i),'tmp;'); eval(str1);for j=1:n2tmp=str2num(fgetl(fid);eval(str2); % 讀方案層的判斷矩陣end致性指標(biāo)endri=0,0,0.58,0.90,1.12,1.24,1.32,1.41,1.45; % x,y=eig(a);lamda=max(diag(y);num=find(diag(y)=lamda);w0=x(:,num)/sum(x(:,num);cr0=(lamda-n1)/(n1-1)/ri(n1)for i=1:n1x,y=eig(eval(char('b',int2str(i);lamda=max(diag(y); num=find(diag(y)=lamda); w1(:,i)=x