QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2864|回复: 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看懂是什么意思  求大神指点!!/ [' I$ ~. `  _

    % f+ Y1 _* l; |, Nglobal v: I" |6 W! N& u7 t
    v = 630589;                % 求解v为630589的收入水平情况
    9 Z4 F$ J0 l* @8 ]7 ^DEBUG = 0;4 X2 i& k& A+ [' @* p( n. X
    rand('state',sum(100*clock))        % 初始化随机数生成器0 E2 Z& P! o/ |, w
    format long g3 n& `$ {  R* N
    5 R6 m5 m# q! n# r. ^
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%& Y  V$ r5 f: w* E  T+ Y/ n6 j
    % 求解开始1 y" x( `; P2 H3 N8 j# }- {$ t2 }
    % 对于K1型
    9 X% Y' H$ m6 e- g: U, z- pp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];0 S& L6 D* `* ~" B# Z; ?
    % 6个奖项6个变量* g$ y- Q. X6 @; _; W3 s5 I) \; b
    Aeq=[1,1,1,0,0,0];
    . ~* |; H9 L) X; ~% f8 sbeq=1;
    ( W2 L: ?; m2 j! g& aa_lb=[10,4,3,4,2];
    ( L8 G' n; ^4 r) H" ]' u  ub_ub=[233,54,17,20,10];: h, a- s; |3 O7 r. x  k
    A= [0,0,0,-1,a_lb(4),0;
    " {& ?1 F$ Z) Y' d6 H: {) l$ l( k    0,0,0,1,-b_ub(4),0;
    1 m& L6 }+ D4 X- L4 v2 Y    0,0,0,0,-1,a_lb(5);
    0 N: G, J" Y: P    0,0,0,0,1,-b_ub(5)];
    * F0 Y2 N% C- }% f2 Tb= [0;0;0;0];
    . C. e* H% I/ E8 {& C$ h6 \9 a: Rlb=[0.5;0;0;0;0;0];
    ' Z! `+ }3 N  F* d7 D! Q$ X3 K% g- Kub=[0.8;1;1;inf;inf;inf];
    ; L( d5 K' d9 p( m+ j8 yp_test = p_k1;
    % g0 o# x2 n# D$ g5 W+ i5 grx0_tmp = zeros(6,1);
    . N0 J1 c) i6 @* ?8 S" vrx_meta_result = zeros(6,1);3 @! P9 t2 A( e" ?5 e1 y7 ]
    fval_meta_result = inf;
    4 S& T: M1 @5 f& pflag_meta_result = nan; %用以判断有没有得到过可行解4 s9 Q8 S9 ]+ {# `  Q0 v- Z9 c& b
    if DEBUG == 1, ?- J0 t+ k- O% l3 D# P
        output_meta_result = [];
    . A5 x% o3 e, |2 L3 P% i: Zend
    $ u" {8 C- j  b: {6 ~3 P" f" J6 Yfor j = 1:nums_test_of_initial_value  S3 T( q- B8 E' q
        %随机生成多个初始值rx0_tmp,以避免局部最优
    / Y, [3 n+ E6 }9 s    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;8 F4 q' G0 \, O2 ?: {" Q
        rx0_tmp(2) = rand*(1-rx0_tmp(1));& J' L9 V. q3 Q0 B. {1 [% |
        rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);3 G2 X  @* D6 t
        rx0_tmp(4) = rand*1000;! }# h0 v6 j7 W. y
        rx0_tmp(5) = rand*100;
      Z+ A& T- H, ^! n: o, ^    rx0_tmp(6) = rand*50;
    * V) z- _) v% x0 `4 p) o: k' C    % 寻优
    1 f! X  W  M% _* e    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    , |$ l5 R0 G' O6 v1 }9 c& @            fmincon('cpiao',rx0_tmp,A,b,...
    4 Z: O9 E1 J+ D/ O* G2 s2 M: o                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);) b1 j7 p: ^2 R$ m
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
      P- f  x, J* g3 Z# r: q        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。) k- R  a3 P  e$ |2 }  d+ K% ^
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)( g# o) `* p# B9 N! _
            fval_meta_result = fval_tmp;, l4 j9 v) I# ^; e% e. K" x$ \
            rx_meta_result = rx_tmp;% B$ U$ p4 s+ k/ h0 `
            flag_meta_result = 1;
    * M4 d( S/ L. |2 q# A* A        if DEBUG == 1# ]7 Z! g( K4 r) u
                output_meta_result = output_tmp;
    . N$ c  \1 m* z9 G. I. _1 C        end
    # t( {% I4 C8 z8 i! C( c; z0 R% G& F    end; }0 h9 f/ M0 Z% }, ?7 Z
    end, T. p# p- n+ F  m; {$ @
    % 把求得的最好结果保存下来: I6 z6 E* x9 a! D9 Y0 J0 e. ~
    if ~isnan(flag_meta_result)
    ; b0 _+ g- z8 e& x" j/ ], `( |    rx_k1 = rx_meta_result;; m) z* b% b8 b; Z/ c( {
        fval_k1 = fval_meta_result;
    ! u- y2 d! n  F* m! {; n; Q7 @# t    flag_k1 = flag_meta_result;9 p% f! O8 g% D
        if DEBUG == 1
    , O6 ^, b2 k+ o4 m" K        output = output_meta_result;0 Y- L: F7 k' m/ v/ Y6 t
        end
    ( q& m9 K2 }  p7 o- Melse
    ' I  A9 Z8 \+ b, D! q& C    if DEBUG == 1$ T, M, K! ^9 B& C( ?- @% n( J' v3 {
            rx_k1 = rx_tmp;% A! p2 N3 T7 n4 s$ `, @1 J" W, W/ c
            fval_k1 = fval_tmp;
    5 z* {0 y% B$ P: a0 t        flag_k1 = flag_tmp;
    7 |3 o& `& {# f+ t$ E/ a1 [& |        output = output_tmp;
    % S% H! a/ v  P5 F3 T& h    end
    0 p0 ]6 e& E) }# w1 O0 K1 Fend; m9 n7 u  `. W7 I$ s. }) 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, 2025-12-4 14:49 , Processed in 1.764878 second(s), 58 queries .

    回顶部