可持续的中国城乡居民养老保险体系的数学模型研究研究生数学建模F题

1、参赛密码 （由组委会填写）全第十届华为杯全国研究生数学建模竞赛学 校参赛队号队员姓名123- 3 -参赛密码 （由组委会填写） 第十届华为杯全国研究生数学建模竞赛题 目可持续的中国城乡居民养老保险体系的数学模型研究摘 要：在人口老龄化程度日益加深的背景下，我国的养老保险体系遇到了极大的挑战。养老金缺口已经出现并且缺口可能还会继续扩大。如果不及时解决养老金缺口问题，它将会给社会经济发展带来沉重的负担，如何保障养老保险体系持续健康稳定的发展，成为一个亟待研究的重要问题。为此，本文立足现实，探讨适合国情的养老保险体系，探寻促进我国城乡养老保险体系可持续发展的路径，同时为国家制定城乡居民养老保险相关政

2、策提供参考。 问题一 (建立中国城乡居民养老保险体系数学模型)，首先确定影响养老保险体系的相关因素，此处一方面利用题目给出的影响因素；另一方面借鉴相关专家学者的研究以及我国经济社会发展的现实情况，最终确定影响中国城乡居民养老保险体系的主要因素，包括养老金收入、养老金支出、养老金替代率、缴费率、养老金结余量、投资收益率、投资收益额、政府财政投入、gdp增长率、gdp增长额、分年龄段死亡率、养老金领取人数、cpi增长率、年平均工资水平、个人缴费总额、缴费率、缴费人数、劳动人口数、就业率。其次，确立各因素之间的假设关系。然后，构建基本养老保险体系的系统动力学模型，并进行系统仿真。模型建立要紧密结合社

3、会我国的经济发展状况。问题二，(1) 本文利用问题一中构建的基本养老保险系统动力学模型，设置基期为2010年的实际数据，预测今年至2035年我国养老金缺口，从2021年养老金开始出现缺口，并呈现扩大趋势。选取2000-2010年的数据来预测2011-2035年的养老金结余，其中2011年我国实际养老金结余为19496.6亿元，经过预测得到2011年的养老金结余为20672.1亿元，测量误差仅为6%，可见对养老金缺口的预测非常合理。(2) 如果国家现行养老保险政策保持不变，根据本文构建的数学模型，我国城乡居民养老保险收支矛盾最尖锐的情况将在2035年出现，基金缺口将高达2.00102e+011元

4、 。(3) 根据党的十八大提出的收入倍增计划，城乡居民人均收入比2010年翻一番，本文把模型中的工资提高一倍，则个人缴费总额将上涨，进而基金收入也会提高，在基金支出的增长不变情况下，养老保险基金结余量就会增多。问题三，尝试运用系统动力学模型进行仿真，调试替代率与缴费率，探寻最佳组合，经过不断调试，结合现实情况，最终确定替代率水平在40%-50%区间，缴费率在6%-8%之间，为合理区间。这不仅符合我国的实际情况，也与整个养老金系统发展的趋势较为匹配。 问题四，由于人口结构、死亡率、经济增速、投资收益等因素几乎无法人为较大幅度改动，另外，替代率与缴费率在第三问中已经讨论过。因此本文尝试将

5、增加政府财政投入作为可调节变量的数学模型，经过反复调试，发现增加政府财政投入有利于延缓养老金缺口的出现。本文创新点是利用系统动力学模型对养老保险体系进行仿真，综合多项因素，构成一个综合模型。有效说明我国养老保险体系的运行现状并且对未来养老金结余量进行科学合理地预测。关键词：城乡居民；养老保险体系；系统动力学；回归分析；1.问题的提出党的十八大报告提出了“统筹推进城乡社会保障体系建设”的任务：“社会保障是保障人民生活、调节社会分配的一项基本制度。要坚持全覆盖、保基本、多层次、可持续方针，以增强公平性、适应流动性、保证可持续性为重点，全面建成覆盖城乡居民的社会保障体系”。近年来，我国城镇职工养老保

6、险金制度基本建立，参保覆盖面日益扩大，退休人员待遇逐年提高，养老保险基金规模不断上升，在制度层面上提前完成了城市和农村居民社会保障全覆盖，成为我们社会保障事业发展的重要里程碑。然而我国的社会保障事业也面临着巨大的挑战。我国的养老保险体系建设远远滞后于社会经济发展水平。一方面资金规模相对较少，保障能力有限，地区差距、城乡差距明显；另一方面资金管理分散，运营效率不高，保值增值的压力较大。因此，进一步完善养老保险保障体系仍是我国政府未来最重要的任务之一。面临人口老龄化以及经济结构性减速、财政收入增长速度下降，社会保障体系可持续问题亟待解决，必须进行前瞻性研究，否则将形成社会问题。在研究的过程中，应该

