一、部分题目
2011高教社杯全国大学生数学建模竞赛题目
C题 企业退休职工养老金制度的改革
养老金也称退休金,是一种根据劳动者对社会所作贡献及其所具备享受养老保险的资格,以货币形式支付的保险待遇,用于保障职工退休后的基本生活需要。
我国企业职工基本养老保险实行“社会统筹”与“个人账户”相结合的模式,即企业把职工工资总额按一定比例(20%)缴纳到社会统筹基金账户,再把职工个人工资按一定比例(8%)缴纳到个人账户。这两个账户我们合称为养老保险基金。退休后,按职工在职期间每月(或年)的缴费工资与社会平均工资之比(缴费指数),再考虑到退休前一年的社会平均工资等因素,从社会统筹账户中拨出资金(基础养老金),加上个人工资账户中一定比例的资金(个人账户养老金),作为退休后每个月的养老金。养老金会随着社会平均工资的调整而调整。如果职工死亡,社会统筹账户中的资金不退给职工,个人账户中的余额可继承。个人账户储存额以银行当时公布的一年期存款利率计息,为简单起见,利率统一设定为3%。
养老金的发放与职工在职时的工资及社会平均工资有着密切关系;工资的增长又与经济增长相关。近30年来我国经济发展迅速,工资增长率也较高;而发达国家的经济和工资增长率都较低。我国经济发展的战略目标,是要在21世纪中叶使我国人均国民生产总值达到中等发达国家水平。
现在我国养老保险改革正处于过渡期。养老保险管理的一个重要的目标是养老保险基金的收支平衡,它关系到社会稳定和老龄化社会的顺利过渡。影响养老保险基金收支平衡的一个重要因素是替代率。替代率是指职工刚退休时的养老金占退休前工资的比例。按照国家对基本养老保险制度的总体思路,未来基本养老保险的目标替代率确定为58.5%. 替代率较低,退休职工的生活水准低,养老保险基金收支平衡容易维持;替代率较高,退休职工的生活水准就高,养老保险基金收支平衡较难维持,可能出现缺口。所谓缺口,是指当养老保险基金入不敷出时出现的收支之差。
附件1是山东省职工历年平均工资数据;附件2是2009年山东省某企业各年龄段职工的工资分布情况,附件3是养老金的计算办法。请建立数学模型,解决如下问题:
问题1 对未来中国经济发展和工资增长的形势做出你认为是简化、合理的假设,并参考附件1,预测从2011年至2035年的山东省职工的年平均工资。
问题2 根据附件2计算2009年该企业各年龄段职工工资与该企业平均工资之比。如果把这些比值看作职工缴费指数的参考值,考虑该企业职工自2000年起分别从30岁、40岁开始缴养老保险,一直缴费到退休(55岁,60岁,65岁),计算各种情况下的养老金替代率。
问题3 假设该企业某职工自 2000年起从30岁开始缴养老保险,一直缴费到退休(55岁,60岁,65岁),并从退休后一直领取养老金,至75岁死亡。计算养老保险基金的缺口情况,并计算该职工领取养老金到多少岁时,其缴存的养老保险基金与其领取的养老金之间达到收支平衡。
问题4 如果既要达到目标替代率,又要维持养老保险基金收支平衡,你认为可以采取什么措施。请给出你的理由。![请添加图片描述](https://i-blog.csdnimg.cn/direct/4deb1d50343344c19ce0740fc16750f3.png)
二、部分论文
三、部分源代码
程序1
clear all;clc;
%年平均工资
y0=[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];
%人均国民生产总值
GNP0=[318 419 463 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 29667];
X=[ones(33,1),GNP0'];
[b,bint,r,rint,stats]=regress(y0',X)
GNP=300:30000;
y=1.12*GNP+89.07;
plot(GNP0,y0,'o',GNP,y);
grid on;
yw=1.12*51109+89.07
程序2
clear all;clc;
ymax=1.12*8000*6.3886+89.07;
a=1/(ymax);
syms y c x
t1=33;y1=32074;
b1=(1-a*y1)/(y1*exp(-c*t1))
t2=1;y2=566;
b2=(1-a*y2)/(y2*exp(-c*t2))
%运行后的b1,b2值
%b1=248005506055589/(18056056806035161088*exp(-33*c));
%b2=278696112876881/(159314836818231296*exp(-c));
a=248005506055589;b=18056056806035161088;
c=278696112876881;d=159314836818231296;
C0=-log(a*d/(b*c))/32
b0=248005506055589/(18056056806035161088*exp(-33*C0))
程序3
%Logistic regression analysis model
clear all;clc;
t=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 ];
b0=[0.0020;0.1515];
[b,R,J]=nlinfit(t,y,'fun',b0)
%参数b估计的置信区间
ci=nlparci(b,R,J);
ymax=1.12*8000*6.3886+89.07;
a=1/(ymax);
t1=1:68;
f=1./(a+b(1)*exp(-b(2)*t1));
plot(t,y,'o',t1,f);grid on;
fw=f(34:end)
M函数文件:
function f=fun(b,t)
ymax=1.12*8000*6.3886+89.07;
a=1/(ymax);
f=1./(a+b(1)*exp(-b(2)*t));
程序4
function y=PA(c,m,r, ) %个人账户储存额
p=m:-1:1;
R=(1+ ).^p;
R=R';
y=c*r*r2;
程序5
clear all;clc;
%月收入平均值
ygz=[ 1249.5,1749.5,2249.5,2749.5,3249.5,3749.5,4499.5,6500.0];
%各年龄段人数
rs=[74 165 26 16 1 0 0 0
36 82 94 42 6 3 0 0
0 32 83 95 24 6 2 0
0 11 74 83 36 16 4 2
0 0 43 86 55 21 13 3
0 3 32 32 64 41 18 4
0 7 23 29 44 21 8 3
0 6 17 27 37 7 7 0];
%企业总人数
zrs=sum(sum(rs));
A=zeros(8,8);
for i=1:8
A(:,i)=ygz(i)*rs(:,i);
end
%企业总工资
zgz=sum(sum(A));
%企业月平均工资
ypjgz=zgz/zrs;
%各年龄段总人数
C=[282 263 242 226 221 194 135 101];
%各年龄段各类工资
D=[ 92463 288667.5 58487 43992 3249.5 0 0 0
44982 143459 211453 115479 19497 1124.8 0 0
0 55984 186708.5 261202.5 77988 22497 8999 0
0 1924.5 166460 228210 116980 59990 18000 13000
0 0 96730 236460 178720 78740 58490 19500
0 5250 71980 87980 207970 153730 80990 26000
0 12250 51740 79740 142980 78740 36000 19500
0 10500 38240 74240 120230 26250 31500 0];
%各年龄段总工资
E=zeros(1,8);
for i=1:8
E(i)=sum(D(i,:));
end
%各年龄段职工工资
for i=1:8
F(i)=E(i)/C(i);
end
%各年龄段职工平均工资与企业平均工资之比
lamda=F./ypjgz