一种将测试性虚拟试验数据转换为实物试验数据的方法
【专利摘要】本发明公开了一种将测试性虚拟试验数据转换为实物试验数据的方法。目的是提供一种将虚拟试验数据转换为实物试验数据的方法，将虚拟试验样本量转换为等效的实物试验样本量；技术方案是：利用测试性试验的HPD可信区间平均长度和试验的可信度计算试验的设计效应，通过方程求解得到虚拟试验样本量对应的等效的实物试验样本量，从而将测试性虚拟试验样本量转换为等效的实物试验样本量；采用本发明获得的等效实物试验样本量贴近工程实际；该方法可以应用于各类产品的测试性试验，对加快试验进度、降低试验成本有意义。
【专利说明】一种将测试性虚拟试验数据转换为实物试验数据的方法
【技术领域】
[0001] 本发明涉及测试性试验的方法，尤其是一种将测试性虚拟试验数据转换为实物试 验数据的方法。
【背景技术】
[0002] 测试性是指"产品能及时准确地确定其状态（可工作、不可工作或性能下降）并有 效地隔离其内部故障的一种设计特性"。测试性试验是指通过故障模拟和故障注入等手段， 检查产品的机内测试系统（设备）和外部测试系统（设备）对故障检测和隔离能力，从而 暴露产品测试性设计方面的缺陷，确认是否符合测试性定量要求所进行的试验。
[0003] 测试性试验分为虚拟试验、实物试验。实物试验是指在模拟实际使用条件或在规 定的工作及环境条件下的试验，其试验对象和测试系统都是真实的，故障样本则是人工制 造出来的，其中有的与自然故障在物理、功能表现等各方面几乎完全一致，有的仅在功能表 现方面等效。虚拟试验指在人工构建的虚拟对象上进行测试性试验，理论上它可以虚拟测 试性试验的全部要素。
[0004] 目前，在测试性核查和验证等测试性工作中所依托的测试性试验都是实物试验， 其特点是必须向产品中注入故障，这就不可避免地存在一些故障不能注入或无法有效注 入，比如危害性大、不可撤销的故障，比如在封装产品内部的故障，等等。按照统计学原理， 在一次试验中，如果达到一定数量的试验用例无法获取，将会造成一次试验无效。
[0005] 对此，人们在积极寻找解决办法，比如采取故障等效注入的方法。这些方法不能从 根本上解决问题。依照产品发展趋势来看，产品集成和封装越来越来紧密，光机电一体化产 品越来越多，故障机理更加复杂，故障注入的风险越来越大，且后果难以预料。
[0006] 如果采取虚拟试验，那么上述问题将得到根本性地解决，而且还有更多的好处。首 先，人们可以利用测试性虚拟样机开展故障机理研究，从而帮助人们找到避免或延缓故障 产生的方法；其次，虚拟试验不受场地、时间和次数的限制，并且可以对试验过程进行回放， 从而大大提髙试验效率，降低试验成本。
[0007] 对于虚拟试验数据的最简单处理方法是直接将其等效为实物试验数据，但是在实 际应用中，考虑到测试性虚拟试验模型可信度等问题，虚拟试验并不能直接等效为实物试 验，虚拟试验数据也不能直接等效为实物试验数据。虚拟试验数据这种与实物试验数据相 关性极强，但又不完全等效的特性，给虚实结合的测试性一体化试验方案设计和试验评估 等带来新问题和新挑战。解决此类问题的一个思路是，将虚拟试验数据转换为等效实物试 验数据，这样可以将新问题转化为解决方案比较成熟的传统测试性问题，因此，探索将虚拟 试验数据转换为等效实物试验数据的方法，是新兴的虚实结合的测试性一体化试验技术中 一个重要研究方向。
[0008] 在测试性试验数据中，测试性试验样本量的大小关系到测试性试验方案设计、指 标评估各个方面。虚实结合一体化试验常遇到2个问题，一个问题是已经开展了一定数量 的虚拟试验如何确定后续实物试验样本量；另一个是测试性实物试验开展不充分，指标评 估精度不局，如何利用虚拟试验数据提尚评估精度。这2个问题的核心都与虚拟试验样本 量如何合理利用有关。如何将虚拟试验样本量转化为等效的实物试验样本量是测试性虚拟 试验数据转换中的关键问题。因此，本发明要解决的关键技术问题是，如何将虚拟试验样本 量转换为等效的实物试验样本量。这一问题尚未见到相关解决方案。
【发明内容】
[0009] 本发明要解决的技术问题是：提供一种方法，将虚拟试验数据转换为实物试验数 据，即将虚拟试验样本量转换为等效的实物试验样本量。
[0010] 本发明的技术方案如下：
[0011] 第一步，根据建模资料和专家评判信息，采用基于专家评判的灰色综合方法（具 体内容见参考文献（杨丽娜等，《基于专家评判的灰色综合方法在可信度评估中的应用研 究》，航天控制，2007, 25 (4) :63-67))计算虚拟试验模型的可信度Cs。
[0012] 第二步，从被测产品的测试性试验大纲或者历史试验报告中获得测试性虚拟试 验样本量ns ;取故障检测率p的无信息先验分布为贝塔分布：π。= Beta( α。，β。)，贝塔分 布的两个参数aQ和β。的值可根据经验给定（一般的，可取a。= β。= 1)，或者按照 Jeffreys方法计算获得（Jeffreys方法参考文献（韦来生等编著，《贝叶斯分析》，中国科学 技术大学出版社，2013,第53-55页））；确定故障检测率p落在可信区间之外的概率a，取 值a彡〇. 1。
[0013] 第三步，故障成功检测次数rs取值0, 1，…，ns，按公式（1)计算虚拟试验样本量为 ns、故障成功检测次数为rs时，故障检测率p的后验分布，公式（1)中的r赋值为rs，公式 (1)中的η赋值为ns，得到后验分布函数族{ n (p | (ns，rs)) | rs = 〇,；[，…，。
[0014] π (p I (n，r)) = Beta (r+a 0，（n_r) + β 0) (1)
[0015] 上式中，n表示试验样本量；r表示试验中故障成功检测次数，其取值范围为 0, 1，…，n, Beta(r+a。，（η-ι·) + β。）表示参数为r+a。和（η-Γ) + β。的贝塔分布（贝塔分布 的定义见参考文献（李萍等编著，《应用概率统计》，科学出版社，2013,第159-161页））。
[0016] 第四步，采用最大后验密度（Highest Posterior Density，HPD)可信区间逼近算 法（算法步骤见参考文献（王国玉等编著，《电子系统小子样试验理论方法》，国防工业出版 社，2003,第95-96页））计算后验分布函数族W (p| (ns，rs)) |rs = 0，l，…，ns}中各后验分 布函数的（1-a )HPD可信区间随机长度，得到一组随机长度{La (rs|ns) |rs = 0, 1,…，ns}。
[0017] 第五步，把rj武值给r，ns赋值给n，按照公式（2)所示故障成功检测次数匕的边 缘密度函数表达式计算试验样本量为1时、不同故障成功检测次数r s对应的边缘密度函数 值{m(rs | ns) | rs = 0, 1，... ns}。
[0018] Q . m (2)
[0019] 上式中，:C|表示n选r的组合数，B(r+a。，（n-r) + β。）表示参数为r+ a。和 (η-：τ) + β。的贝塔函数（贝塔函数的定义见参考文献（李萍等编著，《应用概率统计》，科学 出版社，2013,第160页）），Β( α。，β。）表示参数为a。和β。的贝塔函数。
[0020] 第六步，把rs赋值给r，ns赋值给n，La (rs|ns)赋值给La (r|n)，m(rs|ns)赋值 给m(r|n)，按照公式⑶计算样本量为ns的测试性虚拟试验（l-a)HPD区间平均长度 E[La(rs|ns)]。
[0021]
【权利要求】
1. 一种将测试性虚拟试验数据转换为实物试验数据的方法，其特征在于包括以下步 骤： 第一步，根据建模资料和专家评判信息，采用基于专家评判的灰色综合方法计算虚拟 试验模型的可信度Cs; 第二步，从被测产品的测试性试验大纲或者历史试验报告中获得测试性虚拟试验样本 量ns ;取故障检测率P的无信息先验分布为贝塔分布^ =Beta(ameJ，贝塔分布的两个 参数a〇和的值根据经验给定或者按照Jeffreys方法计算获得；确定故障检测率p落 在可信区间之外的概率a; 第三步，故障成功检测次数rs取值0, 1，…，ns，按公式（1)计算虚拟试验样本量为ns、 故障成功检测次数为&时，故障检测率p的后验分布，公式（1)中的r赋值为rs，公式（1) 中的n赋值为ns，得到后验分布函数族{>(pI(ns,rs)) |rs = 0, 1，…，ns}; (pI(n，r)) =Beta(r+a0,（n-r) + 3 0) (I) 上式中,n表示试验样本量；r表示试验中故障成功检测次数，其取值范围为0, 1，一,n,Beta(r+a。，（n-r) + @。）表示参数为r+a。和（n-r) + @。的贝塔分布； 第四步，采用最大后验密度即HH)可信区间逼近算法计算后验分布函数族 {> (p| (ns，rs)) |rs = 0, 1，…，ns}中各后验分布函数的（l-a)HPD可信区间随机长度， 得到一组随机长度ILa (rs|ns) |rs = 0, 1，...，ns}，HPD的含义是最大后验密度即Highest PosteriorDensity的缩写； 第五步，把rs赋值给r，ns赋值给n，按照公式（2)所示故障成功检测次数rs的边缘密 度函数表达式计算试验样本量为ns时、不同故障成功检测次数rs对应的边缘密度函数值 {m(rs|ns) |rs = 0,I,---nj：
上式中，G表示n选r的组合数，B(r+aQ, (n-r) +P。）表示参数为r+a。和（n-r) + 3Q 的贝塔函数，BUtl,Ptl)表示参数为a^和的贝塔函数； 第六步，把rs赋值给r，ns赋值给n，La (rs|ns)赋值给La(r|n)，m(rs|ns)赋值给m(r|n)， 按照公式⑶计算样本量为ns的测试性虚拟试验（l-a)HPD区间平均长度E[La(rs|ns)]:
上式中，La(r|n)表示试验样本量为n，故障成功检测次数为r所对应的后验分布函数 的（l-ci)HPD可信区间长度，此时，故障检测率p落在可信区间内的概率为（1-a) ;m(r|n) 表示试验样本量为n，故障成功检测次数为r所对应的边缘密度函数值； 第七步，把rs赋值给r，ns赋值给n，E[La (rs |ns)]赋值给E[La(r|n)]，按照公式⑷计 算样本量为ns的测试性虚拟试验的设计效应Ed (nsIcs): Ed (nIc) = c?exp (-E[La(r In) ]) (4) 上式中，c表示测试性虚拟试验可信度，exp( ?)表示自然指数，对于实物试验，c= 1 ; 对于虚拟试验，c=Cs，在（0, 1)的区间内取值； 第八步，建立方程（5)，解方程（5)，得到虚拟试验样本量ns的等效实物试验样本量I， 解算步骤如下： ED(ns|cs) =ED(ne|c=I) (5) 8.Dn6初始值设为I; 8. 2)1^取值0, 1，…，Iiej，把;Tej赋值给r，Iiej赋值给n，按公式（1)计算等效实物试验 样本量为I、故障成功检测次数为&时，故障检测率p的后验分布，得到后验分布函数族 {>(pI(ne，re))Ire = 0, 1，…，ne}; 8. 3)根据HH)可信区间逼近算法计算后验分布函数族中各后验分布函数的HH)可信区 间随机长度，得到一组随机长度ILa (re |r〇 |rs = 0, 1，…，rU; 8.4)把rJ武值给r，nJ武值给n，按照公式（2)所示边缘密度函数表达式计算等效实 物试验样本量为时、不同故障成功检测次数&对应的边缘密度函数值{m(r」r〇 |re = 0, 1，…ne}; 8. 5)把 & 赋值给r，:? 赋值给n，LaCrejIr〇 赋值给La(rIn)，mCrejIr〇 赋值给m(rIn), 按照公式⑶计算等效实物试验的（l-a)HPD可信区间平均长度E[La(r」r〇]; 8. 6)c取值为1，把re赋值给r，ne赋值给n,E[La (reIr〇 ]赋值给E[La(rIn)]，按照公式 (4)计算等效试验的设计效应EdOieJc= 1); 8. 7)如果Ed (neIc=I) <Ed (nsIcs)，则ne 更新为ne+l，转步骤 8. 2);若Ed (neIc= 1)彡Ed (nsIc=Cci)，则结束，此时得到的即为虚拟试验样本量ns的等效实物试验样本量。
2. 如权利要求1所述的将测试性虚拟试验数据转换为实物试验数据的方法，其特征在 于所述贝塔分布的两个参数a^和的值取= 1。
3. 如权利要求1所述的将测试性虚拟试验数据转换为实物试验数据的方法，其特征在 于所述故障检测率P落在可信区间之外的概率a<0.1。
【文档编号】G06F19/00GK104268386SQ201410477812
【公开日】2015年1月7日 申请日期:2014年9月18日 优先权日:2014年9月18日
【发明者】杨鹏, 张勇, 刘瑛, 邱静, 刘冠军, 吕克洪, 王超, 沈亲沐, 王刚, 赵晨旭, 赵志傲, 吴超, 谢皓宇, 李乾, 林辰龙, 李华康, 季明江, 彭磊 申请人:中国人民解放军国防科学技术大学