查看: 4754|回复: 6

层次分析法5

[复制链接]

字体大小: 正常放大

OLS

60 主题	6 听众	1594 积分

独孤九剑

升级 59.4%

TA的每日心情

	郁闷 2012-4-30 23:06

签到天数: 6 天

[LV.2]偶尔看看I

自我介绍: 　　“独孤九剑”的要旨，在于一眼见到对方招式中的破绽，便即乘虚而入，后发先至，一招制胜，

群组: 南京邮电大学数模协会

群组: 数学散仙

群组: 数模应用

群组: 数学趣味、游戏、IQ等

群组: 2011年第一期数学建模

电梯直达

1^#

发表于 2009-2-3 21:49 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

Matlab语言编写程序如下:
clc
a=[1,  2,  4,  1,  1/2, 1
1/2,1,  1,  1, 5, 1
1/4,1,  1,  2, 1, 1
1,  1,1/2,  1, 1, 1
2,1/5,  1,  1, 1, 1
1,  1,  1,  1, 1, 1];
[x,y]=eig(a);eigenvalue=diag(y);lamda=eigenvalue(1);
ci1=(lamda-6)/5;cr1=ci1/1.24
w1=x(:,1)/sum(x(:,1))
b1=[1,1,2,1; 1,1,1/3,4;1/2,3,1,1;1,1/4,1,1];
[x,y]=eig(b1);eigenvalue=diag(y);lamda=eigenvalue(1);
ci21=(lamda-4)/3;cr21=ci21/0.90
w21=x(:,1)/sum(x(:,1))
b2=[1 1 4 1;1 1 1/5 4;1/4 5 1 1; 1 1/4 1 1];
[x,y]=eig(b2);eigenvalue=diag(y);lamda=eigenvalue(1);
ci22=(lamda-4)/3;cr22=ci22/0.90
w22=x(:,1)/sum(x(:,1))
b3=[1 7 2 1;1/7 1 1/5 4;1/2 5 1 1;1 1/4 1 1];
[x,y]=eig(b3);eigenvalue=diag(y);lamda=eigenvalue(1);
ci23=(lamda-4)/3;cr23=ci23/0.90
w23=x(:,1)/sum(x(:,1))
b4=[1 5 2 1;1/5 1 1/5 2;1/2 5 1 1 ;1 1/2 1 1];
[x,y]=eig(b4);eigenvalue=diag(y);lamda=eigenvalue(1);
ci24=(lamda-4)/3;cr24=ci24/0.90
w24=x(:,1)/sum(x(:,1))
b5=[1 1 9 4 ; 1 1 1/4 2; 1/9 1/2 1 1 ;1/4 1/2 1 1];
[x,y]=eig(b5);eigenvalue=diag(y);lamda=eigenvalue(2);
ci25=(lamda-4)/3;cr25=ci25/0.90
w25=x(:,2)/sum(x(:,2))
b6=[1 1/2 9 5 ;2 1 1/4 3; 1/9 4 1 1 ; 1/5 1/3 1 1];
[x,y]=eig(b6);eigenvalue=diag(y);lamda=eigenvalue(1);
ci26=(lamda-4)/3;cr26=ci26/0.90
w26=x(:,1)/sum(x(:,1))
w_sum=[w21,w22,w23,w24,w25,w26]*w1
ci=[ci21,ci22,ci23,ci24,ci25,ci26];
cr=ci*w1/sum(0.90*w1)

zan