人口模型预测——数学建模作业

上传是为了分析数学的乐趣，请粘贴复制的时候也多思考哈。为了更多的学子们。2014年数学建模论文第二套题目人口增长模型的确定专业、姓名土木135提交日期2015/7/2晚上题目人口增长模型的确定摘要对美国人口数据的变化进行拟合，并进行未来人口预测，在第一个模型中，考虑到人口连续变化的规律，用微分方程的方法解出其数量随时间变化的方程，用MATLAB里的CFTOOL工具箱求出参数，即人口净增长率R002222,对该模型与实际数据进行对比，并计算了从1980年后每隔10年的人口数据，与实际对比，有很大出入。因此又改进出更为符合实际的阻滞增长模型，应用微分方程里的分离变量法和积分法解出其数量随时间变化的方程，求出参数人口增长率R002858和人口所能容纳最大值2589,MX与实际数据对比，拟合得很好，并预测出1980年后每隔10年的人口数据，与实际对比，比较符合。为了便于比较两个模型与实际数据的描述情况作对比，又做出了两个模型与实际数据的对比图，以及两个模型的误差图。关键词人口预测微分方程马尔萨斯人口增长模型阻滞增长模型一、问题重述17901980年间美国每隔10年的人口记录如下表所示。表1人口记录表试用以上数据建立马尔萨斯MALTHUS人口指数增长模型，并对接下来的每隔十年预测五次人口数量，并查阅实际数据进行比对分析。如果数据不相符，再对以上模型进行改进，寻找更为合适的模型进行预测。二、问题分析由于题目已经说明首先用马尔萨斯人口增长模型来刻划，列出人口增长指数增长方程并求解，并进行未来50年内人口数据预测，但发现与实际数据有较大出入。考虑到实际的人口增长率是受实际情况制约的，因此，使人口增长率为一变化的线性递减函数，列出人口增长微分方程，求出其方程解，并预测未来五十年内人口实际数据。三、问题假设1假设所给的数据真实可靠2各个年龄段的性别比例大致保持不变3人口变化不受外界大的因素的影响；年份1790180018101820183018401850186018701880人口10639537296129171232314386502年份1890190019101920193019401950196019701980人口10662976092010651232131715071793204022654马尔萨斯人口模型（1）单位时间的人口增长率R为常数；（2）将视为T的连续可微函数。XT5改进后的模型（阻滞增长模型）1人口净增长率R为变化量。四、变量说明T时刻的人口数量XT初始时刻的人口数量1790R人口净增长率环境所能容纳的最大人口数量，即MX0MRX五、模型建立1马尔萨斯人口增长模型T1790时的人口数为,在T到TT这一时间间隔内，人口的增长为1790XTTR由于0LITXXTTA则得到可建立含初始条件的微分方程,39（省略106）TRX1790X其解为1790RTTE2阻滞增长模型假设人口增长生长率为人口的线性递减函数，即。XTM假设自然资源和环境条件所能承受的最大人口容量为，显然，当时，XMX。所以。因此有0MRXR/SR。于是建立下列微分方程，/M1TTR。把上式化为。分离常数并积1793X1790MDXX分得到。17901790MRTXXE六、模型求解1马尔萨斯模型求解参数估计R可以用实际数据的线性最小二乘法求解，对于，1790RTXTE直接求解是比较麻烦的，因此在两边取对数，即，记1790LNLT，A。则原方程化为X39EXPRTLNXTY1790LL36X1790。利用17901900年的数据进行拟合，得到R002142所以也能求出方程程序见附录1。但本题还可以应用MATLAB里的CFTOOL工具箱求参数，在命令行中输入得到更精确的解GENERALMODELFX39EXPRT1790COEFFICIENTSWITH95CONFIDENCEBOUNDSR002222002163,002281得到如图所示结果，其中蓝线表示马尔萨斯人口模型预测人口数据，正方形黑点表示实际人口数据。图1马尔萨斯人口模型与实际人口数据则每隔10年预测人口为,19032X04128X20157X,，然而查阅相关年份美国实际人口数据，1990年为20645X20379X2487百万，2000年为2814百万，2010年为3070百万。对于2020年和2030年实际还没有统计，因为没有发生，但通过前三个数据就可以看出马尔萨斯模型预测人口与实际有很大出入，所以必须对该模型做出改进，得到更符合实际的预测模型。2阻滞增长模型求解通过对求导得拐点在时，人口增长速度最大。在问题分析已经得到XT/2MX该模型的表达式，运用MATLAB里的CFTOOL工具箱拟合求出参数GENERALMODELFXA39/39A39EXPRT1790COEFFICIENTSWITH95CONFIDENCEBOUNDSA28592574,3144R002858002763,002953因此。并得到如下图，蓝线表示组织增028517902913TXE长模型预测数据，黑点表示实际人口数据。图2组织增长模型预测数据与实际人口数据根据该方程预测得到23092，24251，25202，25967，2657190X20X201X20X203X1其中1990，2000,2010年这三年的预测人口数斗鱼实际人口数据很接近。但还是有一定的误差，模型也存在一定的改进程度才能更符合实际情况。但从图形看，与实际拟合的很好。3为了便于比较两种模型与实际数据的直观对比，编出程序附录2把他们放在一个坐标系里。图3两个模型与实际人口数据的对比17801800182018401860188019001920194019601980050100150200250300T/万万X/万万万万万万万万万万万万万万万万万万万万万万万万万万图形虽然直观，但不具体，因此应算出两种模型与实际的误差值比较，程序见附录3得到下图。图4马尔萨斯模型与阻滞增长模型误差的比较17801800182018401860188019001920194019601980020100102030405万万万万万万万万万万万万万万万万万万万万万万从图中可以看出阻滞增长模型的误差更小。七、结果分析1马尔萨斯模型结果分析则每隔10年预测人口为,19032X04128X20157X,，然而查阅相关年份美国实际人口数据，1990年为20645X2372487百万，2000年为2814百万，2010年为3070百万。对于2020年和2030年实际还没有统计，因为没有发生，但通过前三个数据就可以看出马尔萨斯模型预测人口与实际有很大出入，所以必须对该模型做出改进，得到更符合实际的预测模型。2阻滞增长模型结果分析根据该方程预测得到23092，24251，25202，25967，2657190X20X201X20X203X1其中1990年实际人口为2487百万，2000年为2814百万，2010年为3070百万，这三年的预测人口数与实际人口数据很接近。但还是有一定的误差，模型也存在一定的改进程度才能更符合实际情况。但从图形看，与实际拟合的比较好。八，模型的评价与推广MALTHUS数学模型在短期内具有较好的准确度，简易易行，但是不能准确的预测处人口长期的发展趋势，不具有预测人口长期增长数量的能力。为此，结合资料，考虑到一些实际因素，建立了能很好滴预测人口数量增长的LOGSTIC模型。在人口增长的整个过程中LOGISTIC模型预测的数据与题中所给数据能很好地在误差范围内，几乎一致。但由于也存在误差，因此也可以通过相关多项式拟合出其方程，也是可以的，比如二次多项式，但次数不一定越高越好，应使模型所预的数据与实际数据更接近，才是比较好的模型。LOGISTIC模型在人口预测中，在医疗卫生中可以预测寻找某一疾病的危险因素（以及疾病的发展趋势），预测自然界中种群数量的增长等都发挥着巨大的作用。九、参考文献1王玉英王建国史加荣鲁萍数学建模及其软件实现北京清华大学出版社,20152赵凤群戴芳王小侠肖艳婷数学实验基础西安理工大学理学院2013十、附录程序1马尔萨斯模型的线性解法T01790101980X0395372961291712323143865026297609201065123213171507179320402265PLOTT,X,ON1APOLYFITT0,X0,NYLOGXPPOLY2SYMA程序2人口数量实际值与两种模型预测数据对比图CLEART1790101980X395372961291712323143865026297609201065123213171507179320402265NX139EXP002222T1790NX22859/（17231EXP002858T1790PLOTT,X,R,T,NX1,B,T,NX2,GLEGEND实际值,马尔萨斯模型,阻滞增长模