| 24 | 1/1 | 返回列表 |
| 查看: 2525 | 回復(fù): 23 | ||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)22次,作者pgx457600增加金幣 16.8 個(gè) | ||
pgx457600金蟲(chóng) (初入文壇)
|
[資源]
模糊層次分析法FAHP中的模糊判斷矩陣的求解IDL程序
|
|
|
本人在利用模糊層次分析法FAHP進(jìn)行多元信息成礦預(yù)測(cè)時(shí)候,需要對(duì)模糊互補(bǔ)判斷矩進(jìn)行求解,在網(wǎng)上沒(méi)有搜到合適的程序,于是自己根據(jù)文獻(xiàn)資料(附件里面有幾篇經(jīng)典文獻(xiàn))用IDL寫(xiě)了個(gè)程序,在此共享出來(lái),供大家參考使用。 ;參考文獻(xiàn) ;徐澤水. 模糊互補(bǔ)判斷矩陣排序的一種算法[J].系統(tǒng)工程學(xué)報(bào),2001, 16(4): 311 314. ;陳華友, 趙佳寶. 模糊判斷矩陣的相容性研究[J].運(yùn)籌與管理,2004, 13(1): 44 47. ;姬東朝,宋筆鋒,喻天翔.模糊層次分析法及其在設(shè)計(jì)方案選優(yōu)中的應(yīng)用[J].系統(tǒng)工程與電子技術(shù).2006,28(11):1692-1696. PRO FAHPMATRIX,event ;輸入矩陣A的階數(shù) N=3 ;定義輸入矩陣A A=fltARR(n,n) ;定義權(quán)重向量矩陣W W=fltARR(n) ;定義判斷矩陣A 的特征矩陣Wij Wij=fltARR(n,n) ;輸入矩陣A ;A=[[0.5,0.3,0.6,0.7],[0.7,0.5,0.7,0.5],[0.4,0.3,0.5,0.4],[0.3,0.5,0.6,0.5]];4階矩陣 ;A=[[0.5,0.6,0.7],[0.4,0.5,0.6],[0.3,0.4,0.5]];3階矩陣 A=[[0.5,0.7,0.8],[0.3,0.5,0.7],[0.2,0.3,0.5]];3階矩陣 print,'A=',A start_time = systime(1) ;1.求權(quán)重向量W的值 for j=0,N-1 DO BEGIN SUMAj=TOTAL(A[0:N-1,j:j]) ;print, 'SUMAj=', SUMAj Wj=(SUMAj+0.5*n-1)/(n*(n-1)) W[j]= Wj;權(quán)重向量 ; print,'j=',j+1,' ;Wj=',Wj ; print,A[0:N-1,j:j] endfor print,'權(quán)重W=',W ;2.求判斷矩陣A 的特征矩陣Wij的值 for i=0,N-1 DO BEGIN Wi=W ;print,'i=',i+1,wi for j=0,N-1 DO BEGIN Wj=W[j] ;print,'j=',j+1,wj Wij[i,j]=Wi/(Wi+Wj) endfor endfor ;矩陣轉(zhuǎn)置 Wij=TRANSPOSE(Wij) print,'特征矩陣Wij=' print,Wij ;3.求判斷矩陣A 與的特征矩陣Wij的相容性指標(biāo) ;矩陣轉(zhuǎn)置 TWij=TRANSPOSE(Wij) ;定義累積指數(shù)相容度FC FC=0.0 for i=0,N-1 DO BEGIN for J=0,N-1 DO BEGIN FC=FC+ abs(A[i,j]+Wij[J,I]-1) endfor endfor print,'FC=',FC ;相容性指數(shù)IAB IAB=1.0*FC/(N*N) print,'相容性指數(shù)IAB=',IAB IF IAB LT 0.10 THEN BEGIN PRINT, '通過(guò)模糊互補(bǔ)判斷矩陣的一致性檢驗(yàn) ' ENDIF ELSE BEGIN PRINT, '不能通過(guò)模糊互補(bǔ)判斷矩陣的一致性檢驗(yàn) ' ENDELSE print,'Process Time :', systime(1) - start_time END 程序運(yùn)行結(jié)果: A= 0.500000 0.700000 0.800000 0.300000 0.500000 0.700000 0.200000 0.300000 0.500000 權(quán)重W= 0.416667 0.333333 0.250000 特征矩陣Wij= 0.500000 0.555556 0.625000 0.444444 0.500000 0.571429 0.375000 0.428571 0.500000 FC= 0.896032 相容性指數(shù)IAB= 0.0995591 通過(guò)模糊互補(bǔ)判斷矩陣的一致性檢驗(yàn) ************ 根據(jù)陳華友(2004)的數(shù)據(jù)計(jì)算結(jié)果如下: ;陳華友, 趙佳寶. 模糊判斷矩陣的相容性研究[J].運(yùn)籌與管理,2004, 13(1): 44 47. A= 0.500000 0.700000 0.800000 0.300000 0.500000 0.700000 0.200000 0.300000 0.500000 權(quán)重W= 0.416667 0.333333 0.250000 特征矩陣Wij= 0.500000 0.555556 0.625000 0.444444 0.500000 0.571429 0.375000 0.428571 0.500000 FC= 0.896032 相容性指數(shù)IAB= 0.0995591 通過(guò)模糊互補(bǔ)判斷矩陣的一致性檢驗(yàn) 文獻(xiàn)中的另外一個(gè)計(jì)算案例如下: A= 0.500000 0.600000 0.700000 0.400000 0.500000 0.600000 0.300000 0.400000 0.500000 權(quán)重W= 0.383333 0.333333 0.283333 特征矩陣Wij= 0.500000 0.534884 0.575000 0.465116 0.500000 0.540541 0.425000 0.459459 0.500000 FC= 0.499151 相容性指數(shù)IAB= 0.0554613 通過(guò)模糊互補(bǔ)判斷矩陣的一致性檢驗(yàn)[ Last edited by pgx457600 on 2013-5-12 at 08:46 ] |
鐵蟲(chóng) (初入文壇)
| 24 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|