一种将测试性虚拟试验数据转换为实物试验数据的方法技术

本发明专利技术公开了一种将测试性虚拟试验数据转换为实物试验数据的方法。目的是提供一种将虚拟试验数据转换为实物试验数据的方法，将虚拟试验样本量转换为等效的实物试验样本量；技术方案是：利用测试性试验的HPD可信区间平均长度和试验的可信度计算试验的设计效应，通过方程求解得到虚拟试验样本量对应的等效的实物试验样本量，从而将测试性虚拟试验样本量转换为等效的实物试验样本量；采用本发明专利技术获得的等效实物试验样本量贴近工程实际；该方法可以应用于各类产品的测试性试验，对加快试验进度、降低试验成本有意义。

【技术实现步骤摘要】
一种将测试性虚拟试验数据转换为实物试验数据的方法
本专利技术涉及测试性试验的方法，尤其是一种将测试性虚拟试验数据转换为实物试验数据的方法。
技术介绍
测试性是指“产品能及时准确地确定其状态(可工作、不可工作或性能下降)并有效地隔离其内部故障的一种设计特性”。测试性试验是指通过故障模拟和故障注入等手段，检查产品的机内测试系统(设备)和外部测试系统(设备)对故障检测和隔离能力，从而暴露产品测试性设计方面的缺陷，确认是否符合测试性定量要求所进行的试验。测试性试验分为虚拟试验、实物试验。实物试验是指在模拟实际使用条件或在规定的工作及环境条件下的试验，其试验对象和测试系统都是真实的，故障样本则是人工制造出来的，其中有的与自然故障在物理、功能表现等各方面几乎完全一致，有的仅在功能表现方面等效。虚拟试验指在人工构建的虚拟对象上进行测试性试验，理论上它可以虚拟测试性试验的全部要素。目前，在测试性核查和验证等测试性工作中所依托的测试性试验都是实物试验，其特点是必须向产品中注入故障，这就不可避免地存在一些故障不能注入或无法有效注入，比如危害性大、不可撤销的故障，比如在封装产品内部的故障，等等。按照统计学原理，在一次试验中，如果达到一定数量的试验用例无法获取，将会造成一次试验无效。对此，人们在积极寻找解决办法，比如采取故障等效注入的方法。这些方法不能从根本上解决问题。依照产品发展趋势来看，产品集成和封装越来越来紧密，光机电一体化产品越来越多，故障机理更加复杂，故障注入的风险越来越大，且后果难以预料。如果采取虚拟试验，那么上述问题将得到根本性地解决，而且还有更多的好处。首先，人们可以利用测试性虚拟样机开展故障机理研究，从而帮助人们找到避免或延缓故障产生的方法；其次，虚拟试验不受场地、时间和次数的限制，并且可以对试验过程进行回放，从而大大提高试验效率，降低试验成本。对于虚拟试验数据的最简单处理方法是直接将其等效为实物试验数据，但是在实际应用中，考虑到测试性虚拟试验模型可信度等问题，虚拟试验并不能直接等效为实物试验，虚拟试验数据也不能直接等效为实物试验数据。虚拟试验数据这种与实物试验数据相关性极强，但又不完全等效的特性，给虚实结合的测试性一体化试验方案设计和试验评估等带来新问题和新挑战。解决此类问题的一个思路是，将虚拟试验数据转换为等效实物试验数据，这样可以将新问题转化为解决方案比较成熟的传统测试性问题，因此，探索将虚拟试验数据转换为等效实物试验数据的方法，是新兴的虚实结合的测试性一体化试验技术中一个重要研究方向。在测试性试验数据中，测试性试验样本量的大小关系到测试性试验方案设计、指标评估各个方面。虚实结合一体化试验常遇到2个问题，一个问题是已经开展了一定数量的虚拟试验如何确定后续实物试验样本量；另一个是测试性实物试验开展不充分，指标评估精度不高，如何利用虚拟试验数据提高评估精度。这2个问题的核心都与虚拟试验样本量如何合理利用有关。如何将虚拟试验样本量转化为等效的实物试验样本量是测试性虚拟试验数据转换中的关键问题。因此，本专利技术要解决的关键技术问题是，如何将虚拟试验样本量转换为等效的实物试验样本量。这一问题尚未见到相关解决方案。
技术实现思路
本专利技术要解决的技术问题是：提供一种方法，将虚拟试验数据转换为实物试验数据，即将虚拟试验样本量转换为等效的实物试验样本量。本专利技术的技术方案如下：第一步，根据建模资料和专家评判信息，采用基于专家评判的灰色综合方法(具体内容见参考文献(杨丽娜等，《基于专家评判的灰色综合方法在可信度评估中的应用研究》，航天控制，2007，25(4)：63-67))计算虚拟试验模型的可信度cs。第二步，从被测产品的测试性试验大纲或者历史试验报告中获得测试性虚拟试验样本量ns；取故障检测率p的无信息先验分布为贝塔分布π0＝Beta(α0,β0)，贝塔分布的两个参数α0和β0的值可根据经验给定(一般的，可取α0＝β0＝1)，或者按照Jeffreys方法计算获得(Jeffreys方法参考文献(韦来生等编著，《贝叶斯分析》，中国科学技术大学出版社，2013，第53-55页))；确定故障检测率p落在可信区间之外的概率α，取值α≤0.1。第三步，故障成功检测次数rs取值0,1,…,ns，按公式(1)计算虚拟试验样本量为ns、故障成功检测次数为rs时，故障检测率p的后验分布，公式(1)中的r赋值为rs，公式(1)中的n赋值为ns，得到后验分布函数族{π(p|(ns,rs))|rs＝0,1,…,ns}。π(p|(n,r))＝Beta(r+α0,(n-r)+β0)(1)上式中，n表示试验样本量；r表示试验中故障成功检测次数，其取值范围为0,1,…,n，Beta(r+α0,(n-r)+β0)表示参数为r+α0和(n-r)+β0的贝塔分布(贝塔分布的定义见参考文献(李萍等编著，《应用概率统计》，科学出版社，2013，第159-161页))。第四步，采用最大后验密度(HighestPosteriorDensity，HPD)可信区间逼近算法(算法步骤见参考文献(王国玉等编著，《电子系统小子样试验理论方法》，国防工业出版社，2003，第95-96页))计算后验分布函数族{π(p|(ns,rs))|rs＝0,1,…,ns}中各后验分布函数的(1-α)HPD可信区间随机长度，得到一组随机长度{Lα(rs|ns)|rs＝0,1,…,ns}。第五步，把rs赋值给r，ns赋值给n，按照公式(2)所示故障成功检测次数rs的边缘密度函数表达式计算试验样本量为ns时、不同故障成功检测次数rs对应的边缘密度函数值{m(rs|ns)|rs＝0,1,…ns}。上式中，表示n选r的组合数，B(r+α0,(n-r)+β0)表示参数为r+α0和(n-r)+β0的贝塔函数(贝塔函数的定义见参考文献(李萍等编著，《应用概率统计》，科学出版社，2013，第160页))，B(α0,β0)表示参数为α0和β0的贝塔函数。第六步，把rs赋值给r，ns赋值给n，Lα(rs|ns)赋值给Lα(r|n)，m(rs|ns)赋值给m(r|n)，按照公式(3)计算样本量为ns的测试性虚拟试验(1-α)HPD区间平均长度E[La(rs|ns)]。上式中，La(r|n)表示试验样本量为n，故障成功检测次数为r所对应的后验分布函数的(1-α)HPD可信区间长度，此时，故障检测率p落在可信区间内的概率为(1-α)；m(r|n)表示试验样本量为n，故障成功检测次数为r所对应的边缘密度函数值。第七步，把rs赋值给r，ns赋值给n，E[La(rs|ns)]赋值给E[La(r|n)]，按照公式(4)计算样本量为ns的测试性虚拟试验的设计效应ED(ns|cs)。ED(n|c)＝c·exp(-E[La(r|n)])(4)上式中，c表示测试性虚拟试验可信度，exp(·)表示自然指数。对于实物试验，c＝1；对于虚拟试验，c＝cs，在(0,1)的区间内取值。第八步，建立方程(5)，解方程(5)，得到虚拟试验样本量ns的等效实物试验样本量ne，解算步骤如下：ED(ns|cs)＝ED(ne|c＝1)(5)8.1)ne初始值设为1。8.2)re取值0,1,…,ne，把re赋值给r，ne赋值给n，按公式(1)计算等效实物试本文档来自技高网...

