数学建模-人口增长模型

1、人口增长模型摘要本文根据某地区的人口统计数据，建立模型估计该地区2010年的人口数量。首先，通过直观观察人口的变化规律后，我们假设该地区的人口数量是时间的二次函数，建立了一个二次函数模型，并用最小二乘法对已有数据进行拟合得到模型的具体参数，从而可以预测2010年的人口数为333.8668百万。 然后，我们发现从1980年开始该地区的人口增长明显变慢，于是我们假设人口增长率是人口数的线性减函数，即随着人口数的增加，人口的增长速度会慢慢下降，从而我们建立了阻滞增长模型，利用此模型我们最后求出2010年的人口预报数为296.3865。关键字：人口预报，二次函数模型，阻滞增长模型问题重述：根据某地区人

2、口从1800年到2000年的人口数据（如下表），建立模型估计出该地区2010年的人口 ，同时画出拟合效果的图形。表1 该地区人口统计数据年 份1800181018201830184018501860人口7.213.817.217.624.733.636.2年 份1870188018901900191019201930人口48.658.173.389.8105.6125.9149.1年 份1940195019601970198019902000人口172.2189.8230.5246.7262.1271.2280.3符号说明 时刻的人口数量 初始时刻的人口数量 人口增长率 环境所能容纳的最大人口

3、数量，即问题分析首先，我们运用软件1编程（见附件1），绘制出1800年到2000年的人口数据图，如图1。图1 1800年到2000年的人口数据图从图1我们可以看出1800年到2000年的人口数是呈现增长的趋势的，而且类似二次函数增长。 所以我们可以建立了一个二次函数模型，并用最小二乘法对已有数据进行拟合得到模型的具体参数。 于是我们假设人口增长率是人口数的线性减函数，即随着人口数的增加，人口的增长速度会慢慢下降,从而我们可以建立一个阻滞增长模型。模型建立模型一：二次函数模型我们假设该地区时刻的人口数量的人口数量是时间的二次函数，即：我们可以根据最小二乘法，利用已有数据拟合得到具体参数。即，要求

4、、和，使得以下函数达到最小值：其中是时刻该地区的人口数，即有：令，可以得到三个关于、和的一次方程，从而可解得、和。我们用编程（见附件2），解得0.006018，,，即：从而我们可以预测2010年的人口数为百万。图2 二次函数模型的拟合效果图图2是所得到的二次函数模型和原数据点的拟合效果图。 从图2可以看出拟合的效果在1950年之前还可以，但是对后期的数据拟合的不好。模型二：阻滞增长模型我们假设人口增长率是人口数的线性减函数，即随着人口数的增加，人口增长速度会慢慢下降：人口数量最终会达到饱和，且趋于一个常数，当时，增长率为0：由上面的关系式可得出：把上式代进指数增长模型的方程中，并利用初始条件，

5、可以得到：解得：我们可以利用已有数据拟合求解得（程序见附件4）:, -0.027958。可以预测2010年的人口数为百万。图4 阻滞增长模型的拟合效果图图4是阻滞增长模型的拟合效果图。 从图4我们可以看出我们的模型对该地区的人口数据拟合得很好。可以看出阻滞增长模型更客观地反映人口的增长规律，基本上都在拟合曲线上，拟合效果好，特别是后期的数据非常的吻合，所以次模型对未来的人口数预测是很适合的，结果更准确，对未来的预测比指数增长模型更为优越。参考文献1 刘卫国, 陈昭平, 张颖. MATLAB程序设计与应用M, 北京:高等教育出版社, 2002年。2 姜启源, 谢金星, 叶俊. 数学建模(第三版)

6、M, 北京：高等教育出版社, 2004年。附录附件1：1800年到2000年的人口数据图x=1800:10:2000;y=7.2 13.8 17.2 17.6 24.7 33.6 36.2 48.6 58.1 73.3 89.8 105.6 125.9 149.1 172.2 189.8 230.5 246.7 262.1 271.2 280.3;figure;plot(x,y,'r*');附件2:线性增长模型的拟合代码x=1800:10:2000;y=7.2 13.8 17.2 17.6 24.7 33.6 36.2 48.6 58.1 73.3 89.8 105.6 125

7、.9 149.1 172.2 189.8 230.5 246.7 262.1 271.2 280.3;plot(x,y,'r*'); % 画点，红色hold on; % 使得以下图形画在同一个窗口p = polyfit(x,y,2) % 多项式拟合，返回系数pxn = 1800:5:2010; % 定义新的横坐标yn = polyval(p,xn); % 估计多项式p的函数值plot(xn,yn) % 把(x,yn)定义的数据点依次连起来% 给图形加上图例xlabel('年份');ylabel('人口数');legend('原始数据

8、9;,'拟合函数',2);box on; grid on;x1=2010;y1 = polyval(p,x1) % 估计多项式p在未知点的函数值附件3：阻滞增长模型的拟合代码clc; % 清屏幕clear; % 清除以前的变量% 数据点(t,y)t=1800:10:2000;y=7.2 13.8 17.2 17.6 24.7 33.6 36.2 48.6 58.1 73.3 89.8 105.6 125.9 149.1 172.2 189.8 230.5 246.7 262.1 271.2 280.3;plot(t,y,'b*');% 定义需要拟合的函数类型my

9、fun(a,t)，a是参数列表,t是变量 myfun = (a,t)a(1)./(1+(a(1)./7.2-1)*exp(a(2)*(t-1800); a0=500,1; % 初始值 % 非线性拟合.最重要的函数，第1个参数是以上定义的函数名，第2个参数是初值，第3、4个参数是已知数据点a=lsqcurvefit(myfun,a0,t,y);disp('a=' num2str(a); % 显示得到的参数% 画出拟合得到的函数的图形ti=1800:10:2010;yi=myfun(a,ti);hold on;plot(ti,yi,'r');% 给图形加上图例xlabel('年份');ylabel('