穷举法举例

552433160

function l=fx000(n,m)
a=196.6;b=600;
! l+ _" N7 E$ nj=10;f=200;K=1000;d=3000;
: i$ P* y) w/ v$ l- Sq=1/11400;
* |) J. l! h) |+ b8 s3 z, @f2=0;f11=0;
/ Z9 I* |: Y+ k0 d; D' yf1=(j*m+K)*(1-normcdf(m*n,b,a));
! _8 j1 @+ w% z, t* Pfor k=1:1:m
    for i=1:1:n
    i1(i)=i*normpdf((k-1)*n+i,b,a);
4 H  C$ l5 f& Y0 L# H& F    f11=f11+i1(i);
    end
5 a0 `% c' e! b- P; u, o    f11;
0 T% X' Z0 l( l    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));
0 `- H" K# g5 g1 B    f2=f2+k1(k);
end
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;
! S  a6 [* ~4 b" ]f3=m*n-fn;
l=(f1+f2)/f3;
l=1000;
u=0;
for n=1:50
3 a6 @3 I7 V7 l+ e# ?    for m=1:50   
             min=fx000(n,m);
# g6 f% j) {; v& L3 N, i  `1 R            if l>min
                l=min;
2 ~+ [' h7 W# B; a! d3 i% n/ ~                u(1)=n;u(2)=m;
9 K" \7 _/ w. e3 b           
            end
    end
end
l
% b9 Z7 t. M6 J) U( U( \u

una

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

wanlun

好程序，值得一读

dukaibing

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

wangjian31791

支持啊   谢谢

何挺

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

wuhao2277

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

炎炎

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

天之叶子

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