岳玫君，王啟智

(1．解放軍電子工程學(xué)院，合肥 230037；2．中國電子科技集團公司第三十八研究所，合肥 230088)

空間目標(biāo)檢測與實現(xiàn)方法研究

岳玫君1，王啟智2

(1．解放軍電子工程學(xué)院，合肥 230037；2．中國電子科技集團公司第三十八研究所，合肥 230088)

空間目標(biāo)由于其高速高機動等特點在觀察時間內(nèi)通?？缭蕉鄠€距離單元，使得為了目標(biāo)檢測進行能量相參積累變得很困難。通過對大時寬高速運動目標(biāo)雷達回波信號的建模，分析了拉伸效應(yīng)、脈沖間距離遷徙等因素對積累的影響。提出了一種基于欠采樣下Keystone變換的高速目標(biāo)檢測方法。該方法能有效補償距離走動，去除多普勒模糊，提高信噪比，改善雷達的檢測性能。通過仿真和實際數(shù)據(jù)分析，驗證了該算法的有效性，同時結(jié)合硬件平臺給出了工程實現(xiàn)方法。

空間目標(biāo)；Keystone變換；距離遷徙；工程實現(xiàn)

0 引 言

空間目標(biāo)具有探測距離遠、目標(biāo)RCS小、高速高機動性等特點[1]，對其檢測和跟蹤是當(dāng)今雷達技術(shù)領(lǐng)域的熱點與難點問題之一。對于此類信號的檢測通過增大天線孔徑、增強發(fā)射功率能達到很好的效果，但這些通常會大大增加系統(tǒng)成本，不符合綠色設(shè)計理念。通?？梢匝娱L積累時間以增加實際應(yīng)用的能量，達到降低系統(tǒng)信噪比的要求。對于高速的空間目標(biāo)進行長時積累，不可避免地會發(fā)生目標(biāo)包絡(luò)跨距離單元走動問題[2-4]。這不僅造成目標(biāo)的距離和速度信息失真，更使得目標(biāo)能量分散而達不到檢測門限，造成雷達作用距離縮短。Keystone變換[5]是在距離向脈沖壓縮的頻域?qū)β龝r間軸進行尺度變換，是回波距離校正的有效方法。

1 高速目標(biāo)回波分析

在窄帶情況下，由于分辨率的原因，目標(biāo)都被認為是一個散射點，同時由于在一定的時間和空間內(nèi)只是對某一個目標(biāo)感興趣，因此這里只對高速運動的單點目標(biāo)回波進行分析。

假設(shè)雷達發(fā)射的基帶信號為線性調(diào)頻脈沖信號：

式中，T為信號時寬，B為帶寬。

雷達在觀測衛(wèi)星、導(dǎo)彈之類的空間高速目標(biāo)時，一般發(fā)射時寬較寬、目標(biāo)速度較快。因此，回波信號存在拉伸效應(yīng)，不能只看作發(fā)射信號的遲延，僅考慮多普勒影響，其時間尺度因子也必須考慮。但是，在|v|/c?1(v是目標(biāo)的徑向速度，c是光速)的條件下，回波信號的基帶信號可以寫為[6]

式中，n為發(fā)射的脈沖個數(shù)，fc為載波頻率，A0是散射點回波的幅度，τn為第n個脈沖的延遲時間，R0為0時刻的距離，Tr為脈沖重復(fù)時間。

為了減小運算量，一般利用快速傅里葉變換(FFT)替代時域卷積來進行脈沖壓縮，脈壓后的信號頻譜為

對應(yīng)的時域信號為

由上式可以看出，脈沖壓縮時域信號的包絡(luò)是一個sinc函數(shù)。由sinc函數(shù)的性質(zhì)知，信號峰值的位置位于τn-Tfd/B。由于目標(biāo)的高速運動，目標(biāo)的包絡(luò)會發(fā)生與多普勒頻率成正比的偏移，其偏移的距離可以表示為cTfd/(2B)，若忽略相參積累時間里fd的變化，則這個偏移量對于每個脈沖回波相同，可以不予考慮。但是，脈壓以后目標(biāo)的包絡(luò)τn分量隨著每個脈沖數(shù)變化。這對于信號的能量積累不利，從而不利于對微弱目標(biāo)的檢測。后續(xù)處理若直接用FFT變換進行多普勒濾波，信號能量積累會有損失，不利于信噪比的提高；而且由于距離走動的影響，目標(biāo)的主瓣明顯展寬，其距離分辨力和多普勒分辨力也會有所降低。

2 欠采樣下Keystone變換及計算仿真

2.1 欠采樣下Keystone變換

所謂Keystone變換就是進行一種變量代換，將f-τn平面的矩形支撐域在f-τm平面變成倒梯形。它是對τn軸的伸縮變換，伸縮幅度與頻率有關(guān)。

當(dāng)fd?fc時，經(jīng)過Keystone變換后，可以寫為[7]

