Matlab求解層次分析法程序代碼

層次分析法1）建立層次結(jié)構(gòu)模型：決策目標(biāo)準(zhǔn)則C1準(zhǔn)則C2準(zhǔn)則C3準(zhǔn)則C4P1 P2 P3 P4 P5 P62）構(gòu)造判斷矩陣標(biāo)度aj含義1C與C的影響相同i j3C比C的影響稍強(qiáng)i j5C比C的影響強(qiáng)i j7C比C的影響明顯的強(qiáng)i j9C比C的影響絕對(duì)的強(qiáng)ijC與C的影響之比在上述兩個(gè)相鄰的等2,4，6，8ij級(jí)之間112'9C與C的影響之比為上面a的復(fù)反數(shù)ijij判斷矩陣A=Cj）應(yīng)為正互反矩陣’而且寫的判斷如下A9尺度法）:（3）單層排序及一致性檢驗(yàn)1、單層排序求解判斷矩陣A的最大特征值入max，再由最大特征值求出對(duì)應(yīng)的特征向量①3=^”'

并將①標(biāo)準(zhǔn)化，即為同一層相對(duì)于上一層某一因素的權(quán)重，根據(jù)此權(quán)重的大小，便可確定該層

因素的排序。2、一致性檢驗(yàn)九一n取一致性指標(biāo)CI=max, (n為A的階數(shù))n—1取隨機(jī)性指標(biāo)RI如下:n123456789RI00令CR=，若CR<0.1，則認(rèn)為A具有一致性。RI否則，需要對(duì)A進(jìn)行調(diào)整，直到具有滿意的一致性為止。(4)層次總排序及一致性檢驗(yàn)假定準(zhǔn)則層C,C,L,C排序完成，其權(quán)重分別為a,a,L,a，方案層P包含m個(gè)方案：TOC\o"1-5"\h\z1 2n 12nP,P,L,P。其相對(duì)于上一層的C.(j=1,2,L,n)對(duì)方案層P中的m個(gè)方案進(jìn)行單層排序，12m j其排序權(quán)重記為b,b,L,b (j=1,2,L,n),則方案層P中第i個(gè)方案Pi的總排序權(quán)重為1j2j mj￡ab.，見下表:從而確定P層的排序。例：純文本文件中的數(shù)據(jù)格式如下11411/212411/211/21531/21/41/41/511/31/3111/33112233111/41/24131/3111/41/5411/2521131/31/311/77111/35171/51/711171171/71/711791/7111/911matlab程序：>>fid=fopen('','r');n1=6;n2=3;a=[];fori=1:n1tmp=str2num(fgetl(fid));a=[a;tmp];%讀準(zhǔn)則層判斷矩陣endfori=1:n1str1=char(['b',int2str(i),'=[];']);str2=char(['b',int2str(i),'=[b',int2str(i),';tmp];']);eval(str1);forj=1:n2tmp=str2num(fgetl(fid));eval(str2);%讀方案層的判斷矩陣endendri=[0,0,,,,,,,];%一致性指標(biāo)[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)fori=1:n1[x,y]=eig(eval(char(['b',int2str(i)])));lamda=max(diag(y));num=find(diag(y)==lamda);w1(:,i)=x(:,num)/sum(x(:,num));cr1(i)=(lamda-n2)/(n2-1)/ri(n2);endcr1,ts=w1*w0,cr=cr1*w0層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法層次分析法實(shí)例與步驟結(jié)合一個(gè)具體例子，說明層次分析法的基本步驟和要點(diǎn)。【案例分析】市政工程項(xiàng)目建設(shè)決策：層次分析法問題提出市政部門管理人員需要對(duì)修建一項(xiàng)市政工程項(xiàng)目進(jìn)行決策，可選擇的方案是修建通往旅游區(qū)的高速路（簡稱建高速路）或修建城區(qū)地鐵（簡稱建地鐵）。除了考慮經(jīng)濟(jì)效益外，還要考慮社會(huì)效益、環(huán)境效益等因素，即是多準(zhǔn)則決策問題，考慮運(yùn)用層次分析法解決。建立遞階層次結(jié)構(gòu)應(yīng)用AHP解決實(shí)際問題，首先明確要分析決策的問題，并把它條理化、層次化，理出遞階層次結(jié)構(gòu)。AHP要求的遞階層次結(jié)構(gòu)一般由以下三個(gè)層次組成：目標(biāo)層（最高層）：指問題的預(yù)定目標(biāo)；準(zhǔn)則層（中間層）：指影響目標(biāo)實(shí)現(xiàn)的準(zhǔn)則；措施層（最低層）：指促使目標(biāo)實(shí)現(xiàn)的措施；通過對(duì)復(fù)雜問題的分析，首先明確決策的目標(biāo)，將該目標(biāo)作為目標(biāo)層（最高層）的元素，這個(gè)目標(biāo)要求是唯一的，即目標(biāo)層只有一個(gè)元素。然后找出影響目標(biāo)實(shí)現(xiàn)的準(zhǔn)則，作為目標(biāo)層下的準(zhǔn)則層因素，在復(fù)雜問題中，影響目標(biāo)實(shí)現(xiàn)的準(zhǔn)則可能有很多，這時(shí)要詳細(xì)分析各準(zhǔn)則因素間的相互關(guān)系，即有些是主要的準(zhǔn)則，有些是隸屬于主要準(zhǔn)則的次準(zhǔn)則，然后根據(jù)這些關(guān)系將準(zhǔn)則元素分成不同的層次和組，不同層次元素間一般存在隸屬關(guān)系，即上一層元素由下一層元素構(gòu)成并對(duì)下一層元素起支配作用，同一層元素形成若干組，同組元素性質(zhì)相近，一般隸屬于同一個(gè)上一層元素（受上一層元素支配），不同組元素性質(zhì)不同，一般隸屬于不同的上一層元素。在關(guān)系復(fù)雜的遞階層次結(jié)構(gòu)中，有時(shí)組的關(guān)系不明顯，即上一層的若干元素同時(shí)對(duì)下一層的若干元素起支配作用，形成相互交叉的層次關(guān)系，但無論怎樣，上下層的隸屬關(guān)系應(yīng)該是明顯的。最后分析為了解決決策問題（實(shí)現(xiàn)決策目標(biāo)）、在上述準(zhǔn)則下，有哪些最終解決方案（措施），并將它們作為措施層因素，放在遞階層次結(jié)構(gòu)的最下面（最低層）。明確各個(gè)層次的因素及其位置，并將它們之間的關(guān)系用連線連接起來，就構(gòu)成了遞階層次結(jié)構(gòu)?！景咐治觥渴姓こ添?xiàng)目進(jìn)行決策：建立遞階層次結(jié)構(gòu)在市政工程項(xiàng)目決策問題中，市政管理人員希望通過選擇不同的市政工程項(xiàng)目，使綜合效益最高，即決策目標(biāo)是“合理建設(shè)市政工稈，使綜合效益最高”。為了實(shí)現(xiàn)這一目標(biāo)，需要考慮的主要準(zhǔn)則有三個(gè)，即經(jīng)濟(jì)效益、社會(huì)效益和環(huán)境效益。但問題絕不這么簡單。通過深入思考，決策人員認(rèn)為還必須考慮直接經(jīng)濟(jì)效益、間接經(jīng)濟(jì)效益、方便日常出行、方便假日出行、減少環(huán)境污染、改善城市面貌等因素（準(zhǔn)則），從相互關(guān)系上分析，這些因素隸屬于主要準(zhǔn)則，因此放在下一層次考慮，并且分屬于不同準(zhǔn)則。假設(shè)本問題只考慮這些準(zhǔn)則，接下來需要明確為了實(shí)現(xiàn)決策目標(biāo)、在上述準(zhǔn)則下可以有哪些方案。根據(jù)題中所述，本問題有兩個(gè)解決方案，即建高速路或建地鐵，這兩個(gè)因素作為措施層元素放在遞階層次結(jié)構(gòu)的最下層。很明顯，這兩個(gè)方案于所有準(zhǔn)則都相關(guān)。將各個(gè)層次的因素按其上下關(guān)系擺放好位置，并將它們之間的關(guān)系用連線連接起來。同時(shí)，為了方便后面的定量表示，一般從上到下用A、B、C、D。。。代表不同層次，同一層次從左到右用1、2、3、4。。。代表不同因素。這樣構(gòu)成的遞階層次結(jié)構(gòu)如下圖。目標(biāo)層A準(zhǔn)則層B準(zhǔn)則層C措施層D合理建設(shè)市政工程，使綜合效益最高(A)經(jīng)濟(jì)效益(B1)社會(huì)效益(B2)準(zhǔn)則層B準(zhǔn)則層C措施層D合理建設(shè)市政工程，使綜合效益最高(A)經(jīng)濟(jì)效益(B1)社會(huì)效益(B2)環(huán)境效益(B3)直接經(jīng)濟(jì)效益(C1)|間接帶動(dòng)效益(C2)方便日常出行(C3)方便假日出行(C4)減少環(huán)境污染(C5)改善城市面貌(C6)|圖1遞階層次結(jié)構(gòu)示意圖構(gòu)造判斷矩陣并賦值根據(jù)遞階層次結(jié)構(gòu)就能很容易地構(gòu)造判斷矩陣。構(gòu)造判斷矩陣的方法是：每一個(gè)具有向下隸屬關(guān)系的元素(被稱作準(zhǔn)則)作為判斷矩陣的第一個(gè)元素(位于左上角)，隸屬于它的各個(gè)元素依次排列在其后的第一行和第一列。重要的是填寫判斷矩陣。填寫判斷矩陣的方法有：大多采取的方法是：向填寫人(專家)反復(fù)詢問：針對(duì)判斷矩陣的準(zhǔn)則，其中兩個(gè)元素兩兩比較哪個(gè)重要，重要多少，對(duì)重要性程度按1-9賦值(重要性標(biāo)度值見下表)。表1重要性標(biāo)度含義表重要性標(biāo)度含 義1表示兩個(gè)兀素相比，具有冋等重要性3表示兩個(gè)兀素相比，刖者比后者稍重要5表示兩個(gè)兀素相比，前者比后者明顯重要7表示兩個(gè)兀素相比，刖者比后者強(qiáng)烈重要9表示兩個(gè)兀素相比，刖者比后者極端重要2,4,6,8表示上述判斷的中間值倒數(shù)若兀素1與兀素j的重要性之比為aij,則兀素j與兀素1的重要性之比為aii=1/aii設(shè)填寫后的判斷矩陣為A=(aj)nxn，判斷矩陣具有如下性質(zhì)：aij〉0aji=1/ajiaii=1根據(jù)上面性質(zhì)，判斷矩陣具有對(duì)稱性，因此在填寫時(shí)，通常先填寫aii=1部分，然后再僅需判斷及填寫上三角形或下三角形的n(n-1)/2個(gè)元素就可以了。在特殊情況下，判斷矩陣可以具有傳遞性，即滿足等式：aij*ajk=aik當(dāng)上式對(duì)判斷矩陣所有元素都成立時(shí)，則稱該判斷矩陣為一致性矩陣?！景咐治觥渴姓こ添?xiàng)目建設(shè)決策：構(gòu)造判斷矩陣并請(qǐng)專家填寫接前例，征求專家意見，填寫后的判斷矩陣如下：

