李 濤，陳 麗，譚 晨

(空軍勤務(wù)學(xué)院 航空彈藥保障系，江蘇 徐州 221003)

對于“長期儲(chǔ)存、一次使用”的制導(dǎo)彈藥而言，儲(chǔ)存在整個(gè)壽命周期中占據(jù)極大比重。針對通過性能檢測的方法來檢查制導(dǎo)彈藥的狀況，檢測通過，即合格這種把制導(dǎo)彈藥狀態(tài)一分為二“合格或故障”的方法，亟需一種制導(dǎo)彈藥的性能評估方法合理區(qū)分制導(dǎo)彈藥狀態(tài)情況。范志峰等[1]通過一些退化敏感參數(shù)等來評估彈藥控制系統(tǒng)的壽命，采用了退化量分布法進(jìn)行評估，但并未對整個(gè)彈藥進(jìn)行評估；王金柱[2]利用環(huán)境因子和混合貝塔分布的貝葉斯法對制導(dǎo)彈藥進(jìn)行評估，此方法針對的是成熟型產(chǎn)品。在保障過程中產(chǎn)生的大量檢測數(shù)據(jù)、維護(hù)數(shù)據(jù)等蘊(yùn)含著豐富的制導(dǎo)彈藥狀態(tài)信息，應(yīng)用KPCA與LSSVM[3-5]根據(jù)檢測數(shù)據(jù)、維護(hù)數(shù)據(jù)等狀態(tài)信息來進(jìn)行制導(dǎo)彈藥性能評估的理論和方法是可行的[6]，不僅能克服上述兩種方法的缺點(diǎn)，還能實(shí)現(xiàn)動(dòng)態(tài)地評估制導(dǎo)彈藥性能，具有重要理論意義和軍事應(yīng)用價(jià)值。

1 KPCA特征提取方法

為了判斷制導(dǎo)彈藥的性能狀況，制導(dǎo)彈藥檢測過程中對電氣系統(tǒng)、引信控制器、綜控機(jī)等各個(gè)子系統(tǒng)都進(jìn)行檢測，每一個(gè)檢測的子系統(tǒng)中所檢測的指標(biāo)少則幾十，多則上百；再加上保障維護(hù)中的事件數(shù)據(jù)，若把這些數(shù)據(jù)用于制導(dǎo)彈藥的性能評估，將導(dǎo)致模型癱瘓[7-8]。因此，需要通過一定的參數(shù)提取方法將各種數(shù)據(jù)進(jìn)行降維，提取出體現(xiàn)制導(dǎo)彈藥儲(chǔ)存性能退化狀態(tài)的參數(shù)為后續(xù)性能退化評估打下良好的基礎(chǔ)。

核主元方法[9](Kernel Principal Component Analysis,KPCA)是近年來比較流行的非線性主元分析方法，其核心思想是通過非線性映射將輸入轉(zhuǎn)換至高維空間，再處理高維空間中對應(yīng)的線性問題，利用核函數(shù)避開繁復(fù)的映射函數(shù)選取。

(1)

設(shè)λ為矩陣C的特征值，μ為空間H的特征向量，μ滿足正交條件，將μ歸一化，設(shè)μj為矩陣C的第j個(gè)特征向量，則φ(ak)在μj投影為

(2)

則gj(ak)為測試樣本ak對應(yīng)的第j個(gè)非線性主元分量。將所有的投影值形成矢量g(ak)作為新特征向量。為簡化計(jì)算，本文使用核函數(shù)Kji=φ(aj)·φ(ai)代替H空間中的點(diǎn)積運(yùn)算，式(1)、式(2)變?yōu)?/p>

(3)

2 檢測維護(hù)數(shù)據(jù)的制導(dǎo)彈藥儲(chǔ)存性能評估方法

制導(dǎo)彈藥在儲(chǔ)存過程其性能隨時(shí)間不斷退化下降直至故障失效。儲(chǔ)存性能越優(yōu)，制導(dǎo)彈藥的技術(shù)狀態(tài)就越好，檢測維護(hù)數(shù)據(jù)大多落在標(biāo)準(zhǔn)值附近；儲(chǔ)存性能越差，制導(dǎo)彈藥越接近失效狀態(tài)，檢測維護(hù)數(shù)據(jù)大多分布在失效閾值邊緣或之外。故可以用檢測維護(hù)參數(shù)落在失效閾值內(nèi)的概率Pi(t)來衡量制導(dǎo)彈藥的儲(chǔ)存性能[6，10]。

2.1 模型構(gòu)建

