QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2967|回复: 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看懂是什么意思  求大神指点!!
    % n& G) A$ ]' Y6 V2 j8 I: o1 E: s% }
    8 E# l( a% k% m3 X2 Gglobal v4 r2 I" o% ^, U& p& ?- ~7 l
    v = 630589;                % 求解v为630589的收入水平情况
    " l& ~% _: ~8 b0 xDEBUG = 0;$ I1 Y% T( m: F2 @2 c- C& R. d+ `
    rand('state',sum(100*clock))        % 初始化随机数生成器
    7 }$ f3 a/ K# Fformat long g
    ) J$ K) x% k$ U: N
    3 j* ^/ z1 N/ Y6 h: S( T%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    5 ]2 w: S. U& b5 J% d3 v5 O% 求解开始
    / B2 v- W/ x0 J, q% 对于K1型
    ; k$ H1 ~2 ~4 p6 lp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
    - d# q! y6 l3 }9 o' R# C0 e% 6个奖项6个变量
    % K/ b+ S3 {7 V5 `9 _5 Z# wAeq=[1,1,1,0,0,0];7 p5 @, v: w( ]% ^  R) W5 L  i
    beq=1;# H& w( z& n9 L  O
    a_lb=[10,4,3,4,2];( c3 W( C# T1 q! m8 q
    b_ub=[233,54,17,20,10];
    0 }5 v6 z+ w& _* }/ [/ HA= [0,0,0,-1,a_lb(4),0;* U% R' J4 w% p" [
        0,0,0,1,-b_ub(4),0;
    5 h$ i4 v' t# g/ ~  V    0,0,0,0,-1,a_lb(5);  ]. \1 e& L& Z! I
        0,0,0,0,1,-b_ub(5)];
    . P& W7 Z5 J: s' N7 G; Hb= [0;0;0;0];# k6 e( U3 H8 k& Z) p1 O
    lb=[0.5;0;0;0;0;0];4 V7 [" }) V" P" U: ^! M: |7 v
    ub=[0.8;1;1;inf;inf;inf];; x) T% q; }. A0 b
    p_test = p_k1;1 L' b* _4 O. n: P
    rx0_tmp = zeros(6,1);( i' n6 Z$ S2 b
    rx_meta_result = zeros(6,1);
    " |5 w/ K- W0 [$ i. m3 M; J* Mfval_meta_result = inf;; s  Z; q5 R: t) Y; R( \
    flag_meta_result = nan; %用以判断有没有得到过可行解1 S3 ~' [% H' c% v9 F
    if DEBUG == 1
    : U% B# q) L: ~5 W0 f    output_meta_result = [];/ q3 `+ U& K: R" j8 q' B$ ?
    end. j! t& [% }! D" A0 g
    for j = 1:nums_test_of_initial_value+ u) ]! Y# Q  E) B0 b
        %随机生成多个初始值rx0_tmp,以避免局部最优
    * L5 \+ F5 q, f& M  ~2 Z$ d    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    3 z) M5 [8 ]  O% x' a! X    rx0_tmp(2) = rand*(1-rx0_tmp(1));4 c1 S0 ^- k+ k1 w: y1 B. W! I/ z
        rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);1 r6 x8 M5 `6 C* Z$ y% L1 f
        rx0_tmp(4) = rand*1000;
    : W, \" F: o7 {" r5 u    rx0_tmp(5) = rand*100;# o% d8 z$ F# ]# G6 e9 R
        rx0_tmp(6) = rand*50;3 R' h7 w' J  o- @1 ~! \
        % 寻优
    : n9 o' c; s; {+ ^8 e6 ?' l    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    # l) B3 b" a5 x/ K5 Q! G            fmincon('cpiao',rx0_tmp,A,b,...% B3 E- h3 w+ {$ w$ F
                        Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
    : h9 R, _2 @: Q3 E$ z        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)! w' C" `: w1 e$ R/ Q) m- R
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。( i$ o% U: Z% W1 j7 M- t- L& F
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
    + D0 h$ O9 U. |        fval_meta_result = fval_tmp;& @' `' L, u6 c
            rx_meta_result = rx_tmp;
    : b3 i- i! N% z7 J+ E$ E" P        flag_meta_result = 1;! B3 N# i/ i3 }5 b8 \' k9 H, C
            if DEBUG == 1; v% g8 H- j, |/ S9 M: b$ {( V
                output_meta_result = output_tmp;
    1 I5 |* X; ^0 M3 r        end
    4 l& c7 D4 t: v& z1 I! D, E0 s; v    end  w8 g" t! j8 X+ ~; w2 \1 G3 d
    end- l9 `$ _6 I3 g+ g* t% G
    % 把求得的最好结果保存下来; ~' p: Z$ y, F% \% ]
    if ~isnan(flag_meta_result). n" Y0 L7 v  T& h
        rx_k1 = rx_meta_result;2 B. k- b  h0 O* v, |
        fval_k1 = fval_meta_result;4 R6 f1 \) \# W& c% e. h
        flag_k1 = flag_meta_result;
    6 F- s1 [) c( b: E/ d    if DEBUG == 1
    % [  h3 n: y; E1 ~7 r  ]' \* m        output = output_meta_result;7 p/ m+ w* I. Z7 U& S
        end3 o* z, n' p# U- w
    else
    7 v4 M% v$ `' D+ C/ i; {- N9 y    if DEBUG == 1$ D! |$ q' r& {+ R9 U: F: u! v
            rx_k1 = rx_tmp;
    , q( t: X/ o3 z! W        fval_k1 = fval_tmp;
    ( }3 |% o- }( [5 \% G/ P        flag_k1 = flag_tmp;
      B" S, M1 Y  H5 g/ i        output = output_tmp;8 Y# r# Z  R& n* E7 P# ?
        end1 `6 z  y' }$ ^3 ~- I
    end+ j# y; w4 l5 a" I. U7 S
    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-14 20:17 , Processed in 0.435931 second(s), 59 queries .

    回顶部