数学建模人口老龄化问题论文

..数学建模人口老龄化问题论文

摘要我国是世界第一人口大国,地球上每5个人中就有一个中国人.在20世纪的一段时间,我国人口的增长速度过快,有效地控制我国人口的增长,解决养老金制度面临一胎化政策、人口老龄化、及通膨加剧、社保基金收益低等一系列问题，不仅是深入贯彻科学开展观的需要,而且对于全人类的美好理想来说,也是我们义不容辞的责任。本文过查询国家统计局，中国统计年鉴-2011；早报，南方都市报，等各个的数据，并对所查信息进展分析，通过模型建立与分析，以解决以下几个问题；针对问题一我们利用Leslie模型来预测的，首先我们将人按年龄大小等间距的分隔，建立Leslie矩阵，建立Leslie模型，通过Matlab的程序运算，预测中国未来40年的各个年龄段的人口数量，通过Excel分析预测中国未来的40年人口的构造：青少年与中年人的数目相对平衡，而老龄人占大多数，老龄化速度加快；针对问题二我们对所搜集到的数据用最小二乘法来分析用残差平方和来度量测量值与回归直线的接近或偏差程度。利用MatlAB程序求解，并作出养老金累计结余金额实际测量值与回归值的拟合曲线图作出判断，并且运用MatlAB程序预测未来养老金数额，得到结论：未来养老金数额呈直线上升趋势，也就是说未来养老金规模将逐渐变大。针对问题三我们采用拟合曲线来分析未来四十年的人均收入，再根据资料取通货膨胀率k=3%10%15%20%条件下判断未来养老金是否能保障退休水平。针对问题四我们从生育问题、国民经济总值、个人税收三个方面提出未来"养老难〞的可行性方案。针对问题五问题重述一、分析中国未来40年的人口构造和中国老龄化的速度；二、对未来的养老金规模进展预测；三、分析在不同水平的通货膨胀率下未来养老金是否能真正保障退休水平〔收入为当时居民收入的平均收入水平或一些保证生活所需的收入水平〕；四、提出解决未来"养老难〞的可行性方法；五、如果将退休年龄从60岁推迟到65岁将会对中国有什么利弊；问题一一、模型假设1、假设人口的出生率与死亡率不随时间段k变化，只与年龄组有关。2、时间以年为单位，年龄按周岁计算，最大年龄。3、较短时间，平均年龄变化较小，可以认为不变。4、不考虑移民引起人口总人数的变动。二、符号说明1、表示年份〔表示规定的初始年〕2、用表示第年岁〔指满岁但不到岁〕的人数，,

