实验5 动态规划.docVIP

年级、专业 2011应数2班 姓名 学号 名单序号 实验时间 2013年 5 月 16 日 MATLAB版本： R2010b 注：实验报告的最后一部分是实验小结与收获 实验五 动态规划 PAGE 2013春 数学实验 实验五 动态规划 PAGE \* Arabic 8 / NUMPAGES \* Arabic 8 1. 某公司打算向它的四个营业区增设4个销售点，从各区赚取的利润与增设的销售店个数有关，其数据为：? 第1题数据表（数据文件：动态规划 Prob_1.xls ） 销售店 增加数 第一区利润 （万元） 第二区利润 （万元） 第三区利润 （万元） 第四区利润 （万元） 0 160 190 200 250 1 310 225 298 308 2 541 445 399 487 3 600 517 601 655 4 705 632 721 674 ? 试求各区应分配几个增设的销售店，才能使总利润最大？其值是多少？ %以下各题均要应用此逆序算法函数： dynprog函数 function [p_opt,fval]=dynprog(x,DecisFun,SubObjFun,TransFun,ObjFun) % x为状态变量，一列代表一个阶段的状态 % M_函数DecisFun（k,x)表示由阶段k的状态值x求出相应的允许决策集合 % M_函数SubObjFun（k,x,u)表示阶段k的指标函数 % M_函数TransFun（k,x,u)是状态转移函数，其中x是阶段k的状态值，u是其决策集合 % M_函数ObjFun（v,f)是第k阶段到最后阶段的指标函数，当ObjFun（v,f)=v+f时，输入ObjFun（v,f)可以省略 % 输出p_opt由4列组成，p_opt=[序号组，最优轨线组，最优策略组，指标函数值组]; % 输出fval是列向量，各元素分别表示p_opt各最优策略组对应始端状态x的最优函数值 k=length(x(1,:)); % k为阶段数 x_isnan=~isnan(x); t_vub=inf; t_vubm=inf*ones(size(x)); % t_vubm为指标函数值的上限 f_opt=nan*ones(size(x)); % f_opt为不同阶段、状态下的最优值矩阵，初值为非数 d_opt=f_opt; % d_opt为不同阶段不同状态下的决策矩阵，初值为非数 tmp1=find(x_isnan(:,k)); % 找出第k阶段状态值（不是非数）的下标 tmp2=length(tmp1); for i=1:tmp2 u=feval(DecisFun,k,x(tmp1(i),k));% 求出相应的允许决策向量 tmp3=length(u); for j=1:tmp3 % 该for语句是为了求出相应的最有函数值以及最优决策 tmp=feval(SubObjFun,k,x(tmp1(i),k),u(j)); if tmp=t_vub f_opt(i,k)=tmp; d_opt(i,k)=u(j); t_vub=tmp; end end end for ii=k-1:-1:1 % 从后往前面递推求出f_opt 以及d_opt tmp10=find(x_isnan(:,ii));tmp20=length(tmp10); for i=1:tmp20 u=feval(DecisFun,ii,x(i,ii)); tmp30=length(u); for j=1:tmp30 tmp00=feval(SubObjFun,ii,x(tmp10(i),ii),u(j)); tmp40=feval(TransFun,ii,x(tmp10(i),ii),u(j)); % 由该状态值及相应的决策值求出下一阶段的状态值 tmp50=x(:,ii+1)-tmp40; tmp60=find(tmp50==0); % 找出下一阶段的状态值在x(:,ii+1)的下标 if ~isempty(tmp60) if nargin5,tmp00=tmp00+f_opt(tmp60(1),ii+1); else