广义指数分布下循环序加试验模型的参数估计.pdfVIP

第41卷 第 1期 曲 阜 师 范 大 学 Vo1．41 N0．1 2015年 1月 Journal of Qufu Normal Jan．2015 广义指数分布下循环序加试验模型的参数估计 邱雷颦 (① 闽江师范高等专科学校计算机系，350002，福建省福州市) 摘要 ：针对产品寿命服从广义指数分布的情况，讨论了循环序进应力加速寿命试验的模型，给出相应参 数的3种估计方法：极大似然估计法、极大似然估计与最小二乘估计相结合的估计方法，以及最小距离估计和 最小二乘估计相结合的估计方法 ，并利用数据模拟 比较 了它们的优劣性． 关键词 ：广义指数分布；加速寿命试验；循环序进应力；最小距离估计 中图分类号：O213．2 文献标识码 ：A 文章编号 ：1001—5337(2015)01—0033—06 在进行序加试验操作过程中，当试验进行到应力达到的限度时，往往 尚未有足够多的产品失效，给统计 分析带来 了困难．文献[1]提出了循环序加试验，让 尚未失效的产品再在相 同的序进应力下继续重复试验 ，直 到有足够多的产品失效时试验停止．从数学上看，如果应力函数是周期函数，则序加试验就是循环序加寿命 试验 ，应力为循环序进应力．文献[1]和文献[2]分别讨论了Weibul1分布，对数正态分布场合下的循环序加 寿命试验的模型和数据分析．然而，在很多情况下，广义指数分布更符合产品寿命数据的分布，它是指数分布 的一种推广 ]．近年来对这一分布的研究已有了一些成果，如文献1-43讨论了其次序统计量的矩，文献 [5]给 出了其贝叶斯统计推断．以下针对产品寿命服从广义指数分布的循环序加寿命试验 ，从几个基本假定 出发 ， 给出了数学模型中相应参数 的几种估计方法 ，并利用数据模拟技术 比较了这几种估计方法 的优劣． 1 数学模型 1．1 基本假定 (i)在恒定应力 S下，产品的寿命分布服从广义指数分布，分布函数为F(￡)一 [1一exp(一 ￡)] ，其 中r， 0，且都是应力水平 5的函数． 与S之间的加速方程为 —exp{一 (n+6 (s))}． (i1)循环序进应力 s(￡)是时间t的周期函数 (周期为 T)，且在一个周期内，s(￡)是非负递增有界的连 续函数． (iii)在一定的应力范围内，产 品的失效机理保持不变_6]． (jv)产品满足累积失效模型，即产品的剩余寿命仅依赖于 已累积失效的部分和当时的应力水平 ，而与 累积的方式无关_7]． 1．2 建立模型 根据假定 (i)，可将时刻 t时应力S(￡)视为恒定应力 ，此时 与时间t之间的关系为 (￡)===exp{一(n+ 6 (S()))}． 首先，对任意给定 的循环序进应力函数 S=S(￡)，可以构造出相应的步进应力函数．分割区间Eo，T]，0 一t。t …t t …t一T ，令 s()：垒=S(￡)，tH tt，i一1，2，…， ．由基本假定(iij)， r 三 z-，i= 1，2，… ，n ． 引理 1嘲 循环序进应力 s(t)下产品的寿命分布函数近似为 收稿 日期 ：2014—09—23 作者简介：邱雷颦，女，1980一，硕士，讲师；研究方向：概率与数理统计．E—mail：qiuleipin@163．corn 34 曲阜师范大学学报(自然科学版) 2O15年 F(z)一l1一exp(--I(“)du)I， 其 中，r 0， — (S(f))是时间的函数． 引理 2 假设产品满足四个基本假定 ：(i)一(iv)，则产品在任何循环应力下所承受的循环次数 的概 率分布为 P 一P (一k)===(1一 q ) 一 (1一 q ) ，走一0，1，2，… ， r下 其中q—exp{--I ()du)． J 0 证明 产品可以承受 k次 (k一0，1，2，… )循环的概率为 P 一