趙 杰，郭德明，陳 果，倪國新

(南京電子技術研究所, 南京 210039)

·信號處理·

基于改進Keystone變換的高速目標相參積累方法

趙 杰，郭德明，陳 果，倪國新

(南京電子技術研究所, 南京 210039)

Keystone變換常常用于校正運動目標在脈沖多普勒雷達相參積累時間內的距離走動，提高回波信號的信噪比。但是常規(guī)Keystone變換應用于高速目標會存在兩方面問題，首先常規(guī)Keystone變換引入插值增加了運算負擔，其次若存在多普勒模糊，常規(guī)Keystone變換實現方法會出現“半盲速點”效應。文中提出了基于尺度變換的Keystone變換方法，新方法可以明顯減少運算量，減小了帶寬對多普勒模糊的影響，抑制了“半盲速點”效應，提高了相參積累效果。

距離走動； Keystone變換；尺度變換；“半盲速點”效應

0 引 言

隨著雷達目標的雷達散射截面積(RCS)不斷減小，超高速飛行器的出現，這些都對雷達的檢測提出了新的挑戰(zhàn)，成為了雷達急需解決的技術問題之一。文獻[1-2]提出了通過增加信號積累時間提高雷達的檢測能力，但是高速目標的寬帶回波信號包絡在積累時間內存在距離單元走動現象，增加了目標檢測的復雜度。文獻[3]提出了一種基于Hough變換的長時間積累算法，但是算法復雜且是非相參積累，效率較低。相參積累可以改善積累效率，但是必須解決積累時間內的距離走動問題。

針對高速目標在相參積累時間內出現的距離走動現象和高速目標存在的多普勒模糊問題，文獻[4]提出的Keystone變換有效地校正了脈沖多普勒(PD)雷達的距離走動。該Keystone變換方法引入了sinc插值，運算量大，不利于工程實現。文獻[5]和文獻[6]提出了基于Chirp_z算法的Keystone變換，運用快速傅里葉變換和復乘運算，優(yōu)化了運算，但文獻[6]指出該方法在回波信號的信噪比較低時，相參積累效果明顯下降。針對相參積累受限于目標運動，文獻[7]提出一種基于Keystone變換的運動補償方案，可以在沒有目標運動速度信息條件下校正距離走動，從而使積累時間不再受目標運動的限制，文獻[8]提出了“半盲速點”效應抑制和消除的方法，還需進一步優(yōu)化引入插值運算帶來的問題。針對Keystone變換在工程實踐中遇到的計算量問題和“半盲速點”效應，本文將文獻[9]提出的尺度變換應用于Keystone變換，該方法通過在慢時間域的復乘和快速傅里葉變換運算，消除了插值運算，解決了常規(guī)Keystone變換中運算量的問題，并結合“半盲速點”效應抑制算法，減小了帶寬對多普勒模糊倍數區(qū)間的影響，能夠快速的實現對距離走動的校正和對多普勒模糊的補償，有效提高相參積累處理后信號的信噪比。

1 Keystone原理分析

本節(jié)將闡述Keystone變換意義，針對多普勒模糊問題改進Keystone變換，結合常用Keystone變換方法，分析基于尺度變換的Keystone變換的意義。

1.1 Keystone變換算法

假設雷達發(fā)射信號為線性調頻信號，則第n個子脈沖發(fā)射信號表達式為

exp(j2πf0(t-nTr))

(1)

式中：t為時間；T0為線性調頻信號的脈沖寬度；b為調頻斜率；Tr為脈沖重復周期；f0為發(fā)射信號載頻。令tn=nTr，稱為慢時間，t′=t-tn，稱為快時間，速度為v的點目標的初始距離為R0，回波信號經下變頻處理后為

s(t′,tn)=Aa(t′-tn)exp(-j2πf0tn)

(2)

其中

為第n個脈沖的延遲時間。

對式(2)在快時間域做傅里葉變換,得到線性調頻信號的頻譜

s(f,tn)=Aa(f)exp(-j2π(f0+f)τn)=

Aa(f)exp(-j4π(f0+f)R0/C)·

exp(j4π(f0+f)vtn/C)=

Ap(f)exp(j4π(f0+f)vtn/C)

(3)

其中

p(f)=a(f)exp(-j4π(f0+f)R0/C)

(4)

式中：f為距離頻率；指數項exp(j4πfvtn/C)表示距離頻率和多普勒頻率之間存在耦合，造成各個脈沖周期內的脈壓結果不同，引起各個點目標回波信號的線性距離走動。通過對慢時間軸作變尺度變換，可以消除這一耦合，即Keystone變換。則

