人口增长模型.docxVIP

　　篇一：数学建模logistic人口增长模型 　　Logistic人口发展模型  　　一、题目描述  　　建立Logistic人口阻滞增长模型 ，利用表1中的数据分别根据从1954年、1963年、1980年到2005年三组总人口数据建立模型，进行预测我国未来50年的人口情况.并把预测结果与《国家人口发展战略研究报告》中提供的预测值进行分析比较。分析那个时间段数据预测  　　的效果好?并结合中国实情分析原因。 　　表1各年份全国总人口数（单位：千万） 　　二、建立模型  　　阻滞增长模型（Logistic模型）阻滞增长模型的原理：阻滞增长模型是考虑到自然资源、环境条件等因素对人口增长的阻滞作用，对指数增长模型的基本假设进行修改后得到的。阻滞作用体现在对人口增长率r的影响上，使得r随着人口数量x的增加而下降。若将r表示为x的函数r(x)。则它应是减函数。于是有：  　　dx  　　?r(x)x,x(0)?x0  　　dt  　　对r(x)的一个最简单的假定是，设r(x)为x的线性函数，即r(x)?r?sx  　　（1）  　　(r?0,s?0)  　　（2）  　　设自然资源和环境条件所能容纳的最大人口数量长率  　　xm，当x?xm时人口不再增长，即增  　　r(xm)?0，代入（2）式得  　　s?  　　r  　　xm，于是（2）式为  　　x)xm  　　（3）  　　r(x)?r(1?将（3）代入方程（1）得： 　　x?dx 　　??rx(1?)  　　xm?dt  　　?x(0)?x0  　　?  　　解得：  　　（4）  　　x(t)?  　　1?( 　　xmxm 　　?1)e?rtx0  　　（5）  　　三、模型求解  　　用Matlab求解，程序如下： t=1954:1:2005;  　　x=[60.2,61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988,130.756];  　　x1=[60.2,61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988];  　　x2=[61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988,130.756];  　　dx=(x2-x1)./x2; a=polyfit(x2,dx,1);  　　r=a(2),xm=-r/a(1)%求出xm和r  　　x0=61.5;  　　f=inline(xm./(1+(xm/x0-1)*exp(-r*(t-1954))),t,xm,r,x0);%定义函数 plot(t,f(t,xm,r,x0),-r,t,x,+b);  　　title(1954-2005年实际人口与理