魏京超，豆俊奇，孫 策

(1 青島興道新特科技有限公司，山東 青島 266400；2 中國科學(xué)院西安光學(xué)精密機(jī)械研究所，陜西 西安 710119)

0 引言

末敏彈作為現(xiàn)代戰(zhàn)場的一種智能彈藥,它通過邊旋轉(zhuǎn)下落邊掃描探測的方式搜尋和打擊目標(biāo)[1]。前端敏感器利用目標(biāo)和背景的輻射差異探測目標(biāo)[2]。當(dāng)發(fā)現(xiàn)目標(biāo)時,引爆戰(zhàn)斗部朝向目標(biāo)發(fā)射[3]。因此,目標(biāo)識別和定位對于成功打擊至關(guān)重要。

單元探測掃描只能獲取一條掃描信號,末敏彈無法判定掃描線所掃到目標(biāo)的位置,無法判定打擊落點是否在中心區(qū)域[4]。多元紅外敏感器具有同時獲取多條掃描線的能力,可理解為對地面目標(biāo)“成二維像”。通過該紅外二維像計算目標(biāo)熱區(qū)尺寸和中心[5]。

和單元目標(biāo)識別只能利用信號響應(yīng)幅值判別相比[6],多元掃描探測可利用尺寸匹配方法[7]進(jìn)一步識別目標(biāo),提高排偽能力,同時又可以獲取戰(zhàn)斗部指向地面目標(biāo)的位置,判斷是否對準(zhǔn)中心區(qū)域。文獻(xiàn)[8]提出了一種多元紅外探測識別的檢測算法,未具體給出尺寸計算方法。在該文獻(xiàn)基礎(chǔ)之上,提出了一種熱區(qū)尺寸計算方法和定位方法,為多元紅外探測提供技術(shù)支撐。

1 多元紅外旋轉(zhuǎn)掃描原理

末敏子彈從拋射出開始邊下落邊旋轉(zhuǎn)進(jìn)入穩(wěn)態(tài)掃描階段[9],該階段原理如圖1所示。

圖1 末敏子彈多元探測穩(wěn)態(tài)掃描Fig.1 Steady-state scanning of multi-element detection of terminal sensitive projectile

末敏彈以掃描角θ、轉(zhuǎn)速ω和下落速度v對地面掃描,掃描間距為ΔR。初始掃描位置為A0(與X軸成角度α),經(jīng)時間t后到達(dá)掃描螺線上的A1點。該點坐標(biāo)表示為：

(1)

相鄰螺線的距離ΔR為：

(2)

式中n=ω/2π為末敏子彈繞鉛垂軸的角速度。多元探測器每個探測單元具有獨立的視場角,能夠探測到目標(biāo)的不同區(qū)域[10]。和單元掃描不同,多元探測器在地面上某個瞬時掃描的是一個區(qū)域,而不是一條線。

2 尺寸解算和定位方法

多元探測掃描到目標(biāo)時產(chǎn)生多路信號。末敏彈需要對該多路信號作尺寸解算和定位判定,才能完成后續(xù)的模版匹配和打擊決策。

2.1 尺寸解算

目標(biāo)熱區(qū)尺寸解算的基本思路是：在經(jīng)過對實時信號數(shù)據(jù)預(yù)處理后,通過對多路信號進(jìn)行邊緣點提取,獲取各個探測單元掃描目標(biāo)的脈沖寬度、探測到的單元數(shù)目和時序關(guān)系。根據(jù)各個探測單元的時序關(guān)系以及高度計的高度信息,計算目標(biāo)位置姿態(tài)角。然后結(jié)合姿態(tài)角,用上述多元的脈寬分布分別計算熱區(qū)的兩個邊長。

2.1.1 多路信號邊緣和脈寬提取

