养老金建模.doc

企业退休职工养老金制度的改革模型评价与改进摘要本文主要通过建立保持替代率下的养老金收支平衡模型来探讨企业退休职工养老金制度的改革方案，我们首先建立logistic模型预测出2011年到2035年山东省职工的年平均工资，然后建立替代率、收支平衡模型来分析各种情况下的替代率和缺口情况，最后通过灵敏度分析对养老金制度的改革方案提出合理的建议。针对问题一：我们根据题中给出的数据和条件，结合中国最近几年经济发展增长趋势和未来经济战略目标，假设未来中国经济发展呈现阻滞指数增长模式，工资增长形势与国家经济发展形势相同，建立符合未来经济发展趋势的增长阻滞logistic模型，用Matlab编程求解得到平均工资与时间的关系：，并进行非线性二乘拟合方法预测2011到2035山东省职工平均工资（见表1）。针对问题二：我们首先计算出该各年龄段职工工资与企业平均工资之比作为职工缴费指数的参考值，然后建立养老金替代率模型进而得出不同情况下的养老金替代率。 针对问题三：我们建立养老金缺口模型，然后分别对（55岁 ，60岁，65岁）退休并从退休后一直领取养老金到75岁3种情况下的养老金缺口，运用 Excel软件进行求解，并分别对3种退休情况的收支平衡计算得到： 时间30-55岁30-60岁30-65岁缺口平衡时间 关键字：养老金 Logistic模型 替代率 收支平衡模型 灵敏度分析 Matlab软件1、问题的重述目前，我国企业职工基本养老保险实行“社会统筹”与“个人账户”相结合的模式。养老金的发放与职工在职时的工资及社会平均工资有着密切关系，工资的增长又与经济增长相关。我国经济发展的战略目标，是要在21世纪中叶使我国人均国民生产总值达到中等发达国家水平。在此条件下，我国企业退休职工养老金制度也要相应的改革和完善。附件中给出了山东省职工历年平均工资数据 、2009年山东省某企业各年龄段职工的工资分布情况以及养老金的计算办法。现在我国养老保险改革正处于过渡期。养老保险管理的一个重要的目标是养老保险基金的收支平衡，它关系到社会稳定和老龄化社会的顺利过渡。本文尝试解决以下问题。问题一：对未来中国经济发展和工资增长的形势做出简化、合理的假设，并参考附件1，预测从2011年至2035年的山东省职工的年平均工资。问题二：根据附件2计算2009年该企业各年龄段职工工资与该企业平均工资之比。把这些比值看作职工缴费指数的参考值，考虑该企业职工自2000年起分别从30岁、40岁开始缴养老保险，一直缴费到退休（55岁，60岁，65岁），计算各种情况下的养老金替代率。问题三：假设该企业某职工自 2000年起从30岁开始缴养老保险，一直缴费到退休（55岁，60岁，65岁），并从退休后一直领取养老金，至75岁死亡。计算养老保险基金的缺口情况，并计算该职工领取养老金到多少岁时，其缴存的养老保险基金与其领取的养老金之间达到收支平衡。问题四：如果既要达到目标替代率，又要维持养老保险基金收支平衡，提出可行的措施，并给出理由。2、模型假设假设1：题目附件中给出的数据真实可靠。假设2：养老金会随着社会平均工资的调整而调整。如果职工死亡，社会统筹账户中的资金不退给职工，个人账户中的余额可继承。假设3：个人账户储存额以银行当时公布的一年期存款利率计息，为简单起见，利率统一设定为3%。 3、符号说明：年份：年份山东省职工的年平均工资：各月收入范围的中间值，：每个月收入范围的总的职工人数，：企业的平均工资：各个年龄段，：工资为、年龄段为的企业人数，：各年龄段职工平均工资与企业平均工资的比值，：职工平均缴费指数：替代率：企业和职工实际缴纳基本养老保险费的年限：参保人员i年度的本人缴费工资，：各年利率：计发月数4、问题分析4.1 问题一的分析对于该问题我们首先根据附录一用Matlab画出1978到2010年的散点图， 得出经济增长的大体趋势。又通过查找文献未来中国经济增长预测【1】得知在未来经济增长率会下降，但总体经济仍保持上升趋势。这与人口增长的Logistic模型相似。建立未来经济预测的Logistic模型,先用最小二乘拟合方法拟合出改革开放100年内山东省职工的平均水平，分析是否符合Logistic模型然后再拟合1978年到2035年内的山东省职工工资水平。4.2问题二的分析因为各年龄段职工的具体工资未明确给出，所以本文假设每个月收入范围内的工资呈均匀分布，取每个工资段中间值作为每个年龄段中每个工资段的职工工资。再求每个年龄段的平均工资作为各年龄段职工工资。通过计算2009年该企业各年龄段职工工资与该企业平均工资之比并把这些比值看作职工缴费指数的参考值，再求其平均值，利用附录三进一步求得替代率。4.3问题三的分析 对于该问题我们只需求出职工从30岁到退休时所缴纳的基金累计包括社会统筹与个人账户的基金。再求出职工从退休到其75岁时所领的养老金总数，两者之差即为养老保险的缺口。令两者相等即可得到收支平衡时的年龄。4.4问题四的分析我们首先参考问题二和问题三 ，然后根据题目所给数据得到2009年30岁开始缴费到60岁，领取养老金到75时的替代率和养老金缺口，然后对各因素对替代率和养老金缺口的影响进行逐步回归分析。5、模型的建立与求解5.1问题一的模型建立与求解2011年前三季度，经济增长率逐季下降，全年总增长率与2010年10.3%的增长率相比，已现下降趋势。现假设未来25年内，经济增长率总体趋势下降，但总体经济仍保持上升趋势，即呈现“S型”增长。这与人口增长的Logistic模型的增长性相似。Logistic模型又称为Logistic回归，主要在流行病学中应用较多，比较常用的情形是探索某疾病的危险因素，在预测人口的数量发展变化规律可以得到很好的效果。同样的经济预测领域也同样有广泛的应用。类比得到2011年至2035年的山东省职工的年平均工资的预测模型，根据附录一中已知数据运用Matlab软件画出散点图利用非线性最小二乘拟和的方法，求解得到参数即得到方程： （1）由此可得改革开放100年内山东省职工年平均工资的变化趋势图1并得到前57年（19782035）平均工资的变化趋势图2根据公式1预测得到2011年到2035年山东省职工年平均工资如表1所示表1 20112035山东省职工年平均工资年份2011201220132014201520162017工资/元37562.742754.148630.455272.762768.571212.480704.7年份2018201920202021202220232024工资/元91351.1103261.3116546.8131318.9147685.3165746.4185590.8年份2025202620272028202920302031工资/元207290.1230893.6256422.3283863.2313164.8344232.3376926.1年份2032203320342035工资/元411060.6446406.7482696.5519630.25.2问题二的模型建立与求解5.2.1计算2009年该企业各年龄段职工工资与该企业平均工资之比由于题目中未明确给出职工的具体工资，所以我们假设每个月收入范围内的工资符合均匀分布，因此取每个工资段的中间值作为每个年龄段中的个工资段的职工工资。根据附件二，求得各月收入范围的中间值分别为 1250，1750，2250，2750，3250，3750，4500，6500，令它们分别为，。根据Excel中的求和函数求得每个月收入范围的总的职工人数分别为110，306，392，410，267，115，52，12，令它们分别是，。所以该企业的平均工资， （2）用表示各个年龄段，即分别代表年龄段20-24，25-29，30-34，35-39，40-44，45-49，50-54，55-59。用表示工资为、年龄段为的企业人数，可以求得各年龄段职工平均工资与企业平均工资的比值 （3）利用附录二和公式（2）和（3），得到表2表2 职工平均工资和企业平均工资的关系年龄段/岁20-2425-2930-3435-3940-4445-4950-5455-59平均工资/元172720772535.12392.3302632683118.52980.2企业平均工资2580.22580.22580.22580.22580.22580.22580.22580.2职工平均工资/企业平均工资0.66930.80500.98250.92721.17281.26661.20861.15505.2.2计算各种情况下的养老金替代率 现在我国养老保险改革正处于过渡期。养老保险管理的一个重要的目标是养老保险基金的收支平衡，它关系到社会稳定和老龄化社会的顺利过渡。影响养老保险基金收支平衡的一个重要因素是替代率。替代率是指职工刚退休时的养老金占退休前工资的比例。用来表示替代率，表示企业和职工实际缴纳基本养老保险费的年限，表示参保人员i年度的本人缴费工资，表示各年利率，表示计发月数。假设把2009年该企业各年龄段职工工资与该企业平均工资的比值看作职工缴费指数的参考值，未来职工平均缴费指数同2009年相比产生的变化可忽略不计，用表示。那么可得公式其中利用Matlab软件求的6种情况的替代率如表3所示缴费年龄段/岁30-5530-6030-6540-5540-6040-650.12290.21870.37110.03970.06900.13465.3 问题三的模型建立和求解假设养老金随社会平均工资调整而调整的变化忽略不计（1）养老保险基金缴费积累（社会统筹个人缴纳）由题目已知2条件，带入数据统筹基金比例，个人账户比例，得到个人缴纳与社会统筹的比值其中用表示前i年社会统筹基金积累，表示前i年职工个人缴费工资额得到养老保险基金缴费积累 ：利用附录一和问题二的预测，由Matlab软件求得总共缴费积累为（）（2）领取的养老金（基础养老金个人账户养老金）由附录三，得到职工退休后领取的养老金为：资金缺口为时员工缴存的养老保险金与其领取的养老保险金之间达到收支平衡。综合附录一和附录二，以及上述（1）（2）分析，编写算法计算得到领取到75岁时的养老金的缺口情况和达到平衡的时间如表4所示：表4时间30-55岁30-60岁30-65岁缺口平衡时间5.4问题四的求解 6、模型的评价6.1 模型的优点（1）运用的模型简单易懂，有很好的实际指导意义。（2）运用表格和图像相结合，对于结果的分析更加清晰。（3）数学软件Matlab和Excel软件的运用提高了结果的可行度，数据更加精确。（4）对于题目中的问题做出了合理的假设，从而使计算变得更加容易。 7.2 模型的缺点（1）本题对数据依赖性比较大，只是根据题中所给数据做了一个理想化的模型可能与实际不相吻合。（2）我们建立的Logistic模型也只是一个预测模型，在实际生活中会有各方面的影响，因此通过我们所建立的模型只能大致体现一个发展趋势而无法精确地描述其发展趋势。参考文献王建.未来五年中国经济增长预测.人民论坛:5052,2012.02.百度百科 /view/2294104.htm 2011-9-12 附录程序1：建立M文件：function f=curvefun1(x,tdata)f=1./(x(1)+x(2)*exp(-x(3).*tdata)指令行：cleartdata=0:1:32;y=566,632,745,755,769,789,985,1110,1313,1428,1782,1920,2150,2292,2601,3149,4338,5145,5809,6241,6854,7656,8772,10007,11374,12567,14332,16614,19228,22844,26404,29688,32074;ydata=y/1000;plot(tdata,ydata,r*)x0=0.001,0.001,0.3;x=lsqcurvefit(curvefun1,x0,tdata,ydata) f=curvefun1(x,tdata) 程序2：x=1727 2077 2535.1 2392.3 3026 3268 3118.5 2980.2;y=2580.2;x/y程序3：y=0,8.772,10.007,11.374,12.567,14.332,16.614,19.228,22.844,26.404,29.688,32.074,37.5627,42.7541,48.6304,55.2727,62.7685,71.2124,80.7047,91.3511,103.2613,116.5468,131.3189,147.6853,165.7464,185.5908,207.2901;x=0.08*y;K=170;P=0.6693,0.8050,0.9825,0.9272,1.1728,1.2666,1.2086,1.1550;ave=0;sum=0;for j=1:8 ave=ave+0.125*P(j);endfor i=1:25 sum=sum+(x(i+1)-x(i).*1.03.i;endR=(1/K*sum+0.005*y(26)/(x(2)/0.08)程序4：y=0,8.772,10.007,11.374,12.567,14.332,16.614,19.228,22.844,26.404,29.688,32.074,37.5627,42.7541,48.6304,55.2727,62.7685,71.2124,80.7047,91.3511,103.2613,116.5468,131.3189,147.6853,165.7464,185.5908,207.2901,230.8936,256.4223,283.8632,313.1648,344.2323;x=0.08*y;K=139;P=0.6693,0.8050,0.9825,0.9272,1.1728,1.2666,1.2086,1.1550;ave=0;sum=0;for j=1:8 ave=ave+0.125*P(j);endfor i=1:30 sum=sum+(x(i+1)-x(i).*1.03.i;endR=(1/K*sum+0.005*y(31)/(x(2)/0.08)程序5：y=0,8.772,10.007,11.374,12.567,14.332,16.614,19.228,22.844,26.404,29.688,32.074,37.5627,42.7541,48.6304,55.2727,62.7685,71.2124,80.7047,91.3511,103.2613,116.5468,131.3189,147.6853,165.7464,185.5908,207.2901,230.8936,256.4223,283.8632,313.1648,344.2323,376.9261,411.0606,446.4067,482.6965,519.6302;x=0.08*y;K=101;P=0.6693,0.8050,0.9825,0.9272,1.1728,1.2666,1.2086,1.1550;ave=0;sum=0;for j=1:8 ave=ave+0.125*P(j);endfor i=1:35 sum=sum+(x(i+1)-x(i).*1.03.i;endR=(1/K*sum+0.005*y(36)/(x(2)/0.08)程序6：y=0,8.772,10.007,11.374,12.567,14.332,16.614,19.228,22.844,26.404,29.688,32.074,37.5627,42.7541,48.6304,55.2727,62.7685;x=0.08*y;K=170;P=0.6693,0.8050,0.9825,0.9272,1.1728,1.2666,1.2086,1.1550;ave=0;sum=0;for j=1:8 ave=ave+0.125*P(j);endfor i=1:15 sum=sum+(x(i+1)-x(i).*1.03.i;endR=(1/K*sum+0.005*y(16)/(x(2)/0.08)程序7：y=0,8.772