我用matlab编了好久没编出来，有没有人变出来了，希望指教：
N=3200
R = lognrnd (2,2,[1,N]);
target=0;
N1=32;
for i=1:100
    r1=R((i-1)*N1+1:(i-1)*N1+N1);
    mu=1/N1*sum(log(r1));
    sigma=(1/(N1-1))*sum((log(r1)-mu).^2);
    F=normcdf(((log(r1)-mu)/sqrt(sigma)),0,1);
   F1=logncdf(r1,2,2).';
    F=F.';
    j=1:N1;
    S = -sum((((2*j)-1)/N1)*(log(F1)+log(1-F1(N1+1-j))))-N1;
            if S<0.852
                target=target+1;
            end
end


S是我自己算出来的，好大呀，参考值才0.8，我哪里错了，有高手有源程序吗，651903249@qq.com，谢谢