www.gusucode.com > 司守奎《Lingo软件及应用习题解答》的源程序 > 司守奎《Lingo软件及应用习题解答》的源程序/《Lingo软件及应用习题解答》程序和数据/12第12章/Lt12_5_5.m

    clc, clear, close all
r=0.03; r2=0.08;
gn=[170	164	158	152	145	139	132	125	117	109	101]; %55岁到65岁退休时，个人账户养老金计发月数
m=[25:35]; %30岁开始缴费时的11种缴费年限;
c=[8772   10007 	11374     12567 	 14332    16614 	 19228 ...    
22844   26404    29688      32074     36996	    41699    46871	 ...
52526     58670  65300	   72402    79949	    87903    96211   ...  
104810	 113624   122572   131564   140513	    149329    157930 ...  
166241   174198   181750   188856   195490   201639	 207298]; %从2000年到2034年35年的平均工资;
L=[0.9825	 1.0667	 1.1728	 1.2666	 1.2087	 1.1550   1.1002]; %从30岁到64岁段的缴费指数;
L=repmat(L,5,1); L=L(:); %从30岁到64岁每一年的缴费指数
LL=cumsum(L); LL=LL(25:end)./[25:35]' %从55岁退休到65岁退休的总平均缴费指数
for j=1:length(c)
    x(j)=c(j)*L(j); %根据缴费指数和平均工资，计算职工各年度的工资
end
for i=1:length(LL)
    sa(i)=(c(m(i))*(1+LL(i))/24*m(i)*0.01+sum(x(1:m(i))*r2.*(1+r).^([1:m(i)])./gn(i)))*12/x(m(i));
end
plot([55:65],sa,'.-k'), hold on
plot([55,65],[0.585,0.585],'--k')
text(62.5,0.6,'0.585')
xlabel('退休年龄'), ylabel('养老金替代率')