QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3033|回复: 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看懂是什么意思  求大神指点!!
    8 D0 h% l6 {5 e* {1 p( w
    : Q; D) I2 ^0 Q' g! h7 ~! Eglobal v
    9 s- s  O+ W( s( i  o+ z6 H0 b$ u$ dv = 630589;                % 求解v为630589的收入水平情况& e- r/ p# @( l, t1 ^% h
    DEBUG = 0;$ U/ L% `. V& o( N1 e$ G7 o
    rand('state',sum(100*clock))        % 初始化随机数生成器
    ! E) y1 G/ X/ |# a# zformat long g
    : U" O5 b( T" x7 J7 @3 y) I  Q
    ! y! Z8 ?5 l" M7 k1 }7 I- A%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%- B  \5 z5 u5 x2 S
    % 求解开始
    # b5 b( d# B) V6 H" ~% 对于K1型
    ( l* H- Q& E! @! ?; ]7 ~p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];( A9 u/ r* Y, b" F+ d4 v; o
    % 6个奖项6个变量5 R  S7 m1 p- h9 b
    Aeq=[1,1,1,0,0,0];$ Z6 B  m4 E$ K3 c' O( s
    beq=1;
      v2 v4 z. d: U1 D, [a_lb=[10,4,3,4,2];% e) ?0 r/ e* |  s% K3 s
    b_ub=[233,54,17,20,10];2 F. P7 k# L5 Z+ |3 i- Y
    A= [0,0,0,-1,a_lb(4),0;: |1 F# U/ u  [
        0,0,0,1,-b_ub(4),0;2 q, J0 U! B3 h6 ~/ W& g& \" C
        0,0,0,0,-1,a_lb(5);6 q$ T/ n; d" a! n: M5 o
        0,0,0,0,1,-b_ub(5)];
    . C" x( x1 l$ f, T0 |b= [0;0;0;0];2 ?% I: w- I2 n' B7 z
    lb=[0.5;0;0;0;0;0];$ Z# W% b* J$ Q+ k6 m; x- x
    ub=[0.8;1;1;inf;inf;inf];
    % R- c- k- Z3 m4 H  o3 B! c2 W5 L; r7 ip_test = p_k1;8 i; H) G- x& G9 W; R8 y" v
    rx0_tmp = zeros(6,1);
    2 @. q& B' _5 L% H- Rrx_meta_result = zeros(6,1);. O" A9 L$ ~+ V- d0 V
    fval_meta_result = inf;
    5 P6 p) ~3 Y1 Y4 f* Iflag_meta_result = nan; %用以判断有没有得到过可行解5 v* W, I: x% m' Z
    if DEBUG == 1
    1 X: @5 @5 l' c# U) f    output_meta_result = [];
    9 I( \; h8 q4 @0 a) V# Y# `end1 |  h/ w% ?7 o( ^: h; n# {' r9 M+ b# c" Y
    for j = 1:nums_test_of_initial_value7 d/ g( Z) d4 `- h
        %随机生成多个初始值rx0_tmp,以避免局部最优( M. ]5 G# S& C$ @/ ~
        rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    4 R; N; _- M  H* b2 m7 W    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    4 C" i/ d, p" h5 O5 K- i, K    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);! i2 v& g2 r( o* d
        rx0_tmp(4) = rand*1000;
    $ y# M+ X/ s. t    rx0_tmp(5) = rand*100;7 a! |! c4 a) f' E* u$ N
        rx0_tmp(6) = rand*50;8 h  }8 c/ }7 D& \2 i( O( G" T+ k) K
        % 寻优
    $ ?- b5 `% I0 Q" \    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ..., \* I4 R8 P0 R
                fmincon('cpiao',rx0_tmp,A,b,...
    % E5 n* C8 p" a                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);9 {5 r3 }4 S' g& d
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)1 w  K. K7 q( Z
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    ! G9 p$ [( G; x$ }. b    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
    6 }% I- t4 Q) s- e" ^' r% I        fval_meta_result = fval_tmp;" \+ z3 c, U% G: E
            rx_meta_result = rx_tmp;  }; I6 a3 g6 v( i) R7 E0 h7 ?
            flag_meta_result = 1;
    3 V0 e& _' z7 x7 Z        if DEBUG == 1/ W8 d) _, m' Z, N; \6 {
                output_meta_result = output_tmp;
    / J( M$ L7 F7 W        end
    , `% ]- {7 }8 ]" [8 j* ^) f8 J6 h5 x    end' j  ^% j% U0 ^" K6 Y
    end. R2 y1 [6 _/ r' V" q8 h$ v! l
    % 把求得的最好结果保存下来5 I9 q+ c7 K: U
    if ~isnan(flag_meta_result)
    ) S  o9 P3 o8 m; b! _. b    rx_k1 = rx_meta_result;! `# y6 i& F# S
        fval_k1 = fval_meta_result;: A! C5 C  \+ o, `. p$ u
        flag_k1 = flag_meta_result;
    2 C7 `6 C) Z8 q    if DEBUG == 1
    9 p6 ]' P! g- I, k        output = output_meta_result;
    1 i: F) k# }. v8 w4 W    end
    ; m" b$ p% i4 Z+ F; S1 F4 {else: C. G! v7 J, \% Q0 @& [$ z. Q; a
        if DEBUG == 1
    0 `; A3 @0 s( F        rx_k1 = rx_tmp;
    8 {: L( Q" t, ?- v! c; i$ @" p        fval_k1 = fval_tmp;- Z% C1 j  i$ F/ k
            flag_k1 = flag_tmp;
    4 S: Y4 K+ x0 T8 l# O' X! C        output = output_tmp;
    5 o6 H( W1 O- ^( w    end. O* ?6 e% A% O, J: }: ~
    end
    9 O" B/ m# G( O; a" C; {
    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-8-6 22:48 , Processed in 0.559975 second(s), 57 queries .

    回顶部