于 航,黃 奕,李豪杰,馬少杰,魏亞偉

(1.南京理工大學(xué) 機械工程學(xué)院,江蘇 南京 210094;2.陸軍裝備科研與訂購部駐南京地區(qū)第四軍代室,江蘇 南京 210094)

定距空炸引信可控制彈丸在預(yù)定的目標(biāo)區(qū)起爆,攻擊具有高機動性與高隱蔽性目標(biāo)[1]。具有定距空炸功能的電子時間引信,其定距精度是影響定距空炸彈藥作用效果的主要因素。而彈丸炮口初速跳動導(dǎo)致定距空炸引信的定距精度存在很大波動[2],因此,為了解決此問題提出了一種具有自測速功能的電子時間引信,該引信對彈丸實際炮口初速進行測量并依據(jù)實測彈丸炮口初速修正引信實際作用時間,以實現(xiàn)精確起爆。因此,準(zhǔn)確、快速地進行炮口初速測量,對保證自測速電子時間引信作用精度至關(guān)重要。

現(xiàn)有引信自測速方法主要分為兩類:一種是磁鐵-線圈法,適用于微旋或不旋彈藥;另一種是適用于旋轉(zhuǎn)彈丸的地磁測速法、章動法、光學(xué)法等[3]。對于利用地磁測速法的旋轉(zhuǎn)彈丸,為了提高引信的測速精度,前人做了多方面工作。文獻[4]研究了火炮膛線纏角有9′的誤差時,引信測速值與雷達測速值相比存在1.75%左右的偏差。文獻[5]分析認為,引信后坐開關(guān)閉合閾值的選擇會影響測速起點,進而影響測速精度,且針對某型彈給出了閉合閾值選取范圍。文獻[6]為了同步武器平臺與引信的時鐘,將雙向時間戳交換方法應(yīng)用于引信時鐘同步,以提高引信定時精度。對于基于地磁測速法的彈丸,其自測速度數(shù)據(jù)來自引信電路所測得的彈丸轉(zhuǎn)數(shù)信號,而引信軟件對轉(zhuǎn)數(shù)數(shù)據(jù)的提取與篩選對初速精度計算非常關(guān)鍵,然而現(xiàn)有文獻中未發(fā)現(xiàn)對轉(zhuǎn)數(shù)數(shù)據(jù)提取策略的研究。

彈丸在出炮口計轉(zhuǎn)數(shù)時容易受到炮口火焰等離子體、磁場突變、低頻小信號、彈丸章動等外界因素影響。在這些外部干擾存在的情況下,引信控制器捕捉到的轉(zhuǎn)數(shù)數(shù)據(jù)中難免會存在一定量的干擾數(shù)據(jù)。為了在數(shù)據(jù)提取時去除干擾數(shù)據(jù),通常設(shè)置固定的數(shù)據(jù)篩選區(qū)間進行轉(zhuǎn)數(shù)數(shù)據(jù)剔除,而為了保證測速區(qū)間的普適性,通常選擇較大范圍的測速區(qū)間。但常規(guī)測速方法中數(shù)據(jù)提取多按時序順序進行提取及計算,即使存在數(shù)據(jù)篩選區(qū)間仍會提取到離真實轉(zhuǎn)數(shù)數(shù)據(jù)偏差較大的干擾數(shù)據(jù)。針對該問題本文提出了一種基于中位數(shù)法的引信轉(zhuǎn)數(shù)數(shù)據(jù)提取策略,當(dāng)測速數(shù)據(jù)數(shù)組中存在干擾數(shù)據(jù)問題時,仍可最大程度精確測速。

1 計轉(zhuǎn)數(shù)測速原理及引信測速模型

1.1 地磁傳感器的工作原理

對于地磁測速法,地磁傳感器是利用被測物體在空間中的不同運動狀態(tài),感知地磁場變換而產(chǎn)生與運動狀態(tài)相關(guān)的感應(yīng)電動勢。引信中的地磁傳感器是用一定直徑的漆包線繞制在特定的線圈骨架上,當(dāng)彈丸在空中旋轉(zhuǎn)飛行時,線圈切割地磁磁感線,產(chǎn)生一個周期性變換的感應(yīng)電動勢。引信控制器識別并捕捉該感應(yīng)電動勢,即可得到彈丸轉(zhuǎn)數(shù)[3,7]。

1.2 計轉(zhuǎn)數(shù)測速模型

