收稿日期: 2022-04-23
网络出版日期: 2022-07-13
基金资助
国家科技重大资助项目(J2019-V-0009-0103);中国博士后科学基金资助项目(2021T140098);中央高校基础研究项目(N2103009)
Three-Parameter Estimation of the Weibull Distribution Based on Least Squares Iteration
Received date: 2022-04-23
Online published: 2022-07-13
Supported by
the National Science and Technology Major Project(J2019-V-0009-0103);China Postdoctoral Science Foundation(2021T140098);Fundamental Research Funds for the Central Universities(N2103009)
威布尔分布因其数学处理的便利性和适应性而被广泛用于描述产品寿命,引入了位置参数的三参数威布尔分布是研究机械零部件可靠性最适合的模型之一,尤其适用于长寿命、高可靠性的产品。三参数威布尔分布的参数估计一直是关注的焦点,本文中提出了一种基于最小二乘的迭代方法对其进行参数估计,将初始位置参数设置为0,使用最小二乘法得到初始形状参数和尺度参数,将其代入位置参数的无偏估计得到新的位置参数,进行多次迭代,在此过程中形状参数和尺度参数逐渐变小,位置参数逐渐变大,最终获得稳定的形状参数、尺度参数和位置参数,即为最终的参数估计值,并计算可靠度为99%的寿命。通过蒙特卡洛(Monte Carlo)仿真证明此方法是收敛的,并在不同的威布尔模型、不同中小样本量(样本量为10、15、20、25和30)下,使用偏差(Bias)和均方差根(RMSE)两个指标与相关系数法进行对比,此方法估计的3个参数及可靠度为99%的寿命更加准确。通过两个实例分析,表明该方法具有可行性和有效性,估计结果与相关系数法相比更加保守,更适于工程应用。
杨小玉, 宋佳昕, 谢里阳, 等 . 基于最小二乘迭代威布尔分布的三参数估计[J]. 华南理工大学学报(自然科学版), 2023 , 51(2) : 20 -26 . DOI: 10.12141/j.issn.1000-565X.220231
The three-parameter Weibull distribution is widely used to describe product longevity because of the convenience and adaptability of its mathematical processing. The three-parameter Weibull distribution with location parameter is one of the most suitable models for studying the reliability of mechanical components, especially for long-life and high-reliability products. Parameter estimation of three-parameter Weibull distribution has always been the focus of attention. This paper proposed an iterative method based on least squares to estimate the parameters of the three-parameter Weibull distribution. The initial location parameter was set to 0, the initial shape parameter and scale parameter were obtained by using least squares, and the new location parameter was obtained by substituting them into the unbiased estimation of the location parameter, and multiple iterations were performed. In this process, the shape parameters and scale parameters gradually become smaller and the location parameters gradually become larger, and finally the stable shape parameters, scale parameters and location parameters were obtained, which are the final parameter estimates, and the lifetime of 99% reliability was calculated. The method was proved to be convergent by Monte Carlo simulation. Compared with the correlation coefficient method by two metrics including Bias and Root Mean Square Error (RMSE) for different Weibull models with different small and medium sample sizes (10, 15, 20, 25 and 30), the three estimated parameters and the 99% reliability of the lifetime of the proposed method are more accurate. The analysis of two examples shows that the method is feasible and valid. Compared with the correlation coefficient method, the estimation results are more conservative and more suitable for engineering application.
| 1 | SILVIA L P F, MICHEL F,DA S,Adjusted profile likelihoods for the Weibull shape parameter [J].Journal of Statistical Computation and Simulation. 2007,77(7):531-548. |
| 2 | 王晓峰,申桂香,张英芝,等 .可靠性模型参数估计方法的对比[J].华南理工大学学报(自然科学版),2011,39(6):47-52. |
| 2 | WANG Xiao-feng, SHEN Gui-xiang, ZHANG Ying-zhi,et al .Comparison of parameter estimation methods for reliability model [J].Journal of South China University of Technology(Natural Science Edition),2011,39(6):47-52. |
| 3 | 张英芝,翟粉莉,郑玉彬,等 .基于累积误差平方和最小的参数估计方法[J].华南理工大学学报(自然科学版),2020,48(11):49-54. |
| 3 | ZHANG Yingzhi, ZHAI Fenli, ZHENG Yubin,et al .Parameter estimation method based on minimum sum of squares of cumulative errors [J] Journal of South China University of Technology(Natural Science Edition),2020,48(11):49-54. |
| 4 | JACQUELIN J .Generalization of the method of maximum Likelihood [J].IEEE Transactions on Electrical Insulation,1993,28(1):65-72. |
| 5 | KANTAR Y M .Generalized least squares and weighted least squares estimation methods for distributional parameters [J].REVSTAT-Statistical Journal,2015,13(3):263-282. |
| 6 | DODSON B .The Weibull analysis handbook [M].Milwaukee,Wisconsin:ASQ Quality Press,2006. |
| 7 | HOSKING J R M .L-moments:analysis and estimation of distributions using linear combinations of order statistics [J].Journal of the Royal Statistical Society(Series B:Methodological),1990,52(1):105-124. |
| 8 | COHEN C A, WHITTEN B .Modified maximum likelihood and modified moment estimators for the three-parameter weibull distribution [J].Communications in Statistics,1982,13(1):47-68. |
| 9 | COUSINEAU D .Nearly unbiased estimators for the three-parameter Weibull distribution with greater efficiency than the iterative likelihood method [J].British Journal of Mathematical and Statistical Psychology,2009,62:167-191. |
| 10 | TEIMOURI M, HOSEINI S M, NADARAJA S .Comparison of estimation methods for the Weibull distribution [J].Statistics,2013,47(1):93-109. |
| 11 | CHENG R C H, AMIN N A K .Estimating parameters in continuous univariate distributions with a shifted origin [J].Journal of the Royal Statistical Society,1983,45(3):394-403. |
| 12 | 严晓东,马翔,郑荣跃,等 .三参数威布尔分布参数估计方法比较[J].宁波大学学报,2005,18(3):301-305. |
| 12 | YAN Xiao-dong, MA Xiang, ZHENG Rong-Yue,et al .Comparison of parameter estimation methods of three parameter Weibull distribution [J] Journal of Ningbo University,2005,18(3):301-305. |
| 13 | NAGATSUKA H, KAMAKURA T, BALAKRISHNAN N .A consistent method of estimation for the three-parameter Weibull distribution [J].Computational Statistics & Data Analysis,2013,58:210-226. |
| 14 | 许伟,程刚,黄林,等 .基于混沌模拟退火PSO算法的威布尔分布参数估计应用研究[J].振动与冲击,2017,36(12):134-139. |
| 14 | XU Wei, CHENG Gang, HUANG Lin,et al .A chaotic simulated POS algorithm application for Weibull distribution parameter estimation [J].Journal of Vibration and Shock,2017,36(12):134-139. |
| 15 | 高镇同 .航空金属材料疲劳性能手册[R].北京:北京航空材料研究所,1981. |
/
| 〈 |
|
〉 |