QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3029|回复: 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看懂是什么意思  求大神指点!!
    ; @: r+ C$ \6 j; r2 s$ e0 w2 T) r+ Q8 p6 q' \
    global v
    ' i/ g  d% k6 U% H. T8 Z/ d/ I5 @v = 630589;                % 求解v为630589的收入水平情况% Y7 `& E- j# B
    DEBUG = 0;- |7 K) j( ~0 [1 y& f- M8 C
    rand('state',sum(100*clock))        % 初始化随机数生成器
    5 v6 K$ ^7 Q7 k/ X1 P) D1 y& x" B! T1 D$ @format long g
    2 T4 ]* p/ c" n* [* A% p1 D
    " Z( h) l. F" T# H/ F  K%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    - [; r" K) ^$ G. j8 J* ?# q% 求解开始# t* l9 L  |- E6 X
    % 对于K1型
    : _4 \0 P" o- y6 zp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];) }1 }; n. C0 c3 a, `* I
    % 6个奖项6个变量/ }7 n3 n$ v  q3 B1 P' v' S, q
    Aeq=[1,1,1,0,0,0];
    / k# S) A  G: N0 zbeq=1;
    4 J, ]$ O2 w' h# r. sa_lb=[10,4,3,4,2];; ^2 J* ^& k% L) X' p  s1 L1 [$ W
    b_ub=[233,54,17,20,10];5 r- O% i$ k7 m/ O
    A= [0,0,0,-1,a_lb(4),0;% @/ S+ ]9 b% l6 n3 c" L/ y
        0,0,0,1,-b_ub(4),0;" f4 V* t' ~+ j) @& i9 J. d% n
        0,0,0,0,-1,a_lb(5);
    , J. a7 \0 e4 Z! n    0,0,0,0,1,-b_ub(5)];# B2 _% K7 C, `0 o
    b= [0;0;0;0];
    9 Z! s' }+ s' v+ v6 M5 plb=[0.5;0;0;0;0;0];
    6 {. b- K3 a% d- W' Pub=[0.8;1;1;inf;inf;inf];
      g, m+ Q( X, w, u$ i, e& dp_test = p_k1;+ Q1 c2 N' x7 R
    rx0_tmp = zeros(6,1);
    & G3 j* Q( I, b4 h  p6 Erx_meta_result = zeros(6,1);
    8 w! {3 N8 I8 t$ Qfval_meta_result = inf;' v& i+ Z' t5 i8 D
    flag_meta_result = nan; %用以判断有没有得到过可行解3 v; n- N/ y, _& b+ h( A  F
    if DEBUG == 1
    ' b$ Q, A5 ?1 x+ w  N" u    output_meta_result = [];
    # X! i+ d5 [% G3 V; lend
    & |) Y. p5 d; s* L) p* g8 a/ Dfor j = 1:nums_test_of_initial_value' |9 _" _9 h5 \
        %随机生成多个初始值rx0_tmp,以避免局部最优
    ( ]* A4 s2 Z5 M( G+ r    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    . Q1 @8 I* c+ k: G& Z  g    rx0_tmp(2) = rand*(1-rx0_tmp(1));/ X% c8 A% Q6 @9 ^3 v4 ]
        rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    * r+ ]& W! D2 E6 a    rx0_tmp(4) = rand*1000;- u6 R& H* \9 _& W
        rx0_tmp(5) = rand*100;9 ^" [$ `& V$ y4 {" P0 N
        rx0_tmp(6) = rand*50;
    2 _* A! t/ E1 y" y    % 寻优
    ' B8 @+ s/ K7 j) ?. P6 M    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    9 _( T) ?" U7 m+ J+ x            fmincon('cpiao',rx0_tmp,A,b,...1 V& c6 j* H7 w6 Z+ [' Q) v6 t3 n
                        Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
    1 F; J  J) M( j" j3 Q5 U- |        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
    # S/ g! o6 P9 Q. G" ?        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    . D1 j, m6 c2 u2 k1 H+ h8 j  `    if (flag_tmp == 1) && (fval_meta_result > fval_tmp). a* _( w* n4 l3 l7 G" V$ S
            fval_meta_result = fval_tmp;
    8 \( T. z7 X. j        rx_meta_result = rx_tmp;3 }4 P) u' T8 e" Y$ v
            flag_meta_result = 1;
    ; j! i# Z. o2 J0 G' {: `        if DEBUG == 1
    ( `) ]8 g0 B2 Q; ^- f5 {' f            output_meta_result = output_tmp;
    / x/ n; u# v" c0 Z. m# ^7 v1 h- w        end
    : Y5 X% f- J4 t    end
    7 n) }! B1 _- \; r/ Jend
    8 I: }  z9 J) R. U; o& w+ A' k. t% 把求得的最好结果保存下来
    * b) R- Z. d4 {if ~isnan(flag_meta_result)& o1 z+ L3 s) x' B2 `/ H- l! @
        rx_k1 = rx_meta_result;# ?- b/ N! E5 {: Q8 S# G$ k' S0 D
        fval_k1 = fval_meta_result;
    9 P7 b( F9 n8 U+ I) N, p& B    flag_k1 = flag_meta_result;
    ; Z( X* ]& [9 x2 P" o  q  d    if DEBUG == 1
    - Y3 A, Q- R; z# i        output = output_meta_result;
    5 h* w- |. M) L    end
    2 h& i) L# m1 V2 s9 C3 a2 J2 Telse
    ( ]/ E( c& g- v8 E( s2 f6 O    if DEBUG == 1
    7 m2 N( ?4 W0 P' ~8 L; {        rx_k1 = rx_tmp;
    * w1 r9 r1 T7 A        fval_k1 = fval_tmp;' E% e7 g/ V9 u( g- K% S
            flag_k1 = flag_tmp;
      v, Y8 `: s& I: S/ Y) y! R        output = output_tmp;# S) P4 Z9 v0 d: {+ _& R
        end" ?3 S& }. Q; o1 |: I
    end
    % `. z, n. h0 C$ @
    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 17:03 , Processed in 0.549009 second(s), 59 queries .

    回顶部