QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2995|回复: 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看懂是什么意思  求大神指点!!6 q' A; G9 e9 W& H) U+ ?

    9 m% q! s0 k( G" O- nglobal v" F0 A- L/ @8 n/ l( q8 R3 ~1 ^0 `) h
    v = 630589;                % 求解v为630589的收入水平情况
    & |. y- ]' K1 V+ s5 l6 B9 z/ i; iDEBUG = 0;4 p) k; v9 B! e5 G. x
    rand('state',sum(100*clock))        % 初始化随机数生成器& ^9 }) D$ B7 E: {
    format long g4 P$ P8 U) x; j6 ?

    4 V1 P: D9 X$ k- r! K7 M5 e6 J( k%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%; O$ _5 m7 V8 X; ]/ y- {
    % 求解开始/ u, {4 i3 s' m
    % 对于K1型
    , D# Y( u0 I2 W5 q0 cp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];9 |/ m: V! q5 o0 w
    % 6个奖项6个变量' S8 D* T: y  g; t4 G
    Aeq=[1,1,1,0,0,0];' R$ `& A" G9 j7 h% c
    beq=1;
    * W* f' b8 m5 {+ B: aa_lb=[10,4,3,4,2];* _. q  s/ S. K5 B" a- m, e+ q
    b_ub=[233,54,17,20,10];
    5 g; ?8 G( \* G; D0 R, a4 i4 ^A= [0,0,0,-1,a_lb(4),0;
    ) K- j0 ^& e' v6 i  Z7 j    0,0,0,1,-b_ub(4),0;( |# H' \0 B. c! `- Z
        0,0,0,0,-1,a_lb(5);5 U; E$ R9 P4 d8 Q4 a
        0,0,0,0,1,-b_ub(5)];
    1 y' v% _6 j5 M% Z' B. kb= [0;0;0;0];' D. @% {0 c# Q; ]5 h! P
    lb=[0.5;0;0;0;0;0];# n) [' G( ]/ m
    ub=[0.8;1;1;inf;inf;inf];$ d$ y! l6 W; [- Y
    p_test = p_k1;
    . u  `  h( ]9 G1 q* o& Y3 zrx0_tmp = zeros(6,1);
    ) V& w  |  m. q& srx_meta_result = zeros(6,1);
    9 d* ?' p/ _. H5 x7 d" k# _. wfval_meta_result = inf;
    8 T- m. _8 n. Fflag_meta_result = nan; %用以判断有没有得到过可行解
    6 }2 g- a) X6 g$ W9 eif DEBUG == 1
    $ w1 M" \2 x4 {1 d    output_meta_result = [];4 E7 _# s) t& s' R0 d
    end
    ; S  {5 y% C7 ?  Gfor j = 1:nums_test_of_initial_value* }6 V- u# X3 v
        %随机生成多个初始值rx0_tmp,以避免局部最优0 h- I% J3 p5 ^1 T3 G& a, m# b
        rx0_tmp(1) = rand*(0.8-0.5) + 0.5;- |* m* ?6 i% `: ]5 p  X" H9 K' y
        rx0_tmp(2) = rand*(1-rx0_tmp(1));
    / m; s6 }1 q! F9 V$ L8 e    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);7 x' t$ b7 B) _
        rx0_tmp(4) = rand*1000;5 ?6 Z/ n) o" @& t( R% d0 S- d
        rx0_tmp(5) = rand*100;) _2 I) M$ k# N# [/ c( T
        rx0_tmp(6) = rand*50;" G' u: \0 E& Y! C8 A1 @
        % 寻优
    ; s  P2 g3 _7 C( R" b    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...3 `: e7 j7 B4 {
                fmincon('cpiao',rx0_tmp,A,b,...
    + j2 ]9 ]+ b5 N, p& L* c                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
    * P5 m5 F7 W# z8 X        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
    : U- g- G/ ~- |; }3 U        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    . e$ ~7 }4 z2 Q' o+ H' [" w- g    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)" N) m2 m2 X3 v( A* f3 Z
            fval_meta_result = fval_tmp;
    + B0 ~7 f$ _% N5 L! o; v        rx_meta_result = rx_tmp;0 c1 f$ H! o0 B' z) b
            flag_meta_result = 1;
    1 s" X' [% U1 M5 M. x* ?9 P        if DEBUG == 1
    & K: `) [1 J* V$ s2 T8 {' ?- m            output_meta_result = output_tmp;
      V) N8 F* g- y7 ]" ]$ A3 w        end
    ) ]2 @8 |; x5 `! O    end9 b& c# O$ c# ^% Z8 m5 S+ P
    end* ~- ?3 ^) b% w
    % 把求得的最好结果保存下来
    ; C- E) ?$ G$ Y4 u0 [if ~isnan(flag_meta_result)
    % H  [* q& ^4 o/ N    rx_k1 = rx_meta_result;
    4 l% C/ B! B" i& \- D" U    fval_k1 = fval_meta_result;
    5 ~. D( Q. [8 _. I; z7 H7 d    flag_k1 = flag_meta_result;
    , j! m( C. u: D! D' h4 ^: E! }    if DEBUG == 1
    ' f! p$ w) B* W        output = output_meta_result;- T1 g! p  l. G9 w2 N5 A
        end6 ?. B. P: t3 T& m! y
    else+ `- K( y2 }2 g1 x! Y5 O; k7 T
        if DEBUG == 11 j* B. O; a" a$ K5 M+ \' t
            rx_k1 = rx_tmp;
    6 o$ o' a  h0 u. E        fval_k1 = fval_tmp;4 [) [7 w' d# b( G, g0 \' S
            flag_k1 = flag_tmp;
    * }2 \2 ?  |" _8 Q* W6 S        output = output_tmp;
    ' \* ]0 e! f" y; M, B* S1 J3 D% |+ V    end
    ; _8 @3 C, T. U, M; z3 s+ |end4 |" H0 r0 V2 {; e- q
    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-6-14 03:18 , Processed in 0.609460 second(s), 59 queries .

    回顶部