关于河南人口的预测模型.docVIP

关于河南人口的预测模型 新乡学院 数学与信息科学系实验报告 实验名称 河南省人口预测模型 所属课程 数学软件与实验 实验类型 开放性实验 专 业 信息与计算科学 班 级 学 号 姓 名 指导教师 一、实验概述 【实验目的】 [1] 学习由实际问题去建立数学模型的全过程; [2] 训练综合应用数学模型、微分方程、函数拟合和预测的知识分析和解决实际问题; [3] 应用matlab软件求解微分方程、作图、函数拟合等功能，设计matlab程序来求解其中的数学模型; [4] 提高论文写作、文字处理、排版等方面的能力; 通过完成该实验，学习和实践由简单到复杂，逐步求精的建模思想，学习如何建立反映人口增长规律的数学模型，学习在求解最小二乘拟合问题不收敛时，如何调整初值，变换函数和数据使优化迭代过程收敛。 【实验原理】 非线性拟合函数用lsqcurvefit格式。 人口预测模型的建立过程中,主要考虑两方面的问题,一方面是人口预 测参数的选择和假定。 【实验环境】 MatlabR2010b 二、实验内容 问题 查阅1980-2010 年河南省人口数据，用查得数据检验马尔萨斯(Malthus)人口指数增长模型，根据检验结果进一步讨论马尔萨斯人口模型的改进，并利用这两种模型来预测河南省2013年的人口数量. 1. 基本假设及符号说明 我们已知Malthus模型的假设是基于人口增长率为常数,而后设符号。 符号假设, r 人口增长率,常数,; x(t) t时刻的人口数; x0 为初始时刻的人口, 1 dx,,,rx得到微分方程如下, ,dt ,x(0),x0, 改进后的模型为Logistic模型,充分的的考虑到资源环境的容纳能力,人口到 达一定数量后,增长率会下降,不再为常数,于是得到改进的微分方程, 符号变量, r 为固有增长率, xm 为在资源环境有限的条件下人口最大容量, r(x)=r(1-x/xm) 为改进后的人口增长率 改进微分方程 2. 问题求解 根据Malthus 模型的基本假设, 我们可以分别求得微分方程的解析解:p1=x0*exp(r*x); 编写拟合函数的文件fun1.m function y=fun1(x,t) p1=x(1).*exp(x(2).*t); 运行的脚本文件 clear all clc x0=[7285,0.002]; x=[1980:1:2010]; p=[7285 7397 7519 7632 7732 7737 7847 7985 8148 8317 8491 8649 8763 8861 8946 9027 9100 9243 9315 9387 9488 9555 9613 9667 9717 9768 9820 9869 9918 9967 10437] a=lsqcurvefit(fun1,x0,(x-1979),p) x1=[1980:1:2010]; p1=fun1(a,(x-1979)); plot(x,p,ro,x1,p1,b) p2013=fun1(a,(2013-1979)) 根据Logistic模型，我们可以分别求得微分方程的解析解:p2= xm/(1+x0*exp(-r*x))编写拟合函数的文件fun2.m function y=fun(x,t) y=17000./(1+x(1).*exp(-x(2).*)); 2 运行的脚本文件 clear all clc x0=[50,0,02]; x=[1980:1:2010]; p=[7285 7397 7519 7632 7732 7737 7847 7985 8148 8317 8491 8649 8763 8861 8946 9027 9100 9243 9315 9387 9488 9555 9613 9667 9717 9768 9820 9869 9918 9967 10437] plot(x,p,ro) hold on a=lsqcurvefit(fun2,x0,(x-1979),p) x1=[1980:1:2010]; p2=fun2(a,(x1-1979)); plot(x1,p2) p2013=fun(a,(2013-1979)) 3. 结论及分析 Malthus模型的最小二程曲线拟合,在2013年,p2013 = 1.0752e+04 10500 10000 9500 9000 8500 8000 7500 70001980198519901995200020052010 3 Logistic模型的最小二乘曲线拟合,在2013年p2013 =1.0633e+04 10500 1