(5)

式中：tk為Keystone變換之后的虛擬慢時間。將式(5)代入式(4)，可得

s(f,tk)=Ap(f)exp(j4πf0vtk/C)

(6)

由式(6)可以看出，經過Keystone變換，回波信號的距離頻率與目標速度的耦合相位得到了補償。

高速目標常常會引入多普勒模糊問題，通常回波信號在多普勒模糊倍數已知的情況下，對式(6)乘以多普勒模糊補償因子即可

(7)

但Keystone變換在長時間相參積累過程中存在“半盲速點”效應，除去插值影響Keystone性能外，信號頻率對多普勒模糊影響如式(8)

(8)

式中：L為多普勒模糊倍數，若目標速度恰處于“半盲速點”附近，多普勒模糊倍數會出現突變現象，嚴重影響目標回波信號的積累效果。

為了減小帶寬對多普勒模糊的影響，抑制“半盲速點”效應，Keystone變換前應先進行多普勒模糊補償，即

p(f)exp(-j4π(f0+f)v0tn/C)

(9)

1.2Keystone變換實現的改進

本節(jié)將介紹sinc插值算法和Chirp_z算法，分析基于尺度變換的Keystone變換的原理。

常規(guī)的Keystone變換利用sinc插值算法實現信號在慢時間軸作變尺度變換，實現距離頻率和多普勒頻率的解耦合，其中m,n分別為處理前后慢時間的采樣序列值，即

(10)

在多普勒模糊的情況下，對高速目標回波信號進行多普勒模糊補償

(11)

Chirp_z算法應用于采樣間隔不均勻時快速求取各采樣值。文獻[3]提出的Chirp_z算法實現了信號在虛擬慢時間tm對應的頻域上進行變尺度離散傅里葉變換(DFT)變換，用循環(huán)卷積計算線性卷積，進而可以在頻域使用快速傅里葉變換(FFT)實現快速計算，降低了插值方式下的Keystone變換的運算量。

文獻[10]提出在聚束式SAR成像處理中，在距離域和方位域的重采樣中應用尺度變換算法，有效地解決了散射點在距離單元和多普勒單元的運動補償，該算法能夠有效地實現時間維的尺度變換。本文將尺度變換引入PD雷達檢測中，實現Keystone變換，為了抑制“半盲速點”效應，需先對回波信號進行多普勒模糊補償，改進基于Keystone變換的多普勒模糊處理。

圖1所示，基于尺度變換算法的Keystone變換可在時域完成變尺度變換，有效地處理“半盲速點”效應，其中n表示回波信號的慢時間采樣序列，f為相應于快時間t的距離頻率的采樣序列。

圖1 基于尺度變換的Keystone變換流程圖

回波信號需要先進行多普勒模糊補償，減小多普勒模糊倍數突變區(qū)間，減少帶寬對多普勒模糊的不利影響，較好地抑制“半盲速點”效應，提高距離走動校正后的相參積累效果。即

Scomp(f,n)=S(f,n)·Φcon(f,n)

(12)

式中：Φcon(f,n)為多普勒模糊補償因子，即

(13)

此外，流程圖中的兩次復乘運算的算子，分別稱為尺度因子和反尺度因子，其中fn為慢時間tn對應的多普勒頻率，即

(14)

(15)

如流程圖所示，在完成時間域的尺度變換過程中，系統(tǒng)還需要進行兩次相位濾波，此過程可以通過快速傅里葉變換實現，濾波器的系統(tǒng)函數h1(f,fn)和h2(f,fn)定義如下

h1(f,n)=wa(n)exp(jnkαn2)

(16)

h2(f,n)=exp(-jπδKαn2)

(17)

式中：wa(t)是一個可選的權值因子，通常使用泰勒窗或漢明窗。系統(tǒng)函數h1(f,n)和h2(f,n)需滿足大時間帶寬積，其中式(16)滿足大時間帶寬積，且尺度變換因子應近似一。

參數δ取值由尺度變換因子決定，即

(18)

參數Kα取值受限于式(16)和式(17)的時帶寬積，此外，尺度變換的中間結果頻譜會對處理效果產生不利影響，故使其位于相位濾波的系統(tǒng)函數的帶寬范圍內。

本文用復數乘法次數表征計算量，設距離頻率f的采樣點數為N，相參積累時間內的脈沖個數為M，則一次sinc插值的計算量為NM2,一次相參積累的計算量為NM·(lbM)/2,一次模糊校正的計算量為NM，則上述三種Keystone變換實現算法的所需計算量如表1所示。

