王澤玉,李 明,盧云龍,陳洪猛,曹潤清,左 磊,張 鵬
(西安電子科技大學(xué) 雷達(dá)信號處理國家重點實驗室,西安 710071)
?
基于頻域校正的快速長時間積累算法
王澤玉,李明,盧云龍,陳洪猛,曹潤清,左磊,張鵬
(西安電子科技大學(xué) 雷達(dá)信號處理國家重點實驗室,西安 710071)
摘要:針對長時間積累過程中回波信號出現(xiàn)的距離徙動和多普勒走動問題,提出一種基于 頻域校正的快速長時間積累算法。該算法首先利用Radon-Ambiguity 變換對脈壓以后的回波信號 進(jìn)行加速度估計,用得到的估計值構(gòu)造加速度補償因子消除距離彎曲和多普勒走動,然后選取回 波包絡(luò)之間的相關(guān)系數(shù)作為代價函數(shù)對目標(biāo)的速度進(jìn)行估計,用得到的估計值構(gòu)造速度補償因子 消除距離走動。仿真結(jié)果驗證了該算法的有效性。
關(guān)鍵詞:長時間積累; 距離徙動; 多普勒走動; 補償因子; 頻域校正
0引言
長時間積累技術(shù)能改善回波的信噪比,是現(xiàn)代雷達(dá)提高檢測性能的一種有效途徑[1]。 但隨著積累時間的增加,運動目標(biāo)回波容易產(chǎn)生距離徙動(包括距離走動和距離彎曲)和多普勒走動現(xiàn)象,使得目標(biāo)能量擴(kuò)散,回波得不到有效積累[2-5]。 因此,如何進(jìn)行有效的包絡(luò)補償和多普勒補償是長時間積累技術(shù)亟待解決的問題[6-8]。
針對運動目標(biāo)長時間積累過程中出現(xiàn)的距離走動問題,文獻(xiàn)[9]提出一種基于約當(dāng)傅立葉變換(Radon-Fourier Transform,RFT)的運動目標(biāo)檢測方法; 文獻(xiàn)[10]提出一種基于分時處理的相干積累及Keystone變換校正相結(jié)合的距離徙動實時相干積累算法,實現(xiàn)高速目標(biāo)的有效檢測; 文獻(xiàn)[11]提出一種改進(jìn)Keystone變換算法。 但以上算法只適用于目標(biāo)作勻速運動時的情況,當(dāng)目標(biāo)作勻加速運動并在長時間積累過程中產(chǎn)生距離徙動和多普勒走動時,這些算法均失效。 針對這個問題,文獻(xiàn)[12]將Keystone變換與Wigner-Hough變換進(jìn)行結(jié)合,提出一種用于機載相控陣預(yù)警雷達(dá)實現(xiàn)高速空中機動目標(biāo)檢測的新方法。 然而當(dāng)目標(biāo)速度較大,產(chǎn)生多普勒模糊時,該算法需要對模糊數(shù)進(jìn)行搜索,計算量較大,并且采用Wigner-Hough變換對加速度進(jìn)行搜索,計算量進(jìn)一步加大。
為減小計算量,本文提出一種基于頻域校正的快速長時間積累算法,該算法首先采用約當(dāng)模糊變換(Radon Ambiguity Transform, RAT)估計目標(biāo)的加速度,構(gòu)造加速度補償因子,補償加速度引起的距離彎曲和多普勒走動,再選取計算量小的代價函數(shù)對目標(biāo)速度進(jìn)行搜索,構(gòu)造速度補償因子,補償速度引起的距離走動。 實驗結(jié)果表明本文算法能有效實現(xiàn)包絡(luò)補償和多普勒補償,且與文獻(xiàn)[12]相比,計算量明顯減少。
1回波信號模型
(1)
(2)
(3)
2基于頻域校正的快速長時間積累算法
(4)
由以上分析可知,距離徙動和多普勒走動是由目標(biāo)速度和加速度引起的回波相位變化,因此,本文構(gòu)造相位補償因子對加速度相位項和速度相位項進(jìn)行補償,從而消除距離徙動和多普勒走動。
2.1加速度相位項補償
信號s(t)的模糊函數(shù)定義為
AFs(τ,ξ)=∫-∞∞s(t+τ/2)s*(t-τ/2)e-j2πξtdt
(5)
式中: []*為共軛運算; (τ,ξ)平面為模糊域。 當(dāng)信號s(t)為線性調(diào)頻信號時,s(t)的模糊函數(shù)為
AFs(τ,ξ)=A2δ(ξ-μ0τ)exp(j2πf0τ)
(6)
式中: f0為中心頻率; μ0為調(diào)頻率。 對線性調(diào)頻信號的模糊函數(shù)取模以后作Radon變換,得到
(7)
由式(7)可知,當(dāng)觀測信號中含有調(diào)頻率為μ0的線性調(diào)頻信號時,檢測器在μ=μ0點出現(xiàn)峰值,可通過峰值的位置估計線性調(diào)頻信號的調(diào)頻率。
由式(4)可知,當(dāng)頻率f為一個定值時(f=fconst,即選定一個距離單元),信號Sr(fconst,tm)是關(guān)于方位慢時間的線性調(diào)頻信號。 該線性調(diào)頻信號的調(diào)頻率是目標(biāo)加速度的函數(shù)。 因此,可以對Sr(fconst,tm)進(jìn)行RAT來檢測該線性調(diào)頻信號的調(diào)頻率,從而估計出運動目標(biāo)的加速度。
設(shè)對信號Sr(fconst,tm)進(jìn)行RAT以后,搜索得到的峰值對應(yīng)的弧度為θ0,則θ0對應(yīng)的線性調(diào)頻信號的歸一化調(diào)頻率可以表示為[14]
(8)
式中:M為脈沖數(shù);Tr為脈沖重復(fù)周期。 同時,由式(4)可計算得理論上線性調(diào)頻信號Sr(f,tm)對應(yīng)的歸一化調(diào)頻率為
(9)
聯(lián)立式(8)~(9)可以求得運動目標(biāo)加速度的估計值為
(10)
(11)
由式(11)可以明顯看出,加速度引起的距離彎曲和多普勒走動均得到補償,而且由于RAT只是模糊平面上的一維搜索,與Wigner-Hough變換法的二維搜索相比,計算量大大減少[14-15]。 但是回波距離走動的指數(shù)相位項仍然存在,接下來需要對目標(biāo)速度引起的距離走動進(jìn)行補償。
2.2速度相位項補償
為了對距離走動進(jìn)行校正,本文選取第一個脈沖和最后一個脈沖的相關(guān)系數(shù)為代價函數(shù),對目標(biāo)速度進(jìn)行搜索,然后構(gòu)造關(guān)于速度的補償因子項,來完成距離走動校正。
將速度范圍(目標(biāo)速度一般不超過1 000 m/s)以相等的步長量化為N個值(步長Δv≤ΔR/2T。 其中,ΔR為距離分辨單元,T為積累時間。 同時為了減少運算量,步長不應(yīng)該選取的太小,即步長Δv的選取是量化誤差和計算量的折衷),然后根據(jù)N個速度值確定N個相位補償因子進(jìn)行搜索補償,具體步驟如下所示:
本文通過構(gòu)造兩個補償因子來實現(xiàn)包絡(luò)補償和多普勒補償,為獲得補償因子,對運動目標(biāo)的加速度和速度進(jìn)行估計,具體實現(xiàn)流程如圖1所示。
圖1基于頻域校正的勻加速目標(biāo)長時間積累算法流程圖
3仿真實驗
雷達(dá)仿真參數(shù)的設(shè)置如表1所示。
表1 雷達(dá)仿真參數(shù)
目標(biāo)參數(shù): 設(shè)目標(biāo)速度為690 m/s; 加速度為50 m/s2; 初始距離為10.5 km; 信噪比SNR(脈壓前)為0 dB。 目標(biāo)的速度搜索范圍選取為(0 m/s, 1 000 m/s),速度的搜索步長選為18 m/s。
將1 024個脈沖回波信號進(jìn)行頻域脈沖壓縮,并變換到時域以后得到的結(jié)果如圖2所示,由圖可以看出,由于在長時間積累過程中有速度和加速度的存在,回波存在嚴(yán)重的距離徙動。 對圖2的結(jié)果直接進(jìn)行相參積累,得到的結(jié)果見圖3。 由圖3可以看出,直接進(jìn)行相參積累時,能量在多普勒域展寬嚴(yán)重,積累效果很差。 因此,必須進(jìn)行加速度和速度的估計并對回波信號進(jìn)行補償。
圖2 原始脈壓回波數(shù)據(jù)
圖3直接相參積累結(jié)果
頻域脈壓信號在一個距離單元上進(jìn)行RAT得到的結(jié)果如圖4(a)所示。 圖中出現(xiàn)一個明顯的峰值,將峰值對應(yīng)的角度代入式(10),計算得到加速度的估計值。 對頻域信號進(jìn)行加速度相位項補償后,利用N個量化速度進(jìn)行速度補償,得到N個量化速度對應(yīng)的代價函數(shù)如圖4(b)所示。 圖中相關(guān)系數(shù)的峰值對應(yīng)的速度即為速度的估計值。
圖4頻域補償因子搜索示意圖
文獻(xiàn)[12]的算法與本文算法對比的結(jié)果如圖5(a)~(d)所示。 圖5(a)~(b)分別為利用文獻(xiàn)[12]算法和利用本文算法補償后的時域回波結(jié)果,由圖可知,兩種算法補償后包絡(luò)都基本對齊。 此時分別對兩種算法校正后的回波信號作相參積累,得到的結(jié)果如圖5(c)~(d)所示。 與圖3對比可知,兩種算法校正后能量聚集到一個距離單元,多普勒擴(kuò)散也已基本消除,積累效果明顯改善。 實驗表明,兩種算法均能有效地校正加速度和速度引起的距離徙動和多普勒走動。
圖5 長時間積累結(jié)果對比
表中RAT為RAT算法;WH為Winger-Hough算法;N為快時間采樣點數(shù);M為脈沖積累數(shù);L為2的整數(shù)冪次,且有L>2M-1;S為搜索次數(shù)。 Winger-Hough算法是時頻平面的二維搜索,而RAT算法是模糊平面的一維搜索。 因此,由仿真結(jié)果和算法復(fù)雜度分析可知,本文算法在有效校正距離徙動和多普勒走動的情況下,明顯減少計算量。
圖6長時間積累結(jié)果對比
由圖6(a)可知,隨著信噪比的增加,目標(biāo)加速度的均方誤差呈下降趨勢,且RMSE均不大于0.12 m/s2。 由圖6(b)可知,隨著信噪比的增加,目標(biāo)速度的RMSE下降,且當(dāng)信噪比大于 -7 dB時,目標(biāo)速度的RMSE為6 m/s并保持穩(wěn)定。 說明本文算法在低SNR條件下具有一定的穩(wěn)健性。
4結(jié)論
本文針對長時間積累算法計算量較大的問題,提出一種不需要考慮多普勒模糊的頻域校正快速算法。 該算法通過RAT對加速度進(jìn)行一維搜索,構(gòu)造加速度補償因子項消除加速度引起的距離彎曲和多普勒走動; 以回波之間的相關(guān)系數(shù)為代價函數(shù)對速度進(jìn)行搜索,構(gòu)造速度補償因子項消除速度引起的距離走動。 仿真結(jié)果表明,本文算法能有效地實現(xiàn)包絡(luò)補償和多普勒補償,并具有較小的計算量,因而更加適合工程應(yīng)用。
參考文獻(xiàn):
[1] 趙永波, 周曉佩, 王娟. 一種用于弱信號檢測的廣義Keystone變換算法[J]. 西安電子科技大學(xué)學(xué)報, 2013, 40(2): 98-102.
[2] 葉春茂, 魯耀兵, 宋建社, 等. 一種ISAR目標(biāo)轉(zhuǎn)動補償及圖像尺寸標(biāo)定方法[J]. 電子學(xué)報, 2014, 42(3): 611-615.
[3] Su Junhai, Xing Mengdao, Wang Genyuan, et al. High-Speed Multi-Target Detection with Narrowband Radar[J]. IET Radar Sonar and Navigation, 2010, 4(4): 595-603.
[4] 張衛(wèi)杰, 高昭昭, 許博, 等. 基于Keystone 變換的警戒雷達(dá)信號處理[J]. 系統(tǒng)工程與電子技術(shù), 2011, 33(9): 2007-2011.
[5] 寧夏, 葉春茂, 楊健. 一種基于時頻匹配的長時間相干積累器設(shè)計[J]. 系統(tǒng)工程與電子技術(shù), 2013, 35(11): 2275-2282.
[6] Yang Jun, Sun Guangcai, Xing Mengdao, et al. Squinted TOPS SAR Imaging Based on Modified Range Migration Algorithm and Spectral Analysis[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(10): 1707-1711.
[7] 鄭紀(jì)彬, 符渭波, 蘇濤, 等. 一種新的高速多目標(biāo)檢測及參數(shù)估計方法[J]. 西安電子科技大學(xué)學(xué)報, 2013, 40(2): 82-88.
[8] Qian Jiang, Lv Xiaolei, Xing Mengdao, et al. Motion Parameter Estimation of Multiple Ground Fast-Moving Targets with a Three-Channel Synthetic Aperture Radar[J]. IET Radar Sonar and Navigation, 2011, 5(5): 582-592.
[9] Xu Jia, Yu Ji, Peng Yingning, et al. Radon-Fourier Transform for Radar Target Detection, I: Generalized Doppler Filter Bank[J]. IEEE Transactions on Aerospace and Electronic Systems, 2011, 47(2): 1183-1202.
[10] 關(guān)欣, 胡東輝, 仲利華, 等. 一種高效的外輻射源雷達(dá)高徑向速度目標(biāo)實時檢測方法[J]. 電子與信息學(xué)報, 2013, 35(3): 581-588.
[11] 王娟, 趙永波. 一種改進(jìn)的Keystone 變換算法及其在微弱信號檢測中的應(yīng)用[J]. 航空兵器, 2011(5): 3-6.
[12] 吳仁彪, 賈瓊瓊, 李海. 機載雷達(dá)高速空中機動目標(biāo)檢測新方法[J]. 電子學(xué)報, 2013, 41(1): 86-90.
[13] 保錚, 邢孟道, 王彤. 雷達(dá)成像技術(shù)[M]. 北京: 電子工業(yè)出版社, 2005.
[14] 蘇軍海, 李亞超, 邢孟道, 等. 采用Radon模糊變換的寬帶雷達(dá)多目標(biāo)檢測方法[J]. 西安交通大學(xué)學(xué)報, 2009, 43(4): 85-89.
[15] Sharif R, Abeysekera S S. Efficient Wideband Signal Parameter Estimation Using a Radon-Ambiguity Transform Slice[J]. IEEE Transactions on Aerospace and Electronic Systems, 2007, 43(2): 673-688.
[16] Grossi E, Lops M, Venturino L. A Novel Dynamic Programming Algorithm for Track-before-Detect in Radar Systems[J]. IEEE Transactions on Signal Processing, 2013, 61(10): 2608-2619.
A Fast Long-Time Accumulation Algorithm Based on Frequency-Domain Correction
Wang Zeyu, Li Ming, Lu Yunlong, Chen Hongmeng, Cao Runqing, Zuo Lei, Zhang Peng
(National Key Laboratory of Radar Signal Processing, Xidian University,Xi’an 710071,China)
Abstract:In order to solve range migration and Doppler walk caused by target motion during long-time coherent integration period, a novel fast algorithm based on frequency-domain correction is proposed. The Radon-Ambiguity transformation is utilized to estimate the acceleration and the acceleration compensation factor is constructed to correct range bend and Doppler walk. Then the relevance coefficient between the two echoes is used to estimate the velocity,and the velocity compensation factor is constructed to correct range walk. The simulation results prove the effectiveness of the proposed algorithm.
Key words:long-time coherent integration; range migration; Doppler walk; compensation factor; frequency-domain correction
中圖分類號:TN958. 5
文獻(xiàn)標(biāo)識碼:A
文章編號:1673-5048( 2016) 02-0061-05
作者簡介:王澤玉(1990-),女,河南許昌人,博士研究生,研究方向為雷達(dá)微弱目標(biāo)檢測。
基金項目:國家自然科學(xué)基金項目(61271297; 61272281; 61301284); 國防預(yù)研基金項目(9140A07020913DZ01001)
收稿日期:2015-08-06
DOI:10.19297/j.cnki.41-1228/tj.2016.02.012