【技术保护点】
一种将测试性虚拟试验数据转换为实物试验数据的方法，其特征在于包括以下步骤：第一步，根据建模资料和专家评判信息，采用基于专家评判的灰色综合方法计算虚拟试验模型的可信度cs；第二步，从被测产品的测试性试验大纲或者历史试验报告中获得测试性虚拟试验样本量ns；取故障检测率p的无信息先验分布为贝塔分布π0＝Beta(α0,β0)，贝塔分布的两个参数α0和β0的值根据经验给定或者按照Jeffreys方法计算获得；确定故障检测率p落在可信区间之外的概率α；第三步，故障成功检测次数rs取值0,1,…,ns，按公式(1)计算虚拟试验样本量为ns、故障成功检测次数为rs时，故障检测率p的后验分布，公式(1)中的r赋值为rs，公式(1)中的n赋值为ns，得到后验分布函数族{π(p|(ns,rs))|rs＝0,1,…,ns}；π(p|(n,r))＝Beta(r+α0,(n‑r)+β0)    (1)上式中，n表示试验样本量；r表示试验中故障成功检测次数，其取值范围为0,1,…,n，Beta(r+α0,(n‑r)+β0)表示参数为r+α0和(n‑r)+β0的贝塔分布；第四步，采用最大后验密度即HPD可信区间逼近算法计算后验分布函数族{π(p|(ns,rs))|rs＝0,1,…,ns}中各后验分布函数的(1‑α)HPD可信区间随机长度，得到一组随机长度{Lα(rs|ns)|rs＝0,1,…,ns}，HPD的含义是最大后验密度即Highest Posterior Density的缩写；第五步，把rs赋值给r，ns赋值给n，按照公式(2)所示故障成功检测次数rs的边缘密度函数表达式计算试验样本量为ns时、不同故障成功检测次数rs对应的边缘密度函数值{m(rs|ns)|rs＝0,1,…ns}：m(r|n)=CnrB(r+α0,(n-r)+β0)B(α0,β0),r=0,1...n---(2)]]>上式中，表示n选r的组合数，B(r+α0,(n‑r)+β0)表示参数为r+α0和(n‑r)+β0的贝塔函数，B(α0,β0)表示参数为α0和β0的贝塔函数；第六步，把rs赋值给r，ns赋值给n，Lα(rs|ns)赋值给Lα(r|n)，m(rs|ns)赋值给m(r|n)，按照公式(3)计算样本量为ns的测试性虚拟试验(1‑α)HPD区间平均长度E[La(rs|ns)]：E[La(r|n)]=Σr=0nLa(r|n)×m(r|n)---(3)]]>上式中，La(r|n)表示试验样本量为n，故障成功检测次数为r所对应的后验分布函数的(1‑α)HPD可信区间长度，此时，故障检测率p落在可信区间内的概率为(1‑α)；m(r|n)表示试验样本量为n，故障成功检测次数为r所对应的边缘密度函数值；第七步，把rs赋值给r，ns赋值给n，E[La(rs|ns)]赋值给E[La(r|n)]，按照公式(4)计算样本量为ns的测试性虚拟试验的设计效应ED(ns|cs)：ED(n|c)＝c·exp(‑E[La(r|n)])    (4)上式中，c表示测试性虚拟试验可信度，exp(·)表示自然指数，对于实物试验，c＝1；对于虚拟试验，c＝cs，在(0,1)的区间内取值；第八步，建立方程(5)，解方程(5)，得到虚拟试验样本量ns的等效实物试验样本量ne，解算步骤如下：ED(ns|cs)＝ED(ne|c＝1)    (5)8.1)ne初始值设为1；8.2)re取值0,1,…,ne，把re赋值给r，ne赋值给n，按公式(1)计算等效实物试验样本量为ne、故障成功检测次数为re时，故障检测率p的后验分布，得到后验分布函数族{π(p|(ne,re))|re＝0,1,…,ne}；8.3)根据HPD可信区间逼近算法计算后验分布函数族中各后验分布函数的HPD可信区间随机长度，得到一组随机长度{Lα(re|ne)|rs＝0,1,…,ne}；8.4)把re赋值给r，ne赋值给n，按照公式(2)所示边缘密度函数表达式计算等效实物试验样本量为ne时、不同故障成功检测次数re对应的边缘密度函数值{m(re|ne)|re＝0,1,…ne}；8.5)把re赋值给r，ne赋值给n，Lα(re|ne)赋值给Lα(r|n)，m(re|ne)赋值给m(r|n)，按照公式(3)计算等效实物试验的(1‑α)HPD可信区间平均长度E[La(re|ne)]；8.6)c取值为1，把re赋值给r，ne赋值给n，E[La(re|ne)]赋值给E[La(r|n)]，按照公式(4)计算等效试验的设计效应ED(ne|c＝1)；8.7)如果ED(ne|c＝1)＜ED(ns|cs)，则ne更新为ne+1，转步骤8.2)；若ED(ne|c＝1)≥ED(ns|c＝c0)，则结束，此时得到的ne...

