QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2682|回复: 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看懂是什么意思  求大神指点!!
    3 w7 y: J3 w+ E0 a: c" [+ E; o% t# q1 _# t. H
    global v1 d- W1 U+ h9 V  M5 K9 ?/ p4 a3 r) i
    v = 630589;                % 求解v为630589的收入水平情况
    9 o9 P8 o2 r  G% C2 PDEBUG = 0;5 M$ v1 O* j; t! S( U
    rand('state',sum(100*clock))        % 初始化随机数生成器2 \( ^7 q* h' \
    format long g! i9 u4 }& t1 X$ P  E* q  _

    4 h5 `2 y3 }( R6 ^%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    % n! y! [' o' R; V, R: Z% 求解开始$ P3 E! s/ R6 B9 G
    % 对于K1型) O$ o7 v3 e' `# K" m
    p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];2 j5 t) H+ i/ R& ^
    % 6个奖项6个变量4 g7 O2 {0 f5 m8 [( q% I, e# X
    Aeq=[1,1,1,0,0,0];
    ( P0 x# O  B- ?5 W0 Nbeq=1;7 ^2 V1 e  o% E; k) i" I
    a_lb=[10,4,3,4,2];
    4 Q8 G1 }# F4 d% B0 n: Rb_ub=[233,54,17,20,10];% C- ?) _3 e/ l0 q  W0 Q* _
    A= [0,0,0,-1,a_lb(4),0;4 x/ Q( C# J! ?# C( T
        0,0,0,1,-b_ub(4),0;: w" V* e* G5 i% r) G% d
        0,0,0,0,-1,a_lb(5);
    , I" o1 F4 l1 r$ n/ U    0,0,0,0,1,-b_ub(5)];
    " ?" p- @; D" ]. e; I: ub= [0;0;0;0];% y9 h" E) L: R, C1 {  H5 l+ e
    lb=[0.5;0;0;0;0;0];! q: O5 G0 ?/ f: h
    ub=[0.8;1;1;inf;inf;inf];) N6 _2 [3 t, y. n9 h
    p_test = p_k1;
    ) m  C+ _1 ^% a# Prx0_tmp = zeros(6,1);
    ! M: @6 s: K$ |$ T7 a% Erx_meta_result = zeros(6,1);
    $ h( S2 T- w$ V/ O0 u4 l* E% cfval_meta_result = inf;) b0 H+ y' H' m3 Z4 P' u
    flag_meta_result = nan; %用以判断有没有得到过可行解' e9 ~3 _) \% h: k. {* ?
    if DEBUG == 1
    7 g3 w7 m/ c, x" f8 i* q. a2 i    output_meta_result = [];: r/ b! P/ A) A
    end4 ?, n( }* n0 v$ {- z
    for j = 1:nums_test_of_initial_value9 e8 M8 F. U5 H, F, x
        %随机生成多个初始值rx0_tmp,以避免局部最优: B9 O- p& D3 V% \
        rx0_tmp(1) = rand*(0.8-0.5) + 0.5;6 l2 h* n$ H( q  c4 k
        rx0_tmp(2) = rand*(1-rx0_tmp(1));$ W: a8 m! s. J2 d# G
        rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    ) a$ G! {  y4 c, P6 O+ w- Q+ `    rx0_tmp(4) = rand*1000;
    : V6 N3 O8 ?6 t! x5 C    rx0_tmp(5) = rand*100;; g1 @" F2 w7 i  W0 ]
        rx0_tmp(6) = rand*50;6 `4 Y- Z# }& p" }0 y
        % 寻优
    6 @5 G6 N7 ?7 Q8 ~    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...) V" x, o1 D% C+ ?; r: N% {
                fmincon('cpiao',rx0_tmp,A,b,...
    $ {% n2 G9 U' N- w9 ^$ w$ r' u( a                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
    & P# q1 Z! f0 \5 q: ]1 n+ x        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)% g5 {8 D* Z. c( ?9 Q& {: e7 p, Z
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。7 b1 J  L  B. i) a1 F
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)( i6 z& N, c6 x7 i9 K) ?. ^
            fval_meta_result = fval_tmp;0 P  B5 o5 p- _/ p
            rx_meta_result = rx_tmp;
    + A1 R4 [& i* X- y. w. `        flag_meta_result = 1;" }. [% Z5 B9 y  z% o* m" L
            if DEBUG == 14 U  h1 {* E6 L2 G1 p7 i& A' j
                output_meta_result = output_tmp;- N- }6 p; o# t' K8 d1 W8 w+ S; R5 G
            end
    4 J# `+ `! ^' P6 p    end1 ~; c; [8 E8 D% a
    end
    ( U# ?( f- [5 H* U% 把求得的最好结果保存下来2 B/ ^' j3 r3 i: W
    if ~isnan(flag_meta_result)
    ( v1 e3 v2 R; Y+ X! ]7 c    rx_k1 = rx_meta_result;
    9 B  p" u6 M6 f    fval_k1 = fval_meta_result;) D% m  ]- U6 T. k
        flag_k1 = flag_meta_result;
    ' U2 e) g5 t# {9 d    if DEBUG == 1) W! I: V3 C# o, `" W
            output = output_meta_result;% @+ I& d; k. |9 g6 a
        end
    ' r+ a( \4 \/ |" D  b% U, Oelse
    * l3 u4 ], e4 i    if DEBUG == 1" v1 @& }8 G! Z$ ?4 H" Q% P3 S, c! w  N  C
            rx_k1 = rx_tmp;
    6 T. U% G# ?9 u4 i  k        fval_k1 = fval_tmp;
    ; S5 U9 u1 k% V: v8 z+ Y) h        flag_k1 = flag_tmp;
    1 B; A' @0 R; M7 L8 T        output = output_tmp;6 q/ f! ?1 C* N# r: W
        end, [; u' K0 o+ ?  Z
    end
    9 V" g9 d, G  M1 j; g/ 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 22:18 , Processed in 0.424472 second(s), 58 queries .

    回顶部