多元探測器掃描目標(biāo)時,掃描到目標(biāo)上的探測單元產(chǎn)生幅值響應(yīng),即波形信號[11]。掃入目標(biāo)和掃出目標(biāo)時產(chǎn)生幅值變化。對應(yīng)在信號曲線上則為邊緣點。響應(yīng)的探測信號包含掃入目標(biāo)時的左邊緣點和掃出目標(biāo)時的右邊緣點。掃描目標(biāo)上的時間即為脈寬。結(jié)合高度和掃描速度,可以計算出目標(biāo)在兩點間的橫向距離。多元探測器每個單元之間具有一定的視場角。掃描時,統(tǒng)計產(chǎn)生響應(yīng)的信號元數(shù)便可獲得縱向的目標(biāo)邊緣點之間的距離。

如圖2(a)所示,探測器第1路信號掃描線離開A′的時刻為T1,探測器第2路信號掃描線離開B′的時刻為T2,T2和T1的時間差為Δτ,由時間差Δτ與探測元掃描的線速度和高度信息可得到目標(biāo)所處視場中的角度,角度計算公式為：

圖2 目標(biāo)尺寸計算Fig.2 Target size calculation

L=H/cosθ×sinβ

(3)

α=arctan(Δτ×vc/L)

(4)

式中：L為兩個探測元視場間距離;H為末敏子彈離地高度;θ為紅外探測器傾斜角;β為探測元之間視場角;α為目標(biāo)所處紅外視場中的角度;vc為中心波元掃描過目標(biāo)的線速度。

2.1.2 掃描切向目標(biāo)的尺寸計算

如圖2(a)所示,掃描切向目標(biāo)的尺寸a可通過式(5)計算得出：

a=BE=BB′·cosα

(5)

其中：BB′=vr·tBB′,vr為末敏子彈旋轉(zhuǎn)速度,tBB′為第2路信號探測元掃描過目標(biāo)的時間;α為目標(biāo)所處紅外視場中的角度。

2.1.3 掃描徑向目標(biāo)的尺寸計算

如圖2(b)所示,掃描徑向目標(biāo)的尺寸b可通過式(6)計算得出：

(6)

其中,

(7)

2.2 目標(biāo)定位

多元掃描探測利用響應(yīng)信號的時空分布判斷目標(biāo)位置[12]。當(dāng)只有邊緣探測單元產(chǎn)生信號響應(yīng)時,并不起爆戰(zhàn)斗部,而是讓末敏彈繼續(xù)旋轉(zhuǎn)。當(dāng)多路信號時空呈現(xiàn)對稱分布時,說明目標(biāo)已進(jìn)入掃描的中心區(qū)域,這時起爆戰(zhàn)斗部,原理如圖3所示。

圖3 實時掃描示意圖Fig.3 Schematic diagram of real-time scanning

如圖3(a)所示,掃描區(qū)域未掃到目標(biāo)熱區(qū)時,多路信號沒有幅值響應(yīng)。如圖3(b)所示,多路信號邊緣探測元掃描到目標(biāo),紅點為掃到的目標(biāo)邊緣點。4路信號掃描目標(biāo)的時刻分別為(t4a,t4b),同理5路信號掃描對應(yīng)的時刻分別為(t5a,t5b)。如圖3(c)所示,目標(biāo)完全進(jìn)入到掃描區(qū)域時,多路信號產(chǎn)生的幅值影響。2、3、4路信號上升沿和下降沿對應(yīng)的時刻分別為(t2a,t2b),(t3a,t3b),(t4a,t4b)。

假設(shè)n為多元信號的總元數(shù),第i路信號產(chǎn)生的目標(biāo)響應(yīng)為Vi,信號產(chǎn)生的脈寬為di=v(tib-tia)。實時定位算法如圖4所示。

圖4 實時定位算法流程圖Fig.4 Flow chart of real-time positioning algorithm

多路信號實時掃描目標(biāo)時,首先進(jìn)入第一個判決條件,即對稱元信號響應(yīng)判別。假設(shè)n=5,i=4時,相對應(yīng)產(chǎn)生的目標(biāo)幅值響應(yīng)為V4,脈寬為d4=v(t4b-t4a)。當(dāng)|V4-V5+1-4|<ΔV時(ΔV為極小數(shù),可忽略為0),第4條信號與第2條信號目標(biāo)響應(yīng)呈對稱分布,此時目標(biāo)位于掃描中心;接下來進(jìn)入第二個判決條件,則存在max(di)=d(5+1)/2=d3,即當(dāng)?shù)?路信號與第2路信號呈對稱分布時,第3路信號存在最大脈寬。最后進(jìn)行目標(biāo)熱區(qū)尺寸判別。判別均通過后計算熱區(qū)中心,實現(xiàn)定位。

