林春风 黄春琳 粟 毅
(国防科学技术大学电子科学与工程学院 长沙 410073)
双基地雷达Radon-Fourier变换弱目标积累检测
林春风*黄春琳 粟 毅
(国防科学技术大学电子科学与工程学院 长沙 410073)
增加信号相参积累时间能够提高弱目标检测能力。双基地雷达长时间相参积累的关键是解决非线性相位回波的目标运动补偿问题。在波束影响可忽略的情况下,该文提出了双基地雷达距离-速度域Radon-Fourier变换(RFT)处理方法。该方法通过频域联合搜索双基地距离和速度参数空间,完成距离走动校正与脉冲积累。理论分析和仿真实验证明了其有效性。
双基地雷达;相参积累;RFT变换
低RCS的弱目标检测是当前雷达领域面临的一个技术难题[1]。增加信号相参积累时间能够提高弱目标检测信噪比,但目标运动引起的距离走动会严重影响目标能量的积累和检测。因此,必须减小或消除距离走动对目标能量积累的影响。研究目标包络距离走动问题的雷达信号积累方法有很多[2],其中积分类方法是一个研究热点。积分类方法包括相参Hough变换[3]、相参Radon变换[4,5]等,文献[6–8]已将这类方法统一为Radon-Fourier变换(RFT),其本质是对目标运动参数决定的直线进行积分的同时增加对目标的多普勒滤波。
然而,现有的相参积累算法多数以单基地雷达目标检测为应用背景,单基地雷达主要接收目标后向散射回波进行相参积累,除时间滤波外很难增加其他空间信息,对弱目标的检测效果和前景有限。双基地雷达能够利用目标侧向或前向散射进行探测[9,10],有优于单基地雷达的目标信息获取能力。此外,增加接收站组成双基地雷达系统还能够增大探测范围和作用距离,增强系统抗干扰和战场生存能力。
本文提出了双基地雷达RFT弱目标积累检测方法,其中第2节介绍了双基地雷达回波模型,第3节分析了双基地雷达RFT算法的特点,给出了在距离-速度域的具体实现方法,第4节为仿真验证结果。
本文主要讨论T/R-R类型的双基地雷达,这种双基地雷达容易在现有雷达基础上,通过增加一个独立接收站构成。其模型如图1所示,雷达T/R位置(–L, 0),独立接收站R位置(L, 0), β为双基地角,δ为目标方向角。
图1 T/R-R双基地雷达几何模型Fig. 1 T/R-R bistatic radar geometric model
设系统工作在侧向区,雷达发射线性调频信号
目标以恒速vT运动,初始时刻与发射站、接收站的斜距分别为RT0, RR0。将接收站接收的直达波作为参考信号[11],那么瞬时双基地距离和可定义为:
为便于分析,忽略目标起伏特性影响,接收站接收到的目标基带回波信号可表示为
其中,τ=nTs为快时间,Ts是快时间采样周期,n=0, 1, ..., Nr–1, Nr为距离窗采样数。t=mTr为慢时间,Tr为脉冲重复周期,m=0, 1, ..., M–1, M为相参积累时间内发射的脉冲数。f0为发射信号载频,c为光速,σ0为目标双基地后向散射系数,w(t)为白噪声。
忽略目标距离的2阶以上分量,式(2)可近似为:
双基地雷达的距离分辨单元和多普勒分辨单元分别为:
若积累时间T内不发生多普勒扩展和距离弯曲,即满足
则目标与雷达间的瞬时距离和可表示为瞬时斜距rs(t)=r0+v0t。忽略指数常数项,回波快时间域的匹配滤波输出为:
令r=cτ为目标距离搜索参数,上式可重写为:
由式(10)可知,经过脉冲压缩处理后目标回波包络走动是一条由目标运动参数(r0, v0)决定的斜线,如图2所示。当目标在Δts时间内移动距离和Δr>ρr时,即产生跨距离单元问题。
图2 双基地距离轨迹Fig. 2 Bistatic distance trajectory
3.1 双基地RFT变换原理
标准RFT用于信号检测的基本思想是将距离-脉冲域空间带有相位信息的直线信号轨迹映射成某参数空间的一个点。由回波模型推导过程可知,双基地雷达可在规定的积累时间内采用标准RFT变换实现相参积累。利用不同的速度搜索值构造多普勒滤波器
与标准RFT算法参数(fd, v, θ)在信号空间中构成一组映射关系相似,双基地系统中,如果目标近似具有1阶运动特性,直线的截距r0与极距ρ的关系仍可表述为ρ=r0sinθ,斜率v0与极角θ的关系表述为v0=–cotθ。因此,可将式(12)从(ρ, θ)域变换到(r, v)域,即
3.2 双基地RFT变换距离-速度域算法实现
双基地RFT变换距离-速度域离散表达式为:
在式(14)中RFT变换通过联合搜索参数空间中目标参数,解决了距离走动(弯曲)与相位调制耦合的问题。由于雷达数据为离散值,参数空间搜索步长通常为整数,这种方法会引入由量化误差带来的能量损失。本节结合频域匹配滤波技术,提出一种计算式(14)的方法,如图3所示。该方法在不采用插值运算的条件下也能显著减小由量化误差引起的积累损失。
图3 双基地RFT变换频域算法流程Fig. 3 Process of bistatic RFT algorithm in frequency domain
步骤1 对接收到的M个基带回波脉冲进行FFT处理;
步骤2 产生距离匹配滤波器的频率响应函数P*(f)并与基带回波脉冲相乘;
步骤3 在频域乘以引起距离走动的指数项的共轭exp(j2πfnvmTr/c),将斜率为v的直线运动轨迹校正为垂直于t轴的直线;
步骤4 通过FFT变换对脉间相位差进行补偿并沿方位维对每一距离频率单元进行累加;
步骤5 对步骤4的结果进行IFFT处理可得G(r, v),峰值检测得到目标结果。
该算法基于FFT和IFFT实现,运算复杂度较低。设距离域采样点数为Nr,速度域采样点数为Nv,相参积累脉冲数为M。
按图3频域算法需要复乘次数分别为:步骤1: (1/2)MNrlog2(Nr),步骤2:MNr,步骤3、步骤4:[2(M+Nv)+4(M+Nv)log2(M+Nv)]Nr,步骤5:(1/2)[MNrlog2(Nr)+NrNv]。
3.3 运动参数空间搜索范围
距离搜索范围。双基地雷达工作于侧向区,由参考文献[10]可知,此时双基地角范围
图5给出了不同SNR条件下的RFT算法检测性能曲线(Pd1),仿真中虚警概率为Pfa=10–6。作为比较,同时给出了MTD相参积累时的检测性能曲线(Pd2)、跨距离单元非相参积累时的检测性能曲线(Pd3)和单个脉冲检测性能曲线(Pd4)。可以看出,Pd=0.9时,本文方法所需SNR最低,为–38.7 dB,(30°, 135°),基线长度范围为单基地雷达最大作用距离。由双基地雷达三角关系,最大双基地搜索距离为:
速度搜索范围。由起始时刻目标双基地径向速度可知,对所有的β值,当时,多普勒频移为正;当时,多普勒频移为负。双基地搜索速度最大值为:
其中,vTmax为观测目标的最大速度。
雷达载频fc=3 GHz,脉冲重复频率fr=1500 Hz,最大作用距离RM=100 km,发射LFM脉冲时宽Tp=10 μs,带宽B=30 MHz,采样频率fs=30 MHz,独立接收机有和单基地雷达接收机相同的参数。如图1所示,设发射站位置(–50 km, 0),接收站位置(50 km, 0),观测区域内有一目标位于(20 km, 50 km)处,其双基地后向散射系数σ0=1,以恒速vT=200 m/s运动,目标方向角δ=45°。忽略波束影响,此时初始双基地距离r0=44.33 km,初始双基地速度v0=–208 m/s,最大可相参积累脉冲数为573个,积累时间内目标运动跨7个搜索步长。
图4 双基地RFT变换仿真Fig. 4 Simulation of bistatic RFT
图5 双基地RFT算法检测性能Fig. 5 Detection performance of bistatic RFT
本文提出了基于RFT变换的双基地雷达相参积累方法,该算法能够在频域实现,且运算复杂度较低。在目标相参积累角内,增加接收机数量,雷达系统能够获得更多的积累脉冲数,从而为弱目标检测提供一条有益的途径。
[1]许稼, 彭应宁, 夏香根, 等. 空时频检测前聚焦雷达信号处理方法[J]. 雷达学报, 2014, 3(2): 129–141. Xu Jia, Peng Ying-ning, Xia Xiang-gen, et al.. Radar signal processing method of space-time-frequency focus-beforedetects[J]. Journal of Radars, 2014, 3(2): 129–141.
[2]田静. 雷达机动目标长时间积累信号处理算法研究[D]. [博士论文], 北京理工大学, 2014: 9–20. Tian Jing. Study on long-time coherent integration for radar maneuvering targets[D]. [Ph.D. dissertation], Beijing Institute of Technology, 2014: 9–20.
[3]Zeng Jian-kui, He Zi-shu, Mathini S, et al.. Modified Hough transform for searching radar detection[J]. IEEE Geoscience and Remote Sensing Letters, 2008, 5(4): 683–686.
[4]Carretero-Moya J, Gismero-Menoyo J, Asensio-Lopez A, et al.. A coherent radon transform for small target detection[C]. IEEE Radar Conference, California, USA, 2009: 1–4.
[5]战立晓, 刘润华, 汤子跃, 等. 基于相参Radon变换的米波相控阵雷达弱目标检测[J]. 现代雷达, 2013, 35(9): 30–35. Zhan Li-xiao, Liu Run-hua, Tang Zi-yue, et al.. Target detection for metric-band phased array radar based on coherent radon transform[J]. Modern Radar, 2013, 35(9): 30–35.
[6]Xu J, Yu J, Peng Y, 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): 1186–1202.
[7]Xu J, Yu J, Peng Y, et al.. Radon-Fourier transform for radar target detection II: Blind speed sidelobe suppression[J]. IEEE Transactions on Aerospace and Electronic Systems, 2011, 47(4): 2473–2489.
[8]Yu J, Xu J, Peng Y, et al.. Radon-Fourier transform for radar target detection III: Optimality and fast implementations[J]. IEEE Transactions on Aerospace and Electronic Systems, 2012, 48(2): 991–1004.
[9]Baker C J, Lhume A, and Qineti Q. Netted radar sensing[J]. IEEE Aerospace & Electronic Systems Magazine, 2003, 18(2): 3–6.
[10]杨振起, 张永顺, 骆永军. 双(多)基地雷达系统[M]. 北京: 国防工业出版社, 1998: 20–22. Yang Zhen-qi, Zhang Yong-shun, and Luo Yong-jun. Bistatic (Multistatic) Radar Systems[M]. Beijing: National Defence Industry Press, 1998: 20–22.
[11]Guan Xin, Zhong Li-hua, Hu Dong-hui, et al.. An extended processing scheme for coherent integration and parameter estimation based on matched filtering in passive radar[J]. Journal of Zhejiang University-Science C (Computers & Electronics), 2014, 15(11): 1071–1085.
[12]Richards M A. 雷达信号处理基础[M]. 北京: 电子工业出版社, 2008: 229–248. Richards M A. Fundamentals of Radar Signal Processing[M]. Beijing: Publishing House of Electronics Industry, 2008: 229–248.
[13]Kay S M. 统计信号处理基础[M]. 北京: 电子工业出版社, 2014: 448–452. Kay S M. Fundamentals of Statical Signal Processing[M]. Beijing: Publishing House of Electronics Industry, 2014: 448–452.
Target Integration and Detection with the Radon-Fourier Transform for Bistatic Radar
Lin Chunfeng Huang Chunlin Su Yi
(School of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China)
Increasing the coherent integration time of the signal can improve the detection of weak targets. The key to the long-time coherent integration for bistatic radar is to solve the target motion compensation from echoes with nonlinear phases. Neglecting the effect of the crossing beam, the compensation of the Radon-Fourier Transform (RFT) was studied in the range and velocity fields. In addition, the method can be operated in the frequency domain. The range walk correction and pulse accumulation were completed by jointly searching along the bistatic range and velocity parameter space. The proposed method was verified by theoretical analysis and simulations.
Bistatic radar; Coherent integration; Radon-Fourier Transform (RFT)
TN957
A
2095-283X(2016)05-0526-05
10.12000/JR16049
林春风, 黄春琳, 粟毅. 双基地雷达Radon-Fourier变换弱目标积累检测[J]. 雷达学报, 2016, 5(5): 526–530.
10.12000/JR16049.
Reference format: Lin Chunfeng, Huang Chunlin, and Su Yi. Target integration and detection with the Radon-Fourier transform for bistatic radar[J]. Journal of Radars, 2016, 5(5): 526–530. DOI: 10.12000/JR16049.
林春风(1984–),男,浙江绍兴人,现为国防科技大学电子科学与工程学院硕士研究生,研究方向为雷达信号处理。
E-mail: 297041767@qq.com
黄春琳(1973–),男,江西赣州人,现为国防科技大学电子科学与工程学院副教授,研究方向为表层穿透雷达系统及应用、雷达与通信信号处理、遥感信息处理等。
E-mail: hclg@163.com
粟 毅(1961–),男,现为国防科技大学电子科学与工程学院教授,研究方向为信号处理、雷达系统、遥感信息处理等。
E-mail: yi.su@yeah.net
2016-03-02;改回日期:2016-05-10;
2016-05-31
*通信作者:林春风 297041767@qq.com
国家自然科学基金(61372160),国防科技大学创新基金(2100040316001)
Foundation Items: The National Natural Science Foundation of China (61372160); The Innovation Foundation of National University of Defense Technology (2100040316001)