我用matlab拟合了一下,得到最好结果是:
y= 1/(2*pi)^0.5/sig0 *exp (- (x-x0)^2/2/sig0^2)
x0=626, sig=180 (第一组)对应的百分比为: 8.87, 44.86,42.87, 3.05 0.36
x0=435,sig=210 (第二组)对应的百分比为:2.95 26.12 31.19 18.12 2.69 7.70 2.50 8.72
这已经是最好的结果了。
源程序为:
x=[830,550,380,270,250,180,150];
pr=[0.05, 0.251,0.316,0.125,0.127,0.088,0.047,0.04];
dpp0=10;
for x0=210:0.5:630
for sig=10:0.5:590
xx= (x -x0)/sig;
p=normcdf(xx);
pp(1)=1-p(1);
for k=1:6
pp(k+1)=p(k)-p(k+1);
end
pp(8)=p(7);
dpp= sum(pp.*(pr-pp).^2);
if(dpp
dpp0=dpp;
xx0=x0;
sig0=sig;
pp0=pp;
end
end
end
X=[0:1200];
Y=pdf('norm',X,xx0,sig0);
plot(X,Y)
xlabel('粒径(微米)');
ylabel('分布');
xx0,sig0,