关于debug的问题 。。。。

Bojay

这是02年CUMCU 彩票问题  debug看懂是什么意思  求大神指点！！
2 b8 T  A/ j7 N' m4 y
& o6 c4 e7 P1 \" Y4 A' z% Fglobal v
( q# i7 o4 n* \4 X% c8 Wv = 630589;                % 求解v为630589的收入水平情况
* q7 o$ G! a' K. I1 Z8 IDEBUG = 0;
0 ]" V* Z  L! v$ Nrand('state',sum(100*clock))        % 初始化随机数生成器
$ B0 }: e& }# T  e) I1 l, t, Pformat long g
. u$ c$ U( K" c
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
' {4 d6 I4 V# M' V% 求解开始
% 对于K1型
, t4 A' Z& t  fp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
* o. |- f" p! m% 6个奖项6个变量
2 Q3 v3 @# P- X( o, wAeq=[1,1,1,0,0,0];
. h( \0 d+ l- K/ y$ Rbeq=1;
, B# F$ c3 d* C  C9 J) Y* ja_lb=[10,4,3,4,2];
b_ub=[233,54,17,20,10];
$ c# y. N- R6 x- d' S& |$ QA= [0,0,0,-1,a_lb(4),0;
    0,0,0,1,-b_ub(4),0;
  Z: {5 |; B9 ~6 E    0,0,0,0,-1,a_lb(5);
    0,0,0,0,1,-b_ub(5)];
4 Z) v" H8 \- [& Q3 `% }b= [0;0;0;0];
  v& C# J, C+ v" o$ @7 G0 T/ W. Xlb=[0.5;0;0;0;0;0];
ub=[0.8;1;1;inf;inf;inf];
p_test = p_k1;
* E" |- K4 l5 U, b( J+ l3 Erx0_tmp = zeros(6,1);
- X# U3 U1 @8 rrx_meta_result = zeros(6,1);
( n, v/ _* g4 s' i8 n6 s! x0 |fval_meta_result = inf;
flag_meta_result = nan; %用以判断有没有得到过可行解
6 I; C  s' s2 i* S' h5 tif DEBUG == 1
, g5 Q8 p% Y* W' P% S7 p) g. |  J* @1 k    output_meta_result = [];
9 w3 |3 O# {' l) L4 h& uend
3 ^/ q; M7 O' {+ S' Qfor j = 1:nums_test_of_initial_value
    %随机生成多个初始值rx0_tmp，以避免局部最优
    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
* Y* \; d" Q# r) X  r    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
/ l6 }& @0 Y9 C4 W: {' K. P5 P    rx0_tmp(4) = rand*1000;
    rx0_tmp(5) = rand*100;
5 y1 l) v5 T5 Z4 X+ w' X    rx0_tmp(6) = rand*50;
/ n9 w7 \# s$ z5 P% b2 D) C# w    % 寻优
8 O; |( l( ~* m1 H4 C    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
            fmincon('cpiao',rx0_tmp,A,b,...
                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
+ x# b1 G! e& i, k8 R: C! @. F    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
  e5 }! p$ {! @0 {4 i        fval_meta_result = fval_tmp;
. p& r- ~4 C/ f) x+ w7 E* W        rx_meta_result = rx_tmp;
        flag_meta_result = 1;
. [0 n/ t4 ~* X& P        if DEBUG == 1
, E1 `1 S4 `  S: c6 y7 Q2 i            output_meta_result = output_tmp;
        end
" _( ~( [# y+ `1 [    end
end
3 ~' I( R4 |8 Z, x- T% 把求得的最好结果保存下来
* Q1 t- `: S) z" zif ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
  V' d7 O3 K0 t    fval_k1 = fval_meta_result;
    flag_k1 = flag_meta_result;
    if DEBUG == 1
8 |6 F* u7 A9 c5 a1 x' s        output = output_meta_result;
. ]7 M4 G5 T0 n$ N  B    end
0 U; P; X7 k; M: w+ ~else
    if DEBUG == 1
        rx_k1 = rx_tmp;
9 I: H) _+ V0 H+ q9 z% Y1 V        fval_k1 = fval_tmp;
  Y0 u8 P9 P. J! O8 T        flag_k1 = flag_tmp;
5 d9 t+ O3 U+ J" }1 k0 l: G  R: g        output = output_tmp;
    end
end
, a3 {# X6 ?+ {9 j, ~& V

Bojay

没人回答啊