文章編號:1003-6199(2011)04-0001-04
摘 要:糾正對1-9標度和指數標度之間關系的錯誤理解,并指出現有一致性改進方法的不合理之處。將調整1-9標度判斷矩陣一致性的方法推廣至指數標度判斷矩陣;并提出基于感覺判斷矩陣的新的一致性調整方法。通過實例說明如何使用這兩種方法,并證明這兩種方法的有效性。
關鍵詞:層次分析法;指數標度; 一致性
中圖分類號: O223 文獻標識碼:A
Research on Methods for Improving Consistency of Judgement Matrix Based on Exponential Scale
DONG Shuai, WANG Yongji, LIU Lei
(Depontment of Control Science and Engineering, Huazhong University of Science and Technology,
Key Laboratory of Image Rrocessing and Intelligent Contral, Wuhan 430074, China)
Abstract:First point out some mistakes about the relation between 1-9 scale and exponential scale in some papers. Then apply the method of improving consistency of 1-9 scale judgment matrix to exponential scale matrix, and propose a new method based on sensory matrix. At last, illustrate how to use these two methods and justify the validation of them through an example.
Key words:analytic hierarchy process; scale; consistency
1 引 言
層次分析法(Analytic hierarchy process, AHP)自上世紀80年代引入我國以來,引起了學者們的廣泛關注。因為1-9標度存在與人的思維不一致等問題,所以標度問題一直是學者研究的焦點之一。由于指數標度能克服1-9標度的缺點[1-3],舒康等[4]提出的基于指數標度的AHP法得到了呂躍進等的推崇[5,6]。有人曾給出了指數標度判斷矩陣一致性的檢驗標準和改進方法[7-9],但這些改進方法過于繁瑣,且意義不明確,所得結果也不合理,因而沒有引起其他學者的關注。針對1-9標度判斷矩陣的一致性,國外學者提出過多種調整方法[10-12],本文將Saaty提出的方法推廣到指數標度,并提出了基于感覺判斷矩陣的調整方法。這兩種方法均能有效地改進指數標度判斷矩陣的一致性,且意義明確,計算量小。
2 指數標度AHP法
與1-9標度AHP的不同,基于指數標度的AHP需要先構造感覺判斷矩陣,再將感覺判斷矩陣轉化為指數標度判斷矩陣。設某屬性支配的下層指標數量為n,則計算這些指標相對于該屬性的重要性排序的步驟為[4]:
1)根據表1所示的指數標度,填寫感覺判斷矩陣B=(bij)n×n,其中,bij為指標i與指標j的重要性之差;
2)將感覺判斷矩陣轉化為指數標度判斷矩陣=(ij)(n×n)=kbij,其中,k有多種取值方法,本文取Saaty較為認可的k=89;
3)計算的最大特征根λmax ;
4)按式C.R.=C.I.R.I.計算的一致性指標,其中,C.I.=λmax -nn-1,RI按表2取值[7];
5)如C.R.≤0.1,則一致性滿足要求,λmax對應的特征向量即為所求的相對重要性排序;否則需對矩陣B進行調整。
可知,矩陣,B滿足如下條件:
ii=1,i=1,2,…,n(1)
bii=0,i=1,2,…,n(2)
ij=-ji,i,j=1,2,…,n(3)
bij=-bij,i,j=1,2,…,n(4)
如果矩陣完全一致,則存在:
ij=kjki,k=1,2,…,n(5)
bij=bkj-bki,k=1,2,…,n(6)
可相應地認為滿足(2)、(4)和(6)式的感覺判斷矩陣完全一致。
值得注意的是,文獻[6]認為1-9標度判斷矩陣(后文將“1-9標度判斷矩陣”直接簡稱“判斷矩陣”)A=(aij)n×n和指數標度判斷矩陣可以按式
ij=(89)aij-1(7)
直接轉換,并在此基礎上證明了兩種標度的一致性是相互矛盾的。但從指數標度判斷矩陣構造過程可可以看出,兩種標度之間不能直接轉換。例如,假設3個指標的重要性等級分別為1,2,4,則1-9標度判斷矩陣為
A=10.50.25210.5421,
感覺判斷矩陣和指數標度判斷矩陣分別為
B=0-1-310-2320,
=10.75980.43871.316110.57742.27951.73241≠9A-18。
可知,CIA=CI=0。故矩陣A和不存在直接的轉換關系,但他們的一致性是相容的。
另外,文獻[1-3]中例2的判斷矩陣為
A=19149149189-1419-389-129-1493819-189-189129181,
此時,CI=0.0287;調整后,
A=19149189189-1419-149-149-18914119-1891411,
CI=0.0578>0.0287,即該方法并未改進判斷矩陣的一致性。
3 一致性調整方法
31 1-9標度調整方法的推廣
由(1)、(3)式可知,也是正互反矩陣,適用于1-9標度的一致性調整方法[10,11],稍加修改后也可用于指數標度。其主要步驟為:
1)計算的主特征向量;
2)計算擾動量εij=ij#8226;j/i(或εij=ij-i/j),得到擾動矩陣E;
3)根據擾動矩陣E確定中擾動量較大的1個或多個元素pq為待調整元素,并依據具體問題來確定調整的元素;
4)令bpq=round(8#8226;log (p/q)9),bqp=-bpq得到新的感覺判斷矩陣B1,若bqp>8,則bqp=8#8226;sgn (bqp);(round為取整運算,sgn為取符號運算);
5)根據B1求取新的1;
6)判斷1的一致性,如不滿足要求,重復1)~4)步。
32 基于感覺判斷矩陣的方法
3.1節中的方法是根據來確定待修正的元素,那是否可以直接由B來確定呢?的一致性決定于B的一致性,且的正互反特性也對應于(2)、(4)式,所以可以直接由B來確定待修正元素。
設指標i的重要性等級為wi,則
=w1-w1w1-w2…w1-wnw2-w1w2-w2…w2-wnwn-w1wn-w2…wn-wn(8)
為完全一致的感覺判斷矩陣。加入擾動量后的感覺判斷矩陣為
B=+ε=(wi-wj+εij)n×n(9)
可得
∑nj=1bij=n#8226;wi-∑ni=1wi+∑nj=1εij(10)
對(10)式兩兩相減可得
wi-wj=-εi-εjn+bi-bjn(11)
其中,εi=∑nj=1εij, bi=∑nj=1bij。當B不一致時,可以認為其所有元素都存在擾動量,即擾動量近似均勻分布于各個元素。則
εij=bij+εi-εjn-bi-bjn
≈bij-bi-bjn(12)
由(12)式可得擾動矩陣E為
E=(εij)n×n≈0ε12ε13…ε1nε210ε23…ε2nε31ε320ε3nεn1εn2εn3…0
=0Δ0-Δ0(13)
E是關于主對角線對稱的,故可只考慮其上三角陣Δ。若Δ中絕對值較大的一個或幾個元素為εpq,則可選擇bpq和bqp為待調整元素。令bpq=bpq-εpq,bqp=-bpq,即減去擾動量,得到新的感覺判斷矩陣B。
由上述推理過程可得到基于感覺判斷矩陣的一致性改進方法:
(1)計算擾動矩陣E=(εij)n×n,其中,εij=bij-∑nk=1bik-∑nk=1bjkn;
(2)選取擾動量較大的一個或多個元素bpq作為確定待修正的元素;
(3)令bpq=bpq-εpq,bqp=-bpq得到新的感覺判斷矩陣B,若bqp>8,則bqp=8#8226;sgn (bqp);
(4)判斷調整后的一致性指標是否滿足要求,如不滿足,重復(1)~(3)步。
4 實例分析
隨機產生感覺判斷矩陣
B=0230-2-20503-3-508000-808230-80,
則
=11.7322.28010.5770.57713.94810.4390.4390.253191110.111191.7322.28010.1111,
=0.3640.4250.6190.4630.297′。
此時,CI=0.9095,CR=1.2632>0.1,原始判斷矩陣的一致性非常差,需要進行調整。
方法一:
按照εij=ij#8226;j/i計算擾動矩陣
E1=12.0223.8751.2700.4710.49515.7491.0880.3070.2580.17416.7260.4800.7870.9190.14915.7832.1223.2602.0830.1731。
ε13,ε23,ε34,ε45明顯較其他元素大,故初步選定13,23,34,45及其對稱位置的元素為待調整元素。
按3.1節方法先單獨調整34和43,有
B1=0230-2-20503-3-501000-108230-80。
此時,CI1=0.4099,CR1=0.5693>0.1。判斷矩陣一致性指標得到改善,但仍不滿足要求。現對4組元素都進行調整,有
B2=02-20-2-20-10-32101000-102230-20
此時,CI2=0.0535,CR1=0.0743<0.1,一致性滿足要求。
方法二:
根據(8)~(11)式,計算3.2中算例擾動矩陣得
E2=01.42.4-0.6-3.2-1.4050-3.6-2.4-508-0.60.60-807.43.23.60.6-7.40
同樣可初步確定b23,b25,b34,b45及其對稱元素為待修正元素。若只調整擾動量最大的b34和b43,可得調整后的CI3=0.3749,CR3=0.5207>0.1;若對4組元素均進行調整,則CI5=0.0590,CR5=0.0819<0.1,一致性滿足要求。
由計算結果可知,兩種方法均能有效地改善指數標度判斷矩陣的一致性。但由于計算過程存在近似、取整和限幅等運算,會使誤差積累,可能導致兩種方法均不能達到要求。因此在初步確定待調整元素后,可以依次選擇調整元素和調整幅度進行嘗試,選取效果最好的調整方案。此外,還要考慮待決策問題的實際情況和決策者的判斷。例如,當εpq相對較大時,如果決策者對于指標p和指標q的關系非常確定,則可能是在構造矩陣B時出現邏輯錯誤。由(12)式可知,矩陣B(或)的一致性越高,對于擾動越敏感,即所得到的擾動量越精確,調整效果會越好。
5 結 論
本文所討論的兩種方法,均是根據擾動量來選擇調整的元素。二者區別在于前者是利用指數標度判斷矩陣來選擇,后者則是直接根據感覺判斷矩陣來選擇。兩種方法均能有效改善一致性指標,但調整元素和幅度的確定沒有嚴格的標準,需要決策者根據待決策問題背景進行具體分析。
參考文獻
[1] 呂躍進, 張維, 曾雪蘭. 指數標度與1-9標度互不相容及其比較研究[J]. 工程數學學報, 2003, 20(8): 77-81.
[2] 侯岳衡, 沈德家. 指數標度及其與幾種標度的比較[J]. 系統工程理論與實踐, 1995, 15(10): 43-46.
[3] 駱正清, 楊善林. 層次分析法中幾種標度的比較[J]. 系統工程理論與實踐, 2004, 24(9):51-60.
[4] 舒康, 梁鎮韓. AHP中的指數標度法[J]. 系統工程理論與實踐, 1990, 10(1): 6-8.
[5] 劉樹林, 邱菀華, 張瑞清. 再論AHP中的指數標度法[J]. 系統工程理論與實踐, 1995, 15(10): 78-80.
[6] 呂躍進, 張維. 指數標度在AHP標度系統中的重要作用[J]. 系統工程學報, 2003, 18(5): 452-456.
[7] 呂躍進. 指數標度判斷矩陣的一致性檢驗方法[J]. 統計與決策, 2006, 22(18): 31-32.
[8] 呂躍進, 占濟舟, 吳小歡. 基于指數標度判斷矩陣的一致性改進方法[A]. 2005中國控制與決策學術年會論文集(下), 2005. 1899-1902.
[9] 劉穎芬, 占濟舟. 基于指數標度下判斷矩陣一致性調整的新算法[J]. 統計與決策 2007, 23(9): 18-19.
[10]SAATY T L. Decision-making with the AHP: Why is the principal eigenvector necessary[J]. European journal of operational research, 2003, 145(1):85-91.
[11]SAATY T L. Fundamentals of decision making and priority[M]. Pittsburgh, PA: RWS Publications, 1994.
[12]LAMATA M T, PELAEZ J I. A method for improving the consistency of judgements[J]. International journal of uncertainty, fuzziness and knowledgebased systems, 2002:10(6): 677-686.
收稿日期:2011-10-18
基金項目:國家自然科學基金項目(60975058,61075095);康復機器人863項目(2008AA04Z207);湖北省自然科學基金項目(2010CDB01904);2010航天支撐技術基金項目;華中科技大學國防自主創新研究基金項目;教育部重點實驗室基金項目
作者簡介:董 帥(1986—),男,湖北天門人,博士研究生,研究方向:飛行器制導控制與軌跡優化(E-mail:dongshuai1986@smail.hust.edu.cn);王永驥(1955—),男,江西吉安人,教授,博士生導師,研究方向:飛行器制導控制與軌跡優化,康復機器人。