一种基于二维局部插值的雷达目标距离联合估计方法技术

本发明专利技术涉及一种基于二维局部插值的雷达目标距离联合估计方法，包括以下步骤：在雷达目标的距离‑多普勒图像上进行检索，得到峰值点的位置、沿距离轴的次强点位置以及沿多普勒频率轴的次强点位置；距离‑多普勒图像的峰值局部进行二维插值运算，得到图像峰值位置与距离轴的偏差值和与多普勒轴的偏差值，最终确定目标的实际距离‑多普勒位置。因此本发明专利技术中通过搜索距离‑多普勒的峰值点和仅次于峰值点的次强点的位置坐标，对峰值点位置进行局部插值，提高了对目标多普勒频率的估计精度以及对目标的距离位置的估算精度，减少了数据采集的工作量，对峰值点位置进行局部插值与对全局进行内插计算相比，计算量小且精度高。

A Joint Range Estimation Method for Radar Targets Based on Two-Dimensional Local Interpolation

The present invention relates to a joint range estimation method for radar target based on two-dimensional local interpolation, which includes the following steps: retrieving the range Doppler image of radar target, obtaining the position of the peak point, the position of the sub-strong point along the range axis and the position of the sub-strong point along the Doppler frequency axis; and the range Doppler image. Two-dimensional interpolation is performed to get the deviation between the peak position of the image and the distance axis and the deviation from the Doppler axis. Finally, the actual range and Doppler position of the target are determined. Therefore, by searching the position coordinates of the peak point of range Doppler and the sub-strong point next to the peak point, the position of the peak point is interpolated locally, which improves the estimation accuracy of the target Doppler frequency and the distance position of the target, reduces the workload of data acquisition and the peak point position. Compared with global interpolation, local interpolation is less computational and more accurate.

