QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3028|回复: 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看懂是什么意思  求大神指点!!
    4 J# {  ?/ T8 e/ e; }. K3 x6 D; d5 K  _- a. k- k
    global v
    # v/ w; {& ]1 q; H! `v = 630589;                % 求解v为630589的收入水平情况7 i( s3 c; G' p1 u  I
    DEBUG = 0;9 f' x9 d7 C: D) p' W% `, ~
    rand('state',sum(100*clock))        % 初始化随机数生成器
    - D7 B3 ]. p! J6 Bformat long g2 m$ O( c* c# w6 E

    6 K- L3 b  j! i) {+ O6 G3 o& l7 y%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2 j3 y: Y7 b) v3 ~- B. e% ^% q3 E
    % 求解开始, Q+ P5 i! }+ L+ Z" |  H( {& W
    % 对于K1型
    1 r  X/ |! s0 z# a6 Y& Mp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];( b+ m& w9 _% X$ V. A
    % 6个奖项6个变量3 b0 e' {5 F" F1 V9 I# u- a
    Aeq=[1,1,1,0,0,0];8 `5 K# c: S: g' n+ k  F  m
    beq=1;) h! O2 N) Q" z8 k; \
    a_lb=[10,4,3,4,2];1 K, j$ L$ x, T% @5 I
    b_ub=[233,54,17,20,10];$ q* `# _1 V& x& v; F, w7 g, N
    A= [0,0,0,-1,a_lb(4),0;$ h' i, O3 c7 ?, S& ~
        0,0,0,1,-b_ub(4),0;' A! g- O! g8 P* R' Z  J
        0,0,0,0,-1,a_lb(5);
    9 M2 z: W! f, v* Y( k    0,0,0,0,1,-b_ub(5)];+ y( s* d) P% P: x5 V
    b= [0;0;0;0];" u' N. j8 j5 J4 R; A0 d
    lb=[0.5;0;0;0;0;0];) \- E2 a* R: G6 `& X# ]
    ub=[0.8;1;1;inf;inf;inf];* c& f- {2 U/ x0 @5 s2 B
    p_test = p_k1;' Y+ y* @$ ?$ Z0 S- n: f5 d
    rx0_tmp = zeros(6,1);4 s9 Z* ~$ [: Q7 d3 U2 t- p
    rx_meta_result = zeros(6,1);5 h' u' u) K( F9 K
    fval_meta_result = inf;/ t5 {3 U. j( }) |5 S
    flag_meta_result = nan; %用以判断有没有得到过可行解  |6 F) w- D# k  z4 i
    if DEBUG == 1: u, A% u  R, C* q
        output_meta_result = [];
    ' A1 Z) y! I' I2 j8 u3 aend2 O- I! R& v. q# z- c
    for j = 1:nums_test_of_initial_value
    5 C4 Q4 Y$ F! |    %随机生成多个初始值rx0_tmp,以避免局部最优
    4 m7 O6 u8 M# V" n    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    % k$ m2 t; o8 w# X    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    # u8 }: M1 b0 j6 [+ Q    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);1 E# a$ l/ S4 p; j6 t
        rx0_tmp(4) = rand*1000;
    7 t" a- v/ a: @, z- W* m    rx0_tmp(5) = rand*100;7 _  y! A( U; F; K
        rx0_tmp(6) = rand*50;
    2 r4 o9 ?8 p0 `8 y$ u    % 寻优* @! o; d0 \) S2 B. {* t
        [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    8 u% Z7 [3 \* i, g* O            fmincon('cpiao',rx0_tmp,A,b,...
    + J4 l7 b* o6 K  B                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);# f/ x5 ]% t# t$ N$ T7 z  l# N" b5 n
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)4 R) X* }- a, Z; j9 x  O
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。% Q0 a" y- }2 `& H% ]' O
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)- r6 O4 y. C" X9 W7 H7 }0 o. n
            fval_meta_result = fval_tmp;' y0 {2 w* a3 u; y. J; U
            rx_meta_result = rx_tmp;! R( b: k: Z( o" }$ _
            flag_meta_result = 1;3 L" `6 O* F5 Y( N
            if DEBUG == 1
    . p* T" B  X5 x- b2 @' Z( w            output_meta_result = output_tmp;7 w$ j  o. M7 M
            end1 h7 R. o) S! t6 J( `: T3 m' j4 u* B
        end
    % _+ W- G9 s! Y# y6 M9 Y  _) L2 pend
    % `) V7 ~8 M# d0 h0 y! O) M% 把求得的最好结果保存下来
    % T' i# l7 s# z. B+ r' zif ~isnan(flag_meta_result)
    3 P7 ~5 A7 @* }    rx_k1 = rx_meta_result;
    % S, c: W& E; [) s9 {- ]& S$ @/ F% D1 V    fval_k1 = fval_meta_result;
    ( B& w2 W4 W3 u* \3 c    flag_k1 = flag_meta_result;3 Z3 l$ j! y9 z$ d) K' D5 q- G1 g
        if DEBUG == 18 I; `8 U8 M0 [/ v; F/ k8 O
            output = output_meta_result;
    * m5 S: |  F( R3 A6 u: R1 ^; f    end9 H9 y2 q: g: G; t& p
    else/ ?4 U% E, W, K* C
        if DEBUG == 12 e; F. w' r. y! \4 x. }7 Z3 s0 O4 L1 l
            rx_k1 = rx_tmp;, x  |& G/ Q3 f6 K" p
            fval_k1 = fval_tmp;
    % N$ [9 Y! q8 n* l$ u        flag_k1 = flag_tmp;
    3 K1 W7 V( X) ~& |        output = output_tmp;, y: x0 c7 x4 n! D" `5 A
        end
    # Q8 D- o( z0 m* L( zend
    5 A9 z+ ?# Y, E3 b+ ?0 Y+ ~4 ?
    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 15:05 , Processed in 0.323756 second(s), 59 queries .

    回顶部