表2判斷矩陣表AB1B2B3B1C1C2B2C3C4B3C5C6B111/31/3C111C313C513B211C21C41C61B31C1D1D2C2D1D2C3D1D2C4D1D2D115D113D111佔(zhàn)D117D21D21D21D21C5D1D2C6D1D2D111/5D111/3D21D21層次單排序(計(jì)算權(quán)向量)與檢驗(yàn)對(duì)于專家填寫后的判斷矩陣，利用一定數(shù)學(xué)方法進(jìn)行層次排序。層次單排序是指每一個(gè)判斷矩陣各因素針對(duì)其準(zhǔn)則的相對(duì)權(quán)重，所以本質(zhì)上是計(jì)算權(quán)向量。計(jì)算權(quán)向量有特征根法、和法、根法、冪法等，這里簡要介紹和法。和法的原理是，對(duì)于一致性判斷矩陣，每一列歸一化后就是相應(yīng)的權(quán)重。對(duì)于非一致性判斷矩陣，每一列歸一化后近似其相應(yīng)的權(quán)重，在對(duì)這n個(gè)列向量求取算術(shù)平均值作為最后的權(quán)重。具體的公式是：1na》aklw=—》akl1nj=1k=1需要注意的是，在層層排序中，要對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn)。在特殊情況下，判斷矩陣可以具有傳遞性和一致性。一般情況下，并不要求判斷矩陣嚴(yán)格滿足這一性質(zhì)。但從人類認(rèn)識(shí)規(guī)律看，一個(gè)正確的判斷矩陣重要性排序是有一定邏輯規(guī)律的，例如若A比B重要，B又比C重要，則從邏輯上講，A應(yīng)該比C明顯重要，若兩兩比較時(shí)出現(xiàn)A比C重要的結(jié)果，則該判斷矩陣違反了一致性準(zhǔn)則，在邏輯上是不合理的。因此在實(shí)際中要求判斷矩陣滿足大體上的一致性，需進(jìn)行一致性檢驗(yàn)。只有通過檢驗(yàn)，才能說明判斷矩陣在邏輯上是合理的，才能繼續(xù)對(duì)結(jié)果進(jìn)行分析。一致性檢驗(yàn)的步驟如下。第一步，計(jì)算一致性指標(biāo).(consistencyindex)第二步，查表確定相應(yīng)的平均隨機(jī)一致性指標(biāo)?(randomindex)據(jù)判斷矩陣不同階數(shù)查下表，得到平均隨機(jī)一致性指標(biāo).。例如，對(duì)于5階的判斷矩陣，查表得到.=表3平均隨機(jī)一致性指標(biāo).表(1000次正互反矩陣計(jì)算結(jié)果)矩陣階數(shù)123 4 5 678■00矩陣階數(shù)91011 12 131415■第三步，計(jì)算一致性比例.(consistencyratio)并進(jìn)行判斷R.I.當(dāng).＜時(shí)，認(rèn)為判斷矩陣的一致性是可以接受的，.＞時(shí)，認(rèn)為判斷矩陣不符合一致性要求，需要對(duì)該判斷矩陣進(jìn)行重新修正。【案例分析】市政工程項(xiàng)目建設(shè)決策：計(jì)算權(quán)向量及檢驗(yàn)上例計(jì)算所得的權(quán)向量及檢驗(yàn)結(jié)果見下：表4層次計(jì)算權(quán)向量及檢驗(yàn)結(jié)果表A單(總A單(總)排序權(quán)值B1B2B3CRC1單排序權(quán)值D1D2CRB1單排序權(quán)值C1C2CRC2單排序權(quán)值D1D2CRB2單排序權(quán)值C3C4CRC3單排序權(quán)值D1D2CRB3單排序權(quán)值C5C6CRC4單排序權(quán)值D1D2CRC5 單排序權(quán)值C5 單排序權(quán)值D1D2CRC6單排序權(quán)值D1D2CR可以看出，所有單排序的.＜，認(rèn)為每個(gè)判斷矩陣的一致性都是可以接受的。層次總排序與檢驗(yàn)總排序是指每一個(gè)判斷矩陣各因素針對(duì)目標(biāo)層(最上層)的相對(duì)權(quán)重。這一權(quán)重的計(jì)算采用從上而下的方法，逐層合成。很明顯，第二層的單排序結(jié)果就是總排序結(jié)果。假定已經(jīng)算出第k-1層m個(gè)元素相對(duì)于總目標(biāo)的權(quán)重W(k-1)=(W](k-1),W2(k-1),...,Wm(k-1))T,第k層n個(gè)元素對(duì)于上一層(第k層)第j個(gè)元素的單排序權(quán)重是呼)=(卩嚴(yán)卩瀘,...,卩詐))丁，其中不受j支配的元素的權(quán)重為零。令P(k)=(p1(k),p2(k),^,pn(k)),表示第k層元素對(duì)第k-1層個(gè)元素的排序，則第k層元素對(duì)于總目標(biāo)的總排序?yàn)椋簃Wi(k)=乙W(k)=(W1(k),W2(k),.,WnmWi(k)=乙Pj.(k)Wj(k-1)I=1,2,.,nj=1

同樣，也需要對(duì)總排序結(jié)果進(jìn)行一致性檢驗(yàn)。假定已經(jīng)算出針對(duì)第k-1層第j個(gè)元素為準(zhǔn)則的、和j=1,2,...,m,則第k層的綜合檢驗(yàn)指標(biāo)（,,…,）w（k-1）（,,.,）w（k-1）C.C.R.(k)C.I.(k)R.I.(k)當(dāng)?（k）＜