3 驗證與分析

通過實驗室內(nèi)旋平臺采集的數(shù)據(jù)驗證文中算法。旋轉(zhuǎn)平臺參數(shù)如表1所示,多元紅外探測器參數(shù)如表2所示。

表1 實驗室旋轉(zhuǎn)平臺參數(shù)Table 1 Parameters of rotating platform in laboratory

表2 多元紅外探測器參數(shù)Table 2 Parameters of multi-element infrared detector

3.1 尺寸解算算法驗證

下面從3個代表性的角度(0°、45°和90°)進(jìn)行驗證,波形信號如圖5所示。圖5(a)為目標(biāo)所處0°時的波形圖,3～9路單元探測到目標(biāo)信號;圖5(b)為目標(biāo)所處45°時的波形圖,2～9路單元探測到目標(biāo)信號;圖5(c)為目標(biāo)所處90°時的波形圖,4～8路單元探測到目標(biāo)信號。

圖5 不同角度波形圖Fig.5 Wave forms at different angles

針對圖5的3種情況進(jìn)行目標(biāo)尺寸解算,解算結(jié)果如表3所示。

表3 不同角度目標(biāo)尺寸解算結(jié)果Table 3 Calculation results of target size at different angles

由表3可知,根據(jù)各路信號提取的邊緣點進(jìn)行尺寸解算,求出各路平行于掃描方向的目標(biāo)距離L1和垂直掃描方向的目標(biāo)距離L2。利用不同路信號的時序關(guān)系解算出姿態(tài)角α。最后根據(jù)角度和高度信息解算出邊長a和b。

3.2 定位解算算法驗證

目標(biāo)以45°姿態(tài)角在紅外視場中不同縱向深度時的波形信號如圖6所示。圖6(a)為目標(biāo)剛進(jìn)入掃描區(qū)域時的波形圖,1和2路單元探測到目標(biāo)信號;圖6(b)為目標(biāo)一半進(jìn)入掃描區(qū)域時的波形圖,1～5路單元探測到目標(biāo)信號;圖6(c)為目標(biāo)完全進(jìn)入掃描區(qū)域時的波形圖,2～9路單元探測到目標(biāo)信號。

圖6 不同縱向位置波形圖Fig.6 Wave forms at different longitudinal positions

由表4可知,根據(jù)各路信號提取的邊緣點進(jìn)行響應(yīng)元數(shù)、是否中心對稱和最大脈寬判別。位置1只有1和2路單元探測到信號,說明目標(biāo)在掃描區(qū)域邊緣。位置2處有1～5路單元探測到目標(biāo)信號,說明目標(biāo)一部分進(jìn)入掃描區(qū)域。位置3處有2～9路單元探測到目標(biāo)信號,說明目標(biāo)完全進(jìn)入掃描區(qū)域。中心探測單元5和6脈寬最大,說明目標(biāo)在掃描區(qū)域中心。

表4 不同位置目標(biāo)尺寸解算結(jié)果Table 4 Calculation results of target size at different positions

4 結(jié)論

提出了一種多元紅外掃描探測目標(biāo)尺寸解算和定位方法,通過獲取多路信號的邊緣點和時序關(guān)系解算目標(biāo)尺寸,并通過響應(yīng)元數(shù)和最大脈寬對稱性定位目標(biāo)。對文中所提方法進(jìn)行了實驗驗證分析。實驗結(jié)果表明與單元紅外掃描探測相比,采用線陣紅外掃描探測按照所述方法可以獲取目標(biāo)熱區(qū)尺寸,實現(xiàn)定位,提高目標(biāo)識別能力和定位精度。