基于Matlab的随机截尾数据下的Weibull分布参数估计.pdfVIP

第28 卷 第5 期 河 南 科 学 Vol.28 No.5 2010 年5 月 HENAN SCIENCE May 2010 文章编号：1004-3918 （2010）05-0584-04 基于Matlab 的随机截尾数据下的 Weibull 分布参数估计 史景钊，张 峰，陈新昌 （河南农业大学机电工程学院，郑州 450002） 摘 要：介绍了随机截尾情况下计算样本失效概率的两种方法，并编写了Matlab 函数. 提出了利用Matlab 的非线 性最小二乘曲线拟合函数对服从Weibull 分布的随机截尾数据进行曲线拟合和参数估计的方法，并编写了相应程序. 计算结果表明，在Matlab 中只需少量代码即可获得较好的拟合效果和估计精度，比用其他计算机语言编程更简单而 实用. 关键词：可靠性；Weibull 分布；参数估什；Matlab 中图分类号：TB 114.3 文献标识码：A 在产品的寿命试验中有完全寿命试验和截尾寿命试验两种类型. 其中截尾寿命试验又分为定时截尾、 [1] 定数截尾和随机截尾等 . 参加试验的部分产品由于某种原因（如人为因素造成产品损坏、统计数据丢失、 试验设备失效、根据试验计划有意撤出等）还没有失效就中途退出试验，这样得到的数据即为随机截尾数据 （也称为右删失数据）；还有现场可靠性数据，是对实际运行的设备进行寿命考察，由于设备使用情况相差较 大，造成观察的结束时间各不相同，也会出现随机截尾的问题. 例如，考察一个运输公司某种汽车的首次大 修里程，得到的就是随机截尾数据，即有些已经大修过，有些尚未大修，如果只记录已经大修的数据，不记录 尚未大修的数据，就会使评估的结果偏离实际. 随机截尾寿命试验是可靠性寿命试验中最一般的情况，其 他寿命试验都可看作它的一个特例. Weibull 分布模型能够根据形状参数的变化表现为各种不同的形状，较好地适用于各类寿命试验数据， 因而在可靠性分析中应用十分广泛. 对于服从Weibull 分布的随机截尾寿命数据的参数估计，国内外学者 [1- 4] [5-8] [9] 进行了大量的研究，提出了一些参数估计方法，主要有极大似然估计法 、贝叶斯估计法 、最小二乘法 、 [10] 图估计法 等. 本文根据文献［11-12］介绍的计算样本失效概率的方法，编写了Matlab 函数，利用Matlab 强 大的数值计算功能，实现了用非线性最小二乘法对随机截尾寿命数据进行分布拟合和参数估计. 1 样本失效概率及Matlab 实现 假设投入寿命试验的产品数量（即样本容量）为n ，产品的寿命为随机变量T，其分布函数为F （t），相应 的样本失效概率为F （t），在试验结束时，其中有r 个产品发生了失效，其失效时间为x ≤x ≤…≤x ，有k＝n－r 1 2 r 个产品由于各种原因中途撤出了试验，其撤出时间分别为y ≤y ≤…≤y ，则观察到的随机截尾寿命数据其 1 2 k 时间按从小到大排序后