鄭宗貴 王繼平 林紅斌 解 靜

?

一種基于雷達(dá)導(dǎo)引頭的彈道導(dǎo)彈爆高解算方法

鄭宗貴①王繼平*①②林紅斌①解 靜①

①(第二炮兵裝備研究院 北京 100085)②(96271部隊(duì) 洛陽 471000)

目前彈道導(dǎo)彈主要通過無線電引信測高來進(jìn)行戰(zhàn)斗部爆高控制，但彈頭的再入環(huán)境對無線電引信測高裝置的要求非常高，且引信存在著易被干擾的問題。為此，該文針對采用雷達(dá)景象匹配末制導(dǎo)的彈道導(dǎo)彈，提出一種基于雷達(dá)導(dǎo)引頭的彈道導(dǎo)彈爆高控制方法，通過雷達(dá)導(dǎo)引頭測高數(shù)據(jù)與慣性導(dǎo)航數(shù)據(jù)的融合處理，實(shí)時(shí)計(jì)算彈頭飛行高程，來實(shí)現(xiàn)爆高控制，既保證了精度，又增強(qiáng)了抗干擾能力。仿真算例證明了該方法的有效性和爆高解算精度。

雷達(dá)導(dǎo)引頭；彈道導(dǎo)彈；爆高

1 引言

近些年出現(xiàn)了一些新的測高設(shè)備與方法。點(diǎn)頭式測高雷達(dá)容易受到欺騙[10]；機(jī)載預(yù)警雷達(dá)測高精度易受地-海雜波、大氣折射、電波傳播多徑效應(yīng)等因素影響[11]；InSAR系統(tǒng)高程測量受信噪比、體散射、時(shí)間去相干等非理想因素的影響[12]；杜曉輝等人[13,14]提出利用已發(fā)展完備的地面移動通信基站作為氣壓差分測量基準(zhǔn)點(diǎn)的差分氣壓測高方法，雖然精度高，但移動通信基站易受干擾，這些測高設(shè)備或方法均不適合在彈道導(dǎo)彈上使用。針對采用景象匹配末制導(dǎo)的彈道導(dǎo)彈，溫求遒等人[15]提出了利用圖像導(dǎo)引頭角量測信息進(jìn)行彈目相對距離估計(jì)，還未用于爆高的控制；蘇昂[8]提出了單幀多點(diǎn)匹配和幀間單點(diǎn)匹配測高方法，但是該方法限制在景象匹配區(qū)測高，不適用于爆高的控制。為此，本文提出了一種基于雷達(dá)導(dǎo)引頭的彈道導(dǎo)彈爆高解算方法，利用雷達(dá)導(dǎo)引頭測高數(shù)據(jù)，與慣性導(dǎo)航數(shù)據(jù)進(jìn)行融合處理，實(shí)時(shí)計(jì)算出滿足戰(zhàn)斗部爆高精度要求的高程，既能增強(qiáng)抗干擾能力[16]，又能保證精度。

由于彈道導(dǎo)彈慣導(dǎo)工具誤差的存在，其導(dǎo)航誤差會隨著工作時(shí)間逐漸增大，由導(dǎo)航定位數(shù)據(jù)確定導(dǎo)彈飛行中的高程，很難滿足爆高精度要求。然而對于采用雷達(dá)景象匹配末制導(dǎo)導(dǎo)彈，其給出的雷達(dá)某一時(shí)間間隔的連續(xù)測高數(shù)據(jù)，精度非常高，可將該數(shù)據(jù)和慣性導(dǎo)航高程解算數(shù)據(jù)進(jìn)行融合處理，實(shí)時(shí)計(jì)算高程以滿足爆高控制精度要求。

雷達(dá)景象匹配導(dǎo)引頭測高時(shí)，所得的測高數(shù)據(jù)是雷達(dá)導(dǎo)引頭天線相對地面的高程，而爆高數(shù)據(jù)要求的是彈頭質(zhì)心相對目標(biāo)點(diǎn)的高程，因而，需要將雷達(dá)導(dǎo)引頭的測高數(shù)據(jù)由相對天線位置轉(zhuǎn)換到相對彈頭質(zhì)心位置。另外，由于慣性導(dǎo)航數(shù)據(jù)是相對于慣組的安裝位置，將雷達(dá)數(shù)據(jù)與慣性導(dǎo)航高程解算數(shù)據(jù)融合時(shí)，需將慣性導(dǎo)航數(shù)據(jù)轉(zhuǎn)換到相對于彈頭的質(zhì)心位置。

