高社杯大学生数学建模竞赛C题企业退休职工养老金制度的改革论

1、2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： C 我们的参赛报名号为（如果赛区设置

2、报名号的话）： 162C02 所属学校（请填写完整的全名）： 浙江树人大学 参赛队员 (打印并签名) ：1. 王 瑶 2. 林海蔚 3. 指导教师或指导教师组负责人 (打印并签名)： 王芬 日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：企业退休职工养老金制度的改革摘要为了适应社会主义市场经济体制改革的需要，同时也是为了应对老

3、龄化社会的风险，我国企业职工基本养老保险实行“社会统筹”与“个人账户”相结合的模式，即“统账结合”式基本养老金制度，但是这个制度还是存在不少问题，如人口老龄化问题、隐性债务问题、个人账户空帐问题等都没用提出有效解决方案，基本养老保险的基金收支平衡仍是个难度很大的问题。本文通过反正切模型对职工未来年平均工资进行预测，并在此基础上对企业退休职工养老金的相关数据进行计算。对于问题一，首先证明了职工年平均工资与人均GDP值成线性关系，相关系数为0.9985。所以经济的发展趋势与工资的增长趋势一致。而经济的发展将会从快速增长向缓慢增长转变，所以工资增长也会呈现这样的规律。为此，采用符合S型增长特征的反正

4、切函数作为预测模型，利用山东省职工年平均工资的历史数据，采用最小二乘方法，拟合出模型的参数。通过这个模型预测从2011年到2036年的年平均工资。预计到2035年职工年平均工资为80063元。对于问题二，我们利用附件中的数据，加上问题一的预测数据计算出30岁、40岁开始缴养老保险到退休时的养老金替代率，结果如下表所示年龄段（岁）30-5530-6030-6540-5540-6040-65养老金（元）1461.412031.342988.51729.641265.832023.03替代率55%76%112%27%48%76%对于问题三，综合利用问题一、问题二的结果，发现该职工自2000年起从30

5、岁开始缴养老保险，一直缴费到退休（55岁，60岁，65岁），并从退休后一直领取养老金到65岁、70岁、74岁时，其缴存的养老保险基金与其领取的养老金之间达到收支平衡，并且出现缺口。由于该职工是75岁死亡，假设他在75岁那年还能领到养老金，那么缺口就在75岁的时候最大，分别为-190965.54元、-133786.04元、-29789.18元。对于问题四，既要达到目标替代率，又要维持养老保险基金收支平衡，我们提出了以下建议措施：1、延长退休年龄2、延长养老金缴费年限。3、改革现有的按工资指数化的管理方法，以相应地减轻养老金支付负担和转制成本，从而缓解养老基金收支不平衡的状况。4、建立专门部门管理

6、职工上缴个人账户养老金，为基金做投资增加基金收益率。关键词：S型增长模式；反正切模型；最小二乘拟合；替代率；收支平衡一、问题的重述企业退休职工养老金过低的问题现已成为社会关注的焦点，现在我国养老保险改革正处于过渡期，企业职工养老保险制度改革要处理好国家利益、集体利益和个人利益，目前利益和长远利益，整体利益和局部利益的关系，主要是对现行的制度办法进行调整、完善。我国企业职工基本养老保险实行“社会统筹”与“个人账户”相结合的模式，“社会统筹”基金它采取的是“现收现付”的筹资方式，即从现有在职人员那里筹集资金，支付已退休职工的养老金，从而形成“下一代供养上一代”的社会契约；“个人账户”基金它选择“储

7、备积累”筹资方式，通过个人账户的积累为将来退休的职工预先筹集养老保险基金。而养老金会随着社会平均工资的调整而调整。如果职工死亡，社会统筹账户中的资金不退给职工，个人账户中的余额可继承。个人账户储存额以银行当时公布的一年期存款利率计息，为简单起见，利率统一设定为3%.养老金的发放与职工在职时的工资及社会平均工资有着密切关系；工资的增长又与经济增长相关。养老保险管理的一个重要的目标是养老保险基金的收支平衡，它关系到社会稳定和老龄化社会的顺利过渡。影响养老保险基金收支平衡的一个重要因素是替代率。照国家对基本养老保险制度的总体思路，未来基本养老保险的目标替代率确定为58.5%.替代率与退休职工的生活水

