关于debug的问题 。。。。

Bojay

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

- W3 }5 u5 u9 N5 o7 ?global v
v = 630589;                % 求解v为630589的收入水平情况
DEBUG = 0;
rand('state',sum(100*clock))        % 初始化随机数生成器
  q$ G/ a3 W" K" {& e/ W8 e( s% kformat long g

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 求解开始
% q5 G; V9 l% j* D* k, F# T/ G% 对于K1型
p_k1 = [2e-7;8e-7;1.8e-5;2.61e-4;3.42e-3;4.1995e-2];
2 @6 |5 R( v+ e1 k3 |: z. z% 6个奖项6个变量
7 q- s6 h! B# @$ R# Y" S! ^Aeq=[1,1,1,0,0,0];
8 O) I1 Z: Z$ U- ~beq=1;
a_lb=[10,4,3,4,2];
* N* g1 _. j- S5 _3 ^b_ub=[233,54,17,20,10];
# _" F! T# F2 F- }1 x0 D- A3 {A= [0,0,0,-1,a_lb(4),0;
/ e& t* R0 V, g    0,0,0,1,-b_ub(4),0;
    0,0,0,0,-1,a_lb(5);
! h, B% T# S5 _: q# Q; ~) q) I; d    0,0,0,0,1,-b_ub(5)];
b= [0;0;0;0];
lb=[0.5;0;0;0;0;0];
4 d1 x7 p+ D- O5 E; Q+ w6 `$ iub=[0.8;1;1;inf;inf;inf];
p_test = p_k1;
rx0_tmp = zeros(6,1);
rx_meta_result = zeros(6,1);
6 P8 ?' \; ~" ?9 W+ Nfval_meta_result = inf;
/ z/ [( T  z$ k0 x+ S& Yflag_meta_result = nan; %用以判断有没有得到过可行解
if DEBUG == 1
; ?! H0 k8 _) Z8 C& x. R: x    output_meta_result = [];
end
  v0 `" D! y  ^  B& V4 P) Q$ Hfor j = 1:nums_test_of_initial_value
  c, t% ?8 i) d- F0 m) t5 ?    %随机生成多个初始值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);
    rx0_tmp(4) = rand*1000;
    rx0_tmp(5) = rand*100;
    rx0_tmp(6) = rand*50;
    % 寻优
$ m  m1 Y9 y8 O0 x! {5 }    [rx_tmp,fval_tmp,flag_tmp,output_tmp]= ...
            fmincon('cpiao',rx0_tmp,A,b,...
                    Aeq,beq,lb,ub,'nonlcon',[],1,p_test,a_lb,b_ub);
        % 上式倒数第四个参数是为了区分彩票的类型(K1/K2/K3/K4)
        % 最后三个是函数cpiao和nonlcon计算中可能要用到的量。
4 u- W* v5 K* F: z  w- h    if (flag_tmp == 1) && (fval_meta_result > fval_tmp)
        fval_meta_result = fval_tmp;
0 Y+ v3 A/ Q6 V, a& C        rx_meta_result = rx_tmp;
        flag_meta_result = 1;
% T' ~1 o% K* y$ s/ k        if DEBUG == 1
            output_meta_result = output_tmp;
! p4 S6 s& m. R. A: ?        end
    end
5 D+ G) c& `  W3 q2 j  X# wend
. d$ Z- U. K: a* O/ B; z% 把求得的最好结果保存下来
if ~isnan(flag_meta_result)
    rx_k1 = rx_meta_result;
$ _, T6 c& c3 H1 y$ P    fval_k1 = fval_meta_result;
6 ?$ a7 X9 w- K( L8 Q6 E) S: G    flag_k1 = flag_meta_result;
    if DEBUG == 1
        output = output_meta_result;
    end
else
    if DEBUG == 1
& f! s8 O& g* d; I1 N        rx_k1 = rx_tmp;
        fval_k1 = fval_tmp;
        flag_k1 = flag_tmp;
        output = output_tmp;
$ o8 y. o0 [; v    end
end
2 t) o7 V* u' {$ C  E( u: B

Bojay

没人回答啊