QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2974|回复: 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 `0 x' X- M  _* P, w  ]3 X* R+ M4 p
    $ Q3 B$ N3 ], I  aglobal v5 x+ _2 ~* [! E0 J
    v = 630589;                % 求解v为630589的收入水平情况2 ^8 ]% }* `# Y2 m( F) x
    DEBUG = 0;# |) B+ L5 \2 W- W) N7 F3 D
    rand('state',sum(100*clock))        % 初始化随机数生成器
    ! A) d7 S; ]$ [( Zformat long g7 Y) @; f. d8 [* |
    ' y1 T* _$ k$ k
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%& f9 k5 \5 F) D: J  k+ g5 B  P# X
    % 求解开始
    5 a0 C8 M( M! U% 对于K1型+ x3 W$ R+ l* V6 V! w
    p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
    6 [& X  B- `3 j! ?% 6个奖项6个变量
    8 r7 M' O$ R8 O0 ^, ]. sAeq=[1,1,1,0,0,0];" X* J( I2 O4 I% i8 r; N2 l, J7 E! Y
    beq=1;
    - A: B+ C, p5 B1 r1 I! fa_lb=[10,4,3,4,2];+ B" K& o1 l- |0 P- h! y# u0 [
    b_ub=[233,54,17,20,10];
    $ |7 z8 H! V0 E. x% ZA= [0,0,0,-1,a_lb(4),0;
    % k! R! _  t/ `    0,0,0,1,-b_ub(4),0;0 a" U: A- d' Q
        0,0,0,0,-1,a_lb(5);
    ( O+ Z+ H- ~& W5 R    0,0,0,0,1,-b_ub(5)];
    + Z8 m+ G6 P3 d8 q1 f4 {1 X7 ib= [0;0;0;0];
    5 u. Q: F8 Q0 E2 i* l# {lb=[0.5;0;0;0;0;0];
    0 J6 R  O2 x9 v; H' W( ^ub=[0.8;1;1;inf;inf;inf];, q& S) V3 `, ~! _
    p_test = p_k1;
    % A) T% ?* r3 Crx0_tmp = zeros(6,1);; [8 r7 D/ k# `9 x0 G/ f6 b9 X* d6 w
    rx_meta_result = zeros(6,1);6 M  E% Y2 ^$ @& x  o; c) q# ?/ `
    fval_meta_result = inf;) O! N* ~* u) _/ X8 ~$ A- R' N( S
    flag_meta_result = nan; %用以判断有没有得到过可行解
    3 q. [. j+ n7 N9 ~) i9 p+ Pif DEBUG == 1
    5 q1 D8 s) V7 H+ f# W    output_meta_result = [];% k( K6 x4 `- U4 q" q, N6 G
    end
    8 {: M* d0 g" o6 Hfor j = 1:nums_test_of_initial_value
    ; {3 f7 F# S, D7 f    %随机生成多个初始值rx0_tmp,以避免局部最优
    - T7 Q7 v3 U5 l: N" g    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    $ @& C1 F8 O6 X9 L7 N    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    8 o6 k4 ~5 ^8 r7 C) v7 \1 @    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    ' u' ]5 V0 k: g0 M  G    rx0_tmp(4) = rand*1000;
    7 r% D6 q8 Q$ T8 _9 e    rx0_tmp(5) = rand*100;
    * c* V1 K8 w% W* ?& A! v    rx0_tmp(6) = rand*50;
    & n9 u) A( w  g2 b$ X    % 寻优& l4 l9 l+ @4 l
        [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...4 @. _* k6 C+ W9 j3 U
                fmincon('cpiao',rx0_tmp,A,b,...
    & P. W7 ^( v& y. j) O0 r+ c                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);  \2 b' C" ~$ E& \4 P9 y7 p
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
    4 {0 Q0 A  M1 ?; E  B        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。9 C5 }4 g/ X, Z, G: Q) }% w
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
    ! b/ }9 m7 y7 A" j        fval_meta_result = fval_tmp;6 ~) K, C9 ^& b. S: p! F
            rx_meta_result = rx_tmp;( R' d  j' @% s* K" Z
            flag_meta_result = 1;
    ) @* `1 f# ]2 T        if DEBUG == 1- j, S2 e. V8 @9 u9 c
                output_meta_result = output_tmp;' N) w9 R8 X0 _' ]$ Q  x
            end
    : }6 }( ]' N4 n9 }/ ^. a5 J3 Q    end& E" r: _4 g7 u- e. z( b' ^
    end
    1 o) u& a. u$ m% 把求得的最好结果保存下来" _* \" w* l+ @* D8 l
    if ~isnan(flag_meta_result)8 m" b7 _6 d% b) R! k7 j8 e
        rx_k1 = rx_meta_result;% ^: L7 M$ E- |  I; i( A! J- i" u* m
        fval_k1 = fval_meta_result;+ z  K0 {8 c; u: q' B) {
        flag_k1 = flag_meta_result;
    1 u! q/ F8 j" J    if DEBUG == 13 Z5 `5 U. u4 q# I( J* l; `2 Z
            output = output_meta_result;
    7 v8 X% ^% \1 o- T9 U! k    end
    4 K, o+ @# r* Gelse. _! x/ Q( o3 N1 M% \7 @
        if DEBUG == 19 n9 k; n' v$ f6 A
            rx_k1 = rx_tmp;. t, E# C5 r3 A7 C
            fval_k1 = fval_tmp;1 ?' Z, X0 O0 _; y8 K# E+ u
            flag_k1 = flag_tmp;
    - B8 R$ |  @4 C& c4 J        output = output_tmp;7 z4 E% s, Y; C- o5 X
        end8 ]; U, c+ C' W4 b3 `
    end
    % `! S# D- f$ h
    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-4-20 10:35 , Processed in 0.435767 second(s), 59 queries .

    回顶部