关于debug的问题 。。。。

Bojay

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

+ A+ l: d/ y( x# p- _; }, Bglobal v
3 J" F/ _: n2 v3 zv = 630589;                % 求解v为630589的收入水平情况
DEBUG = 0;
rand('state',sum(100*clock))        % 初始化随机数生成器
3 M4 K! F" Z) Gformat long g

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 求解开始
, w& P! P( s. O% 对于K1型
5 Q' \+ d' I- W! o& y, @p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
' m/ |4 v( w8 i3 a; T% 6个奖项6个变量
Aeq=[1,1,1,0,0,0];
5 x& F' a; e; B/ w5 w* [- ?beq=1;
4 g& q# v2 _: ra_lb=[10,4,3,4,2];
$ |) f( _  x+ T" W) \3 Mb_ub=[233,54,17,20,10];
A= [0,0,0,-1,a_lb(4),0;
: v: L/ _7 g5 p8 R$ r8 C    0,0,0,1,-b_ub(4),0;
8 V7 N9 K6 |2 @  E+ ~: @0 M8 f1 S. @    0,0,0,0,-1,a_lb(5);
    0,0,0,0,1,-b_ub(5)];
b= [0;0;0;0];
lb=[0.5;0;0;0;0;0];
ub=[0.8;1;1;inf;inf;inf];
p_test = p_k1;
6 N, E  V" p$ p, frx0_tmp = zeros(6,1);
rx_meta_result = zeros(6,1);
& E3 m" Q$ V' k; O% vfval_meta_result = inf;
: ?8 c% x9 P3 \, j7 uflag_meta_result = nan; %用以判断有没有得到过可行解
if DEBUG == 1
    output_meta_result = [];
1 \! m5 N8 ^; _. g" p7 _end
for j = 1:nums_test_of_initial_value
" F; g) S  x, P$ o    %随机生成多个初始值rx0_tmp，以避免局部最优
- ]! l8 h$ E8 ?" \    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
3 S+ {3 F% ^/ \    rx0_tmp(2) = rand*(1-rx0_tmp(1));
    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
0 K. ~' `, I: }2 i* |5 t% _: P+ Y6 t    rx0_tmp(4) = rand*1000;
: r" }! c" f5 q0 {/ n/ }7 H1 K* ]    rx0_tmp(5) = rand*100;
    rx0_tmp(6) = rand*50;
    % 寻优
% z% J+ B( J* Y8 s' U    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
            fmincon('cpiao',rx0_tmp,A,b,...
* S# b) n1 s! X                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
# m5 C& d% Z3 y4 l3 X2 ?        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
4 c% c5 o, e% f6 \8 H; e7 }    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
9 s1 L0 @2 \" Q' d* o        fval_meta_result = fval_tmp;
        rx_meta_result = rx_tmp;
        flag_meta_result = 1;
        if DEBUG == 1
3 J! c, r: p) G- ~+ L) i            output_meta_result = output_tmp;
! X8 }9 `! I4 ^0 v        end
    end
" D3 G  t( E. _: t5 w8 O; _2 Qend
- E, m5 M. I* p! W; q0 `% 把求得的最好结果保存下来
if ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
    fval_k1 = fval_meta_result;
    flag_k1 = flag_meta_result;
    if DEBUG == 1
) V, l4 y; q0 P7 w6 f/ p7 ~        output = output_meta_result;
    end
else
    if DEBUG == 1
        rx_k1 = rx_tmp;
        fval_k1 = fval_tmp;
* A$ E% m, e& L9 _% G! E        flag_k1 = flag_tmp;
        output = output_tmp;
& i5 n) v/ y" }# i9 c, y% u& N    end
0 n) q# a9 v7 H) P) d* q- s# Rend
2 ]0 U% }! ?3 c1 h- G7 }  _

Bojay

没人回答啊