【技术实现步骤摘要】
一种基于二维局部插值的雷达目标距离联合估计方法
本专利技术涉及雷达目标探测
，尤其涉及一种基于二维局部插值的雷达目标距离联合估计方法。
技术介绍
雷达具有全天候、全天时工作能力，能够在雨、雪、雾等恶劣天气以及黑夜等特殊环境下持续探测感兴趣的运动目标，对运动目标的位置信息和多普勒频率进行估计。目前广泛使用的运动目标检测方法是距离-多普勒处理，利用该方法可以同时获取目标的距离和多普勒信息，具有性能稳健，运算量低的特点。但是，传统的距离-多普勒处理方法只是对相干处理时间内的回波信号进行二维离散傅里叶变换(DFT)，然后依据二维DFT变换后的距离-多普勒图像的峰值点坐标换算得到目标的距离和多普勒频率，参数估计的精度受限于二维坐标轴的数据采样点数，且只有在目标的实际运动参数为二维信号处理精度的整数倍时，才能对目标距离、多普勒信息进行精确估计。只依靠更为密集地数据采样来提高二维DFT的估计精度，将极大地增加二维数据矩阵的整体内插运算量，计算费时，且不能实现无偏估计。
技术实现思路
本专利技术所要解决的技术问题是提供一种数据采样点和运算量少的基于二维局部插值的雷达目标距离联合估计方法。为解决上述技术问题，本专利技术的技术方案是：一种基于二维局部插值的雷达目标距离联合估计方法，包括以下步骤：在雷达目标的距离-多普勒图像上进行峰值检索，得到峰值点的位置坐标X(R0，D0)，再分别以峰值点所在位置为原点沿距离和多普勒坐标轴的正负方向进行搜索，确定沿距离轴的仅次于峰值点的次强点位置坐标X(R0+i，D0)，i＝1或-1，以及沿多普勒频率轴的仅次于峰值点的次强点位置坐标X(R0，D0+j)，j＝1或-1；通过所述峰值点位置坐标X(R0，D0)、仅次于峰值点的次强点位置坐标X(R0+i，D0)和X(R0，D0+j)，对所述距离-多普勒图像的峰值局部进行二维插值运算，得到所述图像峰值位置与距离轴的偏差值δR和与多普勒轴的偏差值δD，最终确定目标的实际距离-多普勒位置X(R0+δR，D0+δD)。作为优选的技术方案，所述雷达目标的距离-多普勒图像是通过将接收到的雷达回波进行脉冲压缩和杂波抑制，得到回波的距离-慢时间矩阵，然后，对各距离单元内的脉冲压缩信号沿慢时间轴进行离散傅里叶变换而得到的。作为优选的技术方案，沿距离轴的仅次于峰值点的次强点位置坐标和沿多普勒频率轴的仅次于峰值点的次强点位置坐标的搜索是依据所述雷达回波的理想模糊函数的相应轴线方向搜索。由于采用了上述技术方案，本专利技术的有益效果是：传统的距离-多普勒处理方法，为保证计算精度，需要密集地对多普勒频率轴进行采样，力图在最大程度上减小样本频率错过信号真实多普勒频率的可能性，但这种方法需要对整个二维数据矩阵进行内插运算，计算量大，计算时间长。由于目标的真实多普勒频率通常在DFT的峰值点和次强点之间，因此本专利技术中通过搜索距离-多普勒的峰值点和仅次于峰值点的次强点的位置坐标，对峰值点位置进行局部插值，提高了对目标的多普勒频率的估计精度以及对目标的距离位置的估算精度，减少了数据采集的工作量，对峰值点位置进行局部插值与对全局进行内插计算相比，计算量小且精度高。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术实施例的流程图；图2是本专利技术实施例中雷达回波的距离-慢时间图像；图3是本专利技术实施例中目标距离-多普勒的亮度图；图4是本专利技术实施例中目标距离-多普勒的局部网线图；图5是本专利技术实施例的单脉冲模糊函数峰值周边次强点的取值示意图；图6是本专利技术实施例的线性调频信号模糊函数峰值周边次强点的取值示意图；图7是本专利技术实施例中目标峰值多普勒切面下，理想采样与实际采样的对比示意图。具体实施方式一种基于二维局部插值的雷达目标距离联合估计方法，包括以下步骤：在雷达目标的距离-多普勒图像上进行峰值检索，得到峰值点的位置坐标X(R0，D0)，再分别以峰值点所在位置为原点沿距离和多普勒坐标轴的正负方向进行搜索，确定沿距离轴的仅次于峰值点的次强点位置坐标X(R0+i，D0)，i＝1或-1，以及沿多普勒频率轴的仅次于峰值点的次强点位置坐标X(R0，D0+j)，j＝1或-1；通过峰值点位置坐标X(R0，D0)、仅次于峰值点的次强点位置坐标X(R0+i，D0)和X(R0，D0+j)，对距离-多普勒图像的峰值局部进行二维插值运算，得到图像峰值位置与距离轴的偏差值δR和与多普勒轴的偏差值δD，最终确定目标的实际距离-多普勒位置X(R0+δR，D0+δD)。优选的，雷达目标的距离-多普勒图像是通过将接收到的雷达回波进行脉冲压缩和杂波抑制，得到回波的距离-慢时间矩阵，然后，对各距离单元内的脉冲压缩信号沿慢时间轴进行离散傅里叶变换而得到的。优选的，沿距离轴的仅次于峰值点的次强点位置坐标和沿多普勒频率轴的仅次于峰值点的次强点位置坐标的搜索是依据雷达回波的理想模糊函数的相应轴线方向搜索。如图1至图7共同所示，雷达发射电磁波利用来自运动目标的电磁反射回波对目标进行探测。接收到的回波中包含有目标的距离和运动速度信息，可通过距离-多普勒处理进行提取，并通过局部插值运算提高距离和多普勒频率的估计精度。步骤1：雷达回波预处理。雷达系统发射载频fc＝2.3GHz，雷达带宽为B＝10MHz，脉冲宽度为Tp＝2μs，脉冲重复频率为PRF＝10kHz的雷达信号。目标实际位置为Rref＝2.530km，运动速度为v＝18.5m/s，相干处理时间为100个脉冲，即CPI＝0.10s。接收信号经过解调处理，信号的数学模型为：其中，t为快时间，tm为慢时间，Ar为回波信号的幅度，c为光速，R(tm)为目标随时间的实时距离，λ为信号的波长。回波信号经过脉冲压缩处理：其中，Arm＝TpBAr，R(tm)＝Rref-vtm，tm＝m/PRF将脉冲压缩后的信号表示为离散形式，令快时间的采样间隔为Tft，则t＝2Rref/c+l·Tft，l＝0，...，L-1回波信号的离散形式为：预处理后的回波信号为距离-慢时间坐标下的离散信号矩阵，信号图像如图2所示。步骤2：获得目标的距离-多普勒图像。在一个相干处理间隔(CoherentProcessingInterval，CPI)内，目标近似为匀速运动，对同一相干处理时间内的距离-慢时间信号，沿慢时间轴分别对各距离单元的信号进行离散傅里叶变换，将回波沿慢时间维变换得到多普勒频谱信号，距离-多普勒图像如图3和图4所示。步骤3：确定距离-多普勒图像的峰值位置。对步骤2中得到的距离-多普勒图像进行峰值检索，得到图像峰值点位置。再分别以峰值所在位置为原点沿距离和多普勒坐标轴的正负方向进行搜索，确定仅次于峰值点的次强点的坐标位置。次强点沿距离、多普勒维度进行搜索是由于，对于理想模糊函数，多普勒和距离信息相互独立，故可以沿各自坐标轴分别进行处理。当考虑雷达波形的模糊函数时，次强点的搜索方向需根据雷达波形的模糊函数沿相应的轴线进行搜索。图5和图6分别为两类波形对应的搜索方式，图5所示出的为单脉冲的模糊函数等值线，该本文档来自技高网...