(设为最高年龄)3、表示第年岁女性生育率〔每位女性平均生育的婴儿数〕，育龄区4、表示岁人口的女性比5、表示第年所有育龄女性平均生育的婴儿数6、表示岁人口的死亡率7、表示岁人口的存活率三、模型建立与求解1、模型准备模型主要考虑女性人口数,鉴于男女人口数有一定的比值(可以预测出来),那么由女性人口可以得之总人口数,将女性按年龄顺序分成假设干组，假设每一组的妇女有一样的生育率和死亡率,并且假设最大寿命为90岁(90岁以上按90岁计算)。以年为区间,其中15~49岁年龄段的女性有生育能力。下面建立模型来描述在时刻的性态。第0年龄组在时间的人数为这1年出生并存活的人数，即................................................(1)设在年龄组的人数均正比于前一年龄组在时间的人数,即对．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．(2)其中,表示第个年龄组中至时存活下来并进入年龄组的的人的比例.用矩阵记号,人口增长模型可以表示为……………………(3)其中此模型可以用递推的方式预测各年龄组的人数，进而对人口总量及其年龄构造进展预测，由于人群中只有岁的妇女有生育能力，而最大寿命为90岁，那么矩阵中，显然有．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．〔4〕下面确定动态矩阵的表达式即确定向量，；1，确定经分析计算得．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．〔5〕其中所以将各个年龄值带入即可求得.2，确定分析计算得综合以上的1,2即可求得矩阵关于的函数。用Matlab编程〔附录1-1〕计算2010到2051年的人口总数及各年龄段的总人数〔附录1-2〕，由计算得到的数据通过excel工作表分析整理数据得到2010-2052年的不同年龄阶段〔包括男女〕总人数预测折线图〔图1〕与人口构造预测图〔图2〕及高龄人口占总人口比例趋势预测图〔图3〕图1：图2通过不同年龄阶段〔包括男女〕总人数预测折线图〔图1〕可看出青年人数从2010年的4.57亿增加到2032年3.52亿的人口数量此时到达最顶峰，之后人口数量开场下降，到2050年数量到达3.10亿，总体来说到达平稳状态。中年人的人口数量开展趋势与青少年的类似从2010年的5.59亿增加到2030年的6.86亿到达最顶峰此后开场下降，到2050年数量到达5.47亿，说到达平稳状态。总体来看青年人与老年人的数量在中国未来的40年总体保持平衡状态。总的看青年人与中年人的人数总量开展平稳先增后减最终平稳。通过人口构造预测图〔图2〕可看出2010年各个年龄阶段的人口总数相对平衡而在之后的几年里老年人的人口数量将会有明显的增加，2050年人口将到达4.37亿。综上分析可得结论中国未来的40年人口构造为：青少年与中年人的数目相对平衡，而老龄人占大多数，也就是说中国在未来的40年将进入老龄化社会图3：分析高龄人口占总人口比例趋势预测图〔图3〕可得到以下结论：局局部析可看出2010年〔8.0%〕到2032年〔15.2%〕65岁及65岁以上的人口占全国人口总数的比例从呈增长趋势，但2033年〔15.4%〕到2035年〔15.7%〕这几年65岁及65岁以上的人口占全国人口总数的比例增长速度相对前年来说呈稳定状态。2036〔16.2%〕年到2050年〔24.5%〕年增长速度明显加快，2050年以后增长稳定下来，也就是说在未来的40年65岁及65岁以上的老龄人口会越来越多分析60岁及60岁以上的人口占全国人口总数的比例可知从2010年〔11.5%〕到2026年〔18.4%〕呈增长趋势2027年〔18.9%〕到2030年〔19.5%〕呈稳定增长状态，但以后的几年也就是说从2031年〔20.1%〕到2050年〔30.4%〕增长速度明显加快也就是说在未来的40年60岁及60岁以上的老龄人口会越来越多。综上分析可得结论：中国在未来的40年老龄人口占全国人口总数的比例呈增长趋势，也就是说中国的老龄人数量将会持续的增加，中国在未来的40年将会进入老龄化社会。问题二通过查找统计年鉴等各种资料，找出中国已公布的历年养老金规模如下表：年份养老金规模养老金基金支出养老金累计结余1991215.7173.1144.11992365.8321.9220.61993503.5470.6258.61994707.4661.1304.81995950.1847.6429.819961171.8103.9578.619971337.91251.3682.8199814591511.6587.819991965.11924.9733.520002278.52115.5947.1200124892321.31054.120023171.52842.91608200336803122.12205.520044258.43502.1297520055093.34040.3404120066309.84896.75488.920077834.25964.97391.420089740.27389.69931200911490.88894.412526.1201013419.510554.915365.3单位：亿元图1现将1991——2010年的养老金累计结余金额随时间的数据转换成图标形式如下：图4模型建立与求解：模型假设：1、假设在未来的预测中社会安定，2、假设一回归分析：如果把年代作为横坐标，把历年养老金额作为纵坐标，将这些数据放入直角坐标系中，那么得图4，把这个图称为散点图。从图4可以看出，1991——2000年，养老金累计结余金额缓慢上升，变化趋势不是很明显，但从2000年开场，养老金累计结余金额迅速上升，且数据点根本落在一条呈上升趋势的直线附近。对此，我们分析2000年之后的数据，从2000年往后变量X(设X=当时年代-2000)与Y的关系大致可看作是线性关系，即它们之间的相互关系可以用线性关系来描述。但是由于并非所有的数据点完全落在一条直线上，因此X与Y的关系并没有确切到可以唯一地由一个X值确定一个Y值的程度。其他因素，诸如疾病灾荒、经济危机、人口迁徙、人口出生率或增或减、人口死亡率或增或减等都会影响养老金Y的结果。如果我们要研究X与Y的关系，可以作线性拟合……………