7、承认历史、立足现实、尊重规律、借鉴国外、留有余地，既量力而行，又尽力而为地完善我国的养老保险体系。本题要求运用数学模型方法研究、解决该问题。为此，要做好以下几项工作。1分别建立合乎国情、适应国力的中国城乡居民(含新农保)养老金收入、支出的宏观数学模型，至少包括替代率(基本养老保险人均养老金占城镇单位在岗职工平均工资比率)、缴费率(基本养老保险人均缴费占城镇单位在岗职工平均工资比率)、人口结构、分年龄段死亡率、经济增速、财政补贴、工资水平或物价指数、投资效益等主要因素，要做到模型结合现实，分多个层次(含企业基金等)，体现“多缴多得，长缴多得”(不考虑分省、分地区模型)。2根据所建立的数学模型、对

8、养老金缺口的理解和对未来有关情况的合理估计，估计从今年至2035年我国养老金缺口，并说明对养老金缺口分析的合理性。如果全部情况维持不变，根据建立的数学模型我国城乡居民养老保险收支矛盾最尖锐的情况发生在什么时间，严重程度如何？考虑到党的十八大提出的收入倍增计划，建立的数学模型哪些部分需要调整？3养老保险制度也是调节社会分配，分析各国养老保险的不同模式，取其精华，去其糟粕，根据建立的数学模型和中国的实际情况，利用仿真手段寻找替代率和缴费率的合理区间以保证我国养老保险体系的可持续性(因为人口结构、分年龄段死亡率、经济增速、投资效益等主要因素几乎无法人为较大幅度改动)；在步入良性循环之前，在矛盾最尖锐

9、到来前的过渡期内应该采取哪些政策措施实现平稳过渡并仿真预测相关政策的效果。4. 尝试建立第三问增加可调节变量的数学模型。2.问题分析与建模思路本文首先借助文章所提示的影响养老保险体系的主要因素，结合我国的经济社会发展情况以及未来发展的趋势，构建系统动力学模型，进行系统仿真。然后运用该模型预测我国从今年至2035年养老金缺口。通过模型的模拟，探讨养老保险收支矛盾的程度。结合党的十八大提出的收入倍增计划，提出工资翻一倍，应当将此处在模型中体现。最后，分析各国养老保险的不同模式，总结各种模式的优缺点，学习他们成功的经验。以基期的数据为基本数据，在其它因素不变的情况下，利用仿真手段，寻找替代率和缴费率

10、的合理区间。同时，为了解决养老保险缺口问题，根据模型运行结果，提出政府的解决对策。3.基本假设与符号说明3.1模型假设1 假定预测得到的数据与实际中的数据相对误差在10%以内，就认为模型具有效性。2假定经济的发展态势是总体趋向良好，允许经济和平均工资水平出现小波动范围，但不存在大幅度的波动，就业市场能够实现充分自由流动。3 假定基期为2010年，未将2011年作为基期，是因为考虑到将用预测到的2011年的数据与2011年实际的数据进行对比，以检验实际预测数据误差的大小。这样更有利于预测2011年-2035年养老金账户的变化趋势。4 假定劳动年龄人口均参加养老保险 5 假设目前运行的养老保险制度

11、不发生变化3.2符号说明表3-1 符号说明符号意义gdp国民经济发展水平cpi消费物价指数4.第一题求解4.1 系统动力学的基本原理4.1.1系统动力学简介系统动力学是系统科学理论与计算机仿真紧密结合、研究系统反馈结构与行为的一门学科，是系统科学与管理科学的一个重要分支。该理论以系统作为研究对象，初期主要应用于工业企业管理，处理生产和雇员情况的波动、市场股票与市场增长不稳定性等问题1。系统动力学的主要原理是系统的行为模式与特性主要取决于其内部的结构，系统内部各个因素相互影响最终形成反馈链，各反馈链最终影响整个系统的运行。系统动力学研究和处理复杂问题的办法是定性和定量分析相结合、系统综合推理的方

12、法，其建模过程就是学习、调查、研究的过程。模型的功能主要在于为人们提供一个学习和分析的工具。应用系统动力学建模需要充分考虑模型与现实情况的相似性，对现实数据复制的准确度、精准度；对于非线性、高阶次的复杂系统，系统动力学相对于其他研究方法更有优势。4.1.2 建模步骤系统动力学的建模步骤主要包括以下环节：第一步，明确建模的目的，选择需要建模的系统，确定建模时限、系统边界，选择关键变量，这是利用系统动力学建模不可缺少的前期准备工作；第二步，提出动态假设，根据现有的理论或者当前系统的运行情况，提出一个由系统内部反馈结构导致动态变化的假设，根据初始假设、关键变量、参考模式和其他可用数据绘制出因果回路图