表1 Keystone變換實現算法比較

尺度變換算法在計算量方面相較于sinc插值算法優(yōu)化明顯，但較于Chirp_z算法優(yōu)化不明顯，但是回波信號的信噪比低于-16 dB時，Chirp_z算法相參積累出現明顯下降，尺度變換算法較Chirp_z算法提高了在低信噪比情況下的處理效果，其可以在對距離走動作有效校正的前提下優(yōu)化計算，有利于工程應用。

2 仿真分析

2.1 常規(guī)Keystone變換實現距離校正

仿真所需的雷達參數如下，載頻f0=1 GHz，信號帶寬B=30 MHz，脈沖寬度T0=100 μs，脈沖重復周期T=1 ms，一個CPI內發(fā)射的脈沖數np=128，目標相對雷達初始距離R0=35 km，信號的初始信噪比設為-16 dB，目標的徑向運動速度v=310 m/s。

圖2 Keystone變換校正前后回波信號的等高線對比圖

目標在128個脈沖周期，即積累時間12.8 ms內的距離走動量δr=39.68 m，距離分辨率ΔR=5 m，由圖2顯示回波信號在第1個脈沖周期和第 128個脈沖周期之間跨距離單元數約為8個，與理論分析一致，圖3顯示經Keystone變換校正后，信號相參積累后的幅值得到很大提高，相參積累后信號的信噪比達18.7 dB,較回波信號直接進行相參積累的信噪比提高大約6.5 dB。

圖3 Keystone變換校正前后回波信號相參積累效果對比圖

2.2 “半盲速點”效應分析

仿真所需的雷達參數如下，載頻f0=1 GHz，信號帶寬B=10 MHz，脈沖寬度T0=100 μs，脈沖重復周期T=1 ms，一個CPI內發(fā)射的脈沖數np=128，目標相對雷達初始距離R0=35 km，信號的初始信噪比設為-16 dB，目標的徑向運動速度v=370 m/s。

圖4比較了“半盲速點”效應對高速目標回波信號相參積累性能的影響，從圖4中可知“半盲速點”效應嚴重影響了Keystone變換的效果，通過對“半盲速點”效應做抑制處理后，信號相參積累后的信噪比達到20.4 dB，較未做抑制處理的結果提高了2 dB。因此，在校正距離走動的過程中，Keystone變換處理多普勒模糊問題應采用上述抑制“半盲速點”效應的方法，即先對多普勒模糊做補償，再進行慢時間域的變換，可以提高Keystone變換的魯棒性。

圖4 “半盲速點”效應處理前后相參積累效果

2.3 Keystone變換實現方式比較

仿真所需的雷達參數如下，載頻f0=1 GHz，信號帶寬B=10 MHz，脈沖寬度T0=100 μs，脈沖重復周期T=1 ms，一個CPI內發(fā)射的脈沖數np=128，目標相對雷達初始距離R0=35 km，信號的初始信噪比設為-16 dB，目標的徑向運動速度v=310 m/s。

圖5比較了采用sinc插值和尺度變換兩種算法實現Keystone變換后信號相參積累的效果。圖5的上圖是采用文獻[3]提出的插值算法，未處理 “半盲速點”效應，下圖采用了尺度變換算法，其處理后信號相參積累的信噪比大約為19.6 dB，運行時間為29.92 s，較插值提高了1 dB，節(jié)約了110 s左右。

圖5 Sinc插值和尺度變換實現Keystone變換的效果

圖6比較了采用Chirp_z算法和尺度變換算法實現Keystone變換后信號相參積累的效果。圖6的上圖是采用文獻[4]提出的Chirp_z算法，未處理 “半盲速點”效應，下圖采用了尺度變換算法，其處理后信號相參積累的信噪比大約為19.63 dB，信噪比較Chirp_z算法提高了3 dB。

圖6 Chirp_z和尺度變換實現Keystone變換的效果

3 結束語

雷達目標檢測的效果取決于回波信號信噪比，相參積累可以有效提高低速目標回波信號的信噪比，但高速目標的回波信號在長時間積累過程中會出現距離走動現象，造成回波信號相參積累后信噪比下降，為了提高信號相參積累性能，本文應用Keystone變換校正高速目標回波信號處理的距離走動問題，但是傳統(tǒng)Keystone變換常常存在計算量大、“半盲速點”效應引起多普勒模糊補償性能下降等問題，本文針對上述問題，應用尺度變換算法實現Keystone變換，優(yōu)化了計算，提高了對微弱回波信號的相參積累性能，并有效結合抑制“半盲速點”效應處理算法，有效補償了在半盲速點附近的多普勒模糊，經仿真分析，改進后的Keystone變換算法可以有效地解決高速目標信號的距離走動問題，提高信號經相參積累處理后的信噪比，為目標檢測奠定了基礎。

