求动态规划的MATLAB程序，要通用性强的

majiancumt

求动态规划的MATLAB程序，要通用性强的如dynprog

无语的昭昭

额额额，其他地方很多啊
1 D8 X, M7 W) V8 t- z

张坤14

hugoczx 发表于 2009-12-21 23:20
有是有                              
1 V/ A$ Y$ ^. b4 p% d& `- u不过不能白给！

同求动态规划的源程序，现在其他的算法基本上都可以集成，就是动态规划没有了。希望大神赐教，我的qq邮箱1240258213@qq.com，谢谢。关于其他算法的集成， 可以与大神共享。

m存在

madio

1. function [p_opt,fval]=dynprog(x,DecisFun,ObjFun,TransFun)      
   ! j- Q& U. A7 }0 \
2. % [p_opt,fval]=dynprog(x,DecisFun,ObjFun,TransFun)9 C  a5 l: P% ^8 S# S) x1 d
   
3. % 自由始端和终端的动态规划,求指标函数最小值的逆序算法递归
   4 i0 ]/ O% |8 i# p; v0 v+ X
4. % 计算程序。x是状态变量，一列代表一个阶段状态；M-函数$ w8 ]  L8 ]( K: `8 v7 \6 z
   
5. % DecisFun(k,x)由阶段k的状态变量x求出相应的允许决策变量;6 B( \9 T6 c0 M9 s
   
6. % M-函数ObjFun(k,x,u)是阶段指标函数，M-函数TransFun(k,x,u)3 `- L3 T, y4 H$ I
   
7. % 是状态转移函数,其中x是阶段k的某状态变量，u是相应的决策变量；
   % a5 k\" C: J+ ?
8. % 输出p_opt由4列构成，p_opt=[序号组;最优策略组;最优轨线组;
   - L2 S/ ~# p+ f, d( N+ F2 Z
9. % 指标函数值组]；fval是一个列向量，各元素分别表示p_opt各2 H: b3 t/ I4 E3 M
   
10. % 最优策略组对应始端状态x的最优函数值；$ d, o6 b; p6 K8 J/ l: L
    
11. %' _6 \\" w$ V. [+ j/ ^
    
12. %例（参看胡良剑等编《数学实验--使用MATLAB》P180+ E2 k! L4 e5 p3 [( c/ {2 D
    
13. %先写3个函数
    5 {+ {) Y: s5 m7 \. T
14. %                 eg13f1_2.m/ |* H- I; L5 h8 _& O$ Z
    
15. %    function u=DecisF_1(k,x)
    ' h! A, z, H+ {# {7 Z; T
16. %    在阶段k由状态变量x的值求出其相应的决策变量所有的取值5 F: O6 U# @# A6 B' F
    
17. %    c=[70,72,80,76];q=10*[6,7,12,6];
      _4 z! y$ g5 T9 R- Q+ i. ~6 ?
18. %    if q(k)-x<0,u=0:100;       %决策变量不能取为负值
    . K/ \% P5 x  W5 B
19. %    else,u=q(k)-x:100;end;     %产量满足需求且不超过100/ {\" b* A2 K8 t/ ~8 @- R( E' K4 V/ m3 h
    
20. %    u=u(:);
    , M+ m' J: `0 R/ J) E( p
21. %                 eg13f2_2.m- q: Z: e# ^6 U6 E( K9 O  @9 w. u
    
22. %    function v=ObjF_1(k,x,u)
    % N  k  W! I- b# g
23. %    阶段k的指标函数
    ( T2 |$ q5 `5 S3 {
24. %    c=[70,72,80,76];v=c(k)*u+2*x;/ @5 c; ^5 y0 M) {& x/ O
    
25. %                 eg13f3_2.m0 Q3 ?/ o# k3 K0 r- S
    
26. %    function y=TransF_1(k,x,u)8 `- g+ b) Y& C; K; z7 R7 j
    
27. %     状态转移方程
    / E# j4 `+ W8 R# D4 W\" ]2 ^# }0 V
28. %     q=10*[6,7,12,6];y=x+u-q(k);0 ~9 |: ]! {8 e) t+ y7 T
    
29. %调用DynProg.m计算如下：
    ' L5 e% s- n0 V
30. %    clear;x=nan*ones(14,4);% x是10的倍数，最大范围0≤x≤130,
    - O* R6 I3 n4 l3 t1 u
31. %       %因此x=0,1,...13，所以x初始化取14行，nan表示无意义元素
    + M9 w4 Y# u! |# c
32. %    x(1:7,1)=10*(0:6)';     % 按月定义x的可能取值8 h$ Q& v; S/ C
    
33. %    x(1:11,2)=10*(0:10)';x(1:12,3)=10*(2:13)';$ e. Z4 f& p) y6 p- y\" m; `
    
34. %    x(1:7,4)=10*(0:6)';
    ( }) x. b' |9 l. e
35. %    [p,f]=dynprog(x,'eg13f1_2','eg13f2_2','eg13f3_2')
    ' g( N: a) u$ E/ ]\" Q
36. 
    ; S\" l+ p, M% _0 ?& l! Y' r
37. % By X.D. Ding June 2000
    1 K6 p5 D9 I, w+ x
38. 
    9 N; W! t, I; q* g8 L
39. k=length(x(1,:));f_opt=nan*ones(size(x));d_opt=f_opt;/ L9 T! }8 d1 H% A9 u: I) e- i' m
    
40. t_vubm=inf*ones(size(x));x_isnan=~isnan(x);t_vub=inf;
    5 l- h! N\" A/ K2 N! }5 M$ ^  n
