function [u,I]=s173(a,b,m,N,d,T)
% N 为阶段数
% a 为初始库存量,非负整数
% b 为终止库存量,非负整数
% m 为最大生产能力向量,非负整数向量
% d 为社会需求向量,非负整数向量
% T 为最大库存向量,非负整数向量
x(1)=a;
x(N+1)=b;
if d(N)+x(N+1)<T(N)
s=d(N)+x(N+1);
else
s=T(N);
end
for i=1:s+1
V(i,N)=cost(x(N+1)+d(N)-i+1,N,m)+stochar(i-1,N);
W(i,N)=x(N+1)+d(N)-i+1;
end
for k=N-1:-1:2
z=x(N+1);
for i=k:N
z=z+d(i);
end
if z>T(k)
s(k)=T(k);
else
s(k)=z;
end
for j=1:s(k)+1
if z-j+1>m(k)
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余2页未读,立即下载