1罗应泽 实验五 层次分析法11.15.docVIP

数学建模实验报告 实验序号：5 日期：2012年 12月 3日 班级 10级应数B班 姓名 罗庭佑 学号 201005050255 实验名称 层次分析法的程序设计 实验所用软件及版本 Matlab软件 1、实验目的 （1）熟悉使用Matlab软件完成层次分析法中的特征值、权向量、一致性比率的计算。 （2）进一步加深对层次分析法的流程和具体步骤的理解。 2、实验内容 以教材8.1节旅游问题为例 （1）编写程序代码用于计算第229页对比较阵对应的最大特征值、对应的特征向量以及一致性比率，并将这些结果列表（类似于表3） （2）（选做题）使用MATLAB中求特征值函数命令eig,编写通用性强的求特征值和特征向量的程序代码，即输入任意成对比较阵B，输出对应的最大特征值、特征向量，并判断B是否通过一致性检验。运用定理2（教材P238）判断B是否为一致阵。 （3）（选做题）编写根法的程序代码，即输入任意成对比较阵B，输出对应的最大特征值、特征向量。 3、详细设计（包括算法描述和程序） （1）求比较阵对应的最大特征值、对应的特征向量(将其归一化处理)的程序： format short; %输入待求的矩阵B1 B1=[1 2 5;... 1/2 1 2 ;... 1/5 1/2 1]; [v,d]=eigs(B1); a=max(d(:));%最大特征值 [m,n]=size(v); %得到行数和列数 sum = 0; for i=1:m sum = sum + v(i,1); end tbvector = v(:,1); for i=1:m tbvector(i,1)= v(i,1)/sum; end disp(最大的特征值为：); a disp(最大的特征值对应的特征向量为（标准化后的）：); tbvector format short; %输入待求的矩阵B2 B2=[1 1/3 1/8;... 3 1 1/3;... 8 3 1]; [v,d]=eigs(B2); a=max(d(:));%最大特征值 [m,n]=size(v); %得到行数和列数 sum = 0; for i=1:m sum = sum + v(i,1); end tbvector = v(:,1); for i=1:m tbvector(i,1)= v(i,1)/sum; end disp(最大的特征值为：); a disp(最大的特征值对应的特征向量为（标准化后的）：); tbvector format short; %输入待求的矩阵B3 B3=[1 1 3;... 1 1 3;... 1/3 1/3 1]; [v,d]=eigs(B3); a=max(d(:));%最大特征值 [m,n]=size(v); %得到行数和列数 sum = 0; for i=1:m sum = sum + v(i,1); end tbvector = v(:,1); for i=1:m tbvector(i,1)= v(i,1)/sum; end disp(最大的特征值为：); a disp(最大的特征值对应的特征向量为（标准化后的）：); tbvector format short; %输入待求的矩阵B3 B3=[1 1 3;... 1 1 3;... 1/3 1/3 1]; [v,d]=eigs(B3); a=max(d(:));%最大特征值 [m,n]=size(v); %得到行数和列数 sum = 0; for i=1:m sum = sum + v(i,1); end tbvector = v(:,1); for i=1:m tbvector(i,1)= v(i,1)/sum; end disp(最大的特征值为：); a disp(最大的特征值对应的特征向量为（标准化后的）：); tbvector format short; %输入待求的矩阵B4 B4=[1 3 4;... 1/3 1 1;... 1/4 1 1]; [v,d]=eigs(B4); a=max(d(:));%最大特征值 [m,n]=size(v); %得到行数和列数 sum = 0; for i=1:m