微 分 方 程 建 模

,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,模型,1,当我们描述实际对象的某些特性随时间（空间）而演变的过程、分析它的变化规律、预测它的未来形态、研究它的控制手段时,，,通常要建立对象的动态模型。,在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立,微分方程模型,的方法来研究该问题。,2,我们来建立如下的一些问题的模型：,我们通过一些最简单的实例来说明微分方程建模的一般方法。在连续变量问题的研究中，微分方程是十分常用的数学工具之一。,1,、,Malthus,模型,2,、,Logistic,模型,3,为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。首先我们建立两个简单的单种群增长模型，以简略分析一下这方面的问题。,种群的数量本应取离散值，但由于种群数量一般较大，为建立微分方程模型，可将种群数量看作连续变量,由此引起的误差将是十分微小的。,4,世界人口数量统计数据：,年,1625,1830,1930,1960,1974,1987,1999,人口亿,5,10,20,30,40,50,60,中国人口数量统计数据：,年,1908,1933,1953,1964,1982,1990,2000,人口,3.0,4.7,6.0,7.2,10.3,11.3,12.95,5,1,马尔萨斯（,Malthus,）模型,马尔萨斯在分析人口出生与死亡情况的资料后发现，人口净增长率,r,基本上是一常数，（,r,=,b,-,d,b,为出生率，,d,为死亡率），因而提出了著名的,人口指数增长模型,。,分析与建模：,人口的净增长率是一个常数，也就是单位时间内人口增长量与当时人口数成正比。,设,t,时刻人口数为,N(t),，,t=t,0,时，,N(t,0,)=N,0,，,则,6,这个方程的解为：,马尔萨斯模型的一个显著特点,：,种群数量翻,一番所需的时间是固定的,。,令种群数量翻一番所需的时间为,T,，则有：,故,即,Malthus,模型,7,模型检验,比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如，,1961,年世界人口数为,30.6,（即,3.0610,9,），人口增长率约为,2%,，人口数大约每,35,年增加一倍。检查,1700,年至,1961,的,260,年人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每,34.6,年增加一倍，两者也几乎相同。,模型预测,假如人口数真能保持每,34.6,年增加一倍，那么人口数将以几何级数的方式增长。例如，到,2510,年，人口达,210,14,个，即使海洋全部变成陆地，每人也只有,9.3,平方英尺的活动范围，而到,2670,年，人口达,3610,15,个，只好一个人站在另一人的肩上排成二层了。 故,马尔萨斯模型是不完善的。,几何级数的增长,Malthus,模型实际上只有在群体总数不太大时才合理，到总数增大时，生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。,所以,Malthus,模型假设的人口,净增长率不可能始终保持常数，它应当与人口数量有关。,8,2,Logistic,模型,人口净增长率应当与人口数量有关，即：,r,=,r,(,N,),从而有：,（,1,）,r,(,N,),是未知函数，但根据实际背景，它无法用拟合方法来求 。,为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。,r,(,N,),最简单的形式是常数，此时得到的就是马尔萨斯模型。对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）,此时得到微分方程：,或,（,2,）,（,2,）,被称为,Logistic,模型或生物总数增长的统计筹算律，是由荷兰数学生物学家弗赫斯特（,Verhulst,）首先提出的。一次项系数是负的，因为当种群数量很大时，会对自身增大产生抑制性，故一次项又被称为竞争项。,（,2,）,可改写成：,（,3,）,(3),式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为,K,（近似地将,K,看成常数），,N,表示当前的种群数量，,K,-,N,恰为环境还能供养的种群数量，（,3,）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（,3,）也被称为统计筹算律的原因。,9,对,（,3,）,分离变量：,两边积分并整理得：,令,N,(0)=,N,0,，求得：,故,（,3,）,的满足初始条件,N,(0)=,N,0,的解为：,（,4,）,易见：,N,(0)=,N,0,，,N,(,t,),的图形请看右图,10,模型检验,用,Logistic,模型来描述种群增长的规律效果如何呢？,1945,年克朗比克（,Crombic,）做了一个人工饲养小谷虫的实验，数学生物学家高斯（,EFGauss,）也做了一个原生物草履虫实验，实验结果都和,Logistic,曲线十分吻合。,大量实验资料表明用,Logistic,模型来描述种群的增长，效果还是相当不错的。例如，高斯,把,5,只草履虫放进一个盛有,0.5cm,3,营养液的小试管，他发现，开始时草履虫以每天,230.9%,的速率增长，此后增长速度不断减慢，到第五天达到最大量,375,个，实验数据与,r,=2.309,，,a,=0.006157,，,N,(0)=5,的,Logistic,曲线：,几乎完全吻合，见右图,11,Malthus,模型和,Logistic,模型的总结,Malthus,模型和,Logistic,模型,均为对微分方程（,1,）所作的模拟近似方程。前一模型假设了种群增长率,r,为一常数，（,r,被称为该种群的内禀增长率）。后一模型则假设环境只能供养一定数量的种群，从而引入了一个竞争项。,用模拟近似法建立微分方程来研究实际问题时必须对求得的解进行检验，看其是否与实际情况相符或基本相符。相符性越好则模拟得越好，否则就得找出不相符的主要原因，对模型进行修改。,Malthus,模型与,Logistic,模型虽然都是为了研究种群数量的增长情况而建立的，但它们也可用来研究其他实际问题，只要这些实际问题的数学模型有相同的微分方程即可。,12,