13、、存量流量图；第三步，写入公式，确定参数、行为关系和初始化条件。根据现有数据资料进行整理分析，对系统中涉及的变量进行预测，用系统动力学语言写入存量流量图中；第四步，对于得出的预测结果进行测试，检验模型得出的行为结果是否符合现实，对不符合现实部分结论的产生条件进行修正，反复调试，对模型中各个参数、初始化条件、模型边界和概括程度进行灵敏度测试；第五步，对于模型产生结果进行分析，得出具体化的方案设计。4.2 模型构建4.2.1建模思想和数据来源(1) 建模思想通过系统动力学软件仿真养老保险基金收支的过程，测算养老保险基金收支以及影响因素，力求找到针对未来养老保险基金运营可能存在的“空账”问题，找到现

14、实解决办法，以保证我国养老保险制度稳定，为我国离退休群体提供更好的制度保护和防范措施。(2) 数据来源历年的中国统计年鉴、中国人力资源和社会保障发展公报、世界卫生组织网站、中国劳动统计年鉴、国家统计局和社会保障基金理事会等网站。4.2.2 参数设定及参数方程(1) 参数设定根据养老保险基金系统理论，考虑到基金流入和流出的动态性，将各个参数因子的选取和定性如下表格所示，其中结余量、年平均工资水平和人均gdp为存量，即状态变量，基金收入、支出、年平均工资增长额和gdp增长额为流量，其他的参数设定为辅助变量。表4-1 参数选取变量名性质变量名性质结余量存量领取人数辅助变量年平均工资水平存量死亡率辅助

15、变量人均gdp存量个人领取养老金辅助变量基金支出流量缴费率辅助变量基金收入流量投资收益辅助变量工资年平均增长额流量投资收益率辅助变量gdp增长额流量缴费人数辅助变量政府财政投入辅助变量劳动人口数辅助变量个人缴费总额辅助变量就业率辅助变量劳动人口总数辅助变量人均gdp增长率辅助变量工资增长率辅助变量替代率辅助变量cpi增长率辅助变量领取人数总表辅助变量死亡率总表辅助变量(2) 参数方程根据参数因子的选取与各因子之间的关系，设定参数方程如下：年平均工资水平=年平均工资水平+年平均工资增长额年平均增长额=年平均工资水平*工资年平均增长率工资年平均增长率=1.2*cpi增长率+人均gdp增长率结余量=

16、基金收入-基金支出基金收入=个人缴费总额+政府财政投入+投资收益额基金支出=领取人数*(1-死亡率)*个人领取养老金投资收益额=结余量*投资收益率政府财政投入=-4456.35+0.2101*人均gdp个人缴费总额=年平均工资水平*缴费率*缴费人数缴费人数=劳动人口数*就业率个人领取养老金=年平均工资水平*替代率人均gdp=人均gdp+gdp增长额gdp增长额=人均gdp*人均gdp增长率劳动人口数=劳动人口总数（time）= (2010, 8.891e+008)( 2010,9.9938e+008), (2000,8.891e+008),(2001,8.9849e+008),(2002,9.

17、0302e+008),(2003,9.0976e+008),(2004,9.2184e+008),(2005,9.4197e+008),(2006,9.5068e+008),(2007,9.5833e+008),(2008,9.668e+008),(2009,9.7484e+008),(2010,9.9938e+008)领取人数=领取人数总表 (time)= (2000,0)-(2100,8e+008),(2010,1.18e+008),(2011,1.22e+008),(2012,1.28e+008),(2013,1.35e+008),(2014,1.43e+008),(2015,1.52e

18、+008),(2016,1.62e+008),(2017,1.74e+008),(2018,1.88e+008),(2019,2.03e+008),(2020,2.2e+008),(2021,2.39e+008),(2022,2.6e+008),(2023,2.83e+008),(2024,3.09e+008),(2025,3.37e+008),(2026,3.67e+008),(2027,4e+008),(2028,4.36e+008),(2029,4.75e+008),(2030,5.17e+008),(2031,5.62e+008),(2032,6.11e+008),(2033,6.62

19、e+008),(2034,7.18e+008),(2035,7.77e+008)死亡率=死亡率总表 (time)= (2000,0.00645), ( 2010, 0.00711), (2000,0.00645),(2001,0.00643),(2002,0.00641),(2003,0.0064),(2004,0.00642),(2005,0.00651),(2006,0.00681),(2007,0.00693),(2008,0.00706),(2009,0.00708),(2010,0.00711)4.2.3 主要参数因子解释(1) 年平均工资增长率年平均工资增长率=（当年的工资-上一年

