毕业论文 下载的matlab程序 但是mat文件下载不到a01-05.mat，谁能通过下面程序帮把a01-05.mat
算出来，谢谢了

load a01-05.mat
for r=1:91
    p1(r,1)=(a2001(r,2)+a2001(r,4));
    p2(r,1)=(a2001(r,6)+a2001(r,8));
    p3(r,1)=(a2001(r,10)+a2001(r,12));
    p(r,1)=(p1(r,1)*(z01_05(2,1)+z01_05(3,1))+p2(r,1)*(z01_05(4,1)+z01_05(5,1))+p3(r,1)*(z01_05(6,1)+z01_05(7,1)))/(z01_05(2,1)+z01_05(3,1)+z01_05(4,1)+z01_05(5,1)+z01_05(6,1)+z01_05(7,1));
   
    p1(r,2)=(a2002(r,2)+a2002(r,4));
    p2(r,2)=(a2002(r,6)+a2002(r,8));
    p3(r,2)=(a2002(r,10)+a2002(r,12));
    p(r,2)=(p1(r,2)*(z01_05(2,2)+z01_05(3,2))+p2(r,2)*(z01_05(4,2)+z01_05(5,2))+p3(r,2)*(z01_05(6,2)+z01_05(7,2)))/(z01_05(2,2)+z01_05(3,2)+z01_05(4,2)+z01_05(5,2)+z01_05(6,2)+z01_05(7,2));
   
    p1(r,3)=(a2003(r,2)+a2003(r,4));
    p2(r,3)=(a2003(r,6)+a2003(r,8));
    p3(r,3)=(a2003(r,10)+a2003(r,12));
    p(r,3)=(p1(r,3)*(z01_05(2,3)+z01_05(3,3))+p2(r,3)*(z01_05(4,3)+z01_05(5,3))+p3(r,3)*(z01_05(6,3)+z01_05(7,3)))/(z01_05(2,3)+z01_05(3,3)+z01_05(4,3)+z01_05(5,3)+z01_05(6,3)+z01_05(7,3));
   
    p1(r,4)=(a2004(r,2)+a2004(r,4));
    p2(r,4)=(a2004(r,6)+a2004(r,8));
    p3(r,4)=(a2004(r,10)+a2004(r,12));
    p(r,4)=(p1(r,4)*(z01_05(2,4)+z01_05(3,4))+p2(r,4)*(z01_05(4,4)+z01_05(5,4))+p3(r,4)*(z01_05(6,4)+z01_05(7,4)))/(z01_05(2,4)+z01_05(3,4)+z01_05(4,4)+z01_05(5,4)+z01_05(6,4)+z01_05(7,4));
   
    p1(r,5)=(a2005(r,2)+a2005(r,4));
    p2(r,5)=(a2005(r,6)+a2005(r,8));
    p3(r,5)=(a2005(r,10)+a2005(r,12));
    p(r,5)=(p1(r,5)*(z01_05(2,5)+z01_05(3,5))+p2(r,5)*(z01_05(4,5)+z01_05(5,5))+p3(r,5)*(z01_05(6,5)+z01_05(7,5)))/(z01_05(2,5)+z01_05(3,5)+z01_05(4,5)+z01_05(5,5)+z01_05(6,5)+z01_05(7,5));
end;

krt=0.460012;
bata=1.8;
swl=sum(a2005(:,2).*a2005(:,3))+sum(a2005(:,4).*a2005(:,5));

%×ÛºÏ
t=0:50;
%•Ö²¼Òò×Ó
io1=0.002*t+0.25;
io2=0.003*t+0.15;
io3=-0.007*t+0.6;
%ÀÏÁ仯µ¼ÖÂËÀÍöÂÊÔö¼ÓÒò×Ó
lp=0.000019212*t;
%ËÀÍöÂʽµµÍÒò×Ó
sw=0.000034523*t;

%×ÛºÏ
for t=0:50;
    f(t+1)=0;
    for r=0:90
        if t<=r
             jf1=((sum(a2005((r-t+1)r+1),2).*a2005((r-t+1)r+1),3))+sum(a2005((r-t+1)r+1),4).*a2005((r-t+1)r+1),5)))*io1(t+1)+(sum(a2005((r-t+1)r+1),6).*a2005((r-t+1)r+1),7))+sum(a2005((r-t+1)r+1),8).*a2005((r-t+1)r+1),9)))*io2(t+1)+(sum(a2005((r-t+1)r+1),10).*a2005((r-t+1):(r+1),11))+sum(a2005((r-t+1):(r+1),12).*a2005((r-t+1):(r+1),13)))*io3(t+1))/100000-lp(t+1)+sw(t+1);
            if t>4
                 p(r+1,t+1)=p(r-t+1,1)*exp(-jf1);
            end;
        end;
        if t>r
            jf2=((sum(a2005(1:(r+1),2).*a2005(1:(r+1),3))+sum(a2005(1:(r+1),4).*a2005(1:(r+1),5)))*io1(t+1)+(sum(a2005(1:(r+1),6).*a2005(1:(r+1),7))+sum(a2005(1:(r+1),8).*a2005(1:(r+1),9)))*io2(t+1)+(sum(a2005(1:(r+1),10).*a2005(1:(r+1),11))+sum(a2005(1:(r+1),12).*a2005(1:(r+1),13)))*io3(t+1))/100000-lp(t+1)+sw(t+1);
            if t>4
                p(r+1,t+1)=f(t-r+1)*exp(-jf2);
            end;
        end;
    end;
    for j=15:49
        f(t+1)=f(t+1)+(a2005(j+1,15))*krt*p(j+1,t+1)/100000;
    end;
    f(t+1)=bata*f(t+1);
    if f(t+1)<0.002
        f(t+1)=0.002;
    end;
    zrzzl(t+1)=f(t+1)-swl/100000-lp(t+1)+sw(t+1);
end;

rk(1)=127627;
rk(2)=128453;
rk(3)=129227;
rk(4)=129988;
rk(5)=130756;
rk(6)=131448;

for i=6:50
    rk(i+1)=rk(i)*(1+zrzzl(i+1));
end;

t=0:50;
plot(t,rk,'*');
grid on;
hold on;