一种三参数Weibull分布极大似然估计的求解方法.pdfVIP

第27卷 第7期 河 南 科 学 Vo1．27 No．7 2009年 7月 HENAN SCIENCE Ju1．2009 文章编号：1004—3918(2009)07—0832—03 一 种三参数Weibul1分布极大似然估计的求解方法 史景钊 1，任学军 ， 陈新昌t， 李祥付 (1．河南农业大学机电工程学院，郑州 450002； 2．商丘职业技术学院，河南 商丘 476000) 摘 要：提出了一种求解极大似然估计的新算法．根据极大似然估计原理求出尺度参数的表达式，把该表达式代 入对数似然函数，使对数似然函数中只包含位置参数和形状参数，把求解非线性方程组问题变成了求解满足约束条 件的最优4Ef*q题，使问题得到了简化．该法具有计算精度高、运算速度快的优点，利用EXCEL的规划求解即可求 解，便于工程应用． 关键词：可靠性；威布尔分布；参数估什；极大似然估计 中图分类号：TB114．3 文献标识码：A 威布尔分布是瑞典物理学家W．Weibul1分析材料强度时在实际经验的基础上推导出来的分布形式[1】， 国内外大量研究表明，用三参数威布尔分布比用对数正态分布往往能更准确地描述结构疲劳寿命或腐蚀损 伤的概率分布[2]，物理意义更加合理；在 以损耗为特征的机械零件寿命评估中，采用三参数威布尔分布比采 用二参数威布尔分布拟合精度更高．因此，三参数威布尔分布在强度与环境研究领域及机械零件磨损寿命 评价中得到越来越广泛的应用． 、 威布尔分布参数估计方法有很多，国内外一直有人在进行相关研究[3-8]，但是三参数Weibul1分布的参数 估计比较复杂，人们也提出了很多方法，如有极大似然估计法[71、最大相关系数优化法恻、概率加权矩法[101、灰 色估计法 1【l】、图估计法、矩估计法r玎等，除图估计法外，其他方法大都计算复杂．图参数估计法虽有使用方便、 直观易懂、截尾寿命试验时也适用等特点，但这种仅靠 目力估计的方法，对作图质量要求较高，估计结果偏 差大】2【1，其他方法大都应用不便，只能通过Matlab或其他计算机语言编程计算，不利于一般工程技术人员的 应用 ． 极大似然估计是一种十分有效和通用的参数估计方法，在参数估计问题中占有重要地位．威布尔分布 的极大似然估计精度较高，但一般要用迭代法求解联立的3个超越方程，相当复杂．而且威布尔分布的极 大似然估计可能不存在或有多个解[134】．为此人们一直在探讨极大似然法的改进，提出各种更加合理简便 的方法．Qiao等先用Newton—Raphson法求二参数Weibul1分布的极大似然估计，在此基础上，逐步求解三 参数Weibul1分布的参数估计值堋【．Gove采用非线性最优化方法求解 ．曲延碌等提出一种叫做降阶的似 然方法求解方程幺且切【．王华胜等以似然函数值最大化为 目标，以一定步长寻找最佳位置参数，当位置参数确 定以后，似然方程组退化为二维方程组 l8【1．杨某存等提出的方法将三元方程组转化为二元方程组，先求解二 元方程组的解，然后将尺度参数和位置参数表示成形状参数的函数，使极大似然函数变成形状参数的单变 量函数，再估计出形状参数}l圳．这些方法不同程度地简化了极大似然估计法的求解问题，但仍不够理想． 本文根据极大似然估计原理先求出尺度参数的表达式，然后把该表达式代入对数似然函数，使对数似 然函数中只包含位置参数和形状参数，把求解三元非线性方程组问题变成了求解满足约束条件的极值问题， 使问题得到了简化，估计精度也得到了进一步提高． 1 威布尔分布的极大似然估计 威布尔分布的概率密度函数由下式给出 f(t)= ．exp(-( ， (1) 收稿 日期：2009—02—02 基金项 目：河南农业大学博士基金项 目 作者简介：史景钊 (1963一)，男，河南商丘人，副教授，研究方向为农业装备及其可靠性． 2009年7月 史景钊等：一种三参数Weibul1分布极大似然估计的求解方法