對于線膛火炮發(fā)射的旋轉(zhuǎn)彈丸,彈丸出炮口后,在一定射擊范圍內(nèi),其每旋轉(zhuǎn)一圈彈丸沿發(fā)射方向前進一個火炮膛線導(dǎo)程,且該特性與彈丸出炮口速度無關(guān)[8]。導(dǎo)程與火炮身管口徑的關(guān)系為

l=ηD

(1)

式中:l為火炮的導(dǎo)程,η為身管膛線纏度,D為火炮身管直徑。對于等齊膛線火炮身管,其導(dǎo)程與纏角之間的關(guān)系為

(2)

式中:α為身管纏角。聯(lián)立式(1)、式(2),可得到tanα=π/η。

(3)

圖1為理想狀態(tài)下引信轉(zhuǎn)數(shù)波形,根據(jù)引信實際測速過程,引信控制器分別捕捉轉(zhuǎn)數(shù)波形上升沿、下降沿與閾值電壓范圍[U1,U2]交點對應(yīng)的時刻,其中設(shè)定閾值電壓的作用是剔除幅值過低的干擾波形。將相鄰2個上升沿、下降沿時刻點差值作為彈丸旋轉(zhuǎn)周期,對該周期進行提取與篩選,最終用來計算彈丸炮口初速。

例如,在圖1捕捉到的上升沿時刻t1,t3,t5,t7,…,下降沿時刻t2,t4,t6,t8,…,定義ΔTi=ti+2-ti,(1≤i≤s,i∈Z*),其中,s為捕捉到的時刻點數(shù),Z*表示正整數(shù)集。ΔTi為所測得的彈丸旋轉(zhuǎn)信號對應(yīng)的計時周期值。在上述分析的基礎(chǔ)上,可計算獲得彈丸炮口速度。

圖1 彈丸理想轉(zhuǎn)數(shù)波形

2 彈丸自測速偏差分析

在引信測速過程中,引信控制器接收來自地磁傳感器并經(jīng)過信號處理電路處理后的轉(zhuǎn)數(shù)波形電壓信號,采用軟件捕捉該電壓信號與設(shè)定閾值電壓交點對應(yīng)的時刻并加以存儲。待測速結(jié)束后,讀取部分存儲的數(shù)據(jù)用于初速計算。

彈丸自測速時偏差來源有很多,主要可分為兩類誤差:一類是測量系統(tǒng)誤差,包含計時起點誤差、測速原理誤差、計時器基準(zhǔn)誤差及計時器漂移誤差等[5-6];另一類為被測狀態(tài)誤差,包含火炮膛線纏角誤差[4],彈丸章動引起的誤差等。在已有文獻中除測速原理誤差外,上述文獻均有分析,因此本文對測速原理中的數(shù)據(jù)提取誤差進行分析。

根據(jù)式(3),可得彈丸旋轉(zhuǎn)周期的處理誤差對初速的影響,即

(4)

可見測速精度與彈丸旋轉(zhuǎn)周期值呈負相關(guān)關(guān)系。當(dāng)實測周期值較標(biāo)準(zhǔn)值周期值小時,測得的初速比實際初速偏大;當(dāng)實測周期值較標(biāo)準(zhǔn)值周期值大時,測得的初速比實際初速偏小。

現(xiàn)有轉(zhuǎn)數(shù)數(shù)據(jù)提取算法如圖2所示。圖中,D為火炮口徑;η為火炮纏度;v1,v2分別為允許引信開始測速的速度上限、下限;Tall為符合篩選區(qū)間內(nèi)的周期值總和;j為符合測速區(qū)間周期值個數(shù);m為設(shè)定的總周期提取個數(shù);T′all為平均周期值;v為獲得的測速值。將處理后的地磁傳感器信號送至引信控制器進行閾值判別并保存波形上升沿與下降沿對應(yīng)的時刻點,將相鄰上升沿、下降沿2個時刻點之間的差值作為彈丸旋轉(zhuǎn)周期值。但地磁傳感器及引信電路受外界與自身狀態(tài)影響,不可避免地會產(chǎn)生干擾信號,而這些干擾信號就混合在已經(jīng)保存的時刻點中,因此在對周期值進行提取時,要用篩選區(qū)間[Dη/v2,Dη/v1]進行篩選,將不滿足該周期區(qū)間的周期值進行剔除。將滿足篩選區(qū)間范圍的旋轉(zhuǎn)周期值進行進一步保存,將不滿足者剔除。最終提取m個有效周期,并將計算平均值T′all作為彈丸真實旋轉(zhuǎn)周期用于初速計算。

圖2 引信測速流程圖

