QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2680|回复: 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看懂是什么意思  求大神指点!!! e8 D: |4 u, q
    - i0 ~( f. s# U) O; W6 _
    global v
    % _% v$ w( P! O/ Q# j) _. I7 b" Ev = 630589;                % 求解v为630589的收入水平情况
    3 Y' H# \' \' C: w0 MDEBUG = 0;5 V; J6 ^4 Z; h5 N! D, Z% _
    rand('state',sum(100*clock))        % 初始化随机数生成器
    5 I; e* V: v6 ^: J2 L( K) uformat long g
    5 {# P, ~7 i' b, m; H! _1 r  d9 Y/ C/ N9 H
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 s1 X+ s) J1 I# ^- }" W
    % 求解开始
    . C6 @# L  [/ L( ~1 H& }% 对于K1型- N0 L, `* _' ?! o8 v
    p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
    $ N5 H, @2 _: ~6 ]; |7 d% 6个奖项6个变量2 V, E- g# F$ G5 M5 H( s  y
    Aeq=[1,1,1,0,0,0];
    % Y& l8 k! n( \% V- Dbeq=1;% u- F% y8 s% Y* ^: v$ z
    a_lb=[10,4,3,4,2];
    % }, \2 R4 Y) U5 S! v9 [% nb_ub=[233,54,17,20,10];
    - s2 V! Q2 c* \; {1 b# m5 ?A= [0,0,0,-1,a_lb(4),0;# G( i7 t2 j# j& O+ z/ f" k: B
        0,0,0,1,-b_ub(4),0;
    $ J- Y* ]" {5 Y" a    0,0,0,0,-1,a_lb(5);
    9 C6 e' ?" d7 b$ ^) y; {5 v2 O% j3 @    0,0,0,0,1,-b_ub(5)];" H- s6 [) }" s6 |+ W2 b
    b= [0;0;0;0];+ v/ }7 ^% J# I5 i2 F' ^' M( \# u
    lb=[0.5;0;0;0;0;0];& n' n% O2 Y3 \+ o: i# t
    ub=[0.8;1;1;inf;inf;inf];* ~( Y3 F9 k8 w6 {3 E1 [+ ^7 {- l
    p_test = p_k1;" ^# g# C- S& @: `3 l/ z( Q: X) w
    rx0_tmp = zeros(6,1);
    $ D% e  e5 ]7 b8 s+ n/ e) ^. Brx_meta_result = zeros(6,1);( D4 v: p: Y! ]0 n9 d
    fval_meta_result = inf;
    : F0 W% O7 B& q" gflag_meta_result = nan; %用以判断有没有得到过可行解( o  b. u) \: O- N# Q
    if DEBUG == 1# G9 r% e% I6 T  |+ L! d! D
        output_meta_result = [];1 R0 h7 F% {: B5 b$ t6 ^
    end
    $ y! n1 ?. W. a# pfor j = 1:nums_test_of_initial_value4 P2 U8 F# ]( Y2 B2 W9 C
        %随机生成多个初始值rx0_tmp,以避免局部最优
    , I& P* \% m8 d( D6 n: J    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    $ z5 ?$ h) r! U1 L- M    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    ) H7 q& G9 N# `5 n    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    + @' f& V5 Y. V9 O2 E    rx0_tmp(4) = rand*1000;% ^$ p+ V8 T% K$ Q) }/ G. ?
        rx0_tmp(5) = rand*100;5 _* `! T3 z" D" [8 U, G* P) n" Z$ |+ [
        rx0_tmp(6) = rand*50;
    0 w- ~! C7 i/ X3 Y  J    % 寻优0 D0 c1 K2 b+ O  |7 Q) N
        [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    7 t# `6 z) C6 O# a+ p8 s# y# I            fmincon('cpiao',rx0_tmp,A,b,...  r& X7 W+ v1 d+ m7 B
                        Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);: H9 l% y0 {4 }* N; ]
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)& R' q9 n& I! b" f
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    . c$ h' J$ L/ I& l, K" S    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)2 i6 z4 m6 X3 i( \
            fval_meta_result = fval_tmp;' V  O" |" o$ x; w0 u
            rx_meta_result = rx_tmp;& l5 g0 _. _: T' F$ ^- |2 N
            flag_meta_result = 1;
    5 v! a4 P" g7 S" B3 |0 P        if DEBUG == 1! K. ]' @$ `  C
                output_meta_result = output_tmp;
    - E0 m0 a  S/ C3 I7 H( T4 Y        end
    5 U5 Z$ Z) y9 g6 n- B8 U    end
    6 F# r3 P, Z* B0 o" f0 l  e3 U0 oend% t' R7 t; O3 U4 R
    % 把求得的最好结果保存下来
    7 e* ^+ E5 ~. `7 Q/ r5 A& F: vif ~isnan(flag_meta_result)6 O' ^( R) B" q  X, n
        rx_k1 = rx_meta_result;
    / u: f. d$ g8 N8 v    fval_k1 = fval_meta_result;
    + J1 S5 O  g3 ?4 W2 H: x! P4 M    flag_k1 = flag_meta_result;
    ' E9 T( A8 v4 O, g$ o: |    if DEBUG == 1; Z7 q" N/ i( C# b. C9 {* ~/ G* ^
            output = output_meta_result;# ?/ r/ Z9 G* Q9 S* U! F& r
        end
    7 U) ^9 O6 L" [1 j4 L+ z* Ielse2 T! e9 h, {/ r5 b, }" ^0 b
        if DEBUG == 1
    3 A8 ^( o% t% C  s7 N        rx_k1 = rx_tmp;3 Q8 R' ]* {1 A! D0 a" I; f& Z  m
            fval_k1 = fval_tmp;/ c: L$ z3 D& ]% f2 T+ {: i4 Z  L
            flag_k1 = flag_tmp;
    # _/ _8 C! ^0 x5 p        output = output_tmp;
    6 a" w, {" S0 b# k$ a% d    end, ?4 z0 f( f6 m- r
    end/ u" k0 n1 N0 z* X+ u$ m" @
    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, 2025-7-24 04:01 , Processed in 0.625348 second(s), 58 queries .

    回顶部