綜上所述，基于雷達(dá)導(dǎo)引頭的彈道導(dǎo)彈爆高解算思路是首先將慣性導(dǎo)航數(shù)據(jù)轉(zhuǎn)換到相對于導(dǎo)彈彈頭的質(zhì)心位置，并解算對地高程；接著將景象匹配雷達(dá)最后一次成像后的連續(xù)測高數(shù)據(jù)轉(zhuǎn)換到導(dǎo)彈質(zhì)心位置；然后將兩者之差(即慣導(dǎo)解算高程誤差值)擬合成時(shí)間的函數(shù)，根據(jù)導(dǎo)彈飛行時(shí)間，修正慣導(dǎo)解算高程誤差，獲得高精度的高程，滿足爆高控制精度要求。

2 基于雷達(dá)導(dǎo)引頭的爆高解算方法

首先，給出慣組安裝位置偏離彈頭質(zhì)心引起的速度位置偏差模型和雷達(dá)導(dǎo)引頭天線位置偏離彈頭質(zhì)心引起的高程偏差模型；在此基礎(chǔ)上給出兩者的數(shù)據(jù)融合模型及爆高解算方法。

2.1慣測組合安裝位置偏離彈頭質(zhì)心引起的速度位置偏差模型

一般情況下，彈載慣性測量裝置安裝在導(dǎo)彈儀器艙內(nèi)，安裝位置與彈頭質(zhì)心不重合，因此再入段加速度表測量的視加速度為相對安裝位置處的值，與相對質(zhì)心處的值存在差異。然而，爆高控制是按彈頭質(zhì)心運(yùn)動的高程進(jìn)行控制，因而需將加速度表測量的視加速度歸算到質(zhì)心，再計(jì)算彈頭的質(zhì)心運(yùn)動。

2.1.1基本假設(shè) 為便于將加速度表測量的視加速度歸算到質(zhì)心，作如下基本假設(shè)：

(1)忽略彈頭飛行時(shí)的動態(tài)過程，即滿足瞬時(shí)平衡假設(shè)；

(2)在標(biāo)準(zhǔn)條件下飛行，不計(jì)各種干擾對測量裝置輸出信號的影響；

(3)忽略彈頭的彈性振動。

圖1 慣組測量裝置安裝偏差矢量圖

令

根據(jù)彈體坐標(biāo)系與慣性坐標(biāo)系間的關(guān)系矩陣，可得慣測裝置偏離質(zhì)心安裝時(shí)的視加速度偏差在慣性坐標(biāo)系各軸上的分量式為

彈頭質(zhì)心視加速度為

用式(8)修正慣測視速度增量，并代入彈頭導(dǎo)航計(jì)算式(9)，可計(jì)算彈頭質(zhì)心的速度和位置。

式中

2.2 雷達(dá)天線偏離彈頭質(zhì)心引起的高程偏差模型

圖2 目標(biāo)坐標(biāo)系

高程偏差為

2.3 高精度爆高解算方法

雷達(dá)測高方程為

若目標(biāo)區(qū)平坦，雷達(dá)天線對地高可近似為對目標(biāo)點(diǎn)高；若已知目標(biāo)區(qū)地形高程數(shù)據(jù)，可根據(jù)慣導(dǎo)數(shù)據(jù)計(jì)算彈頭經(jīng)緯坐標(biāo)，插值獲取當(dāng)前地形高程值，將雷達(dá)對地高歸算到雷達(dá)對目標(biāo)點(diǎn)高。則雷達(dá)測高與慣性導(dǎo)航計(jì)算高之差為

令

3 引爆時(shí)間預(yù)報(bào)

4 仿真算例

為了驗(yàn)證本文爆高解算方法的有效性，分別在地形平坦區(qū)和地形起伏區(qū)，采用最小二乘二階濾波模型，進(jìn)行爆高解算。

4.1 仿真條件

對標(biāo)準(zhǔn)彈道添加不同的速度誤差和位置誤差，形成500條慣組測量誤差彈道，進(jìn)行爆高解算精度仿真計(jì)算。

4.2 仿真結(jié)果

針對圖3的目標(biāo)區(qū)平坦地形情況，最小二乘濾波解算爆高仿真結(jié)果如圖4，圖5所示。

由仿真結(jié)果可以看出：濾波前后，測高誤差均值差異不大，都約為0 m；濾波后測高誤差標(biāo)準(zhǔn)差約為1.2 m，而單點(diǎn)測高的標(biāo)準(zhǔn)差約為2.2 m，可見濾波后，測高精度提高。

