关于debug的问题 。。。。

Bojay

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

& [0 n1 x# n8 H  b3 t$ p, E4 C8 }; t6 nglobal v
( H* L3 J, [/ _4 hv = 630589;                % 求解v为630589的收入水平情况
5 }0 b% F$ Y( {1 [7 f; ?DEBUG = 0;
  U. B; Q/ ]+ p6 |( Grand('state',sum(100*clock))        % 初始化随机数生成器
format long g

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
1 z# X. f1 \! A# \% 求解开始
% 对于K1型
; n; m' {; p7 ep_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];
8 h' C! K' M" x1 |" f6 Bbeq=1;
a_lb=[10,4,3,4,2];
b_ub=[233,54,17,20,10];
A= [0,0,0,-1,a_lb(4),0;
' L; c, W2 i; w3 h    0,0,0,1,-b_ub(4),0;
; u( p6 u) C8 r: H' d4 V( K3 f6 t4 T    0,0,0,0,-1,a_lb(5);
2 Y8 k6 J% ]4 @4 {    0,0,0,0,1,-b_ub(5)];
b= [0;0;0;0];
* w3 b9 _+ W4 b' Nlb=[0.5;0;0;0;0;0];
! I( w4 B! m* _: hub=[0.8;1;1;inf;inf;inf];
1 H* _. v- g8 i$ r3 }# h' i$ M6 Ip_test = p_k1;
7 X% k/ v6 |9 n$ \; a/ L5 U' Grx0_tmp = zeros(6,1);
+ G& \, @2 A, `6 c- U3 }7 frx_meta_result = zeros(6,1);
fval_meta_result = inf;
: N" M$ H* w6 Xflag_meta_result = nan; %用以判断有没有得到过可行解
if DEBUG == 1
( Q4 r+ G: _6 N! f. L, M& l    output_meta_result = [];
end
for j = 1:nums_test_of_initial_value
    %随机生成多个初始值rx0_tmp，以避免局部最优
% \* c, w6 z5 q3 H) Q    rx0_tmp(1) = rand*(0.8-0.5) + 0.5;
1 e% f. r5 _2 {1 V2 j    rx0_tmp(2) = rand*(1-rx0_tmp(1));
" [4 o3 z9 w% h  C$ E: `7 o    rx0_tmp(3) = 1 - rx0_tmp(1) - rx0_tmp(2);
8 t; Q9 V2 V! `! `4 s    rx0_tmp(4) = rand*1000;
    rx0_tmp(5) = rand*100;
    rx0_tmp(6) = rand*50;
    % 寻优
    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
, Q" C9 m) z) X4 z( ~* k            fmincon('cpiao',rx0_tmp,A,b,...
1 E) Q' v1 Z8 Z3 g" N                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
$ j1 L7 W# Y& Y, b# E" }+ O* A        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
0 B$ X& E3 ?2 f$ D7 T5 r1 |        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
5 i, s% `) k/ u* Z# D    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
        fval_meta_result = fval_tmp;
: L+ y* a- l2 Q' x        rx_meta_result = rx_tmp;
1 P  x0 [  z$ |        flag_meta_result = 1;
        if DEBUG == 1
$ e1 E  l: [" y7 i0 Q            output_meta_result = output_tmp;
        end
    end
end
2 @! H$ H" V3 I7 Q; g* C% 把求得的最好结果保存下来
- ]2 p5 `/ Q0 A7 X, Q, r, Fif ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
    fval_k1 = fval_meta_result;
( j1 ?- o& `# Z$ j# W* `) T    flag_k1 = flag_meta_result;
) |" w9 c* z" q/ G" p0 @% g    if DEBUG == 1
        output = output_meta_result;
    end
else
    if DEBUG == 1
( n' \9 M: C2 c  q; w* A& B        rx_k1 = rx_tmp;
' e* J3 J" u0 K) ?: U. ]        fval_k1 = fval_tmp;
        flag_k1 = flag_tmp;
        output = output_tmp;
8 ^( E# a. L$ e1 X3 F    end
end
$ b% G: f! c( u) e6 |2 E

Bojay

没人回答啊