matlab解动态规划,数太大解不出来
function u=DecisF1(k,x) %在阶段k由状态变量x的值求出其相应的决策变量所有的取值?
if k==16,u=x;
else u=1:x-1;
end
function v=SubObjF1(k,x,u) %阶段k的指标函数?
c=[345.78;358;571.5;196.5;345.78;358;524.5;19.5;345.78;358;396.5;28.5;345.78;358;1185.5;18.5]
v=c(k).*u;
v=-v;
function y=TransF1(k,x,u) %状态转移方程?
y=x-u;
function y=ObjF1(v,f) %基本方程中的函数g?
y=v+f;
clear;
n=37350;
x1=[n;nan*ones(n-1,1)];
x2=1:n;
x2=x2';
x=[x1,x2,x2,x2,x2,x2,x2,x2,x2,x2,x2,x2,x2,x2,x2,x2,];
[p,f]=dynprog(x,'DecisF1','SubObjF1','TransF1','ObjF1')