层次分析法-数学建模实验报告.docVIP

《数学建模与数学实验》实验报告 实验1 离散模型 专业、班级 信息 学号 姓名 课程编号 实验类型 验证性 学时 2 实验（上机）地点 教七楼数学实验中心 完成时间 任课教师 评分 一、实验目的及要求 1．掌握数学软件Matlab，c++的基本用法和一些常用的规则，能用该软件进行编程； 2．能够借助数学软件进行离散数学模型问题的求解和分析； 3．理解层次分析数学模型和效益的合理分配数学模型的数学原理，并能够分别利用这两种离散数学模型进行实际问题的建模。 二、借助数学软件，研究、解答以下问题 一位同学准备购买一部手机，他考虑的因素有质量，颜色，价格，外形，实用，品牌等因素，比较中意的手机有苹果，三星，HTC。但不知选择哪一款为好，请你建立数学模型给他一个好的建议。 （一）建立层次结构模型 构造成对比较矩阵 （三）最大特征值及其对应的特征向量的求解（调用MATLAB函数、和法、根法、幂法） 1.调用MATLAB函数 部分代码如下： %调用函数eig求特征值和特征向量 function [cr,ci,ri,b,lambda]=EIG(A) N=size(A); n=N(1); [x,y]=eig(A); eigenvalue=diag(y); lambda=eigenvalue(1); b=x(:,1); ci=(lambda-n)/(n-1); RI=[0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45]; ri=RI(n); cr=ci/ri; end 和法 部分代码如下： %特征值和特征向量的简化计算——和法（取列向量的算术平均） function [cr,ci,ri,w,lambda]=HEFA(A) [M,N]=size(A); B=zeros(M,N); w=ones(M,1); Aw=ones(M,1); %列向量归一化 for i=1:N B(:,i)=A(:,i)./(sum(A(:,i))); end %算术平均 for j=1:M w(j,1)=sum(B(j,:))/M; end %A*w for i=1:M for j=1:N A(i,j)=A(i,j)*w(j,1); end Aw(i,1)=sum(A(i,:)); end %求lambda lambda=sum(Aw./w)/M; %求cr,ci,ri ci=(lambda-M)/(M-1); RI=[0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45]; ri=RI(M); cr=ci/ri; end 根法 根法计算思想与和法类似，不同之处在于列向量归一化后，和法取其算术平均，根法取其几何平均。 部分代码如下（与和法不同之处）： %几何平均 for i=1:M for j=1:N w(i,1)=w(i,1)*B(i,j); end w(i,1)=w(i,1)^(1/M); end 幂法 部分代码如下： %特征值和特征向量的简化计算——幂法（迭代算法） function [cr,ci,ri,b,lambda]=MIFA(A,w0) %设置精度 epsilon=10^(-4); MAXi=1; N=size(A); n=N(1); W=zeros(50,n); w=zeros(50,n); w(1,:)=w0; k=1; while MAXiepsilon %计算A*w(k,:) for i=1:n for j=1:n A(i,j)=A(i,j)*w(k,j); end W(k,i)=sum(A(i,:)); end %归一化 w(k+1,:)=W(k,:)./sum(W(k,:)); %计算MAXi MAXi=max(abs(w(k+1,:)-w(k,:))); k=k+1; end k=k-1; lambda=(sum(W(k,:)./w(k,:))/n); b=(w(k,:)); ci=(lambda-n)/(n-1); RI=[0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45]; ri=RI(n); cr=ci/ri; end （四）层次单排序及一致性检验 调用不同的求特征值和特征向量的函数来求解层次单排序的方法类似，这里以调用和法的函数为例： 部分代码如下： %成对比较矩阵 A=[1,1,1,4,1,1/2;1,1,2,4,1,1/2;1,1/2,1,5,3,1/2;... 1/4,1/4