一种基于频域切变的相参积累方法
【专利摘要】本发明涉及一种基于频域切变的相参积累方法，当目标回波包络在相参处理间隔内出现跨距离单元走动时，将距离-慢时间维信号空间的包络走动轨迹看作图像，再通过频域切变处理方法进行一定角度的旋转平移，每一次旋转平移过程通过慢时间维快速傅里叶变换得到相参积累峰值；这样每一个旋转角度对应一个积累峰值，以最大峰值位置得到对应的目标轨迹的斜率、及起始距离，从而得到目标的运动轨迹。本发明有效实现距离走动校正和目标能量的积累，适用于雷达微弱目标信号检测、及雷达成像技术中距离走动校正，其精度高、运算速度快、易于工程实现、鲁棒性好。
【专利说明】一种基于频域切变的相参积累方法
【技术领域】
[0001] 本发明属于雷达系统中目标检测和跟踪【技术领域】，特别涉及一种基于频域切变的 相参积累方法。
【背景技术】
[0002] 隐身目标、无人机以及巡航导弹等微弱目标，对传统雷达系统构成严重的威胁，而 长时间相参积累技术能有效的对隐身等微弱目标进行检测。但在雷达信号长时间相参积 累处理中，目标运动引起的跨距离走动严重影响目标能量的积累和检测。并且伴随着雷达 对高分辨率的需求，距离单元越来越小，跨距离单元现象更加严重。对于跨距离单元补偿问 题，需要目标运动参数的精确估计，而一般情况下低信噪比（SNR)时，目标运动参数的精确 估计是难以做到的。因此未知目标运动参数情况下，如何进行积累检测和改善积累检测的 性能是长时间信号积累检测要研究的主要问题。
[0003] 在成像系统中，逆合成孔径雷达（ISAR)的平动补偿包络对齐方法，如互相关法、 模-1或模-2法、最小熵法等，在回波SNR较高时才有较好的包络对齐效果。对目标运动 参数搜索类的方法，其时域移位或寻址量为整数，在未插值时有能量积累损失，并且运算量 大。Keystone变换能够在一定速度范围内的跨距离单元进行统一补偿（详见"张顺生，曾 涛，基于Keystone变换的微弱目标检测，电子学报，2005, 33 (6) : 1675-1678. "），但是该 算法需要对多普勒模糊数进行估计，否则性能严重下降。图像处理中，直线检测的Hough 变换和Radon变换应用于"检测前跟踪（TBD) "技术中，但是由于非相参积累存在SNR阈值 现象，对低SNR积累效率不高。相参Radon变换对微弱目标进行检测，其采用旋转和双线 性插值实现Radon变换并通过多普勒滤波器进行相位补偿实现相参积累（详见"JavierC M,JavierGM,etal. ,ACoherentRadonTransformforSmallTargetDetection,IEEE RaderConference, 2009:1-4. "），但是插值存在损失且容易导致旁瓣的产生。
【发明内容】
[0004] 本发明的目的是提供一种适用于隐身等微弱目标回波信号，采用长时间相参积累 技术过程出现的跨距离单元走动问题的解决方法。按照本发明提供的方法对微弱目标回波 信号进行距离走动校正和相参积累，得到低信噪比下目标回波信号的能量积累，同时得到 运动目标的轨迹信息。
[0005] 为了实现上述目的，本发明提供一种基于频域切变的相参积累方法，其中包含如 下步骤：
[0006]S1、雷达对接收信号处理后得到基带信号，进一步得到包含目标的包络轨迹信号 距离-慢时间维信号空间；
[0007]S2、根据目标的速度范围确定所要进行频域切变的角度范围，并确定角度搜索过 程中的角度步进；
[0008]S3、对每一个搜索角度进行频域切变，来进行包络轨迹的搜索；
[0009]S4、沿慢时间维快速傅里叶变换完成一次角度搜索结果的多普勒滤波，得到快速 傅里叶变换之后该角度对应的距离-多普勒参数空间；
[0010]S5、判断距离-多普勒参数空间峰值是否为最大峰值，如果是最大峰值，则返回最 大峰值位置和对应的旋转角度；如果不是最大峰值，重复S3?S5,直至搜索角度遍历完毕， 输出最大峰值位置及其对应的旋转角度；
[0011]S6、对返回的最大峰值所在的距离-多普勒参数空间进行恒虚警检测，如果过门 限，则输出此时参数空间的距离信息，并根据此时的角度值得到目标的速度信息；
[0012] 通过S1?S6,完成距离走动情况下的距离走动补偿，并得到相参积累结果，从而 得到目标的运动轨迹。
[0013] 优选示例的一种基于频域切变的相参积累方法中，包含如下步骤：
[0014] (1)雷达对接收信号处理后得到基带信号，进一步得到包含目标的包络轨迹信号 距离-慢时间维信号空间。
[0015] (2)根据目标的速度范围确定所要进行频域切变的角度范围，并确定角度搜索过 程中的角度步进A0。
[0016] (3)对包络轨迹s(t，〇进行频域切变进行包络轨迹的搜索。首先，初始化切变旋 转角度。
[0017] (4)整个切变过程始终围绕信号空间的中心进行旋转平移。为了表示方便，将 f(x，y)表示包络轨迹s(t，tm)，u表示距离维x的频域，v表示慢时间维y的频域。则行切 变过程表示为
[0018]Mx (x,y) =IFFTU{Q: (u, y)FFTX{f(x,y)}}
[0019] 其中，FFTX表示对x进行快速傅里叶变换；IFFTU表示对u进行快速傅里叶逆变换； Qjuj)为
【权利要求】
1. 一种基于频域切变的相参积累方法，其特征在于，包含如下过程： 51、 雷达对接收信号处理后得到基带信号，进一步得到包含目标的包络轨迹信号距 离-慢时间维信号空间； 52、 根据目标的速度范围确定所要进行频域切变的角度范围，并确定角度搜索过程中 的角度步进； 53、 对每一个搜索角度进行频域切变，来进行包络轨迹的搜索； 54、 沿慢时间维快速傅里叶变换完成一次角度搜索结果的多普勒滤波，得到快速傅里 叶变换之后该角度对应的距离-多普勒参数空间； 55、 判断距离-多普勒参数空间峰值是否为最大峰值，如果是最大峰值，则返回最大峰 值位置和对应的旋转角度；如果不是最大峰值，重复S3?S5,直至搜索角度遍历完毕，输出 最大峰值位置及其对应的旋转角度； 56、 对返回的最大峰值所在的距离-多普勒参数空间进行恒虚警检测，如果过门限，则 输出此时参数空间的距离信息，并根据此时的角度值得到目标的速度信息； 通过Sl?S6,完成距离走动情况下的距离走动补偿，并得到相参积累结果，从而得到 目标的运动轨迹。
2. 如权利要求1所述的相参积累方法，其特征在于， 所述步骤1中，设远场、单目标回波信号在一个相参积累间隔内雷达接收到的第m个回 波基带信号为 sr (t, tm) = Arp (t-2R (tm) /c) exp (-j4 n fcR (tm) /c) 其中，t和^分别为快时间和慢时间，^ = ml；，且I；为脉冲重复间隔； f。为载波频率；c为光速;4为点目标的回波幅度；p (t)表示发射波形基带信号；R(tm) =Rc^vtm其中Rtl和R(tJ分别为h和tm时刻目标相对于雷达的径向距离。
3. 如权利要求2所述的相参积累方法，其特征在于， 如果发射波形基带信号P (t)为方波，则回波的包络轨迹表示为 s (t，tm) Sr(t，tm); 如果发射波形基带信号为线性调频信号，则包络轨迹表示为 s (t, tm) = A1Bsinc [B (t~2R (tm) /c) ] exp (-j4 n fcvtm/c) 其中，A1 = Arexp(-j4 Ji fcRQ/c) ;B 为信号带宽。
4. 如权利要求3所述的相参积累方法，其特征在于， 所述步骤2中，目标的速度范围为Vmin < V < Vmax时，所需要搜索的角度范围为
其中，AR为距离分辨率。
5. 如权利要求4所述的相参积累方法，其特征在于， 所述步骤2中，角度搜索过程中的角度步进A 0 :
其中，P为当前速度，速度搜索步长为Av= ARAMI；)。
6.如权利要求1或5所述的相参积累方法，其特征在于， 所述S3中进一步包含如下过程： S3-1、对包络轨迹s (t，tm)进行频域切变进行包络轨迹的搜索；首先，初始化切变旋转 角度； S3-2、整个切变过程始终围绕信号空间的中心进行旋转平移；将f(x，y)表示包络轨迹 s (t，〇，u表示距离维X的频域，V表示慢时间维y的频域；则进行行切变的过程表示为 Mx (x, y) = IFFTu (Q1 (u, y) FFTx {f (x, y)}} 其中，FFTx表示对x进行快速傅里叶变换；IFFTu表示对u进行快速傅里叶逆变换； Qju, y)由下式得到：
其中，9为切变旋转角度；M为相参处理间隔内的脉冲数；N为距离维采样点数； S3-3、对第一次行切变之后的信号空间进行列切变变换；即对Mx(x，y)进行列方向的变 换： Myx (x, y) = IFFTv (Q2 (x, v) FFTy {Mx (x, y)}} 其中，FFTy表示对y进行快速傅里叶变换；IFFTv表示对v进行快速傅里叶逆变换； Q2 (X，V)由下式得到：
S3-4、对Myx(x，y)进行行方向的变换： Mxyx (X，y) = IFFTu (Q1 (u，y) FFTx {Myx (X，y)}} 所述S3-1?S3-6完成一次角度切变旋转变换。
【文档编号】G01S13/64GK104330791SQ201410577050
【公开日】2015年2月4日 申请日期:2014年10月24日 优先权日:2014年10月24日
【发明者】刘俊豪, 陈潜, 付朝伟, 王海涛 申请人:上海无线电设备研究所