leslie人口增长模型模型重点讲解论文

-.z.人口增长预测模型摘要本文建立了我国人口增长的预测模型，对各年份全国人口总量增长的中短期和长期趋势作出了预测，并对人口老龄化、人口抚养比等一系列评价指标进行了预测。最后提出了有关人口控制与管理的措施。模型Ⅰ：建立了Logistic人口阻滞增长模型，利用2中数据，结合网上查找补充的数据，分别根据从1954年、1963年、1980年到2005年三组总人口数据建立模型，进行预测，把预测结果与1《国家人口发展战略研究报告》中提供的预测值进行分析比较。得出运用1980年到2005年的总人口数建立模型预测效果好，拟合的曲线的可决系数为0.9987。运用1980年到2005年总人口数据预测得到2010年、2020年、2033年我国的总人口数分别为13.55357亿、14.18440亿、14.70172亿。模型Ⅱ：考虑到人口年龄结构对人口增长的影响，建立了按年龄分布的女性模型（Leslie模型）：以2中提供的2001年的有关数据，构造Leslie矩阵，建立相应Leslie模型；然后，根据中外专家给出的人口更替率1.8，构造Leslie矩阵，建立相应的Leslie模型。首先，分别预测2002年到2050年我国总人口数、劳动年龄人口数、老年人口数（见附录8），然后再用预测求得的数据分别对全国总人口数、劳动年龄人口数的发展情况进行分析，得出：我国总人口在2010年达到14.2609亿人，在2020年达到14.9513亿人，在2023年达到峰值14.985亿人；预测我国在短期劳动力不缺，但须加强劳动力结构方面的调整。其次，对人口老龄化问题、人口抚养比进行分析。得到我国老龄化在加速，预计本世纪40年代中后期形成老龄人口高峰平台，60岁以上老年人口达4.45亿人，比重达33.277%；65岁以上老年人口达3.51亿人，比重达25.53%；人口抚养呈现增加的趋势。再次，讨论我国人口的控制，预测出将来我国育龄妇女人数与生育旺盛期育龄妇女人数，得到育龄妇女人数在短期将达到高峰，随后又下降的趋势的结论。最后，分别对模型Ⅰ与模型Ⅱ进行残差分析、优缺点评价与推广。关键词Logistic人口模型Leslie人口模型人口增长预测MATLAB软件§1、问题重述一、背景知识：中国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。我国人口发展经历了多个阶段，近年来中国的人口发展出现了一些新的特点，例如，老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，这些都影响着中国人口的增长。全面建设小康社会时期是我国社会快速转型期，人口发展面临着前所未有的复杂局面，人口安全面临的风险依然存在二、相关数据：1《国家人口发展战略研究报告》2人口数据（《中国人口统计年鉴》中的部分数据）及其说明根据已有数据三、要解决的问题：1、试从中国的实际情况和人口增长的上述特点出发，参考2中的相关数据（也可以搜索相关文献和补充新的数据），建立中国人口增长的数学模型，并由此对中国人口增长的中短期和长期趋势做出预测；特别要指出你们模型中的优点与不足之处。2、利用所建立模型的预测结果，参照1的相关叙述对反映中国人口增长特点的一系列指标如人口老龄化、人口抚养比等进行分析预测。3、根据模型的计算结果，对未来人口发展高峰进行预测并针对中国人口的调控和管理进行分析。§2、问题分析人口的变化受到众多方面因素的影响，因此对人口的预测与控制也就十分复杂，很难在一个模型中综合考虑到各个因素的影响。为了更好的解决此问题，我们分析了题目以及附录1中所给的相关信息，考虑到可以根据对人口增长不同的评价指标及不同的时期建立多个模型分别加以讨论。一、从1中，我们看到过去一些专家对中国的总人口数做出了2010年、2020年分别达到13.6亿人和14.5亿人，2033年前后达到峰值15亿人左右的预测。因而，我们也可以先对总人口的增长趋势做出自己的预测与专家预测数据进行比较，对于预测所要用到的一些相关数据，我们作了相应的补充，由此我们建立了模型Ⅰ：阻滞增长模型。二、模型Ⅰ只考虑了人口总数，对人口总数进行了预测分析。但实际中在对人口进行分析时，按年龄段分布的人口结构是非常重要的。在人口总数一定时，不同年龄段的人的生育率和死亡率是不同的，它们对人口未来发展的影响也是很不一样的。为了讨论不同年龄段的人口分布对人口增长的影响，我们依据2建立了模型Ⅱ：按年龄分布的Leslie模型。三、由模型Ⅰ和模型Ⅱ的结果我们预测了人口总数的发展趋势，由模型Ⅱ的计算结果我们还能够得到各年份处在各年龄段的人口数量、男女比率的预测值。根据这些预测值我们可以计算出反映人口增长特点的其他指标，由此我们可以对模型的计算结果进行进一步的分析。§3、合理的假设1、社会稳定，不会发生重大自然灾害和战争不随时间而变化2、超过90岁的妇女（老寿星）都按90岁年龄计算3、在较短的时间，平均年龄变化较小，可以认为不变4、不考虑移民对人口总数的影响§4、名词解释与符号说明一、名词解释1、总和生育率——指一定时期（如*一年）各年龄组妇女生育率的合计数，说明每名妇女按照*一年的各年龄组生育率度过育龄期，平均可能生育的子女数，是衡量生育水平最常用的指标之一。2、更替水平——指这样一个生育水平，同一批妇女生育女儿的数量恰好能替代她们本身。一旦达到生育更替水平，出生和死亡将逐渐趋于均衡，在没有国际迁入与迁出的情况下，人口将最终停止增长，保持稳定状态。3、人口抚养比——指人口总体中非劳动年龄人口数与劳动年龄人口数之比。通常用百分比表示。说明每100名劳动年龄人口大致要负担多少名非劳动年龄人口。用于从人口角度反映人口与经济发展的基本关系。根据劳动年龄人口的两种不同定义（15-59岁人口或15-64岁人口），计算总抚养有两种方式4、人口老龄化——指人口中老年人比重日益上升的现象。促使人口老龄化的直接原因是生育率和死亡率降低，主要是生育率降低。一般认为，如果人口中65岁及以上老年人口比重超过7%，或60岁及以上老年人口比重超过10%，则该人口就属于老年型。5、出生人口性别比——是活产男婴数与活产女婴数的比值，通常用女婴数量为100时所对应的男婴数来表示。正常情况下，出生性别比是由生物学规律决定的，保持在103～107之间。二、符号说明序号符号意义1：表示年份（选定初始年份的）2人口增长率3：人口数量4：自然资源和环境条件所能容纳的最大人口数量5：可决系数6：在时间段第年龄组的人口总数7：第年龄组的生育率8：第年龄组的死亡率9：第年龄组的存活率10：Leslie矩阵11：2001年全国人口总数12：2001年城市总人口13：2001年镇总人口14：2001年乡总人口15：2001年第年龄段的人口总数16：时分别表示市、镇、乡的女孩出生率17：时段具有劳动能力的人口18：社会的抚养比指数19:总和生育率20：时段年龄组中女性所占的百分比§5、模型的建立与求解模型Ⅰ：阻滞增长模型（Logistic模型）[1]一、模型的准备阻滞增长模型的原理：阻滞增长模型是考虑到自然资源、环境条件等因素对人口增长的阻滞作用，对指数增长模型的基本假设进行修改后得到的。阻滞作用体现在对人口增长率的影响上，使得随着人口数量的增加而下降。若将表示为的函数。则它应是减函数。于是有：（1）对的一个最简单的假定是，设为的线性函数，即（2）设自然资源和环境条件所能容纳的最大人口数量，当时人口不再增长，即增长率，代入（2）式得，于是（2）式为（3）将（3）代入方程（1）得：（4）解方程（4）可得：（5）二、模型的建立为了对以后一定时期的人口数做出预测，我们首先从中国经济统计数据库（211.86.245.155/inde*.asp*）上查到我国从1954年到2005年全国总人口的数据如表1。表1各年份全国总人口数（单位：千万）年份195419551956195719581959196019611962总人口60.261.562.864.666.067.266.265.967.3年份196319641965196619671968196919701971总人口69.170.472.574.576.378.580.783.085.2年份197219731974197519761977197819791980总人口87.189.290.992.493.795.096.25997.598.705年份198119821983198419851986198719881989总人口100.1101.654103.008104.357105.851107.5109.3111.026112.704年份199019911992199319941995199619971998总人口114.333115.823117.171118.517119.850121.121122.389123.626124.761年份1999200020012002200320042005总人口125.786126.743127.627128.453129.227129.988130.7561、将1954年看成初始时刻即，则1955为，以次类推，以2005年为作为终时刻。用函数（5）对表1中的数据进行非线性拟合，运用Matlab编程（程序见附录1）得到相关的参数，可以算出可决系数（可决系数是判别曲线拟合效果的一个指标）：由可决系数来看拟合的效果比较理想。所以得到中国各年份人口变化趋势的拟合曲线：（6）根据曲线（6）我们可以对2010年（）、2020年（）、及2033年（）进行预测得（单位：千万）：结果分析：从附录1所给信息可知从1951年至1958年为我国第一次出生人口高峰，形成了中国人口规模“由缓到快”的增长基础；因此这段时期人口波动较大，可能影响模型结果的准确性。1959、1960、1961年为三年自然灾害时期，这段时期人口的增长受到很大影响，1962年处于这种影响的滞后期，人口的增长也受到很大影响。总的来说1951-1962年的人口增长的随机误差不是服从正态分布，由于上面的曲线拟合是用最小二乘法，所以很难保证拟合的准确性。因此我们再选择1963年作为初始年份对表1中的数据进行拟合。2、将1963年看成初始时刻即，以2005年为作为终时刻。运用Matlab编程（程序见附录2）得到相关的参数，可以算出可决系数得到中国各年份人口变化趋势的另一拟合曲线：（7）根据曲线（7）我们可以对2010年（）、2020年（）、及2033年（）进行预测得（单位：千万）：结果分析：1963年-1979年其间，人口的增长基本上是按照自然的规律增长，特别是在农村是这样，城市受到收入的影响，生育率较低，但都有规律可寻。总的来说，人口增长的外界大的干扰因素基本上没有，可以认为这一阶段随机误差服从正态分布；1980-2005年这一时间段，虽然人口的增长受到国家计划生育政策的控制，但计划生育的政策是基本稳定的，这一阶段随机误差也应服从正态分布（当然均值与方差可能不同）因此用最小二乘法拟合所得到的结果应有较大的可信度。3、从1980-2005年，国家计划生育政策逐渐得到完善及贯彻落实，这个时期的人口增长受到国家计划生育政策的控制，人口的增长方式与上述的两个阶段都不同。因此我们进一步选择1980年作为初始年份2005年作为终时刻进行拟合。运用Matlab编程（程序见附录3）得到相关的参数，可以算出可决系数得到中国各年份人口变化趋势的第三条拟合曲线：（8）根据曲线（7）我们可以对2010年（）、2020年（）、及2033年（）进行预测得（单位：千万）：结果分析：这一时期，国家虽然对人口大增长进行了干预，但国家的计划生育的政策是基本稳定的，在此其间没有其他大的干扰，所以人口增长的随机误差应服从正态分布。所以我们的结果应是比较可信的。我们分别根据拟合曲线（6）、（7）、（8）对各年份中国总人口进行预测得到结果如表2：表2各年份全国总人口用不同拟合曲线预测数（单位：千万）年份全国总人口预测（单位：千万）预测曲线（6）预测曲线（7）预测曲线（8）2000126.7649126.3338126.4732003130.5141129.2303129.51682006134.1131.8447132.27582009137.516134.1926134.76382012140.7577136.2917136.99712015143.8231138.1607138.99332018146.7117139.819140.7712021149.4251141.2856142.34892024151.9662142.579143.74522027154.3392143.7168144.97782030156.5494144.7157146.06322033158.6028145.5908147.01722036160.5063146.3562147.85412039162.267147.0247148.58712042163.8924147.6077149.22842045165.3903148.1158149.78862048166.7683148.558150.2775由上表可以看出：用拟合曲线（6）预测得到的数据比较大，在2024年总人口就已经超过了151.9662千万，而且一直以比较快的速度增长到2048年达到了166.7683千万。用拟合曲线（7）预测得到的数据偏小，到2048年人口只有148.558千万。相比较而言用拟合曲线（8）预测的数据比较接近1中的预测。画出图形如图1：图1：对各年份全国总人口数的预测模型Ⅱ：按年龄分布的Leslie模型[2]一、模型的准备将人口按年龄大小等间隔地划分成个年龄组（譬如每10岁一组），模型要讨论在不同时间人口的年龄分布，对时间也加以离散化，其单位与年龄组的间隔相同。时间离散化为.设在时间段第年龄组的人口总数为，定义向量，模型要研究的是女性的人口分布随的变化规律，从而进一步研究总人口数等指标的变化规律。设第年龄组的生育率为，即是单位时间第年龄组的每个女性平均生育女儿的人数；第年龄组的死亡率为，即是单位时间第年龄组女性死亡人数与总人数之比，称为存活率。设、不随时间变化，根据、和的定义写出与应满足关系：（9）在（9）式中我们假设中已经扣除婴儿死亡率，即扣除了在时段以后出生而活不到的那些婴儿。若记矩阵（10）则（9）式可写作（11）当、已知时，对任意的有（12）若（10）中的元素满足（ⅰ）；（ⅱ），且至少一个。则矩阵称为Leslie矩阵。只要我们求出Leslie矩阵并根据人口分布的初始向量，我们就可以求出时段的人口分布向量。二、模型的建立我们以2001年为初始年份对以后各年的女性总数及总人口数进行预测，根据2中所给数据，以一岁为间距对女性分组。(1)计算2001年处在各个年龄上的妇女人数的分布向量：2给了2001年中国人口抽样调查数据，提取为表3表3城市男147907城市女147465镇男80279镇女77976乡男394690乡女372242根据抽样调查的结果，可以算出2001年城市、镇、乡人口占2001年全国总人口的比率分别为：我们由表1数据知2001年全国总人口（单位：千万），因此可以算出2001年城市、镇、乡的总人口分别为（单位：千万）：、、根据2给的2001年城市、镇、乡各个年龄段的女性比率，可以分别算出2001年城市、镇、乡处在第年龄段的女性的总数分别为。以城市为例，设2001年城市中处在年龄段妇女占城市总人口比率分别为，则（镇、乡类似）。于是可以算出2001年处在第年龄段上的妇女总人数（见附录7）。（2）计算处在第年龄段的每个女性平均生育女儿的人数。2中分别给出了2001年城市、镇、乡育龄妇女（15岁—49岁）的生育率（此处应该是包含男孩和女孩）（或时都为0），则可以分别算出2001年处在第年龄段的城市、镇、乡育龄妇女总共生育的小孩数（包含男孩和女孩），记为：。以城市为例计算：`（镇、乡类似）。2中还分别给出了2001年市、镇、乡的男女出生人口性别比（女100计），据此可以分别计算出城市、镇、乡女孩的出生率。由此就可以求出2001年处在第年龄段的每个女性平均生育女儿的人数：，由于总和生育率：经计算得到总和生育率小于1.8，误差很大，我们对生育率进行修正：具体计算结果见附录7。(3)计算第年龄段的女性总存活率率：记第年龄段的女性的死亡率为。2中分别给出了城市、镇、乡处在第年龄段的女性死亡率，则处在第年龄段的女性总死亡率为：，于是总存活率为：见附录4。用E*CEL对计算出来的数据进行整理，然后运用MATLAB软件进行编程，计算出Leslie矩阵，于是可以用上面(12)式进行预测。三、对模型结果作进一步讨论我国人口发展形势复杂，目前人口的低生育水平面临着严峻的挑战，下面我们分别从如下方面分析预测我国人口发展将要面临的复杂局面。（1）人口总量与劳动力人口的发展变化根据考虑种群结构的Leslie离散模型，利用2001年的数据建立人口预测模型。通过分析，计算出我国人口的预测值，对应作出的我国劳动年龄人口与总人口的折线图如下：图2我国全国总人口与劳动年龄人口折线图根据图2可以知道从2001年到2023年预测我国全国总人口是呈现上升趋势的，随后几年呈现缓慢下降的趋势。总人口在2010年、2020年分别达到14.2609亿人和14.9513亿人，在2023年达到峰值14.985亿人，在2033年达到14.7455亿人。把预测数值与2中所提供的预测数值进行比较，发现我们预测的未来人口的高峰期提前10年。这一方面可能由我国男女的出生性别比例中女性所占的比例较小的原因；另一方面，我们计算出人口更替率仅为1.42（此为5年的均值），而中外专家对我国90年代中期以来的人口更替率的计算结果为1.8（见附录10），两者相差甚远，这说明附录---提供的数据可能不够真实，从而导致了我国人口峰值的预测年份提前。根据图2，我国劳动年龄人口庞大，15-64岁的劳动年龄人口2010年为10.4421亿人，2013年将达到高峰10.4852亿人，随后劳动年龄人口呈现下降的趋势。由此，可知在相当长的时间，我国不缺劳动力，但需要加强劳动力结构性的调整，同时由于我国计划生育等宏观政策的影响，近几年总和生育率已降低到1.8，并将稳定在1.8的水平上，所以经过较长的时期，我国的劳动年龄人口将有所降低。（2）人口老龄化与人口抚养比通过计算分析人口结构持续老龄化，运用Leslie离散模型，通过MATLAB软件计算出我国60岁以上与65岁以上的老龄人口数，做出散点图如下：图3我国老年人口预测值的折线图从图3可以直观的看出我国老龄人口在持续增加，说明我国老龄化进程在加速。同时做出未来我国老龄人口占总人口的比例的折线图如下：图4我国老龄人口占总人口预测比例的折线图从图3，图4得到：2001年我国60岁以上老年人口已达到1.5538亿人，占总人口的11.5693%。到2020年，60岁以上老年人口将达到2.907亿人比重为19.443%；65岁以上老年人口将达到2.0628亿人比重从2000年的8.009%增长到13.797%。预计本世纪40年代中后期形成老龄人口高峰平台，60岁以上老年人口达4.45亿人，比重达33.277%；65岁以上老年人口达3.51亿人，比重达25.53%。综上可知我国老龄人口数量大，老龄化速度快，高龄趋势明显，加上我国人口基数大，所以我国是个老龄人口多的国家。老龄化也在一定程度上导致了我国人口抚养比的不断增高。下面计算人口抚养比指数：设与分别为男性与女性中具有劳动能力的年龄组，则时段具有劳动能力的人口为,而为时段由社会抚养的失去劳动能力与老人或尚未具有劳动能力的为成年人的数量。定义社会的抚养比指数，即平均每一劳动者抚养的无劳动能力的人数。我们以0—14岁为没有劳动能力的儿童，以15-64岁为具有劳动能力的年龄劳动人口，以65岁及以上的为老龄人口。首先，通过MATLAB编程计算出2002到2051年0-14岁、15-64岁、65岁及5以上三段的人数；其次，根据人口抚养比的含义，计算出每一年份的人口抚养比得出人口抚养比。得出的每年人口抚养比的折线图如下：图5预测人口抚养比从图5可以看出预测的以后各年的人口抚养比呈增长的趋势。人口抚养比比较高主要原因有：每年新生婴儿数目在增加；老龄化的加剧，老龄人口数量大；15-64岁年龄段中的人的残疾、生病而无劳动能力等。（3）人口调控与管理现阶段我国生育水平的不稳定性，根据建立的Leslie模型，运用MATLAB软件计算出2000年到2050年我国育龄妇女（15-49岁）人口，并做出的散点图如下：图6未来我国育龄妇女（15-49岁）人口预测从图6中可以看出我国育龄妇女（15-49岁）人口在2010年左右到达到高峰，图7未来我国生育旺盛期育龄妇女（20-29）人数预测从图7我们发现，我国生育旺盛期育龄妇女（20-29）人数在2012年将达到高峰，到2025年左右有进入一个小低谷，然后再2037年左右有达到一个小高峰。第二个我国生育旺盛期育龄妇女（20-29）人数小高峰的原因在于在2012年人口出生高峰期的女婴到2037年时达到生育旺盛期，因此，在2025年生育旺盛期育龄妇女（20-29）人数达到低谷时有回升的形势。§6、误差分析与灵敏度分析一、模型的残差分析：1、运用Matlab软件计算出用1954年到2005年的总人口数进行拟合产生的残差,再利用E*CEL作出残差的散点图如下：图8残差分析从图8可以看出残差在坐标轴上下波动，但是，不是呈现正态分布，并且残差绝对值之和为57.9992，是比较大，因此拟合的效果不太好。2、利用1963年到2005年的总人口数，根据Logistic模型的形式，用Matlab软件进行拟合，并求出残差序列，再利用E*CEL进行处理，并作出残差散点图如下：图9残差分析图通过图9，可以看出残差值大致分布在坐标轴的上下，呈现对称分布，又有Matlab软件计算出拟合的残差绝对值之和为27.8046，因此效果较好。3、利用1980年到2005年的人口总数居，同样运用Matlab、E*CEL软件进行分析、处理，作出散点图如下：图10残差分析图通过Matlab软件计算，得出拟合的残差绝对值之和为10.1699，从图10可以看出，图形基本关于坐标轴对称，所以你和效果比较好。二、灵敏度分析：1、在不同的总合生育率下按照前面的方法分别计算从2001年到2050年全国人口总数的预测值（程序见附录6），并画出图形如图11图11：在不同的k值下对各年份全国总人口数的预测由图11可以看出当值很小时人口增长比较缓慢，达到峰值后人口数量很快下降出现严重负增长；当值很大时人口增长速度很快，达到峰值后下降的速度缓慢，在此情况下人口数量急剧膨胀。只有当值适中时，总人口增长才比较稳定。2、再在不同的总和生育率下按照前面的方法分别计算从2001年到2050年全国老龄化变化趋势（程序见附录6），并画出图形如图12图12：在不同的k值下对各年份老龄化变化趋势由图12可以看出值越小，老龄化增大的速度越快；值越大老龄化指数增长平缓年龄结构稳定，有利于社会发展。由以上分析可知国家在制定人口政策时要多方面考虑，如果只看重对人口总数的控制可能导致社会老龄化严重、劳动力不足这显然是不利于社会经济发展的；相反如果为了防止社会老龄化加快而放任人口的增长，也会导致社会人口过多对资源和环境带来巨大压力。因此只有掌握好一个“平衡点”正确制定政策才能使国民经济持续增长，人民生活水平不断提高。§7、模型的评价与推广一、模型的优点：1、在用模型Ⅰ对各年全国人口总数预测时结合实际情况，分别用不同时间段的数据拟合确定了三个预测函数。并对三个函数预测的数据进行了对比分析，使模型的计算结果更加准确。2、利用E*CEL软件对数据进行处理并作出各种平面图，简便，直观、快捷；3、运用多种数学软件进行计算，取长补短，使计算结果更加准确;4、在模型Ⅱ中我们充分考虑到不同年龄的个体具有不同的生育能力和死亡率，采用leslie模型，建立年龄结构的离散模型，并通过合理假设，在时间跨度不大的前提下，对人口数量仅此进行了预测，得到人口数量变化趋势图2与<<国家人口发展战略研究：人口发展预测>>课题中未来我国总人口，劳动人口及人口扶养比预测及未来我国人口老龄化预测趋势图基本一致。因为原始数据得到的人口总和生育率跟实际情况不符，我们对此进行了合理修正，使预测更为准确。在模型Ⅰ中我们还进行了参差分析，在模型Ⅱ中我们对不同的平均妇女生育胎数下人口总数及老龄化趋势进行了分析，得到适合平均生育胎数的最佳值。二、模型的缺点：在模型假设中我们及不随时段的变迁而改变这一理想状态下，但出生率及死亡率会随时间的变化而有所该变，本模型没有建立与死亡率随时间变化的动态模型，因而存在一定的误差；三、模型的改进：随着人民的生活水平的提高和医疗卫生的改善，各年龄的死亡率不断下降，存活率不断提高。因此我们可以对Leslie模型进行进一步改变：记时段年龄组中女性所占的百分比为，并设为育龄女性的年龄组，则时段新生儿为我们引入控制变量，使得=1，这里，，称为女性生育模式，我们将lestie矩阵变成：其中在一定时期（这里j从0到90），为平均生育胎数，和可视为与无关的常数，我们可以通过控制结婚年龄和生育两胎间的年龄差来求的最佳值，从而达到控制人口数量和年龄结构的目的。四、模型的推广：本文首先不考虑年龄结构对人口增长的影响，建立Logistic人口预测模型；然后，逐步改进，考虑年龄结构对人口增长的影响，建立Leslie模型，对人口增长进行预测，这种由简到繁，逐步加深的思路，可以应用到较复杂问题的处理上。参考文献[1]启源,金星,叶俊.数学模型[M].:.2003年8月第三版；[2]启源.数学模型[M].:高等教育.1987年4月第一版；[3]于洪彦.E*cel统计分析与决策[M].:高等教育.2006年4月；[4]胡守信,柏年.基于MATLAB的数学实验[M].:科学.2004年6月；[5]扬启帆,康旭升,等.数学建模[M].:高等教育.2006年5月；[6]于学军.《中国人口科学》2000年第2期,时间:2000-4-6,中国人口信息网.附录附录1：t=0:51;%令1954年为初始年*=[60.261.562.864.66667.266.265.967.369.170.472.574.576.378.580.78385.287.189.290.992.493.79596.25997.598.705100.1101.654103.008104.357105.851107.5109.3111.026112.704114.333115.823117.171118.517119.85121.121122.389123.626124.761125.786126.743127.627128.453129.227129.988130.756];[c,d]=solve('c/(1+(c/60.2-1)*e*p(-5*d))=67.2','c/(1+(c/60.2-1)*e*p(-20*d))=90.9','c','d');%求初始参数b0=[241.9598,0.02985];%初始参数值fun=inline('b(1)./(1+(b(1)/60.2-1).*e*p(-b(2).*t))','b','t');[b1,r1,j1]=nlinfit(t,*,fun,b0)y=180.9871./(1+(180.9871/60.2-1).*e*p(-0.0336.*t));%非线性拟合的方程plot(t,*,'*',t,y,'-or')%对原始数据与曲线拟合后的值作图R1=r1.^2;R2=(*-mean(*)).^2;R=1-R1/R2%可决系数W=sum(abs(r1))%残差绝对值之和附录2：t=46:3:94y=180.9871./(1+(180.9871/60.2-1).*e*p(-0.0336.*t))%对总人口进行预测t=0:42;%令1963年为初始年*=[69.170.472.574.576.378.580.78385.287.189.290.992.493.79596.25997.598.705100.1101.654103.008104.357105.851107.5109.3111.026112.704114.333115.823117.171118.517119.85121.121122.389123.626124.761125.786126.743127.627128.453129.227129.988130.756];[c,d]=solve('c/(1+(c/69.1-1)*e*p(-5*d))=78.5','c/(1+(c/69.1-1)*e*p(-20*d))=103.008','c','d');%求初始参数b0=[134.368,0.056610];%初始参数值fun=inline('b(1)./(1+(b(1)/69.1-1).*e*p(-b(2).*t))','b','t');[b1,r1,j1]=nlinfit(t,*,fun,b0)y=151.4513./(1+(151.4513/69.1-1).*e*p(-0.0484.*t));%非线性拟合的方程plot(t,*,'*',t,y,'-or')%对原始数据与曲线拟合后的值作图R1=r1.^2;R2=(*-mean(*)).^2;R=1-R1/R2%可决系数W=sum(abs(r1))%残差绝对值之和附录3：t=37:3:85y=151.4513./(1+(151.4513/69.1-1).*e*p(-0.0484.*t))%对总人口进行预测t=0:25;%令1980年为初始年*=[98.705100.1101.654103.008104.357105.851107.5109.3111.026112.704114.333115.823117.171118.517119.85121.121122.389123.626124.761125.786126.743127.627128.453129.227129.988130.756];[c,d]=solve('c/(1+(c/98.705-1)*e*p(-5*d))=105.851','c/(1+(c/98.705-1)*e*p(-8*d))=111.026','c','d');%求初始参数b0=[109.8216,-0.19157];%初始参数值fun=inline('b(1)./(1+(b(1)/98.705-1).*e*p(-b(2).*t))','b','t');[b1,r1,j1]=nlinfit(t,*,fun,b0)y=153.5351./(1+(153.5351/98.705-1).*e*p(-0.0477.*t));%非线性拟合的方程plot(t,*,'*',t,y,'-or')%对原始数据与曲线拟合后的值作图R1=r1.^2;R2=(*-mean(*)).^2;R=1-R1/R2%可决系数W=sum(abs(r1))%残差绝对值之和t=20:3:53y=153.5351./(1+(153.5351/98.705-1).*e*p(-0.0477.*t))%对总人口进行预测附录4：计算0-14岁，15-64岁，65岁及以上的程序、绘画出未来我国育龄人数的程序N=[0.6808912720.584591720.5845582070.6922202170.724110210.7755360410.8473689180.8344187030.9179220420.9514668191.0700157171.2492560631.1992639881.2021985251.2742189171.1110508390.9923144250.8937975440.8746573470.9843568770.8595767780.852153460.908644180.8979448070.8805393231.0190867241.042186671.1148237311.1928671991.2035665721.2729739951.3285135761.2549924031.3338194451.1031861231.224703071.2206434421.2367363191.3907264150.9807651110.6466840690.7856606230.7016275920.9104201120.9601576460.9142587130.9539805680.9274299560.8510077590.8254823590.8079428230.7365520020.690432040.605802950.6155106240.5547856630.503701350.4800517620.4687228170.4553640590.4843865410.4473446810.4201644980.442380330.4265290910.4281838750.391329530.3804091290.3853399670.3279245740.3346977110.3073300120.2628648340.2706631830.2358721650.2087254950.2120015490.1784567720.1642603160.1498428330.1387349160.1098999490.0973582770.07657620.06381350.0557941230.0493960160.03828810.0335447770.0238706160.070211606];N0=N';%第0年（2001年）的女性个年龄段的人口数A=eye(90);b=[0.9749069660.9993212310.997724330.9992476160.9995674180.9991806630.9998879480.9993875960.9996185860.9999856720.9993894340.9997243540.9998017960.9996276260.9997047950.9996396860.9997284620.9999745330.9991733270.9989541180.9994410670.9993573920.9992906750.9989991760.9998816040.9988963470.9983559390.9991353390.9990745270.9988726520.9991807940.9989181590.9990461120.9990423540.9993960270.9986249720.9982527160.9995978550.9987109450.9990032740.9994434440.9991414150.9987721010.9989405050.9979050050.9983745620.9977837740.9975966660.9973449060.9969544990.9966697840.9960307590.9950066390.9961574880.9946477440.9957794350.9956523130.995777130.9924778060.9949695640.9881305370.9892848680.9887039610.9883025630.984208240.9844954160.9852987350.9800620890.9789283070.9773584460.9711269890.9693038990.9699798180.964050590.9617403120.967297060.9483023460.9465715590.9496413870.9359493910.9124894820.92618050.9237578630.9287579060.9182303330.8877613890.8853068580.8751780860.8824957520.824428701];fori=1:90A(i,:)=A(i,:)*b(1,i);endA;c=[0000000000000004.478E-050.0003221690.0003582460.0010046040.0046833670.0110111650.0336164920.0578753940.0748717270.0691820060.0760391410.067248950.0524294060.0437324640.0343505020.0246327330.0232525320.0183438470.0147012750.0110399610.0071175570.0050948430.003592910.0025148580.0024847810.0017647090.0014716440.0006769530.0002654760.0004014740.0004087790.0001104470.0001924010.0003894210.0002240690000000000000000000000000000000000000000];c1=1.295274487*c;M=sum(c1');%总合生育率d=zeros(91,1);B=[c1;A];L=[B,d];%构造的lestie矩阵fori=0:49H=L^i*N0;%第i年人口总数Q(1,i+1)=sum(H([16:50],:));%第i年15-49育龄妇女总数P(1,i+1)=sum(H([21:30],:));%第i年20-29生育旺盛期妇女总数end*=2001:2050;y1=Q*10;y2=P*10;plot(*,y1,'*')plot(*,y2,'-or')gridon附录5：计算2001到2051年的人口总数程序p=0.464429182;%女性占总人口的比例N=[0.6808912720.584591720.5845582070.6922202170.724110210.7755360410.8473689180.8344187030.9179220420.9514668191.0700157171.2492560631.1992639881.2021985251.2742189171.1110508390.9923144250.8937975440.8746573470.9843568770.8595767780.852153460.908644180.8979448070.8805393231.0190867241.042186671.1148237311.1928671991.2035665721.2729739951.3285135761.2549924031.3338194451.1031861231.224703071.2206434421.2367363191.3907264150.9807651110.6466840690.7856606230.7016275920.9104201120.9601576460.9142587130.9539805680.9274299560.8510077590.8254823590.8079428230.7365520020.690432040.605802950.6155106240.5547856630.503701350.4800517620.4687228170.4553640590.4843865410.4473446810.4201644980.442380330.4265290910.4281838750.391329530.3804091290.3853399670.3279245740.3346977110.3073300120.2628648340.2706631830.2358721650.2087254950.2120015490.1784567720.1642603160.1498428330.1387349160.1098999490.0973582770.07657620.06381350.0557941230.0493960160.03828810.0335447770.0238706160.070211606];N0=N'/10;%第0年（2001年）的女性各个年龄段的人口数(千万）N00=N0/10%把单位化成亿（人）A=eye(90);b=[0.9749069660.9993212310.997724330.9992476160.9995674180.9991806630.9998879480.9993875960.9996185860.9999856720.9993894340.9997243540.9998017960.9996276260.9997047950.9996396860.9997284620.9999745330.9991733270.9989541180.9994410670.9993573920.9992906750.9989991760.9998816040.9988963470.9983559390.9991353390.9990745270.9988726520.9991807940.9989181590.9990461120.9990423540.9993960270.9986249720.9982527160.9995978550.9987109450.9990032740.9994434440.9991414150.9987721010.9989405050.9979050050.9983745620.9977837740.9975966660.9973449060.9969544990.9966697840.9960307590.9950066390.9961574880.9946477440.9957794350.9956523130.995777130.9924778060.9949695640.9881305370.9892848680.9887039610.9883025630.984208240.9844954160.9852987350.9800620890.9789283070.9773584460.9711269890.9693038990.9699798180.964050590.9617403120.967297060.9483023460.9465715590.9496413870.9359493910.9124894820.92618050.9237578630.9287579060.9182303330.8877613890.8853068580.8751780860.8824957520.824428701];fori=1:90A(i,:)=A(i,:)*b(1,i);endA;c=[0000000000000004.478E-050.0003221690.0003582460.0010046040.0046833670.0110111650.0336164920.0578753940.0748717270.0691820060.0760391410.067248950.0524294060.0437324640.0343505020.0246327330.0232525320.0183438470.0147012750.0110399610.0071175570.0050948430.003592910.0025148580.0024847810.0017647090.0014716440.0006769530.0002654760.0004014740.0004087790.0001104470.0001924010.0003894210.0002240690000000000000000000000000000000000000000];%由2001年原始数据得到的生育率c1=1.295274487*c;%修正后的生育率M=sum(c1');%总合生育率d=zeros(91,1);B=[c1;A];L=[B,d];%构造的lestie矩阵fori=0:1:50*=L^i*N0;%第i年后女性各个年龄段的人口数(千万）Z=*./p;%第i年在各个年龄段的人口总数预测K(i+1,1)=sum(Z);S1=sum(Z([1:15],:));%第i年0-14岁的总人数D(1,i+1)=S1;S2=sum(Z([16:65],:));%第i年15-64岁的总人数S3=sum(Z([61:91],:));%第i年60-90岁人数G(1,i+1)=S3;E(1,i+1)=S2;S4=sum(Z([66:91],:));%第i年65-90岁人数F(1,i+1)=S4;endK%2001-2051的人口总数D%年龄在0-14岁总人数（包括男女）E%年龄在15-64岁总人数（包括男女）F%年龄在65岁及65岁以上总人数（包括男女）G%年龄在60岁及60岁以上总人数（包括男女）附录6：functionW=pare(*)p=0.464429182;%女性占总人口的比例N=[0.6808912720.584591720.5845582070.6922202170.724110210.7755360410.8473689180.8344187030.9179220420.9514668191.0700157171.2492560631.1992639881.2021985251.2742189171.1110508390.9923144250.8937975440.8746573470.9843568770.8595767780.852153460.908644180.8979448070.8805393231.0190867241.042186671.1148237311.1928671991.2035665721.2729739951.3285135761.2549924031.3338194451.1031861231.224703071.2206434421.2367363191.3907264150.9807651110.6466840690.7856606230.7016275920.9104201120.9601576460.9142587130.9539805680.9274299560.8510077590.8254823590.8079428230.7365520020.690432040.605802950.6155106240.5547856630.503701350.4800517620.4687228170.4553640590.4843865410.4473446810.4201644980.442380330.4265290910.4281838750.391329530.3804091290.3853399670.3279245740.3346977110.3073300120.2628648340.2706631830.2358721650.2087254950.2120015490.1784567720.1642603160.1498428330.1387349160.1098999490.0973582770.07657620.06381350.0557941230.0493960160.03828810.0335447770.0238706160.070211606];N0=N';A=eye(90);b=[0.9749069660.9993212310.997724330.9992476160.9995674180.9991806630.9998879480.9993875960.9996185860.9999856720.9993894340.9997243540.9998017960.9996276260.9997047950.9996396860.9997284620.9999745330.9991733270.9989541180.9994410670.9993573920.9992906750.9989991760.9998816040.9988963470.9983559390.9991353390.9990745270.9988726520.9991807940.9989181590.9990461120.9990423540.9993960270.9986249720.9982527160.9995978550.9987109450.9990032740.9994434440.9991414150.9987721010.9989405050.9979050050.9983745620.9977837740.9975966660.9973449060.9969544990.9966697840.9960307590.9950066390.9961574880.9946477440.9957794350.9956523130.995777130.9924778060.9949695640.9881305370.9892848680.9887039610.9883025630.984208240.9844954160.9852987350.9800620890.9789283070.9773584460.9711269890.9693038990.9699798180.964050590.9617403120.967297060.9483023460.9465715590.9496413870.9359493910.9124894820.92618050.9237578630.9287579060.9182303330.8877613890.8853068580.8751780860.8824957520.824428701];b1=[0.9749069660.9993212310.997724330.9992476160.9995674180.9991806630.9998879480.9993875960.9996185860.9999856720.9993894340.9997243540.9998017960.9996276260.9997047950.9996396860.9997284620.9999745330.9991733270.9989541180.9994410670.9993573920.9992906750.9989991760.9998816040.9988963470.9983559390.9991353390.9990745270.9988726520.9991807940.9989181590.9990461120.9990423540.9993960270.9986249720.9982527160.9995978550.9987109450.9990032740.9994434440.9991414150.9987721010.9989405050.9979050050.9983745620.9977837740.9975966660.9973449060.9969544990.9966697840.9960307590.9950066390.9961574880.9946477440.9957794350.9956523130.995777130.9924778060.9949695640.9881305370.9892848680.9887039610.9883025630.984208240.9844954160.9852987350.9800620890.9789283070.9773584460.9711269890.9693038990.9699798180.964050590.9617403120.967297060.9483023460.9465715590.9496413870.9359493910.9124894820.92618050.9237578630.9287579060.9182303330.8877613890.8853068580.8751780860.8824957520.8244287010.7717624];fori=1:90A(i,:)=A(i,:)*b(1,i);endA;c1=[0000000000000004.478E-050.0003221690.0003582460.0010046040.0046833670.0110111650.0336164920.0578753940.0748717270.0691820060.0760391410.067248950.0524294060.0437324640.0343505020.0246327330.0232525320.0183438470.0147012750.0110399610.0071175570.0050948430.003592910.0025148580.0024847810.0017647090.0014716440.0006769530.0002654760.0004014740.0004087790.0001104470.0001924010.0003894210.0002240690000000000000000000000000000000000000000];%由2001年原始数据得到的生育率t=sum(c1);c=((**p-t)/t+1)*c1%修正后的生育率M=sum(c');%总合生育率d=zeros(91,1);B=[c;A];L=[B,d];%构造的lestie矩阵[V,d]=eig(L);%求特征根与特征向量p=d(42,42);%特征根Q=-V(:,42);%对应的正特征向量fori=0:49D=L^i*N0;%第i年女性人口分布E(i+1,1)=sum(D)/p%第i年总人口（2001为第0年）forj=0:90%大于90岁的按90岁算F(j+1,1)=j*D(j+1,1)/p;T(j+1,1)=e*p(-b1(1,j+1));endY(i+1)=sum(F)/E(i+1,1);%平均年龄endY%输出01-50年平均年龄矩阵T=0;s=0;fori=0:90%大于90岁的按90岁算T=T+e*p(b1(1,j+1)-1);%求平均寿命，不随年份而变化endTW=Y/T;%社会老龄化指数*=2001:2050;W1=pare(1.6);W2=pare(1.8);W3=pare(2.0);W4=pare(2.2);plot(*,W1,'-r')holdonplot(*,W2,'-G')plot(*,W3,'-B')plot(*,W4,'-Y')附录7：年龄女性总人数（万人）出生婴儿数平均有孩子数女性死亡率存活率0680.8912721000.0250930340.9749069661584.5917197000.0006787690.9993212312584.5582068000.002275670.997724333692.2202166000.0007523840.9992476164724.1102102000.0004325820.9995674185775.5360408000.0008193370.9991806636847.368918000.0001120520.9998879487834.4187027000.0006124040.9993875968917.9220422000.0003814140.9996185869951.4668194001.43283E-050.999985672101070.015717000.0006105660.999389434111249.256063000.0002756460.999724354121199.263988000.0001982040.999801796131202.198525000.0003723740.999627626141274.218917000.0002952050.999704795151111.0508390.049752824.478E-050.0003603140.99963968616992.31442550.319692670.0003221690.0002715380.99972846217893.79754350.320199030.0003582462.54672E-050.99997453318874.65734670.8786842150.0010046040.0008266730.99917332719984.35687724.6101047520.0046833670.0010458820.99895411820859.57677859.4649417930.0110111650.0005589330.99944106721852.153460128.64641020.0336164920.0006426080.99935739222908.644180352.588139950.0578753940.0007093250.99929067523897.944807567.230678280.0748717270.0010008240.99899917624880.539323360.917476650.0691820060.0001183960.999881604251019.08672477.490478720.0760391410.0011036530.998896347261042.186