数模实验报告-计算机1班

1、河北大学数学模型实验实验报告 班级专业 计算机1班 姓名 张建华 学号 20144138021 实验地点 C1-229 指导老师 司建辉 成绩 实验项目 层次分析模型 实验9-1层次分析模型一一正互反阵最大特征根和特征向量的实用算法 一、实验目的 1、了解层次分析模型； 2 、了解幕法函数及和法函数； 3、掌握根法函数； 二、实验要求 1. 用给出的幕法函数求A的最大特征根和特征向量。调用命令和运行结果。 2. 用给出的和法函数求A的最大特征根和特征向量。调用命令和运行结果。 3. 编写根法函数，用该函数求 A的最大特征根和特征向量。程序、调用命令和运行结果。 三、实验内容 1、幕法： %俞入

2、参数：A为正互反方阵，精度e %输出参数：lamda为最大特征根，w为其归一化特征列向量 fun ctio n lamda w=Mlfu n(A,e) if(nargin=1)%如果只输入一个变量(即A),则e取0.000001 e=0.000001; end n=size(A,1);% 取A的行数 w0=ones(n,1)/n;% 归一向量 while 1 ww=A*w0; w=ww/sum(ww);归一化 if all(abs(w-w0) 5876 0- 323a 根法： S9 实验9-2层次分析模型一一旅游决策问题 一、实验目的 掌握层次分析法； 二、实验要求 1. 显示第2层的以下数据

3、： 最大特征根入；特征向量(权向量)w; 致性指标CI ; 一致性比率CF。 2. 显示第3层的以下数据： 特征向量(权向量)W(3)k ;最大特征根入k ; 一致性指标CIk。 3. 显示最下层(第3层)对目标(第1层)的组合权向量。 4. 显示第2层和第3层的组合一致性比率，以及最下层对第1层的组合一致性比率。 三、实验内容 1、脚本式M文件： format compact; A=1 1/2 4 33 ;2 17 5 5 ;1/4 1/7 1 1/2 1/3 ;1/3 1/5 2 1 1;1/3 1/5 3 1 1; B1=125 ;1/2 12 ; 1/5 1/2 1; B2=1 1/3

4、 1/8; 3 1 1/3; 8 3 1; B3=113 ;1 13 ; 1/3 1/3 1; B4=134 ;1/3 11 ; 1/4 1 1; B5=1 1 1/4; 1 1 1/4; 4 4 1; Iamda2 w2 CI2 CR2=TRfun(A); if CR2=0.1 disp( CR2= ,n um2str(CR2), 0.1,A? ? o D 1 y o ?D? i 2 e ?) return ; end lamda3_1 w3_1 CI3_1 CR3_1=TRfu n( B1); if CR3_10.1 disp( CR3_1= ,num2str(CR3_1), 0.1,B1

5、? ? o D 1 y o ? ? ? D ? ? i 2 e ?) return ; end lamda3_2 w3_2 CI3_2 CR3_2=TRfu n( B2); if CR3_20.1 disp( CR3_2= ,num2str(CR3_2), 0.1,B2? ? o D 1 y o ? ? ? D ? ? i 2 e ?) return end lamda3_3 w3_3 CI3_3 CR3_3=TRfun(B3); if CR3_30.1 disp( CR3_3= , num2str(CR3_3), 0.1,B3? ? o D 1 y o ?D? i 2 e ?) return

6、 ; end lamda3_4 w3_4 CI3_4 CR3_4=TRfun(B4); if CR3_40.1 disp( CR3_4= ,num2str(CR3_4), 0.1,B4? ? o D 1 y o ? ? ? D ? ? i 2 e ?) return ; end lamda3_5 w3_5 CI3_5 CR3_5=TRfun(B5); if CR3_50.1 disp( CR3_5= ,num2str(CR3_5), 0.1,B5? ? o D 1 y o ? ? ? D ? ? i 2 e ?) return ; end w3=w3_1,w3_2,w3_3,w3_4,w3_5

7、*w2; CI3=CI3_1 CI3_2 CI3_3 CI3_4 CI3_5*w2; RI=0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51; RI3=RI(3) RI(3) RI(3) RI(3) RI(3)*w2; CR3=CI3/RI3; if CR30.1 disp( CR3= ,num2str(CR3), 0.1,卩。3 2 ? ? ? o D i - 1 y xe o ? o ? ? ? D ? ? i 2 e ?) return ; end CR=CR2+CR3; if CR0.1 disp( CR= ,n um2str(CR),

8、0.1, ? ? o D i 1 y xe o ? o ? ? ? D ? ? i 2 e ?) return ; end lamda2 w2 CI2 CR2=TRfun(A) lamda3_1 w3_1 CI3_1=TRfun(B1) lamda3_2 w3_2 CI3_2=TRfun(B2) lamda3_3 w3_3 CI3_3=TRfun(B3) lamda3_4 w3_4 CI3_4=TRfun(B4) lamda3_5 w3_5 CI3_5=TRfun(B5) w3 CR2 CR3 CR 2、函数式M文件： % 输入参数：A为成对比较阵（正互反方阵） %俞出参数： % lamda为

9、最大特征根值 % w为A的归一化特征列向量(权向量) % CI为一致性指标值 % CF为一致性比率值 fun ctio n lamda w CI CR=TRfu n(A) lamda w=HEfun(A); % 求A的最大特征根及归一化特征列向量 %随机一致性指标RI的数值(下标对应成对比较方阵的阶数)： RI=0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51; n=size(A,1); CI=(lamda-n)/(n-1); %一致性指标，CI=0时A为一致阵；CI越大A的不一致程度越严 CR=CI/RI(n); % 一致性比率，CR0.1时认为A的不一致程度在容许范围之内 四、实验结果及分析 UCL 01 磧 R lajriiioi2 丨 &!B D7SS Q. 7-14 O. 口命护 0. 085 , 1 18 GIS = , 132 CR2 = 0- O 163 X aon.d.a.3_1= 3. DOBS wJ_l = 5019 CL 27tSU D* 1S8G 1 0. rifljes 1 sundsi 3_2 = 3. nn i e T一 2 = a. 020 a, 0, SIS CT3_ 2 = 九? 1 04-0口A X akmdzi! 3_ 3 = 3 心 _.3= . 4:2君