数学建模系列-常用模型课件

数学建模常用模型数学建模常用模型模型Ⅰ：层次分析法模型Ⅰ：层次分析法问题1选择旅游地现有三个旅游胜地可供选择，分别为苏杭、黄山、桂林，下面将作出旅游地的选择。问题1选择旅游地现有三个旅游胜地可供选择，分别为

面临各种各样的方案，要进行比较、判断、评价、最后作出决策。这个过程主观因素占有相当的比重给用数学方法解决问题带来不便。T.L.saaty等人20世纪在七十年代提出了一种能有效处理这类问题的实用方法。

层次分析法（AnalyticHierarchyProcess,AHP)这是一种定性和定量相结合的、系统化的、层次化的分析方法。过去研究自然和社会现象主要有机理分析法和统计分析法两种方法，前者用经典的数学工具分析现象的因果关系，后者以随机数学为工具，通过大量的观察数据寻求统计规律。近年发展的系统分析是又一种方法，而层次分析法是系统分析的数学工具之一。面临各种各样的方案，要进行比较、判断、评价、最后层次分析法（AHP）具体步骤：明确问题递阶层次结构的建立建立两两比较的判断矩阵层次单排序层次综合排序层次分析法（AHP）具体步骤：层次分析法的基本步骤1建立层次结构模型一般分为三层，最上面为目标层，最下面为方案层，中间是准则层或指标层。若上层的每个因素都支配着下一层的所有因素，或被下一层所有因素影响，称为完全层次结构，否则称为不完全层次结构。层次分析法的基本步骤1建立层次结构模型若上层的每选择旅游地景色费用居住饮食旅途苏杭、黄山、桂林建立选择旅游地层次结构

准则层A

方案层B目标层Z选择景色费用居住饮食旅途苏杭、建立选择旅游地层次结构准则分别分别表示景色、费用、居住、饮食、旅途。分别表示苏杭、黄山、桂林。分别分别表示景色、费用、分别表示苏杭、黄山、桂林。设某层有个因素，2构造成对比较矩阵要比较它们对上一层某一准则（或目标）的影响程度，确定在该层中相对于某一准则所占的比重。（即把个因素对上层某一目标的影响程度排序）用表示第个因素相对于第个因素的比较结果，则则称为成对比较矩阵。上述比较是两两因素之间进行的比较，比较时取1~9尺度。设某层有个因素，2构造成对比较矩阵要比较它们对上一层某标度定义与说明1两个元素对某个属性具有同样重要性3两个元素比较，一元素比另一元素稍微重要5两个元素比较，一元素比另一元素明显重要7两个元素比较，一元素比另一元素重要得多9两个元素比较，一元素比另一元素极端重要2,4,6,8表示需要在上述两个标准之间拆衷时的标度1/bij两个元素的反比较标度定义与说明1两个元素对某个属性具有同样重要性3两个元旅游问题中，第二层A的各因素对目标层Z的影响两两比较结果如下：ZA1A2A3A4A5A1A2A3A4A511/2433217551/41/711/21/31/31/52111/31/5311分别表示景色、费用、居住、饮食、旅途。旅游问题中，第二层A的各因素对目标层Z的影响两两比较结由上表，可得成对比较矩阵旅游问题的成对比较矩阵共有6个（一个5阶，5个3阶）。问题：两两进行比较后，怎样才能知道，下层各因素对上层某因素的影响程度的排序结果呢？由上表，可得成对比较矩阵旅游问题的成对比较矩阵共有6个（一个3层次单排序及一致性检验层次单排序：确定下层各因素对上层某因素影响程度的过程。用权值表示影响程度，先从一个简单的例子看如何确定权值。例如一块石头重量记为1，打碎分成各小块，各块的重量分别记为：则可得成对比较矩阵由右面矩阵可以看出，3层次单排序及一致性检验层次单排序：确定下层各因素即，但在例2的成对比较矩阵中，在正互反矩阵中，若则称为一致阵。一致阵的性质：4.的任一列(行)都是对应于特征根的特征向量。即，但在例2的成对比较矩阵中，在正互反矩阵中，若若成对比较矩阵是一致阵，则我们自然会取对应于最大特征根的归一化特征向量定理：阶互反阵的最大特征根，当且仅当时，为一致阵。表示下层第个因素对上层某因素影响程度的权值。