由上式可以看出，經(jīng)Keystone變換后脈沖壓縮信號峰值位置只與初始時刻有關(guān)，存在一個固定的速度距離耦合偏差，與雷達各脈沖周期回波無關(guān)，脈沖間的距離走動得以消除。

當(dāng)目標(biāo)的速度較快或者雷達的脈沖重復(fù)頻率較低時，常常會導(dǎo)致欠采樣，這時多普勒頻率是模糊的，可表示為

fd=fd0+nkPRF

式中，fd0是折疊以后的多普勒頻率，nk是折疊因子，PRF是脈沖重復(fù)頻率。

如果出現(xiàn)多普勒頻率模糊，就需要對Y(f,m)進行相位補償：

在處理時，模糊次數(shù)將是算法的一個限制條件。但是，通常雷達模糊數(shù)不會特別嚴(yán)重，或者通過對目標(biāo)特性的解析或搜索雷達的引導(dǎo)，這一限制條件相對寬松。

2.2 計算機仿真

用大時寬LFM脈沖串信號來驗證Keystone變換 在高速目標(biāo)相參積累中校正包絡(luò)徙動的作用，取雷達參數(shù)為：載頻fc=3.2 GHz，目標(biāo)速度v=5 km/s，信號帶寬B=10 MHz， 信號時寬T=100 μs，采樣率fs=12 MHz，脈沖個數(shù)N=128，脈沖重復(fù)頻率PRF=2 kHz，假設(shè)回波信噪比為-10 dB。圖1為仿真結(jié)果?？梢姡?jīng)過脈沖壓縮后，由于目標(biāo)的運動，在一個相參處理間隔內(nèi)，目標(biāo)跨越了幾十個距離單元，如圖1(a)所示。若如圖1(c)直接方位FFT相參積累，目標(biāo)能量將被分散至多個距離和多普勒單元，不利于目標(biāo)的檢測及測距、測速。對回波經(jīng)Keystone變換后，目標(biāo)回波包絡(luò)被拉回至同一距離，如圖1(b)，然后再進行濾波處理。比較圖1(d)和圖1(c)，不難發(fā)現(xiàn)，經(jīng)過本文的方法處理后，檢測信噪比得到顯著提升。

(a)脈沖壓縮 (b)Keystone變換包絡(luò)校平

(c)直接方位FFT的檢測結(jié)果 (d)欠采樣Keystone變換檢測結(jié)果

圖1 計算機仿真結(jié)果

3 實測數(shù)據(jù)處理及工程實現(xiàn)

3.1 實測數(shù)據(jù)處理

為了進一步驗證算法，利用實測數(shù)據(jù)也進行了處理。針對某S波段雷達回波，信號帶寬5 MHz，回波脈沖數(shù)為64，處理結(jié)果如圖2所示。可見，Keystone變化前后信噪比有3.7 dB改善，目標(biāo)的測距、測速精度均有所提高，表明了方法的有效性。

(1)脈沖壓縮

(2)Keystone變換包絡(luò)校平

(3)直接方位FFT的檢測結(jié)果

(4)欠采樣Keystone變換檢測結(jié)果

3.2 工程實現(xiàn)

ADSP-TS201是一款高性能的靜態(tài)超標(biāo)量處理器，主要結(jié)構(gòu)特點包括并行運算、內(nèi)部存儲器、四指令執(zhí)行、可升級性及多核處理器。運行在600 MHz時，ADSP-TS201內(nèi)核的指令周期為1.67 ns，內(nèi)部存儲器大小為24 Mbit。ADSP-TS201靜態(tài)超標(biāo)量結(jié)構(gòu)使DSP每周期能夠執(zhí)行多達4條指令，執(zhí)行24個16位定點或者6個浮點運算。

通用信號處理板卡是一款自主研發(fā)的通用信號處理硬件平臺。板卡以DSP(ADSP-TS201)為主處理器，F(xiàn)PGA提供時序、控制及數(shù)據(jù)輸入輸出接口，板卡上有8片TS201，并行工作峰值運算能力高達8×3.6GFLOPS。板卡采用的是4片/簇的并行處理結(jié)構(gòu)，每片TS201均通過一路鏈路口與FPGA相連，簇內(nèi)4片TS201彼此通過鏈路互連，單向數(shù)據(jù)率可達500Mb/s。處理結(jié)果通過鏈路口或總線輸出至FPGA，再通過CPCI總線傳送給主控計算機。板卡框圖如圖3～圖4所示。

Keystone變換的步驟如下：

步驟1 對回波信號x(t,n)進行快時間頻率維度的FFT變換得到快時間頻域數(shù)據(jù)X(f,n)；

步驟2 頻域脈壓匹配處理，即Spc(f,n)=X(f,n)P*(f)，其中P*(f)為頻域脈壓匹配權(quán)系數(shù)；

步驟3 對信號Spc(f,n)進行脈沖維度的Keystone變換得到變換后信號Yk(f,m)；

圖3 通用信號處理板卡框圖

