QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2993|回复: 1
打印 上一主题 下一主题

[问题求助] 关于debug的问题 。。。。

[复制链接]
字体大小: 正常 放大
Bojay        

1

主题

3

听众

17

积分

升级  12.63%

  • TA的每日心情
    开心
    2012-5-16 15:39
  • 签到天数: 1 天

    [LV.1]初来乍到

    群组MCM优秀论文解析专题

    群组数学建模培训课堂2

    群组数学建模培训课堂1

    群组2011年第一期数学建模

    群组第二届数模基础实训

    跳转到指定楼层
    1#
    发表于 2012-7-22 13:27 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    这是02年CUMCU 彩票问题  debug看懂是什么意思  求大神指点!!1 e/ \* @+ g5 J* E, j  o* q3 Z8 h# J
    ( f4 s# z  r3 ?. a& ?
    global v( b1 B6 {6 S8 ?: ~8 I4 R8 T3 S
    v = 630589;                % 求解v为630589的收入水平情况. [2 B5 R8 J! _
    DEBUG = 0;0 a& M# L3 O9 ~) O1 h7 i# w
    rand('state',sum(100*clock))        % 初始化随机数生成器' k% p/ P  P9 \  I2 U+ T' ?% k! f7 U* A+ a
    format long g
    : C# s* @# y  x# Y
    ; f2 q+ |1 d0 \, W. k+ ^%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    5 i, u  B8 D3 H: b+ z2 w! I, i% 求解开始
    / M3 D1 @  K$ x7 M4 l% 对于K1型. J5 L+ g9 D  \; t
    p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];/ d4 ^, N3 D8 J/ a! j( s* `
    % 6个奖项6个变量9 M) D: G" g7 R$ A6 d
    Aeq=[1,1,1,0,0,0];& l9 f, t+ m) d& w7 f+ V  }
    beq=1;
    & ]) v- L4 s% P1 Qa_lb=[10,4,3,4,2];6 ], ^2 y7 ?* a1 g" q5 }, h- ?4 U
    b_ub=[233,54,17,20,10];
    9 ~) E2 @, u* O0 ^% I" dA= [0,0,0,-1,a_lb(4),0;
    8 ^% [' w+ o4 p, V* y    0,0,0,1,-b_ub(4),0;
    ! P. t" L+ x' v    0,0,0,0,-1,a_lb(5);% m6 t5 T% A% G0 H2 M$ L, s4 H
        0,0,0,0,1,-b_ub(5)];+ Y( M# q0 i" o( S# t& L
    b= [0;0;0;0];) n9 R! n( v, q- x& h
    lb=[0.5;0;0;0;0;0];9 u, u: A  U6 ?& i  h3 H4 {# s
    ub=[0.8;1;1;inf;inf;inf];
    9 q8 l. w3 q, D2 I+ Sp_test = p_k1;
    . O- k5 f& \2 z' Z6 \5 M+ a+ S# K3 Z0 crx0_tmp = zeros(6,1);! S( L# ?' [" w, t7 `( Z( J
    rx_meta_result = zeros(6,1);  e& Y# t1 a5 a$ m7 Y* X7 `
    fval_meta_result = inf;) \3 m3 t- F( r. c4 d* I+ C8 D
    flag_meta_result = nan; %用以判断有没有得到过可行解
    9 c, L1 A' b1 a  gif DEBUG == 1
    ' f4 N# y8 q' i2 Z% ^# J    output_meta_result = [];9 N* i* q' G9 y: \  {2 A* |/ e
    end
    ! m; A% L$ A1 R. I* }for j = 1:nums_test_of_initial_value) e0 Z! z3 |1 K
        %随机生成多个初始值rx0_tmp,以避免局部最优
    # r1 W# K) `5 G( J    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    ) Z9 ]: `# q4 b5 Y6 y    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    ; u) E8 S- y# l  N. {: o    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    5 |8 T6 ~/ ]& P$ w% |6 ^    rx0_tmp(4) = rand*1000;
    , a  D0 a4 ~; d; N    rx0_tmp(5) = rand*100;
    + E. K" L- d3 R* z: ]    rx0_tmp(6) = rand*50;
    + R) m6 a  G+ K# _    % 寻优
    $ Q" ]/ q" _9 B% j$ w    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    9 ]2 F( Z9 T1 q9 K& Q0 X            fmincon('cpiao',rx0_tmp,A,b,...
    ( a/ ]) O5 k6 K7 |' J                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);! p! o' G3 \) u- `$ _9 _
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
    ; t9 o5 v; b7 _8 Y5 M, P        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
      L9 n, Z) d, c* c/ _3 k# B$ D    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
    " ], }- O4 F( ^9 V! @. T        fval_meta_result = fval_tmp;
    / w4 N2 G. ~. N% ^3 u9 H        rx_meta_result = rx_tmp;
    ( G4 e& P) c, ?; M- B- ]        flag_meta_result = 1;6 C9 B, v+ B  O4 \
            if DEBUG == 1
    7 I3 G1 C9 ]" E9 q* U; t8 {            output_meta_result = output_tmp;
    $ b' S# n; Q& g, L) {. Q( C        end# D6 Q% c( h  }) V* v; ^
        end
    / W% V+ n5 T9 E1 N$ Uend
    - W( d2 L7 G) a1 g6 E( o+ p* C& X% 把求得的最好结果保存下来" e* I) w6 Z% l; Q6 K
    if ~isnan(flag_meta_result)2 O0 ~+ \$ x0 U! x2 U
        rx_k1 = rx_meta_result;
    ' h4 b1 t* v; Z7 D* K    fval_k1 = fval_meta_result;
    : l- d, u! m& q, x/ k    flag_k1 = flag_meta_result;8 O. v  k/ Y5 o/ F
        if DEBUG == 1
    ; B2 ^1 N6 x0 J% {        output = output_meta_result;% c# A. B: c: `1 Q: c/ K: Q* y8 |. q
        end  e! c- ]: x) F4 X3 w8 ~& o7 h
    else/ S6 p. ^% ^3 M
        if DEBUG == 19 o/ ~  J+ F. y1 ^9 @; u
            rx_k1 = rx_tmp;% O5 l. v% h% l3 Q+ T* l  X+ b% w
            fval_k1 = fval_tmp;
      ~) U  {8 I7 k$ a        flag_k1 = flag_tmp;
    / }' z$ C' G8 K6 b0 |2 h3 P* ^% j        output = output_tmp;# u+ L$ g, O; Y8 v, t+ S. M0 C
        end
    ; c# x2 m3 ]7 G# h) J+ zend
    7 D- \" x! q* X2 S& y+ K
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    我那崇高的理想 有时卑鄙的幻想
    Bojay        

    1

    主题

    3

    听众

    17

    积分

    升级  12.63%

  • TA的每日心情
    开心
    2012-5-16 15:39
  • 签到天数: 1 天

    [LV.1]初来乍到

    群组MCM优秀论文解析专题

    群组数学建模培训课堂2

    群组数学建模培训课堂1

    群组2011年第一期数学建模

    群组第二届数模基础实训

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-11 02:01 , Processed in 0.393580 second(s), 59 queries .

    回顶部