41. % 计算终端相关值, u. S3 L/ [# ~4 c. ?- L; Y
    
42. tmp1=find(x_isnan(:,k));tmp2=length(tmp1);
    & A% \: U' Z. |
43. for i=1:tmp2
    . e! L0 p0 g* o6 Y' V4 ]
44.    u=feval(DecisFun,k,x(i,k));tmp3=length(u);1 U8 b. O. R: C
    
45.    for j=1:tmp3
    * F* j& w\" `7 J& |\" @3 n
46.          tmp=feval(ObjFun,k,x(tmp1(i),k),u(j));7 m) I* y& z& m3 X+ I, w
    
47.          if tmp<=t_vub, 6 K9 H# ^/ j* Z' H4 X\" `
    
48.             f_opt(i,k)=tmp;d_opt(i,k)=u(j);t_vub=tmp;
    # }( Z% S7 O# i
49. end;end;end
    # |5 Z( `+ O! W. f$ i\" C( L
50. % 逆推计算各阶段的递归调用程序+ m: P/ B0 K  c6 [* l5 `+ _9 ~$ H
    
51. for ii=k-1:-1:1/ U  J3 C$ w. q: d
    
52.    tmp10=find(x_isnan(:,ii));tmp20=length(tmp10);  {+ u5 R& Z2 _\" G. n# _
    
53.    for i=1:tmp20\" \; y3 M/ C, A: L0 s
    
54.       u=feval(DecisFun,ii,x(i,ii));tmp30=length(u);# R- Y+ r4 m: i: W' d6 {
    
55.       for j=1:tmp30
    & @* M( t  L9 ?' f1 A8 w# z
56.          tmp00=feval(ObjFun,ii,x(tmp10(i),ii),u(j));
    - R- D) S5 E5 ?/ [' D\" _
57.          tmp40=feval(TransFun,ii,x(tmp10(i),ii),u(j));
    ! @( h\" T7 B  m6 P8 Y* T
58.          tmp50=x(:,ii+1)-tmp40;5 Z& E& m9 s9 c) n\" J
    
59.          tmp60=find(tmp50==0);5 a' ]# t' j4 p' ]' H, o7 H
    
60.          if ~isempty(tmp60),
    & y5 M; u8 h7 Z4 f\" C8 w# C
61.             tmp00=tmp00+f_opt(tmp60(1),ii+1);    1 f% z2 \\" |1 W( @: s
    
62.             if tmp00<=t_vubm(i,ii), x\" x, u$ G$ }4 D8 i) X( d( m
    
63.                f_opt(i,ii)=tmp00;d_opt(i,ii)=u(j);# s4 F( o9 P2 Z3 M) s1 |
    
64.                t_vubm(i,ii)=tmp00;% E$ R3 _' @0 ~6 U) f4 @6 u
    
65. end;end;end;end;end;
    ( G/ H6 D- C; C2 \0 c, z9 C
66. fval=f_opt(tmp1,1);
    , M9 y$ C; _. N, d. w: i4 h
67. % 记录最优决策、最优轨线和相应指标函数值$ q0 W& S' `. @4 }
    
68. p_opt=[];tmpx=[];tmpd=[];tmpf=[];& O, c' {- [. G8 t8 w0 J8 [
    
69. tmp0=find(x_isnan(:,1));tmp01=length(tmp0);1 r& A) ?( x+ o! F4 `
    
70. for i=1:tmp01,
    5 R9 x! `( ^% p) U0 O3 ]$ m$ s
71.   tmpd(i)=d_opt(tmp0(i),1); $ y3 {( h5 N: O\" z3 c/ `. S$ |. @
    
72.   tmpx(i)=x(tmp0(i),1);
    ( r$ j( X8 ^7 T4 }' x
73.   tmpf(i)=feval(ObjFun,1,tmpx(i),tmpd(i));
    . E  m1 q. R6 x4 L# N; o8 M
74.   p_opt(k*(i-1)+1,[1,2,3,4])=[1,tmpx(i),...
    - [* z( [\" R- o( i* N2 n+ Q. I: E( H
75. tmpd(i),tmpf(i)];5 s6 ?- V8 |( C8 X% H
    
76.   for ii=2:k' ]0 a  G& [# E1 c
    
77.      tmpx(i)=feval(TransFun,ii-1,tmpx(i),tmpd(i));
    9 e3 @. v; G& o- {
78.      tmp1=x(:,ii)-tmpx(i);tmp2=find(tmp1==0);5 W' L! c! }8 c7 O
    
79.      if ~isempty(tmp2)
    * ]5 N1 @- ~) y, K
80.         tmpd(i)=d_opt(tmp2(1),ii);0 b4 E7 J2 U. e- b) a+ {\" T
    
81.      end;
    4 L7 B8 w! c# g1 ~
82.      tmpf(i)=feval(ObjFun,ii,tmpx(i),tmpd(i));
    4 Y% Z7 T1 L+ v  R
83.      p_opt(k*(i-1)+ii,[1,2,3,4])=[ii,tmpx(i),...  D. a\" i; }  b9 ]& V& o
    
84. tmpd(i),tmpf(i)];! A) n2 C* W* Q0 |& }
    
85. end;end;
    6 U* v4 Q( ?- @. E

复制代码

danieldu

看看看看那

sam_yat

同学习下，希望有人解答。  

花开终为败

harry0326

呵呵，动态规划是一门艺术不是技术，不可能有通解的，这句话是数学大师说的。具体哪个忘了，你们有？可能吗

jmedwardliu1205