s174.m

来自「该文件主要包括了matlab的一些主要应用例程!对学习matlab很有帮助」· M 代码 · 共 33 行

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