QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3000|回复: 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看懂是什么意思  求大神指点!!
    5 Q1 S* H" \: r5 W4 Z( s
    " V% W5 F/ g# e1 Tglobal v( b2 c# `$ ?* D( d' j
    v = 630589;                % 求解v为630589的收入水平情况
    " h2 R1 p0 z5 EDEBUG = 0;- W/ \4 e0 C9 {2 G( |9 |+ ]. H
    rand('state',sum(100*clock))        % 初始化随机数生成器, p1 d: w4 Z4 a9 ^+ C
    format long g
    : w9 l  S: {5 T- |/ u" _& t/ h9 G+ u" v7 ?& ?' ^& `: ]
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1 n  z: h2 x* F& O& t3 S1 g" m: n
    % 求解开始
    5 u! A& `; z$ B% 对于K1型
    ! Q  z% E( V/ i# o6 qp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
    3 ~  U; o5 A+ l, w# _, b- G3 \9 @% 6个奖项6个变量
    % I  g; k  R9 ^3 c9 s5 g- F# SAeq=[1,1,1,0,0,0];
    ) `% Z# R! \, M, X9 z0 o. V+ @. tbeq=1;0 w4 P: d9 m/ t/ y/ V* L" p
    a_lb=[10,4,3,4,2];3 f) I! m( J8 E8 C- ]
    b_ub=[233,54,17,20,10];
    ; ~, ?, x0 y+ m$ cA= [0,0,0,-1,a_lb(4),0;
      ~; P- X$ Q  e3 l  H. p    0,0,0,1,-b_ub(4),0;
    2 p) G9 F! Y3 V    0,0,0,0,-1,a_lb(5);8 r# ]+ T) |3 q
        0,0,0,0,1,-b_ub(5)];
    7 p( |, Q( e: [( h/ xb= [0;0;0;0];6 G( y5 W7 J  w4 w3 ^; E; m$ C
    lb=[0.5;0;0;0;0;0];$ \: M+ _% R, |( Y+ E+ _
    ub=[0.8;1;1;inf;inf;inf];, s& J7 N% |6 @# T* J/ l) Q1 D
    p_test = p_k1;
    9 R# S, t& l) k  P1 V" w6 u, Y5 ~; Rrx0_tmp = zeros(6,1);
      b6 p. F# ^( ?" ~rx_meta_result = zeros(6,1);7 T/ T: }# I# w; T
    fval_meta_result = inf;  A- w1 e' w+ G
    flag_meta_result = nan; %用以判断有没有得到过可行解* Q9 R& q; C6 G
    if DEBUG == 1- ^1 ]6 \5 B% o. E% _8 s" m
        output_meta_result = [];9 }7 Z% Z- m! f; G, l$ L5 r
    end
    ( G( H( a6 |1 D4 l: rfor j = 1:nums_test_of_initial_value' i* t, w* d! V) v+ ?* a6 d
        %随机生成多个初始值rx0_tmp,以避免局部最优
    ) R; v% B2 @$ G, q( t+ n, v6 u7 e    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    / O0 B+ @9 N+ j' d& u  H! S! b    rx0_tmp(2) = rand*(1-rx0_tmp(1));4 S. `2 t) i6 M* ~
        rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    + Z4 ]0 m. C0 x$ a  m" F- X    rx0_tmp(4) = rand*1000;
    # x2 M1 L: N! _* m: ]    rx0_tmp(5) = rand*100;
    8 w# R, e. x8 A7 R( b$ O    rx0_tmp(6) = rand*50;
    ( v) \0 [5 V8 L' T    % 寻优
    % T+ X* t; R( H' O" S5 ~    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    ' ?: D7 p4 b3 z0 N* O* D7 w1 i            fmincon('cpiao',rx0_tmp,A,b,...
    ' B3 K) e, _! y# d! o' E                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);: o. W: I$ c4 W/ J; S# J
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)* F5 [+ E9 J* I* k% b
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    , c0 g+ U1 M5 A+ _* T    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)3 u; ~  n7 r2 r! B
            fval_meta_result = fval_tmp;
    : x) H* g. B8 O7 a        rx_meta_result = rx_tmp;
    3 D* S' i( r5 I+ V# g$ Q) L. D6 ^4 ?        flag_meta_result = 1;4 \& z2 V( q% L( P/ ^5 p
            if DEBUG == 1
    1 q. F& F4 k+ m; D: ~3 c$ l( |% ?/ B            output_meta_result = output_tmp;
    4 e, `. Y. L+ m! E        end/ ^8 L' T0 ]8 Z
        end
    - x! ]% [, v* s1 {end4 V( B4 z; U" ]; R3 R& T0 n5 I
    % 把求得的最好结果保存下来
    ! z% G- g, j1 E  p$ G: m7 i% Cif ~isnan(flag_meta_result)* R7 ^  U( A$ w5 p  S) l5 u# U
        rx_k1 = rx_meta_result;
    ) _) v, d; d: ?( }( X    fval_k1 = fval_meta_result;
    8 T+ O& S1 k) X) o4 Z( Q    flag_k1 = flag_meta_result;, C% K# x4 x4 S4 j7 x" A
        if DEBUG == 17 R7 j0 F, p! ^( N5 ], Q  \/ J
            output = output_meta_result;. Y& \- w( u# U1 e
        end
    ; C- Q$ W$ q/ S2 K% i. Selse" \5 @/ ?% }. `+ A
        if DEBUG == 1% ?% q3 d% b  t& s, {& Q/ w, ?" T
            rx_k1 = rx_tmp;
    ; B! u- h5 Z, @+ c, q2 o        fval_k1 = fval_tmp;5 E6 U/ g% @/ F, w- Z: M4 t+ t
            flag_k1 = flag_tmp;
    ! T  q' F9 W: p: C5 c1 V7 i        output = output_tmp;
    4 }0 T/ z- H8 Y% [% e7 N" v2 |  |    end
    0 P8 f2 M4 u, zend: a, j) w; z1 J/ A6 P. r
    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-15 13:47 , Processed in 0.425518 second(s), 58 queries .

    回顶部