QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2999|回复: 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看懂是什么意思  求大神指点!!* T  ~- L5 }- d; V2 w4 o

    4 N( h9 k) ]1 H3 s, e& l! E1 l; |global v
    5 I, P) @0 K- u; w0 uv = 630589;                % 求解v为630589的收入水平情况/ A* J- w* n! ^! [; b3 K
    DEBUG = 0;
    ) U$ H% V4 t7 ^7 K1 \7 F% Trand('state',sum(100*clock))        % 初始化随机数生成器; w9 O' p2 }; x* \& h. n0 o' W
    format long g& b0 }# t. Q8 k: J
    8 a! p$ [# r( g; [. N
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    9 @' k0 y- x" {8 ~% 求解开始
    ) l3 }9 s  J4 z- @% 对于K1型
    / A4 U( X( C9 T6 x% tp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
    / y+ P6 K$ T6 z0 g) D% 6个奖项6个变量$ P, @" z% I. H% t" V8 N
    Aeq=[1,1,1,0,0,0];3 g: ]' ?& G  ~) ?0 a5 q6 ^
    beq=1;/ g* I; r$ S0 n, S
    a_lb=[10,4,3,4,2];
    ; E9 i# G& `" f3 Q* x  V; mb_ub=[233,54,17,20,10];
    9 o0 U& F1 E1 n# y# HA= [0,0,0,-1,a_lb(4),0;- `( F  v* D9 C2 q# N! d
        0,0,0,1,-b_ub(4),0;, ^4 J1 v( a( N5 A. M
        0,0,0,0,-1,a_lb(5);+ H+ g( _$ ?1 Z; H, i
        0,0,0,0,1,-b_ub(5)];
    6 O  V0 z+ X9 F9 wb= [0;0;0;0];, w, X0 Q( P9 X
    lb=[0.5;0;0;0;0;0];
    ; T$ _" t  t  G% Z" {7 ]# bub=[0.8;1;1;inf;inf;inf];
    , V+ W: I% @4 g! ^0 bp_test = p_k1;" N) ^6 W& P( B, @$ P; J- [! u* P
    rx0_tmp = zeros(6,1);
    , C6 e4 l5 W/ P. W6 xrx_meta_result = zeros(6,1);
    2 A% v8 w# ^7 e6 C2 z) yfval_meta_result = inf;
    - }7 x2 p' C) X, {. nflag_meta_result = nan; %用以判断有没有得到过可行解: {  ^9 f2 Z8 \5 z2 @2 ^% Y
    if DEBUG == 1- F4 |" H% y# y7 z3 l& _
        output_meta_result = [];
    2 [3 _/ u5 f1 v2 \0 S) d9 y# Oend  s; L3 ?9 ?. g3 W
    for j = 1:nums_test_of_initial_value
    + Y0 y( Q) o, k    %随机生成多个初始值rx0_tmp,以避免局部最优
    7 x8 |7 R. ~8 Q9 ?0 Y* ^& q/ t    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;# p) x) I- l7 [& x% Q
        rx0_tmp(2) = rand*(1-rx0_tmp(1));3 j. M  ~1 w0 |( {5 j: ?$ ]4 J
        rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);. ~% }5 C8 O' R( ^0 m9 ]
        rx0_tmp(4) = rand*1000;+ U/ u* `6 X+ G' {7 j
        rx0_tmp(5) = rand*100;& C( a1 S4 L, _  v0 \5 ^
        rx0_tmp(6) = rand*50;
    . M/ H% ?) ^4 m. @1 u    % 寻优( r6 q4 o% u& l/ k9 X3 L
        [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
    - g  O, w6 w/ Y            fmincon('cpiao',rx0_tmp,A,b,...; o$ O/ ?5 A; H+ W9 c; d: [9 H) C/ B
                        Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
    ( H' U  }6 }1 D' A) h( b. V1 N        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)4 m; R% J( u$ o1 ?& H+ s) [
            % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。& V# g" m; X8 x. ~( C- S8 f
        if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
    ) Z' T2 d. Y; @: A  w: w        fval_meta_result = fval_tmp;+ v1 [6 F; m* d( q
            rx_meta_result = rx_tmp;
    , `5 e6 v: s! l" {/ p7 W5 P        flag_meta_result = 1;
    " g5 i+ P8 }2 L        if DEBUG == 1& B8 |5 h% I1 @* P4 F. l5 F
                output_meta_result = output_tmp;* V4 E" G) c+ G2 C
            end
    3 g/ A; `8 u; Z- A7 m3 Q    end7 ?% v, A2 O- U( i5 l7 E6 c: h9 c0 ]
    end
    3 b$ c7 ~% D% e% 把求得的最好结果保存下来9 v* K9 D! q  w( D/ r
    if ~isnan(flag_meta_result)
    * W4 {* U) |' s% D' [% K& O    rx_k1 = rx_meta_result;
    6 w# K8 m: i4 U$ N  K3 F; E3 r    fval_k1 = fval_meta_result;
    : Z5 e) P) V6 h! ?; V5 h1 `    flag_k1 = flag_meta_result;
    ) C) K- A+ c" ~; C( g    if DEBUG == 11 v9 {# h2 f; d% X9 `
            output = output_meta_result;3 Z# K  r2 i- n/ P2 Y
        end  d/ `& j! j5 y
    else  N; i, q  v- R. H( y( \/ c/ k
        if DEBUG == 1: L' Q% N* b6 U( M7 ~7 [; I7 O
            rx_k1 = rx_tmp;
    " b: g- ^2 U' V        fval_k1 = fval_tmp;
    / [/ h7 U9 w; y7 g; d0 B        flag_k1 = flag_tmp;
    % \2 ]9 S8 m$ i' V+ S% _# N$ ?        output = output_tmp;& O' a/ z9 v7 @% s
        end: N' l* |0 A  R
    end% ~: q% q. Y, Z5 s
    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 17:34 , Processed in 0.417772 second(s), 59 queries .

    回顶部