若成对比较矩阵不是一致阵，Saaty等人建议用其最大特征根对应的归一化特征向量作为权向量，则这样确定权向量的方法称为特征根法.若成对比较矩阵是一致阵，则我们自然会取对应于最定理：由于连续的依赖于，则比大的越多，的不一致性越严重。用最大特征值对应的特征向量作为被比较因素对上层某因素影响程度的权向量，其不一致程度越大，引起的判断误差越大。因而可以用数值的大小来衡量的不一致程度。定义一致性指标其中为的对角线元素之和，也为的特征根之和。由于连续的依赖于，则比大的越多，则可得一致性指标定义随机一致性指标随机构造500个成对比较矩阵随机一致性指标RI的数值：n1234567891011RI000.580.901.121.241.321.411.451.491.51则可得一致性指标一致性检验：利用一致性指标和一致性比率<0.1及随机一致性指标的数值表，对进行检验的过程。

一般，当一致性比率的不一致程度在容许范围之内，可用其归一化特征向量作为权向量，否则要重新构造成对比较矩阵，对加以调整。时，认为一致性检验：利用一致性指标和一致性比率<0.1一般，当一致4层次总排序及其一致性检验确定某层所有因素对于总目标相对重要性的排序权值过程，称为层次总排序从最高层到最低层逐层进行。设：对总目标Z的排序为的层次单排序为4层次总排序及其一致性检验对总目标Z的排序为的层次单排序为即层第个因素对总目标的权值为：层的层次总排序为：B层的层次总排序AB即层第个因素对层的层次总排序为：B层的层次AB层次总排序的一致性检验设层对上层(层)中因素的层次单排序一致性指标为，随机一致性指为，则层次总排序的一致性比率为：当时，认为层次总排序通过一致性检验。到此，根据最下层（决策层）的层次总排序做出最后决策。层次总排序的一致性检验设层对上层1.建立层次结构模型该结构图包括目标层，准则层，方案层。层次分析法的基本步骤归纳如下3.计算单排序权向量并做一致性检验2.构造成对比较矩阵从第二层开始用成对比较矩阵和1~9尺度。对每个成对比较矩阵计算最大特征值及其对应的特征向量，利用一致性指标、随机一致性指标和一致性比率做一致性检验。若检验通过，特征向量（归一化后）即为权向量；若不通过，需要重新构造成对比较矩阵。1.建立层次结构模型层次分析法的基本步骤归纳如下3.计算单排计算最下层对最上层总排序的权向量。4.计算总排序权向量并做一致性检验进行检验。若通过，则可按照总排序权向量表示的结果进行决策，否则需要重新考虑模型或重新构造那些一致性比率较大的成对比较矩阵。利用总排序一致性比率计算最下层对最上层总排序的权向量。4.计算总排序权向量并做一三层次分析法建模举例旅游问题(1)建模分别分别表示景色、费用、居住、饮食、旅途。分别表示苏杭、黄山、桂林。三层次分析法建模举例旅游问题分别分别表示景色、费用、分别（2）构造成对比较矩阵（2）构造成对比较矩阵(3)计算层次单排序的权向量和一致性检验成对比较矩阵的最大特征值表明通过了一致性验证。故则该特征值对应的归一化特征向量(3)计算层次单排序的权向量和一致性检验成对比较矩阵

