关于debug的问题 。。。。

Bojay

这是02年CUMCU 彩票问题  debug看懂是什么意思  求大神指点！！
( D# W( u/ {; @7 G" P- o
global v
v = 630589;                % 求解v为630589的收入水平情况
* z, \$ s9 P9 j1 RDEBUG = 0;
+ D7 o, W; ^0 X' urand('state',sum(100*clock))        % 初始化随机数生成器
0 S8 N2 {3 N4 Y4 g8 J: I% d4 b1 kformat long g
0 d4 w6 i* R. Z
9 P0 ?. U3 C: R: ^%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 求解开始
) Z: n8 r; n/ ^; V1 j: d; v% 对于K1型
p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
% 6个奖项6个变量
! o) b" P/ X" Y* A# ^Aeq=[1,1,1,0,0,0];
) a% R  _6 C2 }+ sbeq=1;
a_lb=[10,4,3,4,2];
& ?9 s4 ^( y# xb_ub=[233,54,17,20,10];
7 K5 t! \$ B- C1 G$ x6 t- N+ UA= [0,0,0,-1,a_lb(4),0;
" B; m2 d6 o3 B" o* w. I9 b    0,0,0,1,-b_ub(4),0;
' b7 P) g" B9 Y5 ~' y- w+ Z# R    0,0,0,0,-1,a_lb(5);
    0,0,0,0,1,-b_ub(5)];
; W  E- J; q9 }" ?" wb= [0;0;0;0];
lb=[0.5;0;0;0;0;0];
. B. T8 L) S% }9 v1 R: i6 Y) Sub=[0.8;1;1;inf;inf;inf];
p_test = p_k1;
rx0_tmp = zeros(6,1);
9 B) Z; P0 h4 V5 arx_meta_result = zeros(6,1);
* ?3 y! v+ d3 |' Y4 V$ jfval_meta_result = inf;
flag_meta_result = nan; %用以判断有没有得到过可行解
( o3 S. e3 b  t9 nif DEBUG == 1
    output_meta_result = [];
end
* `% u0 Y3 d! X7 u7 bfor j = 1:nums_test_of_initial_value
6 ?- Y1 }" D3 [  R3 m3 J    %随机生成多个初始值rx0_tmp，以避免局部最优
7 @9 w. _; M1 G7 S1 K7 O3 K    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
: ?4 p. O# c: ^) ?8 X% u( ^$ w    rx0_tmp(2) = rand*(1-rx0_tmp(1));
6 r: n4 b8 W% m% p/ Q2 y2 W    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
    rx0_tmp(4) = rand*1000;
    rx0_tmp(5) = rand*100;
! B- V' o* A! ]- T9 m2 N    rx0_tmp(6) = rand*50;
7 W8 e( w9 }5 N' Y* G6 t# Q    % 寻优
4 f" A  H- w9 R0 n8 u' E5 k$ F    [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)
  O3 i. V+ H. N% G8 }: Z        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
$ D( |, u6 D7 s3 O) f+ j# F    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
        fval_meta_result = fval_tmp;
; z6 N/ s8 f+ y) f$ O        rx_meta_result = rx_tmp;
2 c( h- T, a9 U* y4 \        flag_meta_result = 1;
8 S2 X% K8 |' P        if DEBUG == 1
" K" f7 T  ?# o. `' y            output_meta_result = output_tmp;
        end
* F6 s% V% }' z. z5 |) t    end
  i# X$ k, F; ^% Pend
/ Z/ T% h5 M3 n8 z% 把求得的最好结果保存下来
4 G( o+ \  ]5 P6 L4 W5 Xif ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
    fval_k1 = fval_meta_result;
    flag_k1 = flag_meta_result;
    if DEBUG == 1
+ k  k4 r$ G( E  k  ~        output = output_meta_result;
4 x: ^9 M" w7 J    end
else
& T- L* w2 ?: P4 R& r! r    if DEBUG == 1
* ~, v2 Q( y; C+ Z8 U        rx_k1 = rx_tmp;
        fval_k1 = fval_tmp;
6 P. T4 A5 I5 N+ W& B$ g1 u        flag_k1 = flag_tmp;
        output = output_tmp;
    end
end
* R9 l, S$ A8 v5 }

Bojay

没人回答啊