在上述分析過程中,轉(zhuǎn)數(shù)數(shù)據(jù)的提取按照其保存順序進行,而干擾數(shù)據(jù)的周期值也可能滿足周期篩選區(qū)間。因此存在干擾數(shù)據(jù)也被選入有效周期序列中求均值的可能。然而在求平均值時,受到極端干擾數(shù)據(jù)影響,所求的轉(zhuǎn)數(shù)平均值會偏離整體數(shù)據(jù)的特性,使按順序篩選后直接求均值的方法存在初速計算不準(zhǔn)確的情況。

3 提高測速精度的引信轉(zhuǎn)數(shù)數(shù)據(jù)處理算法

由于干擾數(shù)據(jù)數(shù)值大小與正常轉(zhuǎn)數(shù)數(shù)據(jù)不一致,考慮到中位數(shù)的優(yōu)點是不會受到偏大或偏小的數(shù)據(jù)的影響,具有良好的主體數(shù)據(jù)特性[9],因此在轉(zhuǎn)數(shù)數(shù)據(jù)提取過程中,選擇提取中位數(shù)周圍的數(shù)據(jù),可對干擾數(shù)據(jù)進行剔除,獲得具有整體數(shù)據(jù)特性的轉(zhuǎn)數(shù)數(shù)據(jù)。

中位數(shù)求取過程的核心是數(shù)據(jù)的排序,而數(shù)據(jù)排序方法非常多,但要適用于引信,必須具有算法高效、快捷、易編程、占用空間小等特點。

傳統(tǒng)冒泡排序法具有算法簡單、穩(wěn)定、易編程等優(yōu)點,其操作機理是對任意兩相鄰的數(shù)據(jù)進行比較大小并進行排序,冒泡排序法的數(shù)據(jù)處理執(zhí)行次數(shù)與數(shù)組元素個數(shù)關(guān)系式為

(5)

式中:S1(a)為算法數(shù)據(jù)處理執(zhí)行次數(shù),a為數(shù)組中元素個數(shù),其時間復(fù)雜度為O(a2)[10]。

目前最為常用的排序方法是Tony Hoare于1962年提出的快速排序法。該方法采用分治思想,首先在數(shù)組中任意選擇一個基準(zhǔn)數(shù),然后對剩余數(shù)據(jù)進行分組,小于等于基準(zhǔn)數(shù)的在左側(cè),大于等于基準(zhǔn)數(shù)的在右側(cè)。這樣會得到2個新數(shù)組,再在新數(shù)組中選擇一個基準(zhǔn)數(shù),按以上過程進行遞歸,得到一個有序數(shù)列?？焖倥判蚍〝?shù)據(jù)處理執(zhí)行次數(shù)與數(shù)組元素個數(shù)關(guān)系式為

(6)

式中:k為計數(shù)常數(shù)。如式(6)所示,快速排序法的數(shù)據(jù)處理執(zhí)行次數(shù),其時間復(fù)雜度為O(alog2a)[10]。

由于引信為弱硬件,數(shù)據(jù)運算能力一般,在數(shù)據(jù)量一定的情況下,需要選擇對數(shù)據(jù)排序執(zhí)行次數(shù)較少的排序方法,以減少軟件執(zhí)行次數(shù),提高數(shù)據(jù)處理速度。因此,在相同數(shù)據(jù)量下對快速排序法與冒泡排序法進行排序,對軟件需執(zhí)行的次數(shù)進行比較,如圖3所示。

可見快速排序法的優(yōu)勢是,當(dāng)排序數(shù)據(jù)數(shù)量大于7時,軟件執(zhí)行次數(shù)比冒泡排序法少,并且數(shù)據(jù)越多時優(yōu)勢越明顯。對于引信而言,其控制器捕捉到的轉(zhuǎn)數(shù)數(shù)據(jù)量一般小于100,因此可選擇快速排序法作為引信轉(zhuǎn)數(shù)數(shù)據(jù)排序算法。

圖3 2種排序算法排序數(shù)據(jù)量與軟件執(zhí)行次數(shù)的關(guān)系

在以上分析的基礎(chǔ)上,本文提出了一種利用中位數(shù)法原理進行彈丸旋轉(zhuǎn)周期值篩選的方法。該方法是對引信控制器所存儲的彈丸旋轉(zhuǎn)周期值數(shù)組求取中位數(shù),并提取中位數(shù)左右兩側(cè)多組數(shù)據(jù)。數(shù)據(jù)提取總數(shù)為控制程序中預(yù)先設(shè)定的提取數(shù)量,并將這些數(shù)據(jù)求均值,再將得到的旋轉(zhuǎn)周期值進行篩選區(qū)間判斷,通過后加以運算即可得到彈丸炮口初速。