設(shè)某型制導(dǎo)彈藥有n個(gè)檢測維護(hù)參數(shù)，第i(i=1,2，…,n)個(gè)檢測維護(hù)數(shù)據(jù)值、失效閾值上限、失效閾值下限分別記為xi、Ui和Li。在t時(shí)刻對制導(dǎo)彈藥進(jìn)行一次測試，其退化量xi(t)的分布函數(shù)為G(x,βi)，其中βi=(β1,β2,…,βk)為該分布的參數(shù)向量，對應(yīng)的概率密度函數(shù)g(x,βi)可表示為

(4)

記Ti為測出第i個(gè)檢測參數(shù)不合格的時(shí)間，第i個(gè)檢測參數(shù)在t時(shí)刻未失效的概率為

(5)

考慮某型制導(dǎo)彈藥具有n個(gè)檢測維護(hù)參數(shù)，任一檢測維護(hù)參數(shù)的失效均會(huì)引起制導(dǎo)彈藥故障，因此一枚正常的制導(dǎo)彈藥需要保證所有參數(shù)均處在閾值規(guī)定范圍內(nèi)。假設(shè)制導(dǎo)彈藥參數(shù)相互獨(dú)立，則制導(dǎo)彈藥在t時(shí)刻未失效的概率為

(6)

因?qū)崪y數(shù)據(jù)需進(jìn)行主成分提取，提取出的主成分已不是原有意義上的參數(shù)，故需對提取出的主成分參數(shù)重新確定閾值。基于“3σ原則”的閾值法是一種常用、性能良好的確定方法。假定制導(dǎo)彈藥的性能檢測維護(hù)數(shù)據(jù)服從正態(tài)分布，分布由均值μ和σ方差決定，故該方法適合用于確定制導(dǎo)彈藥失效閾值的上下限Ui和Li的確定。對于正態(tài)分布檢測數(shù)據(jù)落到區(qū)間(μ-3σ,μ+3σ)之外的概率為0.3%，越靠近分布邊緣概率越小。當(dāng)概率足夠小時(shí)，可以認(rèn)為該數(shù)據(jù)異常。故可據(jù)此制定健康閾值區(qū)間(μ-nσ,μ+nσ)。

2.2 數(shù)據(jù)分布類型的確定和檢驗(yàn)

假定制導(dǎo)彈藥的檢測維護(hù)數(shù)據(jù)服從正態(tài)分布[11]，對其進(jìn)行的擬合檢驗(yàn)為正態(tài)性檢驗(yàn)，可以利用樣本{x1,x2,…,xn}來檢驗(yàn)總體分布是否服從N(μ,σ2)，含有未知參數(shù)μ和σ2，則K-S檢驗(yàn)中的假設(shè)變?yōu)?/p>

(7)

檢驗(yàn)統(tǒng)計(jì)量為

(8)

(9)

2.3 數(shù)據(jù)分布參數(shù)預(yù)測

制導(dǎo)彈藥檢測參數(shù)與時(shí)間緊密相關(guān)，相應(yīng)的分布參數(shù)也隨時(shí)間變化，因此可將其看作時(shí)間序列。通過上述分析，制導(dǎo)彈藥的Pi(t)值的預(yù)測即是對檢測維護(hù)數(shù)據(jù)分布參數(shù)的預(yù)測。制導(dǎo)彈藥檢測維護(hù)數(shù)據(jù)分布參數(shù)具有小樣本、短時(shí)序、非線性等特點(diǎn)，而在處理小樣本預(yù)測和估計(jì)問題時(shí)，統(tǒng)計(jì)學(xué)習(xí)理論有較好的效果。最小二乘支持向量機(jī)(least squares support vector machine，LS-SVM)是一種建立在統(tǒng)計(jì)學(xué)習(xí)理論基礎(chǔ)上的機(jī)器學(xué)習(xí)方法，是標(biāo)準(zhǔn)支持向量機(jī)的擴(kuò)展，在解決分布參數(shù)時(shí)序小樣本的問題的同時(shí)，避免了標(biāo)準(zhǔn)SVM用于預(yù)測時(shí)計(jì)算負(fù)擔(dān)過大的問題。因此，在預(yù)測檢測維護(hù)數(shù)據(jù)分布參數(shù)時(shí)間序列時(shí)[12]，本文中采用LS-SVM預(yù)測模型。

1)模型描述

f(x)=WT(x)+b

(10)

其中:xi為i個(gè)輸入向量；yi為對應(yīng)的輸出量；l為樣本容量；Wi∈Rn,b∈R,φ(·)為核函數(shù)。

