数学建模讲座玩具、照片…~实物模型风洞中的飞机…~物理模型地图、电路图…~符号模型模型是为了一定目的,对客观事物的一部分进行简缩、抽象、提炼出来的原型的替代物。模型集中反映了原型中人们需要的那一部分特征。我们常见的模型什么是数学模型第一章建立数学模型你碰到过的数学模型——“航行问题”甲乙两地相距750公里,船从甲到乙顺水航行需30小时,从乙到甲逆水航行需50小时,问船的速度是多少。用x表示船速,y表示水速,列出方程:75050)(75030)(yxyx求解得到x=20,y=5,答:船速每小时20公里航行问题建立数学模型的基本步骤•作出简化假设(船速、水速为常数);•用符号表示有关量(x,y表示船速和水速);•用物理定律(匀速运动的距离等于速度乘以时间)列出数学式子(二元一次方程);•求解得到数学解答(x=20,y=5);•回答原问题(船速每小时20公里)。数学模型(MathematicalModel)和数学建模(MathematicalModeling)数学模型:对于一个现实对象,为了一个特定目的,根据其内在规律,作出必要的简化假设,运用适当的数学工具,得到的一个数学结构。数学建模:建立数学模型的全过程(包括建立、求解、分析、检验)。数学建模的重要意义•电子计算机的出现及飞速发展•数学以空前的广度和深度向一切领域渗透数学建模作为用数学方法解决实际问题的第一步,越来越受到人们的重视。数学建模计算机技术如虎添翼知识经济建模示例椅子能在不平的地面上放稳吗?问题椅子能在不平的地面上放稳吗?1.椅子四条腿一样长,椅脚与地面接触处可视为一人点,四脚的连线呈正方形;2.地面高度是连续变化的,沿任何方向都不会出现间断(没有像台阶那样的情况),即地面可视为数学上的连续曲面;3.对于椅脚的间距和椅腿的长度而言,地面是相对平坦的,使椅子的任何位置至少有三只脚同时着地。模型假设ABCDtA‘B‘C‘D‘Ox模型构成椅脚连线为正方形ABCD(如右图)。t~椅子绕中心点O旋转角度f(t)~A,C两脚与地面距离之和g(t)~A,C两脚与地面距离之和f(t),g(t)0模型构成由假设1,f和g都是连续函数由假设3,椅子在任何位置至少有三只脚同时着地:对任意t,f(t)和g(t)中至少有一个为0。当t=0时,不妨设g(t)=0,f(t)0,原题归结为证明如下的数学命题:已知f(t)和g(t)是t的连续函数,对任意t,f(t)•g(t)=0,且g(0)=0,f(0)0。则存在t0,使f(t0)=g(t0)=0模型求解OxA‘B‘C‘D‘ABCDt最后,因为f(t)•g(t)=0,所以f(t0)=g(t0)=0。令h(t)=f(t)-g(t),则h(0)0和h()0,由f和g的连续性知h也是连续函数。根据连续函数的基本性质,必存在t0(0t0),使h(t0)=0,即f(t0)=g(t0)。22将椅子旋转90º,对角线AC与BD互换。由g(0)=0,f(0)0可知g()0,f()=022建模示例商人们怎样安全过河问题(智力游戏)3名商人3名随从河小船(至多2人)随从们密约,在河的任一岸,一旦随从的人数比商人多,就杀人越货.但是乘船渡河的方案由商人决定.商人们怎样才能安全过河?问题分析多步决策过程决策~每一步(此岸到彼岸或彼岸到此岸)船上的人员要求~在安全的前提下(两岸的随从数不比商人多),经有限步使全体人员过河模型构成xk~第k次渡河前此岸的商人数yk~第k次渡河前此岸的随从数xk,yk=0,1,2,3;k=1,2,sk=(xk,yk)~过程的状态S={(x,y)x=0,y=0,1,2,3;x=3,y=0,1,2,3;x=y=1,2}S~允许状态集合uk~第k次渡船上的商人数vk~第k次渡船上的随从数dk=(uk,vk)~决策D={(u,v)u+v=1,2}~允许决策集合uk,vk=0,1,2;k=1,2,sk+1=skdk+(-1)k~状态转移律求dkD(k=1,2,n),使skS按转移律由s1=(3,3)到达sn+1=(0,0).多步决策问题模型求解xy3322110穷举法~编程上机图解法状态s=(x,y)~16个格点~10个点允许决策D~移动1或2格;k奇,左下移;k偶,右上移.s1sn+1d1,d11给出安全渡河方案评注和思考规格化方法,易于推广考虑4名商人各带一随从的情况d1d11允许状态SS={(x,y)x=0,y=0,1,2,3;x=3,y=0,1,2,3;x=y=1,2}D={(u,v)u+v=1,2}习题•模仿这一案例,作下面一题:人带着猫、鸡、米过河,船除需要人划之外,至多能载猫、鸡、米三者之一,而当人不在场时猫要吃鸡、鸡要吃米。试设计一安全过河方案,并使渡河次数尽量地少。背景年1625183019301960197419871999人口(亿)5102030405060世界人口增长概况中国人口增长概况年1908193319531964198219901995人口(亿)34.76710.111.312研究人口变化规律控制人口过快增长建模示例如何预报人口的增长指数增长模型常用的计算公式kkrxx)1(0马尔萨斯(1788--1834)提出的指数增长模型(1798)x(t)~时刻t人口r~人口(相对)增长率(常数)ttrxtxttx)()()(今年人口x0,年增长率rk年后人口0)0(,xxrxdtdxrtextx0)(trextx)()(0trx)1(0随着时间增加人口按指数规律无限增长指数增长模型的应用及局限性•与19世纪以前欧洲一些地区人口统计数据吻合•适用于19世纪后迁往加拿大的欧洲移民后代•可用于短期人口增长预测•不符合19世纪后多数地区人口增长规律•不能预测较长期的人口增长过程19世纪后人口数据人口增长率r不是常数(逐渐下降)阻滞增长模型(Logistic模型)人口增长到一定数量后,增长率下降的原因:资源、环境等因素对人口增长的阻滞作用且阻滞作用随人口数量增加而变大假定:)0,()(srsxrxrr~固有增长率(x很小时)xm~人口容量(资源、环境能容纳的最大数量))1()(mxxrxrr是x的减函数mxrs0)(mxr阻滞增长模型(Logistic模型)rxdtdx)1()(mxxrxxxrdtdxdx/dtx0xmxm/2xmxtxxxemmrt()()110tx0x(t)~S形曲线,x增加先快后慢x0xm/2模型的参数估计用指数增长模型或阻滞增长模型作人口预报,必须先估计模型参数r或r,xm•利用统计数据用最小二乘法作拟合例:美国人口数据(单位~百万)17901800181018201830……19501960197019803.95.37.29.612.9……150.7179.3204.0226.5r=0.2072,xm=464•专家估计模型检验用模型预报1990年美国人口,与实际数据比较]/)1980(1)[1980()1980()1980()1990(mxxrxxxxx实际为251.4(百万)5.250)1990(x模型应用——人口预报用美国1790~1990年人口数据重新估计参数r=0.2083,xm=457.6x(2000)=275.0x(2010)=297.9Logistic模型在经济领域中的应用(如耐用消费品的售量)基本方法•机理分析•测试分析根据对客观事物特性的认识,找出反映内部机理的数量规律将研究对象看作“黑箱”,通过对量测数据的统计分析,找出与数据拟合最好的模型机理分析没有统一的方法,主要通过实例研究(CaseStudies)来学习。以下建模主要指机理分析•二者结合机理分析建立模型结构,测试分析确定模型参数数学建模的方法和步骤数学建模的一般步骤模型准备模型假设模型构成模型求解模型分析模型检验模型应用怎样学习数学建模数学建模与其说是一门技术,不如说是一门艺术技术大致有章可循艺术无法归纳成普遍适用的准则想象力洞察力判断力•学习、分析、评价、改进别人作过的模型•亲自动手,认真作几个实际题目创新意识看谁答得快1、某甲早8时从山下旅店出发沿一路径上山,下午5时到达山顶并留宿。次日早8时沿同一路径下山,下午5时回到旅店。某乙说,甲必在两天中的同一时刻经过路径中的同一地点,为什么?2、某人家住T市在他乡工作,每天下班后乘火车于6时抵达T市车站,他的妻子驾车准时到车站接他回家。一日他提前下班搭早一班火车于5时半抵T市车站,随即步行回家,他的妻子像往常一样驾车前来,在路上遇到他接回家时,发现比往常提前了10分钟,问他步行了多长时间?3、两兄妹分别在离家2千米和1千米且方向相反的两所学校上学,每天同时放学后分别以4千米/小时和2千米/小时的速度步行回家,一小狗以6千米/小时的速度从哥哥处奔向妹妹,又从妹妹处奔向哥哥,如此往返直至回家中,问小狗奔波了多少路程?录象机计数器的用途问题经试验,一盘录象带从头走到尾,时间用了183分30秒,计数器读数从0000变到6152。在一次使用中录象带已经转过大半,计数器读数为4580,问剩下的一段还能否录下1小时的节目?要求不仅回答问题,而且建立计数器读数与录象带转过时间的关系。思考计数器读数是均匀增长的吗?第二章初等模型问题分析录象机计数器的工作原理0000左轮盘右轮盘磁头主动轮压轮计数器录象带录象带运动方向录象带运动右轮盘半径增大右轮转速不是常数录象带运动速度是常数计数器读数增长变慢观察计数器读数增长越来越慢!模型假设•录象带的运动速度是常数v;•计数器读数n与右轮转数m成正比,记m=kn;•录象带厚度(加两圈间空隙)为常数w;•空右轮盘半径记作r;•时间t=0时读数n=0.建模目的建立时间t与读数n之间的关系(设V,k,w,r为已知参数)模型建立建立t与n的函数关系有多种方法1.右轮盘转第i圈的半径为r+wi,m圈的总长度等于录象带在时间t内移动的长度vt,所以knmvtwirmi1)(2nvrknvwkt222模型建立2.考察右轮盘面积的变化,等于录象带厚度乘以转过的长度,即wvtrwknr])[(22nvrknvwkt2223.考察t到t+dt录象带在右轮盘缠绕的长度,有vdtkdnwknr2)(思考vtwirmi1)(2wvtrwknr])[(22vdtkdnwknr2)(nvrknvwkt2221.3种建模方法得到同一结果但仔细推算会发现稍有差别,请解释。2.模型中有待定参数,,,,kvwr确定参数的一种办法是测量或调查,试设计测量方法。参数估计确定参数的另一种方法——测试分析将模型改记作,2bnant只需估计,,ba理论上,已知t=183.5,n=6152,再有一组(t,n)数据即可;实际上,由于测试有误差,最好用足够多的数据作拟合。现有一批测试数据:t020406080n00001153204528003466t100120140160183.5n40684621513556196152用最小二乘法可得.1044.1,1051.226ba模型检验应该另外测试一批数据检验模型:bnant2)1044.1,1051.2(26ba模型应用1.回答提出的问题:由模型算得n=4580时t=118.5分,剩下的录象带能录183.5-118.5=65分钟的节目。2.揭示了“t与n之间呈二次函数关系”这一普遍规律,当录象带的状态改变时,只需重新估计a,b即可。第2章模糊聚类分析§2.1模糊矩阵定义1设R=(rij)m×n,若0≤rij≤1,则称R为模糊矩阵.当rij只取0或1时,称R为布尔(Boole)矩阵.当模糊方阵R=(rij)n×n的对角线上的元素rii都为1时,称R为模糊自反矩阵.定义2设A=(aij)m×n,B=(bij)m×n都是模糊矩阵,相等:A=Baij=bij;包含:A≤Baij≤bij;并:A∪B=(aij∨bi