关于debug的问题 。。。。

Bojay

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

global v
v = 630589;                % 求解v为630589的收入水平情况
DEBUG = 0;
rand('state',sum(100*clock))        % 初始化随机数生成器
format long g

9 B. I4 O1 ^1 ^: z. o9 [%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% L7 T$ k- n& _% 求解开始
% 对于K1型
* V% ~* ]" r4 ^3 n% F& rp_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
% 6个奖项6个变量
Aeq=[1,1,1,0,0,0];
) {% F0 V' x+ `5 V  r- z$ \beq=1;
a_lb=[10,4,3,4,2];
b_ub=[233,54,17,20,10];
, z2 V( f+ S, S" n( p# D% EA= [0,0,0,-1,a_lb(4),0;
    0,0,0,1,-b_ub(4),0;
" I5 ^3 E  w$ _* V5 N    0,0,0,0,-1,a_lb(5);
  _- }& r- P% K! @$ r3 |1 ^  u, W    0,0,0,0,1,-b_ub(5)];
4 ^" }  l+ @' T* i% zb= [0;0;0;0];
lb=[0.5;0;0;0;0;0];
, S$ k" }2 l, w- e6 P3 `ub=[0.8;1;1;inf;inf;inf];
. \6 v' W. x2 j' I' qp_test = p_k1;
rx0_tmp = zeros(6,1);
3 q" S- }* d4 @+ b) O! C& Prx_meta_result = zeros(6,1);
3 x( v% n+ `+ B1 `fval_meta_result = inf;
% X1 Z! f8 ]+ C2 n# [flag_meta_result = nan; %用以判断有没有得到过可行解
; S2 g" R( u8 I8 @/ aif DEBUG == 1
    output_meta_result = [];
/ Y1 S. k' Z  ]9 u; yend
for j = 1:nums_test_of_initial_value
    %随机生成多个初始值rx0_tmp，以避免局部最优
    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
' W" i3 q  [: S+ `$ J5 x    rx0_tmp(2) = rand*(1-rx0_tmp(1));
# a: ^- O' T  I- Y    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
6 D: Q( g3 I" f    rx0_tmp(4) = rand*1000;
    rx0_tmp(5) = rand*100;
4 K: `- ?- W* K, s    rx0_tmp(6) = rand*50;
    % 寻优
; \1 E6 v5 F3 @2 \$ ^  x# Z    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
. f2 v  B1 d, \/ F( U            fmincon('cpiao',rx0_tmp,A,b,...
- N- h* X- g4 ?: ~0 F$ j0 X  c/ {$ l                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
$ o" ~* A# r7 Q8 B5 r, F% e8 d        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
6 |! b8 A# O+ }7 ?! Y* u3 u' v* C    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
: R7 C' B2 \. S+ L4 D' r" K        fval_meta_result = fval_tmp;
        rx_meta_result = rx_tmp;
% x( J# h5 a* m        flag_meta_result = 1;
        if DEBUG == 1
            output_meta_result = output_tmp;
        end
0 [. {4 w$ r9 A+ X' U9 h    end
- C0 s1 A; g: s  k* J# b3 r: h, dend
( f/ t4 v4 b  c+ ^7 r( o% 把求得的最好结果保存下来
if ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
: {+ k% |0 y7 @- N- N8 e5 w    fval_k1 = fval_meta_result;
    flag_k1 = flag_meta_result;
) p) U6 f2 J. _; \; c4 A$ T    if DEBUG == 1
* Y- N5 m9 n  T! r' ~2 L        output = output_meta_result;
    end
else
: {5 u+ Q3 x: p# V7 {# j9 ?    if DEBUG == 1
, Y/ ?$ }$ J6 e- C' ^        rx_k1 = rx_tmp;
        fval_k1 = fval_tmp;
        flag_k1 = flag_tmp;
3 X. k' {/ t+ j+ f/ ~& A        output = output_tmp;
1 o; V1 r3 o) X2 E/ a    end
end

Bojay

没人回答啊