s174.m

来自「matlab 设备更新问题」· M 代码 · 共 36 行

function[u,I]=s174(m,N,r,a,c)
%m为第一年所所使用的设备在年初时的役龄
%N为所考虑的年数
%r为具有不同役龄的设备使用一年可得到的经济收入向量
%a为具有不同役龄的设备使用一年所需的维修费用向量
%c为卖掉具有不同役龄的设备同时买进一台新设备的纯支出费用向量
for i=1:m+N
    V(i,N+1)=0;
end
for k=N:-1:2
    for i=1:m+k-1
       if r(i+1)-a(i+1)+V(i+1,k+1)<r(1)-a(1)-c(i+1)+V(1,k+1)
           V(i,k)=r(1)-a(1)-c(i+1)+V(1,k+1);
           w(i,k)=0;
       else
           V(i,k)=r(i+1)-a(i+1)+V(i+1,k+1);
           w(i,k)=1;
       end
    end
end
if r(m+1)-a(m+1)+V(m+1,2)<r(1)-a(1)-c(m+1)+V(1,2)
    I =r(1)-a(1)-c(m+1)+V(1,2);
    u(1)=0;
else
     I = r(m+1)-a(m+1)+V(m+1,2);
     u(1)=1;
end
x(1)=m;
for k=2:N
    x(k)=x(k-1)*u(k-1)+1;
    u(k)=w(x(k),k);
end