在LS-SVM中，上述回歸問題對應(yīng)的優(yōu)化問題為

(11)

約束條件為

yi=wTφ(xi)+b+ei

(12)

其中：ei為第i個(gè)估計(jì)值與真實(shí)值之間的誤差，γ為正則化參數(shù)。得到LS-SVM預(yù)測函數(shù)為

(13)

式中，核函數(shù)K(x,xi)仍使用徑向高斯核：

為了檢驗(yàn)?zāi)Ｐ偷念A(yù)測性能，用最優(yōu)參數(shù)構(gòu)成的LS-SVM模型對測試樣本進(jìn)行訓(xùn)練，得到預(yù)測輸出和均方誤差的值。比較訓(xùn)練集與測試集的均方誤差，若兩者接近則說明該模型具有較好的預(yù)測性能；若相差較大，則要對訓(xùn)練樣本重新進(jìn)行訓(xùn)練，直到找到最優(yōu)的參數(shù)[13]。

2)實(shí)現(xiàn)步驟

① 訓(xùn)練LS-SVM

給定時(shí)間序列{xt}(t=1,2,…,N)，采用前n(n

(14)

式中：K(xi,xt)為核函數(shù)；αi∈R為拉格朗日因子；b?R為偏置項(xiàng)。

② 分布參數(shù)預(yù)測

由于xn-m+1={xn-m+1,xn-m+2, …,xn}，則第1步預(yù)測模型可表示為

(15)

將其加入時(shí)間序列得第二次預(yù)測模型。以此類推，第k步LS-SVM預(yù)測模型表示為

(16)

3 制導(dǎo)彈藥儲(chǔ)存性能評估

對制導(dǎo)彈藥進(jìn)行儲(chǔ)存性能評估，首先要在一定置信度下對檢測維護(hù)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析；其次利用LS-SVM預(yù)測未來某一時(shí)刻的分布函數(shù)，再計(jì)算未來時(shí)序的庫存制導(dǎo)彈藥不發(fā)生失效的概率；最后將所得值P(t)與評估閾值相比較，評估該批制導(dǎo)彈藥的庫存性能狀況并預(yù)測該批制導(dǎo)彈藥的壽命。制導(dǎo)彈藥儲(chǔ)存性能評估流程如圖1所示。

圖1 制導(dǎo)彈藥儲(chǔ)存性能評估流程框圖

4 制導(dǎo)彈藥儲(chǔ)存性能評估應(yīng)用

以假設(shè)為例，原始參數(shù)為90個(gè)，采集6年時(shí)間的150組樣本數(shù)據(jù)。

4.1 檢測維護(hù)數(shù)據(jù)的特征提取

1)粗篩選。90種參數(shù)應(yīng)用到模型中會(huì)產(chǎn)生巨大的“尾巴”導(dǎo)致評估結(jié)果的不精確，所以，先采用方差剔除法去除一部分對彈藥性能影響不大的參數(shù)，或者說參數(shù)變動(dòng)很小，基本處于穩(wěn)定狀態(tài)。計(jì)算結(jié)果見表1。

表1 取樣數(shù)據(jù)波動(dòng)情況(部分)

將方差值大于10的參數(shù)認(rèn)定為波動(dòng)參數(shù)，選取為輸入?yún)?shù)，使用KPCA進(jìn)行降維處理。計(jì)算其在各維投影，結(jié)果如表2所示。

表2 貢獻(xiàn)率大于90%特征向量

4.2 儲(chǔ)存性能評估及壽命預(yù)測

根據(jù)提取出的2個(gè)主成分KPCA1、KPCA2對制導(dǎo)彈藥儲(chǔ)存性能進(jìn)行評估。

1)分布類型假設(shè)。假設(shè)某批次制導(dǎo)彈藥的綜合參數(shù)KPCA1、KPCA2均符合正態(tài)分布。對其檢驗(yàn)，結(jié)果顯示均為：H=0，P=0.500，通過檢驗(yàn)，數(shù)據(jù)服從正態(tài)分布。

2)參數(shù)估計(jì)。按檢測的時(shí)間序列將數(shù)據(jù)分組，調(diào)用nomfit函數(shù)，計(jì)算出每年綜合數(shù)據(jù)分布的正態(tài)參數(shù)，對于部分彈藥的缺失數(shù)據(jù)采用擬合法進(jìn)行補(bǔ)充，得到正態(tài)參數(shù)的估計(jì)值見表3。

表3 正態(tài)參數(shù)估計(jì)值

