谷廣宇 劉建敏 喬新勇
摘要:??針對(duì)小樣本條件下單純采用使用期無(wú)法準(zhǔn)確衡量裝甲裝備技術(shù)狀況的問(wèn)題,本文以裝甲車輛發(fā)動(dòng)機(jī)狀態(tài)評(píng)估為例,分析了小樣本條件下,裝備各技術(shù)狀況等級(jí)特征參數(shù)分布及其有效性和適用性。同時(shí),采用Bayes決策法,從裝備屬于各技術(shù)狀況等級(jí)的概率和錯(cuò)判風(fēng)險(xiǎn)兩方面出發(fā),對(duì)裝備狀態(tài)進(jìn)行評(píng)估。評(píng)估結(jié)果表明,利用兩種Bayes方法均能準(zhǔn)確評(píng)估發(fā)動(dòng)機(jī)技術(shù)狀況,評(píng)估準(zhǔn)確率達(dá)90%以上,并且基于最小風(fēng)險(xiǎn)的Bayes評(píng)估方法具有更好評(píng)估效果,評(píng)估準(zhǔn)確率較高,符合裝甲車輛技術(shù)狀況隨使用時(shí)間增加而逐漸變差的一般規(guī)律。該方法在裝備狀態(tài)評(píng)估中具有實(shí)用性及可行性。
關(guān)鍵詞:??狀態(tài)評(píng)估;?小樣本;?Bayes決策法;?發(fā)動(dòng)機(jī)
中圖分類號(hào):?TK428;?E923.1?文獻(xiàn)標(biāo)識(shí)碼:?A
收稿日期:?2019-03-25;?修回日期:?2019-06-27
基金項(xiàng)目:??國(guó)家自然科學(xué)基金資助項(xiàng)目(U1836101)
作者簡(jiǎn)介:??谷廣宇(1989-),男,漢族,博士研究生,主要研究方向?yàn)榘l(fā)動(dòng)機(jī)狀態(tài)評(píng)估與預(yù)測(cè)。
裝甲裝備在使用過(guò)程中,其技術(shù)狀況將隨著使用時(shí)間的增長(zhǎng)發(fā)生變化,且呈現(xiàn)劣化趨勢(shì)。雖然其技術(shù)狀況在一定程度上可以采用運(yùn)行時(shí)間來(lái)衡量[1],但由于使用環(huán)境、使用強(qiáng)度及工作地點(diǎn)等因素不盡相同,該方法有時(shí)不能完全真實(shí)地反映裝備的實(shí)際技術(shù)狀況。因此,需要運(yùn)用科學(xué)的方法對(duì)裝甲裝備進(jìn)行技術(shù)狀況評(píng)估[2]。通常情況下,在裝備的狀態(tài)檢測(cè)評(píng)估過(guò)程中,首先要采集裝備狀態(tài)信息,并對(duì)數(shù)據(jù)樣本進(jìn)行統(tǒng)計(jì)分析,才能得到相關(guān)結(jié)論。目前,基于數(shù)據(jù)驅(qū)動(dòng)的評(píng)估方法主要有模糊聚類[3]、支持向量機(jī)[4]、神經(jīng)網(wǎng)絡(luò)[5]、Bayes決策法[6]等。這些方法就統(tǒng)計(jì)分析而言,均需要盡可能大的樣本量,然而在實(shí)際評(píng)估過(guò)程中,受試驗(yàn)成本、試驗(yàn)周期等條件影響,無(wú)法進(jìn)行大量試驗(yàn)來(lái)獲取大樣本。在此情況下,仍采用傳統(tǒng)大樣本條件下的統(tǒng)計(jì)方法,難以保證其最終得到評(píng)估結(jié)果的可靠性。因此,研究在小樣本條件下的狀態(tài)檢測(cè)評(píng)估方法具有實(shí)際意義[7]。另外由于在狀態(tài)信息采集測(cè)量過(guò)程中,用以描述各個(gè)技術(shù)狀況的特征參數(shù)都是隨機(jī)變量,同一技術(shù)狀況下所有采集測(cè)量的樣本只可能是相似的,不可能完全一樣,而且某一特征值可能不只出現(xiàn)在某一特定狀況下,在其他狀況中也可能出現(xiàn)。鑒于Bayes決策理論既能參考各類狀態(tài)的分類錯(cuò)誤率,又能考慮分類錯(cuò)誤的風(fēng)險(xiǎn)程度,具有較強(qiáng)的分類決策能力。因此,本文采用Bayes決策法對(duì)裝備技術(shù)狀況進(jìn)行評(píng)估,并利用小樣本統(tǒng)計(jì)方法構(gòu)造各技術(shù)狀況下特征參數(shù)的分布,以獲取Bayes決策中的類條件概率密度。該研究可以對(duì)即將變“差”的車輛提前進(jìn)行預(yù)防性維修,降低故障風(fēng)險(xiǎn)。
1?小樣本統(tǒng)計(jì)原理與方法
為提高統(tǒng)計(jì)分析精度,在小樣本條件下,最有效的解決方法是設(shè)法增加信息量。在此類方法中,以美國(guó)斯坦福大學(xué)B.Efron教授提出的Bootstrap法[8]及B.R.Donald提出的Bayes?Bootstrap法[9]最為簡(jiǎn)單有效,而且在工程應(yīng)用中取得了豐碩成果[10-13]。其基本思路都是依據(jù)現(xiàn)有小樣本模仿未知分布,進(jìn)而獲取再生樣本,將小樣本問(wèn)題轉(zhuǎn)化成大樣本。因此,該方法適用于本文對(duì)不同技術(shù)狀況下小樣本狀態(tài)信息的統(tǒng)計(jì)分析。
1.1?Bootstrap方法
Bootstrap方法又稱自助估計(jì)法,其原理主要根據(jù)觀測(cè)到來(lái)自于未知總體分布F的隨機(jī)子樣X(jué)=(X1,…,Xn),估計(jì)總體分布F的某一分布特征R(X,F(xiàn)),如均值、方差等,從而推測(cè)總體分布F。具體方法如下:
設(shè)總體分布F的某個(gè)分布特征θ=θ(F)(如均值,方差等),由觀測(cè)子樣X(jué)=(X1,…,Xn)構(gòu)造經(jīng)驗(yàn)分布Fn,則有對(duì)θ的估計(jì)=(Fn),估計(jì)誤差為
Tn=(Fn)-(F)(1)
根據(jù)經(jīng)驗(yàn)分布Fn,重新抽取再生子樣X(jué)(1)=(X(1)1,…,X(1)n),進(jìn)而構(gòu)造經(jīng)驗(yàn)分布函數(shù)F(1)n。于是由X(1)又可得到θ的估計(jì)(F(1)n)。此時(shí),可得到估計(jì)誤差Tn的Bootstrap統(tǒng)計(jì)量R(1)n,即
R(1)n=(F(1)n)-(Fn)(2)
重復(fù)抽取多組再生子樣X(jué)(i),i=1,2,…,m,可計(jì)算相應(yīng)R(i)n,i=1,2,…,m,并利用R(i)n的分布去逼近Tn的分布,并根據(jù)式(1)得到θ(F)的樣本,即
θ(i)(F)=(Fn)-Tn≌(Fn)-R(i)n(3)
該方法與傳統(tǒng)統(tǒng)計(jì)方法相比,在小樣本估計(jì)中具有較高精度。
1.2?Bayes?Bootstrap方法
Bayes?Bootstrap方法,即隨機(jī)加權(quán)自助法,也是將子樣信息“提攜”的方法,是對(duì)Bootstrap方法的改進(jìn)。其原理主要利用Dirichilet隨機(jī)向量對(duì)子樣X(jué)=(X1,…,Xn)進(jìn)行加權(quán),從而計(jì)算隨機(jī)加權(quán)統(tǒng)計(jì)量Dn,進(jìn)而利用Dn模仿估計(jì)誤差Tn的分布。具體方法如下:
設(shè)子樣X(jué)=(X1,…,Xn)來(lái)自于未知總體分布F,μ和σ2為F的未知期望和方差,和S2為X的統(tǒng)計(jì)期望和方差,則估計(jì)誤差為
T1=-μ(4)
T2=nn-1S2-σ2(5)
式中,T1和T2為期望的估計(jì)誤差和方差的估計(jì)誤差。生成N個(gè)參數(shù)為(1,…,1)的Dirichilet隨機(jī)向量Vi=(vi1,vi2,…,vin),i=1,2,…,N,可獲取隨機(jī)加權(quán)統(tǒng)計(jì)量為
Di1=∑nk=1vikxk-(6)
Di2=nn-1∑nk=1vik(xk-)2-nn-1S2(7)
利用Di1和Di2模仿估計(jì)誤差T1和T2的分布,根據(jù)式(4)和式(5),可得
μ=-T1≈-1(8)
S2=nn-1S2-T1≈nn-1S2-2(9)
已有成果的實(shí)際應(yīng)用及仿真計(jì)算表明,該方法分析精度較高于Bootstrap方法[14]。
2?Bayes決策理論
Bayes決策理論是模式識(shí)別問(wèn)題中經(jīng)典的分類決策理論,其基本思想是通過(guò)各狀態(tài)下的先驗(yàn)概率分布推定后驗(yàn)概率[15]。根據(jù)決策規(guī)則,Bayes決策又分為基于最小錯(cuò)誤率[16]和基于最小風(fēng)險(xiǎn)的Bayes決策[17]。
2.1?基于最小錯(cuò)誤率的貝葉斯決策
假設(shè)共有m個(gè)決策狀態(tài)總體W=(w1,w2,…,wm),已知其先驗(yàn)概率P(wi),i=1,2,…,m,依據(jù)前文小樣本統(tǒng)計(jì)理論可推定各狀態(tài)下特征參數(shù)x的分布模型,并獲得其類條件概率密度P(x|wi)。則根據(jù)Bayes公式,可得各狀態(tài)后驗(yàn)概率為
P(wi|x)=P(x|wi)P(wi)∑mi=1P(x|wi)P(wi)(10)
并有決策規(guī)則:若P(wi|x)=max\[P(wi|x)\],則x∈wi。此類決策規(guī)則稱為基于最小錯(cuò)誤率的貝葉斯決策。
2.2?基于最小風(fēng)險(xiǎn)的貝葉斯決策
實(shí)際決策過(guò)程中除了錯(cuò)誤率外,有時(shí)也需要考慮決策錯(cuò)誤所帶來(lái)的風(fēng)險(xiǎn),在這種情況下,通??梢赃x擇基于最小風(fēng)險(xiǎn)的貝葉斯決策方法[18]。
基于最小風(fēng)險(xiǎn)的貝葉斯方法在基于最小錯(cuò)誤方法上引入了“損失”的概念,即把屬于wi的狀態(tài)錯(cuò)判為aj狀態(tài)的損失,可記為cji,一般決策損失表如表1所示。
根據(jù)表1及前文所得后驗(yàn)概率,可計(jì)算決策為wi的狀態(tài)條件風(fēng)險(xiǎn)為
R(wi|x)=E\[c(i|j)\]=∑mj=1c(j|i)P(wj|x)(11)
并有決策規(guī)則:若R(wi|x)=min\[R(wi|x)\],則x∈wi。
3?應(yīng)用實(shí)例
基于貝葉斯方法的裝備狀態(tài)評(píng)估,實(shí)質(zhì)上就是將裝備的技術(shù)狀況劃分為從優(yōu)到劣幾個(gè)等級(jí)的問(wèn)題,本文以裝甲車輛發(fā)動(dòng)機(jī)狀態(tài)評(píng)估為例,并依據(jù)摩托小時(shí)(使用時(shí)間)將其技術(shù)狀況分為“優(yōu)”、“良”、“一般”、“較差”、“差”這5個(gè)等級(jí),分別以“0~200”,“200~300”,“300~400”,“400~500”,“500~600”這5個(gè)階段統(tǒng)計(jì)特征表示,即有W=(w1,w2,w3,w4,w5)={“優(yōu)”,“良”,“一般”,“較差”,“差”}。
3.1?特征提取
評(píng)估發(fā)動(dòng)機(jī)技術(shù)狀況的參數(shù)很多,但實(shí)踐證明,通過(guò)減速時(shí)間評(píng)估技術(shù)狀況是一種簡(jiǎn)單有效的辦法[19],因此本文選取柴油機(jī)減速時(shí)間作為技術(shù)狀況評(píng)估的特征參數(shù)。減速時(shí)間是指發(fā)動(dòng)機(jī)在最高空轉(zhuǎn)轉(zhuǎn)速下運(yùn)轉(zhuǎn)時(shí)迅速停止供油,使其自行減速,計(jì)算從最高空轉(zhuǎn)轉(zhuǎn)速減到轉(zhuǎn)速為0的時(shí)間。隨著使用時(shí)間增長(zhǎng),發(fā)動(dòng)機(jī)活塞缸套磨損增加,阻力減小,減速時(shí)間將變長(zhǎng)[20]。
發(fā)動(dòng)機(jī)瞬時(shí)轉(zhuǎn)速如圖1所示,將實(shí)際計(jì)算瞬時(shí)轉(zhuǎn)速?gòu)膎3減速到n4的時(shí)間作為發(fā)動(dòng)機(jī)減速時(shí)間,此時(shí)取n3=2?200?r/min,n4=800?r/min。
依此方法在5種技術(shù)狀況等級(jí)下,分別對(duì)同一型號(hào)6輛不同摩托小時(shí)的裝甲車輛進(jìn)行測(cè)試,提取減速時(shí)間,并作為統(tǒng)計(jì)樣本,用以計(jì)算該等級(jí)下分布特征,以檢驗(yàn)本文評(píng)估方法的準(zhǔn)確性。另外,隨機(jī)抽取5輛未知狀況的裝甲車,提取修正減速時(shí)間作為測(cè)試樣本,利用本文方法對(duì)其進(jìn)行評(píng)估。統(tǒng)計(jì)樣本減速時(shí)間如表2所示,測(cè)試樣本減速時(shí)間如表3所示。
3.2?小樣本統(tǒng)計(jì)
假設(shè)發(fā)動(dòng)機(jī)減速時(shí)間在各個(gè)使用階段均服從正態(tài)分布,利用Bayes?Bootstrap小樣本統(tǒng)計(jì)方法,對(duì)各個(gè)階段減速時(shí)間進(jìn)行統(tǒng)計(jì),根據(jù)式(4)~式(9),得各階段分布的均值和標(biāo)準(zhǔn)差,Bayes?Bootstrap法下各技術(shù)狀況分布特征如表4所示,各技術(shù)狀況等級(jí)分布特征如圖2所示。由表4和圖2可以看出,減速時(shí)間隨技術(shù)狀況等級(jí)變化較為明顯,因此能夠較好地反映發(fā)動(dòng)機(jī)技術(shù)狀況的變化。
3.3?基于Bayes決策的狀態(tài)評(píng)估
由于測(cè)試車輛摩托小時(shí)在0~600?h內(nèi)呈均勻分布,因此各狀態(tài)的先驗(yàn)概率為
P(W)=\[1/3,1/6,1/6,1/6,1/6\]
考慮評(píng)估過(guò)程中錯(cuò)誤率的存在,根據(jù)經(jīng)驗(yàn),引入發(fā)動(dòng)機(jī)狀態(tài)評(píng)估決策損失C,發(fā)動(dòng)機(jī)狀態(tài)評(píng)估決策損失如表5所示。
3.4?評(píng)估準(zhǔn)確性對(duì)比
根據(jù)表4及式(10)和式(11),分別利用兩種Bayes方法對(duì)統(tǒng)計(jì)樣本進(jìn)行評(píng)估,統(tǒng)計(jì)樣本評(píng)估結(jié)果如表6所示。
由表6可以看出,利用兩種Bayes方法均能夠準(zhǔn)確評(píng)估發(fā)動(dòng)機(jī)技術(shù)狀況,評(píng)估準(zhǔn)確率達(dá)90%以上,并且基于最小風(fēng)險(xiǎn)的Bayes評(píng)估方法具有更好評(píng)估效果。
3.5?測(cè)試樣本狀態(tài)評(píng)估
采用基于最小風(fēng)險(xiǎn)的Bayes評(píng)估方法,對(duì)表3測(cè)試樣本進(jìn)行評(píng)估,在5種狀態(tài)下,各測(cè)試樣本的條件風(fēng)險(xiǎn)及評(píng)估結(jié)果如表7所示。
由表7可以看出,評(píng)估結(jié)果符合裝甲車輛技術(shù)狀況隨使用時(shí)間增加而逐漸變差的一般規(guī)律,證明該方法在裝備狀態(tài)評(píng)估中的實(shí)用性及可行性。
該方法還能判斷測(cè)試樣本的變化趨勢(shì),例如測(cè)試樣本1根據(jù)表7被評(píng)估為“優(yōu)”狀態(tài),但其“良”狀態(tài)的風(fēng)險(xiǎn)與“優(yōu)”狀態(tài)只有很小差別。因此,可判斷該樣本接近“優(yōu)”狀態(tài)與“良”狀態(tài)的邊緣,即將退化為“良”狀態(tài)。依據(jù)該方法可以對(duì)即將變“差”的車輛提前進(jìn)行預(yù)防性維修,以降低故障風(fēng)險(xiǎn)。
4?結(jié)束語(yǔ)
本文通過(guò)結(jié)合Bootstrap小子樣統(tǒng)計(jì)方法及Bayes決策法,有效地解決了小樣本條件下的裝甲裝備狀態(tài)評(píng)估的問(wèn)題,并以裝甲車輛發(fā)動(dòng)機(jī)為例,驗(yàn)證了該方法的有效性。采用Bootstrap小子樣統(tǒng)計(jì)方法,能夠有效擴(kuò)充發(fā)動(dòng)機(jī)狀態(tài)參數(shù)樣本,實(shí)現(xiàn)小樣本條件下的裝備狀態(tài)評(píng)估的先驗(yàn)統(tǒng)計(jì);Bayes決策法在裝備狀態(tài)評(píng)估過(guò)程中,不但考慮了各技術(shù)狀況出現(xiàn)的概率,還考慮了錯(cuò)誤評(píng)估的風(fēng)險(xiǎn),比其他技術(shù)狀況評(píng)估方法更具合理性;采用本文方法進(jìn)行裝甲車輛發(fā)動(dòng)機(jī)的技術(shù)狀況等級(jí)評(píng)估,評(píng)估結(jié)果較符合實(shí)際規(guī)律,正確率可達(dá)90%以上。該研究具有一定工程實(shí)用性及可行性。
參考文獻(xiàn):
[1]?李國(guó)錦.?基于層次分析及模糊綜合評(píng)判的裝甲車輛發(fā)動(dòng)機(jī)狀態(tài)評(píng)估[J].?車輛與動(dòng)力技術(shù),?2017(4):?40-44.
[2]?張會(huì)奇,?陳春良,?劉峻巖,?等.?基于BP神經(jīng)網(wǎng)絡(luò)的裝甲車輛發(fā)動(dòng)機(jī)使用狀態(tài)評(píng)價(jià)[J].?兵工自動(dòng)化,?2014,?33(1):?32-34.
[3]?辛江慧,?鄒政耀,?任成龍.??基于灰色定權(quán)聚類的混合動(dòng)力汽車發(fā)動(dòng)機(jī)運(yùn)行狀態(tài)評(píng)估[J].??重慶理工大學(xué)學(xué)報(bào):自然科學(xué),?2014,?28(7):?26-29.
[4]?汪瑾,?陳果,?王洪偉,?等.??基于后驗(yàn)概率SVM的航空發(fā)動(dòng)機(jī)滾動(dòng)軸承狀態(tài)評(píng)估[J].??軸承,?2015(5):?55-59.
[5]?劉伯運(yùn),?歐陽(yáng)光耀,?常漢寶.??基于神經(jīng)網(wǎng)絡(luò)和DS證據(jù)理論的柴油機(jī)狀態(tài)評(píng)估[J].??車用發(fā)動(dòng)機(jī),?2005(5):?14-17.
[6]?朱偉明.??基于貝葉斯融合的柴油機(jī)性能狀態(tài)評(píng)估方法研究[J].??內(nèi)燃機(jī),?2010(2):?18-19,?36.
[7]?高攀東,?陳曉陽(yáng),?沈雪瑾,?等.??小樣本下軸承可靠性評(píng)估方法對(duì)比分析[J].??軸承,?2014(12):??44-48.
[8]?Efron?B,?Tibshirani?R?J.??An?introduction?to?the?bootstrap[M].??London:?Chapman?and?Hall,?1993.
[9]?Donald?B?R.??The?bayesian?bootstrap[J].??The?Annals?of?Statistics,?1981,?9(1):?130-134.
[10]?李會(huì)鵬,?明廷鋒,?賀國(guó).??小樣本條件下一種旋轉(zhuǎn)機(jī)械振動(dòng)狀態(tài)閾值的確定方法及應(yīng)用[J].??武漢理工大學(xué)學(xué)報(bào):?交通科學(xué)與工程版,?2012,?36(5):?1059-1063.
[11]?喬新勇,?劉建敏,?張小明.??小子樣統(tǒng)計(jì)理論在裝備狀態(tài)檢測(cè)中的應(yīng)用[J].??裝甲兵工程學(xué)院學(xué)報(bào),?2009,?23(2):?41-43,?87.
[12]?唐雪梅,?周伯昭,?李榮.??武器裝備小子樣綜合試驗(yàn)設(shè)計(jì)與鑒定技術(shù)[J].??戰(zhàn)術(shù)導(dǎo)彈技術(shù),?2007(2):?51-56,?79.
[13]?袁?鵬,?張衛(wèi)兵,?韓自剛.?銀川地區(qū)建筑地基沉降計(jì)算驗(yàn)系數(shù)取值[J].?科學(xué)技術(shù)與工程,?2017,?17(20):??236-240.
[14]?樊鶴紅,?劉盼,?趙興群,?等.??小樣本條件下行波管可靠性評(píng)估方法的研究[J].??電子學(xué)報(bào),?2010,?38(6):?1394-1398.
[15]?張湘平,?張金槐,?謝紅衛(wèi),?等.??基于隨機(jī)加權(quán)法的BAYES精度評(píng)定[J].??國(guó)防科技大學(xué)學(xué)報(bào),?2001,?23(3):?98-102.
[16]?全星慧,?牟海維,?呂秀麗,?等.??基于PCA與貝葉斯決策的人臉識(shí)別算法[J].?光學(xué)儀器,?2014,?36(2):?122-125.
[17]?徐智康,?李旸,?李德玉.?基于可變最小貝葉斯風(fēng)險(xiǎn)的層次多標(biāo)簽分類方法[J].?南京大學(xué)學(xué)報(bào):?自然科學(xué)版,?2017,?53(6):?1023-1032.
[18]?毛思敏.??基于最小風(fēng)險(xiǎn)貝葉斯準(zhǔn)則的工業(yè)過(guò)程故障診斷[D].?武漢:華中科技大學(xué),?2016.
[19]?張更云,?何來(lái)龍.??一種評(píng)價(jià)發(fā)動(dòng)機(jī)實(shí)際技術(shù)狀況的方法[J].??裝甲兵工程學(xué)院學(xué)報(bào),?2005,?19(4):?67-70.
[20]?張小明,?劉建敏,?喬新勇,?等.?柴油機(jī)健康狀況評(píng)估及其方法[J].?裝甲兵工程學(xué)院學(xué)報(bào),?2008,?22(2):?30-34.
Equipment?States?Evaluation?Based?on?Bayes?Method?with?Small?Sample
GU?Guangyu,?LIU?Jianmin,?QIAO?Xinyong
(Vehicle?Engineering?Department,?Army?Academy?of?Armored?Forces,?Beijing?100072,?China)
Abstract:??In?order?to?solve?the?problem?that?the?technical?condition?of?armored?equipment?can?not?be?accurately?measured?by?using?the?service?life?alone?under?the?condition?of?small?sample,?this?paper?takes?the?engine?state?evaluation?of?armored?vehicle?as?an?example,?and?analyzes?the?effectiveness?and?applicability?of?the?characteristic?parameters?distribution?in?each?technical?condition?of?the?equipment?under?the?condition?of?small?sample.??At?the?same?time,?the?Bayes?decision?method?is?used?to?evaluate?the?equipment?status?from?two?aspects:??the?probability?of?the?equipment?belonging?to?each?technical?condition?grade?and?the?misjudgment?risk.??The?evaluation?results?show?that?both?Bayes?methods?can?accurately?evaluate?the?engine?technical?status,?and?the?evaluation?accuracy?is?larger?than?90%.??The?Bayes?evaluation?method?based?on?minimum?risk?has?better?evaluation?effect?and?higher?evaluation?accuracy,?which?accords?with?the?general?law?that?the?technical?condition?of?armored?vehicles?becomes?worse?with?the?increase?of?service?time.??This?method?is?practical?and?feasible?in?equipment?condition?evaluation.
Key?words:??states?evaluation;?small?sample;?Bayes?decision;?engine