2. 中国人民解放军32181部队, 西安 710000
2. Unit 32181 of PLA, Xi'an 710000, China
测试性验证是根据装备的技术规范和指标要求,在承制方和使用方风险约束下,综合考虑试验成本、研制周期等因素,采用单次抽样、多次抽样或序贯检验等方法制定验证试验方案,进而判断是否达到规定的测试性水平[1]。
由于现有中国军标规定的测试性试验方法所需样本量较大,相关学者对基于Bayes融合的试验方法开展了研究,文献[2]运用Bayes理论建立了融合先验信息的决策模型,通过定义不确定性和支持度得到加权因子,然后对先验信息进行融合,明显降低了样本量;文献[3]针对测试性试验数据小子样变总体的特点,建立了指标动态增长模型,并结合最大熵原理得到先验分布,依据Bayes最大后验风险准则确定试验方案,提高了评估结果置信度。文献[4]研究了先验分布参数计算方法,提出混合Beta分布的小子样测试性评估方法。然而,上述方法中的系统级先验信息很难直接获取,先验分布权值确定大都依据数据关系或主观经验,没有考虑先验信息来源,导致现有方法的准确性和适用性存在一定问题。
基于此,本文提出了基于分系统先验信息的测试性试验优化方法,通过引入分系统贡献率对先验信息进行折算,得到系统级先验数据;在此基础上提出了基于理想点排序-层次分析法(Technique for order preference by similarity to ideal solution-analytic hierarchy process, TOPSIS-AHP)的先验分布可信度计算方法,确定了先验分布权值;最后,研究了基于截尾序贯后验加权检验(Sequential posterior odd test, SPOT)的试验方法,并给出了截尾阈值计算方法。
1 基于贡献率的数据折算方法 1.1 测试性多源先验信息分析在装备研制过程中,不同的寿命周期、不同的层次结构、不同的试验手段都会积累多种试验信息,这些信息都可以作为测试性验证试验的先验信息[5]。测试性先验信息的种类较多,从层次结构的角度分为部组件级、分系统级和系统级试验信息;从寿命周期各阶段的角度分为设计阶段信息、生产阶段信息和试验阶段信息;从信息来源的角度分为专家经验信息、虚拟试验信息、摸底试验信息等。测试性多源先验信息分类如图 1所示。
本文从信息来源的角度选取用于Bayes融合的测试性先验信息,主要包括以下5种典型的先验信息:
1) 专家经验信息。由设计人员或领域专家根据自身经验知识给出的测试性指标评估值,通常是连续区间估计或单侧置信下限估计,此类信息主观性较强。
2) 虚拟试验信息。通过搭建装备分系统虚拟样机,利用计算机仿真技术进行测试性验证故障注入试验,代替部分实装故障注入试验,试验结果通常为成败型数据,此类信息的可信度需要考虑。
3) 增长试验信息。部组件或分系统在不用阶段开展的有计划、有继承的测试性试验,各阶段的试验数据服从不同总体,属于“异总体”情况,通常为成败型试验数据。
4) 单元摸底试验信息。在设计研制阶段进行的针对某些部组件或分系统的摸底试验,可信度较高,可以向上综合为系统信息,多为成败型试验数据。
5) 相似装备历史信息。许多装备的设计研制具有继承性,相似装备测试性水平是十分重要的历史信息,可以作为新型号装备的先验信息,通常表示为指标区间估计或单侧置信下限估计。
不同先验信息具有不同的数据形式,为有效利用多源先验信息,需要将不同的数据形式折合为成败型试验数据,相关文献已经给出了详细的成败型试验数据折合方法[6],这里不再赘述。通常采用Beta分布作为成败型试验数据的先验分布[7],分布形式为
$ \pi(p)=\frac{P^{a-1}(1-P)^{b-1}}{B(a, b)} $ | (1) |
式中:p为检测/隔离率,a,b为分布超参数。
1.2 分系统贡献率确定不同分系统在装备对测试性指标评估的影响是不同的,借鉴装备体系贡献率评估中的计算方法[8-9],依据不同分系统的测试性指标分配值计算分系统贡献率,通过指标分配的不同表示在测试性验证评估中,分系统对指标评估的作用大小。
定义1 贡献率ci,在测试性指标评估中,分系统指标作为贡献者,系统指标作为受益者,待测系统的分系统指标的增加或减少对系统指标变化的贡献程度为分系统贡献率,表达式为
$ 分系统A的概率 = \frac{{{{系统指标}_{包括分系统A}} - {{系统指标}_{不包括分系统A}}}}{{{{系统指标}_{包括分系统A}}}} \times 100\% $ | (2) |
第i个分系统的测试性指标分配值为
$ \gamma_i=1-\frac{\lambda_{\mathrm{s}}\left(1-\gamma_{\mathrm{r}}\right)}{k \lambda_i} $ | (3) |
式中:γi为测试性指标分配值,λs为系统的总故障率,由故障模式及危害性分析(Failure mode, effects and criticality analysis, FMECA)确定[10],γr为系统测试性指标要求值,λi为第i个分系统的故障率,k为分系统数量。
由分系统指标计算系统指标时,主要考虑诊断能力问题,即发生和诊断的故障次数或相对故障率,则由测试性指标分配值得到的系统指标值为
$ {\gamma _{\rm{s}}} = \frac{{\sum\limits_{i = 1}^k {{\lambda _i}} {\gamma _i}}}{{\sum\limits_{i = 1}^k {{\lambda _i}} }} $ | (4) |
则第i个分系统的贡献率为
$ {c_i} = \frac{{{\lambda _i}{\gamma _i}}}{{\sum\limits_{i = 1}^k {{\lambda _i}} {\gamma _i}}} $ | (5) |
在可靠性综合评估中,由信息量等效原则和信息量可加性可知,分系统可靠性试验提供的信息量与系统的总信息量相等[11]。在测试性验证中,信息量等效的方法也应用于虚拟试验数据向实装试验数据的折合[12]。基于此,结合分系统贡献率,提出基于分系统贡献率和信息论的数据折算方法。
设系统由k(k≥2)个分系统组成,第i个分系统每次试验的成功概率为pi,则基于分系统贡献率修正后的第i个分系统每次试验的平均信息量[13]为
$ I_i=-c_i\left[\bar{p}_i \ln \bar{p}_i+\left(1-\bar{p}_i\right) \ln \left(1-\bar{p}_i\right)\right] $ | (6) |
设折合后的第i个分系统的成败型试验数据为(ni , fi ),其中ni 为试验总次数,fi为试验失败次数。根据信息量可加性,k个分系统全部试验的总信息量为
$ {I_{{\rm{ss}}}} = - \sum\limits_{i = 1}^k {{c_i}} {n_i}\left[ {{{\bar p}_i}\ln {{\bar p}_i} + \left( {1 - {{\bar p}_i}} \right)\ln \left( {1 - {{\bar p}_i}} \right)} \right] $ | (7) |
设系统等效折算后的试验信息为(n, f),系统每次试验检测/隔离成功的概率为
$ I_{\mathrm{st}}=-n[\tilde{p} \ln \tilde{p}+(1-\tilde{p}) \ln (1-\tilde{p})] $ | (8) |
取pi为对应分系统的极大似然估计,即
$ \bar{p}_i=\frac{n_i-f_i}{n_i} $ | (9) |
结合各分系统贡献率和测试性指标的极大似然估计值,可得各先验信息对应的系统测试性指标估计值为
$ \tilde p = \sum\limits_{i = 1}^k {{c_i}} {{\bar p}_i} $ | (10) |
根据信息量等效原则,即Iss=Ist,则可得k个分系统的成败型试验数据向系统级成败型数据的折算公式为
$ \left\{ {\begin{array}{*{20}{l}} {n = \frac{{\sum\limits_{i = 1}^k {{c_i}} \left[ {{f_i}\ln {f_i} + \left( {{n_i} - {f_i}} \right)\ln \left( {{n_i} - {f_i}} \right) - {n_i}\ln {n_i}} \right]}}{{\tilde p\ln \tilde p + (1 - \tilde p)\ln (1 - \tilde p)}}}\\ {f = n(1 - \tilde p)} \end{array}} \right. $ | (11) |
为保证有效利用先验信息,需要检验先验信息和实装试验数据的相容性。已知先验分布π(p)的类型为Beta分布,采用贝叶斯置信区间估计的方法进行相容性检验。设π(p) 的贝叶斯置信区间为[p1, p2](置信度为1-αh),其中,p1,p2由式(12)确定。
$ \left\{\begin{array}{l} \int_0^{p_1} \pi(p) \mathrm{d} p=\frac{\alpha_{\mathrm{h}}}{2} \\ \int_{p_2}^1 \pi(p) \mathrm{d} p=\frac{\alpha_{\mathrm{h}}}{2} \end{array}\right. $ | (12) |
由实装试验数据得到的指标点估计值为
即使先验信息满足相容性检验条件,其获取方式和数据处理过程的可信程度仍不明确,在用于测试性试验前需要评估其可信度。通常,实装试验数据的可信度最高,以实装数据为基准,用先验分布和实装试验数据分布函数的相似性度量表示先验信息的可信度。设X为实装试验样本,Y为先验信息试验样本,给出如下假设:
$ \left\{\begin{array}{l} H_0: X \text { 与 } Y \text { 是相同的总体 } \\ H_1: X \text { 与 } Y \text { 不是相同的总体 } \end{array}\right. $ | (13) |
记A为接受H0,拒绝H1的事件;A为拒绝H0,拒绝H1的事件。设犯第1类错误的概率为αh,犯第1类错误的概率为βh,则
$ \left\{\begin{array}{l} P\left(\bar{A} \mid H_0\right)=\alpha_{\mathrm{h}} \\ P\left(A \mid H_1\right)=\beta_{\mathrm{h}} \end{array}\right. $ | (14) |
定义2 先验分布可信度ρ,指接受H0的条件下H0为真的概率,表示X与Y来自相同总体的概率,即X与Y的相似程度,其表达式为
$ \begin{aligned} & \rho=\frac{P\left(A \mid H_0\right) P\left(H_0\right)}{P\left(A \mid H_0\right) P\left(H_0\right)+\left(1-P\left(H_0\right)\right) P\left(A \mid H_1\right)}= \\ & \frac{\left(1-\alpha_{\mathrm{h}}\right) P\left(H_0\right)}{\left(1-\alpha_{\mathrm{h}}\right) P\left(H_0\right)+\beta_{\mathrm{h}}\left(1-P\left(H_0\right)\right)} \\ & \end{aligned} $ | (15) |
式(15)中,两类错误概率αh、βh由承制方和使用方协商确定。P(H0)为先验相似度,指在获得先验信息前,先验信息与试验信息属于相同总体的概率,与信息获取途径、数据处理方法、试验方法等多种因素有关,实质上反应的是对先验信息来源和类型的信任程度,属于多准则综合评判问题,很难直接通过显式定量计算,在无法获得先验相似度信息时,通常取P(H0)=0.5。
为降低可信度计算的主观性,提出逼近理想解排序法(TOPSIS)与层次分析法(AHP)相结合的方法确定P(H0)。首先由m名试验参与人员对P(H0)进行初始打分,得到TOPSIS决策矩阵V=
$ \left\{ {\begin{array}{*{20}{l}} {d_i^ + = \sqrt {\sum\limits_{j = 1}^m {{{\left( {{v_{ij}} - v_j^ + } \right)}^2}} } , i = 1, 2, \cdots , n}\\ {d_i^ - = \sqrt {\sum\limits_{j = 1}^m {{{\left( {{v_{ij}} - v_j^ - } \right)}^2}} } , i = 1, 2, \cdots , n} \end{array}} \right. $ | (16) |
则先验信息与理想解的贴近度为
$ \boldsymbol{A}=\left[{a}_{i j}\right]_{n \times n} $ | (17) |
依据先验信息贴近度,采用三标度层次分析法[15]对先验信息信任程度进行两两比较,建立比较矩阵
$ {u_i} = \sum\limits_{j = 1}^n {{a_{ij}}} $ | (18) |
设信任程度判断矩阵为
$ s_{i j}= \begin{cases}\frac{u_i-u_j}{u_{\max }-u_{\min }}\left(b_{\mathrm{m}}-1\right)+1, & u_i>u_j \\ 1, & u_i=u_j \\ {\left[\frac{u_j-u_i}{u_{\max }-u_{\min }}\left(b_{\mathrm{m}}-1\right)\right]^{-1}, } & u_i<u_j\end{cases} $ | (19) |
式中,bm为基点比较标度。
对判断矩阵S进行按行求和,并进行标准化后可得各先验分布的P(H0)值,代入式(15)可得先验分布可信度ρi。
设有r种先验信息通过了相容性检验,其可信度为ρi, i=1, 2, …, r,可得混合先验分布为
$ \left\{ {\begin{array}{*{20}{l}} {\pi (p) = \sum\limits_{i = 1}^r {{\omega _i}} {\pi _i}(p)}\\ {{\omega _i} = {\rho _i}/\sum\limits_{i = 1}^r {{\rho _i}} } \end{array}} \right. $ | (20) |
基于分系统先验信息的测试性试验方案设计流程如图 2所示。
首先,由分系统先验信息和贡献率,通过分系统数据折算得到系统级先验信息;然后,对先验信息进行相容性检验,并通过计算先验分布可信度得到先验分布权值;最后,基于混合先验分布,采用序贯验后加权检验(SPOT)法进行试验方案设计。
3.2 序贯验后加权检验试验方案相比于单次抽样方法和序贯概率比检验(SPRT)方法[16],融合先验信息SPOT法可以有效减少试验样本量[17],故采用SPOT方法进行试验方案设计。设试验方案约束条件为:最低可接受值p0,设计要求值p1,承制方风险α,使用方风险β。
设测试性指标p的检验假设为:H0:p∈Θ0,H1:p∈Θ1。其中,
则对于实装试验结果为(n, c)的样本X=(X1, X2, …, Xn),验后加权比[18]为
$ O_n=\frac{\int_{p_1}^1 p^{n+a-c+1}(1-p)^{n+a-c+1} \mathrm{~d} p}{\int_0^{p_0} p^{n+a-c+1}(1-p)^{n+a-c+1} \mathrm{~d} p} $ | (21) |
检验准则为:当On≤A时,采纳假设H0,判为拒收;当On≥B时,采纳假设H1,判为接收;当A < On < B时,不做判决,继续试验。A,B为常数,由式(22)、(23)得到:
$ \left\{\begin{array}{l} A=\frac{\beta_{\pi_1}}{\int_0^{p_0} \pi(p) \mathrm{d} p-\alpha_{\pi_0}} \\ B=\frac{\int_{p_1}^1 \pi(p) \mathrm{d} p-\beta_{\pi_1}}{\alpha_{\pi_0}} \end{array}\right. $ | (22) |
$ \left\{\begin{array}{l} \alpha_{\pi_0}=\int_{\Theta_0} \alpha \cdot \pi(p) \mathrm{d} p \\ \beta_{\pi_1}=\int_{\Theta_1} \beta \cdot \pi(p) \mathrm{d} p \end{array}\right. $ | (23) |
上述判决准则可能会一直出现A < On < B的情况,造成试验样本量很大。因此,需要制定截尾试验方案,在A,B间插入合理的检验点C。当A < On < C时,判为拒收;当C < On < B时,判为接收。设截尾样本量为ns,成功次数为xs,记
$ h\left(x_{\mathrm{s}}\right)=\frac{1-\int_0^{p_1} \frac{\varGamma\left(a+b+n_{\mathrm{s}}\right)}{\varGamma\left(a+x_{\mathrm{s}}\right) \varGamma\left(b+n_{\mathrm{s}}-x_{\mathrm{s}}\right)} p^{a+x_{\mathrm{s}}-1}(1-p)^{b+n_{\mathrm{s}}-x_{\mathrm{s}}-1} \mathrm{~d} p}{\int_0^{p_0} \frac{\varGamma\left(a+b+n_{\mathrm{s}}\right)}{\varGamma\left(a+x_{\mathrm{s}}\right) \varGamma\left(b+n_{\mathrm{s}}-x_{\mathrm{s}}\right)} p^{a+x_{\mathrm{s}}-1}(1-p)^{b+n_{\mathrm{s}}-x_{\mathrm{s}}-1} \mathrm{~d} p} $ | (24) |
令Ma=h-1(A),Mc=h-1(C),Mb=h-1(B),则由截尾判决造成的承制方和使用方风险增量分别为:
$ \Delta {\alpha _{{{\rm{s}}_0}}} = \frac{{\mathit{\Gamma }(a + b)}}{{\mathit{\Gamma }(a)\mathit{\Gamma }{{(b)}_{k = {M_c} + 1}}}}\sum\limits_{k = {M_a} + 1}^{{M_c}} {C_{{n_s}}^k} {(1 - p)^{b + {n_{\rm{s}}} - k - 1}}{\rm{d}}p $ | (25) |
$ \Delta {\beta _{{\rm{s}}{\pi _1}}} = \frac{{\mathit{\Gamma }(a + b)}}{{\mathit{\Gamma }(a)\mathit{\Gamma }(b)}}\sum\limits_{k = {M_a} + 1}^{{M_c}} {C_{{n_s}}^k} \int_{{p_1}}^1 {{p^{a + k - 1}}} {(1 - p)^{b + {n_{\rm{s}}} - k - 1}}{\rm{d}}p $ | (26) |
由Bayes后验风险准则[19],对于给定的风险增量上界αsπ0和βsπ1,将确定截尾阈值转化为以下优化问题:
$ \begin{array}{*{20}{c}} {\mathop {\min }\limits_{{n_{\rm{s}}}, {x_{\rm{s}}} \in {N^ + }} }&{{n_{\rm{s}}}, {x_{\rm{s}}}}\\ {{\rm{ s}}{\rm{. t}}{\rm{. }}}&{\Delta {\alpha _{{\rm{s}}{\pi _0}}} \le {\alpha _{{\rm{s}}\pi 0}}}\\ {}&{\Delta {\beta _{{\rm{s}}{\pi _1}}} \le {\beta _{{{\rm{s}}^2}{\pi _1}}}}\\ {}&{{M_a} < {M_c} < {M_b}} \end{array} $ | (27) |
将式(27)的优化结果代入式(24)即可得到检验点C。
4 案例分析某航空装备控制系统由4个分系统组成,分别为控制功能分系统、发射功能分系统、通信功能分系统和导航功能分系统。以故障检测率(Fault detection rate, FDR)为例进行测试性验证试验方案设计。根据装备测试性设计要求确定SPOT试验方案的约束条件为:最低可接受值p0=0.85,设计要求值p1=0.95,双方风险α=β=0.1,截尾方案中风险增量上界αsπ0=βsπ1=0.01。
4.1 分系统数据折算由文献[6]的数据等效模型可得各分系统先验信息的成败型试验数据见表 1。
根据试验双方确定的控制系统FMECA,由式(3)~(5) 可得分系统贡献率矩阵为c=[0.282 6, 0.399 1, 0.190 4, 0.127 9]。将表 1数据和分系统贡献率代入式(9)、(10)可得各先验信息对应的指标估计值为
在置信度0.9条件下对先验分布进行相容性检验,根据式(12)可得相容性检验区间分别为[0.748 4, 0.987 7],[0.747 2, 0.919 5],[0.821 5, 0.957 9],而由实装试验数据得到的点估计值
依据基于TOPSIS-AHP的先验分布融合方法可得先验分布的可信度矩阵ρ=[0.930 6, 0.667 2, 0.686 6],进而可得先验分布取值向量为[0.402 8, 0.290 7, 0.301 2]。
4.3 试验方案设计由式(20)可得混合先验分布为
$ \pi (p) = \frac{{\mathit{\Gamma }(34.71)}}{{\mathit{\Gamma }(30.42)\mathit{\Gamma }(4.29)}}{p^{30.42 - 1}}{(1 - p)^{4.29 - 1}} $ | (28) |
图 3为先验分布函数,以及考虑先验分布可信度和未考虑先验分布可信度的分布函数曲线。
将试验约束条件和混合先验分布代入式(22)、(23),可得判决阈值A=0.024 0,B=1.950 1。根据风险增量约束,可得截尾SPOT方案的截尾次数为43,此时检验点C=0.253 5。
对控制系统进行实装故障注入试验,设当前试验结果为(ni, ci),代入式(21)可得当前试验方案的验后加权比值,与判决阈值进行对比,即可得出判决结果,判决过程如图 4所示。
在相同试验约束条件下,采用同为序贯类试验方案的SPRT方法的判决过程如图 5所示。
由图 4可以看出,对于基于可信度融合的判决过程,当试验结果为(N, C)=(36, 2)时,做出接收判决;对于直接融合的判决过程,当试验结果为(N, C)=(38, 2)时,做出接收判决。由图 5可以看出,当试验的成功失败数为(N, C)=(42, 2)时,做出接收判决,此时试验样本量为42。
4.4 试验结果对比分析对控制系统进行多次试验,分别采用SPRT、SPOT方法进行判决,均做出接收判决,不同失败次数下的试验次数对比见表 3。
为进一步验证该方法的有效性,在相同的试验约束条件下,分别采用经典定数抽样试验方案、文献[2]传统Bayes试验方案和文献[20]的改进Bayes试验方案,各方案的对比分析结果见表 4。
根据上述试验结果,从先验信息获取、先验分布权值、试验样本量、双方风险4个方面对比分析本文试验方案的总体性能:
1) 先验信息获取。由表 1、2可知,本文方法可以有效实现分系统试验数据折算,与文献[20]方法相比,考虑了分系统贡献率,使得数据折算更为合理。
2) 先验分布权值。基于可信度确定权值为[0.402 8, 0.290 7, 0.301 2],文献[2]方法确定权值为[0.293 5, 0.313 6, 0.392 8],文献[2]确定的虚拟试验信息和专家经验信息权值较大,而本文方法确定的单元试验信息权值较大,相对不准确的专家经验信息权值较小,更符合试验实际情况。由图 3也可以看出,直接融合先验分布(先验分布权值相等)得到的指标点估计值小于基于可信度的方法,这是因为专家经验对指标的估计比较保守,但其可信度较低,若直接融合先验分布,导致可信度较低的专家经验数据一定程度上湮没可信度较高的单元试验数据,使得专家经验对评估结果产生的影响较大,导致融合结果过于保守。由图 4(b)和表 3也可以看出,直接融合先验分布的试验方案在相同失败次数下所需的样本量更多,导致评估结果过于保守。
3) 试验样本量。由表 3可以看出,同为截尾序贯试验方案,在相同判决结果下,基于可信度融合先验分布和直接融合先验分布下的SPOT试验方案与SPRT试验方案相比,样本量分别平均减少18.6%和13.1%;由表 4可以看出,传统定数抽样方案得到的样本量普遍比较大,经典试验方案为(60, 5),文献[2]的方案为(38, 0),文献[20]的方案为(34, 0),在相同约束条件下,与Bayes试验方案相比,本文试验方案所需样本量平均减少61.1%。
4) 双方风险。对于基于可信度的SPOT方案,试验方案(14, 0)即可在满足双方风险的条件下做出接收判决,而对于表 4中的传统定数抽样方案,若采用(14, 0)方案,其使用方风险分别高达0.51、0.36、0.27,均无法满足双方风险要求。
5 结论1) 在分析多源先验信息的基础上,定义了分系统贡献率,并基于信息量等效原则进行数据折算,为获取系统级先验分布提供了有效途径。
2) 采用TOPSIS-AHP方法计算先验相似度P(H0)能从先验信息来源的角度给出合理的计算可信度的方法,得到的混合先验分布更加准确。
3) 从先验信息获取、先验分布权值、试验样本量、双方风险4个方面进行对比分析,表明该试验方案能够有效解决系统级先验信息获取难的问题,确定的先验分布权值更合理,而且可以有效降低试验样本量和双方风险。
[1] |
黄晓晴, 缪永生, 于盛林. 可测试性在光电雷达检测设备中的应用[J]. 哈尔滨工业大学学报, 2009, 41(3): 161. HUANG Xiaoqing, MIAO Yongsheng, YU Shenglin. Application of testability design in the equipment for detection of synthetic photovoltaic radar parameters[J]. Journal of Harbin Institute of Technology, 2009, 41(3): 161. DOI:10.3321/j.issn:0367-6234.2009.03.039 |
[2] |
张西山, 黄考利, 闫鹏程, 等. 基于验前信息的测试性验证试验方案确定方法[J]. 北京航空航天大学学报, 2015, 41(8): 1505. ZHANG Xishan, HUANG Kaoli, YAN Pengcheng, et al. Method of confirming testability verification test scheme based on prior information[J]. Journal of Beijing University of Aeronautics and Astronautics, 2015, 41(8): 1505. DOI:10.13700/j.bh.1001-5965.2014.0628 |
[3] |
汤巍, 景博, 黄以锋. 小子样变总体下的Bayes测试性验证方法[J]. 系统工程与电子技术, 2014, 36(12): 2566. TANG Wei, JING Bo, HUANG Yifeng. Testability verification test method based on Bayes theory under small sample and varying population circumstance[J]. Systems Engineering and Electronics, 2014, 36(12): 2566. DOI:10.3969/j.issn.1001-506X.2014.12.37 |
[4] |
张西山, 黄考利, 闫鹏程, 等. 基于混合验前分布的复杂装备测试性评估[J]. 振动、测试与诊断, 2015, 35(4): 697. ZHANG Xishan, HUANG Kaoli, YAN Pengchenng, et al. Evaluation of complex equipment testability based on mixed prior distribution[J]. Journal of Vibration, Measurement & Diagnosis, 2015, 35(4): 697. DOI:10.16450/j.cnki.issn.1004-6801.2015.04.016 |
[5] |
KLEYNER A, ELMORE D, BOUKAI B. A Bayesian approach to determine test sample size requirements for reliability demonstration retesting after product design change[J]. Quality Engineering, 2015, 27(3): 290. DOI:10.1080/08982112.2014.990035 |
[6] |
王京, 李天梅, 何华锋, 等. 多源测试性综合评估数据等效折合模型与方法研究[J]. 兵工学报, 2017, 38(1): 151. WANG Jing, LI Tianmei, HE Huafeng, et al. Research on multi-source data equivalent methods for testability integrated evaluation[J]. Acta Armamentarii, 2017, 38(1): 151. DOI:10.3969/j.issn.1000-1093.2017.01.020 |
[7] |
MYHRE J, JESKE D R, LI Jun, et al. Combining binomial test data via two-stage solutions[J]. Applied Stochastic Models in Business and Industry, 2018, 34(1): 22. DOI:10.1002/asmb.2255 |
[8] |
张先超, 马亚辉. 体系能力模型与装备体系贡献率测度方法[J]. 系统工程与电子技术, 2019, 41(4): 843. ZHANG Xianchao, MA Yahui. Capability model of combat system of systems and measurement method of armament contribution to combat system of systems[J]. Systems Engineering and Electronics, 2019, 41(4): 844. DOI:10.3969/j.issn.1001-506X.2019.04.20 |
[9] |
罗承昆, 陈云翔, 张洋铭, 等. 基于混合参数证据网络的装备体系贡献率评估[J]. 兵工学报, 2018, 39(12): 2490. LUO Chengkun, CHEN Yunxiang, ZHANG Yangming, et al. Evaluation of equipment contribution rate to system-of-systems based on hybrid parameter evidential network[J]. Acta Armamentarii, 2018, 39(12): 2490. DOI:10.3969/j.issn.1000-1093.2018.12.023 |
[10] |
QIU Wenhao, LIAN Guangyao, XUE Mingxi, et al. Physics of failure-based failure mode effects, and critically analysis for Integrated Circuits[J]. Systems Engineering, 2018, 21(6): 518. DOI:10.1002/sys.21451 |
[11] |
孙有朝. 指数型单元混联系统可靠性综合的信息论方法[J]. 航空动力学报, 2005, 20(5): 784. SUN Youchao. Information theory method of reliability synthesis for series-parallel or parallel-series system with hierarchical exponential units[J]. Journal of Aerospace Power, 2005, 20(5): 784. DOI:10.13224/j.cnki.jasp.2005.05.014 |
[12] |
WANG Chao, QIU Jing, LIU Guanjun, et al. Testability demonstration with component level data from virtual and physical tests[J]. Proceedings of the Institution of Mechanical Engineering, Part C: Journal of Mechanical Engineering and Science, 2015, 229(2): 267. DOI:10.1177/0954406214532909 |
[13] |
XU Weihua, LI Mengmeng, WANG Xizhao. Information fusion based on information entropy in fuzzy multi-source incomplete information system[J]. International Journal of Fuzzy Systems, 2016, 19(4): 4. DOI:10.1007/s40815-016-0230-9 |
[14] |
李想, 黄胜, 张修远. 应用改进TOPSIS法的舰船气流场方案评价[J]. 哈尔滨工业大学学报, 2016, 48(4): 137. LI Xiang, HUANG Sheng, ZHANG Xiuyuan. Evaluation of airflow on warship deck based on improved TOPSIS method[J]. Journal of Harbin Institute of Technology, 2016, 48(4): 137. DOI:10.11918/j.issn.0367-6234.2016.04.022 |
[15] |
胡仕成, 徐永东. QFD和AHP相结合的产品质量多维评估模型[J]. 哈尔滨工业大学学报, 2014, 46(11): 66. HU Shicheng, XU Yongdong. A QFD and AHP combined multiple dimensional evaluation model for product quality[J]. Journal of Harbin Institute of Technology, 2014, 46(11): 66. DOI:10.11918/j.issn.0367-6234.2014.11.011 |
[16] |
AFSHARI R, CILDEH B S. Modified sequential sampling plan using fuzzy SPRT[C]//Proceedings of the 5th Iranian Joint Congress on Fuzzy & Intelligent Systems (CFIS). Qazvin, Iran: IEEE, 2017: 116. DOI: 10.1109/CFIS.2017.8003668
|
[17] |
WANG Chao, QIU Jing, LIU Guanjun, et al. Testability demonstration test planning based on sequential posterior odds test method[J]. Journal of Risk and Reliability, 2014, 228(2): 198. DOI:10.1177/1748006X13508108 |
[18] |
MAO Zhaoyong, SONG Baowei. The sampling inspection method based on sequential posterior odd test[C]//Proceedings of the International Conference on Computational and Information Science. Chengdu: IEEE, 2010: 164. DOI: 10.1109/ICCIS.2010.46
|
[19] |
VANDERHORN E, MAHADEVAN S. Bayesian model updating with summarized statistical and reliability data[J]. Reliability Engineering & System Safety, 2018, 172: 14. DOI:10.1016/j.ress.2017.11.023 |
[20] |
周奎, 孙世岩, 严平. 基于后验风险确定故障样本量的Bayes方法[J]. 系统工程与电子技术, 2019, 41(7): 1672. ZHOU Kui, SUN Shiyan, YAN Ping. Bayes method for determining fault sample size based on posterior risk[J]. Systems Engineering and Electronics, 2019, 41. DOI:10.3969/j.issn.1001-506X.2019.07.32 |