24小时热门版块排行榜

返回列表

本帖产生 1 个博学EPI ，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

shuiliushui

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 807
帖子: 695
在线: 232.7小时
虫号: 271855

[交流] 正态分布图的绘制，横坐标为一范围，高手帮忙

数据如下，要求表示出多少范围占到多少的意思，非常感谢！！！
粒径μm          质量百分含量%       粒径μm                质量百分含量%
850μm以上          5             830μm以上                      0.5
850μm~600μm    45             830μm~550μm             25.1
600~300μm          40             550~380μm                   31.6
300~150μm          8                380~270μm                   12.5
150μm以下             2                270~250μm                12.7
                                       250~180μm                8.8
                                       180~150μm          4.7
                                       <150μm                   4

回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-02-19 22:02:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

leedobb

金虫 (正式写手)

博学EPI: 3
应助: 77 (初中生)
金币: 152.4
帖子: 856
在线: 270.4小时
虫号: 1199535

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