[1] 李 海, 吳嗣亮, 莫 力, 等. 微弱信號長時間積累的檢測方法[J]. 北京理工大學學報, 2001, 21(5): 614-617. LI Hai, WU Siliang, MO Li, et al. A method for long-term signal integral detection of weak targets[J]. Journal of Beijing Institute of Technology, 2001, 21(5): 614-617.

[2] 劉國靖, 王宏軍, 秋 實. 短脈沖雷達模擬實驗系統(tǒng)及目標探測初步實驗研究[J]. 電子機械工程, 2001(3): 2-12. LIU Guojing, WANG Hongjun, QIU Shi. Short pulse radar simulation experimental system and basic experimental research on the target detection[J]. Electro-Mechanical Engineering, 2001(3): 2-12.

[3] MO Li. Radar detection of range migrated weak target through long-term integration[J]. Acta Electronica Sinica, 2003(4): 539-544.

[4] 張順生, 曾 濤. 基于keystone變換的微弱目標檢測[J]. 電子學報, 2005, 33(9): 1675-1678. ZHANG Shunsheng, ZENG Tao. Weak target detection based on keystone transform[J]. Acta Electronica Sinica, 2005, 33(9): 1675-1678.

[5] 王 娟, 趙永波. Keystone變換實現方法研究[J]. 火控雷達技術, 2011, 40(1): 45-51. WANG Juan, ZHAO Yongbo. Research on implementation of keystone transform[J]. Fire Control Radar Technology, 2011, 40(1): 45-51.

[6] 王 霞. 高速目標的跨距離門長時間積累研究[D]. 南京: 南京理工大學, 2014. WANG Xia. Research on the integration of high-velocity target across range-gate inlong time[D]. Nanjing: Nanjing University of Science and Technology, 2014.

[7] 張順生, 曾 濤. 基于Keystone變換的微弱目標檢測[J]. 電子學報, 2005, 33(9): 141-144. ZHANG Shunsheng, ZENG Tao. Weak target detection based on keystone translorm[J]. Acta Electronica Sinic, 2005, 33(9): 141-144.

[8] 洪永彬, 高梅國, 王俊嶺, 等. Keystone變換半盲速點效應的抑制和消除[J]. 電子與信息學報, 2014, 36(1): 175-180. HONG Yongbin, GAO Meiguo, WANG Junling, et al. The supression and elimination of half blind velocity effect associated with keystone transform[J]. Journal of Electronics & Information Technology, 2014, 36(1): 175-180.

[9] ZHU D, LI Y, ZHU Z. A keystone transform without interpolation for SAR ground moving-target imaging[J]. IEEE Geoscience & Remote Sensing Letters, 2007, 4(1): 18-22.

[10] ZHU D, YE S, ZHU Z. Polar format agorithm using chirp scaling for spotlight SAR image formation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2008, 44(4):1433-1448.

趙 杰 男，1990年生，碩士研究生。研究方向為雷達目標檢測。

A Method of High-velocity Target Coherent Integration Based on Improved Keystone Transform

ZHAO Jie，GUO Deming，CHEN Guo, NI Guoxin

(Nanjing Research Institute of Electronics Technology, Nanjing 210039, China)

Keystone transform can correct range migration on the moving target in the PD radar coherent integration time, and improve the signal-to-noise ratio of the echo signal. But conventional Keystone transform applied to the high-velocity target will have two problems. Firstly, conventional Keystone transform increases computation burden with interpolation. Secondly, with Doppler ambiguity, 'half-blind-velocity effect' occurs when the conventional implementation method of Keystone transform is applied. The paper presents novel Keystone transform based on chirp scaling. It can obviously reduce computation cost and influence of bandwidth on Doppler ambiguity, suppress 'half-blind-velocity' and improve the performance of coherent integration.

range migration; keystone transform; scaling principle; half-blind-velocity effect

10.16592/ j.cnki.1004-7859.2016.04.008

趙杰 Email：zhaojie657@sina.com

2015-11-13

2016-01-27

TN911.7

A

1004-7859(2016)04-0034-04