圖4 基于中位數(shù)法的引信測速流程圖

如圖4所示,引信控制器捕捉到彈丸轉(zhuǎn)數(shù)波形,分別提取到波形上升沿時刻t1,t3,t5,t7,…;下降沿時刻t2,t4,t6,t8,…。

根據(jù)ΔTi=ti+2-ti,(1≤i≤s,i∈Z*),計算相鄰兩時刻點差值,得到對應(yīng)的彈丸旋轉(zhuǎn)周期值數(shù)組;使用快速排序法對彈丸旋轉(zhuǎn)周期值數(shù)組進行排序,得到新的數(shù)組ΔT′i,并判斷出新數(shù)組中的中位數(shù)ΔT′m,在數(shù)據(jù)提取時,以中位數(shù)ΔT′m為中心分別向左右兩側(cè)提取數(shù)據(jù)(表示為ΔT′m±p,p為0及正整數(shù),并依次增加),并將數(shù)據(jù)進行累加,提取滿程序規(guī)定數(shù)量2q+1個數(shù)據(jù)(q為程序設(shè)定的常數(shù))。對提取的數(shù)據(jù)求得均值T′all,并判斷是否在測速區(qū)間內(nèi)。在滿足條件的情況下,若在區(qū)間內(nèi),則將T′all作為彈丸真實的旋轉(zhuǎn)周期值;若不在區(qū)間內(nèi),說明引信測速值無效,不進行修正。

4 中位數(shù)轉(zhuǎn)數(shù)處理方法有效性驗證

本試驗編制了如圖5所示的可視化模擬引信測速軟件,將前期在靶場回收到的轉(zhuǎn)數(shù)數(shù)據(jù)導(dǎo)入該測速軟件中,輸入幅值篩選區(qū)間與測速區(qū)間進行半實物仿真。該軟件采用2種計算模式分別對應(yīng)平均值算法與中位數(shù)算法,進行運行計算即可得到該組數(shù)據(jù)的彈丸初速。

圖5 模擬引信測速系統(tǒng)

利用模擬引信測速系統(tǒng)運算多組數(shù)據(jù)后,將平均值法與中位數(shù)法的測速結(jié)果匯總,如表1所示。將不同算法得到的炮口測速值v與雷達初速v0進行對比,表中,Δv為雷達初速與不同算法測速偏差的絕對值,δ為相對應(yīng)的不同測速方法的偏差率的絕對值。

由表1中數(shù)據(jù)可看出,相同的原始數(shù)據(jù),基于中位數(shù)法轉(zhuǎn)數(shù)數(shù)據(jù)提取方法明顯較基于平均值法得到的炮口初速值更接近雷達測速值。中位數(shù)法其平均測速偏差率為0.21%,而平均值法平均偏差率為0.975%,因此在測速精度方面,中位數(shù)法要優(yōu)于平均值法。結(jié)合本文第3部分對排序方法的分析,與傳統(tǒng)平均值法相比,選擇快速排序法并結(jié)合中位數(shù)法對引信轉(zhuǎn)數(shù)數(shù)據(jù)進行提取,在相同原始數(shù)據(jù)條件下,獲得的測速值更接近雷達測速值。因此,本文提出的中位數(shù)法轉(zhuǎn)數(shù)數(shù)據(jù)處理方法可提高彈丸測速精度,實現(xiàn)可編程電子時間引信精確打擊目標(biāo)的目的。

表1 試驗數(shù)據(jù)

5 結(jié)論

為了減小可編程電子時間引信測速誤差,提高炸點精度,本文提出了可編程時間引信計轉(zhuǎn)數(shù)測速中位數(shù)轉(zhuǎn)數(shù)數(shù)據(jù)處理方法,該方法考慮到當(dāng)引信轉(zhuǎn)數(shù)數(shù)據(jù)中存在部分干擾數(shù)據(jù)時,現(xiàn)有方法不能可靠剔除干擾的情況,鑒于干擾數(shù)據(jù)對中位數(shù)無影響的特性,結(jié)合快速排序法,使用中位數(shù)法提取數(shù)據(jù)。采用模擬引信測速系統(tǒng)進行驗證,與雷達測速數(shù)據(jù)相比,中位數(shù)法的平均測速偏差率為0.21%,而平均值法平均偏差率為0.975%,因此,采用中位數(shù)法測速精度要優(yōu)于平均值法。半實物仿真試驗驗證結(jié)果表明,該方法提高了引信的測速精度,進而提高了彈丸精確打擊能力。