【技术保护点】
1.一种基于二维局部插值的雷达目标距离联合估计方法，其特征在于，包括以下步骤：在雷达目标的距离‑多普勒图像上进行峰值检索，得到峰值点的位置坐标X(R0,D0)，再分别以峰值点所在位置为原点沿距离和多普勒坐标轴的正负方向进行搜索，确定沿距离轴的仅次于峰值点的次强点位置坐标X(R0+i,D0),i＝1或‑1，以及沿多普勒频率轴的仅次于峰值点的次强点位置坐标X(R0,D0+j)，j＝1或‑1；通过所述峰值点位置坐标X(R0,D0)、仅次于峰值点的次强点位置坐标X(R0+i,D0)和X(R0,D0+j),对所述距离‑多普勒图像的峰值局部进行二维插值运算，得到所述图像峰值位置与距离轴的偏差值δR和与多普勒轴的偏差值δD，最终确定目标的实际距离‑多普勒位置X(R0+δR,D0+δD)。

【技术特征摘要】
1.一种基于二维局部插值的雷达目标距离联合估计方法，其特征在于，包括以下步骤：在雷达目标的距离-多普勒图像上进行峰值检索，得到峰值点的位置坐标X(R0,D0)，再分别以峰值点所在位置为原点沿距离和多普勒坐标轴的正负方向进行搜索，确定沿距离轴的仅次于峰值点的次强点位置坐标X(R0+i,D0),i＝1或-1，以及沿多普勒频率轴的仅次于峰值点的次强点位置坐标X(R0,D0+j)，j＝1或-1；通过所述峰值点位置坐标X(R0,D0)、仅次于峰值点的次强点位置坐标X(R0+i,D0)和X(R0,D0+j),对所述距离-多普勒图像的峰值局部进行二维插值运算，得到所述图像峰值位置与距离轴的...

【专利技术属性】
技术研发人员：李强，
申请(专利权)人：湖北晧天智能科技有限公司，
类型：发明
国别省市：湖北,42