人口增长模型 (2)

上传人:jin****ng 文档编号:122410388 上传时间:2022-07-20 格式:DOC 页数:13 大小:223KB
收藏 版权申诉 举报 下载
人口增长模型 (2)_第1页
第1页 / 共13页
人口增长模型 (2)_第2页
第2页 / 共13页
人口增长模型 (2)_第3页
第3页 / 共13页
资源描述:

《人口增长模型 (2)》由会员分享,可在线阅读,更多相关《人口增长模型 (2)(13页珍藏版)》请在装配图网上搜索。

1、2009高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网 上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的 资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规 则的行为,我们将受到严肃处理。我们参赛选择的题号是(从A/B/C/D中选择一项填写):B我们的参赛报名号为(如果赛区设置报

2、名号的话):所属学校(请填写完整的全名):西安理工大学参赛队员(打印并签名):1.2.3.指导教师或指导教师组负责人(打印并签名):日期:2011年7月3日赛区评阅编号(由赛区组委会评阅前进行编号)2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号(由赛区组委会评阅前进行编号)赛区评阅记录(可供赛区评阅时使用)评阅人评分备注全国统一编号(由赛区组委会送交全国前编号)全国评阅编号(由全国组委会评阅前进行编号)人口增长模型、摘要本文根据某个地区从1982年到1998年间的人口统计数据,对其做出相应的分析, 做出一些合理的假设,利用Mat lab软件画出该地区历年的人口数据图,对数据图进行

3、 仔细观察后,认识人口数量的变化规律通过曲线的插值和拟合,我认为该地区的人口呈 现三次函数增长,并建立了三次函数的数学模型对未来人口的预测。但是,考虑到,自 然资源、环境条件等对人口增长起着阻滞作用,并且随着人口的增加,阻滞作用越来越 大,于是,建立了阻滞增长模型。最后,把两种预测的结果的人口数量与已知的人口数 量进行比较,求出其误差,检验模型的正确性。关键词:人口增长模型 曲线插值 曲线拟合 阻滞增长模型二、问题重述下表列出了中国1982-1998年的人口统计数据,取1982年为起始年(t 0 ), 1982年的人口 101654万人,人口自然增长率为千分之14,以36亿为我国人口的容纳量,

4、 试建立一个较好的人口数学模型并给出相应的算法和程序,并于实际人口进行比较。年198219831984198519861987198819891990人口(万)101654103008104357105851107507109300111026112704114333年19911992199319941995199619971998人口(万)115823117171118517119850121121122389123626124810三、基本假设模型一假设:(1) 假设不存在某抽样年龄段出现0死亡概率(2) 假设人口平稳增长,无大型自然灾害,战争等因素的影响(3) 假设境内外迁移率对我国未来

5、人口影响不计(4) 人口的净增长率(即出生率减去死亡率)为常数(5) 时刻t的人口函数是连续可微的;模型二假设:(1) 假设不存在某抽样年龄段出现0死亡概率;(2) 人口的净增长率(即出生率减去死亡率)随人口的增长而单调递减;(3) 考虑人口的增长受到环境的因素及资源的影响不能无限的增长;(4) 时刻t的人口函数是连续可微的;四、符号与说明为i点的横坐标值;为第i点的纵坐标的值;r增长率xo今年人口xkk年后人口x (t)t时刻人口xm人口容量s确定系数五、问题分析5.1名词解释:曲线插值:根据已知若干点的函数值来求出其他点的函数值的方法,常用的插 值方法有:拉格朗日插值,牛顿插值,分段线性插

6、值,三次样条插值。曲线拟合:根据若干有观察和测试得到的数据进行拟合,给出逼近曲线。最常 用的拟合的方法是最小二乘法。5.2问题的总体分析:根据1982年到1998的人口数据统计,在matlab中进行描点,然后根据这些点 的变化规律选择合适的方法及参数进行曲线插值拟合得到数学模型对未来人口的预测。 又考虑实际的环境因素对人口增长的影响,建立了阻滞增长模型对未来人口的预测。5.3问题的准备:插值问题设函数y=f(x)在区间a,b连续,若给定n+1个点aW x x 0, s 0)这里r称作固有增长率。为了确定系数s的意义,引入自然资源和环境条件所能容纳的最大数量x,称人口容量。当x = x时人口不再

