六讲微分方程模型.ppt

,3.1,小船行走路线问题,3.2,单种群模型与人口问题,3.3,交通管理中的黄灯问题,附：,数学建模常用软件简介,0,图,3.1.1,3.1,小船行走路线问题,一小船从河边某点驶向对岸码头,若考虑水的流速影响,船行走的路线如何,?,如图,3.1.1,建立坐标系,问 题,(1),水流方向为,y,轴正向，速度大小为,a,;,模型假设,(2),船在,A,处,轮船匀速行驶,速度为,b,为了到达码头,总是朝向码头,O,前进；,(3),船行走的路线为,(4),河宽为,l,米,.,在曲线,y=y,(,x,),上任取一点,P,(,x,y,),因为水流方向为,y,轴正向，大小为,a,所以水流速度向量,0,图,3.1.1,划船方向指向原点,O,(0,0),，,大小为,b,.,模型建立,划船速度向量；,水流速度向量；,船行速度向量,所以划船速度向量,因此船行速度向量,0,图,3.1.1,(1),由船行速度向量的方向为船行路线的切线方向，所以有,或,0,图,3.1.1,于是我们的问题是求上述方程满足条件,的解,.,（,2,）,由 ，可求出任意常数,.,方程,(1),是齐次方程，由齐次方程解法得,模型求解,（*）,此外,有,（*）,(1),当,x,=0,时确有,y,=0,，,故小船一定能到达码头,.,思考：这条路线是最优路线吗？若不是，考虑最优路线,应以时间最短为标准,.,最佳路线应是,O,到,A,的直线,.,由,(*),、,(*),可解出,动植物种群本身是离散变量，谈不到可微性，但由于突然增加或减少的只是单一个体或少数几个个体，与全体数量相比，这种增量是很微小的，所以我们可以近似地假设大规模种群随时间是连续地甚至是可微地在变化，进而可以引用微分方程这一数学工具来研究,.,3.2,单种群模型与人口问题,模型,1 Malthus,模型,这个模型的,基本假设,：在人口自然增长过程中，人口的净增长率为常数,即单位时间内人口增量与当时的人口总量成正比,.,英国人口统计学家马尔萨斯,(Malthus,1766,1834),于,1798,年提出了著名的人口指数增长模型,根据,Malthus,的假设,在,t,到,t,+,t,时间内人口的增长量为,设时刻,t,的人口数为,p,(,t,),t,=,t,0,时的人口数为,p,0,人口增长率为,r,如何建立,Malthus,的数学模型呢,?,令,t0,p,(,t,),满足方程,（,3,）,这个问题的解为,(3),是一个线性微分方程，称为,Malthus,模型,.,模型中的参数,p,0,r,可用所给数据用最小二乘法拟合得到（第六章讨论）,如果,r,0,，上式表明人口将以指数规律无限增长特别地，当,t,+,时,p,(,t,)+,这似乎不太可能,另外，此模型虽然与十九世纪以前欧洲一些地区的人口统计数据可以很好地吻合但是当人们用十九世纪以后许多国家的人口统计资料与,Malthus,模型比较时，却发现了很大的差异,附录,列出了美国十九世纪、二十世纪的人口统计数据与这个模型比较的结果,用这个模型预报的结果远远超过了实际人口的增长引起误差的原因是,10,年增长率,估计过高,按照附录中第,2,列给出的实际人口可以算出，,19,世纪,100,年和,20,世纪前,80,年的,10,年增长率分别为,0.266,和,0.137,，远小于,1790,到,1800,年的增长率,0.307,这个事实对,r,是常数的基本假设提出了异议,.,产生上述现象的主要原因是，随着人口的增加，自然资源、环境条件等因素对人口继续增加的抑制作用越来越显著，如果人口较少时,(,相对于资源而言,),人口增长率还可以看作常数的话，那么当人口增加到一定数量后,增长率就会随着人口的继续增加而逐渐减少,.,许多国家人口增长的实际情况完全证实了这一点,.,我们不妨利用附录第,2,列给出的数据计算一下美国人口每,10,年的增长率,可以知道大致是逐渐下降的,.,将增长率,r,表示为人口,p,(,t,),的,函数,r,(,p,),按照前面的分析,r,(,p,),是,p,的减函数,.,一个最简单的假定是：设,r,(,p,),是,p,的线性函数，即,此模型是荷兰生物学家,Verhulst,于十九世纪中叶提出,模型,Logistic,模型,修改假设：,为了使人口预报特别是长期预报更好地符合实际情况，必须修改,Malthus,模型关于人口增长率是常数这个基本假设,.,从而得到,Logistic,模型,.,这里,r,相当于,p,=0,时的增长率，称,固有增长率,为确定系数,s,的意义，引入自然资源和环境条件所能容纳的最,大人口容量,p,m,，称,最大人口容量,当,p,=,p,m,时增长率为零，即,r,(,p,m,)=0,，由此可确定出,s,=,r,/,p,m,.,显然对于任意的,p,0,，增长率,r,(,p,)dsolve(D2y+4*Dy+12*y=0,y(0)=0,Dy(0)=5,x),ans=,5/4*2(1/2)*exp(-2*x)*sin(2*2(1/2)*x),即,THANK YOU,SUCCESS,2025/1/15 周三,28,可编辑,