一种不固定形状参数的威布尔分布参数估计迭代方法

本发明属于三参数威布尔分布参数估计方法，涉及一种不固定形状参数的威布尔分布参数估计迭代方法。

背景技术：

1、威布尔分布成为近年来研究寿命概率最为广泛的参数模型。威布尔分布是由瑞典科学家weibull提出。他按照最弱链模型推算出该分布。威布尔分布因其数学的便利性，尤其是通过累积分布取两次对数变换将其线性化，使得参数拟合较为方便。很多学者更多专注的是两参数威布尔分布的应用和参数估计，但是对于长寿命、高可靠性的产品，三参数威布尔分布更能够逼近实际，更富于弹性，以及更能够充分地描述失效机理然。

2、三参数和两参数威布尔分布的区别在于位置参数，它沿水平坐标移动。然而，包含位置参数会导致参数估计的有效推断尤为困难。当形状参数小于2或等于2时，不满足极大似然估计(mle)渐近性的经典正则性条件，所以此时极大似然方法不适合估计三参数威布尔分布的参数。常用的参数估计方法还有相关系数法、最小偏差法等。

3、当样本量较小时，经典的统计分析可能导致估计结果的不稳定性或不准确性。因此，对于样本较小时，发明一种更精确、更稳健、更有效的三参数威布尔分布数估计方法具有重要意义。

技术实现思路

1、为解决上述技术问题，本发明通过建立形状参数、尺度参数、位置参数与其他参数之间的关系，提出了一种新的三参数威布尔分布的参数估计的方法。

2、本发明的一种不固定形状参数的威布尔分布参数估计迭代方法，包括：

3、步骤1：确定形状参数与位置参数、样本量、样本标准差及样本与对应秩的相关系数之间的关系式；

4、步骤2：根据极大似然法确定尺度参数与样本、位置参数、形状参数之间的关系式；

5、步骤3：根据三参数威布尔分布最小次序统计量的期望，确定位置参数的无偏估计；

6、步骤4：设位置参数的初始估计值为0，代入步骤1的关系式中得到形状参数的初始估计值并代入步骤2的关系式中获得尺度参数的初始估计值，再将形状参数的初始估计值和尺度参数的初始估计值代入步骤3获得新的位置参数的估计值，如此多次迭代最终获得稳定的形状参数、位置参数和尺度参数。

7、在本发明的不固定形状参数的威布尔分布参数估计迭代方法中，所述步骤1确定形状参数与位置参数、样本量、样本标准差及样本和对应秩的相关系数之间的关系式的推导过程为：

8、(1)服从三参数威布尔分布w(β,η,γ)的随机变量t概率密度函数为：

9、

10、式中，β,η,γ其中分布是威布尔分布的形状参数、尺度参数和位置参数，满足γ＞0,β＞0,η＞0,t≥γ；ti(i＝1,2,…,n)为样本，n为样本量；

11、(2)三参数威布尔分布的概率分布函数为：

12、

13、(3)样本ti与其对应的秩ri之间的相关系数ρ为：

14、

15、式中，corr是相关系数，μ为期望μ＝e(t)，σ2为总体方差σ2＝var(t)；

16、(4)根据式(1)和(2)，可得：

17、

18、令t-γ＝x则：

19、

20、

21、

22、所以得到：

23、

24、将式(8)带入式(3)得到：

25、

26、样本变异系数将样本标准差s代替上述的总体标准差σ，则确定形状参数与位置参数、样本量、样本标准差及样本与对应秩的相关系数之间的关系：

27、

28、在本发明的不固定形状参数的威布尔分布参数估计迭代方法中，所述步骤2根据极大似然法确定尺度参数与样本、位置参数、形状参数之间的关系具体为：

29、(1)极大似然法主要是求似然函数的极大值等效为求对数似然函数的极大值，三参数威布尔分布参数对数似然函数的表达式为：

30、

31、(2)求解最大值，使每个参数的偏导等于0，这里对位置参数求偏导为0，则：

32、

33、化简得到：

34、

35、在本发明的不固定形状参数的威布尔分布参数估计迭代方法中，所述步骤3根据三参数威布尔分布最小次序统计量的期望，确定位置参数的无偏估计具体为：

36、(1)根据f(t(1,n))＝n[1-f(t)]n-1f(t)，三参数威布尔分布最小样本值t(1,n)的概率密度函数为：

37、

38、(2)三参数威布尔分布最小样本值t(1,n)的期望为：

39、

40、(3)位置参数的无偏估计为：

41、

42、其中，n为样本量，其中γ(x)为伽马函数。

43、在本发明的不固定形状参数的威布尔分布参数估计迭代方法中，所述步骤4具体为：

44、步骤4.1：用monte carlo随机抽取三参数威布尔分布w(β,η,γ)，获得n个样本，分别为t1,t2,…,tn；

45、步骤4.2：设置迭代次序n取1000，两次迭代结果之差的绝对值取10-5；

46、步骤4.3：将位置参数的初始估计值设置为0，即：

47、步骤4.4：计算样本变异系数cv、样本标准差s、样本和对应秩的相关系数ρ，将样本量n和位置参数的初始估计值代入公式(10)中，计算出形状参数的初始估计值

48、步骤4.5：将形状参数的初始估计值位置参数的初始估计值和样本值代入公式(13)，得到尺度参数的初始估计值

49、步骤4.6：将形状参数的初始估计值和尺度参数的初始估计值和样本最小t(1,n)值代入公式(16)中，可计算出新的位置参数的估计值若结束迭代过程，则否则继续进行下述步骤；

50、步骤4.7：将新的位置参数的估计值代入公式(10)中，可计算出新的形状参数估计值

51、步骤4.8：将新的位置参数的估计值和新的形状参数的估计值代入公式(13)中，可计算出尺度参数的估计值

52、重复步骤4.6-4.8，计算出本次的位置参数的估计值、形状参数的估计值和尺度参数的估计值，与上次计算的位置参数的估计值、形状参数的估计值、尺度参数的估计值分别进行对比，若三个参数的两次迭代结果之差的绝对值都小于10-5，则停止计算，本次结果即为最终估计值，否则继续迭代计算1000次，为最终结果。

53、在本发明的不固定形状参数的威布尔分布参数估计迭代方法中，采用不固定形状参数迭代方法进行三参数威布尔分布参数估计，多次迭代后位置参数逐渐增大到稳定值即为位置参数的估计值，尺度参数逐渐减小到稳定值即为尺度参数的估计值，形状参数逐渐减小到稳定值即为形状参数的估计值。

54、在本发明的不固定形状参数的威布尔分布参数估计迭代方法中，样本量为10-30个。

55、本发明的一种不固定形状参数的威布尔分布参数估计迭代方法，采用不固定形状参数迭代方法进行三参数威布尔分布参数估计，位置参数逐渐增大到稳定值即位置参数估计值，尺度参数逐渐减小到稳定值即尺度参数估计值，形状参数逐渐减小到稳定值即为形状参数估计值，迭代过程简单且快速收敛，可有效估计出威布尔分布的三个参数。在小样本估计时，此不固定形状参数的迭代方法与传统的极大似然方法、矩法及相关系数法相比更加准确、更加稳定，更适用于工程应用。