20、的工资）/上一年的工资，考虑到一个地区的工资水平主要受到当地的经济增长率和通货膨胀率的影响，而通货膨胀率与居民消费价格水平之间有一定的影响，通货膨胀率= (报告期消费者物价指数 - 基期消费者物价指数)/基期消费者物价指数* 100% ，所以在本文中，将影响年平均工资增长率的因素界定为通货膨胀率和gdp增长率。通过折线图，发现gdp增长率和cpi增长率大致上呈现正增长关系，但是增长比例不是特别大，通过对cpi增长率倍数进行赋值1.1、1.2、1.3发现当取1.2时，总体上来看比较接近，为了进一步验证所取值的正确性，对其进行逐步回归分析，对两者进行线性拟合，结果显示所取值有可信性，因此可以假定年

21、平均增长率=gdp增长率+1.2*cpi增长率。图4-1 2001-2010年工资水平-gdp-cpi增长率折线示意图本文的模型假设gdp增长率选用2010年基期的数据，为9.9%， 同样，cpi增长率也选取基期数据，cpi增长率为3.2%。(2) 缴费人数考虑到养老保险参保人数的准确性，将劳动人口（15-64周岁）和就业率纳入到指标构建体系中，假定该年龄段的就业人口均已参加养老保险，即缴费人数等于劳动人口与就业率之间的乘积。个人缴费总额即为个人的缴费比例与缴费人数的乘积。根据国家统计年鉴规定参保人数为报告期末按照国家法律、法规和有关政策规定参加基本养老保险并在社保经办机构已建立缴费记录档案的

22、职工人数，包括中断缴费但未终止养老保险关系的职工人数，不包括只登记未建立缴费记录档案的人数，因此得到2000年-2011年我国的劳动人口数如下表4-2所示。表4-2 2000 -2011年15-64岁劳动人口（单位：万人）年份人口200088,910200189,849200290,302200390,976200492,184200594,197200695,068200795,833200896,680200997,484201099,938（数据来源：中国统计年鉴2011）基于指标构建的可信性，能够更好地反映模型中参保人数与实际参保人数 的一致性，本文采用城镇登记失业率变量来表现就业率的

23、变化趋势。从下表中可以看出失业率在2000年以后都有较大的提升，其中在2005年达到最大值，之后呈现下降趋势，但是总体数值还是较高。为了能够更好地反映该变量的变化趋势，对1992年至2010年的城镇登记失业率取平均值，得到平均登记失业率4.22%，根据国家统计年鉴的相关规定，就业率指就业人口与16岁以上总人口的百分比。凡在制定时期内届满一定下限年龄，有工作并取得报酬或收益的人，或有职位而暂时没有工作（如生病、工商、劳资纠纷、假期等）的人以及家庭企业或农场的无报酬工作者，均算为就业人口。因此本文将就业率假定为95.78%。图4-2 城镇登记失业率（数据来源：根据历年中国统计年鉴整理）(3) 政府

24、财政投入养老保险基金的财政投入能力与国民经济增长水平是密不可分的，因此选用财政投入作为对gdp的函数，分析政府财政对养老保险的投入与gdp之间的关系，以找到两者之间的关系。表4-3 我国2000年-2010年财政支出与gdp数据年份财政支出（亿元）gdp(亿元)200015886.599214.6200118902.58109655.2200222053.15120332.7200324649.95135822.8200428486.89159878.3200533930.28184937.4200640422.73216314.4200749781.35265810.3200862592.6

25、6314045.4200976299.93340902.8201089874.16401512.8(数据来源：中国统计年鉴2010)运用eviews软件对我国2000年-2010年财政支出与gdp进行回归 ，得到以下方程式：caizheng= -4.076290 + 1.197968*gdp -13.16241 47.38784 r2= 0.995567 f= 2245.608 dw= 1.106238从估计检验结果可以看出，回归系数的标准误差非常小，t统计量较大，决定系数 r2= 0.995567 非常高，f= 2245.608，说明gdp对财政支出影响的回归模型非常显著，同时dw= 1.1

26、06238 > r2= 0.995567，说明总体回归模型的随机误差ui之间不存在相关关系，得出的结果是合理的。(4) 投资收益额1997年，国务院关于建立统一的企业职工基本养老保险的决定提出：“基本养老保险基金实行收支两条线管理，要保证专款专用，全部用于职工养老保险，严禁挤占挪用和挥霍浪费基金结余额，除预留相当于2个月的支付费外，全部购买国家债券和存入银行专户，严禁投入其他金融和经营性事业。”可见我国养老保险基金投资实行较为严格的控制，渠道单一，主要是依靠国家债券和银行存款。下表为我国1999-2010年银行存款一年期的年利率，每一年的年利率是取的当年内几次利率调整的平均值。 表4-4

