利用ECE的规划求解进行求解威布尔分布参数讲课稿.pptVIP

利用EXCE的规划求解进行求解威布尔分布参数;采用三参数威布尔分布比采用二参数威布尔分布拟合精度更高。因此，三参数威布尔分布在强度与环境研究领域及机械零件磨损寿命评价中得到越来越广泛的;应用。

在农业机械的强度设计中也经常要用到威布尔分布。

威布尔分布参数估计方法有很多，国内外一直有人在进行相关研究[3-8;]，现有几十种参数估计方法，但多数只能用于形状参数和尺度参数的估计。在众多的估计方法中，能用于三参数估计的并不多，见诸文献的有极大似然估;计法、最大相关系数优化法、概率权重矩法、灰色估计法、图估计法等，除图估计法外，其他方法大都计算复杂，应用不便，即便是计算机水平发达的今天;，也只能通过Matlab或其他计算机语言编程计算。EXCEL提供了超强的数学运算、统计分析等实用程序，利用它的规划求解功能可以快速、;高效地求解三参数威布尔分布的参数估计问题。

2.三参数威布尔分布模型威布尔分布的寿命分布函数由下式给出式中：m称为形状参数;，m0；eta;称为尺度参数，eta;0；gamma;称为位置参数，也称最小寿命，表示产品在gamma;以前不会失效，对;于产品寿命有gamma;≥0,gamma;=0时退化为二参数威布尔分布；t是产品的工作时间，t≥gamma;。;当m1时，由式(3)给出的失效率是递减型的，适合于建模早期失效；当m=1时，失效率为常数，即退化为指数分布，适合于建;模随机失效；当m1时，失效率是递增的，适合于建模磨耗或老化失效。

设有n个产品进行寿命试验数据，按失效时间先后得到的寿命数据;失效时间(顺序统计量)为nt≤t≤Lambda;≤t12，对应的累计失效概率(经验分布函数)为()()();12nFt≤Ft≤Lambda;≤Ft。

其中到第i个产品失效时的累计失效概率F(ti)可用中位秩算;法求得：

Ftii（2）根据失效时间和累计失效概率即可用各种方法对其参数进行估计。

3.最大相关系数优化法对（;1）式做变形处理，并取两次自然对数得到：4.用EXCEL进行参数估计（6）式所表示的方程十分复杂，解该方程一般是通过编程，用??值解法;求出gamma;，然后求再用最小二乘法或其他方法求解形状参数和尺度参数。MSEXCEL具有强大的统计和计算功能，其“规划求解”功;能更是求解最优化问题的强有力工具，（6）式所表示的方程利用EXCEL的“规划求解”功能可很容易解出，然后再利用其散点图的趋势线功能即可求;出形状参数和尺度参数。本文通过实例，就相关系数优化法，用EXCEL进行求解。

例：选取5台某产品进行可靠性试验，失效时间分别是;27,32,36,42,49，已知产品寿命服从威布尔分布，试估计分布参数。

1)准备数据表按图1准备数据表⑴在A2～A6;单元格中输入产品失效的顺序号1～5；⑵在B2～B6单元格中输入产品的失效时间27、32、36、42、49；

⑶在I8单元;格中输入位置参数gamma;的迭代初值，初值可选择接近于第一个失效时间，也可用图估计法的估计值作为初值；⑷在C2单元格中输入公式;“=LN(B2-$I$8)”，用填充柄填充C3～C6单元格，C2～C6单元格的值为ix，即ln(minus;gamma;);it；⑸在D2单元格中输入公式“=C2*C2”，用填充柄填充D3～D6单元格，D2～D6单元格的值为为2ix；⑹在;E2单元格中输入公式“=(A2-0.3)/5.4”，用填充柄填充E3～E6单元格，E2～E6单元格的值为为()iFt，这里;()iFt采用中位值算法，即F(t)=(iminus;0.3)(n+0.4)i；⑺在F2单元格中输入公式“;=LN(LN(1/(1-E2)))”，用填充柄填充F3～F6单元格，F2～F6单元格的值为为iy，即1()lnln1;iminus;Ft；⑻在G2单元格中输入公式“=C2*F2”，用填充柄填充G3～G6单元格，G2～G6单元格的值为i;ixsdot;y；⑼在C7单元格中输入公式“=AVERAGE(C2:C6)”，C7单元格的值为x；⑽在C8单元格中;输入公式“=SUM(D2:D6)”，C8单元格的值为Sigma;=niix`2；⑾在C9单元格中输入公式“=AVER;AGE(F2:F6)”，C9单元格的值为y；⑿在C10单元格中输入公式“=SUM(G2: