[发明专利]基于贝叶斯理论和线性拟合的可靠度置信下限构建方法在审
申请号: | 201811463446.2 | 申请日: | 2018-12-03 |
公开(公告)号: | CN109614586A | 公开(公告)日: | 2019-04-12 |
发明(设计)人: | 贾祥;程志君;郭波 | 申请(专利权)人: | 中国人民解放军国防科技大学 |
主分类号: | G06F17/18 | 分类号: | G06F17/18 |
代理公司: | 长沙国科天河知识产权代理有限公司 43225 | 代理人: | 邱轶 |
地址: | 410073 湖*** | 国省代码: | 湖南;43 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 置信 可靠度 构建 线性拟合 贝叶斯 样本 产品寿命 试验数据 威布尔分布函数 可靠性统计 数学表达式 工程实践 过程计算 理论计算 失效概率 下限曲线 线性函数 转化 应用 | ||
1.基于贝叶斯理论和线性拟合的威布尔分布可靠度置信下限构建方法,其特征在于,包括如下步骤:
(1)对于给定的各个样本产品寿命试验数据即各样本产品的失效时间,首先根据贝叶斯理论计算各个样本产品失效时间处失效概率的置信上限;
(2)将威布尔分布函数转化为线性函数;
(3)基于线性拟合构建威布尔分布可靠度的置信下限曲线,从而构建了可靠度的置信下限。
2.根据权利要求1所述的贝叶斯理论和线性拟合的威布尔分布可靠度置信下限构建方法,其特征在于,步骤(1)的实现方法如下:
设给定的n个样本产品的失效时间分别为t1,…,tn,并假定t1≤…≤tn,称pi=F(ti)为失效时间ti的失效概率,其中
根据贝叶斯理论,引入失效概率pi的验前分布,取该验前分布的分布形式为贝塔分布Beta(pi;a,b),
其中a和b为贝塔分布的2个分布参数,在此取贝塔分布的2个分布参数同时为0.5,即设pi的验前分布为贝塔分布Beta(pi;0.5,0.5);
再确定失效概率pi的似然函数为
根据贝叶斯公式
推得失效概率pi的验后分布为
即贝塔分布Beta(pi;1.5,n-i+0.5);
在给定的置信水平(1-α)下,记为失效概率pi的置信上限,为的估计结果,则即为贝塔分布Beta(pi;1.5,n-i+0.5)的(1-α)分位点,将代入式(7)可得于是在置信水平(1-α)下,可得失效概率pi的置信上限估计结果可通过调用数学软件Matlab中的函数betainv来获得。
3.根据权利要求2所述的贝叶斯理论和线性拟合的威布尔分布可靠度置信下限构建方法,其特征在于,步骤(2)的实现方法如下:
威布尔分布函数为:
将威布尔分布函数通过两次取对数变换,可得
ln[-ln(1-F)]=mlnt-lnη (8)
根据失效概率的定义,可令p=F(t)代表失效概率;再记y=ln[-ln(1-p)],x=lnt,从而可将威布尔分布函数转化为线性函数y=mx-mlnη,其中(x,y)为坐标轴中的点。
从而可将威布尔分布函数转化为线性函数。
4.根据权利要求3所述的贝叶斯理论和线性拟合的威布尔分布可靠度置信下限构建方法,其特征在于,步骤(3)的实现方法如下:
基于n个样本产品的失效时间t1≤…≤tn以及各个失效时间处失效概率pi的置信上限的估计值令xi=lnti,如此可得坐标轴中的n个点
设是真值yi的估计值,与yi之间存在一定的误差;根据线性回归的思想,利用诸点拟合式(8)中的线性函数,当拟合误差和最小时,可得威布尔分布参数m和η的估计为:
再将威布尔分布参数m和η的估计值和代入威布尔可靠度函数中,可得拟合后的可靠度置信下限曲线为
5.根据权利要求2、3或4所述的贝叶斯理论和线性拟合的威布尔分布可靠度置信下限构建方法,其特征在于,步骤(1)中,样本产品的数量n≥2。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国人民解放军国防科技大学,未经中国人民解放军国防科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811463446.2/1.html,转载请声明来源钻瓜专利网。