27、 1999-2010年我国银行存款一年期年利率(%)年份199920002001200220032004利率3.022.252.252.21.982.1年份200520062007200820092010利率2.252.352.333.542.252.25计算得出1999-2010年的平均一年期银行存款年利率为2.4%，短期国债收益率也在2%以上。下表为我国养老保险基金历年投资收益率，计算得出10年来平均投资收益率为10.04%，但是从整个投资收益率发展趋势来看，2006年和2007年投资收益率有一个激增的现象，这与2006国家实行养老基金和社会保险基金投资体系改革有一定关系。这从整体上拉高了

28、投资收益率。因此本文不选用该平均投资收益率。图4-3 2001-2010年投资收益率据新闻媒体报道 ，养老金投资体系改革方案将有望于今年内出台，社会普遍关心的养老金整体投资运营将取得实质性突破。专家表示，未来改革思路面临两种选择。一是延续广东模式，委托全国社保基金理事会投资运营；二是在中央一级新设投资运营机构统一投资管理。这说明我国未来养老保险基金投资运营科技能会有一个不断完善改进的过程，此处，我们根据国家一年期投资收益率和短期国债收益率，将养老基金投资收益率设定为2.5%。那么投资收益额就等于投资收益率与结余量之间的乘积。(5) 养老金收入与支出养老金收入主要是由养老保险缴费、政府转移支付和

29、养老基金投资收益三部分组成。养老保险缴费由于此处是城乡居民养老保险，因此不考虑单位缴费，只考虑个人缴费部分养老保险基金收入=政府财政投入+个人缴费总额+投资收益额养老金支出指按照国家政策规定的开支范围和开支标准从养老保险基金中支付给参加基本养老保险的个人的养老金、丧葬抚恤补助，以及由于保险关系转移、上下级之间调剂资金等原因而发生的支出。养老金支出=领取人数*（1-分年龄段死亡率）*个人领取养老金(6) 人均gdp国内生产总值指按市场价格计算的一个国家（或地区）所有常住单位在一定时期内生产活动的最终成果 国家统计局网站。人均gdp常作为发展经济学中衡量经济发展状况的指标，是重要的宏观经济指标之一

30、，它是人们了解和把握一个国家或地区的宏观经济运行状况的有效工具。将一个国家核算期内（通常是一年）实现的国内生产总值与这个国家的常住人口（目前使用户籍人口）相比进行计算，得到人均国内生产总值。(7) 养老金支出养老金支出是指按照国家政策规定的开支范围和开支标准从养老保险基金中支付给参加基本养老保险的个人的养老金、丧葬抚恤补助，以及由于保险关系转移、上下级之间调剂资金等原因而发生的支出2。此外，本文为计算方便，只考虑退休人数与养老保险待遇水平的乘积。养老金支出=领取人数*（1-死亡率）*个人领取养老金数(8) 领取人数领取人数是考虑达到退休年龄，适龄领取养老金的人口数，实事上就是参保退休人口数，但

31、是在我国目前不是所有的达到退休年龄的人口都能享受养老金待遇，但是国家建立完善社会保障制度，就是为了使所有的达到退休年龄的人口都能享受养老金，假设所有达到退休年龄的人都可以享受养老保险待遇，因此此处将所有退休年龄人口纳入其中。目前我国退休年龄男女不同，男性退休年龄60岁，女职工50岁退休，女干部55岁退休，在延迟退休年龄的相关研究日间成熟，呼吁延迟退休年龄的呼声越来越高的背景之下，本文将退休年龄统一设在65周岁，所有达到65周岁以上的人口数即为养老金领取人数。此处根据2000至2010年10年间的65岁以上人口数，运用spss软件，选取arima模型进行人口预测，预测至2035年的65岁以上人口

32、数。经过调试，选取最佳结果arima(2，3，3)。预测值与实际值之间的相对误差值为0.79和0.72，在可接受的范围之内。表4-5 2011-2035年65岁以上人口预测model201020112012201320142015201620172018人口模型forecast1.18e+081.22e+081.28e+081.35e+081.43e+081.52e+081.62e+081.74e+081.88e+08ucl1.31e+081.37e+081.45e+081.61e+081.78e+081.95e+082.18e+082.44e+082.73e+08lcl1.04e+081.0

33、8e+081.12e+081.09e+081.08e+081.09e+081.07e+081.05e+081.03e+08model201920202021202220232024202520262027人口模型forecast2.03e+082.20e+082.39e+082.60e+082.83e+083.09e+083.37e+083.67e+084.00e+08ucl3.04e+083.39e+083.78e+084.21e+084.69e+085.21e+085.77e+086.38e+087.04e+08lcl1.02e+081.01e+081.00e+08990308819740

