6.1人口模型ppt课件

作业,搜集某个地区(国、省、市)的人口数据（不少于6个数据），1）保留时间上最晚的两个数据做检验，根据其余数据分别用Malthus模型和Logistic模型拟合最佳参数，并检验模型在上述两个用于检验的数据上的误差。2）根据样本数据散点图特点，考虑采用Malthus和Logistic两个经典模型对搜集的数据建模是否合适？如不合适，请给出具体改进的建议方法。,1,常微分方程模型,2,3,6.1人口增长模型,英国人口学家Malthus(1766-1834),模型假设,人口自然增长率r为常数,模型建立,1.指数增长模型,4,人口学原理的基本思想：如没有限制，人口是几何级数（即：2，4，8，16，32，64，128等）增长。而食物供应呈现线性（即：1，2，3，4，5，6，7等）增长。食物为人类生存的最重要之条件。只有自然原因（事故和衰老），灾难（战争，瘟疫，及各类饑荒），道德限制和罪恶（马尔萨斯所指包括杀婴，谋杀，节育和同性恋）能够限制人口的过度增长马尔萨斯注意到许多人误用他的理论，痛苦地阐明他没有仅仅预测未来的大灾难。他辩解道，“周期性灾难持续存在的原因自人类有史以来就已经存在，目前仍然存在，并且将来会继续存在，除非我们的大自然的物理结构发生决定性的变化。”因此，马尔萨斯认为他的人口论是对人类过去和目前状况的解释，以及对我们未来的预测。人口论的论点为粮食增加仅会呈等差数列，而人口的增加却会呈现等比数列。所以他觉得人类必须顾虑食物的缺乏，而减少结婚的预防限制，以及受现实穷困的折磨、对已生人口所加压迫积极限制。,AnEssayonthePrincipleofPopulation,1798,5,模型分析,人口将按指数规律无限增长！,人口将始终保持不变！,人口将按指数规律减少直至绝灭！,人口倍增时间,模型求解,人口以几何级数增加！,6,Malthus模型预测美国人口,7,Malthus模型预测的优缺点,优点,短期预报比较准确,缺点,不适合中长期预报,原因,预报时假设人口增长率r为常数。没有考虑环境对人口增长的制约作用。,8,2.阻滞增长模型,假设人口增长率r是人口x的减函数：,其中，xm为考虑到受自然资源和环境条件限制所能容纳的最大人口数量,（称环境容纳量）,模型假设,模型建立,9,模型分析（定性分析）,人口将递减并趋向于xm！,人口将始终保持xm不变！,人口将递增并趋向于xm！,无论在哪种情况下，人口最终将趋向于最大人口容量！,模型求解,10,人口增长率达到最大值,11,阻滞增长模型预测美国人口,12,阻滞增长模型预测的优缺点,优点,中期预报比较准确,缺点,理论上很好，实用性不强,原因,预报时假设固有人口增长率r以及最大人口容量xm为定值。,13,functionmain()%PopulationofUSclearclcdigits(2);unit=1e6;year=0:10:210;years=year+1790;RealPopulation=3.93,5.31,7.24,9.64,12.87,17.07,23.19,31.44,38.56,.50.16,62.95,76.0,91.97,105.71,122.76,131.67,150.70,179.32,203.21,.226.51,248.71,281.42;%非拟合方法估计模型参数x0=3.93;x1=5.31;r=log(x1/x0)/10.0;ExpPopu=x0.*exp(r.*(year-year(1);errorEXP=100.*(ExpPopu-RealPopulation)./RealPopulation;error1=norm(ExpPopu-RealPopulation)*unitsubplot(1,2,1)plot(years,RealPopulation*unit,b,years,ExpPopu*unit,r);title(Malthus模型预测美国人口);subplot(1,2,2)plot(years,100,r*,years,errorEXP,b);title(Malthus模型预测美国人口相对误差);,14,rl=0.31;xm=197;LogPopu=xm./(1+(xm/x0-1).*exp(-r.*(year-year(1);errorLog=100.*(LogPopu-RealPopulation)./RealPopulation;error2=norm(LogPopu-RealPopulation)*unitsubplot(1,2,1)plot(years,RealPopulation*unit,b,years,LogPopu*unit,r);title(阻滞增长模型预测美国人口);subplot(1,2,2)plot(years,errorLog,b);title(阻滞增长模型预测美国人口相对误差);pause,15,%EstimatethreeparametersinLogisticmodel:p(1)=growthrater,p(2)=carryingcapacityK,%andp(3)=initialpopulationP0.p0=.01,1000,x0;%initialestimatedparametervalueslowerbound=zeros(1,3);options=optimset(MaxFunEvals,10000);p,Error=lsqcurvefit(logistic,p0,year,RealPopulation,lowerbound,options)error3=sqrt(Error)*unitmodelpops=logistic(p,year);errorEstLog=100.*(modelpops-RealPopulation)./RealPopulation;subplot(1,2,1)plot(years,RealPopulation*unit,o,years,modelpops*unit)title(参数估计阻滞增长模型预测美国人口);subplot(1,2,2)plot(years,errorEstLog,b);title(参数估计阻滞增长模型预测美国人口相对误差);functionP=logistic(p,t);%MATLABfunctionfilethattakestimet,growthrat