穷举法举例

552433160

function l=fx000(n,m)
a=196.6;b=600;
$ Z* D) w9 k6 y1 oj=10;f=200;K=1000;d=3000;
! Z0 N6 Z2 K0 w0 j5 l3 z* zq=1/11400;
& O- Q* J# |  Ef2=0;f11=0;
% l2 @0 g1 d! Qf1=(j*m+K)*(1-normcdf(m*n,b,a));
) @! j2 ^' m2 A% Y+ d: f5 q" u8 a: Xfor k=1:1:m
- @. T, H2 `, G) ~    for i=1:1:n
    i1(i)=i*normpdf((k-1)*n+i,b,a);
# ]- c& g: S! y/ F/ f1 C( n' N1 p: R    f11=f11+i1(i);
    end
- G) y5 z, x! S' P1 x9 ?* N    f11;
    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));
9 v) L  |2 h1 l! F% `' b    f2=f2+k1(k);
2 n0 H) y% t5 pend
7 u1 \7 C$ x: e0 F9 ]! mf2;

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;
0 i6 u7 O. X& H$ y! of3=m*n-fn;
8 S9 D% |2 q$ b: D$ Tl=(f1+f2)/f3;
l=1000;
, q3 t1 p. ^# G2 c* ?- f+ ^u=0;
' W2 ^' c2 j( W+ [9 b3 u8 ^for n=1:50
    for m=1:50   
             min=fx000(n,m);
            if l>min
( S' Q9 r+ n+ D" T                l=min;
                u(1)=n;u(2)=m;
+ z" e8 T$ [- I- ~           
) l1 K% e* H: p& k# @- d            end
    end
end
* R7 l0 W, s+ y$ ]* q4 K" k7 T: |l
u

una

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

wanlun

好程序，值得一读

dukaibing

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

wangjian31791

支持啊   谢谢

何挺

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

wuhao2277

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

炎炎

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

天之叶子

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