我用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,谢谢


雷达卡


京公网安备 11010802022788号