8、平恰好成正比养老保，保险基金收支平衡较难维持，可能出现缺口。为了促进企业退休职工养老金制度的改革，下列问题是值得深入研究的。（1）对未来中国经济发展和工资增长的形势做出你认为是简化、合理的假设，并参考附件1，预测从2011年至2035年的山东省职工的年平均工资。（2）根据附件2计算2009年该企业各年龄段职工工资与该企业平均工资之比。如果把这些比值看作职工缴费指数的参考值，考虑该企业职工自2000年起分别从30岁、40岁开始缴养老保险，一直缴费到退休（55岁，60岁，65岁），计算各种情况下的养老金替代率。（3）假设该企业某职工自 2000年起从30岁开始缴养老保险，一直缴费到退休（55岁，6

9、0岁，65岁），并从退休后一直领取养老金，至75岁死亡。计算养老保险基金的缺口情况，并计算该职工领取养老金到多少岁时，其缴存的养老保险基金与其领取的养老金之间达到收支平衡。（4）如果既要达到目标替代率，又要维持养老保险基金收支平衡，你认为可以采取什么措施。请给出你的理由。二、基本假设1、假设国外社会稳定，政府宏观政策方面没有发生大的失误，导致金融风险、财政风险积累，甚至爆发危机等无其它特殊因素。2、假设参保人员从开始缴纳养老保险时每月按时、按额缴纳，中间无间断和拖欠现象。3、假设缴纳到社会统筹基金账户的资金每年以3%的存款利息利率增长。4、假设养老金缴纳完毕后从社会统筹资金中拨出的基础养老金的

10、余额也以每年3%的存款利息利率增长。5、假设不考虑通货膨胀因素。6、假设基金收入不考虑政府财政补贴，基金支出不考虑用于管理方面的费用支出及职工死亡丧葬费用等，基金收入完全用于参保退休职工养老金的发放。三、符号说明为了便于描述问题，我们用一些符号来代替问题中涉及的一些基本变量，如表1所示。其他一些变量将在文中陆续说明。表1 主要变量符号说明一览表符号意义本人指数化月平均缴费工资（Average Indexed Monthly）参保人员实际缴费月数和参保人员退休前1年、2年、m年本人缴费工资额参保人员退休前1年、2年、m年全国/省/地市“职工平均工资”或称“社会平均工资”企业和职工实际缴纳基本养老

11、保险费的月数合计退休前第i年的缴费指数，i =1,m.工资指数参保人员i年度的本人缴费工资通过工资指数得到指数化缴费工资表示第年份在职职工的月平均缴费工资，分别表示社会统筹账户和缴费比例与个人账户的缴费比例四、问题分析4.1 问题一的分析题目中定性地告诉我们养老金的发放与职工在职时的工资及社会平均工资有着密切关系，工资的增长又与经济增长相关。为了更好地说明经济发展与工资增长之间的相关性，我们可以收集我国人均GDP的历史值（如图4.1.1所示），去计算出GDP值和职工年平均工资之间的相关系数，从而证明经济的发展趋势与工资的发展趋势一致。根据我国对21世纪中叶达到中等发达国家水平的预期，我国未来经

12、济增长率将会象现有的中等发达国家一样，保持较低的水平，即经济的发展将会从快速增长向缓慢增长转变，所以工资增长也会呈现这样的规律。为此，可以采用符合S型8增长特征的反正切函数作为预测模型3，利用山东省职工年平均工资的历史数据，采用最小二乘方法，拟合出模型的参数。对原有职工年平均工资数据进行处理，建立模型，进行预测拟合。图4.1.1我国1978年-2010年人均GDP4.2 问题二的分析主要是对数据进行计算以及处理，利用附件2所给的某企业分年龄职工数量及薪酬分布表，把2009年该企业各年龄段职工工资与该企业平均工资之比看作职工缴费指数的参考值，利用附件3提供的养老金计算办法，计算出养老金=个人账户

13、养老金+基础养老金，计算出替代率：4.3 问题三的分析问题三的提问方式与问题二的相似，是一道计算题，只是有给了我们一个具体的例子，只要把握好附件提供的数据资料，我们不难算出职工缴存的养老保险基金与其领取的养老金之间达到收支平衡的年龄。由于退休年龄有（55岁、60岁、65岁）3个阶段，所以要计算3次。利用：算出养老保险基金和养老金以及其达到收支平衡状态的年龄即可。4.4 问题四的分析问题四让我们提出措施，并给定理由。从两个方面出发：（1）既要达到目标替代率（2）又要维持养老保险基金收支平衡替代率的高低是影响养老保险基金收支平衡的一个重要因素，可见二者是有内在联系的，紧扣两个方面，从人口老龄化，职

