穷举法举例

552433160

function l=fx000(n,m)
a=196.6;b=600;
j=10;f=200;K=1000;d=3000;
q=1/11400;
f2=0;f11=0;
( w1 P% t9 n5 R1 J" f2 Sf1=(j*m+K)*(1-normcdf(m*n,b,a));
for k=1:1:m
+ R; c* K! j0 T1 p( r3 I; Z8 R7 f    for i=1:1:n
" c2 k& ~/ b' {& R. [4 W    i1(i)=i*normpdf((k-1)*n+i,b,a);
  n$ _6 {0 \) E. M    f11=f11+i1(i);
5 r  ~' M' W5 m) u( m& d    end
: j2 w, Y3 o/ M) ?. V* |    f11;
1 Y; M0 ?; W% s5 Q& W- f    k1(k)=(j*k+d+f11/(normcdf(k*n,b,a)-normcdf((k-1)*n,b,a)))*(normcdf(k*n,b,a)-normcdf((k-1)*n,b,a));
    f2=f2+k1(k);
end
/ N6 C: O, S& w& Z1 |f2;

fn=1125899906842624/5644425081792261*2^(1/2)*erf(1/2*2^(1/2)*(n*m-b)/a)*pi^(1/2)*n*m-1125899906842624/5644425081792261*2^(1/2)*erf(1/2*2^(1/2)*(n*m-b)/a)*pi^(1/2)*b+2251799813685248/5644425081792261*exp(-1/2*(n*m-b)^2/a^2)*a+1125899906842624/5644425081792261*erf(1/2*b/a*2^(1/2))*2^(1/2)*pi^(1/2)*n*m-1125899906842624/5644425081792261*erf(1/2*b/a*2^(1/2))*2^(1/2)*pi^(1/2)*b-2251799813685248/5644425081792261*exp(-1/2*b^2/a^2)*a;
& F8 v" ~6 ]' A8 d8 `! c% F0 mf3=m*n-fn;
' v6 s( e; Y; P8 {8 h) @; Xl=(f1+f2)/f3;
: i7 T; Y; d' wl=1000;
5 I! p5 Z6 I# K/ i: U# Au=0;
for n=1:50
    for m=1:50   
             min=fx000(n,m);
            if l>min
                l=min;
  i) q7 M% T0 P3 W  b6 d: y                u(1)=n;u(2)=m;
5 @& ~$ \- d$ Z) ^/ j           
9 Q% w& J1 p& x2 N            end
    end
. ^" {% Q, o5 [( ~7 N+ F& Uend
l
' \" L1 T% \* u. M" tu

una

我看不懂呢啊。。。。我要学习

wanlun

好程序，值得一读

dukaibing

1# 552433160    要是能有注释就好了，老大，注释一下吧。。。。3q

wangjian31791

支持啊   谢谢

何挺

????????????????????????????

wuhao2277

顶~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

炎炎

简练啊，好程序，拿去看看~~~~~~~

天之叶子

顶~~~~~~~~~~~~~~~~~~~~~~~~