• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 16
    • 下载费用:30 金币  

    重庆时时彩怎么算大小: 基于相似产品信息的小子样产品贝叶斯可靠性评估方法.pdf

    关 键 词:
    基于 相似 产品信息 小子 产品 贝叶斯 可靠性 评估 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610982338.0

    申请日:

    2016.11.08

    公开号:

    CN106570281A

    公开日:

    2017.04.19

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20161108|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 上海无线电设备研究所
    发明人: 李鸿志; 周郁; 黄勇; 张衡; 李雁斌; 张志俊
    地址: 200090 上海市杨浦区黎平路203号
    优先权:
    专利代理机构: 上海信好专利代理事务所(普通合伙) 31249 代理人: 潘朱慧;朱成之
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610982338.0

    授权公告号:

    |||

    法律状态公告日:

    2017.05.17|||2017.04.19

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明涉及一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,包含:S1、对相似产品进行相似度打分,确定先验样本集;S2、获取先验信息,包括相似产品的失效寿命数据与伪失效寿命数据;S3、确定先验分布形式;S4、相似产品失效机理一致性检验;S5、确定先验分布和后验分布,实现可靠度的贝叶斯估计。本发明有效增加了小子样产品的样本量,先验信息包括失效寿命数据和伪失效寿命数据,并将双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,有效提高了产品贝叶斯可靠性评估的精度,简化了贝叶斯可靠性参数估计的计算方法。

    权利要求书

    1.一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,需要进行可靠性评估
    的产品为本产品,与本产品相似的其他产品为相似产品,其特征在于,包含以下步骤:
    S1、对不同的相似产品进行相似度打分,把具有相同相似度的相似产品归类为一个先
    验样本,确定先验样本集;
    S2、获取先验信息,包括相似产品的失效寿命数据与利用性能退化失效模型得到的伪
    失效寿命数据;
    S3、采用双参数威布尔分布为总体分布函数,构造范·蒙特福特统计量,验证先验样本
    的样品寿命服从双参数威布尔分布,确定先验分布形式;
    S4、利用威布尔分布的形状参数表征产品失效机理,构造巴特利特检验统计量,验证相
    似产品失效机理的一致性;
    S5、利用失效机理一致性原则,结合相似度转化成的先验样本权重,把双参数威布尔分
    布贝叶斯估计问题转化成单参数贝叶斯估计,利用贝叶斯公式获取后验分布,得到本产品
    可靠度的贝叶斯估计。
    2.如权利要求1所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特
    征在于,所述的S1中,具体包含以下步骤:
    S11、将相似产品的相似程度取为nj个离散值,相似度记为sj,j=1,2,...nj,且

    S12、记表示第i个专家zi对相似产品pk的相似度打分,i=1,2,...,ni,k=1,2,...,
    nk,对相似产品pk的相似度计算加权和,得到
    S13、为建立样本集,若则认为相似产品pk的相似度为sj,将具有
    相同相似度sj的产品归类为一个先验样本Xj,从而形成先验样本集。
    3.如权利要求2所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特
    征在于,所述的S2中,具体包含以下步骤:
    S21、记先验样本Xj的样品为若样品在使用期间未发生失效,利用其
    性能退化数据及退化失效模型由
    最小二乘估计法,计算得到:
    <mrow> <msup> <mi>SSD</mi> <mi>j</mi> </msup> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>l</mi> <mi>j</mi> </msub> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>&epsiv;</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>l</mi> <mi>j</mi> </msub> </munderover> <msup> <mrow> <mo>&lsqb;</mo> <msubsup> <mi>y</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>g</mi> <mi>p</mi> </mrow> <mi>j</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <msubsup> <mi>f</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> <mi>j</mi> </msubsup> <mo>|</mo> <msubsup> <mi>&beta;</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>&beta;</mi> <mrow> <mi>g</mi> <mn>2</mn> </mrow> <mi>j</mi> </msubsup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>&beta;</mi> <mrow> <mi>g</mi> <mi>p</mi> </mrow> <mi>j</mi> </msubsup> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>;</mo> </mrow>
    估算出之后,确定退化失效模型;
    S22、记性能失效阈值为利用方程求解t,记为即得到先验样本Xj
    中的样品的伪失效寿命数据;
    S23、若先验样本Xj中的样品在使用期间已发生失效,失效寿命为记则得
    到先验样本Xj中所有样品的失效寿命数据分别为
    S24、对先验样本Xj的所有样品的失效寿命数据进行大小排序,得到顺序统计

    4.如权利要求3所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特
    征在于,所述的S3中,具体包含以下步骤:
    S31、建立假设:H0表示先验样本X的样品寿命服从威布尔分布W(m,η),m表征失效机理的
    形状参数,η表征特征寿命的尺度参数;令:
    <mrow> <msub> <mi>Y</mi> <mi>g</mi> </msub> <mo>=</mo> <mi>l</mi> <mi> </mi> <mi>n</mi> <msub> <mover> <mi>t</mi> <mo>~</mo> </mover> <mi>g</mi> </msub> <mo>;</mo> </mrow>
    Zg=(Yg-μ)/σ;
    μ=l nη;
    σ=1/m;
    在假设H0成立的条件下,是极值分布FY(y)=1-exp(-exp((y-μ)/σ))的
    前ng个顺序统计量,是标准极值分布FZ(z)=1-exp(-exp(z))的前ng个顺序
    统计量;
    S32、构造范·蒙特福特统计量为:
    <mrow> <msub> <mi>G</mi> <mi>g</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Y</mi> <mrow> <mi>g</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>Y</mi> <mi>g</mi> </msub> </mrow> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>Z</mi> <mrow> <mi>g</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>Z</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mi>g</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>n</mi> <mi>g</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>;</mo> </mrow>
    取n′g=[ng/2],则统计量为:
    <mrow> <mi>F</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>g</mi> <mo>=</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mo>&prime;</mo> </msubsup> <mo>+</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>n</mi> <mi>g</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mfrac> <msub> <mi>G</mi> <mi>g</mi> </msub> <mrow> <msub> <mi>n</mi> <mi>g</mi> </msub> <mo>-</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mo>&prime;</mo> </msubsup> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </mrow> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>g</mi> <mo>=</mo> <mn>1</mn> </mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mo>&prime;</mo> </msubsup> </munderover> <mfrac> <msub> <mi>G</mi> <mi>g</mi> </msub> <msubsup> <mi>n</mi> <mi>g</mi> <mo>&prime;</mo> </msubsup> </mfrac> </mrow> </mfrac> <mo>;</mo> </mrow>
    在假设H0成立的条件下,渐近服从分布F(2(ng-n′g-1),2n′g);
    S33、对于给定的显著性水平α,如果有F<Fα/2(2(ng-n′g-1),2n′g),或者有F>F1-α/2(2
    (ng-n′g-1),2n′g),则认为假设H0不成立;否则,接收假设H0成立,即先验样本X的样品寿命服
    从威布尔分布W(m,η)。
    5.如权利要求4所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特
    征在于,所述的S4中,失效机理是指引起故障的物理的、化学的、生物的或其他的过程,是产
    品失效的内在原因;利用相似产品信息评估本产品的可靠性时,要求两个产品的失效机理
    一致,即先验样本Xj所服从的威布尔分布的尺度参数mj是相同的。
    6.如权利要求5所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特
    征在于,所述的S4中,具体包含以下步骤:
    S41、建立假设:H1表示由于σj=1/mj,此假设等价于极值分布G(σj,μj)
    中各σj相等,即等价于假设:H′1表示
    S42、记是σj的线性无偏估计,的方差为:
    <mrow> <mi>var</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>&sigma;</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> </msub> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> <mo>;</mo> </mrow>
    S43、按照近似服从自由度有为的χ2分布,构造巴特利特检验统计量:
    <mrow> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>=</mo> <mn>2</mn> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>j</mi> </msub> </munderover> <msubsup> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>&lsqb;</mo> <mi>ln</mi> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>j</mi> </msub> </munderover> <msubsup> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mover> <mi>&sigma;</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>ln</mi> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>j</mi> </msub> </munderover> <msubsup> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>-</mo> <mn>2</mn> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>j</mi> </msub> </munderover> <msubsup> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mi>ln</mi> <msub> <mover> <mi>&sigma;</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
    <mrow> <mi>C</mi> <mo>=</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <mn>6</mn> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mi>j</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>&lsqb;</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> </msub> <mo>-</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msubsup> <mi>l</mi> <mrow> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> <mo>,</mo> <msubsup> <mi>n</mi> <mi>g</mi> <mi>j</mi> </msubsup> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&rsqb;</mo> <mo>;</mo> </mrow>
    S44、在假设H′1成立的条件下,B2/C近似服从自由度为nj-1的χ2分布;对于给定的显著性
    水平α,当时,认为假设H′0不成立;否则,认为假设H′0成立,即假设H1成立。
    7.如权利要求6所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特
    征在于,所述的S5中,具体包含以下步骤:
    S51、记先验样本Xj的失效分布函数为:
    <mrow> <msup> <mi>F</mi> <mi>j</mi> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>/</mo> <msub> <mi>&eta;</mi> <mi>j</mi> </msub> </mrow> <mo>)</mo> </mrow> <msub> <mi>m</mi> <mi>j</mi> </msub> </msup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
    进行数据变换得到:
    <mrow> <mi>l</mi> <mi> </mi> <mi>n</mi> <mi> </mi> <mi>l</mi> <mi> </mi> <mi>n</mi> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>F</mi> <mi>j</mi> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>=</mo> <msub> <mi>m</mi> <mi>j</mi> </msub> <mi>x</mi> <mo>-</mo> <msub> <mi>m</mi> <mi>j</mi> </msub> <msub> <mi>&mu;</mi> <mi>j</mi> </msub> <mo>;</mo> </mrow>
    μj=l nηj;
    x=l n t;
    令Bj=mjμj;则得到Aj=mjx-Bj;
    利用失效分布函数和最小二乘估计法,根据先验样本Xj的失效寿命数据,获得参
    数估计值
    S52、记为本产品的现场试验信息,其为本产品的现场试验的失效寿命数
    据或伪失效寿命数据的顺序统计量,按S4所述的方法对其进行威布尔分布假设的验证,并
    按S51所述的方法获取其形状参数估计值
    S53、经过S4的威布尔分布假设的验证,相似产品与本产品失效机理一致,即
    则把相似产品的相似度sj转化成形状参数mj的权重,融合先验信息
    和现场试验信息,得到后验分布形状参数m的估计值
    <mrow> <mover> <mi>m</mi> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>j</mi> </msub> </munderover> <msub> <mi>s</mi> <mi>j</mi> </msub> <msub> <mover> <mi>m</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>+</mo> <msub> <mover> <mi>m</mi> <mo>^</mo> </mover> <mn>0</mn> </msub> </mrow> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>n</mi> <mi>j</mi> </msub> </munderover> <msub> <mi>s</mi> <mi>j</mi> </msub> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow>
    其中,为先验样本Xj服从的威布尔分布形状参数mj的估计值,为现场试验信息服
    从的威布尔分布形状参数m0的估计值;从而将双参数贝叶斯可靠性估计转换成单参数贝叶
    斯可靠性估计;
    S54、此时后验分布形状参数参数m已知,则是以为参数
    的指数分布的前g0个顺序统计量,把参数η的贝叶斯可靠性估计转化成指数分布进行处理;
    S55、根据贝叶斯公式,后验分布为:
    <mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>&lambda;</mi> <mo>|</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <munder> <mo>&Integral;</mo> <mi>A</mi> </munder> <mi>p</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&lambda;</mi> </mrow> </mfrac> <mo>;</mo> </mrow>
    采用共轭先验分布,取先验分布为Γ(a,b),则有:
    π(λ)∝λa-1exp(-bλ);
    其中,符号“∝”表示两边仅相差一个不依赖于λ的常数因子;此时后验分布为:
    <mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>&lambda;</mi> <mo>|</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>1</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mover> <mi>m</mi> <mo>^</mo> </mover> </msup> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>2</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mover> <mi>m</mi> <mo>^</mo> </mover> </msup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <msub> <mi>g</mi> <mn>0</mn> </msub> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mover> <mi>m</mi> <mo>^</mo> </mover> </msup> <mo>)</mo> </mrow> <mo>&Proportional;</mo> <msup> <mi>&lambda;</mi> <mrow> <mi>a</mi> <mo>+</mo> <msub> <mi>g</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mo>(</mo> <mrow> <mi>T</mi> <mo>+</mo> <mi>b</mi> </mrow> <mo>)</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
    <mrow> <mi>T</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>g</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>g</mi> <mn>0</mn> </msub> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mi>g</mi> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mover> <mi>m</mi> <mo>^</mo> </mover> </msup> <mo>;</mo> </mrow>
    则参数λ的贝叶斯估计为:
    <mrow> <mover> <mi>&lambda;</mi> <mo>^</mo> </mover> <mo>=</mo> <mi>E</mi> <mrow> <mo>(</mo> <mi>&lambda;</mi> <mo>|</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>1</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>2</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <msub> <mi>g</mi> <mn>0</mn> </msub> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mo>+</mo> <msub> <mi>g</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>b</mi> <mo>+</mo> <mi>T</mi> </mrow> </mfrac> <mo>;</mo> </mrow>
    <mrow> <mover> <mi>R</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>/</mo> <mover> <mi>&eta;</mi> <mo>^</mo> </mover> </mrow> <mo>)</mo> </mrow> <mover> <mi>m</mi> <mo>^</mo> </mover> </msup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
    由于可靠度R(t)=exp(-λt),l n R(t)=-λt,其是λ的线性函数,从而求得可靠度R
    (t)的后验分布为:
    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>&lt;</mo> <mi>q</mi> <mo>|</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>1</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>2</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <msub> <mi>g</mi> <mn>0</mn> </msub> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>)</mo> </mrow> <mo>=</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&lambda;</mi> <mo>&gt;</mo> <mo>-</mo> <mfrac> <mrow> <mi>ln</mi> <mi> </mi> <mi>q</mi> </mrow> <mi>t</mi> </mfrac> <mo>|</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>1</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>2</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <msub> <mi>g</mi> <mn>0</mn> </msub> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>ln</mi> <mi> </mi> <mi>q</mi> </mrow> <mi>t</mi> </mfrac> </mrow> <mi>&infin;</mi> </msubsup> <mfrac> <msup> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>T</mi> <mo>)</mo> </mrow> <mrow> <mi>a</mi> <mo>+</mo> <msub> <mi>g</mi> <mn>0</mn> </msub> </mrow> </msup> <mrow> <mi>&Gamma;</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <msub> <mi>g</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mi>&lambda;</mi> <mrow> <mi>a</mi> <mo>+</mo> <msub> <mi>g</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mo>(</mo> <mrow> <mi>b</mi> <mo>+</mo> <mi>T</mi> </mrow> <mo>)</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&lambda;</mi> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
    进而得到可靠度R(t)的贝叶斯估计,即后验期望为:
    <mrow> <mover> <mi>R</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>E</mi> <mrow> <mo>(</mo> <mi>R</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>|</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>1</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <mn>2</mn> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mover> <mi>t</mi> <mo>~</mo> </mover> <msub> <mi>g</mi> <mn>0</mn> </msub> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mi>m</mi> </msup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
    τ=(T+b)/t。

    说明书

    基于相似产品信息的小子样产品贝叶斯可靠性评估方法

    技术领域

    本发明涉及一种小子样产品贝叶斯可靠性评估方法,具体是指一种基于相似产品
    信息的小子样产品贝叶斯可靠性评估方法,属于可靠性建模、可靠性评估的技术领域。

    背景技术

    卫星在航天领域发挥着非常重要的作用,其可靠性和寿命直接影响卫星运行姿态
    控制、测量精度等指标,是影响航天任务成败的关键。卫星在轨长期服役,这就要求卫星的
    有效载荷必须可靠性高、寿命长。因此,对于卫星有效载荷,在全寿命周期的各个阶段,必须
    开展对应的可靠性工作。然而,由于卫星有效载荷本身具有的造价高昂、系统复杂、小子样
    的特点,不能实现将大量产品用于寿命试验,从而缺乏相关的寿命数据用于可靠性分析和
    验证工作。其次,在轨卫星有效载荷处于失重、真空、强辐射、高温差等复杂、恶劣的环境中,
    地面模拟试验环境很难做到对真实运行环境的全覆盖,传统的可靠性评估技术难以准确反
    映产品失效的本质。因此,利用有限的失效寿命数据统计出来的失效率模型指导建立产品
    可靠性模型,对于产品的可靠性评估会存在较大的误差,无法准确判断产品的可靠性水平。

    对于小子样失效寿命数据缺乏统计推断的问题,常规方法具有很大的局限性,目
    前有许多学者提出了新的小子样产品可靠性分析方法。其中,《基于最大熵试验法的小样本
    产品可靠性验证方法》,通过加大样本载荷或减少样本承载能力,使单个样本的试验熵尽可
    能大,达到用小样本验证产品高可靠性、高置信度水平的目的?!痘谙嗨撇沸畔⒌哪诚?br />统任务可靠性预计方法》,利用环境试验、可靠性研制试验等不同试验的相似产品试验信
    息,扩大了统计样本量,提高了小子样系统可靠性评估的可靠性。该方法可用于同一型号不
    同研制阶段的产品,也适用于具有一定继承性的新旧不同型号的产品。

    对于小子样产品的可靠性统计分析问题,贝叶斯方法已得到越来越多的应用。贝
    叶斯统计推断方法把概率分布中的参数作为未知的服从一定概率分布的随机变量,因此可
    以充分利用大量的先验信息,结合小样本信息进行产品可靠性评估。贝叶斯方法在提高产
    品可靠性评估准确性的同时,可以减少试验样本、缩短试验时间从而节约试验经费、缩短研
    制周期。由于贝叶斯方法的优越性,许多学者在这一领域开展了大量新的研究工作。其中,
    有学者针对高可靠长寿命产品失效寿命数据越来越难以获得的现状,基于加速退化试验实
    时获取的大量的产品性能退化数据,结合贝叶斯方法提出《基于多源退化数据的贝叶斯可
    靠性综合评估方法》。也有学者提出了一种基于混合贝叶斯先验分布的可靠性验证测试方
    法,通过先验矩方法和最大熵方法分别求出共轭先验分布的两组不同参数,按第二类极大
    似然方法确定的两种先验分布的权重融合两组参数,进而得到最终的先验分布,比单纯使
    用其中一种方法得到的先验分布更加准确?;褂醒д呓⒘讼低扯匆端雇绾拖低呈?br />效模式的动态贝叶斯网络模型,并离散化转化成多个单连通或多连通静态贝叶斯网络,运
    用静态贝叶斯网络和桶消元法进行双向推导分别求解的方法求解航空发动机涡轮叶盘系
    统的失效率和各失效模式的故障率。

    但目前贝叶斯可靠性评估方法的研究多集中在利用同型号或同类产品寿命数据
    或性能退化数据作为先验信息评估产品可靠性方面,未将相似产品寿命数据或性能退化数
    据纳入先验信息范畴。对于单批次小批量生产的产品,可以获取的同型号产品样本较少,缺
    乏稳定的寿命数据,利用常规的统计分析方法,可靠度模型参数容易存在较大的误差,不能
    很好的反映产品真实的可靠性。许多单批次小批量产品,如航天、航空类产品,由于自身的
    特殊性,在设计上往往具有很强的继承性,很多新产品是在老产品的基础上继承发展的,因
    此可以利用与本产品相似的已有产品的寿命数据来评估本产品的可靠性。通过对相似产品
    可靠性数据的挖掘获取更多的有效先验信息来构造先验分布,结合现场试验数据,建立产
    品可靠性参数的后验分布,实现融合相似产品可靠性信息的小子样产品贝叶斯可靠性评
    估。

    发明内容

    本发明的目的在于提供一种基于相似产品信息的小子样产品贝叶斯可靠性评估
    方法,有效增加了小子样产品的样本量,先验信息包括失效寿命数据和伪失效寿命数据,并
    将双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,有效提高了产品贝叶斯
    可靠性评估的精度,简化了贝叶斯可靠性参数估计的计算方法。

    为了达到上述目的,本发明提供一种基于相似产品信息的小子样产品贝叶斯可靠
    性评估方法,其中,需要进行可靠性评估的产品为本产品,与本产品相似的其他产品为相似
    产品,包含以下步骤:

    S1、对不同的相似产品进行相似度打分,把具有相同相似度的相似产品归类为一
    个先验样本,确定先验样本集;

    S2、获取先验信息,包括相似产品的失效寿命数据与利用性能退化失效模型得到
    的伪失效寿命数据;

    S3、采用双参数威布尔分布为总体分布函数,构造范·蒙特福特统计量,验证先验
    样本的样品寿命服从双参数威布尔分布,确定先验分布形式;

    S4、利用威布尔分布的形状参数表征产品失效机理,构造巴特利特检验统计量,验
    证相似产品失效机理的一致性;

    S5、利用失效机理一致性原则,结合相似度转化成的先验样本权重,把双参数威布
    尔分布贝叶斯估计问题转化成单参数贝叶斯估计,利用贝叶斯公式获取后验分布,得到本
    产品可靠度的贝叶斯估计。

    所述的S1中,具体包含以下步骤:

    S11、将相似产品的相似程度取为nj个离散值,相似度记为sj,j=1,2,...nj,且

    S12、记表示第i个专家zi对相似产品pk的相似度打分,i=1,2,...,ni,k=1,
    2,...,nk,对相似产品pk的相似度计算加权和,得到

    S13、为建立样本集,若则认为相似产品pk的相似度为sj,
    将具有相同相似度sj的产品归类为一个先验样本Xj,从而形成先验样本集。

    所述的S2中,具体包含以下步骤:

    S21、记先验样本Xj的样品为若样品在使用期间未发生失效,利
    用其性能退化数据及退化失效模型
    由最小二乘估计法,计算得到:


    估算出之后,确定退化失效模型;

    S22、记性能失效阈值为利用方程求解t,记为即得到先验样
    本Xj中的样品的伪失效寿命数据;

    S23、若先验样本Xj中的样品在使用期间已发生失效,失效寿命为记
    则得到先验样本Xj中所有样品的失效寿命数据分别为

    S24、对先验样本Xj的所有样品的失效寿命数据进行大小排序,得到顺序
    统计量

    所述的S3中,具体包含以下步骤:

    S31、建立假设:H0表示先验样本X的样品寿命服从威布尔分布W(m,η),m表征失效
    机理的形状参数,η表征特征寿命的尺度参数;令:


    Zg=(Yg-μ)/σ;

    μ=lnη;

    σ=1/m;

    在假设H0成立的条件下,是极值分布FY(y)=1-exp(-exp((y-μ)/
    σ))的前ng个顺序统计量,是标准极值分布FZ(z)=1-exp(-exp(z))的前ng
    个顺序统计量;

    S32、构造范·蒙特福特统计量为:


    取n′g=[ng/2],则统计量为:


    在假设H0成立的条件下,渐近服从分布F(2(ng-n′g-1),2n′g);

    S33、对于给定的显著性水平α,如果有F<Fα/2(2(ng-n′g-1),2n′g),或者有F>F1-α/2(2
    (ng-n′g-1),2n′g),则认为假设H0不成立;否则,接收假设H0成立,即先验样本X的样品寿命服
    从威布尔分布W(m,η)。

    所述的S4中,失效机理是指引起故障的物理的、化学的、生物的或其他的过程,是
    产品失效的内在原因;利用相似产品信息评估本产品的可靠性时,要求两个产品的失效机
    理一致,即先验样本Xj所服从的威布尔分布的尺度参数mj是相同的。

    所述的S4中,具体包含以下步骤:

    S41、建立假设:H1表示m1=m2=...=mnj;由于σj=1/mj,此假设等价于极值分布G
    (σj,μj)中各σj相等,即等价于假设:H1′表示

    S42、记是σj的线性无偏估计,的方差为:


    S43、按照近似服从自由度有为的χ2分布,构造巴特利特检验统计
    量:



    S44、在假设H′1成立的条件下,B2/C近似服从自由度为nj-1的χ2分布;对于给定的
    显著性水平α,当时,认为假设H′0不成立;否则,认为假设H′0成立,即假设
    H1成立。

    所述的S5中,具体包含以下步骤:

    S51、记先验样本Xj的失效分布函数为:


    进行数据变换得到:


    μj=lnηj;

    x=lnt;

    令Bj=mjμj;则得到Aj=mjx-Bj;

    利用失效分布函数和最小二乘估计法,根据先验样本Xj的失效寿命数据,获
    得参数估计值

    S52、记为本产品的现场试验信息,其为本产品的现场试验的失效寿
    命数据或伪失效寿命数据的顺序统计量,按S4所述的方法对其进行威布尔分布假设的验
    证,并按S51所述的方法获取其形状参数估计值

    S53、经过S4的威布尔分布假设的验证,相似产品与本产品失效机理一致,即
    则把相似产品的相似度sj转化成形状参数mj的权重,融合先验信息
    和现场试验信息,得到后验分布形状参数m的估计值


    其中,为先验样本Xj服从的威布尔分布形状参数mj的估计值,为现场试验信
    息服从的威布尔分布形状参数m0的估计值;从而将双参数贝叶斯可靠性估计转换成单参数
    贝叶斯可靠性估计;

    S54、此时后验分布形状参数参数m已知,则是以为
    参数的指数分布的前g0个顺序统计量,把参数η的贝叶斯可靠性估计转化成指数分布进行
    处理;

    S55、根据贝叶斯公式,后验分布为:


    采用共轭先验分布,取先验分布为Γ(a,b),则有:

    π(λ)∝λa-1exp(-bλ);

    其中,符号“∝”表示两边仅相差一个不依赖于λ的常数因子;此时后验分布为:



    则参数λ的贝叶斯估计为:



    由于可靠度R(t)=exp(-λt),lnR(t)=-λt,其是λ的线性函数,从而求得可靠度R
    (t)的后验分布为:


    进而得到可靠度R(t)的贝叶斯估计,即后验期望为:


    τ=(T+b)/t。

    综上所述,本发明提供的基于相似产品信息的小子样产品贝叶斯可靠性评估方
    法,有效弥补了小子样产品在进行可靠性评估时样本量不足的问题,提高了对产品可靠性
    评估的精度;先验信息不仅包括失效寿命数据,还包括伪失效寿命数据,提高了贝叶斯可靠
    性评估的精度;将双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,简化了
    贝叶斯可靠性参数估计的计算方法。

    附图说明

    图1为本发明中的基于相似产品信息的小子样产品贝叶斯可靠性评估方法的流程
    图。

    具体实施方式

    以下结合图1,详细说明本发明的优选实施例。

    如图1所示,为本发明所述的基于相似产品信息的小子样产品贝叶斯可靠性评估
    方法,其中,需要进行可靠性评估的产品称为本产品,与本产品相似的其他产品称为相似产
    品,包含以下步骤:

    S1、对不同的相似产品进行相似度打分,把具有相同相似度的相似产品归类为一
    个先验样本,确定先验样本集;

    S2、获取先验信息,包括相似产品的失效寿命数据与利用性能退化失效模型得到
    的伪失效寿命数据;

    S3、采用双参数威布尔分布为总体分布函数,构造范·蒙特福特统计量,验证先验
    样本的样品寿命服从双参数威布尔分布,确定先验分布形式;

    S4、利用威布尔分布的形状参数表征产品失效机理,构造巴特利特检验统计量,验
    证相似产品失效机理的一致性;

    S5、利用失效机理一致性原则,结合相似度转化成的先验样本权重,把双参数威布
    尔分布贝叶斯估计问题转化成单参数贝叶斯估计,利用贝叶斯公式获取后验分布,得到本
    产品可靠度的贝叶斯估计。

    所述的S1中,在工程实际中,相关技术领域的专家对产品的设计、制造、功能和使
    用环境具有丰富的经验,对相似产品进行相似度打分具体是指:在充分考虑相似要素数量
    和相似要素相似度的基础上,对相似产品与本产品的相似程度给出一个相对可信性较高的
    评估;具体包含以下步骤:

    S11、将相似产品的相似程度取为nj个离散值,相似度记为sj,j=1,2,...nj,且

    S12、记表示第i个专家zi对相似产品pk的相似度打分,i=1,2,...,ni,k=1,
    2,...,nk,对相似产品pk的相似度计算加权和,得到

    S13、为建立样本集,若则认为相似产品pk的相似度为sj,
    将具有相同相似度sj的产品归类为一个先验样本Xj,从而形成先验样本集。

    所述的S2中,若相似产品在使用期间已发生失效,则可以直接获取该相似产品的
    失效寿命数据;若相似产品在使用期间未发生失效,将关键性能退化信息或与之相关的参
    数作为时间的函数,获取退化轨迹,使用退化失效模型或回归方程进行描述,当性能退化信
    息达到预先设定的失效阈值所需要的时间,即可获取该相似产品的伪失效寿命数据。

    所述的S2中,具体包含以下步骤:

    S21、记先验样本Xj的样品为若样品在使用期间未发生失效,
    利用其性能退化数据及退化失效模型
    由最小二乘估计法,计算得到:


    估算出之后,确定退化失效模型;

    S22、记性能失效阈值为利用方程求解t,记为即得到先验样
    本Xj中的样品的伪失效寿命数据;

    S23、若先验样本Xj中的样品在使用期间已发生失效,失效寿命为记
    则得到先验样本Xj中所有样品的失效寿命数据分别为

    S24、对先验样本Xj的所有样品的失效寿命数据进行大小排序,得到顺序
    统计量

    所述的S3中,采用双参数威布尔分布作为总体分布函数,并进行各样本分布函数
    的假设检验。该双参数威布尔分布是可靠性定量分析中常用的一类典型寿命分布,在可靠
    性工程中得到广泛的应用。

    所述的S3中,具体包含以下步骤:

    S31、不失一般性,建立假设:H0表示先验样本X的样品寿命服从威布尔分布W(m,
    η),m表征失效机理的形状参数,η表征特征寿命的尺度参数;令:


    Zg=(Yg-μ)/σ;

    μ=lnη;

    σ=1/m;

    在假设H0成立的条件下,是极值分布FY(y)=1-exp(-exp((y-μ)/
    σ))的前ng个顺序统计量,是标准极值分布FZ(z)=1-exp(-exp(z))的前ng
    个顺序统计量;

    S32、构造范·蒙特福特统计量为:


    取n′g=[ng/2],则统计量为:


    在假设H0成立的条件下,渐近服从分布F(2(ng-n′g-1),2n′g);

    S33、对于给定的显著性水平α,如果有F<Fα/2(2(ng-n′g-1),2n′g),或者有F>F1-α/2(2
    (ng-n′g-1),2n′g),则认为假设H0不成立;否则,接收假设H0成立,即先验样本X的样品寿命服
    从威布尔分布W(m,η)。

    所述的S4中,失效机理是指引起故障的物理的、化学的、生物的或其他的过程,是
    产品失效的内在原因。利用相似产品信息评估本产品的可靠性时,要求两个产品的失效机
    理一致,即先验样本Xj所服从的威布尔分布的尺度参数mj是相同的。

    所述的S4中,具体包含以下步骤:

    S41、建立假设:H1表示由于σj=1/mj,此假设等价于极值分布G
    (σj,μj)中各σj相等,即等价于假设:H1′表示

    S42、记是σj的线性无偏估计,的方差为:


    S43、按照近似服从自由度有为的χ2分布,构造巴特利特检验统计
    量:



    S44、在假设H1′成立的条件下,B2/C近似服从自由度为nj-1的χ2分布;对于给定的
    显著性水平α,当时,认为假设H′0不成立;否则,认为假设H′0成立,即假设
    H1成立。

    所述的S5中,具体包含以下步骤:

    S51、记先验样本Xj的失效分布函数为:


    进行数据变换得到:


    μj=lnηj;

    x=lnt;

    令Bj=mjμj;则得到Aj=mjx-Bj;

    利用失效分布函数和最小二乘估计法,根据先验样本Xj的失效寿命数据,获
    得参数估计值

    S52、记为本产品的现场试验信息,其为本产品的现场试验的失效寿
    命数据或伪失效寿命数据的顺序统计量,按S4所述的方法对其进行威布尔分布假设的验
    证,并按S51所述的方法获取其形状参数估计值

    S53、经过S4的威布尔分布假设的验证,相似产品与本产品失效机理一致,即
    则把相似产品的相似度sj转化成形状参数mj的权重,融合先验信息
    和现场试验信息,得到后验分布形状参数m的估计值


    其中,为先验样本Xj服从的威布尔分布形状参数mj的估计值,为现场试验信
    息服从的威布尔分布形状参数m0的估计值;从而将双参数贝叶斯可靠性估计转换成单参数
    贝叶斯可靠性估计;

    S54、此时后验分布形状参数参数m已知,则是以为
    参数的指数分布的前g0个顺序统计量,把参数η的贝叶斯可靠性估计转化成指数分布进行
    处理;

    S55、根据贝叶斯公式,后验分布为:


    采用共轭先验分布,取先验分布为Γ(a,b),则有:

    π(λ)∝λa-1exp(-bλ);

    其中,符号“∝”表示两边仅相差一个不依赖于λ的常数因子;此时后验分布为:



    则参数λ的贝叶斯估计为:



    由于可靠度R(t)=exp(-λt),lnR(t)=-λt,其是λ的线性函数,从而求得可靠度R
    (t)的后验分布为:


    进而得到可靠度R(t)的贝叶斯估计,即后验期望为:


    τ=(T+b)/t。

    综上所述,本发明所提供的基于相似产品信息的小子样产品贝叶斯可靠性评估方
    法,通过确定已有的相似产品与需要进行可靠性评估的本产品的相似度;利用性能退化模
    型推导相似产品的伪失效寿命数据,拓展了贝叶斯先验信息的获取,结合已有的失效寿命
    数据,建立先验分布形式;根据失效机理一致性原则,结合由相似度转化的先验样本权重,
    把双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计;最后,根据贝叶斯公式
    融合先验分布和试验信息获取后验分布,获取可靠度的贝叶斯估计。

    因此,本发明所提供的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,
    具有以下优点和有益效果:

    1、充分融合了专家经验、相似产品信息、现场试验信息,弥补了小子样产品在进行
    可靠性评估时样本量不足的问题,提高了对产品可靠性评估的精度。

    2、先验信息不仅包括已有的失效寿命数据,还利用产品性能退化数据建立退化模
    型推导伪失效寿命数据,从而获取了更多的先验信息,提高了贝叶斯可靠性评估的精度。

    3、基于失效机理一致性的原则,结合由相似度转化的先验样本权重,把双参数贝
    叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,简化了贝叶斯可靠性参数估计的计
    算方法。

    尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的
    描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的
    多种修改和替代都将是显而易见的。因此,本发明的?;し段вτ伤降娜ɡ罄聪薅?。

    关于本文
    本文标题:基于相似产品信息的小子样产品贝叶斯可靠性评估方法.pdf
    链接地址://www.4mum.com.cn/p-6092756.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 60彩票群 手机捕鱼之海底捞攻略 pk10牛牛app 快慢高手一波中特 封闭基金 赚钱吗 双色球赚钱 意甲冠军次数排名 英雄联盟隐藏分查询 北京赛车8码 天天捕鱼秘籍 急速赛车漂移 甘肃十一选五走势图解 普通免费单机麻将 摇奖现场视频' 034期特码资料 捕鱼大师pc版