3)分布擬合檢驗(yàn)。利用表3中參數(shù)計(jì)算2008—2013年的檢測維護(hù)數(shù)據(jù)分布概率密度函數(shù)g(x,βi)，并進(jìn)行擬合檢驗(yàn)。分析可知，Hi=0，Pi>0.05，CVi>KSSTATi。即符合假設(shè)接受條件，擬合效果較好。見表4。

4)LS-SVM預(yù)測

① 參數(shù)尋優(yōu)

將每次尋優(yōu)確定的正規(guī)化參數(shù)γ和高斯核寬δ代入式(13)行第一輪預(yù)測。將預(yù)測值與時(shí)序中2013年數(shù)據(jù)作對比，當(dāng)相對誤差小于0.005時(shí)尋優(yōu)結(jié)束，確定最優(yōu)正規(guī)化參數(shù)γ和高斯核寬δ，尋優(yōu)結(jié)果如表5所示。

表4 K-S分布擬合結(jié)果

表5 尋優(yōu)結(jié)果

② 訓(xùn)練LS-SVM

將最優(yōu)參數(shù)和訓(xùn)練樣本輸入LS-SVM進(jìn)行學(xué)習(xí)訓(xùn)練，計(jì)算拉格朗日因子向量和偏置項(xiàng)b，計(jì)算結(jié)果如表6所示。

表6 預(yù)測模型參數(shù)結(jié)果

③ 參數(shù)預(yù)測

已知參數(shù)序列作為樣本，代入式(4)第一步預(yù)測。之后將預(yù)測值加入訓(xùn)練樣本，代入式(16)作第二步預(yù)測，依次類推，預(yù)測未來10年的參數(shù)序列(表7)。

表7 分布參數(shù)預(yù)測時(shí)序

5)壽命預(yù)測及儲(chǔ)存性能評估

① 閾值選取

分別確定KPCA1、KPCA2的失效閾值。結(jié)果見表8。

表8 失效閾值

② 預(yù)測結(jié)果

將表8、表9數(shù)據(jù)代入式(8)得到2008—2023年分布密度函數(shù)序列。再將函數(shù)時(shí)序和閾值代入式(5)分別計(jì)算出兩個(gè)主成分的Pi(t)時(shí)序。從而得出聯(lián)合概率P(t)的時(shí)序。結(jié)果如圖2所示。

圖2 彈藥儲(chǔ)存性能評估

若按照P(t)值低于0.7界定彈藥壽命，該批次出廠時(shí)間為2005年，可初步認(rèn)為該批彈藥2019年左右到壽命。

③ 評估結(jié)果

將儲(chǔ)存性能P(t)分為4個(gè)區(qū)間，以界定彈藥儲(chǔ)存性能等級，如表9所示。按照評估標(biāo)準(zhǔn)，劃分出某批次彈藥性能實(shí)時(shí)狀態(tài)。由圖2可知，截止2017年，該批次彈藥正處在健康和良好的過渡期，技術(shù)狀態(tài)應(yīng)處于優(yōu)良狀態(tài)。

表9 彈藥儲(chǔ)存性能等級

④ 經(jīng)濟(jì)效益

假設(shè)該批次彈藥的定檢周期是每年檢測維護(hù)一次，費(fèi)用為萬元。按照計(jì)劃，在其到壽前，一共需要檢測維護(hù)費(fèi)用15萬元。若按照表10優(yōu)化檢測維護(hù)周期，則共需檢測維護(hù)費(fèi)用13.5萬元，共節(jié)省了10%的檢測維護(hù)經(jīng)費(fèi)(表10)。

表10 優(yōu)化后檢測周期費(fèi)用

5 結(jié)論

1)彈藥在實(shí)際儲(chǔ)存過程中，因儲(chǔ)存環(huán)境、測量誤差等原因與規(guī)定壽命稍有出入，可通過改善庫房儲(chǔ)存條件、增加樣本量進(jìn)一步提高精度。

2)當(dāng)預(yù)測P(t)值較高時(shí)，可大量減少檢測次數(shù)；當(dāng)預(yù)測P(t)值到中期時(shí)，通過增加檢測次數(shù)來評估制導(dǎo)彈藥性能；年末進(jìn)入危險(xiǎn)狀態(tài)，要縮短檢測周期，增大維護(hù)力度，作好延壽或報(bào)廢準(zhǔn)備。

3)根據(jù)制導(dǎo)彈藥的性能評估進(jìn)行合理任務(wù)安排和維護(hù)的方法，有利于克服過檢和漏檢現(xiàn)象，減少任務(wù)工作量，節(jié)約大量的人力物力，可廣泛推廣應(yīng)用。