14、工退休时的工资，退休前的养老金等方面进行周密考虑，而且养老保险基金收支平衡其实质就是政府的财政收支平衡，是一个动态平衡的环节（见图4.4.1），我们可以针对它们的影响因素来提出相应的措施，言之成理即可。20%“老人”的养老金企业缴费社会统筹基金过渡性养老金“中人”的养老金基础性养老金8%个人缴费“新人”的养老金个人账户养老金个人账户基金养老保险基金收入养老保险基金支出图 我国基本养老保险基金收入与支出示意图五、模型的建立与求解5.1 问题一 先将已有的我国人均GDP值6与职工年平均工资计算出相关系数为0.9985，进行线性拟合后如图 所示，证明了工资的增长与经济增长有显著的相关性。图拟合图这说

15、明工资的增长趋势与经济的发展趋势一致。而经济的发展将会从快速增长向缓慢增长转变，即会呈现S型增长的模式。所以工资增长也会呈现这样的规律。为此，采用符合S型增长特征的反正切函数作为预测模型。然后将附件1中的1978-2010年山东省职工年平均工资的数据用matlab做出散点图1，然后建立反正切拟合模型，求出职工年平均工资y。模型为：通过已有数据的散点图的走向分析，加上模型进行拟合，预测2011-2035年职工年平均工资得出2035年工资预测图，见图5.1.2，发现这个模型在前期的和原已有的实测数据拟合起来的效果成平滑的“S型”增长趋势，基本一致，由此可以证明这个模型建立的可行性，为后期预测从20

16、11年到2035年企业在岗职工工资的数据提供了科学依据。图5.1.2 2011-2035年工资预测图为了找出增长率最大的点，要求我们通过对模型求导得出一阶导数：然后再求出二阶导数并令其为0：找出拐点应为：用matlab计算出拐点是在2011年处，这与我国实际政策恰好相符，因为2010年是“十一五”规划的收官之年7，2008年的一场美国引发的次贷危机使得世界格局发生改变，中国是社会主义国家受到危机影响较小，中国开始扩大内需，刺激消费这架马车，我国年均GDP增长指标定是7.5%，实际在“十一五”期间，中国经济年均增长11.2%。在国民经济和社会发展第十二个五年规划纲要（草案）提出了今后五年我国经济

17、社会发展的主要目标：今后五年GDP年均增7%。由此可见2011年是一个增长速度最快的一年，具有事实依据。模型的具体形式是：并预测到2035年职工年平均工资为80063元。5.2 问题二 充分利用附件3数据，由养老金=个人账户养老金+基础养老金总公式出发，个人账户养老金=个人账户储存额计发月数基本养老金基础养老金=(参保人员退休时当地上年度在岗职工月平均工资+本人指数化平均缴费工资)/2缴费年限1%（缴费年限须满15年）个人账户养老金=个人账户储存额国家规定的计发月数（具体数据参考表）即用公式可以表示为：式中表示第年份在职职工的月平均缴费工资，和分别表示社会统筹账户和缴费比例与个人账户的缴费比例

18、。表各年龄段职工工资与该企业平均工资之比职工年龄段（岁）20-2425-2930-3435-39平均月收入（元）1726.9520772533.062732.3各年龄段职工工资与该企业平均工资之比0.670.810.991.07职工年龄段（岁）40-4445-4950-5455-59平均月收入（元）2987.563208.763064.812962.87各年龄段职工工资与该企业平均工资之比1.171.251.21.16其中职工平均工资2559.80元，各年龄段职工工资与该企业平均工资之比的平均数为1.04，根据各年龄段职工工资与该企业平均工资之比的平均数作为参考，再由“本人指数化月平均缴费工资

19、”指标的计算公式：计算出各年龄段的基础养老金。计算出的职工年龄段基础养老金如表.表职工年龄段每月基础养老金职工年龄段（岁）30-5530-6030-6540-5540-6040-65基础养老金(元)855.171060.351258.44382.25656.72897.86用个人账户储蓄额以及附件3表一查到的计发月数求出个人账户养老金，最后算出养老金。表个人账户养老金职工年龄段（岁）30-5530-6030-6540-5540-6040-65个人账户存储额（元）1030611349681747385905784666113642计发月数（个）170139101170139101个人账户养老金（

