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

第28卷 第5期 河 南 科 学 Vol_28 No．5 2010年 5月 HENAN SCIENCE May2010 文章编号：1004—3918(2010)05—0584—04 基于Matlab的随机截尾数据下的 Weibul1分布参数估计 史景钊， 张 峰， 陈新昌 (河南农业大学机 电工程学院，郑州 4500O2) 摘 要：介绍了随机截尾情况下计算样本失效概率的两种方法，并编写了Matlab函数 ．提出了利用Matlab的非线 性最小二乘曲线拟合函数对服从Weibul1分布的随机截尾数据进行曲线拟合和参数估计的方法，并编写了相应程序． 计算结果表明，在 Matlab中只需少量代码即可获得较好的拟合效果和估计精度，比用其他计算机语言编程更简单而 实用 ． 关键词：可靠性；Weibul1分布；参数估什；Matlab 中图分类号：TB114．3 文献标识码：A ． 在产品的寿命试验中有完全寿命试验和截尾寿命试验两种类型．其中截尾寿命试验又分为定时截尾、 定数截尾和随机截尾等ll1．参加试验的部分产品由于某种原因 (如人为因素造成产品损坏、统计数据丢失、 试验设备失效、根据试验计划有意撤出等)还没有失效就中途退出试验，这样得到的数据即为随机截尾数据 (也称为右删失数据)；还有现场可靠性数据，是对实际运行的设备进行寿命考察，由于设备使用情况相差较 大，造成观察的结束时间各不相同，也会出现随机截尾的问题．例如，考察一个运输公司某种汽车的首次大 修里程，得到的就是随机截尾数据，即有些已经大修过，有些尚未大修，如果只记录已经大修的数据，不记录 尚未大修的数据，就会使评估的结果偏离实际．随机截尾寿命试验是可靠性寿命试验中最一般的情况，其 他寿命试验都可看作它的一个特例 ． Wejbul】分布模型能够根据形状参数的变化表现为各种不同的形状，较好地适用于各类寿命试验数据， 因而在可靠性分析中应用十分广泛 ．对于服从Weibul1分布的随机截尾寿命数据的参数估计，国内外学者 进行了大量的研究，提出了一些参数估计方法，主要有极大似然估计 41、贝叶斯估计法5【．81、最小二乘法恻、 图估计法㈣等 ．本文根据文献 [11．123介绍的计算样本失效概率的方法，编写了Matlab函数，利用Matlab强 大的数值计算功能，实现了用非线性最小二乘法对随机截尾寿命数据进行分布拟合和参数估计 ． 1 样本失效概率及Matlab实现 假设投入寿命试验的产品数量 (即样本容量)为 n，产品的寿命为随机变量 ，其分布函数为F(￡)，相应 的样本失效概率为 ()，在试验结束时，其 中有 r个产品发生了失效，其失效时间为 ≤ ≤…≤轴 有k=n—r 个产品由于各种原因中途撤出了试验，其撤出时间分别为Y。≤)，2≤…≤ ，则观察到的随机截尾寿命数据其 时间按从d,N大排序后可表示为 产品失效时 q=l，2，…，r；i=l，2，…，n)， I产品撤出时 =1，2，…，k；1，2，…，凡)． 显然这类数据不能按照完全样本数据的处理方法计算样本失效概率，必须寻找其他合适的方法． 1．1 Johnson的平均次序法 Johnson认为中途撤出试验的产品会造成失效产品的时间次序发生变化，应该计算失效产品的平均次序 号!”1，第r个失效产品的平均次序号为 -．L1+，r， (1) ． 收稿 日期：2010—02—22 基金项目：河南农业大学博士基金项目 作者简介：史景钊 (1963一)，男，河南柘城人，副教授，博士，主要从事农业装备及其可靠性方面的研究． 2010年5月 史景钊等：基于Maflab的随机截尾数据下的Weibul1分布参数估计