系统课设人口问题的仿真

系统仿真课程设计题目专业小组成员用偏微分方程进行人口仿真摘要建立中国人口增长的数学模型，由建立的人口发展的偏微分方程来预测中国未来人口的数量和结构。在预测的基础上，考虑到降低生育率与人口数量和老龄化有着直接的关系，所以在预测人口基础之上，我们进一步拓展对未来人口控制进行研究。即在对人口数量预测的同时对其控制及其优化做出探讨。关键词一、提出背景人口增长预测是随着社会经济发展而提出来的。由于人类社会生产力水平低，生产发展缓慢，人口变动和增长也不明显，生产自给自足或进行简单的以货易货，因而对未来人口发展变化的研究并不重要，根本不用进行人口增长预测。人口是社会经济活动的主体,人口的发展变动趋势,对社会经济发展的影响关系极大,因此人口预测在社会经济实践中占有十分重要的地位。现阶段，中国在享受计划生育政策带来的红利的同时，依然面连着人口结构性失调的严重性问题，而当今社会，经济发展迅速，生产力达到空前水平，这时的生产不仅为了满足个人需求，还要面向社会的需求，所以必须了解供求关系的未来趋势。而人口增长预测是对未来进行预测的各环节中的一个重要方面。准确地预测未来人口的发展趋势，制定合理的人口规划和人口布局方案具有重大的理论意义和实用意义。二、问题重述与分析1）、基本假设1假设本问题所使用的数据均真实有效，具有统计分析价值。2假设本问题所研究的是一个封闭系统，也就是说不考虑我国与其它国家的人口迁移问题。3无重大毁灭性自然灾害和疾病，无战争等暴烈活动，即扰动人口发展的因素只有人。4在对人口进行分段处理时，假设同一年龄段的人死亡率相同，同一年龄段的育龄妇女生育率相同。5生育模式在预测时间内保持不变，并且假设一胎只生一个孩子。6人类的生育观念不发生太大改变，如没有集体不愿生小孩的想法7中国各地各民族的人口政策相同。8人口生存环境为一般常态的自然和社会环境。9中短期内，总和生育率、死亡率和出生性别比不发生大的波动。2）、符号说明R表示人的年龄；RM表示人的最高年龄；T表示时间；FR,T时刻T年龄小于R的人口的分布函数；PR,TDR表示时刻T年龄在区间R,RDR内的人数；PR,T表示人口密度函数；初始密度函数R；R,T表示时刻T年龄R的人的死亡率；T表示总和生育率或生育胎次；FT表示婴儿出生率；HR,T表示年龄为R女性的生育加权因子，称生育模式；NT表示人口总数函数；RT表示平均年龄函数；ST表示平均寿命函数；WT表示老龄化指数。3）、数据来源数据来自。三、整体思路及主要涉及因素本次课设是一个关于人口预测的问题，与其他不同，本问题需要根据中国特殊的国情在特有的现状之下去研究，我们根据对问题的分析并结合实际情况认为对人口产生主要影响的因素有以下四个生育率、死亡率、年龄结构、男女比例。在这里需要说明的是对于人口产生影响的一些因素，如经济发展状况，生态环境情况、已婚夫妇对生育所持的态度、医疗技术的发展等，我们认为它们对人口的增长是通过作用于以上四个指标而间接发挥作用的。而对于诸如战争爆发、疾病流行等突发因素，由于其不可预测性，我们不考虑。在考虑上述主要因素的前提下，我们根据不同的因素所对应的不同情况建立与其最接近的偏微分方程，从而在综合情况下对未来人口的总趋势做出一个预测。以下对四种主要因素做一个简单概述1生育率生育率代表育龄妇女生育人口的能力，从一定意义上讲生育率的高低控制着人口增长率高低，通常来说生育率越高人口增长率越高，所以说生育率是人口增长的源头。生育率的影响因素很多，首先是年龄因素，不同年龄段的育龄妇女的生育率不同，通常20岁至30岁的育龄妇女的生育率最强；此外是地域因素，受政策因素、观念认识、周边环境等影响乡村育龄妇女的生育率高于城市育龄妇女的生育率；还有其它因素的影响，比如大规模疾病会降低育龄妇女的生育率。2死亡率死亡率表示一定时期内一个人口群体中死亡的人数占该人口群体的比值，和生育率一样死亡率的高低同样控制着人口增长率高低，如果说生育率是人口增长的源头，则死亡率是人口增长的汇点。同样影响死亡率的因素很多，首先不同年龄段的死亡率不同，通常老年人和刚出生的婴儿的死亡率较高；从长远来看，随着医疗水平的提高，整个人口群体的死亡率将会成下降趋势；此外一些突发事件，如战争、疾病等，将会使使那一段的人口死亡率大幅度提高。3年龄结构年龄结构反映了总体人口在各年龄段分布情况，年龄结构蕴涵的信息量很大，从其中我们可以实现对很多问题的分析，比如从年龄结构我们可以分析出社会的老年化程度，此外从年龄结构我们可以判断出不同时间段人口出生的情况，比如年龄结构不仅反映了总体人口在各年龄段分布情况，而且考虑到不同年龄段人口生育率、死亡率不同等情况，我们可以在年龄结构中有效反映这些差异。4男女比例男女比例反映了总体人口中男性与女性人数的比较关系，男女比例值能反映出体人口中男性与女性人数是否协调，男女比例主要受男女出生比和男女死亡率的影响，男女出生比正常范围在103107，也就是说出生100个女儿的同时会有103107个男儿出生，但是在现实社会中，女性死亡率低于男性，所以男性与女性人数大致相等，社会维持在一个稳定状态。但目前我国男女出生比超过110，这不仅将导致男女比例失调，还会对人口的预测产生影响，所以在人口预测时必须将男女比例问题考虑进去。四、模型建立1）、人口发展的偏微分方程建立造成人口数量和结构变化的因素不外乎出生、死亡和迁移。为简化，预测期间只考虑自然的出生与死亡，不计迁移等社会因素的影响。为研究任意时刻不同年龄的人口数量，引入人口的分布函数和密度函数。时刻T年龄小于R的人口称为人口分布函数，记FR,T，其中T，（R0），均为连续变量，设F连续、可微，时刻T的人口总数记作NT，最高年龄记作RM，理论推导时设RM，于是对非负非降函数FR,T有1人口密度函数定义2其中，PR,TDR表示时刻T年龄在区间R,RDR内的人数。记R,T为时刻T年龄R的人的死亡率，其含义是，R,TPR,TDR表示时刻T年龄在R,RDR内单位时间死亡的人数。为得到PR,T满足的方程，考察时刻T年龄在R,RDR内的人到时刻TDT的情况。他们中活着的那一部分人的年龄变为RDR1,RDRDR1，这里DR1DT。而在DT这段时间内死亡的人数为R,TPR,TDRDT。于是3变为DR1DT可得4即人口密度函数PR,T的一阶偏微方程，死亡率R,T为已知函数。方程（4）有两个定解条件初始密度函数记为PR,0P0R，是已知函数单位时间出生的婴儿数记为P0,TFT,是婴儿出生率。故所得模型为456方程5求解过程比较复杂，在这儿只给出一种特殊情况求解。社会安定且不太长的时间内，死亡率大致与时间无关，故可近似的设R,TR，可得5的解为7这个解在OTR平面上有一个前线的解释图中对应的对角线RT将OTR面的第一象限分为两部分，TR区域内，PR,T则由未来的生育状况FTR和死亡率S0SR112并取2，N/2，此时有RCR1N213可以看出，提高R1意味着晚婚，而增加N意味着晚育。这样，人口发展方程5和单位时间出生的婴儿数FT的表达式12，构成了我们的连续型人口模型。模型中死亡率函数R,T、性别比函数KR,T和初始密度函数P0R可以有人口统计资料直接得到，或在资料的基础上估计，而生育率T和生育模式HR,T则是可以用于控制人口发展过程的两种手段，T可以控制生育的多少，HR,T可以控制生育的早晚或疏密。在上面的模型分布函数FR,T固然是人口发展过程最完整的描述，但使用起来并不方便。下面定义一些人口指数的及它们与PR,T等数量之间的关系。1）人口总数N（T）152）平均年龄RT153平均寿命ST它表示时刻T出生的人不论活到什么时候，死亡率都是按时刻T的R,T计算，这些人的平均存活时间D16ST实际上是预估寿命。通常说目前平均寿命已达多少岁，是指今年出生婴儿的预估寿命，即S0。根据统计资料得到当前的死亡率R,0后就可以算出S04老龄化指数WT定义17显然，平均年龄RT越大WT越大；对于RT相同的两个国家或地区，平均寿命ST大的，表示健康水平高，一个人能工作的时间一生中占得比例越大，于是老龄化的指数WT较小。五、数据处理与仿真1）参数的确立1初始密度函数R的确定初始密度选择的是2010年的人口数据。之所以没有选择最近的数据，是因为这样更方便对预测数据的检验。人口单位数为千人2死亡率函数R,T的确定在社会和平，不出现大的瘟疫等突发情况下，可认为人口的死亡率基本不变。我们分别对20082012年和世界发达国家的各年龄段分性别的死亡率做加权平均，作为我国未来人口的死亡率。3性别比例函数KR,T的确定在未来一段时期，我们认为各年龄段性别比例基本保持不变。我们根据20082012年的人口统计数据，采取平均计算得出各年龄段的性别比例，将之作为计算所用的性别比例。4生育模式函数HR,T的确定生育模式是对年龄为R的女性的生育加权因子，在社会稳定情况下可认为其基本不变。根据惯例，我们认为生育年龄为1549岁。根据20082012年的统计数据，分别对各年龄段的生育率取平均，然后再做归一化处理，将得到的数据作为计算所需的生育模式函数。2）方程求解先取T15，根据方程（5）和已确定的参数，求出每隔5年的各个年龄段的数据，然后再根据方程（15），（16），（17），（18）分别求出人口总数和老龄化指数。然后改变T的值，重复以上步骤计算。运行RENKOU1M程序分别调整不同时间阶段的T的值，然后再计算求出各个年龄段的数据，人口总数和老龄化指数。将取得的数据作横向比较，选出较优结果。六、计算结果与分析取不同值时的各个时间段的人口总数和老龄化指数（表1）取不同值时的人口总数和老龄化指数最大值（表2）总人口数最大值（千人）老龄化最大值101419E00607088121458E00606759141498E0060648151522E00606346161545E00606216181607E0060602521688E00605856以下是取15和18所对应的人口增长曲线和老龄化指数曲线（由于图形太多，只取其中两个具有代表性的）图1图2根据相关文献，中国目前的值约为15，由表一也可知15时的计算结果与真实的误差是最小的。此时人口的峰值约为1552亿，老龄化指数的峰值为0634618时人口的峰值约为1607亿。也就是说，当每个妇女平均生育的小孩数从15提高到18时，整个国家就会多出生九千万人。所以我们现在不应该也不能放开二胎的限制。由表可知，没有一组数据同时满足人口峰值小于15亿和老龄化峰值小于06，下面我们通过调整生育政策来调整人口结构。七、实验拓展由上述可知，最大老龄化指数是与最大人口数呈反相关的。我们可以分别调整不同时间阶段的T的值，以达到优化人口总数和老龄化指数的目的。下面我们分别计算运行YOUHUA1M2010年T15，2020年T182010年T15，2030年T182010年T15，2030年T20图32010年T152020年T18图42010年T152030年T15图52010年T152030年T20经分析可知，图3和图4均满足要求。它们的人口峰值和老龄化指数峰值分别为9亿、0597和1479亿、0596从峰值上比较，二者并无太大区别。但2100年时，后一幅图的值比前者多出近一个亿，达到125亿。所以第一幅图代表方案更优，即20102020年每个妇女平均生育小孩数控制在15附近，20202100年生育政策可以适度放松，将每个妇女平均生育小孩数控制在17附近，这样人口将逐步优化。八、结论概述由以上仿真结果可以看出，在忽略一些不可控和不可预测的因素条件下，我们可以根据当时所特有的情况确定符合的值，从而预测某年的人口总数和老龄化，对未来人口数量和结构有个了解。而在改变值得取值范围我们也可以对未来人口的生育率有所了解，在此基础之上，可以采取相应的措施对我国人口的发展趋势有所掌控。以上结果表明，在中短期，仍然要坚持严格的人口政策，中长期可以适度放松，但是每个妇女平均生育小孩数仍要控制在20以内，否则人口的峰值会很大，超出环境所能承受的极限。附录程序及其说明1RENKOU1MCLC,CLF,CLEARB15每个妇女一生生的小孩数R75533708827490899889127413101014971381180161247541055957875381312586674111332972238521337356321984初始人口数据X11912501186560116241010817001009540101316010400201047820104019010377901051400102121010347301018840990010897050793690600950495520性别比例SH0005936947840845841871751468000000000出生率SM134035038054068084113163234334562791126320973651582193611395120060男人死亡率SW12802302302603003605107311016028942372512812414402470811107218020女人死亡率NLENGTHRSHSUMSUMSHHBSH/SHSUMFORI1NNVI100/100XI归一化性别比例HIHINVIENDFORI1N计算死亡率SMISMIXIXIXI100SWISWI100SISMISWI/XISISI/1000SISI/4ENDAZEROS19A1,1S1HFORI2N计算以后每年的人口数SI1SI5AI,I1SIENDMZEROS19,20储存未来100年的人口数M,1RFORI220M,IAM,I1ENDZSHSUMM人口总数ZSHMAX,KMAXZSH求最大人口总数及其下标NLZEROS19FORI119NLI,I5I3ENDNNLMZNSUMNFORI120ZNIZNI/ZSHI平均年龄SHMI70I预期寿命WIZNI/SHMI老龄化指数ENDWMAX,PMAXW求W的最大值及其下标FORI120NFI20055IENDSUBPLOT1,2,1NFKZSHMAXPLOTNF,ZSH,NFK,ZSHMAX,RO,SETGCA,FONTSIZE,15TEXTNFK,ZSHMAX,BULLETLEFTARROWTHEMAXVALUEIS,NUM2STRZSHMAX,ATITNF,NUM2STRNFKGRIDONXLABEL年份YLABEL人口总数TITLEIT人口增长曲线B15SUBPLOT1,2,2NFPWMAXPLOTNF,W,NFP,WMAX,RO,SETGCA,FONTSIZE,15TEXTNFP,WMAX,BULLETLEFTARROWTHEMAXVALUEIS,NUM2STRWMAX,ATITNF,NUM2STRNFPGRIDONXLABEL年份YLABEL老龄化指数TITLEIT老龄化指数变化曲线B15ZSH,W2YOUHUA1MCLC,CLF,CLEARB15,B118,每个妇女一生生的小孩数J7从（20105（J1）年改变生育政策R75533708827490899889127413101014971381180161247541055957875381312586674111332972238521337356321984初始人口数据X11912501186560116241010817001009540101316010400201047820104019010377901051400102121010347301018840990010897050793690600950495520性别比例SH0005936947840845841871751468000000000出生率SM134035038054068084113163234334562791126320973651582193611395120060男人死亡率SW12802302302603003605107311016028942372512812414402470811107218020女人死亡率NLENGTHRSHSUMSUMSHHBSH/SHSUMFORI1NNVI100/100XI归一化性别比例HIHINVIENDFORI1N计算死亡率SMISMIXIXIXI100SWISWI100SISMISWI/XISISI/1000SISI/4ENDAZEROS19A1,1S1HFORI2N计算以后每年的人口数SI1SI5AI,I1SIENDMZEROS19,20储存未来100年的人口数M,1RFORI220M,IAM,I1ENDZSHSUMM人口总数ZSHMAX,KMAXZSH求最大人口总数及其下标NLZEROS19FORI119NLI,I5I3ENDNNLMZNSUMNFORI12ZNIZNI/ZSHI平均年龄SHMI70I预期寿命WIZNI/SHMI老龄化指数ENDWMAX,PMAXW求W的最大值及其下标X11912501186560116241010817001009540101316010400201047820104019010377901051400102121010347301018840990010897050793690600950495520性别比例SH0005936947840845841871751468000000000出生率SM134035038054068084113163234334562791126320973651582193611395120060男人死亡率SW12802302302603003605107311016028942372512812414402470811107218020女人死亡率NLENGTHRSHSUMSUMSHHB1SH/SHSUMFORI1NNVI100/100XI归一化性别比例HIHINVIENDFORI1N计算死亡率SMISMIXIXIXI100SWISWI100SISMISWI/XISISI/1000SISI/4ENDAZEROS19A1,1S1HFORI2N计算以后每年的人口数SI1SI5AI,I1SIENDFORIJ20M,IAM,I1ENDZSHSUMM人口总数ZSHMAX,KMAXZSH求最大人口总数及其下标NLZEROS19FORI119NLI,I5I3ENDNNLMZNSUMNFORI120ZNIZNI/ZSHI平均年龄SHMI70I预期寿命WIZNI/SHMI老龄化指数ENDWMAX,PMAXW求W的最大值及其下标FORI120NFI20105IENDSUBPLOT1,2,1NFKZSHMAXPLOTNF,ZSH,NFK,ZSHMAX,RO,SETGCA,FONTSIZE,15TEXTNFK,ZSHMAX,BULLETLEFTARROWTHEMAXVALUEIS,NUM2STRZSHMAX,ATITNF,NUM2STRNFKGRIDONXLABEL年份YLABEL人口总数TITLEIT人口增长曲线SUBPLOT1,2,2NFPWMAXPLOTNF,W,NFP,WMAX,RO,SETGCA,FONTSIZE,15TEXTNFP,WMAX,BULLETLEFTARROWTHEMAXVALUEIS,NUM2STRWMAX,ATITNF,NUM2STRNFPGRIDONXLABEL年份YLABEL老龄化指数TITLEIT老龄化指数变化曲线J20055J1K20055K1达到最大人口总数的年份P20055P1达到最大老龄化指数的年份ZSHMAX,K,WMAX,P,J,ZSH,W3YOUHUA2MCLC,CLF,CLEARB15,B120,B218每个妇女一生生的小孩数J7从（20105（J1）年改变生育政策J19从（20105（J11）年改变生育政策R75533708827490899889127413101014971381180161247541055957875381312586674111332972238521337356321984初始人口数据X11912501186560116241010817001009540101316010400201047820104019010377901051400102121010347301018840990010897050793690600950495520性别比例SH0005936947840845841871751468000000000出生率SM134035038054068084113163234334562791126320973651582193611395120060男人死亡率SW12802302302603003605107311016028942372512812414402470811107218020女人死亡率NLENGTHRSHSUMSUMSHHBSH/SHSUMFORI1NNVI100/100XI归一化性别比例HIHINVIENDFORI1N计算死亡率SMISMIXIXIXI100SWISWI100SISMISWI/XISISI/1000SISI/4ENDAZEROS19A1,1S1HFORI2N计算以后每年的人口数SI1SI5AI,I1SIENDMZEROS19,20储存未来100年的人口数M,1RFORI220M,IAM,I1ENDZSHSUMM人口总数ZSHMAX,KMAXZSH求最大人口总数及其下标NLZEROS19FORI119NLI,I5I3ENDNNLMZNSUMNFORI12ZNIZNI/ZSHI平均年龄SHMI70I预期寿命WIZNI/SHMI老龄化指数ENDWMAX,PMAXW求W的最大值及其下标X11912501186560116241010817001009540101316010400201047820104019010377901051400102121010347301018840990010897050793690600950495520性别比例SH0005936947840845841871751468000000000出生率SM134035038054068084113163234334562791126320973651582193611395120060男人死亡率SW12802302302603003605107311016028942372512812414402470811107218020女人死亡率NLENGTHRSHSUMSUMSHHB1SH/SHSUMFORI1NNVI100/100XI归一化性别比例HIHINVIENDFORI1N计算死亡率SMISMIXIXIXI100SWISWI100SISMISWI/XISISI/1000SISI/4ENDAZEROS19A1,1S1HFORI2N计算以后每年的人口数SI1SI5AI,I1SIENDFORIJ20M,IAM,I1ENDZSHSUMM人口总数ZSHMAX,KMAXZSH求最大人口总数及其下标NLZEROS19FORI119NLI,I5I3ENDNNLMZNSUMNFORI120ZNIZNI/ZSHI平均年龄SHMI70I预期寿命WIZNI/SHMI老龄化指数ENDWMAX,PMAXW求W的最大值及其下标X119125011865601162410108170