1、如何建立数学模型及实例数学建模培训科研处数学建模小组第五章：如何建立数学模型怎样撰写数学建模的论文？1什么是数学模型？数学模型是对于现实世界的一个特定对象，一个特定目的，根据特有的内在规律，做出一些必要的假设，运用适当的数学工具，得到一个数学结构。简单地说：就是系统的某种特征的本质的数学表达式（或是用数学术语对部分现实世界的描述） ，即用数学式子（如函数、图形、代数方程、微分方程、积分方程、差分方程等）来描述（表述、模拟）所研究的客观对象或系统在某一方面的存在规律。2什么是数学建模? 数学建模是利用数学方法解决实际问题的一种实践。即通过抽象、简化、假设、引进变量等处理过程后，将实际问题用数学方
2、式表达，建立起数学模型，然后运用先进的数学方法及计算机技术进行求解。观点：“所谓高科技就是一种数学技术”注 数学建模其实并不是什么新东西，可以说有了数学并需要用数学去解决实际问题，就一定要用数学的语言、方法去近似地刻划该实际问题，这种刻划的数学表述的就是一个数学模型，其过程就是数学建模的过程。数学模型一经提出，就要用一定的技术手段（计算、证明等）来求解并验证，其中大量的计算往往是必不可少的，高性能的计算机的出现使数学建模这一方法如虎添翼似的得到了飞速的发展，掀起一个高潮。注 数学建模 将各种知识综合应用于解决实际问题中，是培养和提高同学们应用所学知识分析问题、解决问题的能力的必备手段之一。3数
3、学建模的一般方法和步骤建立数学模型的方法和步骤并没有一定的模式，但一个理想的模型应能反映系统的全部重要特征：模型的可靠性和模型的使用性建模的一般方法： 机理分析 测试分析方法机理分析：根据对现实对象特性的认识，分析其因果关系，找出反映内部机理的规律，所建立的模型常有明确的物理或现实意义。测试分析方法：将研究对象视为一个“黑箱”系统，内部机理无法直接寻求，通过测量系统的输入输出数据，并以此为基础运用统计分析方法，按照事先确定的准则在某一类模型中选出一个数据拟合得最好的模型。 测试分析方法也叫做系统辩识。将这两种方法结合起来使用，即用机理分析方法建立模型的结构，用系统测试方法来确定模型的参数，也是
4、常用的建模方法。 在实际过程中用那一种方法建模主要是根据我们对研究对象的了解程度和建模目的来决定。机理分析法建模的具体步骤大致可见下图。4、数学模型及其分类 模型 抽 象 模 型具 体 模 型 数 学 模 型符 号 模 型思 维 模 型物 理 模 型直 观 模 型 图 形 模 型数 式 模 型数学模型的分类： 按研究方法和对象的数学特征分：初等模型、几何模型、优化模型、微分方程模型、符合实际不符合实际交付使用，从而可产生经济、社会效益实际问题抽象、简化、假设确定变量、参数建立数学模型并数学、数值地求解、确定参数用实际问题的实测数据等来检验该数学模型图论模型、逻辑模型、稳定性模型、扩散模型等。
5、按研究对象的实际领域（或所属学科）分：人口模型、交通模型、环境模型、生态模型、生理模型、城镇规划模型、水资源模型、污染模型、经济模型、社会模型等。5怎样撰写数学建模的论文？1、 摘要: 摘要是论文内容不加注释和评论的简短陈述，其作用是使读者不阅读论文全文即能获得必要的信息。在数学建模论文中，摘要是非常重要的一部分。数学建模论文的摘要应包含以下内容：所研究的实际问题、建立的模型、求解模型的方法、获得基本结果以及对模型的检验或推广。论文摘要需要用概括、简练的语言反映这些内容，尤其要突出论文的优点，如巧妙的建模方法、快速有效的算法、合理的推广等。一般科技论文的摘要要求不列举例证，不出现图、表和数学公
6、式，不自我评价，且字数应在 200 以内。前几年，全国大学生数学建模竞赛要求摘要字数应在 300 字以内。但从 2001 年开始，为了提高论文评选效率，要求将论文第一页全用作摘要，对字数已无明确限制。故在摘要中也可适当出现反映结果的图、表和数学公式。 2、 问题重述： 数学建模比赛要求解决给定的问题，所以论文中应叙述给定问题。撰写这部分内容时，不要照抄原题，应把握住问题的实质，再用较精练的语言叙述问题。3、 模型假设： 建模时，要根据问题的特征和建模的目的，抓住问题的本质，忽略次要因素，对问题进行必要的简化，做出一些合理的假设。模型假设部分要求用精练、准确的语言列出问题中所给出的假设，以及为了
7、解决问题所做的必要、合理的假设。假设作得不合理或太简单，会导致错误的或无用的模型；假设作得过分详尽，试图把复杂的众多因素都考虑进去，会使工作很难或无法继续下去，因此常常需要在合理与简化之间作出恰当的折中。4、 分析与建立模型： 根据假设，用数学的语言、符号描述对象的内在规律，得到一个数学结构。建模时就尽量采用简单的数学工具，使建立的模型易于被人理解。在撰写这一部分时，对所用的变量、符号、计量单位应作解释，特定的变量和参数应在整篇文章保持一致。为使模型易懂，可借助于适当图形、表格来描述问题或数据。5、 模型求解： 使用各种数学方法或软件包求解数学模型。此部分应包括求解过程的公式推导、算法步骤及计
8、算结果。为求解而编写的计算机程序应放在附录部分。有时需要对结果进行数学上的分析，如结果的误差分析、模型对数据的稳定性或灵敏度分析等。6、 模型检验： 指导求解和分析结果翻译回到实际问题，与实际的现象、数据比较，检验模型的合理性和适用性。如果结果与实际不符，问题常出现在模型假设上，应该修改、补充假设，重新建模。这一步对于模型是否真的有用十分关健。7、 模型推广： 将该问题的模型推广到解决更多的类似问题，或讨论给出该模型的更一般情况下的解法，或指出可能的深化、推广及进一步研究的建议。8、参考文献： 在正文中提及或直接引用的材料或原始数据，应注明出处，并将相应的出版物列举在参考文献中。需标明出版物名
9、称、页码、著者姓名、出版日期、出版单位等、9、附录： 附录是正文的补充，与正文有关而又不方便于编入正文的内容都收集在这里。包括：计算机程序、比较重要但数据量较大的中间结果等。为便于阅读，应在源程序中加入足够多的注释和说明语句。附：2006 高教社杯全国大学生数学建模竞赛论文格式规范 甲组参赛队从 A、B 题中任选一题，乙组参赛队从 C、D 题中任选一题。 论文（答卷）用白色 A4 纸，上下左右各留出 2.5 厘米的页边距。 论文第一页为承诺书，具体内容和格式见本规范第二页。 论文第二页为编号专用页，用于赛区和全国评阅前后对论文进行编号，具体内容和格式见本规范第三页。 论文题目和摘要写在论文第三
10、页上，从第四页开始是论文正文。 论文从第三页开始编写页码，页码必须位于每页页脚中部，用阿拉伯数字从“1”开始连续编号。 论文不能有页眉，论文中不能有任何可能显示答题人身份的标志。 论文题目用 3 号黑体字、一级标题用 4 号黑体字，并居中。论文中其他汉字一律采用小 4 号黑色宋体字，行距用单倍行距。 提请大家注意：摘要在整篇论文评阅中占有重要权重，请认真书写摘要（注意篇幅不能超过一页）。全国评阅时将首先根据摘要和论文整体结构及概貌对论文优劣进行初步筛选。 引用别人的成果或其他公开的资料(包括网上查到的资料) 必须按照规定的参考文献的表述方式在正文引用处和参考文献中均明确列出。正文引用处用方括号
11、标示参考文献的编号，如13等；引用书籍还必须指出页码。参考文献按正文中的引用次序列出，其中书籍的表述方式为：编号 作者，书名，出版地：出版社，出版年。参考文献中期刊杂志论文的表述方式为：编号 作者，论文名，杂志名，卷期号：起止页码，出版年。参考文献中网上资源的表述方式为：编号 作者，资源标题，网址，访问时间（年月日）。 在不违反本规范的前提下，各赛区可以对论文增加其他要求（如在本规范要求的第一页前增加其他页和其他信息，或在论文的最后增加空白页等） ；从承诺书开始到论文正文结束前，各赛区不得有本规范外的其他要求（否则一律无效） 。 本规范的解释权属于全国大学生数学建模竞赛组委会。注 赛区评阅前将
12、论文第一页取下保存，同时在第一页和第二页建立“赛区评阅编号” （由各赛区规定编号方式） ， “赛区评阅纪录”表格可供赛区评阅时使用（各赛区自行决定是否在评阅时使用该表格） 。评阅后，赛区对送全国评阅的论文在第二页建立“全国统一编号”（编号方式由全国组委会规定，与去年格式相同） ，然后送全国评阅。论文第二页（编号页）由全国组委会评阅前取下保存，同时在第二页建立“全国评阅编号” 。全国大学生数学建模竞赛组委会2006 年 5 月修订2006 高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、
13、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D 中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名) ： 日期： 年 月 日赛区
14、评阅编号（由赛区组委会评阅前进行编号）：2006 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：第六章几个例子（一）数学建模实例：人口预报问题1.问题人口问题是当前世界上人们最关心的问题之一.认识人口数量的变化规律，作出较准确的预报，是有效控制人口增长的前提.下面介绍两个最基本的人口模型，并利用表 1 给出的近两百年的美国人口统计数据，对模型做出检验，最后用它预报 2000 年、2010 年美国人口.表 1 美国
15、人口统计数据年（公元）人口（百万）17903.918005.318107.218209.6183012.9184017.1185023.2年（公元）人口（百万）186031.4187038.6188050.2189062.9190076.0191092.01920106.5年（公元）人口（百万）1930123.21940131.71950150.71960179.31970204.01980226.51990251.42.指数增长模型（马尔萨斯人口模型）此模型由英国人口学家马尔萨斯（Malthus17661834）于 1798 年提出.1 假设：人口增长率 r 是常数（或单位时间内人口的增长量
16、与当时的人口成正比）.2 建立模型： 记时刻 t=0 时人口数为 x0, 时刻 t 的人口为 ，由于量大， 可视txtx为连续、可微函数.t 到 时间内人口的增量为：ttrxtx于是 满足微分方程：t（1）00xxrdt3 模型求解： 解微分方程（1）得（2）rtext0表明： 时， （r0）.tx4 模型的参数估计：要用模型的结果（2）来预报人口，必须对其中的参数 r进行估计，这可以用表 1的数据通过拟合得到.拟合的具体方法见本书第 16章或第 18章.通过表中 1790-1980的数据拟合得：r=0.307. 5 模型检验：将 x0=3.9，r=0.307 代入公式（2） ，求出用指数增长
17、模型预测的 1810-1920的人口数，见表 2.表 2 美国实际人口与按指数增长模型计算的人口比较指数增长模型年实际人口（百万） 预测人口（百万）误差（%）1790 3.91800 5.31810 7.2 7.3 1.41820 9.6 10.0 4.21830 12.9 13.7 6.21840 17.1 18.7 9.41850 23.2 25.6 10.31860 31.4 35.0 10.81870 38.6 47.8 23.81880 50.2 65.5 30.51890 62.9 89.6 42.41900 76.0 122.5 61.21910 92.0 167.6 82.11
18、920 106.5 229.3 115.3从表 2可看出，1810-1870 间的预测人口数与实际人口数吻合较好，但 1880年以后的误差越来越大.分析原因，该模型的结果说明人口将以指数规律无限增长.而事实上，随着人口的增加，自然资源、环境条件等因素对人口增长的限制作用越来越显著.如果当人口较少时人口的自然增长率可以看作常数的话，那么当人口增加到一定数量以后，这个增长率就要随着人口增加而减少.于是应该对指数增长模型关于人口净增长率是常数的假设进行修改.下面的模型是在修改的模型中著名的一个. 3. 阻滞增长模型（Logistic 模型）1假设：（a）人口增长率 r为人口 的函数 （减函数） ，最
19、简单假定txxr（线性函数） ，r 叫做固有增长率.0, sxr（b）自然资源和环境条件年容纳的最大人口容量 .m2建立模型：当 时，增长率应为 0，即 =0，于是 ，代入 得：mxxrmxrssxr（3）mxrx1将（3）式代入（1）得：模型为： （4） 01xrdtm3 模型的求解： 解方程组（4）得 （5）rtmext10根据方程（4）作出 曲线图，见图 1-1，由该图可看出人口增长率随人口xdt数的变化规律.根据结果（5）作出 xt曲线，见图 1-2，由该图可看出人口数随时间的变化规律. o xm2mxdtx图 1-1 曲线图 dtxm20xto图 1-2 xt 曲线图4 模型的参数估
20、计：利用表 1中 1790-1980的数据对 r和 xm拟合得：r=0.2072, x m=464.5 模型检验：将 r=0.2072, xm=464代入公式（5） ，求出用指数增长模型预测的 1800-1990的人口数，见表 3第 3、4 列.也可将方程（4）离散化，得t=0,1,2, (6)(1()(1( txrtxttxm用公式（6）预测 1800-1990的人口数，结果见表 3第 5、6 列.表 3 美国实际人口与按阻滞增长模型计算的人口比较阻滞增长模型公式（5） 公式（6）年实际人口（百万）预测人口（百万）误差（%）预测人口（百万）误差（%）1790 3.91800 5.3 5.90
21、25 0.1137 3.9000 0.26421810 7.2 7.2614 0.0085 6.5074 0.09621820 9.6 8.9332 0.0695 8.6810 0.09571830 12.9 10.9899 0.1481 11.4153 0.11511840 17.1 13.5201 0.2094 15.1232 0.11561850 23.2 16.6328 0.2831 19.8197 0.14571860 31.4 20.4621 0.3483 26.5228 0.15531870 38.6 25.1731 0.3478 35.4528 0.08151880 50.2
22、30.9687 0.3831 43.5329 0.13281890 62.9 38.0986 0.3943 56.1884 0.10671900 76.0 46.8699 0.3833 70.1459 0.07701910 92.0 57.6607 0.3733 84.7305 0.07901920 106.5 70.9359 0.3339 102.4626 0.03791930 123.2 0.2917 118.9509 0.034587.26741940 131.7 107.3588 0.1848 137.8810 0.04691950 150.7 132.0759 0.1236 148.
23、7978 0.01261960 179.3 162.4835 0.0938 170.2765 0.05031970 204.0 199.8919 0.0201 201.1772 0.01381980 226.5 245.9127 0.0857 227.5748 0.00471990 251.4 302.5288 0.2034 250.4488 0.00386 模型应用：现应用该模型预测人口.用表 1中 1790-1990年的全部数据重新估计参数，可得r=0.2083, xm=457.6. 用公式（6）作预测得：x(2000)=275; x(2010)=297.9.也可用公式（5）进行预测.（二
24、）椅子能在不平的地面上放稳吗？把椅子往不平的地面上一放，通常只有三只脚着地，放不稳，然而只要稍挪动几次，就可以四脚着地，放稳了。下面用数学语言证明。一、 模型假设对椅子和地面都要作一些必要的假设：1、椅子四条腿一样长，椅脚与地面接触可视为一个点，四脚的连线呈正方形。2、地面高度是连续变化的，沿任何方向都不会出现间断（没有像台阶那样的情况） ，即地面可视为数学上的连续曲面。3、对于椅脚的间距和椅脚的长度而言，地面是相对平坦的，使椅子在任何位置至少有三只脚同时着地。二、模型建立中心问题是数学语言表示四只脚同时着地的条件、结论。首先用变量表示椅子的位置，由于椅脚的连线呈正方形，以中心为对称点，正方形
25、绕中心的旋转正好代表了椅子的位置的改变，于是可以用旋转角度 这一变量来表示椅子的位置。BAC A xDD其次要把椅脚着地用数学符号表示出来，如果用某个变量表示椅脚与地面的竖直距离，当这个距离为 0 时，表示椅脚着地了。椅子要挪动位置说明这个距离是位置变量的函数。由于正方形的中心对称性，只要设两个距离函数就行了，记 A、C 两脚与地面距离之和为 ，B、D 两脚与地面距离之和为 ，显然 、 ，由假设 2 知 f、g 都f gf0g是连续函数，再由假设 3 知 、 至少有一个为 0。当 时，不妨设f ，这样改变椅子的位置使四只脚同时着地，就归结为如下命题：0,fg命题 已知 、 是 的连续函数，对任
26、意 ， * =0，且f f，则存在 ，使 。 ,0f00fg三、模型求解将椅子旋转 ，对角线 AC 和 BD 互换，由 可知09 0,f。令 ，则 ，由 f、g 的连2,fgfgh2h续性知 h 也是连续函数，由零点定理，必存在 使 ，00，由 ，所以 。0f*0fgf四、评 注模型巧妙在于用已元变量 表示椅子的位置，用 的两个函数表示椅子四脚与地面的距离。利用正方形的中心对称性及旋转 并不是本质的，同学们可以考虑四脚呈长方形09的情形。（三）双层玻璃的功效北方城镇的有些建筑物的窗户是双层的，即窗户上装两层厚度为 的玻璃夹着一层d厚度为 的空气，如左图所示，据说这样做是为了保暖，即减少室内向室
27、外的热量流失。l我们要建立一个模型来描述热量通过窗户的热传导（即流失）过程，并将双层玻璃窗与用同样多材料做成的单层玻璃窗（如右图，玻璃厚度为 ）的热量传导进行对比，2对双层玻璃窗能够减少多少热量损失给出定量分析结果。一、 模型假设1、 热量的传播过程只有传导，没有对流。即假定窗户的密封性能很好，两层玻璃之间的空气是不流动的；2、 室内温度 和室外温度 保持不变，热传导过程已处于稳定状态，即沿热传导1T2方向，单位时间通过单位面积的热量是常数；3、 玻璃材料均匀，热传导系数是常数。二、 符号说明室内温度1T室外温度2单层玻璃厚度d两层玻璃之间的空气厚度l内层玻璃的外侧温度aT外层玻璃的内侧温度b
28、热传导系数k热量损失Q三、 模型建立与求解由物理学知道，在上述假设下，热传导过程遵从下面的物理规律：厚度为 的均匀介质，两侧温度差为 ，则单位时间由温度高的一侧向温度低的dT一侧通过单位面积的热量为 ，与 成正比，与 成反比，即Qd（1）dk其中 为热传导系数。k1、双层玻璃的热量流失记双层窗内窗玻璃的外侧温度为 ，外层玻璃的内侧温度为 ，玻璃的热传导系数aTbT为 ，空气的热传导系数为 ，由（1）式单位时间单位面积的热量传导（热量流失）为：1 2k（2）dTkTkdkQbaa 2121 由 及 可得T1b1 121)(kQdba再代入 就将（2）中 、 消去，变形可得：dkba2aTb（3）
29、dlhkssTQ , ,)(2112、单层玻璃的热量流失对于厚度为 的单层玻璃窗户，容易写出热量流失为： d（4）Tk213、 单层玻璃窗和双层玻璃窗热量流失比较比较（3） （4）有： （5）2sQ显然， 。为了获得更具体的结果，我们需要 的数据，从有关资料可知，不流通、干燥空21,k气的热传导系数 （焦耳/厘米.秒.度） ，常用玻璃的热传导系数42105.k（焦耳/厘米.秒.度） ，于是331804k2621k在分析双层玻璃窗比单层玻璃窗可减少多少热量损失时，我们作最保守的估计，即取 ，由（ 3） （5）可得：162k（6）dlhQ 184、模型讨论比值 反映了双层玻璃窗在减少热量损失上的功
30、效，它只与 有关，下图 dlh给出了 的曲线，当 由 0 增加时， 迅速下降，而当 超过一定值（比如hQ）后 下降缓慢，可见 不宜选得过大。hh四、 模型的应用这个模型具有一定的应用价值。制作双层玻璃窗虽然工艺复杂会增加一些费用，但它减少的热量损失却是相当可观的。通常，建筑规范要求 。按照这个模型，4dlh，即双层玻璃窗比用同样多的玻璃材料制成的单层窗节约热量 97%左右。不难%3Q发现，之所以有如此高的功效主要是由于层间空气的极低的热传导系数 ，而这要求空气2k是干燥、不流通的。作为模型假设的这个条件在实际环境下当然不可能完全满足，所以实际上双层玻璃窗的功效会比上述结果差一些。（四）钢管订购
31、和运输优化模型要铺设一条 的输送天然气的主管道, 如图一所示(见反面)。经1521AA筛选后可以生产这种主管道钢管的钢厂有 。图中粗线表示铁路，单细线表示公72,S路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位 km)。为方便计，1km 主管道钢管称为 1 单位钢管。一个钢厂如果承担制造这种钢管，至少需要生产 500 个单位。钢厂 在指定期限内能iS生产该钢管的最大数量为 个单位，钢管出厂销价 1 单位钢管为 万元，如下表：is ipi1 2 3 4 5 6 7is800 800 1000 2000 2
32、000 2000 3000ip160 155 155 160 155 150 1601 单位钢管的铁路运价如下表：里程(km) 300 301350 351400 401450 451500运价(万元) 20 23 26 29 32里程(km) 501600 601700 701800 801900 9011000运价(万元) 37 44 50 55 601000km 以上每增加 1 至 100km 运价增加 5 万元。公路运输费用为 1 单位钢管每公里 0.1 万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点 ，而是管道全线）。1521,A问题：（1）请制定一
33、个主管道钢管的订购和运输计划，使总费用最小（给出总费用)。思考题：（2）请就（1）的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果。（3）如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对图二按（1）的要求给出模型和结果。A13258010103120124270 1088107062703020 2030450104 301750 606194 205201680480 300220 21042050060030601952027206
34、90520170690462160 320160110290115011001200A2A3A4A5A6A11A711A11A8A11A911A11A10A11A12A13A14A15S1S2S3S4S5S6S7图一一基本假设：1 沿铺设的主管道以有公路或者有施工公路。2 在主管道上，每公里卸 1 单位的钢管。3 公路运输费用为 1 单位钢管每公里 0.1 万元（不足整公里部分按整公里计算）4 在计算总费用时，只考虑运输费和购买钢管的费用，而不考虑其他费用。5 在计算钢厂的产量对购运计划影响时，只考虑钢厂的产量足够满足需要的情况，即钢厂的产量不受限制。6 假设钢管在铁路运输路程超过 1000k
35、m 时，铁路每增加 1 至 100km，1 单位钢管的运价增加 5 万元。二符号说明：第 个钢厂； iS 7,21i：第 个钢厂的最大产量； s ：输送管道（主管道）上的第 个点； jAj 5,jA132580101031201242 70 1088107062703020 2030450104 301750 606194 205201680480 3002202104205006003060195202720690520170690462160 320160110290115011001200A19130190260100A2A3A4A5A6A7A8A11A9 A10A11A12A13A14
36、A15S1S2S3S4S5S6S7A16 A17A18A20(A21)图二：第 个钢厂 1 单位钢管的销价； ip 7,21i：钢厂 向点 运输的钢管量； ijxiSjA 15,2j：在点 与点 之间的公路上，运输点 向点 方向铺设的钢管量；tj1 jA1j( )4,32,0t：1 单位钢管从钢厂 运到结点 的最少总费用，即公路运费铁路运费和ijaij钢管销价之和； 7,2 5,2j：与点 相连的公路和铁路的相交点； jbjA 13：相邻点 与 之间的距离； 1.j1j 4,j三模型的建立与求解问题一：讨论如何调整主管道钢管的订购和运输方案使总费用最小由题意可知，钢管从钢厂 到运输结点 的费用
37、 包括钢管的销价钢管的铁路运iSjAija输费用和钢管的公路运输费用。在费用 最小时，对钢管的订购和运输进行分配，可得出ija本问题的最佳方案。1、 求钢管从钢厂 运到运输点 的最小费用iSjA1）将图一转换为一系列以单位钢管的运输费用为权的赋权图。由于钢管从钢厂 运到运输点 要通过铁路和公路运输，而铁路运输费用是分段函i j数，与全程运输总距离有关。又由于钢厂 直接与铁路相连，所以可先求出钢厂 到铁路iSiS与公路相交点 的最短路径。如图三jb图三 铁路网络图依据钢管的铁路运价表，算出钢厂 到铁路与公路相交点 的最小铁路运输费用，并iSjb把费用作为边权赋给从钢厂 到 的边。再将与 相连的公
38、路、运输点 及其与之相连iSjbjbiA的要铺设管道的线路（也是公路）添加到图上，根据单位钢管在公路上的运价规定，得出每一段公路的运费，并把此费用作为边权赋给相应的边。以 为例得图四1S图四 钢管从钢厂 运到各运输点 的铁路运输与公路运输费用权值图1SjA2）计算单位钢管从 到 的最少运输费用1j根据图四，借助图论软件包中求最短路的方法求出单位钢管从 到 的最少运输费1SjA用依次为：170.7，160.3，140.2，98.6，38，20.5，3.1，21.2，64.2，92，96，106，121.2，128，142（单位：万元） 。加上单位钢管的销售价 ，得出从钢厂 购买单位钢管运输到点
39、的最ip1 j小费用 依次为：ja1330.3，320.3，300.2，258.6，198，180.5，163.1，181.2，224.2，252，256，266，281.2，288，302（单位：万元） 。同理，可用同样的方法求出钢厂 到点 的最小费用，从而得出2S3456S7jA钢厂到点的最小总费用（单位：万元）为：表一 到点 最小费用ijAa2 a3 a4 a5 a6 a7 a8 a9 a10 a11 a12 a13 a14 a15s1 320.3 300.2 258.6 198 180.5 163 181.2 224.2 252 256 266 281.2 288 302s2 360.
40、3 345.2 326.6 266 250.5 241 226.2 269.2 297 301 311 326.2 333 347s3 375.3 355.2 336.6 276 260.5 251 241.2 203.2 237 241 251 266.2 273 287s4 410.3 395.2 376.6 316 300.5 291 276.2 244.2 222 211 221 236.2 243 257s5 400.3 380.2 361.6 301 285.5 276 266.2 234.2 212 188 206 226.2 228 242s6 405.3 385.2 366.
41、6 306 290.5 281 271.2 234.2 212 201 195 176.2 161 178s7 425.3 405.2 386.6 326 310.5 301 291.2 259.2 237 226 216 198.2 186 1622、建立模型运输总费用可分为两部分：运输总费用=钢厂到各点的运输费用+铺设费用。运输费用：若运输点 向钢厂 订购 单位钢管，则钢管从钢厂 运到运输点jAiSijxiS所需的费用为 。由于钢管运到 必须经过 ，所以可不考虑 ，那么所有钢管jAijxa12A1A从各钢厂运到各运输点上的总费用为： 。527jiija铺设费用：当钢管从钢厂 运到点 后，钢
42、管就要向运输点 的两边 段和iSj j 1j段运输（铺设）管道。设 向 段铺设的管道长度为 ，则 向 段j1 jA1y的运输费用为 （万元） ；由于相邻运输点 与 之间的20)21(.0jjty jA1j距离为 ，那么 向 段铺设的管道长为 ，所对应的铺设费用为.jAj1j jjtA1.（万元） 。所以，主管道上的铺设费用为：201.1. jj tt 41 120j jjjjj tA总费用为： 715241 1.1.2020ijj jjjjji tAttaxf又因为一个钢厂如果承担制造钢管任务，至少需要生产 500 个单位，钢厂 在指定期iS限内最大生产量为 个单位，故 或 因此本问题可建立如
43、下isijisx52152jix的非线性规划模型：3、模型求解：由于 MATLAB 不能直接处理约束条件： 或 ，我们可先将此ijisx15200152jix 015,2,71 ,1523j . )(0)(min.1527 152714 1.1. jjijiiij jiijj jjjjAtjxsxts atAttf 0,2,71 0 5,1523j . 0)1)(0)(min.127 152714 .1. jji jijiiij jiijj jjjjjAtxxsxts axtAttf 或条件改为 ，得到如下模型：ijisx152用 MATLAB 求解，分析结果后发现购运方案中钢厂 的生产量不足
44、 500 单位，下面我7S们采用不让钢厂 生产和要求钢厂 的产量不小于 500 个单位两种方法计算：7S7S1）不让钢厂 生产计算结果： 1278632（万元） （此时每个钢厂的产量都满足条件） 。1f2）要求钢厂 的产量不小于 500 个单位7计算结果： 1279664 （万元） （此时每个钢厂的产量都满足条件） 。2f比较这两种情况，得最优解为， =1278632（万元） 具121),min(ifff体的购运计划如表二：表二 问题一的订购和调运方案订购量 A2 A3 A4 A5 A6 A7 A8 A9 A10 A11 A12 A13 A14 A15S1 800 0 201 133 200 266 0 0 0 0 0 0 0 0 0S2 800 179 11 14 295 0 0 300 0 0 0 0 0 0 0S3 1000 139 11 186 0 0 0 664 0 0 0 0 0 0 0S4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0S5 1015 0 358 242 0 0 0 0 0 0 415 0 0 0 0S6 1556 0 0 0 0 0 0 0 0 0 351 86 333 621 165S7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
