关于debug的问题 。。。。

Bojay

这是02年CUMCU 彩票问题  debug看懂是什么意思  求大神指点！！

* l% G  y9 H* f6 j6 iglobal v
v = 630589;                % 求解v为630589的收入水平情况
2 ?4 n# x* ?9 i4 H) a, @6 @- HDEBUG = 0;
5 V. m7 p$ d# j3 Nrand('state',sum(100*clock))        % 初始化随机数生成器
format long g

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 P4 ?1 ]) y+ v) Y% 求解开始
% 对于K1型
7 A) o% j4 v. e6 ?8 w7 q; ^p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
* P: p: L1 H( I& _% 6个奖项6个变量
0 Q; p$ T& \) A4 T. YAeq=[1,1,1,0,0,0];
  E7 w- u! y) I0 j: l2 z+ tbeq=1;
a_lb=[10,4,3,4,2];
: p/ t# Y3 G  _+ O) Ub_ub=[233,54,17,20,10];
A= [0,0,0,-1,a_lb(4),0;
    0,0,0,1,-b_ub(4),0;
    0,0,0,0,-1,a_lb(5);
6 Y2 D- s  p/ B/ P, y    0,0,0,0,1,-b_ub(5)];
b= [0;0;0;0];
1 W& A' P1 M: w; ~" Plb=[0.5;0;0;0;0;0];
ub=[0.8;1;1;inf;inf;inf];
8 Q2 y' }* O0 H0 t2 ~: ]p_test = p_k1;
rx0_tmp = zeros(6,1);
) I2 s5 }  O& k+ P4 E  drx_meta_result = zeros(6,1);
' U9 C; S! {- Z  qfval_meta_result = inf;
flag_meta_result = nan; %用以判断有没有得到过可行解
if DEBUG == 1
  V' F" k# m( B+ C/ w$ ~    output_meta_result = [];
end
1 |0 h! t7 u2 `- T) a4 y, Pfor j = 1:nums_test_of_initial_value
$ ]8 V, l1 K. _# c: c    %随机生成多个初始值rx0_tmp，以避免局部最优
6 K1 F7 X, T0 {, z* ]# p    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
+ a* N4 i! Z  y- `; H    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
& ^6 i; f3 ~5 Q$ \% t; H. F    rx0_tmp(4) = rand*1000;
* ^% a) P4 ^, f) P* o3 g' \' K    rx0_tmp(5) = rand*100;
4 Y7 [3 s: n4 v" c' a  e1 c    rx0_tmp(6) = rand*50;
    % 寻优
! y7 j7 q+ P/ Y5 X    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
" y4 I! ?5 @3 S; N' m/ i; c" A            fmincon('cpiao',rx0_tmp,A,b,...
                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
! w% ?  B' o# t5 r1 ~    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
  E; g2 @! P/ R* |2 S        fval_meta_result = fval_tmp;
        rx_meta_result = rx_tmp;
        flag_meta_result = 1;
        if DEBUG == 1
            output_meta_result = output_tmp;
        end
    end
/ G* ~# k/ S- Rend
5 Y0 A6 Z+ d7 ~4 }% 把求得的最好结果保存下来
if ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
4 J" e' l$ z2 ?/ Y7 [1 q    fval_k1 = fval_meta_result;
: _/ g  O1 e6 J    flag_k1 = flag_meta_result;
2 C+ N; ?; ^, l, r" m1 q    if DEBUG == 1
        output = output_meta_result;
1 O& }% `: ~# w- z( k' k4 I    end
% o, p& W" N: P) y4 `0 d1 r  eelse
    if DEBUG == 1
. |* }0 I6 I) t8 r/ r        rx_k1 = rx_tmp;
+ g" ^# E" c8 R: u/ ?) n! x        fval_k1 = fval_tmp;
- [/ j& q4 K1 d9 C' K' J# g- [        flag_k1 = flag_tmp;
        output = output_tmp;
    end
end

Bojay

没人回答啊