針對圖6的目標(biāo)區(qū)地形起伏區(qū)域情況，最小二乘濾波解算爆高仿真結(jié)果如圖7，圖8所示。

由仿真結(jié)果可看出：濾波前后，測高誤差均值差異不大，都約為0 m；濾波后測高誤差標(biāo)準(zhǔn)差約為1.4 m，而單點(diǎn)測高的標(biāo)準(zhǔn)差約為2.1 m，可見濾波后，測高精度提高。

5 結(jié)束語

本文針對采用雷達(dá)景象匹配末制導(dǎo)的彈道導(dǎo)彈，提出了基于雷達(dá)導(dǎo)引頭的彈道導(dǎo)彈爆高解算方法，由于該方法進(jìn)行了信息融合處理，而無線電裝置定高技術(shù)直接利用測高結(jié)果，相比之下，精度提高了2~3倍，但由于基于雷達(dá)導(dǎo)引頭的彈道導(dǎo)彈爆高解算方法采用了最小二乘信息融合計(jì)算，運(yùn)算量相比無線電裝置定高技術(shù)較大，但彈載計(jì)算機(jī)性能在不斷提高，目前計(jì)算能力已能適應(yīng)最小二乘信息融合計(jì)算的要求，工程實(shí)現(xiàn)已成為可能。

圖3 彈道、地面高程和天線波束照射區(qū)域

圖4 測高誤差均值

圖5 測高誤差標(biāo)準(zhǔn)差

圖6 彈道，地面高程和天線波束照射區(qū)域

圖7 測高誤差均值

圖8 測高誤差標(biāo)準(zhǔn)差

[1] 王磊, 許可, 徐曦煜, 等. 一種計(jì)算雷達(dá)高度計(jì)重跟蹤修正表的新方法及其應(yīng)用[J]. 電子與信息學(xué)報(bào), 2013, 35(4): 908-914.

Wang Lei, Xu Ke, Xu Xi-yu,.. A new method for computing radar altimeter look-up correction table and its application[J].&, 2013, 35(4): 908-914.

[2] 鮮勇. 基于攝動理論的爆高確定方法研究[J]. 飛行力學(xué), 2006, 24(4): 50-52.

Xian Yong. On the method of determinating explosion height based on perturbation theory[J]., 2006, 24(4): 50-52.

