24小时热门版块排行榜    

查看: 2780  |  回复: 23
【奖励】 本帖被评价22次,作者pgx457600增加金币 16.8

pgx457600

金虫 (初入文坛)


[资源] 模糊层次分析法FAHP中的模糊判断矩阵的求解IDL程序

本人在利用模糊层次分析法FAHP进行多元信息成矿预测时候,需要对模糊互补判断矩进行求解,在网上没有搜到合适的程序,于是自己根据文献资料(附件里面有几篇经典文献)用IDL写了个程序,在此共享出来,供大家参考使用。



;参考文献
;徐泽水. 模糊互补判断矩阵排序的一种算法[J].系统工程学报,2001, 16(4): 311 314.
;陈华友, 赵佳宝. 模糊判断矩阵的相容性研究[J].运筹与管理,2004, 13(1): 44 47.
;姬东朝,宋笔锋,喻天翔.模糊层次分析法及其在设计方案选优中的应用[J].系统工程与电子技术.2006,28(11):1692-1696.

PRO FAHPMATRIX,event
;输入矩阵A的阶数
N=3
;定义输入矩阵A
A=fltARR(n,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.求权重向量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;权重向量
;  print,'j=',j+1,'    ;Wj=',Wj
;  print,A[0:N-1,j:j]
endfor
    print,'权重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
;矩阵转置
        Wij=TRANSPOSE(Wij)
  print,'特征矩阵Wij='
  print,Wij
;3.求判断矩阵A 与的特征矩阵Wij的相容性指标
;矩阵转置
TWij=TRANSPOSE(Wij)

;定义累积指数相容度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

;相容性指数IAB
IAB=1.0*FC/(N*N)
print,'相容性指数IAB=',IAB

IF IAB LT 0.10  THEN BEGIN

PRINT, '通过模糊互补判断矩阵的一致性检验 '
ENDIF ELSE BEGIN
PRINT, '不能通过模糊互补判断矩阵的一致性检验 '
ENDELSE
print,'Process Time :', systime(1) - start_time
END

程序运行结果:
A=     
     0.500000     0.700000     0.800000
     0.300000     0.500000     0.700000
     0.200000     0.300000     0.500000
权重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
相容性指数IAB=    0.0995591
通过模糊互补判断矩阵的一致性检验

************
根据陈华友(2004)的数据计算结果如下:
;陈华友, 赵佳宝. 模糊判断矩阵的相容性研究[J].运筹与管理,2004, 13(1): 44 47.

A=     0.500000     0.700000     0.800000
     0.300000     0.500000     0.700000
     0.200000     0.300000     0.500000
权重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
相容性指数IAB=    0.0995591
通过模糊互补判断矩阵的一致性检验

文献中的另外一个计算案例如下:
A=     0.500000     0.600000     0.700000
     0.400000     0.500000     0.600000
     0.300000     0.400000     0.500000
权重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
相容性指数IAB=    0.0554613
通过模糊互补判断矩阵的一致性检验[ Last edited by pgx457600 on 2013-5-12 at 08:46 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 模糊互补判断矩阵FAHP求解IDL源代码0416.rar
  • 2013-05-11 16:37:51, 444.42 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (初入文坛)


本帖仅楼主可见
4楼2013-05-18 20:04:38
已阅   申请EPI   回复此楼   编辑   查看我的主页

分水岭

新虫 (小有名气)


★★★★★ 五星级,优秀推荐

楼主太有爱了,
6楼2013-10-02 15:26:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

花开有声lzq

铁虫 (初入文坛)


★★★★★ 五星级,优秀推荐

赞一个
7楼2013-10-07 11:14:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

richdog

银虫 (初入文坛)


11楼2014-04-16 18:36:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mutusas

新虫 (初入文坛)


★★★ 三星级,支持鼓励

????
19楼2015-08-15 20:36:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
liucfcuil2楼
2013-05-12 05:41   回复  
五星好评  顶一下,感谢分享!
cnuhzy3楼
2013-05-13 13:13   回复  
五星好评  顶一下,感谢分享!
noaags5楼
2013-07-15 08:45   回复  
五星好评  顶一下,感谢分享!
2013-12-05 15:52   回复  
五星好评  顶一下,感谢分享!
型格9楼
2014-01-09 14:23   回复  
五星好评  顶一下,感谢分享!
richdog10楼
2014-04-16 18:35   回复  
五星好评  顶一下,感谢分享!
chenxin81212楼
2014-04-17 19:55   回复  
五星好评  顶一下,感谢分享!
liudi2325313楼
2014-06-12 09:55   回复  
五星好评  顶一下,感谢分享!
chenlj051314楼
2014-06-13 15:24   回复  
五星好评  顶一下,感谢分享!
hbj30615楼
2014-06-14 00:25   回复  
五星好评  顶一下,感谢分享!
mathfeng16楼
2014-07-19 10:54   回复  
五星好评  顶一下,感谢分享!
270808010717楼
2015-06-30 12:26   回复  
五星好评  顶一下,感谢分享!
liaotao2518楼
2015-07-17 15:59   回复  
五星好评  顶一下,感谢分享!
wangkexiao20楼
2015-10-05 11:05   回复  
五星好评  顶一下,感谢分享!
anseuji21楼
2016-03-31 11:43   回复  
三星好评  顶一下,感谢分享!
baixue718922楼
2016-04-26 22:36   回复  
五星好评  顶一下,感谢分享!
vinceliu23楼
2016-06-20 14:49   回复  
一般  顶一下,感谢分享!
托帕石24楼
2017-01-04 12:15   回复  
五星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 pgx457600 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见