圖4 通用信號處理板卡實物圖

步驟4 對變換后信號進行快時間頻域的IFFT處理，得到快時間時域信號yk(t,m)；

步驟5 對IFFT處理后的數(shù)據(jù)進行多普勒濾波。

在工程實際應(yīng)用中，若知道目標(biāo)大致位置，在步驟2后可采取脈壓時域開窗，然后再轉(zhuǎn)換至頻域做后續(xù)計算，以降低運算量。

通常有3種方法實現(xiàn)Keystone變換：(1)DFT+FFT，(2)Chirp-z變換算法，(3)sinc內(nèi)插算法。3種算法性能相當(dāng)，主要是運算量的區(qū)別。假設(shè)快時間距離頻率域采樣點數(shù)為N，雷達在一個CPI內(nèi)發(fā)射的脈沖數(shù)為M。根據(jù)以上介紹，3種實現(xiàn)方法需要的復(fù)數(shù)運算量如下[8]：

DFT+FFT算法：

Chirp-z算法：

Sinc內(nèi)插算法：

M2×N

當(dāng)N及M取值較大時，DFT+FFT方法的運算量很大，sinc內(nèi)插方法次之，Chirp-z方法最小。經(jīng)實際工程應(yīng)用，基于Keystone變換進行空間目標(biāo)處理，信噪比改善效果明顯，計算量可接受，具有很強的工程運用價值。在某工程運用中，跟蹤波門3 km，積累脈沖數(shù)為60，單板可實現(xiàn)同時對8個目標(biāo)的連續(xù)跟蹤，跟蹤數(shù)據(jù)率可達10 Hz。

圖5 實際工程處理效果圖

4 結(jié)束語

本文研究了高速空間目標(biāo)相參積累過程中出現(xiàn)跨距離和跨多普勒頻道問題，給出了目標(biāo)回波模型，分析了拉伸效應(yīng)、脈沖間距離遷徙等因素對積累的影響；采用基于欠采樣下Keystone變換的高速目標(biāo)檢測方法，有效補償距離走動，去除多普勒模糊，提高信噪比，改善雷達的檢測性能。通過仿真和實際數(shù)據(jù)分析，驗證了該算法的有效性，同時結(jié)合硬件平臺，給出了工程實現(xiàn)方法。

[1] 張月輝,朱玉鵬,游鵬,黎湘.高速動態(tài)空間目標(biāo)寬帶回波仿真方法研究[J].系統(tǒng)仿真學(xué),2009(5):2721-2729.

[2] 保錚,邢孟道,王彤.雷達成像技術(shù)[M].北京:電子工業(yè)出版社,2005.

[3] Liu Yimin, Meng Huadong, Zhang Hao, et al. Motion Compensation of Moving Targets for High Range Resolution Stepped-Frequency Radar [ J ]. Sensor, 2008 (5) :3429-3437.

[4] Chen J J, Chen J, Wang S L. Detection of ultra-high speed moving target based on matched Fourier transform[C]//Proceedings of the CIE International Conference on Radar, 2006: 1-4.

[5] Zhang S S, Zeng T, Long T, et al. Dim Target Detection Based on Keystone Transform [C]//Proc. of IEEE Int. Radar Conf. Arlington, 2005: 889-894.

[6] 張順生,曾濤. 基于Keystone變換的微弱目標(biāo)檢測[J]. 電子學(xué)報, 2005(9):1675-1678.

[7] 蘇軍海,李亞超,邢孟道. 窄帶雷達高速多目標(biāo)檢測研究[J]. 西安電子科技大學(xué)學(xué)報(自然科學(xué)版), 2009(6):1003-1009.

[8] 王娟,趙永波. Keystone變換實現(xiàn)方法研究[J]. 火控雷達技術(shù),2011(3):45-51.

Research on space target detection and realization

YUE Mei-jun1, WANG Qi-zhi2

(1. Electronic Engineering Institute of PLA, Hefei 230037;2. No.38 Research Institute of CETC, Hefei 230088)

Generally, high-speed and high maneuvering space targets go through several range cells in the observation period, which makes it difficult to perform coherent power accumulation for target detection. Through the modeling of radar echo signals of high-speed moving targets with large time-width, the influences of stretching effect, pulse-to-pulse range migration and other factors on accumulations are analyzed. A high-speed target detection method based on the Keystone transform in under-sampling is proposed, which can effectively compensate range migration, eliminate Doppler ambiguity, and improve the SNR and radar detection performance. It is verified that the algorithm is effective through the simulation and data analysis, and the engineering realization method is also given by combining hardware platform.

space target; Keystone transform; range migration; engineering realization

2014-07-20；

2014-08-01

岳玫君(1981-)，女，講師，碩士，研究方向：目標(biāo)識別及電大尺寸目標(biāo)電磁特性；王啟智(1981-)，男，高級工程師，研究方向：寬帶陣列信號及高速目標(biāo)檢測。

TN957.51

A

1009-0401(2014)03-0001-05