查看: 4378|回复: 8

跪求基于MATLAB的层次分析法程序！！

[复制链接]

字体大小: 正常放大

ninipig

10 主题	4 听众	314 积分

升级 4.67%

TA的每日心情

	无聊 2013-9-6 13:22

签到天数: 7 天

[LV.3]偶尔看看II

电梯直达

1^#

发表于 2009-8-16 11:08 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

哪位高手有基于MATLAB的层次分析法程序，可否贡献出来大家一起分享了~~~
O(∩_∩)O谢谢！

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

yyfanghui010

1 主题	4 听众	45 积分

升级 42.11%

该用户从未签到

2^#

发表于 2009-8-16 11:15 |只看该作者 |招呼Ta 关注Ta

A=[1 2 3 2；1/2 1 2 1;1/2 1 2 1;1/3 1/2 1 1/2;1/2 1 2 1]
[x,y]=eig(A)；%y是特征值且从大到小排列。x是特征矩阵,每一列为相应特征值的一个特征向量。
eigenvalue=diag(y);
lamda=eigenvalue(1);
y_lamda=x(:,1)
v1=y_lamda./norm(x,1) % 单位化
ads(v1) ％绝对化
不知道是不是你要的，A是判断矩阵

我也说一句

发表

使用道具举报

dsshirley

0 主题	4 听众	227 积分

升级 63.5%

该用户从未签到

3^#

发表于 2009-8-17 17:19 |只看该作者 |招呼Ta 关注Ta

没有学习讨论的专区　就没有长期固定的活跃网友

我也说一句

发表

使用道具举报

ninipig

10 主题	4 听众	314 积分

升级 4.67%

TA的每日心情

	无聊 2013-9-6 13:22

签到天数: 7 天

[LV.3]偶尔看看II

4^#

发表于 2009-8-28 08:33 |只看该作者 |招呼Ta 关注Ta

喔，谢谢大家了。O(∩_∩)O~

我也说一句

发表

使用道具举报

chendaolee

0 主题	2 听众	16 积分

升级 11.58%

该用户从未签到

5^#

发表于 2011-6-5 15:40 |只看该作者 |招呼Ta 关注Ta

我也说一句

发表

使用道具举报

muyitianmo

0 主题	7 听众	139 积分

升级 19.5%

TA的每日心情

	开心 2013-4-4 20:07

签到天数: 36 天

[LV.5]常住居民I

自我介绍: up

群组: Matlab讨论组

群组: 数学建模

6^#

发表于 2012-12-7 22:03 |只看该作者 |招呼Ta 关注Ta

好像菏泽个拉丁火把节说的话

我也说一句

发表

使用道具举报

568506561

0 主题	7 听众	30 积分

升级 26.32%

TA的每日心情

	衰 2014-1-10 10:54

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: 我喜欢学习

7^#

发表于 2013-4-27 23:18 |只看该作者 |招呼Ta 关注Ta

yyfanghui010 发表于 2009-8-16 11:15
A=[1 2 3 2；1/2 1 2 1;1/2 1 2 1;1/3 1/2 1 1/2;1/2 1 2 1]
[x,y]=eig(A)；%y是特征值且从大到小排列。x ...

你的标点都错了，不一致，而且矩阵也不是方阵，怎么可能求出来呢？

我也说一句

发表

使用道具举报

568506561

0 主题	7 听众	30 积分

升级 26.32%

TA的每日心情

	衰 2014-1-10 10:54

签到天数: 8 天

[LV.3]偶尔看看II

自我介绍: 我喜欢学习

8^#

发表于 2013-4-27 23:19 |只看该作者 |招呼Ta 关注Ta

%层次分析法的matlab程序
disp('请输入判断矩阵A(n阶)');
A=input('A=');
[n,n]=size(A);
x=ones(n,100);
y=ones(n,100);
m=zeros(1,100);
m(1)=max(x(:,1));
y(:,1)=x(:,1);
x(:,2)=A*y(:,1);
m(2)=max(x(:,2));
y(:,2)=x(:,2)/m(2);
p=0.0001;i=2;k=abs(m(2)-m(1));
while  k>p
  i=i+1;
  x(:,i)=A*y(:,i-1);
  m(i)=max(x(:,i));
  y(:,i)=x(:,i)/m(i);
  k=abs(m(i)-m(i-1));
end
a=sum(y(:,i));
w=y(:,i)/a;
t=m(i);
disp('权向量');disp(w);
disp('最大特征值');disp(t);
      %以下是一致性检验
CI=(t-n)/(n-1);RI=[0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59];
CR=CI/RI(n);
if CR<0.10
disp('此矩阵的一致性可以接受!');
disp('CI=');disp(CI);
disp('CR=');disp(CR);
else
disp('此矩阵的一致性不可以接受!');
end
不知道能否帮到你