动态规划二维分配问题-MATLAB.doc

姓名评分 实验报告 课程名称： 动态规划 实验名称： 二维分配问题 专 业： 信息与计算科学 指导教师： 滕宇 完成日期： 2014年 11月 07日 MATLAB程序： clc;clear; P=[0.4 0.1 0.5; 0.2 0.4 0.2]; A=[3 2 5]; M=[6 10]; sv=@(u,w,pi)A(pi)*(1-((1-P(1,pi))^u)*((1-P(2,pi))^w));% 朝第pi目标发送u个第一种导弹，w个第二种导弹，的价值; lmd=0.7; for ll=1:99 lmd=ll/100; v=@(u,w,pi)sv(u,w,pi)-lmd*w; vu=zeros(7,3);% vk(uk)的值 wi=vu;% 相应的决策 for l=1:3 lv=@(u,w)v(u,w,l); for m=0:6 mv=@(w)lv(m,w); for n=0:10 ff=mv(n); if vu(m+1,l)ff vu(m+1,l)=ff; wi(m+1,l)=n; end end end end fx=zeros(7,3); ui=fx; fx(:,3)=vu(:,3); ui(:,3)=0:6; for i=1:2 l=3-i; for m=0:6 for n=0:m ff=vu(n+1,l)+fx(m-n+1,l+1); if fx(m+1,l)ff fx(m+1,l)=ff; ui(m+1,l)=n; end end end end %% 动态规划 x=zeros(1,3); u=x; w=x; x(1)=find(fx(:,1)==max(fx(:,1)))-1; u(1)=ui(x(1)+1,1); x(2)=x(1)-u(1); u(2)=ui(x(2)+1,2); x(3)=x(2)-u(2); u(3)=ui(x(3)+1,3); w(1)=wi(u(1)+1,1); w(2)=wi(u(2)+1,2); w(3)=wi(u(3)+1,3); %% 判断是否符合 if sum(u)==6sum(w)==10 disp(lmd符合 最大价值为); max(fx(:,1)) disp(方案为); u w lmd else disp(lmd不符合); end end 最后的运行结果为：