对成对比较矩阵可以求层次总排序的权向量并进行一致性检验，结果如下：计算可知通过一致性检验。对成对比较矩阵对总目标的权值为：（4）计算层次总排序权值和一致性检验又决策层对总目标的权向量为：同理得，对总目标的权值分别为：故，层次总排序通过一致性检验。对总目标的权值为：（4）计算层次总排序权值和一致性检验又决策可作为最后的决策依据。故最后的决策应为去桂林。又分别表示苏杭、黄山、桂林，即各方案的权重排序为可作为最后的决策依据。故最后的决策应为去桂林。又模型Ⅱ：线性规划模型Ⅱ：线性规划丁的蛙泳成绩退步到1’15”2；戊的自由泳成绩进步到57”5,组成接力队的方案是否应该调整？如何选拔队员组成4100米混合泳接力队?问题二混合泳接力队的选拔

甲乙丙丁戊蝶泳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种。丁的蛙泳成绩退步到1’15”2；戊的自由泳成绩进步到57”5目标函数若选择队员i参加泳姿j的比赛，记xij=1,否则记xij=0

0-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每种泳姿有且只有1人目标函数若选择队员i参加泳姿j的比赛，记xij=1,否则模型求解最优解：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甲~自由泳、乙~蝶泳、丙~仰泳、丁~蛙泳.模型求解最优解：x14=x21=x32=x43丁蛙泳c43

=69.675.2，戊自由泳c54=62.4

57.5,方案是否调整？乙~蝶泳、丙~仰泳、丁~蛙泳、戊~自由泳最优解：x21=x32=x43=x51=1,成绩为4’17”7c43,c54的新数据重新输入模型，用LINDO求解指派(Assignment)问题：每项任务有且只有一人承担，每人只能承担一项，效益不同，怎样分派使总效益最大.讨论甲~自由泳、乙~蝶泳、丙~仰泳、丁~蛙泳.原方案丁蛙泳c43=69.675.2，戊自由泳c54=62.4模型Ⅲ：微分方程模型（一）Malthus模型（三）传染病模型（房室模型）（二）Logistic模型模型Ⅲ：微分方程模型（一）Malthus模型（三）传染病模型模型1马尔萨斯（Malthus）模型马尔萨斯在分析人口出生与死亡情况的资料后发现，人口净增长率r基本上是一常数，（r=b-d,b为出生率，d为死亡率），即：

或（1）

（2）

（1）的解为：其中N0=N(t0)为初始时刻t0时的种群数。马尔萨斯模型的一个显著特点：种群数量翻一番所需的时间是固定的。令种群数量翻一番所需的时间为T，则有：故模型1马尔萨斯（Malthus）模型马尔萨斯在分模型检验比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如，1961年世界人口数为30.6（即3.06×109），人口增长率约为2%，人口数大约每35年增加一倍。检查1700年至1961的260年人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每34.6年增加一倍，两者也几乎相同。模型预测假如人口数真能保持每34.6年增加一倍，那么人口数将以几何级数的方式增长。例如，到2510年，人口达2×1014个，即使海洋全部变成陆地，每人也只有9.3平方英尺的活动范围，而到2670年，人口达36×1015个，只好一个人站在另一人的肩上排成二层了。故马尔萨斯模型是不完善的。几何级数的增长Malthus模型实际上只有在群体总数不太大时才合理，到总数增大时，生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。所以Malthus模型假设的人口净增长率不可能始终保持常数，它应当与人口数量有关。模型检验比较历年的人口统计资料，可发现人模型2Logistic模型人口净增长率应当与人口数量有关，即：r=r(N)

从而有：（3）r(N)是未知函数，但根据实际背景，它无法用拟合方法来求。为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。r(N)最简单的形式是常数，此时得到的就是马尔萨斯模型。对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）对马尔萨斯模型引入一次项（竞争项），令r(N)=r-aN

此时得到微分方程：或（4）（4）可改写成：

（5）

(5)式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为K（近似地将K看成常数），N表示当前的种群数量，K-N恰为环境还能供养的种群数量，（5）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（5）也被称为统计筹算律的原因。模型2Logistic模型人口净增长率应当与人口数量有对（5）分离变量：两边积分并整理得：令N(0)=N0，求得：故（5）的满足初始条件N(0)=N0的解为：（6）易见：N(0)=N0

，N(t)的图形请看图对（5）分离变量：两边积分并整理得：令N(0)=N0，求得模型检验