34、082496777675967844229650683796629117model20282029203020312032203320342035人口模型forecast4.36e+084.75e+085.17e+085.62e+086.11e+086.62e+087.18e+087.77e+08ucl7.75e+088.50e+089.31e+081.02e+091.11e+091.21e+091.31e+091.42e+09lcl981028771.01e+081.03e+081.07e+081.11e+081.17e+081.24e+081.32e+08人口死亡率是用来衡量一部分人口中、

35、一定规模的人口大小、每单位时间的死亡数目（整体或归因于指定因素）。死亡率通常以每年每一千人为单位来表示。由于65岁以上人口死亡率会影响到领取养老金的人数，因此本文搜集整理2000-2010年10年的65岁以上人口死亡率，如表所示：表4-6 65岁以上人口死亡率（）年份20002001200220032004200520062007200820092010人口死亡率6.456.436.416.46.426.516.816.937.067.087.11运用spss软件进行人口死亡率的预测，经过反复调试，选取最佳结arima(1.0.2)，对2011至2035年的65岁以上人口死亡率进行预测，201

36、1年实际人口死亡率为7.14，预测结果2011年为7.15，误差为0.14%，效果较为理想，预测结果如下表。表4-7 65岁以上人口死亡率预测（）年份65岁以上人口死亡率年份65岁以上人口死亡率20117.15120248.07620127.20920258.14820137.28520268.2220147.35820278.29120157.4320288.36320167.50220298.43520177.57420308.50720187.65%20318.57820197.71820328.6520207.78920338.72220217.86120348.79320227.93

37、320358.86520238.005(9) 养老金替代率养老金替代率是退休人员的养老金待遇水平与在岗职工平均工资的比例，是衡量劳动者退休前后生活保障水平的基本指标之一。柳清瑞(2003)利用精算模型测算我国养老金替代率在60%左右是合适的，可以使退休人口的养老金待遇与中等偏上户的消费额大体相当3。郝勇 (2012)利用结构方程和系统动力学模型对我国养老金替代率进行测算，认为养老金替代率适度水平应该稳定在39%-56%附近4。因此，本文设定养老金替代率为60%，这也符合国际上关于养老金替代率水平的规定。个人领取的养老金数额与退休前年平均工资水平以及养老金替代率有关。个人领取养老金数额=年平均工

38、资水平*替代率4.2.4 因果关系图及流程图通过对各参数之间关系的梳理，画出因果关系图如下所示：图4-4 因果关系图因为因果关系图只能反映各个参数之间的大致关系，不能具体进行操作，这就需要流程图对相关的参数变量进行进一步的表达，因此绘制结构流程图如下所示：图4-5 结构流程图4.2.5 模型有效性检验对模型的有效性检验主要包括：结构确认测试、参数确认测试、单位一致性测试、极端情况测试和行为重现测试。(1) 结构确认测试直接比较所建立的养老保险基金模型符合真实中的系统结构。通过反馈关系图、因果关系图进行检验，发现本模型与真实的失业保险基金政策结构是一致的。a 反馈关系图通过对建立的模型流程图进行

39、分析，得到以下的反馈关系图，通过观察认为该反馈关系图具有科学的合理性，因此认为建立的养老保险基金模型满足反馈关系的检验，具有有效性。b 因果树检验通过模型的因果树检验，得到以下因果树关系图，通过比较分析，认为该因果树关系图具有科学的合理性，因此认为建立的养老保险基金模型满足因果树检验，具有有效性。(2) 参数确认测试直接比较所构建的养老保险基金模型与模型所代表的真实知识是否在概念和数字上符合。本文所用的数据参见历年的国家统计年鉴、社会保障基金理事会网站、中国人力资源与社会保障事业发展公报以及通过所学的统计知识进行数据处理与加工得到，在概念和数字上符合养老保险基金所代表的真实的系统相关的参考文献

40、。(3) 单位一致性测试确认所建立的模型各变量的单位是否一致。通过vensim软件的units check测试所建的模型，得到建立的模型已通过单位一致性测试，显示model is ok。(4) 极端情况测试测试所建立的模型在极端情况下是否仍符合真实情况。当极端情况测试的结果不符合真实情况时，则开始修改系统，再次测试，当极端情况测试出现的结果与真实情况相符合时，停止测试。因为在初始值设置的时候，本文对增长率和比例都进行了设定，所以在调整参数比例变动时，就可以避免发生极端情况。(5) 行为重现测试测试所建构的模型是否会重新出现现实中发生的问题行为。在本文中对参保人数和基金结余量进行了测试，模拟结果

