数学模型(MathematicalModel)是用数学符号、数学式子、程序、图形等对实际课题本质属性的抽象而又简洁的刻划,它或能解释某些客观现象,或能预测未来的发展规律,或能为控制某一现象的发展提供某种意义下的最优策略或较好策略。数学建模(MathematicalModeling)应用知识从实际课题中抽象、提炼出数学模型的过程。§1.1数学模型与数学建模例(万有引力定律的发现)十五世纪中期,哥白尼提出了震惊世界的日心说。丹麦著名的实验天文学家第谷花了二十多年时间观察纪录下了当时已发现的五大行星的运动情况。第谷的学生和助手开普勒对这些资料进行了九年时间的分析计算后得出著名的Kepler三定律。牛顿根据开普勒三定律和牛顿第二定律,利用微积分方法推导出牛顿第三定律即万有引力定律。1.行星轨道是一个椭圆,太太阳位于此椭圆的一个焦点上。2.行星在单位时间内扫过的面积不变。3.行星运行周期的平方正比于椭圆长半轴的三次方,比例系数不随行星而改变(绝对常数)开普勒三大定律这其中必定是某一力学规律的反映,哼哼,我要找出它。。。。如图,有椭圆方程:cos1eprdrdA221矢径所扫过的面积A的微分为:由开普勒第二定律:wrdtdA221常数立即得出:wrwrrwrdtd222)(0即:02wrwr椭圆面积wTrdtdtdAabT2021由此得出Tabwr22常数简单推导如下:行星r太阳我们还需算出行星的加速度,为此需要建立两种不同的坐标架。第一个是固定的,以太阳为坐标原点,沿长轴方向的单位向量记为i,沿短轴方向的单位向量记为j,于是:jirsinrcosr进而有加速度·cossin)(2(sin)(cos()sinr(dtd)cosr(dtd22222))jijijirawrwrrwr以行星为坐标原点建立活动架标,其两个正交的单位向量分别是jiejiecossinsincosθr,因此得出rrwrea)(2由于02wrwr也就是说行星的加速度为rrTaea223214由开普勒第三定律知23/Ta为常数。若记2324MTaG那么就导出著名的万有引力定律:再将椭圆方程)cos1(erp两边微分两次,得0)(1)(2232wrrrprwr将前面得到的结果和焦参数代入,即得22322rTarwr14Tabwr22abp2rrMmGFe21.了解问题的实际背景,明确建模目的,收集掌握必要的数据资料。2.在明确建模目的,掌握必要资料的基础上,通过对资料的分析计算,找出起主要作用的因素,经必要的精炼、简化,提出若干符合客观实际的假设。3.在所作假设的基础上,利用适当的数学工具去刻划各变量之间的关系,建立相应的数学结构——即建立数学模型。4.模型求解。5.模型的分析与检验。在难以得出解析解时,也应当借助计算机求出数值解。§1.2数学建模的一般步骤实体信息(数据)假设建模求解验证应用§1.3数学模型的分类分类标准具体类别对某个实际问题了解的深入程度白箱模型、灰箱模型、黑箱模型模型中变量的特征连续型模型、离散型模型或确定性模型、随机型模型等建模中所用的数学方法初等模型、微分方程模型、差分方程模型、优化模型等研究课题的实际范畴人口模型、生态系统模型、交通流模型、经济模型、基因模型等①数学建模实践的每一步中都蕴含着能力上的锻炼,在调查研究阶段,需要用到观察能力、分析能力和数据处理能力等。在提出假设时,又需要用到想象力和归纳简化能力。②在真正开始自己的研究之前,还应当尽可能先了解一下前人或别人的工作,使自己的工作成为别人研究工作的继续而不是别人工作的重复,你可以把某些已知的研究结果用作你的假设,去探索新的奥秘。因此我们还应当学会在尽可能短的时间内查到并学会我想应用的知识的本领。③还需要你多少要有点创新的能力。这种能力不是生来就有的,建模实践就为你提供了一个培养创新能力的机会。§1.4数学建模与能力的培养开设数学建模课的主要目的为了提高学生的综合素质,增强应用数学知识解决实际问题的本领。仅最近几年里,我校学生都在只参加了半年左右的学习和实践后,就在国际性的竞赛(美国大学生数学建模竞赛)中交出了非常出色的研究论文,夺得了特等奖兼INFORMS奖2项(1999年、2003年各一项)、18项一等奖、10项二等奖的好成绩。•例1某人平时下班总是按预定时间到达某处,然然后他妻子开车接他回家。有一天,他比平时提早了三十分钟到达该处,于是此人就沿着妻子来接他的方向步行回去并在途中遇到了妻子,这一天,他比平时提前了十分钟到家,问此人共步行了多长时间?§1.5一些简单实例似乎条件不够哦。。换一种想法,问题就迎刃而解了。假如他的妻子遇到他后仍载着他开往会合地点,那么这一天他就不会提前回家了。提前的十分钟时间从何而来?显然是由于节省了从相遇点到会合点,又从会合点返回相遇点这一段路的缘故,故由相遇点到会合点需开5分钟。而此人提前了三十分钟到达会合点,故相遇时他已步行了二十五分钟。请思考一下,本题解答中隐含了哪些假设?例2某人第一天由A地去B地,第二天由B地沿原路返回A地。问:在什么条件下,可以保证途中至少存在一地,此人在两天中的同一时间到达该地。分析本题多少有点象数学中解的存在性条件及证明,当然,这里的情况要简单得多。假如我们换一种想法,把第二天的返回改变成另一人在同一天由B去A,问题就化为在什么条件下,两人至少在途中相遇一次,这样结论就很容易得出了:只要任何一人的到达时间晚于另一人的出发时间,两人必会在途中相遇。(请自己据此给出严格证明)•例3交通灯在绿灯转换成红灯时,有一个过渡状态——亮一段时间的黄灯。请分析黄灯应当亮多久。设想一下黄灯的作用是什么,不难看出,黄灯起的是警告的作用,意思是马上要转红灯了,假如你能停住,请立即停车。停车是需要时间的,在这段时间内,车辆仍将向前行驶一段距离L。这就是说,在离街口距离为L处存在着一条停车线(尽管它没被画在地上),见图1-4。对于那些黄灯亮时已过线的车辆,则应当保证它们仍能穿过马路。马路的宽度D是容易测得的,问题的关键在于L的确定。为确定L,还应当将L划分为两段:L1和L2,其中L1是司机在发现黄灯亮及判断应当刹车的反应时间内驶过的路程,L2为刹车制动后车辆驶过的路程。L1较容易计算,交通部门对司机的平均反应时间t1早有测算,反应时间过长将考不出驾照),而此街道的行驶速度v也是交管部门早已定好的,目的是使交通流量最大,可另建模型研究,从而L1=v*t1。刹车距离L2既可用曲线拟合方法得出,也可利用牛顿第二定律计算出来(留作习题)。黄灯究竟应当亮多久现在已经变得清楚多了。第一步,先计算出L应多大才能使看见黄灯的司机停得住车。第二步,黄灯亮的时间应当让已过线的车顺利穿过马路,即T至少应当达到(L+D)/v。DL例4餐馆每天都要洗大量的盘子,为了方便,某餐馆是这样清洗盘子的:先用冷水粗粗洗一下,再放进热水池洗涤,水温不能太高,否则会烫手,但也不能太低,否则不干净。由于想节省开支,餐馆老板想了解一池热水到底可以洗多少盘子,请你帮他建模分析一下这一问题。盘子有大小吗?是什么样的盘子?盘子是怎样洗的?………不妨假设我们了解到:盘子大小相同,均为瓷质菜盘,洗涤时先将一叠盘子浸泡在热水中,然后一清洗。不难看出,是水的温度在决定洗盘子的数量。盘子是先用冷水洗过的,其后可能还会再用清水冲洗,更换热水并非因为水太脏了,而是因为水不够热了。那么热水为什么会变冷呢?假如你想建一个较精细的模型,你当然应当把水池、空气等吸热的因素都考虑进去,但餐馆老板的原意只是想了解一下一池热水平均大约可以洗多少盘子,杀鸡焉用牛刀?不妨可以提出以下简化假设:(1)水池、空气吸热不计,只考虑盘子吸热,盘子的大小、材料相同(2)盘子初始温度与气温相同,洗完后的温度与水温相同(3)水池中的水量为常数,开始温度为T1,最终换水时的温度为T2(4)每个盘子的洗涤时间△T是一个常数。(这一假设甚至可以去掉不要)根据上述简化假设,利用热量守衡定律,餐馆老板的问题就很容易回答了,当然,你还应当调查一下一池水的质量是多少,查一下瓷盘的吸热系数和质量等。可见,假设条件的提出不仅和你研的问题有关,还和你准备利用哪些知识、准备建立什么样的模型以及你准备研究的深入程度有关,即在你提出假设时,你建模的框架已经基本搭好了。例5将形状质量相同的砖块一一向右往外叠放,欲尽可能地延伸到远方,问最远可以延伸多大距离。设砖块是均质的,长度与重量均为1,其重心在中点1/2砖长处,现用归纳法推导。Zn(n-1)n(n+1)由第n块砖受到的两个力的力矩相等,有:1/2-Zn=(n-1)Zn故Zn=1/(2n),从而上面n块砖向右推出的总距离为,nkk121112121,时nnknkn故砖块向右可叠至任意远,这一结果多少有点出人意料。例6某人住在某公交线附近,该公交线路为在A、B两地间运行,每隔10分钟A、B两地各发出一班车,此人常在离家最近的C点等车,他发现了一个令他感到奇怪的现象:在绝大多数情况下,先到站的总是由B去A的车,难道由B去A的车次多些吗?请你帮助他找一下原因AB发出车次显然是一样多的,否则一处的车辆将会越积越多。由于距离不同,设A到C行驶31分钟,B到C要行驶30分钟,考察一个时间长度为10分钟的区间,例如,可以从A方向来的车驶离C站时开始,在其后的9分钟内到达的乘客见到先来的车均为B开往A的,仅有最后1分钟到达的乘客才见到由A来的车先到。由此可见,如果此人到C站等车的时间是随机的,则他先遇上B方向来的车的概率为90%。例4飞机失事时,黑匣子会自动打开,发射出某种射线。为了搞清失事原因,人们必须尽快找回匣子。确定黑匣子的位置,必须确定其所在的方向和距离,试设计一些寻找黑匣子的方法。由于要确定两个参数,至少要用仪器检测两次,除非你事先知道黑匣子发射射线的强度。方法一点光源发出的射线在各点处的照度与其到点光源的距离的平方成反比,即2k/dI黑匣子所在方向很容易确定,关键在于确定距离。设在同一方向不同位置检测了两次,测得的照度分别为I1和I2,两测量点间的距离为a,则有22212)(//dadadkdkII112/IIad方法二在方法一中,两检测点与黑匣子位于一直线上,这一点比较容易做到,主要缺点是结果对照度测量的精度要求较高,很少的误差会造成结果的很大变化,即敏感性很强,现提出另一方法,在A点测得黑匣子方向后,到B点再测方向,AB距离为a,∠BAC=α,∠ABC=β,利用正弦定理得出d=asinα/sin(α+β)。需要指出的是,当黑匣子位于较远处而α又较小时,α+β可能非常接近π(∠ACB接近于0),而sin(α+β)又恰好位于分母上,因而对结果的精确性影响也会很大,为了使结果较好,应使a也相对较大。BACaαβ