1基于层次分析法的大学生志愿选择模型摘要本文主要讨论了大学生毕业后志愿选择的问题。针对问题,利用层次分析法将决策问题分解为目标层(志愿)、准则层(贡献、收入、发展、声誉、人际关系及地理位置)和方案层(工作、学习及其他)。通过成对比较法确定各准则对于目标的权重及各方案对于各准则的权重,构造出准则层对目标层和方案层对准则层的成对比较阵,建立层次结构模型并用MATLAB程序计算各成对比较阵的权向量以及方案层对目标层的组合权向量,得到的结论如下:大学生毕业后志愿选择时工作、学习和其他的权重分别为0.4864、0.2630和0.2506。可见选择工作、学习和其他的大学生分别占总人数的48.64%、26.30%和25.06%。关键词层次分析法;成对比较阵;权重;MATLAB2一、问题分析在日常生活中经常会碰到许多决策问题,在解决这些决策问题时通常会使用离散模型。以是否能够发挥自己的才干为国家做贡献、丰厚的收入、适合个人的兴趣及发展、良好的声誉、人际关系和地理位置六个方面为大学生毕业后志愿选择的主要因素,选择方案有三种,即工作、学习、其他。运用层次分析法得到指标评价体系,建立大学生志愿选择的层次结构模型,利用相对比较矩阵求得各项指标的权向量,给出大学生青年志愿选择得分并进行分析。二、问题假设1.假设调查的数据是合理的;2.假设除已经考虑的因素之外的其他因素对评价模型造成的影响很小,可以不予考虑。三、模型的建立与求解经过讨论,确定大学生青年志愿选择的主要指标为是否能够发挥自己的才干为国家做贡献、丰厚的收入、适合个人的兴趣及发展、良好的声誉、人际关系和地理位置六个方面。利用层次分析法([1])确定大学生志愿选择作为目标层A;以是否能够发挥自己的才干为国家做贡献、丰厚的收入、适合个人的兴趣及发展、良好的声誉、人际关系和地理位置六个方面分别作为准则层1C、2C、3C、4C、5C、6C;以工作、学习、其他分别作为方案层1P、2P、3P。建立以下层次结构(如图1):图1层次结构模型利用成对比较法对准则层、方案层进行列表(表1)3表1-1主要影响因素对比因素贡献收入发展声誉人际关系地理位置贡献11511/23收入113154发展1/51/311/51/31/2声誉115123人际关系21/531/211地理位置1/31/421/311表1-2贡献对比因素工作学习其他工作113/1213/11学习12/13112/11其他11/1311/121表1-3收入对比因素工作学习其他工作132学习1/311/2其他1/221表1-4发展对比因素工作学习其他工作132学习1/311/2其他1/221表1-5声誉对比因素工作学习其他工作125/1625/14学习16/2518/7其他14/257/814表1-6人际关系对比因素工作学习其他工作124学习1/213其他1/41/31表1-7地理环境对比因素工作学习其他工作123学习1/213其他1/31/313.1准则层对目标层的成对比较阵及权向量根据表1数据和层次分析法,给出如下成对比较阵:11511/231131541/51/311/51/31/211512321/531/2111/31/421/311A利用MATLAB可求得最大特征根为5373.6和特征向量)1(w(程序见附录1)。根据公式1161(1)11jjji将特征向量归一化得权向量T0.19750.29470.04940.23070.14460.0830w一致性检验:由一致性指标公式:1nnCI可得6.537360.1074661CI由于6n查矩阵阶数表,如表2所示。5表2矩阵阶数表矩阵阶数1234567891011RI000.580.91.121.241.321.411.451.491.51得到相应的随机一致性指标:1.24RI,则一致性比率为:0.107460.08670.11.24CICRRI故通过一致性检验,即认为A的一致性程度在容许的范围内,可以用归一化后的特征向量w作为排序权重向量。3.2方案层对准则层的成对比较阵及权向量1.1113/1213/1112/13112/1111/1311/121C,最大特征根为3,对应的权向量为:T210.3611,0.3333,0.3056w一致性比率为:2100.1CR;2.21321/311/21/221C,最大特征根为3.0092,对应的权向量为:T220.5396,0.1634,0.2970w一致性比率为:220.00790.1CR;3.31321/311/21/221C,最大特征根为3.0092,对应的权向量为:T220.5396,0.1634,0.2970w一致性比率为:220.00790.1CR;4.4125/1625/1416/2518/714/257/81C,最大特征根为3,对应的权向量为:T210.4545,0.2909,0.2545w一致性比率为:2400.1CR;65.51241/2131/41/31C,最大特征根为3.0183,对应的权向量为:T250.5584,0.3196,0.1220w一致性比率为:250.01580.1CR;6.61231/2131/31/31C,最大特征根为3.0536,对应的权向量为:T260.5278,0.3325,0.1396w一致性比率为:260.04620.1CR。上述结果如表3所示:表3方案层的计算结果k1234562kw0.36110.33330.30560.53960.16340.29700.53960.16340.29700.45450.29090.25450.55840.31960.12200.52780.33250.1396k33.00923.009233.01833.0536kCI00.03950.039500.0790.231kRI0.580.580.580.580.580.58kCR00.00790.007900.01580.0462)6,,2,1(,1.0kCRk,故通过一致性检验,认为)6,,2,1(kCk的一致性程度在容许的范围之内,用归一化后的特征向量作为排序权向量。3.3计算组合权向量并做一致性检验由各准则对目标的权向量w和各方案对每一准则的权向量2kw,计算各方案对目标的权向量,即组合权向量,记作1w,对于方案1P,它在人际关系等6个准则中的权重用1w中的第一个分类表示,所以方案1P在目标中的组合权重应为它们相应项的两两乘积之和,即0.19750.36110.29470.53960.04940.53960.23070.45450.14460.55840.08300.52780.48642P在目标中的组合权重为:0.19750.33330.29470.16340.04940.16340.23070.29090.14460.31960.08300.33250.26303P在目标中的组合权重为:0.19750.30560.29470.29700.04940.29700.23070.25450.14460.12200.08300.13960.25067于是组合权向量T10.4864,0.2630,0.2506w。结果表明:方案1P在大学生志愿选择中所占的权重近于1/2,远大于2P和3P,所以大学生毕业后最有可能选择工作。组合一致性检验可逐层进行。若第p层的一致性指标为)()(1,,pnpCICI,n(是第)1层因素的数目p随机一致性指标为)()(1,,pnpRIRI。)1()()(1)(],,[ppnppCICICIw)1()()(1)(],,[ppnppRIRIRIw则第p层的组合一致性比率为:spRICICRppp,,4,3,)()()(第3层的一致性指标为0,0.0395,0.0395,0,0.079,0.231,随机一致性指标为)6,,2,1(58.0kRIk,所以0.197500.29470.03950.04940.0395CI0.230700.14460.0970.08300.2310.0442;0.19750.580.29470.580.04940.58RI0.23070.580.14460.580.08300.580.5800;此时0.07620.1CR1.0CR,故组合一致性检验通过,组合权向量符合一致性要求。根据组合权向量的分量取值,大学生毕业后最有可能选择工作。参考文献[1]姜启源,谢金星.数学模型[M].北京:高等教育出版社,2001.8附录1.clear;clc;A=[11511/231131541/51/311/51/31/211512321/531/2111/31/421/311];n=length(A);[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;endq=qci=(b(1,1)-n)/(n-1)ri=1.24;cr=ci/ri2.clear;clc;A=[113/1213/1112/13112/1111/1311/121];n=length(A);[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;endq=qci=(b(1,1)-n)/(n-1)ri=0.58;cr=ci/ri3.clear;clc;9A=[1321/311/21/221];n=length(A);[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;endq=qci=(b(1,1)-n)/(n-1)ri=0.58;cr=ci/ri4.clear;clc;A=[1321/311/21/221];n=length(A);[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;endq=qci=(b(1,1)-n)/(n-1)ri=0.58;cr=ci/ri5.clear;clc;A=[125/1625/1416/2518/714/257/81];n=length(A);10[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;endq=qci=(b(1,1)-n)/(n-1)ri=0.58;cr=ci/ri6.clear;clc;A=[1241/2131/41/31];n=length(A);[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;endq=qci=(b(1,1)-n)/(n-1)ri=0.58;cr=ci/ri7.clear;clc;A=[1231/2131/31/31];n=length(A);[a,b]=eig(A)b=b(1,1)s=sum(a(:,1));fori=1:nq(i,1)=a(i,1)/s;end11q=qci=(b(1,1)-n)/(n-1)ri=0.58;cr=ci/ri