20、元）606.24970.991730.08347.39609.111125.17表职工年龄退休前工资职工年龄段（岁）30-5530-6030-6540-5540-6040-65职工退休前工资（元）2942.402832.922709.982942.402844.362709.98最终计算出该企业职工自2000年起分别从30岁、40岁开始缴养老保险，一直缴费到退休（55岁，60岁，65岁），各种情况下的养老金替代率（表）。表5.2.4 30岁、40岁开始缴养老保险到退休时的养老金替代率年龄段（岁）30-5530-6030-6540-5540-6040-65养老金（元）1461.412031.34

21、2988.51729.641265.832023.03替代率55%76%112%27%48%76%5.3 问题三 通过2009年山东省某企业各年龄段工资分布表计算出从30岁开始缴纳养老保险金到退休（55岁，60岁，65岁）的社会统筹基金的总数为：表 各年龄段社会统筹基金总额年龄段（岁）总收入（元）20%社会统筹基金总额（元）30-55871589.5820%174317.9230-601049361.8620%209872.3730-651219337.3420%243867.47然后根据表职工年龄段每月基础养老金中得到的数据求出每年基础养老金，再由：每年社会统筹基金每年养老金=每年剩余统筹基

22、金，最后按利率可以算出每年剩余养老金的本金和利息。三种不同情况如表所示。 表职工年龄段每月基础养老金职工年龄段（岁）30-5530-6030-6540-5540-6040-65基础养老金(元)6213.758486.506647.4244645682.085682.08表5.3.3 不同退休年龄出现的养老金收支平衡及缺口情况年龄(岁）剩余本金和利息（元）55岁退休60岁退休65岁退休55160343.880056145950.610057131125.560058115855.740059100127.84006083928.10192116.6006167242.36173828.16062

23、50056.06154991.0606332354.16135588.8506414121.21115604.57065-4523.0495020.77222426.7766-23167.2973819.45200342.8567-41811.5451982.09177596.4168-60455.7929489.61154167.5769-79100.046322.36130035.8870-97744.29-17029.04105180.2371-116388.54-40380.4479578.9172-135032.79-63731.8453209.5573-153677.04-8708

24、3.2426049.1274-172321.29-110434.64-1870.0375-190965.54-133786.04-29789.18最终得到结果发现，该职工自2000年起从30岁开始缴养老保险，一直缴费到退休（55岁，60岁，65岁），并从退休后一直领取养老金到65岁、70岁、74岁时，其缴存的养老保险基金与其领取的养老金之间达到收支平衡，并且出现缺口。由于该职工是75岁死亡，假设他在75岁那年还能领到养老金，那么缺口就在75岁的时候最大，分别为-190965.54元、-133786.04元、-29789.18元。 5.4 问题四伴随着我国人口老龄化问题的加剧，当下我国养老保险制

25、度主要遇到的难题集中表现在职工基本养老保险基金4供求严重失衡，资金缺口大，养老金支付能力差等难题。养老保险金收支缺口意味着在养老保险确保发放上财政处于赤字状态，基金支付能力不强，此问题决定着养老保险制度是否能持续、有效的运转下去。因此迫切需要提出既要达到目标替代率，又要维持养老保险基金收支平衡建议措施及理由：1、延长退休年龄5我国延长退休年龄的理由：国家对退休年龄的制定关系到养老金的收支情况、劳动者的工作年限等方面，所给合理确定退休年龄是至关重要的。当前我国的退休年龄仍然沿用上世纪50年代职工退休制度的规定，即男职工满60岁，女干部满55岁，女职工满50岁。我国在上世纪50年代规定的退休年龄在

26、当时是合理的，因为当时中国人口的平均预测寿命男性为40岁，女性为42.3岁。但是，到2000年我国人口的平均预期寿命男性和女性分别增长到69.0岁和72.4岁，如今我们依旧沿用以前的退休年龄制度，明显不合适了，因此我国的退休年龄与人均预期寿命延长趋势显得不相称。 我们考虑到人均受教育年限和人均预期寿命都有向后延伸的趋势，并且工作生命周期比下降、退休生命周期比上升对退休工作比存在拓展效应。为了缓解个人自我养老负担加重的趋势，就要提高工作生命周期比，降低退休工作比，也就得延长退休年龄。 退休年龄延长就个人而言：可以使个人有更多的工作时期的积累，以满足退休生命周期比的延长需要，也可以顺应人力资本储存

27、期延长的要求，同时相应延长其释放期，让个人可以更多的参与社会公共事务，为国家多做贡献；就国家而言：可以充分发挥劳动力资源和结构上的优势，将有利于经济增长。 减缓养老基金缺口。根据目前的人口状况，2020年将是人口老龄化的高峰期。延迟退休对于社会保险基金有很大好处：晚退5年，社保基金就可少知府年的养老金，同时又多加收5年的养老保险，一来一去就有10年养老保险基金的差距。从全国来看，我国退休年龄每延迟一年，养老统筹基金可增长40亿元，减支160亿元，减缓基金缺口200亿元。 人均受教育年限的增加带来的实际劳动起始年龄下限上移，同时，人均预期寿命的提高引发退休生涯扩张，实际劳动起始年龄后移、退休生涯

28、延长对工作年限存在挤压效应，这将给个人自我养老造成巨大压力。2、延长养老金缴费年限。理由：在新政策里，养老金的缴费年限与退休后的养老金领取额挂钩体现了激励机制的作用。按照新规定，缴费不满15年的，退休后不享受基础养老金待遇；满15年后，每超过1年，基础养老金的领取额就将增加1。不仅如此，而且养老金调整的幅度高低也与缴费年限的长短直接挂钩，如果缴费年限只有15年，那么在调整养老金待遇时始终都要受到影响。因此，缴费年限越长，缴费额就越大，个人账户的总积累额就越多，退休后领取的基础养老金也就越多。3、改革现有的按工资指数化的管理方法，以相应地减轻养老金支付负担和转制成本，从而缓解养老基金收支不平衡的

29、状况。理由：由于职工缴纳的保险费与领取的保险费均按一定的比例与工资挂钩，因此一旦工资增长后，职工所缴纳与所领取的保险费也会随着增长。所以降低工资增长率可以缓解养老金支付负担。但是工资增长率是很难降低的，因为经济学原理表明在市场经济下生产率的增长最终将导致劳动需求的增长，以及一般工资水平的提高，而工资的增长又富有刚性。所以要通过降低工资的增长率来缩减养老金负担是很难实现的。但若我们从另外一个角度来分析的话，即养老金的发放额不与工资挂钩、不按工资进行指数化的话，我们仍能在工资增长时实现减轻养老金负担的目标。这时，养老金管理就应从与工资增长指数挂钩转向与消费者物价指数挂钩，这种转换，理论分析其实也是

30、可行的。因为养老金管理的首要目标是使其保值，保证退休人员生活水平不下降，所以首先应与消费者物价指数挂钩。过去20年，我国工资增长一直较快，年平均57，预计今后20年仍将以3一5的速度增长。如此高的增长速度，大大增加了与其挂钩的养老金负担和转制成本。既然理论分析表明，养老金应首先与消费者物价指数挂钩，使其保值，那么我们就应该尽快改革现有的按工资指数化的管理方法，以相应地减轻养老金支付负担和转制成本，从而缓解养老基金收支不平衡的状况。4、建立专门部门管理职工上缴个人账户养老金，为基金做投资增加基金收益率。理由：据有关统计资料显示，截至2006年底，我国养老保险基金的收益率达到了9.6。但是在从20

31、03年至2007年的5年问，我国养老保险基金名义投资收益率为2.18，而同期加权通胀率为2.29。可见，养老基金的实际收益率是在贬值缩水。而低的基金收益率是目前我国部分积累基金制度下收支不平衡的重要原因。可以说，决定缴费率的最主要的因素便是养老金投资的实际收益率。因为当养老金支付水平相同时，较高的基金收益率可以要求较低的缴费率。当前，我国由于受到投资政策及风险管理水平的约束，养老基金的80左右必须投到国有银行或购买国家债券，这样便无法保证较高的收益率，相应也就要求较高的缴费率，而高缴费率又造成低覆盖率和低收缴率，从而间接地影响了养老保费的收入总额，大大弱化了部分积累式的自我保障功能，甚至导致支

32、付危机。再者，由于我国的经济运行存在着一定的通货膨胀率，这又使得原本不多的低收益率的养老基金积累显得更少，从而进一步加速了养老金支付的困难。因此，必须加强对基金收益与风险辩证关系的研究，在保证基金安全的前提下，采取科学有效的措施，逐步拓宽养老基金的投资渠道，提高养老基金投资的收益率。六、模型的分析与评价2 我们对题目所给数据进行了充分利用，与模型进行拟合，挖掘其中隐含拐点的各因素变化趋势，对题目要求解决的每一个问题都做了比较彻底、圆满的解答。所给的反正切模型是一个对称模型，也是一个典型的“S”型模型，与实际中国经济发展势头基本吻合，是一个较为完美、经典的模型。其模型的优点：1、思路清晰、简明易

33、懂、彻底全面。模型建立使得预测的中国未来几十年的经济发展趋势形象表达出来。2、对四个问题解答环环相扣。上一问题的解答为下一问题的解答提供依据。3、体现了中国经济发展呈S型增长的基本规律。4、该模型对回归系数统计显著，符合问题一的要求。模型的不足之处依然存在：现实存在诸多不确定的因素，例如政府对每月缴纳养老金的金额和数量或者年龄利率有所改变时，可能会影响假设的不成立或者计算结果不合理。根据我们的预测模型得到的年平均工资的预测结果，在2016年以后增长率就下降到5%以下，到2025年以后甚至低于1%。这表明模型的预测结果在2016年后增长非常缓慢。这说明受到模型本身曲线形态的局限性，参考文献1 周

34、品，赵新芬，MATLAB数学建模与仿真M，北京：国防工业出版社，2009.2 赫孝良，戴永红，周义仓，数学建模竞赛赛题简析与论文点评，西安：西安交通大学出版社，2002.3 刘承平，数学建模方法，北京：高等教育出版社，2002.4 穆怀中，柳清瑞等，中国养老保险制度改革关键问题研究M，北京：中国劳动社会保障出版，2006.5 毕小龙，中国社会养老保险制度:经济转型、人口老龄化与社会养老保险M，广州：暨南大学出版社，2009.6 中华人民共和国国家统计局，7 腾讯财经聚焦十二五，十一五我国GDP年均实际增长11.4% 国力更强，8 姜启源，数学模型M，第3版，北京：高等教育出版社，2003.9

35、新浪财经，养老金空涨超9000亿元实际收益率贬值缩水, 附件表一：个人账户养老金计发月数表退休年龄计发月数退休年龄计发月数402335616441230571584222658152432235914544220601394521661132462126212547208631174820464109491996510150195669351190678452185687553180696554175705655170山东省职工历年平均工资统计表年份平均工资1978566 1979632 1980745 1981755 1982769 1983789 1984985 19851110 1986

36、1313 19871428 19881782 19891920 19902150 19912292 19922601 19933149 19944338 19955145 19965809 19976241 19986854 19997656 20008772 200110007 200211374 200312567 200414332 200516614 200619228 200722844 200826404 200929688 201032074 注：1.2010年前的数据来源于山东省统计年鉴2010 x=1:32;x=1:32;y=381 419 463 492 528 583 6

37、95 858 963 1112 1366 1519 1644 1893 2311 2998 4044 .5046 5846 6420 6796 7159 7858 8622 9398 10542 12336 14185 16500 20169 23708 25575;fun=inline(a(1)*atan(a(2)*x-a(3)+a(4),a,x);a,J=lsqcurvefit(fun,20000,3,30,20000,x,y);z=a(1)*atan(a(2)*x-a(3)+a(4);x1=1:58;z=a(1)*atan(a(2)*x1-a(3)+a(4);plot(x,y,k),ho

38、ld on,plot(x1,z,k:)预测2035年工资预测模型编程x=1:33;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 ;fun=inline(a(1)*atan(a(2)*x-a(3)+a(4),a,x);a,J=lsqcurvefit(fun,20000,1,20,45500,x,y)z=a(

39、1)*atan(a(2)*x-a(3)+a(4);x1=1:58;z=a(1)*atan(a(2)*x1-a(3)+a(4);plot(x,y,k-o),hold on,plot(x1,z,k:)年平均工资与GDP之间的先行关系拟合图a = 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 bb = 381 419 463

40、492 528 583 695 858 963 1112 1366 1519 1644 1893 2311 2998 4044 5046 5846 6420 6796 7159 7858 8622 9398 10542 12336 14185 16500 20169 23708 25575 plot(a,b,o) p=polyfit(a,b,1);q1=p(1)*a+p(2); hold on;plot(a,q1,g)年份预测全省年工资水平退休前一年全省年工资水平c1/cm每个年龄阶段的年工资缴费额水平xm*(c1/cm)20008772778388.8734612431.6821577.41767200110007778387.7783552431.6818914.47065200211374778386.8435032431.6816641.20871200312567778386.1938412431.6815061.43931200414332778385.4310632431.6813206.60814200516614778384.6850852622.7212287.6657920061