2023年数模实验报告计算机

河北大学《教学模型》实龄实验报告班级专业计算机班级专业计算机1班姓名实验地点C1-229指导老师实验项目层次分析模型张建华学号司建辉成绩实验9-1层次分析模型一一正互反阵最大特性根和特性向量的实用算法一、实验目的1、了解层次分析模型；2、了解累法函数及和法函数；3、掌握根法函数；二、实验规定.用给出的幕法函数求A的最大特性根和特性向量。调用命令和运营结果。.用给出的和法函数求A的最大特性根和特性向量。调用命令和运营结果。.编写根法函数,用该函数求A的最大特性根和特性向量。程序、调用命令和运营结果。三、实验内容1、幕法:%输入参数：A为正互反方阵，精度e%输出参数：lamda为最大特性根，w为其归一化特性列向量function[lamdaw]=MIfun(a,e)if(nargin==l)%假如只输入一个变量(即A),则㊀取0.000001e=0.000001；endn=size(A,1);%取人的行数wO=ones(n,1)/n;%归历来量whi1e1ww=A*wO;w=ww/sum(ww)；%归一化ifa11(abs(w-wO)<e)break;endw0=w；end1amda=sum(ww./wO)/n;2、和法：%输入参数:A为正互反方阵%输出参数：lamda为最大特性根,w为其归一化特性列向量function[1amdaw]=HEfun(A)n=size(A,1);fork=l:n%a.将A的每一列向量归一化AA(：,k)=A(：,k)./sum(A(:,k));endww=sum(AA,2);%b.对AA按行求和，ww为列向量w=ww./sum(ww);枇.归一化，得w为近似特性列向量lamda=sum(A*w./w)/n;%d.计算最大特性根的近似值人3、根法：function[1amdaw]=RTfun(A)n=size(A,1);fork=l:nAA(:,k)=A(:,k)./sum(A(：,k))；endww=prod(AA,2)."(1/n)；w=ww./sum(ww);1amda=sum(A*w./w)/n;四、实验结果及其分析lainda=3.0092幕法:

0.58760.32340.0890lamda=3.00920.58690.3238和法:0.0893和法:lamda=3.0092根法:

0.58760.32340.0890实验9-2层次分析模型一一旅游决策问题一、实验目的掌握层次分析法;二、实验规定

.显示第2层的以下数据:最大特性根入；特性向量（权向量）w;一致性指标CI;一致性比率CR。.显示第3层的以下数据：特性向量（权向量）W（3）k；最大特性根入k；一致性指标C工k。.显示最下层（第3层）对目的（第1层）的组合权向量。.显示第2层和第3层的组合一致性比率，以及最下层对第1层的组合一致性比率。三、实验内容1、脚本式M文献：formatcompact;A=[11/2433;21755A=[11/2433;21755A=[11/2433;21755;l/41/711/21/3;1/31/5211；l/31/5311];B1=[125;1/212;1/51/21];A=[11/2433;21755B2=[l1/31/8；B2=[l1/31/8；311/3;831];B3=[113;113;1/31/31];4;1/311;1/411];B5=[111/4；111/4;B5=[111/4；111/4;4411;[lamda2w2CI2CR2]=TRfun(A);ifCR2>=0.1A»ODI"A»ODI"1yd»OADOMi2e£i7])A»ODI"1yd»OADOMi2e£i7])disp(A»ODI"1yd»OADOMi2e£i7])end[1amda3_lw3_lCI3_1CR3_l]=TRfun(Bl);ifCR3_l>0.1disp([,CR3_1=,,num2str(CR3_1),!>0.1,B1A»ODI**1yb»OADO^i2e£j1])return；end[1amda3_2w3_2CI3_2CR3_2]=TRfun(B2)；ifCR32>0.1disp([7CR3_2=1,num2str(CR3_2),!>0.1,B2A»ODI1yb»OADOMi2e£j1])return;end[lamda3_3w3_3CI3_3CR3_3]=TRfun(B3)；ifCR3_3>0.1disp([1CR3_3=!,num2str(CR3_3),'>0.1,B3A»ODI**1yb»OADOMi2e£;)return;end[1amda3_4w3_4CI3_4CR3_4]=TRfun(B4)；ifCR3_4>0.1disp([7CR3_4=1rnum2str(CR3_4),'>0.1,B4A»ODI"1yd»OADOMi2e£j7])return；endLlamda3_5w3_5CI3_5CR3_5]=TRfun(B5);ifCR3_5>0.1disp([1CR3_5=1,num2str(CR3_5),1>0.1,B5A»ODI"1yd»0AD0^i2e£i'])return；endw3=[w3_l,w3_2,w3_3,w3_4,w3_5]*w2；CI3=[CI3_1CI3_2CI3_3CI3_4CI3_5]*w2;Rl=[000.580.901.121.241.321.411.451.491.51]；RI3=[Ri(3)RI(3)RI(3)RI(3)RI(3)]*w2;CR3=CI3/RI3;ifCR3>0.1disp([,CR3=,,num2str(CR3),7>0.1,pU31yxe°IO»dADOMi2e£j1])return；endCr=cR2+Cr3；ifCr>0.1disp([7CR=',num2str(CR),f>0.1,A»6DI"1yxe°i6»OADOMi2e£i1])return;end[1amda2w2CI2CR2]=TRfun(A)[lamda3_lw3_1CI3_1]=TRfun(Bl)[1amda3_2w3_2CI3_2]=TRfun(B2)[lamda3_3w3_3CI3_3]=TRfun(B3)Elamda3_4w3_4Cl3_4]=TRfun(B4)[1amda3_5w3_5CI3_5]=TRfun(B5)w3CR2CR3CR2、函数式M文献：%输入参数:A为成对比较阵(正互反方阵)魅俞出参数：%1amda为最大特性根值%w为A的归一化特性列向量(权向量)%CI为一致性指标值%CR为一致性比率值function[lamdawCICR]=TRfun(A)[lamdaw]=HEfun(A);%求人的最大特性根及归一化特性列向量lamdta.3_4.=3.0092w3_4=O.6327O.1924O.1749CI3_4=O.0046lamd.a.3_5=3w3_5=lamdta.3_4.=3.0092w3_4=O.6327O.1924O.1749CI3_4=O.0046lamd.a.3_5=3w3_5=O.1667O.1667O.6667CI3_5=Ow3=O.2990O.2454O.4556CR2=O.0163CR3=O.0027CRRI=[O00.580.901.121.241.321.411.451.491.51]；n=size(A,1)；CI=(lamda-n)/(n—1);%一致性指标,CI=0时A为一致阵；CI越大A的不一致限度越严重CR=CI/RI(n)；%一致性比率，CR<0.1时认为A的不一致限度在允许范围之内四、实验结果及分析lamd.a2=5.0729w2=O.2623O.4744O.0545O.09