用Logistic模型来描述种群增长的规律效果如何呢？1945年克朗皮克（Crombic）做了一个人工饲养小谷虫的实验，数学生物学家高斯（E·F·Gauss）也做了一个原生物草履虫实验，实验结果都和Logistic曲线十分吻合。

大量实验资料表明用Logistic模型来描述种群的增长，效果还是相当不错的。例如，高斯把5只草履虫放进一个盛有0.5cm3营养液的小试管，他发现，开始时草履虫以每天230.9%的速率增长，此后增长速度不断减慢，到第五天达到最大量375个，实验数据与r=2.309，a=0.006157，N(0)=5的Logistic曲线：

几乎完全吻合，见图模型检验用Logistic模型来描述种群增长的规律效Malthus模型和Logistic模型的总结

Malthus模型和Logistic模型均为对微分方程所作的模拟近似方程。前一模型假设了种群增长率r为一常数，（r被称为该种群的内禀增长率）。后一模型则假设环境只能供养一定数量的种群，从而引入了一个竞争项。

用模拟近似法建立微分方程来研究实际问题时必须对求得的解进行检验，看其是否与实际情况相符或基本相符。相符性越好则模拟得越好，否则就得找出不相符的主要原因，对模型进行修改。Malthus模型与Logistic模型虽然都是为了研究种群数量的增长情况而建立的，但它们也可用来研究其他实际问题，只要这些实际问题的数学模型有相同的微分方程即可。Malthus模型和Logistic模型的总结

模型三

传染病模型传染病是人类的大敌，通过疾病传播过程中若干重要因素之间的联系建立微分方程加以讨论，研究传染病流行的规律并找出控制疾病流行的方法显然是一件十分有意义的工作。在本节中，我们将主要用多房室系统的观点来看待传染病的流行，并建立起相应的多房室模型。医生们发现，在一个民族或地区，当某种传染病流传时，波及到的总人数大体上保持为一个常数。即既非所有人都会得病也非毫无规律，两次流行（同种疾病）的波及人数不会相差太大。如何解释这一现象呢？试用建模方法来加以证明。问题的提出：模型三传染病模型传染病是人类的大敌，通过设某地区共有n+1人，最初时刻共有i人得病，t时刻已感染（infective）的病人数为i(t)，假定每一已感染者在单位时间内将疾病传播给k个人（k称为该疾病的传染强度），且设此疾病既不导致死亡也不会康复模型1此模型即Malthus模型，它大体上反映了传染病流行初期的病人增长情况，在医学上有一定的参考价值，但随着时间的推移，将越来越偏离实际情况。已感染者与尚未感染者之间存在着明显的区别，有必要将人群划分成已感染者与尚未感染的易感染，对每一类中的个体则不加任何区分，来建立两房室系统。则可导出：故可得：

（3.1）

设某地区共有n+1人，最初时刻共有i人得病模型2记t时刻的病人数与易感染人数（susceptible）分别为i(t)与s(t)，初始时刻的病人数为i。根据病人不死也不会康复的假设及（竞争项）统计筹算律，其中：解得：（3.17）可得：（3.16）统计结果显示，(3.17)预报结果比(3.15)更接近实际情况。医学上称曲线为传染病曲线，并称最大值时刻t1为此传染病的流行高峰。令：得：此值与传染病的实际高峰期非常接近，可用作医学上的预报公式。

模型2仍有不足之处，它无法解释医生们发现的现象，且当时间趋与无穷时，模型预测最终所有人都得病，与实际情况不符。为了使模型更精确，有必要再将人群细分，建立多房室系统模型2记t时刻的病人数与易感染人数（susinfectiverecoveredsusceptiblekl

（3.18）

l称为传染病恢复系数求解过程如下：对（3）式求导，由（1）、（2）得：解得：记：则：将人群划分为三类（见右图）：易感染者、已感染者和已恢复者（recovered）。分别记t时刻的三类人数为s(t)、i(t)和r(t)，则可建立下面的三房室模型：模型3infectiverecoveredsusceptibl