41、与统计数据在变化趋势上是一致的，结果比较有信度。经过检验，表明该模型可用，可以进行数学建模。5第二题求解5.1 预算至2035年我国养老金缺口利用第一问对养老保险系统的仿真，以2000-2010年数据为源数据，设置2010年为基期，对2011年到2035年的基金结余量未来进行预测，设定2010年的养老金替代率为50%，gdp增长率为9.9%，cpi增长率为3.2%，投资收益率为4.23%，就业率为95.78%。将这些数据录入模型之中，运行得到结果如下：表5-1 2010-2035年养老保险基金结余量预测time (年)结余量time (年)结余量20101.34e+122023-1.31e+1

42、320112.07e+122024-2.43e+1320122.84e+122025-4.02e+1320133.64e+122026-6.21e+1320144.42e+122027-9.18e+1320155.12e+122028-1.32e+1420165.66e+122029-1.84e+1420175.94e+122030-2.54e+1420185.76e+122031-3.45e+1420194.89e+122032-4.62e+1420203.04e+122033-6.13e+142021-2.00e+112034-8.07e+142022-5.35e+122035-1.05e

43、+15图5-1 2010-2035年养老保险基金预测结余量图根据预测结果，养老保险基金在2018年出现账户亏空，随着人口老龄化程度的加大和青壮年人口的减少，养老保险基金的亏空额度继续加大，在2035年，养老保险收支矛盾最尖锐，达到2.00088e+012元，难以保障退休人员的生活，并影响到养老保险体系的运行。5.2 收入倍增计划根据十八大提出的收入倍增计划，将原始数据年平均工资水平在基期水平之上增加一倍，也就是年平均工资为74294元，则个人缴费总额将上涨，进而基金收入提高，在基金支出的增长不变情况下，养老保险基金结余量增多。考虑到工资的增长幅度与国家的经济发展水平是密不可分的，因此首先要将g

44、dp增长率纳入到收入倍增模型中，虽然cpi增长率会对gdp增长率有一定的影响，但是不会直接构成对年平均工资水平影响。虽然gdp增速减缓，但并不代表的经济发展水平下降，也不意味着政府会减少对社会保障的财政投入，随着我国经济结构的战略性调整和转型升级，必将会进一步促进经济增长，增加政府的社会保障财政收入，减少养老金缺口，以保障人民群众的切身利益。如下表为收入倍增计划下的基金结余量预测。表5-2 收入倍增计划下基金结余量预测 time (年)结余量time (年)结余量20101.34e+122023-2.85e+1320112.74e+122024-5.10e+1320124.23e+122025

45、-8.28e+1320135.76e+122026-1.27e+1420147.25e+122027-1.86e+1420158.58e+122028-2.66e+1420169.60e+122029-3.72e+1420171.01e+132030-5.11e+1420189.66e+122031-6.93e+1420197.82e+122032-9.28e+1420204.04e+122033-1.23e+152021-2.52e+122034-1.62e+152022-1.29e+132035-2.11e+155-2 收入倍增计划下基金结余量预测 从上表中可以看出，养老保险基金在2021

46、年出现账户亏空，收不抵支，亏空额度达到-2.52e+12元，在收入倍增计划下，随着年平均工资水平的上升，个人缴费额度增大，基金的收入量加大，在一定程度上起到缓冲作用，能够较好地防止缺口的扩大，预防养老金支出带来的支付危机6. 第三题求解6.1 各国养老保险模式分析目前世界各国实行的养老保险制度，归纳起来有三种模式，分别为传统型养老保险制度、国家统筹型养老保险制度和强制储蓄型养老保险制度。传统型养老保险制度为德国俾斯麦政府在1889年颁布养老保险法所创设，后来逐渐被美国和日本所采纳，这种模式对不同的社会成员选用不同的保险标准，并以劳动者为核心建立养老保险制度，强调劳动者个人在社会保险方面的责任，

47、社会保险费用由国家、雇主和劳动者三方负担，以劳动者和雇主的缴费为主，国家财政给予适当支持。劳动者享受的养老保险待遇水平与个人缴费多少和个人收入情况相联系。该种模式重视社会保险中的权利和义务的密切联系，强化社会成员的自我保障意识，同时，社会保险基金可以在参保成员间统筹使用，符合保险互助互济的标准，但采取现收现付的方式筹集基金，费率受人口年龄结构与人口就业比例影响较大，在人口老龄化和就业比例下降时，费率过高增加了企业和个人的负担5。国家统筹型的养老保险制度最早为英国所创设，目前以英国、瑞典、挪威和澳大利亚等国家为代表。该制度实行完全的现收现付制，并以“支付确定”的方式来确定养老金水平，该制度的优点

