QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2969|回复: 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看懂是什么意思  求大神指点!!9 y" ~6 ^% _. C6 g4 P  |: j

    . y4 o& ~: T6 h# R0 R/ tglobal v
    # `/ i' G8 [( H$ m$ Q" [5 Jv = 630589;                % 求解v为630589的收入水平情况
    0 s! A2 p5 `& V! G' g# b) DDEBUG = 0;* _& F( F& \7 X% h0 f
    rand('state',sum(100*clock))        % 初始化随机数生成器
    , B2 S9 y# h& s+ R6 N% Z5 u' lformat long g
    2 `$ d9 T, J( }6 u1 F) M
    # o) c: \0 a% Z% m. a%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    + x, ~/ |# m" q8 _0 I1 u+ p# U% 求解开始- h( t3 E- M2 Z) X0 i+ O
    % 对于K1型
    - T9 ^9 [8 ?6 J0 b1 \$ Wp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
    ) o9 {. F& w7 _( n  a- M% 6个奖项6个变量9 y9 R& k; i8 D# Q: R3 K: k! e
    Aeq=[1,1,1,0,0,0];
    ( o2 R9 s* J% H# S4 ?beq=1;
    # G( w# S3 t' s& ca_lb=[10,4,3,4,2];* i! C9 ^: `  S: U( K
    b_ub=[233,54,17,20,10];
    3 a2 Y# {, r0 T* U$ |# U0 H2 zA= [0,0,0,-1,a_lb(4),0;$ m: e! h, [* L  h: V. a+ [! {
        0,0,0,1,-b_ub(4),0;6 c( f- H' u: d3 T& v) c
        0,0,0,0,-1,a_lb(5);) [) C: C; c$ d
        0,0,0,0,1,-b_ub(5)];) [5 R+ g& B- G" [5 x
    b= [0;0;0;0];8 Z6 l1 ?! K7 d% x/ h# Z
    lb=[0.5;0;0;0;0;0];
    * U/ @" w/ t/ L6 W. y2 S7 bub=[0.8;1;1;inf;inf;inf];8 F6 e$ e! }) V
    p_test = p_k1;
    4 a+ D7 S, f" L6 @; ]) xrx0_tmp = zeros(6,1);$ ~  h4 X. ?4 K# D9 |3 `
    rx_meta_result = zeros(6,1);( y1 I' u! i: V# o8 V
    fval_meta_result = inf;
    ; f& z- f( w0 K( s3 u+ f0 }8 wflag_meta_result = nan; %用以判断有没有得到过可行解) W* X5 ?0 O0 D; I  v0 q
    if DEBUG == 1% \& H: t2 C! ~- }2 F: u& B
        output_meta_result = [];
    & N7 d- @+ S' r- @0 Uend  a! q' v. m8 z6 ^* }: T( T% r
    for j = 1:nums_test_of_initial_value
    , g% @2 I: x. l+ g6 ~! n    %随机生成多个初始值rx0_tmp,以避免局部最优
    " |& D' d+ H: H9 I) ~4 R' m    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;' C/ {8 w: @. B
        rx0_tmp(2) = rand*(1-rx0_tmp(1));
    ; M; d  }: {. O; u    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);. a6 I. v; {3 B  a
        rx0_tmp(4) = rand*1000;7 }& l" T. x2 I# Y% A" v
        rx0_tmp(5) = rand*100;4 Z7 w0 k6 u- |. U5 |6 Y9 C' f
        rx0_tmp(6) = rand*50;
    , b& c" p; B# [* Z0 P0 ~; @    % 寻优0 W$ a2 y+ k1 W- M! s
        [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    % H, F3 }- Z2 o6 w7 \            fmincon('cpiao',rx0_tmp,A,b,...
    $ s% Q6 H; B4 C) q# l( D                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);$ d6 o( s$ t5 `
            % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
    & \) \/ ~' D: w- G        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。- O" q4 U4 y7 j6 v
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)5 W4 D5 ~, M. p/ G
            fval_meta_result = fval_tmp;
    * x8 c" Y7 Y2 U3 c        rx_meta_result = rx_tmp;
    & K" I8 l  k$ c7 G7 V, y- w! e        flag_meta_result = 1;
    5 U, X8 v6 v6 L+ [& g" r        if DEBUG == 1
    9 e; V; G9 g" L1 P# K5 W9 \* P  a            output_meta_result = output_tmp;
    - \6 D9 z7 i# y3 f8 L8 M# R9 S        end
    6 i& Y! s, d" Z9 [" m! ^    end
    3 D( u2 |% f: i5 Kend4 M) i. ?3 ], i
    % 把求得的最好结果保存下来
    4 \( d2 b9 g. H5 b" Yif ~isnan(flag_meta_result)
    8 g' C$ x9 z' T+ s6 o; F    rx_k1 = rx_meta_result;5 {6 O/ Q* Z- R( |1 K. P7 P
        fval_k1 = fval_meta_result;
    ; o' {: e/ h8 s9 [( `3 a    flag_k1 = flag_meta_result;
    # i1 A) \& }4 T- \    if DEBUG == 1
    - S# g- M2 J- Z  s  D) s        output = output_meta_result;
    , O' l. o( ?0 ^% H, ?    end
    / c- L1 d2 w7 W* P3 Belse
    1 y1 s( d% |' l3 o; O8 {    if DEBUG == 1
    ( X$ a, V' U! d& ^0 I        rx_k1 = rx_tmp;
    7 |8 P0 ]% T: j        fval_k1 = fval_tmp;
    4 f1 P- @) t: ?8 `, T  d        flag_k1 = flag_tmp;6 A3 g- d- F; ]$ J) I& U7 _  Y! s
            output = output_tmp;& V0 m  l  u2 M# U. K& j
        end$ V; @: B8 T3 q! F; E+ ^9 J3 u
    end
    4 \# I3 q6 T4 c' m+ \# A9 ]# E4 a
    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-15 14:20 , Processed in 0.437872 second(s), 62 queries .

    回顶部