数学建模系列-常用模型

数学建模常用模型模型Ⅰ：层次分析法问题1选择旅游地现有三个旅游胜地可供选择，分别为苏杭、黄山、桂林，下面将作出旅游地的选择。面临各种各样的方案，要进行比较、判断、评价、最后作出决策。这个过程主观因素占有相当的比重给用数学方法解决问题带来不便。T.L.saaty等人20世纪在七十年代提出了一种能有效处理这类问题的实用方法。层次分析法（AnalyticHierarchyProcess,AHP)这是一种定性和定量相结合的、系统化的、层次化的分析方法。过去研究自然和社会现象主要有机理分析法和统计分析法两种方法，前者用经典的数学工具分析现象的因果关系，后者以随机数学为工具，通过大量的观察数据寻求统计规律。近年发展的系统分析是又一种方法，而层次分析法是系统分析的数学工具之一。层次分析法（AHP）具体步骤：明确问题递阶层次结构的建立建立两两比较的判断矩阵层次单排序层次综合排序层次分析法的基本步骤1建立层次结构模型一般分为三层，最上面为目标层，最下面为方案层，中间是准则层或指标层。若上层的每个因素都支配着下一层的所有因素，或被下一层所有因素影响，称为完全层次结构，否则称为不完全层次结构。选择旅游地景色费用居住饮食旅途苏杭、黄山、桂林建立选择旅游地层次结构准则层A方案层B目标层ZZ1A2A3A4A5A1B2B3B54321,,,,AAAAA321,,BBB分别分别表示景色、费用、居住、饮食、旅途。分别表示苏杭、黄山、桂林。设某层有个因素，nnxxxX,,,21ijaijjiijaa1nnnnnnnnijaaaaaaaaaaA212222111211An2构造成对比较矩阵要比较它们对上一层某一准则（或目标）的影响程度，确定在该层中相对于某一准则所占的比重。（即把个因素对上层某一目标的影响程度排序）用表示第个因素相对于第个因素的比较结果，则则称为成对比较矩阵。上述比较是两两因素之间进行的比较，比较时取1~9尺度。标度定义与说明1两个元素对某个属性具有同样重要性3两个元素比较，一元素比另一元素稍微重要5两个元素比较，一元素比另一元素明显重要7两个元素比较，一元素比另一元素重要得多9两个元素比较，一元素比另一元素极端重要2,4,6,8表示需要在上述两个标准之间拆衷时的标度1/bij两个元素的反比较旅游问题中，第二层A的各因素对目标层Z的影响两两比较结果如下：ZA1A2A3A4A5A1A2A3A4A511/2433217551/41/711/21/31/31/52111/31/531154321,,,,AAAAA分别表示景色、费用、居住、饮食、旅途。由上表，可得成对比较矩阵1135131112513131211714155712334211A旅游问题的成对比较矩阵共有6个（一个5阶，5个3阶）。问题：两两进行比较后，怎样才能知道，下层各因素对上层某因素的影响程度的排序结果呢？3层次单排序及一致性检验nn层次单排序：确定下层各因素对上层某因素影响程度的过程。用权值表示影响程度，先从一个简单的例子看如何确定权值。例如一块石头重量记为1，打碎分成各小块，各块的重量分别记为：则可得成对比较矩阵11121212121由右面矩阵可以看出，jkkiji即，nji,,2,1,1321231321234,2,7aaaaaaAijkjikaaaijkjikaaaAnjiaaaiijiij,,2,1,,1,1.11.2ArankA的各行成比例，则但在例2的成对比较矩阵中，在正互反矩阵中，若则称为一致阵。一致阵的性质：。特征根均等于个其余的最大特征根（值）为01,.3n-nλAA4.的任一列(行)都是对应于特征根的特征向量。n若成对比较矩阵是一致阵，则我们自然会取对应于最大特征根的归一化特征向量定理：阶互反阵的最大特征根，当且仅当时，为一致阵。Ann11niiwiwinnnA表示下层第个因素对上层某因素影响程度的权值。若成对比较矩阵不是一致阵，Saaty等人建议用其最大特征根对应的归一化特征向量作为权向量，则Anw这样确定权向量的方法称为特征根法.1nnCI由于连续的依赖于，则比大的越多，的不一致性越严重。用最大特征值对应的特征向量作为被比较因素对上层某因素影响程度的权向量，其不一致程度越大，引起的判断误差越大。因而可以用数值的大小来衡量nijanAA的不一致程度。定义一致性指标其中为的对角线元素之和，也为的特征根之和。AnARI50021,,,AAA50021,,,CICICI15005005002150021nnCICICIRI则可得一致性指标定义随机一致性指标随机构造500个成对比较矩阵随机一致性指标RI的数值：n1234567891011RI000.580.901.121.241.321.411.451.491.511.0RICICRAA一致性检验：利用一致性指标和一致性比率0.1及随机一致性指标的数值表，对进行检验的过程。一般，当一致性比率的不一致程度在容许范围之内，可用其归一化特征向量作为权向量，否则要重新构造成对比较矩阵，对加以调整。时，认为A4层次总排序及其一致性检验确定某层所有因素对于总目标相对重要性的排序权值过程，称为层次总排序从最高层到最低层逐层进行。设：Z1A2AmA1B2BnB,,,,21mAAAmA个因素层对总目标Z的排序为maaa,,,21jAAnB中因素为个因素对上层层的层次单排序为),,2,1(,,,21mjbbbnjjj即层第个因素对总目标的权值为：BnmmnnnmmmmbababaBbababaBbababaB22112222211211221111:::Bimjijjba1层的层次总排序为：B层的层次总排序mAAA,,,21maaa,,,21nBBB2112111nbbb22212nbbbnmmmbbb21AB111bbamjjj212bbamjjjnmjnjjbba1层次总排序的一致性检验设层对上层(层)中因素的层次单排序一致性指标为，随机一致性指为，则层次总排序的一致性比率为：BnBBB,,,21A),,2,1(mjAjjCIjRImmmmRIaRIaRIaCIaCIaCIaCR221122111.0CR当时，认为层次总排序通过一致性检验。到此，根据最下层（决策层）的层次总排序做出最后决策。1.建立层次结构模型该结构图包括目标层，准则层，方案层。层次分析法的基本步骤归纳如下3.计算单排序权向量并做一致性检验2.构造成对比较矩阵从第二层开始用成对比较矩阵和1~9尺度。对每个成对比较矩阵计算最大特征值及其对应的特征向量，利用一致性指标、随机一致性指标和一致性比率做一致性检验。若检验通过，特征向量（归一化后）即为权向量；若不通过，需要重新构造成对比较矩阵。计算最下层对最上层总排序的权向量。4.计算总排序权向量并做一致性检验1.0CRCR进行检验。若通过，则可按照总排序权向量表示的结果进行决策，否则需要重新考虑模型或重新构造那些一致性比率较大的成对比较矩阵。mmmmRIaRIaRIaCIaCIaCIaCR22112211利用总排序一致性比率三层次分析法建模举例Z1A2A3A4A5A1B2B3B54321,,,,AAAAA321,,BBB1旅游问题(1)建模分别分别表示景色、费用、居住、饮食、旅途。分别表示苏杭、黄山、桂林。（2）构造成对比较矩阵1135131112513131211714155712334211A1215121215211B1383113813112B131313113113B114111314314B144411141115B(3)计算层次单排序的权向量和一致性检验A073.5110.0,099.0,055.0,475.0,263.0018.0155073.5CI12.1RI1.0016.012.1018.0CRA成对比较矩阵的最大特征值表明通过了一致性验证。故则该特征值对应的归一化特征向量对成对比较矩阵可以求层次总排序的权向量并进行一致性检验，结果如下：54321,,,,BBBBBk1k2k3kkkCIkRI12345595.0082.0429.0633.0166.0277.0236.0429.0193.0166.0129.0682.0142.0175.0668.0005.3002.33009.33003.0001.000005.058.058.058.058.058.0计算可知通过一致性检验。kCR54321,,,,BBBBB对总目标的权值为：1B3.0110.0166.0099.0633.0055.0429.0475.0082.0263.0595.032,BB,456.0,246.0456.0,246.0,3.01.0015.058.0/)0110.0005.0099.00055.0001.0475.0003.0263.0(CR（4）计算层次总排序权值和一致性检验又决策层对总目标的权向量为：同理得，对总目标的权值分别为：故，层次总排序通过一致性检验。可作为最后的决策依据。456.0,246.0,3.0213BBB321,,BBB故最后的决策应为去桂林。又分别表示苏杭、黄山、桂林，即各方案的权重排序为模型Ⅱ：线性规划丁的蛙泳成绩退步到1’15”2；戊的自由泳成绩进步到57”5,组成接力队的方案是否应该调整？如何选拔队员组成4100米混合泳接力队?问题二混合泳接力队的选拔甲乙丙丁戊蝶泳1’06”857”21’18”1’10”1’07”4仰泳1’15”61’06”1’07”81’14”21’11”蛙泳1’27”1’06”41’24”61’09”61’23”8自由泳58”653”59”457”21’02”45名候选人的百米成绩穷举法：组成接力队的方案共有5!=120种。目标函数若选择队员i参加泳姿j的比赛，记xij=1,否则记xij=00-1规划模型cij(秒)~队员i第j种泳姿的百米成绩约束条件每人最多入选泳姿之一ciji=1i=2i=3i=4i=5j=166.857.2787067.4j=275.66667.874.271j=38766.484.669.683.8j=458.65359.457.262.44151jiijijxcZMin每种泳姿有且只有1人5,1,141ixjij4,1,151jxiij模型求解最优解：x14=x21=x32=x43=1,其它变量为0;成绩为253.2(秒)=4’13”2MIN66.8x11+75.6x12+87x13+58.6x14+……+67.4x51+71x52+83.8x53+62.4x54SUBJECTTOx11+x12+x13+x14=1……x41+x42+x43+x44=1x11+x21+x31+x41+x51=1……x14+x24+x34+x44+x54=1ENDINT20输入LINDO求解甲乙丙丁戊蝶泳1’06”857”21’18”1’10”1’07”4仰泳1’15”61’06”1’07”81’14”21’11”蛙泳1’27”1’06”41’24”61’09”61’23”8自由泳58”653”59”457”21’02”4甲~自由泳