48、在于运作简单运行，通过收入再分配的方式，为老年人提供基本的生活保障，但该制度的缺点是政府财政负担过重，缺乏对个人的激励机制，注重公平而忽视效率。强制储蓄型养老保险制度是一种主要依靠个人缴费的养老保险制度，以新加坡的中央公积金作为代表，该模式主要强调自我保障，由劳动者在职期间与其雇主共同缴纳养老保险费，劳动者退休后从个人账户领取养老金，个人账户的基金在劳动者退休后可以一次性领取，也可以分批分期领取，政府通过中央公积金局对个人账户的基金进行统一管理和运营投资，该种模式的优点是强调效率，调动了劳动者自我保障的积极性，缺点是忽视了社会公平，难以体现社会保险的保障功能。各国养老保险的不同模式各有自己的优

49、点和缺点，我国建立的养老保险制度借鉴了国外模式的成功经验，同时从本国的国情出发建立社会统筹和个人账户相结合的养老保险制度，该制度即吸收了传统型养老保险制度的优点，也借鉴了个人账户的长处，不仅强调了社会成员的自我保障意识，也体现了社会保险互助互济、分散风险的特点。随着我国人口结构的逐渐变化，劳动力供应总量的不断减少，我国人口红利的优势逐渐消失，随之而来的是是沉重的养老压力，由于在现行养老保险制度确立之前，企业员工基本无需缴纳养老保险费用，所以现有的养老保险基金中没有这部分职工的个人账户部分。在建立养老保险制度后，这部分职工退休后却要从养老保险基金中领取养老金，这是我国社会保险制度必须承担的转轨成

50、本的一部分。仅靠统筹账户难以应对当期发放，同时各地财政实力不同，多数地区不得不在实际上采用了“现收现付制”的方法，即挪用个人账户的资金，用正在工作人合计缴纳的28%的月工资来支付现有退休人员的退休金，现在工作人的账户仅仅用来记账，上述2.5万亿的记账额由此形成。为解决这一问题，中国从2000年开始了“做实”个人账户试点。截至到2011年底，参与试点的辽宁、江苏和山东等13个省份共积累个人账户基金达到2703亿元，但其与记账额之间的差额，仍达到2.23万亿元，此“空账”被舆论定义成了养老保险“缺口”概念而广为传播。根据已有的数学模型，对缴费率和缴费比例进行不断调试，假定一变量固定不变，另一变量在

51、某一个范围内按照某一定的规定进行变动，依次调试，得到最理想的缴费率和替代率。6.2 替代率和缴费率合理区间探讨根据中国社科院发布2013年经济蓝皮书指出，稳增长是当前经济工作中的首要任务。据中国社科院世界社保研究中心发布的中国养老金发展报告2012对城镇基本养老保险替代率进行测算，数据显示，养老金替代率由2002年的72.9%下降到2005年的57.7%，此后一直呈下降趋势，到2011年，这一数字更是降至50.3%。中国社科院财贸所学者高培勇和汪德华最近的一项研究指出，在现有社保制度框架下，要确保2020年中国每个退休者都能领取养老金，以替代率(即领取养老金占该年社会平均工资的比重)52.4%

52、。而现阶段企业退休职工的养老金替代率为到42.9%，为了能够更好地涵盖所有可能发生突变的变量，本文将替代率分为40%、50%、60%和62%。缴费率设定为2%、3%、4%、5%、6%、7%、8%。依次对其进行分析，得到合理的区间范围。根据以上基本建模思路，进行构建基于养老金替代率和缴费率的养老保险基金模型。当缴费比例为2%和4%时，预测得到的基金结余量在2011年和2012年出现负值，当缴费比率为5%，6%时，基金结余量可能在2012和2013年出现负值。这与统计年鉴中的数值不相匹配，因此将这四个比例剔除。同样的道理，剔除替代率的60%和62%。如下表为替代率和缴费率的模拟组合表。表6-1 替

53、代率和缴费率的模拟组合表替代率缴费率替代率缴费率40%5%50%5%6%6%7%7%8%8%替代率缴费率替代率缴费率60%5%62%5%6%6%7%7%8%8%当替代率50%，缴费率7% 时，养老金结余发展趋势如下所示。表6-2 替代率50%，缴费率7%组合下养老金结余情况表time (年)结余量time (年)结余量20101.34e+122023-5.74e+1320112.02e+122024-8.54e+1320122.67e+122025-1.24e+1420133.21e+122026-1.75e+1420143.52e+122027-2.43e+1420153.46e+122028-3.33e+1420162.85e+122029-4.50e+1420171.42e+122030-6.02e+142018-1.24e+122031-7.99e+142019-5.68e+122032-1.05e+152020-1.25e+132033-1.37e+152021-2.26e+132034-1.78e+152022-3.71e+132035-2.30e+15在2018年将出现养老金缺口，届时养老金缺口为-1.24e