关于debug的问题 。。。。

Bojay

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

9 |8 d8 k4 e0 Y8 B- Aglobal v
' q  z7 x+ j( |( s* p& I  Mv = 630589;                % 求解v为630589的收入水平情况
& D6 w( A- X4 _% R/ T+ l: ZDEBUG = 0;
rand('state',sum(100*clock))        % 初始化随机数生成器
format long g
, g3 G/ b& i8 F" U# T7 ~: X
) J. ~- f% Y6 [. q, H" a0 z7 f%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 }* [4 \7 z. ]5 c+ y) R8 x3 y% 求解开始
% 对于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];
beq=1;
a_lb=[10,4,3,4,2];
5 g* R) T4 G' h  x( tb_ub=[233,54,17,20,10];
# F; A9 @. n" n, a/ ~3 U2 x- tA= [0,0,0,-1,a_lb(4),0;
* I6 i, u  H: q' e' r4 d: x    0,0,0,1,-b_ub(4),0;
- A4 M- L8 A; |  D1 {) `8 X! x    0,0,0,0,-1,a_lb(5);
; J5 Q. X4 H8 [6 {" E( y5 }+ J    0,0,0,0,1,-b_ub(5)];
+ e' k3 l* G8 D. cb= [0;0;0;0];
lb=[0.5;0;0;0;0;0];
: z' g8 U  m# P  xub=[0.8;1;1;inf;inf;inf];
p_test = p_k1;
rx0_tmp = zeros(6,1);
rx_meta_result = zeros(6,1);
fval_meta_result = inf;
. x! U9 r% S2 {1 ]flag_meta_result = nan; %用以判断有没有得到过可行解
if DEBUG == 1
    output_meta_result = [];
7 J5 }5 q% @- ~0 A1 V# vend
for j = 1:nums_test_of_initial_value
, X6 ]! A7 I5 W5 x. z    %随机生成多个初始值rx0_tmp，以避免局部最优
/ U" I' R- M1 ?3 {2 C% p) R! r    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
    rx0_tmp(2) = rand*(1-rx0_tmp(1));
% V0 L& K4 Z; o& H# A) _6 _    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
3 s, P* X; A. z: S, y# B    rx0_tmp(4) = rand*1000;
, Q# `9 A+ h; Q9 [* b% s+ ^: ^9 G    rx0_tmp(5) = rand*100;
) H+ a0 o& h& D  E    rx0_tmp(6) = rand*50;
    % 寻优
    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
0 H% J4 b6 r! g            fmincon('cpiao',rx0_tmp,A,b,...
                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
& d2 x, n4 J8 `; r* G' A! H        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
5 p" n( D) j1 P5 _9 F        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
" {% a, B' i+ f        fval_meta_result = fval_tmp;
& z/ g) B- ?; k. x+ c        rx_meta_result = rx_tmp;
        flag_meta_result = 1;
        if DEBUG == 1
            output_meta_result = output_tmp;
        end
/ |. y6 j- P; S; {# ]    end
, r: \" V0 v  l* ^& D! Pend
! H( H' n$ i4 w& c& n% 把求得的最好结果保存下来
6 X7 A. P+ M7 l9 p) N6 G4 Mif ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
    fval_k1 = fval_meta_result;
    flag_k1 = flag_meta_result;
% ?' C. r2 d: i" D    if DEBUG == 1
        output = output_meta_result;
8 u& |8 T1 A2 b    end
+ G/ w0 L6 }' x# G/ [& \else
% E* t: G' V$ @4 D5 F5 z: J: _    if DEBUG == 1
5 F1 q3 s4 [6 G0 U1 }        rx_k1 = rx_tmp;
        fval_k1 = fval_tmp;
        flag_k1 = flag_tmp;
5 B6 M; M4 k* d4 A  l! S0 x9 h        output = output_tmp;
9 I  k( z" W1 W2 r- D$ E9 c  o' G    end
end

Bojay

没人回答啊