[3] 廖洪昌, 郭荊燕, 何念念, 等. 一種基于BP網(wǎng)絡(luò)的爆高延時(shí)方法[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2009, 29(4): 119-122.

Liao Hong-chang, Guo Jing-yan, He Nian-nian,.. A method of calculating the high-explosive delay based on back propagation[J].,,, 2009, 29(4): 119-122.

[4] 孟海東, 廖洪昌,陳奇昌, 等. 一種基于均勻設(shè)計(jì)的爆高延時(shí)計(jì)算方法[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2008, 28(6): 103-105.

Meng Hai-dong, Liao Hong-chang, Chen Qi-chang,.. An algorithm for burst height delay based on uniform design journal of projectiles[J].,, 2008, 28(6): 103-105.

[5] 嚴(yán)輝, 熊瑋, 南英, 等. 戰(zhàn)術(shù)地地導(dǎo)彈爆高延時(shí)方案計(jì)算與分析[J]. 導(dǎo)彈與航天運(yùn)載技術(shù), 1996, 24(2): 32-35.

Yan Hui, Xiong Wei, Nan Ying,.. Analysis and computation of timelag control for explosive height of tactical surface-to- surface missile[J]., 1996, 24(2): 32-35.

[6] 陳濤, 梁明. 某型無線電高度表模擬測高試驗(yàn)的結(jié)果分析[J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2011, 31(1): 118-121.

Chen Tao and Liang Ming. Data analysis of static simulating of altitude measuring using radio altimeter[J]., 2011, 31(1): 118-121.

[7] 高善清, 魏新亮. 氣壓測高誤差分析及修正[J]. 電子測量技術(shù), 2012, 35(9): 60-63.

Gao Shan-qing and Wei Xin-liang. Error analysis and correction for the atmospheric pressure measurement of altitude[J]., 2012, 35(9): 60-63.

[8] 蘇昂. 飛行器景象匹配測高測姿方法研究[D]. [碩士論文], 國防科技大學(xué), 2011.

Su Ang . Study on aerial vehicles height and pose estimation method based on scene matching[D]. [Master dissertation], National University of Defense Technology, 2011.

[9] 張麗榮, 馬利華, 王曉嵐. CAPS氣壓測高的仿真分析[J]. 微計(jì)算機(jī)信息, 2011, 27(3): 185-186.

Zhang Li-rong, Ma Li-hua, and Wang Xiao-lan. Simulation analysis of CAPS barometric altimetry[J]., 2011, 27(3): 185-186.

[10] 劉玉良, 肖文杰, 崔海濤. 點(diǎn)頭式測高雷達(dá)高度欺騙研究[J]. 艦船電子對抗, 2012, 35(4): 23-25.

Liu Yu-liang, Xiao Wen-jie, and Cui Hai-tao. Study of height deception for the nodding height-finding radar[J]., 2012, 35(4): 23-25.

[11] 劉波, 劉寶泉, 陳春暉. 機(jī)載預(yù)警雷達(dá)測高精度分析[J]. 雷達(dá)科學(xué)與技術(shù), 2012, 10(2): 133-137.

Liu Bo, Liu Bao-quan, and Chen Chun-hui. Altitude measurement precision analysis of airborne early warning radar[J]., 2012, 10(2): 133-137.

[12] 葛仕奇, 陳亮, 丁澤剛, 等. 相位誤差對InSAR相對測高精度影響[J]. 北京理工大學(xué)學(xué)報(bào), 2012, 32(2): 179-183.

Ge Shi-qi, Chen Liang, Ding Ze-gang,.. Effects of phase error on the relative height accuracy in interferometric synthetic aperture radar[J]., 2012, 32(2): 179-183.

[13] 杜曉輝, 裴軍, 張麗榮, 等. 基于地面移動通信基站的差分氣壓測高方法[J]. 北京航空航天大學(xué)學(xué)報(bào), 2013, 39(1): 82-88.

Du Xiao-hui, Pei Jun, Zhang Li-rong,..Differential barometric altimetry method based on mobile phone base stations[J]., 2012, 39(1): 82-88.

[14] 劉成, 張麗榮, 杜曉輝. 一種氣壓測高輔助下的地面移動網(wǎng)定位方法[J]. 桂林理工大學(xué)學(xué)報(bào), 2012, 32(1): 86-90.

Liu Cheng, Zhang Li-rong, and Du Xiao-hui.Terrestrial mobile network positioning method assisted by barometric altimetry technology[J]., 2012, 32(1): 86-90.

[15] 溫求遒, 夏群利, 蔡春濤. 基于圖像導(dǎo)引頭及慣導(dǎo)相對位置信息的彈目距離估計(jì)算法[J]. 北京理工大學(xué)學(xué)報(bào), 2012, 32(2): 141-145.

Wen Qiu-qiu, Xia Qun-li, and Cai Chun-tao. Method of range estimate based on relative position information of image seeker and INS[J]., 2012, 32(2): 141-145.

[16] 徐春夷. 國外導(dǎo)引頭技術(shù)現(xiàn)狀及發(fā)展趨勢[J]. 制導(dǎo)與引信, 2012, 33(2): 11-15.

Xu Chun-yi. The present state and development trend of forein seeker technology[J].&, 2012, 33(2): 11-15.

鄭宗貴： 男，1970年生，博士，高級工程師，研究方向?yàn)閷?dǎo)航、制導(dǎo)與控制.

王繼平： 男，1977年生，博士，講師，研究方向?yàn)閷?dǎo)航、制導(dǎo)與控制.

林紅斌： 男，1964年生，博士，研究員，研究方向?yàn)閷?dǎo)航、制導(dǎo)與控制.

解 靜： 女，1981年生，博士，助理研究員，研究方向?yàn)閷?dǎo)航、制導(dǎo)與控制.

A Calculation Method of Burst Height of Ballistic Missile Based on Radar Seeker

Zheng Zong-gui①Wang Ji-ping①②Lin Hong-bin①Xie Jing①

①(,100085,)②(96271,471000,)

At present the height of the warhead of ballistic missile is measured by the radio fuze for the burst height control of the warhead, but the reentry environment of the warhead has very high requirement for the measure set of the radio fuze, and there exists the issue of the radio fuze to be jammed. In the paper, for the ballistic missile adopting the radar scene matching terminal guidance, a calculation method of burst height of ballistic missile based on radar seeker is put forward, by fusing the height measure data of the scene matching radar seeker and the inertial navigation data, the height of the warhead is calculated in real time to realize the burst height control, which ensures the precision and improves anti-jamming capability. The validity of the method and the calculation precision of the burst height are proved by simulation and calculation.

Radar seeker; Ballistic missile; Burst height

TJ765.3+31

A

1009-5896(2014)02-0465-06

10.3724/SP.J.1146.2013.01142

王繼平 wjpzj1111@163.com

2013-07-30收到，2013-11-08改回