关于debug的问题 。。。。

Bojay

这是02年CUMCU 彩票问题  debug看懂是什么意思  求大神指点！！
3 w7 y: J3 w+ E0 a: c
global v
v = 630589;                % 求解v为630589的收入水平情况
9 o9 P8 o2 r  G% C2 PDEBUG = 0;
rand('state',sum(100*clock))        % 初始化随机数生成器
format long g

4 h5 `2 y3 }( R6 ^%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% n! y! [' o' R; V, R: Z% 求解开始
% 对于K1型
p_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];
( P0 x# O  B- ?5 W0 Nbeq=1;
a_lb=[10,4,3,4,2];
4 Q8 G1 }# F4 d% B0 n: Rb_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);
, I" o1 F4 l1 r$ n/ U    0,0,0,0,1,-b_ub(5)];
" ?" p- @; D" ]. e; I: ub= [0;0;0;0];
lb=[0.5;0;0;0;0;0];
ub=[0.8;1;1;inf;inf;inf];
p_test = p_k1;
) m  C+ _1 ^% a# Prx0_tmp = zeros(6,1);
! M: @6 s: K$ |$ T7 a% Erx_meta_result = zeros(6,1);
$ h( S2 T- w$ V/ O0 u4 l* E% cfval_meta_result = inf;
flag_meta_result = nan; %用以判断有没有得到过可行解
if DEBUG == 1
7 g3 w7 m/ c, x" f8 i* q. a2 i    output_meta_result = [];
end
for j = 1:nums_test_of_initial_value
    %随机生成多个初始值rx0_tmp，以避免局部最优
    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
) a$ G! {  y4 c, P6 O+ w- Q+ `    rx0_tmp(4) = rand*1000;
: V6 N3 O8 ?6 t! x5 C    rx0_tmp(5) = rand*100;
    rx0_tmp(6) = rand*50;
    % 寻优
6 @5 G6 N7 ?7 Q8 ~    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
            fmincon('cpiao',rx0_tmp,A,b,...
$ {% n2 G9 U' N- w9 ^$ w$ r' u( a                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
& P# q1 Z! f0 \5 q: ]1 n+ x        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
        fval_meta_result = fval_tmp;
        rx_meta_result = rx_tmp;
+ A1 R4 [& i* X- y. w. `        flag_meta_result = 1;
        if DEBUG == 1
            output_meta_result = output_tmp;
        end
4 J# `+ `! ^' P6 p    end
end
( U# ?( f- [5 H* U% 把求得的最好结果保存下来
if ~isnan(flag_meta_result)
( v1 e3 v2 R; Y+ X! ]7 c    rx_k1 = rx_meta_result;
9 B  p" u6 M6 f    fval_k1 = fval_meta_result;
    flag_k1 = flag_meta_result;
' U2 e) g5 t# {9 d    if DEBUG == 1
        output = output_meta_result;
    end
' r+ a( \4 \/ |" D  b% U, Oelse
* l3 u4 ], e4 i    if DEBUG == 1
        rx_k1 = rx_tmp;
6 T. U% G# ?9 u4 i  k        fval_k1 = fval_tmp;
; S5 U9 u1 k% V: v8 z+ Y) h        flag_k1 = flag_tmp;
1 B; A' @0 R; M7 L8 T        output = output_tmp;
    end
end
9 V" g9 d, G  M1 j; g/ m

Bojay

没人回答啊