【技术特征摘要】
1.一种将测试性虚拟试验数据转换为实物试验数据的方法，其特征在于包括以下步骤：第一步，根据建模资料和专家评判信息，采用基于专家评判的灰色综合方法计算虚拟试验模型的可信度cs；第二步，从被测产品的测试性试验大纲或者历史试验报告中获得测试性虚拟试验样本量ns；取故障检测率p的无信息先验分布为贝塔分布π0＝Beta(α0,β0)，贝塔分布的两个参数α0和β0的值根据经验给定或者按照Jeffreys方法计算获得；确定故障检测率p落在可信区间之外的概率α；所述Jeffreys方法参考中国科学技术大学出版社出版的韦来生、张伟平编著2013版《贝叶斯分析》第53‐55页的Jeffreys方法；第三步，故障成功检测次数rs取值0,1,…,ns，按公式(1)计算虚拟试验样本量为ns、故障成功检测次数为rs时，故障检测率p的后验分布，公式(1)中的r赋值为rs，公式(1)中的n赋值为ns，得到后验分布函数族{π(p|(ns,rs))|rs＝0,1,…,ns}；π(p|(n,r))＝Beta(r+α0,(n-r)+β0)(1)上式中，n表示试验样本量；r表示试验中故障成功检测次数，其取值范围为0,1,…,n，Beta(r+α0,(n-r)+β0)表示参数为r+α0和(n-r)+β0的贝塔分布；第四步，采用最大后验密度即HPD可信区间逼近算法计算后验分布函数族{π(p|(ns,rs))|rs＝0,1,…,ns}中各后验分布函数的(1-α)HPD可信区间随机长度，得到一组随机长度{Lα(rs|ns)|rs＝0,1,…,ns}，HPD的含义是最大后验密度即HighestPosteriorDensity的缩写；第五步，把rs赋值给r，ns赋值给n，按照公式(2)所示故障成功检测次数rs的边缘密度函数表达式计算试验样本量为ns时、不同故障成功检测次数rs对应的边缘密度函数值{m(rs|ns)|rs＝0,1,…ns}：上式中，表示n选r的组合数，B(r+α0,(n-r)+β0)表示参数为r+α0和(n-r)+β0的贝塔函数，B(α0,β0)表示参数为α0和β0的贝塔函数；第六步，把rs赋值给r，ns赋值给n，Lα(rs|ns)赋值给Lα(r|n)，m(rs|ns)赋值给m(r|n)，按照公式(3)计算样本量为ns的测试性虚拟试验(1-α)HPD区间平均长度E[La(rs|ns)]：

【专利技术属性】
技术研发人员：杨鹏，张勇，刘瑛，邱静，刘冠军，吕克洪，王超，沈亲沐，王刚，赵晨旭，赵志傲，吴超，谢皓宇，李乾，林辰龙，李华康，季明江，彭磊，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南;43