拆分方程建模分析.ppt

第118页，共173页，星期日，2025年，2月5日模型检验比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如，1961年世界人口数为30.6（即3.06×109），人口增长率约为2%，人口数大约每35年增加一倍。检查1700年至1961的260年人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每34.6年增加一倍，两者也几乎相同。模型预测几何级数的增长Malthus模型实际上只有在群体总数不太大时才合理，到总数增大时，生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。所以Malthus模型假设的人口净增长率不可能始终保持常数，它应当与人口数量有关。第119页，共173页，星期日，2025年，2月5日模型2Logistic模型人口净增长率应当与人口数量有关，即：r=r(N)从而有：（3.7）r(N)是未知函数，但根据实际背景，它无法用拟合方法来求。为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。r(N)最简单的形式是常数，此时得到的就是马尔萨斯模型。对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）此时得到微分方程：或（3.8）（3.8）被称为Logistic模型或生物总数增长的统计筹算律，是由荷兰数学生物学家弗赫斯特（Verhulst）首先提出的。一次项系数是负的，因为当种群数量很大时，会对自身增大产生抑制性，故一次项又被称为竞争项。（3.8）可改写成：（3.9）(3.9)式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为K（近似地将K看成常数），N表示当前的种群数量，K-N恰为环境还能供养的种群数量，（3.9）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（3.9）也被称为统计筹算律的原因。第120页，共173页，星期日，2025年，2月5日图3-5对（3.9）分离变量：两边积分并整理得：令N(0)=N0，求得：故（3.9）的满足初始条件N(0)=N0的解为：（3.10）易见：N(0)=N0，N(t)的图形请看图3.5第121页，共173页，星期日，2025年，2月5日模型检验用Logistic模型来描述种群增长的规律效果如何呢？1945年克朗皮克（Crombic）做了一个人工饲养小谷虫的实验，数学生物学家高斯（E·F·Gauss）也做了一个原生物草履虫实验，实验结果都和Logistic曲线十分吻合。大量实验资料表明用Logistic模型来描述种群的增长，效果还是相当不错的。例如，高斯把5只草履虫放进一个盛有0.5cm3营养液的小试管，他发现，开始时草履虫以每天230.9%的速率增长，此后增长速度不断减慢，到第五天达到最大量375个，实验数据与r=2.309，a=0.006157，N(0)=5的Logistic曲线：几乎完全吻合，见图3.6。图3-6第122页，共173页，星期日，2025年，2月5日第86页，共173页，星期日，2025年，2月5日第87页，共173页，星期日，2025年，2月5日第88页，共173页，星期日，2025年，2月5日第89页，共173页，星期日，2025年，2月5日第90页，共173页，星期日，2025年，2月5日第91页，共173页，星期日，2025年，2月5日第92页，共173页，星期日，2025年，2月5日第93页，共173页，星期日，2025年，2月5日第94页，共173页，星期日，2025年，2月5日第95页，共173页，星期日，2025年，2月5日●Malthus模型设xn是某人类群体在第n个时间段(例如年)末时的总数，若在单位时间段内人口相对增长率为r（出生率与死亡率之差）,那么人口增长数与原人口数成正比,从而xn+1＝xn＋rxn即xn+1=axn其中a=r+1.差分形式的人口增长模型第96页，共173页，星期日，2025年，2月5日这是一个如下线性映射的迭代f(x)=ax从而