QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3081|回复: 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看懂是什么意思  求大神指点!!/ k* ?/ c( n* V8 k/ l- b. r6 o
    1 Y+ Y& g1 M- u8 V0 v9 P& F( ]
    global v
    5 c, D3 Y2 m$ ~9 Zv = 630589;                % 求解v为630589的收入水平情况0 Y4 t9 ~# c% J$ Z5 p
    DEBUG = 0;
    - G& P/ [3 k1 U! |* I0 T. G0 trand('state',sum(100*clock))        % 初始化随机数生成器
    + P" o3 t3 D6 ?. v2 {7 P% Uformat long g
    , [3 P1 D# n3 g# K5 E6 a0 Q# }/ y
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%9 `7 [( t- z1 `) P( Y. S! ~
    % 求解开始  K% j: I3 b# J; j" w' v+ i
    % 对于K1型
    / F0 N) O% O% o* \0 H* g6 Bp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];# c8 V6 t6 n# e, |+ @# `( a$ U
    % 6个奖项6个变量! z0 U- [+ ?1 d& N- z7 I$ E
    Aeq=[1,1,1,0,0,0];
    % f: W$ y/ d. d% C. u% k* _% jbeq=1;3 r  R8 ?! ]' F2 M, g
    a_lb=[10,4,3,4,2];
    " R2 x. R- L/ H. F/ H6 I# hb_ub=[233,54,17,20,10];8 b( A5 a' e1 |0 W% J/ N2 u( n6 ?, {
    A= [0,0,0,-1,a_lb(4),0;
    - x" `- ?0 \3 L2 s6 M    0,0,0,1,-b_ub(4),0;
    . x3 P7 G9 R, O7 S: u5 L7 x    0,0,0,0,-1,a_lb(5);
    4 T7 W9 I" a( c) R2 |) d    0,0,0,0,1,-b_ub(5)];
    / d& x/ i% d1 K. mb= [0;0;0;0];
    4 I0 n& N+ Z, s# ~( o. A* olb=[0.5;0;0;0;0;0];
    + k6 N/ v: }# ^ub=[0.8;1;1;inf;inf;inf];* ]6 E$ j) u' g
    p_test = p_k1;
    : i7 b" O# _1 C) Q! S5 u% y, g" ]) rrx0_tmp = zeros(6,1);
    ; R/ A* b$ X( y' t( [rx_meta_result = zeros(6,1);
    * E$ L9 W) q4 Z/ z. ~fval_meta_result = inf;) ~8 Y4 _( o3 [5 {
    flag_meta_result = nan; %用以判断有没有得到过可行解) I' U4 k- O# F  P) F8 [
    if DEBUG == 1# G" A3 h* P) g/ C- {2 A. |: _
        output_meta_result = [];% U) d  W* p. F
    end2 J3 o0 X6 N/ B9 P' D6 I
    for j = 1:nums_test_of_initial_value  K. T1 I9 k/ w- D1 @- n
        %随机生成多个初始值rx0_tmp,以避免局部最优
    3 S# I; _* Y  H+ R& E    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;3 D/ t3 C7 E' s+ _! ?! K" ]5 }
        rx0_tmp(2) = rand*(1-rx0_tmp(1));
    ; f$ u: ?( P7 N1 _0 b, a7 w. B    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    ) P# ~: \9 g4 d) A0 ~# x    rx0_tmp(4) = rand*1000;
    ) M5 A, H, W0 J+ j+ F' e+ f    rx0_tmp(5) = rand*100;- j3 E( T: l  }# E; F. w0 w
        rx0_tmp(6) = rand*50;
    # H- V/ q  E5 O* Z/ I    % 寻优
    ' J& O' @# v0 |$ W6 ]4 |' e9 H    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    3 n- T; L- \: o* p            fmincon('cpiao',rx0_tmp,A,b,...! C6 m% ]# k9 A9 O9 E. p, N8 V
                        Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
    " x7 H/ d7 M+ m2 w' W* a        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
    : z2 U5 L  h, z6 K5 T5 `        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    / ?9 W/ _' m% v9 y8 P    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)0 f& _; ^3 X! c  \: L' e
            fval_meta_result = fval_tmp;
    7 ^" U$ u: \8 h; V  E        rx_meta_result = rx_tmp;
    ( N: u2 v3 D5 R! N, G        flag_meta_result = 1;6 J5 {+ v0 m9 J2 ?5 F0 v8 k
            if DEBUG == 1
    % ]: m$ B# E, H8 g0 Q+ l5 ^) M# \2 _. o            output_meta_result = output_tmp;
    " d, v$ ]& |3 y; M8 `. k        end( f2 J8 Z1 V) \) l; r
        end6 z+ B! t+ c2 j+ N
    end
    " S" B% ]& A/ Z0 ?" C1 w7 @! ]% 把求得的最好结果保存下来
    ) p# p9 s) _2 k$ Eif ~isnan(flag_meta_result)7 C; ^+ C  L% N. u9 R' ~, ]6 y
        rx_k1 = rx_meta_result;
    + a9 a% P+ R* s( x, k1 a9 f) H! y3 @    fval_k1 = fval_meta_result;
    + \  o% \1 w- @5 [    flag_k1 = flag_meta_result;
    ( i) U5 y" _& H2 [, Z) x& E    if DEBUG == 13 |  H; ^" M3 R2 F8 k' |
            output = output_meta_result;
    : H& C& h# M; L* F    end
    7 K* k, P' m" a9 uelse# i' U. ?! e7 R% }
        if DEBUG == 1; Q7 Q! i0 H0 u: i" }& g% G
            rx_k1 = rx_tmp;
    7 N, q" t$ C1 p8 v6 @9 d, F        fval_k1 = fval_tmp;+ h/ Q" R( {6 H% O
            flag_k1 = flag_tmp;% T! {. u9 ?, s1 D  h  r! |
            output = output_tmp;5 W4 k. l) z2 W0 c$ r( m' e3 Z1 F* {
        end( E' C4 T* R  r0 y4 Q
    end1 h, d2 e, L) R6 l4 a
    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-9-22 22:16 , Processed in 0.367979 second(s), 58 queries .

    回顶部