数学建模之运筹学

第一讲数学建模简介及数学规划模型IntroductionofMMandMathematicalProgrammingModelNetworkProgramming数学建模MathematicalModeling第一讲数学建模简介及数学规划模型Page2of68数学建模简介第一讲数学建模简介及数学规划模型Page3of68一般地，数学模型可以描述为，对于现实世界的一个特定对象，为了一个特定目的，根据特有的内在规律，作出一些必要的简化假设，运用适当的数学工具，得到的一个数学结构。把现实世界中的实际问题加以提炼，抽象为数学模型，求出模型的解，验证模型的合理性，并用该数学模型所提供的解答来解释现实问题，我们把数学知识的这一应用过程称为数学建模。数学模型或者能解释特定现象的现实状态，或者能预测到对象的未来状况，或者能提供处理对象的最优决策或控制。第一讲数学建模简介及数学规划模型Page4of68数学模型的分类1、按模型的应用领域分类：生物数学模型医学数学模型地质数学模型数量经济学模型数学社会学模型2、按是否考虑随机因素分类：确定性模型随机性模型3、按是否考虑模型的变化分类：静态模型动态模型第一讲数学建模简介及数学规划模型Page5of684、按应用离散方法或连续方法分类：离散模型连续模型5、按建立模型的数学方法分类：几何模型微分方程模型图论模型规划论模型马氏链模型第一讲数学建模简介及数学规划模型Page6of686、按人们对是物发展过程的了解程度分类：（1）白箱模型：指那些内部规律比较清楚的模型。如力学、热学、电学以及相关的工程技术问题。（2）灰箱模型：指那些内部规律尚不十分清楚，在建立和改善模型方面都还不同程度地有许多工作要做的问题。如气象学、生态学经济学等领域的模型。（3）黑箱模型：指一些其内部规律还很少为人们所知的现象。如生命科学、社会科学等方面的问题。但由于因素众多、关系复杂，也可简化为灰箱模型来研究。第一讲数学建模简介及数学规划模型Page7of68数学建模的几个过程1、模型准备2、模型假设3、模型建立4、模型构成5、模型求解6、模型分析7、模型检验8、模型应用第一讲数学建模简介及数学规划模型Page8of68模型准备了解实际背景明确建模目的搜集有关信息掌握对象特征形成一个比较清晰的‘问题’模型假设针对问题特点和建模目的作出合理的、简化的假设在合理与简化之间作出折中第一讲数学建模简介及数学规划模型Page9of68模型建立用数学的语言、符号描述问题发挥想像力使用类比法尽量采用简单的数学工具各种数学方法、软件和计算机技术如结果的误差分析、统计分析、模型对数据的稳定性分析模型求解模型分析第一讲数学建模简介及数学规划模型Page10of68模型检验与实际现象、数据比较，检验模型的合理性、适用性模型应用第一讲数学建模简介及数学规划模型Page11of68数学建模有助于培养以下几个方面的素质和能力：•数学素质和能力•计算机应用能力•论文写作能力•团队合作精神和进行协调的组织能力•培养想象能力•发展观察力，形成洞察力•勇于参与的竞争意识和不怕困难、奋力攻关的顽强意志第一讲数学建模简介及数学规划模型Page12of68为培养和选拔优秀的数学人才，世界各国有各种不同形式不同层次的数学竞赛.传统的数学竞赛只局限于演绎、推理等纯数学形式，它不能培养和发展学生运用数学知识解决实际问题的能力，不能满足科学技术飞速发展的时代需要.从1983年起，在美国就有一些有识之士开始探讨组织一项应用数学方面的竞赛的可能性.第一讲数学建模简介及数学规划模型Page13of681985年美国第一届大学生数学建模竞赛（mathematicalcompetitioninmodeling）1988年改为mathematicalcontestinmodeling简称MCM.由美国工业与应用数学会和美国运筹学会联合举办.1985年起每年举行一届，一般在每年的二月下旬或三月初的某个星期五或星期日举行.美国竞赛评出Outstanding,Meritorious,HonorableMention及SuccessfulParticipation等级别.第一讲数学建模简介及数学规划模型Page14of681989年北京的三所大学组队参加美国的MCM竞赛,此后我国的参赛队伍越来越多.1992－1993年中国工业与应用数学学会（CSIAM）举办了两次中国大学生数学建模竞赛.1994年起，由国家教委（教育部）高教司和中国工业与应用数学学会共同于每年9月举办，1999年开始设立大专组的竞赛.第一讲数学建模简介及数学规划模型Page15of68无论是美国还是我国大学本科组数学建模竞赛题每年都是两道，参赛队从中任选一道题目.一般来说其中一道是连续型，另一道是离散型；或者一道是开放型的，另一道是严谨型的.竞赛内容或题目是由工程技术、管理科学中的实际问题简化而成，留有充分余地供参赛者发挥其聪明才智和创造精神.竞赛形式为三名学生组成一队，可以自由地收集资料、调查研究，使用计算机、因特网和任何软件，在三天时间内分工合作完成一篇论文.评奖标准为模型假设的合理性、建模的创造性、结果的准确性和文字表述的清晰程度.第一讲数学建模简介及数学规划模型Page16of68初等模型第一讲数学建模简介及数学规划模型Page17of68一辆汽车在拐弯时急刹车，结果冲到路边的沟里（见下图），交通警察立即赶到了事故现场。司机申辩说，当他进入弯道时刹车失灵，他还一口咬定，进入弯道其车速为每小时４０英里（这是该路的速度上限，约合每秒１７.９２米）。警察验车时证实该车的制动器在事故发生时确实失灵，然而，司机所说的车速是否真实可信呢？第一讲数学建模简介及数学规划模型Page18of68汽车的最终位置刹车痕迹现在，让我们帮警察计算一下司机所报速度的真实性。连接刹车痕迹的初始点和终点，用x表示沿连线汽车横向所走出的距离，用y表示竖直的距离，如下图X0369121516.6182124273033.3Y01.192.152.823.283.533.553.543.312.892.221.290第一讲数学建模简介及数学规划模型Page19of68上面的表中，我们给出了外侧刹车痕迹的有关值，而且，经过测量还发现，该车并没有偏离它所行驶的转弯路线，也就是说，它的车头一直指向切线方向。可以假设，该车的重心是沿一个半径为r的圆做圆周运动。假设磨擦力作用在该车速度的法线方向上，并设汽车的速度v是一个常数。显然，磨擦力提供了向心力，设磨擦系数为μ,则rvmmg2其中m为汽车质量.由上式易得grv如何计算圆周半径r？假设已知弦的长度为c，弓形的高度为h，其图如下所示，由勾股定理知第一讲数学建模简介及数学规划模型Page20of682222)c()hr(r由前面的表中代入近似数据c=33.27,h=3.55后，得r=40.75米根据实际路面与汽车轮胎的情况，可以测量出磨擦系数，经过实际测试得到g=8.175米／秒２将此结果代入我们上面利用第二定律所得到的式子中，得v≈18.25米／秒此结果比司机所报速度（17.92米／秒）略大。但是，我们不得不考虑计算半径r及测试时的误差。如果误差允许在１０％以内，无疑，此计算结果对司机是相当有利的。第一讲数学建模简介及数学规划模型Page21of68椅子能在不平的地面上放稳吗？把四只脚的椅子往不平的地面上一放，通常只有三只脚着地，放不稳，然而有人认为只要稍挪动几次，就可以四脚着地，放稳了，对吗？第一讲数学建模简介及数学规划模型Page22of68问题分析通常三只脚着地放稳的标准:四只脚着地四条腿一样长，椅脚与地面点接触，四脚连线呈正方形;•地面高度连续变化，可视为数学上的连•续曲面;•地面相对平坦，使椅子在任意位置至少三只脚•同时着地。模型假设第一讲数学建模简介及数学规划模型Page23of68建立模型用数学语言把椅子位置和四只脚着地的关系表示出来.•椅子位置利用正方形(椅脚连线)的对称性用(对角线与x轴的夹角)表示椅子位置•四只脚着地椅脚与地面距离为零距离是的函数xBADCOD´C´B´A´四个距离(四只脚)两个距离正方形对称性正方形ABCD绕O点旋转A,C两脚与地面距离之和记为f()B,D两脚与地面距离之和记为g()第一讲数学建模简介及数学规划模型Page24of68用数学语言把椅子位置和四只脚着地的关系表示出来.f(),g()是连续函数对任意,f(),g()至少一个为0数学问题已知：f(),g()是连续函数;对任意，f()•g()=0;且g(0)=0，f(0)0.证明：存在0，使f(0)=g(0)=0.地面为连续曲面椅子在任意位置至少三只脚着地第一讲数学建模简介及数学规划模型Page25of68模型求解将椅子旋转900，对角线AC和BD互换.由g(0)=0，f(0)0，知f(/2)=0,g(/2)0.令h()=f()–g(),则h(0)0和h(/2)0.由f,g的连续性知h为连续函数,据连续函数的基本性质,必存在0,使h(0)=0,即f(0)=g(0).因为f()•g()=0,所以f(0)=g(0)=0.评注和思考建模的关键:和f(),g()的确定.模型假设中四脚呈正方形不是本质的,读者可考虑长方形的情形.第一讲数学建模简介及数学规划模型Page26of68数学规划模型第一讲数学建模简介及数学规划模型Page27of68实际问题中的优化模型mixgtsxxxxfzMaxMiniTn,2,1,0)(..),(),()(1或x~决策变量f(x)~目标函数gi(x)0~约束条件多元函数条件极值决策变量个数n和约束条件个数m较大最优解在可行域的边界上取得重点在模型的建立和结果的分析第一讲数学建模简介及数学规划模型Page28of68•无约束优化•线性规划•非线性规划•整数规划•多目标规划•动态规划等等第一讲数学建模简介及数学规划模型Page29of68线性规划第一讲数学建模简介及数学规划模型Page30of68设每月生产小、中、大型汽车的数量分别为x1,x2,x3321432xxxzMax600535.1..321xxxts60000400250280321xxx0,,321xxx汽车厂生产计划模型建立小型中型大型现有量钢材1.535600时间28025040060000利润234线性规划模型(LP)第一讲数学建模简介及数学规划模型Page31of68模型求解3）模型中增加条件：x1,x2,x3均为整数，重新求解。OBJECTIVEFUNCTIONVALUE1)632.2581VARIABLEVALUEREDUCEDCOSTX164.5161290.000000X2167.7419280.000000X30.0000000.946237ROWSLACKORSURPLUSDUALPRICES2)0.0000000.7311833)0.0000000.0032261）舍去小数：取x1=64，x2=167，算出目标函数值z=629，与LP最优值632.2581相差不大。2）试探：如取x1=65，x2=167；x1=64，x2=168等，计算函数值z，通过比较可能得到更优的解。•但必须检验它们是否满足约束条件。为什么？结果为小数，怎么办？IP可用LINDO直接求解整数规划(IntegerProgramming,简记IP)“gin3”表示“前3个变量为整数”，等价于：ginx1ginx2ginx3IP的最优解x1=64，x2=168，x3=0，最优值z=632max2x1+3x2+4x3st1.5x1+3x2+5x3600280x1+250x2+400x360000endgin3OBJECTIVEFUNCTIONVALUE1)632.0000VARIABLEVALUEREDUCEDCOSTX164.000000-2.000000X2168.000000-3.000000X30.00000