7、增长,即增长率mmr(x )= 0,代入(6)为:mr (x )= r (1 皂)(7)xm(7)式代入方程(5)得:空=rx(l 丄),x (0 )= x (8) dtx0m方程(8)右端的因子rx体现人口自身的增长趋势,因子(1一丄)则体现了资xm源和环境对人口增长的阻滞作用。显然x越大,前一因子越大,后一因子越小,人口增长是两个因子共同作用那个的结果。七、模型求解通过matlab编写程序最后得到合理的数学模型为:模型一求解:77906944*xT-5900627182481437/268435456*x+7632250728293203/524288模型二求解:由得到的数学模型对1982

8、1998年人口计算的结果为:单位(万)年份198219831984198519861987模型预测101470102970104510106090107690109310模型二预测101654102678103709104746105788106837年份198819891990199119921993模型预测110920112520114100115650117160118610模型二预测107892108953110020111092112170113254年份19941995199619971998模型预测119990121290122510123620124630模型二预测114344

9、115439116539117646118757八、结果分析预测结果与统计结果的差值如下:年份198219831984198519861987统计结果101654103008104357105851107507109300模型一与统计差值-184-3815323918310模型二与统计差值0-330-648-1105-1719-2463年份198819891990199119921993统计结果111026112704114333115823117171118517模型一与统计差值-106-184-233-173-1193模型二与统计差值-3134-3751-4313-4731-5001-52

10、63年份19941995199619971998统计结果119850121121122389123626124810模型一与统计差值140169121-6-180模型一与统计差值-5506-5682-5850-5980-6053误差分析:公式为:1 皆 I y1 一 y ID1 =(为)*100%170 y模型一误差:Dl=0.12%模型二误差:D2=3.0768%通过分析可知,该人口增长模型,在误差允许的范围内是比较准确的。九、模型的评价与推广本文是首先将人口增长理想化以后做出的合理简单的人口三次函数增长模型,没有 考虑实际人口的影响因素:出生率,死亡率,迁入人口,迁出人口等,认为人口的增长

11、 率是固定不变的,有一定的局限性。随后,综合考虑各个能影响人口增长的因素,增长 率随人口的增长递减,建立了阻滞增长模型对人口进行预测,通过比较两种模型对未来 人口增长进行预测。事实上,人口的预测是一个相当复杂的问题,影响人口增长的因素除了人口基数与 可利用资源量外,还和医药卫生条件的改善、人们生育观念的变化等因素有关,特别在 做中短期预测时,我们希望得到满足一定预测精度的结果,比如在刚刚经历过战争或是 由于在特定的历史条件下采纳了特殊的人口政策等,这些因素本身以及由此而引起的人 口年龄结构的变动就会变的相当重要,进而需要必须予以考虑。参考文献1秦新强,数学建模,高等学校教材,【M】2010,

12、P45-65送来忠,王志明,数学建模与实验,【M】科学出版社,2003,P75-81。 姜启源等 数学模型(第三版,【M】高等教育出版社2003.08 P10附录程序一:% curvefit02.mt=1982:1998;y=101654 103008 104357 105851 107507 109300 111026 112704 114333115823 117171 118517 119850 121121 122389 123626 124810; plot(t,y,.k,markersize,20) %描点a=polyfit(t,y,3)%求出二次多项式的系数p=poly2sym(a)给出多项式的符号表达式xi=1982:0.1:1998;yi=polyval(a,xi);%求多项式在分点的值hold onplot(xi,yi,b-,markersize,20)%绘制拟合曲线图yii=interp1(t,y,xi,liner); % 求线性插值函数值 plot(xi,yii,r-);%绘制插值函数曲线 yj=polyval(a,t);%求出拟合函数的函数值 yjj=yj-y; %预测结果与统计结果的程序二:y=dsolve(Dy-0.014*(y-yA2/360000),y(0)=101654,t)

展开阅读全文
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!