(1)我们称(1)式为回归方程，与是待定常数，称为回归系数。从理论上讲，〔1〕式有无穷多解，回归分析的任务是求出其最正确的线性拟合。二、最小二乘法所谓最小二乘法，就是选择和使最小，即用最小二乘法得到的回归直线是在所有直线中与测量值残差平方和最小的一条。由〔2〕式可知是关于,的二次函数，所以它的最小值总是存在的。下面讨论和的求法。三、正规方程组根据微分方程中求极值得方法可知,取得最小值应满足〔3〕将〔2〕式与〔3〕式结合，可得〔4〕那么〔4〕式称为正规方程组。解这一方程组可得(5〕其中〔6〕式中，称为的协方差之和，称为的协方差之和。那么改写〔1〕式可得或由Matlab编程〔见附录2-1〕求解得出和的值，分别是并画出了养老金累计结余金额实际测量值与回归值的拟合曲线图：现在我们来建立养老金金额与时间的回归关系式：因此，在未来40年的养老金金额符合上述回归关系式，并顺从拟合曲线图的开展趋势。现根据上面的预测公式将2011——2050年的养老金规模作尝试预测，用Matlab编程〔见附录2-1〕得出未来40年的养老金额，预测结果如下表：2011——2050年养老金规模预测数据年份

金额〔亿元〕年份金额〔亿元〕2011769.920312172.22012840.020322242.32013910.220332312.42014980.320342382.520151050.420352452.620161120.520362522.720171190.620372592.920181260.720382663.020191330.820392733.120201400.920402803.220211471.120412873.320221541.220422943.420231611.320433013.520241681.420443083.620251751.520453153.820261821.620463223.920271891.720473294.020281961.820483364.120292032.020493434.220302102.120503504.3将其未来40年的养老金额数据用Excel绘制成折线图如下：由问题1得出中国老龄化的速度在逐年上升，并根据养老金规模预测图可看出未来养老金数额呈直线上升趋势，也就是说未来养老金规模将逐渐变大。问题三：一、模型的准备：通货膨胀是指商品和劳动的货币价格总水平持续明显上涨的过程，通常用消费物价指数〔CPI〕，来衡量物价上涨指数的幅度。CPI指数能够迅速直接地反映居民生活价格趋势，所以常用来衡量通货膨胀。通货膨胀对经济增长有正效应也，有负效应。通货膨胀率是货币超发局部与实际需要的货币量之比,用以反映通货膨胀、货币贬值的程度;而价格指数那么是反映价格变动趋势和程度的相对数。二、模型建立：近20年来我国经济开展迅速，状态良好，人民生活水平显著提高，但经济快速开展的同时也出现了通货膨胀的现象，根据网上查到的数据，得到近20年的通货膨胀率折线图，在通货膨胀的影响下，根据问题一的人口数据预测和问题二的养老金规模，用Excel表计算得出60岁以上的人均养老保险金额。通过拟合预测2011—2050年的人均收入水平，再在不同通货率水平下根据养老金规模预测与人居收入水平解决该问题。三、模型求解根据附录〔3-1〕给出的1991-2009年的人居收入数据，通过Matlab软件程序附录〔3-2〕计算出拟合的四次函数〔得到拟合数据图像〕：的参数得到拟合函数，，预测出2011-2050年所对应的人均收入〔附录〔3-1〕〕根据查阅网上资料通货膨胀率的上下划分以下三种类型：1、温和的或缓行的通货膨胀2、疾驰的或奔腾的通货膨胀3、恶性通货膨胀故取通货膨胀率k=3%10%15%20%，计算在此情况下的养老保险金额(附录3-3)通过Matlab程序〔附录3-4〕得到，并与人均收入水平比照，得出结论：问题四：通过以上的各个问题的分析，我们已经知道了，中国在将来的40年老龄人口急剧增加，人口构造已成定局，且根据网上资料知道中国养老金个人账户缺口巨大，而中国养老金的财务风险也决定着制度改革，因此解决未来"养老难〞的问题刻不容缓。影响养老金金额的因素很多如生育问题、国民经济总值、个人税收等等。针对生育问题：

由问题一知：在未来似是年无论是10-30岁人口数量还是30-60岁的人口数量都是呈下降趋势的，因此解决"养老难〞问题必须对方案生育政策—一胎化政策做出改动，应适当放宽方案生育政策，局部地区鼓励多生。尽量防止两个劳动者赡养一个老人的局面。针对国民经济总值：由问题三知在通货膨胀条件下，未来养老金在一定程度上不能保障退休水平。因此需要提高国民经济总值来提升未来养老金额，完善医疗保障制度，改善社会基金收益低的现状。针对个人税收问题：

提高个人税收抑制通货膨胀，均衡经济开展解决"养老难〞问题。问题五如果将退休年龄60岁推迟到65岁对国家既有利有弊。首先我们来谈谈他的利处：中国人口急剧增长，尤其是老年人是数量，如果推迟退休的年龄，不仅有可以缓解养老金压力进而可以缓解地方政府的财政压力而且可以充分发挥老一代在企业中的经历作用，挖掘人力资源的利用率。为国家创造更多的财富。其次我们来谈谈他的敝处：中国的情况实在特殊，在人口和劳动力这方面，国情确实是特殊的。因为中国在面临老龄化社会挑战的同时，还面临着另一项挑战，即庞大的新增劳动力就业问题。推迟退休年龄，势必影响到新增劳动力的就业，而后者有没有方法通过迅速增加的新工作岗位予以消化，如果失业成为问题，那就不仅仅只是经济问题，同时也是稳定问题了。我们当前强调的"稳增长〞，虽然稳的是GDP增长速度，但深层次里更是在"稳就业〞参考文献【1】启源数学建模高等教育2000【2】志勇，祖櫻

MATLAB教程R2010a

航空航天大学

2010【3】王培辉，袁薇"我国通货膨胀率的最优目标区间几何？"【4】国家统计局中国统计年鉴——2011[EB/OL]中国人口【5】社科院：中国养老金个人账户缺口高达1.7万亿[N]早报2012【6】燕绥中国"未富先老〞养老困境待破解[N]

南方都市报

2011【7】郭慧，天巍Leslie人口增长模型研究科技信息报第23卷第6期：163-165,2007【8】伟我国大陆省份人均生活水平分析民族师专学报第23卷第2期2003附录〔1-1〕：计算程序p=0.4873034;%女性占总人口的比例N=[0.45930.42700.15780.88440.71740.35910.38940.10100.86410.89451.8661.7961.7941.5471.8661.3311.0060.91080.42741.9610.49240.66130.57300.61640.45431.9691.9341.1361.9801.9361.9321.0101.5781.9101.9961.9251.3761.3901.7651.6990.36750.41130.31100.19661.3720.25221.8270.09830.18230.18230.94350.02840.28650.98690.48390.25610.98700.74190.31690.18230.80070.08720.95460.55990.78770.98500.55990.04570.30380.47080.21470.72890.18430.65700.27540.41200.47270.44240.80270.67010.80270.47270.03040.000000.000000.969630.802660.000000.030370.802660.10221];N0=N'/10;%第0年〔2010年〕的女性各个年龄段的人口数(千万〕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];%由2010年原始数据得到的生育率c1=1.295274487*c;%修正后的生育率M=sum(c1');%总合生育率d=zeros(91,1);B=[c1;A];L=[B,d];%构造的lestie矩阵fori=0:1:42X=L^i*N0;%第i年后女性各个年龄段的人口数(千万〕Z=X./p;%第i年在各个年龄段的人口总数预测K(i+1,1)=sum(Z);S1=sum(Z([1:10],:));%第i年0-10岁的总人数D(1,i+1)=S1;S2=sum(Z([11:20],:));%第i年11-20岁的总人数E(1,i+1)=S2;S3=sum(Z([22:31],:));%第i年21-30岁的总人数F(1,i+1)=S3;S4=sum(Z([32:41],:));%第i年31-40岁的总人数G(1,i+1)=S4;S5=sum(Z([42:51],:));%第i年41-50岁的总人数H(1,i+1)=S5;S6=sum(Z([52:61],:));%第i年51-60岁的总人数M(1,i+1)=S6;S7=sum(Z([61:91],:));%第i年60-90岁人数P(1,i+1)=S7;S8=sum(Z([66:91],:));%第i年65-90岁人数Q(1,i+1)=S8;endK%2001-2051的人口总数D%年龄在0-9岁总人数〔包括男女〕E%年龄在10-20岁总人数〔包括男女〕F%年龄在20-30岁总人数〔包括男女〕G%年龄在30-40岁总人数〔包括男女〕H%年龄在40-50岁总人数〔包括男女〕M%年龄在50-60岁总人数〔包括男女〕P%年龄在60岁及60岁以上总人数〔包括男女〕Q%年龄在65岁及65岁以上总人数〔包括男女〕附录〔1-2〕2010到2051年的人口总数及各年龄段的总人数(单位:亿)年份总人口0-9岁10-20岁20-30岁30-40岁40-50岁50-60岁>=60岁>=65岁201013.4481.6372.3402.2392.5271.8621.2061.5561.077201113.5541.6152.3322.1482.6591.8241.2701.5921.112201213.6621.5932.3402.0992.7051.8551.3441.6381.157201313.7631.5852.3272.0292.6591.9681.4201.6841.194201413.8571.5722.2951.9802.6302.0321.5101.7311.233201513.9451.5712.2221.9682.5892.0871.6001.7861.275201614.0301.5802.1031.9862.5432.1511.6701.8511.310201714.1131.5981.9922.0702.4932.1801.7371.9101.343201814.1961.6431.8592.1342.3972.2681.7961.9841.376201914.2801.6931.7152.1962.3212.3841.7942.0621.410202014.3661.7311.6272.2802.2172.4971.8002.1521.451202114.4531.7351.6032.3242.1272.6271.7622.2411.500202214.5411.7471.5822.3162.0782.6711.7922.3321.544202314.6261.7651.5742.3242.0092.6271.9032.4351.603202414.7061.7841.5602.3101.9602.5981.9662.5401.667202514.7781.8011.5602.2791.9482.5572.0202.6331.740202614.8421.8141.5692.2061.9672.5112.0822.7331.815202714.8951.8181.5872.0872.0502.4612.1102.8181.891202814.9371.8131.6321.9782.1132.3672.1952.8571.977202914.9711.7961.6811.8452.1742.2912.3072.9112.066203014.9931.7661.7191.7032.2572.1892.4152.9322.142203115.0031.7241.7221.6152.3012.1002.5403.0192.223203215.0051.6721.7351.5922.2932.0522.5813.1842.291203314.9971.6151.7521.5712.3011.9842.5393.3132.315203414.9861.5591.7711.5632.2881.9362.5113.4362.354203514.9701.5061.7881.5502.2561.9242.4713.5552.362203614.9491.4591.8011.5492.1841.9432.4263.6452.431203714.9271.4191.8061.5592.0672.0252.3773.7772.574203814.9011.3861.8011.5761.9582.0882.2853.8892.684203914.8731.3611.7841.6211.8272.1482.2124.0122.787204014.8411.3441.7551.6701.6862.2302.1134.1192.885204114.8041.3341.7131.7071.5992.2732.0274.2062.955204214.7651.3321.6611.7101.5762.2651.9824.2883.066204314.7231.3351.6051.7231.5552.2731.9164.3493.156204414.6781.3431.5491.7401.5482.2591.8704.3903.254204514.6321.3521.4971.7591.5342.2281.8604.4223.337204614.5811.3641.4491.7761.5342.1561.8794.4203.400204714.5271.3761.4101.7891.5432.0401.9594.4183.455204814.4711.3881.3771.7941.5611.9332.0204.4153.493204914.4091.3991.3521.7891.6051.8042.0784.3943.510205014.3421.4101.3351.7721.6531.6642.1584.3723.518205114.2701.4181.3251.7421.6901.5792.1994.3703.496205214.1911.4221.3231.7011.6941.5562.1904.3423.473附录〔2-1〕x=0:10;y=[947.11054.116082205.5297540415488.97391.4993112526.115365.3]';p=polyfit(x,y,1)

p=1.0e+003*1.4167

-1.3076**=0:0.5:10;>>yy=polyval(p,**);plot(**,yy,'-b',x,y,'.r','MarkerSize',20)title('拟合曲线图')x=11:50;y=1416.7*x-1.3076;y附录〔3-1〕人居收入水平年份金额〔元〕年份金额〔元〕年份金额〔元〕199134020115194203110802119923702012630120321206341993410201376232033134339199446020149187203414919919955302015110242035165276199665020161