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

danieldu

看看看看那

madio

1. function [p_opt,fval]=dynprog(x,DecisFun,ObjFun,TransFun)       : O7 a: y7 b) Y& `% q. n
   
2. % [p_opt,fval]=dynprog(x,DecisFun,ObjFun,TransFun)
   + ?4 f# o2 ]\" J4 v: s0 l
3. % 自由始端和终端的动态规划,求指标函数最小值的逆序算法递归
   8 S: @! g5 Q) c# n0 C
4. % 计算程序。x是状态变量，一列代表一个阶段状态；M-函数  R3 ~6 f! f0 G# l0 S* h+ a$ b$ v% \
   
5. % DecisFun(k,x)由阶段k的状态变量x求出相应的允许决策变量;
   * R5 v4 T7 b1 X
6. % M-函数ObjFun(k,x,u)是阶段指标函数，M-函数TransFun(k,x,u)2 |' b( R+ k\" `9 ?
   
7. % 是状态转移函数,其中x是阶段k的某状态变量，u是相应的决策变量；
   . D! g( o+ ~# f
8. % 输出p_opt由4列构成，p_opt=[序号组;最优策略组;最优轨线组;; ^$ ^7 V' j& c8 A) K
   
9. % 指标函数值组]；fval是一个列向量，各元素分别表示p_opt各
   . M! c+ D7 k. G. u8 s5 e% B2 I
10. % 最优策略组对应始端状态x的最优函数值；
    + b3 X( K0 d( O  U, b9 ]/ t1 E( N
11. %
    \" M: @: F$ T5 d
12. %例（参看胡良剑等编《数学实验--使用MATLAB》P1804 c! \; ?3 a( z4 F
    
13. %先写3个函数
    7 A- E% S7 b: E\" a% Q6 _
14. %                 eg13f1_2.m  m; c& z& E! J8 W8 z9 `2 r1 @
    
15. %    function u=DecisF_1(k,x)1 Y! I0 W7 I, q! ]\" Q' o* Y- d
    
16. %    在阶段k由状态变量x的值求出其相应的决策变量所有的取值
    # _6 ?: F6 q* ~  q# t7 U  Q/ I5 V
17. %    c=[70,72,80,76];q=10*[6,7,12,6];& E% x; |' u- B: R% ^\" _
    
18. %    if q(k)-x<0,u=0:100;       %决策变量不能取为负值
    / g/ s\" S# J( v. U  }1 S& A1 ^
19. %    else,u=q(k)-x:100;end;     %产量满足需求且不超过100
    % `% J0 A: n2 q' l/ d  ?( N. n: H
20. %    u=u(:);
    % J! k7 }( X* S1 L/ G& F. \
21. %                 eg13f2_2.m\" @) _) D) O% R% I0 e2 H
    
22. %    function v=ObjF_1(k,x,u)9 d, |' @- W% M' f6 \) Y1 D
    
23. %    阶段k的指标函数
    ! z8 T. D' Y) L: P1 n
24. %    c=[70,72,80,76];v=c(k)*u+2*x;! Z) J; F8 |0 `/ u0 g+ n( u
    
25. %                 eg13f3_2.m
    ; g, p* u4 V' p) R$ w+ z& G
26. %    function y=TransF_1(k,x,u)- {( C* I& |% E/ |) l' K
    
27. %     状态转移方程
    ! @! x6 G3 d# r4 b; \
28. %     q=10*[6,7,12,6];y=x+u-q(k);
    ) w4 f3 V) z# b7 r; o9 b
29. %调用DynProg.m计算如下：
    ; U9 Q( [; t6 `: E- J, R
30. %    clear;x=nan*ones(14,4);% x是10的倍数，最大范围0≤x≤130,6 g3 V  p1 M  K  g6 e
    
31. %       %因此x=0,1,...13，所以x初始化取14行，nan表示无意义元素$ N% X& q$ t4 \; F, p1 g
    
32. %    x(1:7,1)=10*(0:6)';     % 按月定义x的可能取值- J) T- w2 |& S: H
    
33. %    x(1:11,2)=10*(0:10)';x(1:12,3)=10*(2:13)';
    5 q8 O- ^7 C# Q4 k% X
34. %    x(1:7,4)=10*(0:6)';
    ! n8 k% K0 O7 |* {
35. %    [p,f]=dynprog(x,'eg13f1_2','eg13f2_2','eg13f3_2')5 i1 E# t/ H# q
    
36. 
    0 [, k; Z- j( q5 I# K* \6 |' |
37. % By X.D. Ding June 2000
    3 U- N7 N& }' I. Y( v' X* a( }* l
38. ; n# d  x, s8 y: G
    
39. k=length(x(1,:));f_opt=nan*ones(size(x));d_opt=f_opt;
    # h8 r6 u3 P+ P/ K) o
40. t_vubm=inf*ones(size(x));x_isnan=~isnan(x);t_vub=inf;\" Z) e8 c9 G0 j
    
41. % 计算终端相关值
    ( x) D& q5 v+ @* \
42. tmp1=find(x_isnan(:,k));tmp2=length(tmp1);' f: T2 M: {8 a1 f+ k; X
    
43. for i=1:tmp2
    & R  _2 u$ _1 f2 |
44.    u=feval(DecisFun,k,x(i,k));tmp3=length(u);
    0 j. R7 s) q) ?, }8 K9 n
45.    for j=1:tmp3
    : _) N! @& F% q: A$ y( ~& C
46.          tmp=feval(ObjFun,k,x(tmp1(i),k),u(j));
    1 f6 S1 u3 U2 e  \
47.          if tmp<=t_vub,
    , `6 b, M0 V2 ?, R
48.             f_opt(i,k)=tmp;d_opt(i,k)=u(j);t_vub=tmp; 8 i) Z  w3 A1 Q/ H$ M
    
49. end;end;end. r5 }% p8 [; ^/ `% s& B
    
50. % 逆推计算各阶段的递归调用程序- V% x7 P$ F6 @0 Z' S8 E  c4 ?
    
51. for ii=k-1:-1:1
    2 K6 N  q( t, f
52.    tmp10=find(x_isnan(:,ii));tmp20=length(tmp10);
    5 Y1 d3 [! x5 w5 l  y) x3 Y
53.    for i=1:tmp209 P  B* X# Y0 i. O' {4 n
    
54.       u=feval(DecisFun,ii,x(i,ii));tmp30=length(u);
    * G; D. s2 ?! N) H
55.       for j=1:tmp30  @# g( N- p& n3 f& ?
    
56.          tmp00=feval(ObjFun,ii,x(tmp10(i),ii),u(j));
    # g\" p' x( J+ \0 E* u
57.          tmp40=feval(TransFun,ii,x(tmp10(i),ii),u(j));$ P$ K2 o( z; |3 {6 J
    
58.          tmp50=x(:,ii+1)-tmp40;
      a  i! ~9 L' g0 _+ ?\" P
59.          tmp60=find(tmp50==0);
    ( N% V/ P3 l  V) ?  ^
60.          if ~isempty(tmp60),* h& H9 s  s; F' p
    
61.             tmp00=tmp00+f_opt(tmp60(1),ii+1);    ! `- j' d6 f4 Q; R+ h4 [
    
62.             if tmp00<=t_vubm(i,ii)! [\" s, [0 W+ H' p. O; K
    
63.                f_opt(i,ii)=tmp00;d_opt(i,ii)=u(j);
    \" p# `! {2 O  Q& r
64.                t_vubm(i,ii)=tmp00;: a0 ]0 i) K( _, M  B
    
65. end;end;end;end;end;\" ~( Q) `! H) @2 Z+ r6 ?
    
66. fval=f_opt(tmp1,1);) G1 X- u3 w/ p* H2 K$ I5 g
    
67. % 记录最优决策、最优轨线和相应指标函数值8 M% F. g0 t: }4 Y  O( C
    
68. p_opt=[];tmpx=[];tmpd=[];tmpf=[];+ m2 N2 Q( H6 Q- u
    
69. tmp0=find(x_isnan(:,1));tmp01=length(tmp0);. ~4 R. i  @% w; |' f$ U1 Q% q* u5 Y
    
70. for i=1:tmp01,, x4 [- v, s) n8 H4 q/ k
    
71.   tmpd(i)=d_opt(tmp0(i),1);
    ) g9 M+ G! g\" d! K; |) |; `: b
72.   tmpx(i)=x(tmp0(i),1);% U: }\" a; B\" t
    
73.   tmpf(i)=feval(ObjFun,1,tmpx(i),tmpd(i));/ W\" V/ o$ B! J. ?$ M* U7 Y\" L4 ?
    
74.   p_opt(k*(i-1)+1,[1,2,3,4])=[1,tmpx(i),...8 W2 l! d) e* [/ K5 [  A
    
75. tmpd(i),tmpf(i)];, o/ Q% e, Z' q! K, Z
    
76.   for ii=2:k
    & ?4 C4 ?  ^0 c7 T3 ?  e, d
77.      tmpx(i)=feval(TransFun,ii-1,tmpx(i),tmpd(i));
      ?, t3 K) O: {0 o9 n$ H' Z
78.      tmp1=x(:,ii)-tmpx(i);tmp2=find(tmp1==0);& j$ t9 p% R. _$ D. b0 _( D
    
79.      if ~isempty(tmp2)# Y% ^3 @# B* `\" H& @3 v2 O$ e\" _& E
    
80.         tmpd(i)=d_opt(tmp2(1),ii);/ |! j' l% v  |( z
    
81.      end;
    7 U. R3 O' o. t
82.      tmpf(i)=feval(ObjFun,ii,tmpx(i),tmpd(i));9 {: V4 Y. K9 T  K* ^
    
83.      p_opt(k*(i-1)+ii,[1,2,3,4])=[ii,tmpx(i),...+ d, e& r* M& P' v( m
    
84. tmpd(i),tmpf(i)];
    ) ~) g' g$ K2 L6 ^0 d
85. end;end;. P3 o$ S+ _/ t% d. o5 ?- I, A
    

复制代码

m存在

张坤14

hugoczx 发表于 2009-12-21 23:20
有是有                              
, C9 ]: w) D7 w  ?不过不能白给！

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

无语的昭昭

额额额，其他地方很多啊
3 l7 J5 G4 r  Z) x9 t