时继庆,王书宁
(清华大学 自动化系,100084 北京 ,shijq04@mails.tsinghua.edu.cn)
基于飞行可靠度的海防战术导弹抽样检验方案
时继庆,王书宁
(清华大学 自动化系,100084 北京 ,shijq04@mails.tsinghua.edu.cn)
为了解决目前海防战术导弹批生产抽样检验中以导弹命中概率为主要检验指标,导致检验指标与检验目的不符、不能对生产质量一致性进行检验的问题,提出了以飞行可靠度为批生产抽样检验主检验指标,研究基于指数分布下飞行可靠度的批生产抽样检验,并给出抽样检验方案、一致性检验和批生产抽样检验的验收准则.模拟计算结果表明,在同等条件下,该方案比原方案检验效率更高,试验样本更少,证实了方法的可行性和有效性.
飞行可靠度;抽样检验;战术导弹;一致性
海防战术导弹批生产抽样检验的2个目的是:1)检验批生产导弹的质量;2)检验批生产导弹质量的一致性.目前,海防战术导弹的批生产抽样检验采用计数抽样检验方案,以导弹命中概率作为主检验指标,以导弹是否命中目标作为判定导弹是否合格的依据.导弹命中概率是导弹命中率和飞行可靠度的积,其中飞行可靠度的指标远大于命中率.导弹飞行可靠度是反映导弹批生产质量的主要指标.
基于飞行可靠度的抽样检验,就是将主检验指标由导弹命中概率换成导弹飞行可靠度,并将导弹飞行可靠度转换为导弹可靠飞行时间,即对导弹飞行寿命进行抽样检验.
对寿命的抽样检验一般采用截尾试验,包括定数截尾和定时截尾寿命试验[1].对于对寿命有明确要求的产品,一般采用定时截尾试验.许多学者[2-9]在这方面做了大量的工作.本文以导弹飞行可靠度为主检验目标,推导出基于可靠飞行时间的抽样检验方案和质量一致性检验方法,提出以抽样数、最大容许失效数和飞行试验截止时间为标志的抽样检验方案,给出包括一致性检验在内的验收标准和确定试验截止时间的原则.最后,通过算例,将本文提出的方法与原抽样检验方法进行了比较,结果表明,本文方法与原方法相比,具有样本量小、检验充分、检验费用低的特点,解决了原方法不能进行批生产质量一致性检验的问题.
设导弹的可靠飞行时间t服从参数为θ的指数分布,其密度函数为
式中θ为导弹的平均可靠飞行寿命(θ>0).其分布函数为
式(2)是导弹在( 0,t]时间内失效的概率,在( 0,t]内不失效的概率是1-F(t).
则飞行可靠度与可靠飞行时间的关系为
在导弹飞行任务中,最大有效航程攻击时对导弹的可靠性要求最高,因此导弹的飞行可靠度是导弹最大有效航程任务剖面内完成作战任务的最低要求.
设导弹最大有效射程下的飞行时间为tmax,导弹飞行可靠度设计指标为R,则
式中θs称为导弹设计可靠飞行时间,也是导弹在飞行可靠度R下可靠飞行时间的期望值.
导弹可靠飞行时间的方差为
抽取n发导弹进行飞行试验,在( 0,t]内失效k发的概率为
在( 0,t]内导弹失效数小于事先给定c的概率为
式中:α为弃真概率(生产方风险);β为采伪概率(使用方风险).
此即为导弹在(n,c)抽样检验方案中被接受的概率.
由方程求得
式(10)确定的(n,c),就是基于可靠飞行寿命的抽样检验方案.
设抽取n枚导弹进行截止时间为t0的飞行试验,试验均在t0时刻结束,在( 0,t0]内飞行失败数为k,这 k枚导弹的飞行时间分别为(x1,x2,…,xk),记(x1,…,xr,k=r)的联合密度为 f(x1,…,xr,r;θ),则
于是,对r≥ 1,由式(11)可得θ的MLE为
式中Tr为n发导弹总的飞行试验时间.由于当r=0时,θ的MLE估计应为=+∞,可采用文献[10]的建议得到:
批生产质量的一致性检验,就是检验试验样本中反映导弹质量特征数据的散布(方差)是否在容许的范围内.因此,进行批生产质量一致性检验,就是进行可靠飞行时间散布(方差)的检验.
文献[3]给出了样本寿命的期望值和方差,对于海防战术导弹,就是可靠飞行时间的期望值和方差为
式中θ为导弹可靠飞行时间的期望值,这里用设计值θs代替,即样本方差的标准值为
式中t0为试验截止时间.
如果 S2≤ var[Sk(t0;θs)],则认为批生产质量的一致性满足要求;否则,为不满足要求.
事先确定:θ0,θ1,α,β,t0.
其中:θ0为可接受平均寿命;θ1为极限平均寿命.
检验规则为:
1)如果 k≤ c,且 S2≤ var[Sk(t0;θs)],接受该批产品;
2)如果 k > c,或 S2> var[Sk(t0;θs)],拒绝该批产品.
制订满足上述检验规则的抽样检验方案,就是在给定的 θ0,θ1,α,β 条件下,求满足方程的(n,c).(n,c)就是抽样检验方案.
设导弹设计的任务飞行可靠度为R,任务飞行时间为 θs,则
基于可靠飞行时间的抽样检验非常关键的一个环节,就是确定t0,t0的取值大小,决定了检验
因此,在海防战术导弹的检验中,一般采用任务飞行时间作为实际试验截止时间.
对飞行可靠度为 0.90、任务飞行时间为1 000 s,命中概率为0.75的海防战术导弹,分别用计数抽样检验和本文的方法进行抽样检验方案设计.
已知命中概率为P=0. 75,取P0=0. 7,P1=0. 50,α=β=0. 20,得到抽样检验方案如表1所示.
表1 计数抽样检验方案
从表1可以看出,序号2、3、4的方案在抽样数、失败数和检验风险上比较适中,如果样本数有严格要求,也可以选择序号 1,但使用方要承担较大风险.
取飞行可靠性R=0. 9,经计算得θs=949 1 s;取 θ0=800 0 s、θ1=4 000 s、t0=1 000 s、α =β=0. 20,得抽样检验方案如表2所示.
表2 基于可靠飞行时间的抽样检验方案
从表2可以看出,序号2、3、4的方案在抽样数、失败数和检验风险上比较适中.
对比2种抽样检验方法,可以看出,本文提出的检验方法与计数抽样检验相比具有的优势为:
1)在同样检验条件下,所需要的试验子样数要少于计数抽样检验的子样数,可以降低检验成本;
2)检验充分.试验飞行过程覆盖了导弹的整个飞行包络线,可以充分考核导弹的飞行可靠性;
3)简化了试验组织指挥的难度,减少了测控系统的负担,降低了试验费用.
1)本文的抽样检验方案与计数抽样检验方案相比,具有检验指标准确,检验内容丰富,检验力强的特点.
2)解决了计数抽样检验中不能进行质量一致性检验的问题,完善了抽样检验内容,达到了检验的目的.
3)基于战术导弹的共同性,该方法可以应用于其他战术导弹的批生产抽样检验中.
[1]LAWLESS J F.Statistical Models and Methods for Lifetime Data[M].New York:Wiley,1982.
[2]王宏.混合截尾寿命试验情形下指数分布平均寿命的置信限[J].应用概率统计, 1996,12(4):348-354.
[3]刘宝友,王燕,茆诗松.指数分布定时截尾试验失效率近似置信区间[J].应用概率统计, 1999,15(3):234-238.
[4]吴启光,吕建华.定时截尾试验下指数分布产品可靠性抽样检验方案[J].系统科学与数学, 2003,23(2):145-154.
[5]LAM Y.An optimal single variable sampling plan with censoring[J].The Statistician, 1990,39(1):53 -66.
[6]LAM Y.Bayesian variable sampling plans for the exponential distribution with type I censoring[J].The Annals of Statistics, 1994,22(2):696 -711.
[7] LAM Y,CHOY S T B.Bayesian variable sampling plans for the exponential distribution with uniformly distributed random censoring[J].Journal of Statistical Planning and Inference, 1995,47(3):277-293.
[8]HUANG Wentao,LIN Yupin.Bayesian sampling plans for exponential distribution based on uniform random censored data[J].Computation Statistics& Data Analysis, 2004,44(4):669 -691.
[9]LIN Chientai,HUANG Yenlung,BALAKRISHNAN N.Exact Bayesian variable sampling plans for the exponential distribution based on type-I and type-II hybrid censored samples [J].Communications in Statiscs——Simulation and Computation, 2008,37:1101 -1116.
[10]BARTHOLOMEW D J.A problem in life testing[J],JournalofAmerican StatisticalAssociation, 1957,52(279):350-355.
Sampling inspection plan for coast defense tactical missile based on flying reliability
SHI Ji-qing,WANG Shu-ning
(Dept.of Automation,Tsinghua University,100084 Beijing,China,shijq04@mails.tsinghua.edu.cn)
The present sampling inspection for the coast defense tactical missile often takes the hitting probability as the main test index,which may result in the mismatch of the test index and test purpose as well as the inability to test the consistency of manufacture quality.To over come this problem,in this paper,the flying reliability of missile is then taken as the main test index.Assume the flying reliability conform to the exponential distribution,the sampling inspection plan is proposed based on this new main test index.Besides,the acceptance rules of sampling inspection and the consistency test are given.Simulation results show that under similar conditions,the proposed plan is more effective and requires fewer samples,which indicats the feasibility and validity of the new method.
flying reliability;sampling inspection;tactical missile;consistency
TP1
A
0367-6234(2011)09-0101-04
2010-01-06.
国家自然科学基金资助项目( 60674025, 60534060,60974008);国家高技术研究发展计划资助项目(2007AA04Z193);高等学校博士学科点专项科研基金资助项目(200900030029).
时继庆(1962—),男,博士研究生;
王书宁(1956—),男,教授,博士生导师.
(编辑 张 红)