一种目标距离走动快速补偿方法技术

针对现有距离徙动校正算法计算量大、需要先验信息、低信噪比情况下估计性能较差等问题，本发明专利技术提出了一种目标距离走动快速补偿方法，该方法充分挖掘目标在距离

【技术实现步骤摘要】
一种目标距离走动快速补偿方法


[0001]本专利技术涉及雷达
，是一种长时间目标能量积累下的目标距离走动快速补偿方法。本专利技术可用于雷达长时间观测内的距离走动量稳健估计及补偿，实现对跨距离单元走动目标能量的有效积累，提升雷达目标检测性能。

技术介绍

[0002]雷达作为一种主动微波目标探测设备，具有全天时、全天候和作用距离远等优势，在目标监视以及预警探测等公共和国防安全领域应用广泛。然而，由于环境的复杂性以及目标的机动性，导致目标能量在距离和多普勒维出现扩散的现象，使得目标回波信号被淹没在强杂波或强噪声中，以致雷达对动目标的探测性能难以满足实际需求。
[0003]长时间积累是实现低信噪比目标检测的有效方法，通过积累时间的增加可增强接收目标的回波能量，即用时间代价换取能量增益。针对目标的运动模型不同，其采用的积累时间也不同。长时间积累检测可分为非相参积累检测、相参积累检测和混合积累检测三种。随着积累时间的增加，目标回波可能出现距离徙动(包括距离走动、距离弯曲)等现象，因此进行有效的包络补偿是实现长时间相参积累的关键。常用的包络补偿方法有：相邻相关法、频域法、Radon变换、一阶/二阶Keystone法等。相邻相关法是通过计算相邻距离单元与参考距离单元之间的相关函数，并通过搜索最大相关值，估计出目标的速度信息，进而根据速度估计值进行回波包络补偿。该方法要求较高的信噪比(Signal to Noise Ratio,SNR)，在信噪比较低时速度估计能力下降，导致检测微弱目标性能下降。频域法是把快时间
‑
慢时间二维回波信号变换到快时间频率
‑
慢时间域，再通过构造频率域距离走动补偿函数对包络距离徙动进行补偿，将回波信号集中到同一距离单元。频域补偿法虽操作简单，但在运动参数未知的情况下，需要对目标的运动速度进行大量的搜索操作。Radon变换是在距离
‑
慢时间二维平面内，利用Radon变换方法将匀速运动目标的轨迹映射至Radon域内的一个点上，通过在Radon域内进行最大值搜索，可对目标运动速度进行估计，从而根据估计的速度行进包络补偿。但该方法在进行参数估计时，需要对Radon的旋转角度进行搜索，搜索过程将会极大提高雷达系统的运算复杂度。Keystone变换可以在目标径向速度未知情况下，通过变量代换实现对慢时间轴的伸缩变换，通过解除快时间和慢时间的线性耦合消除包络的距离徙动。Keystone变换无需进行搜索操作，但一阶和二阶Keystone变换都要进行大量的插值运算，影响了对弱目标的实时检测性能。

技术实现思路

[0004]要解决的技术问题
[0005]针对现有距离徙动校正算法计算量大、需要先验信息、低信噪比情况下估计性能较差等问题，本专利技术提供了一种目标距离走动快速补偿方法。
[0006]技术方案
[0007]一种目标距离走动快速补偿方法，其特征在于步骤如下：
[0008]步骤1：构造原始回波矩阵
[0009]雷达接收的目标回波基带信号可表示为：
[0010][0011]其中，η分别表示快时间和慢时间，A表示幅度系数，f
c
表示载波频率；
[0012]假设目标相对于雷达做匀速运动，τ
n
(η)表示目标时延：
[0013][0014]其中，R0为目标初始距离，v为目标径向速度，c为光速，1＜n≤N，N为方位维脉冲个数；
[0015]经过距离脉压后，其距离
‑
慢时间二维回波数据可表示为：
[0016][0017]其中，A
p
表示脉压后的幅度，ρ(
·
)表示参考信号与回波信号的互相关函数；
[0018]将式(3)沿快时间维进行傅里叶变换，可得距离频域
‑
慢时间二维回波数据：
[0019][0020]其中，P(
·
)为ρ(
·
)的傅里叶变换结果；
[0021]参考信号与目标回波进行互相关处理后得到所有脉冲内的距离脉冲压缩信号离散化后的数据记为D
mn
，其中1＜m≤M,1＜n≤N，M为距离维傅立叶变换后的点数，N为方位维脉冲个数；
[0022]步骤2：阈值处理
[0023]门限因子为T
k
，回波矩阵D
mn
的能量最大值为E
max
，在回波矩阵D
mn
中选取出幅度值大于T
k
×
E
max
的所有点，记为(d1,d2,
…
,d
K
)，其中K<M，K为预处理后选取的点数；(x
i
,y
i
)为点d
i
的坐标，再利用选取出来的K个点，构造新的矩阵，计作P
kn
；
[0024]步骤3：均值处理
[0025]新矩阵P
kn
的平均向量A(K
x
,N
y
)定义为：
[0026][0027]其中，K
x
、N
y
分别表示距离和方位坐标平均值；
[0028]最终数据库矩阵是通过从矩阵P
kn
中减去平均向量A得到：
[0029][0030]步骤4：矩阵分解
[0031]对矩阵利用Godec快速RPCA算法进行矩阵分解，得到包含目标信息的低秩矩阵L和包含噪声的稀疏矩阵S：
[0032][0033]步骤5：求斜率
[0034]RPCA分解出的包含目标信息的低秩矩阵L是一组秩为1的二维矩阵，其大小为(2,K)；因此L矩阵的轨迹斜率为：
[0035][0036]其中，(x1,y1),(x2,y2)为L矩阵中任意两个点；
[0037]估计出的目标速度参数为：
[0038][0039]其中，PRF是脉冲重复频率，f
s
是采样频率；
[0040]步骤6：构造线性补偿函数
[0041][0042]式(4)与补偿函数相乘后的表达式为：
[0043][0044]可以看出，由距离维频域f和速度v引起的耦合已经被完全消除，因此线性距离走动被完全校正；
[0045]沿着距离频域维f进行IFT，得到距离
‑
方位的信号表达式为：
[0046][0047]一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的方法。
[0048]一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现上述的方法。
[0049]有益效果
[0050]本专利技术提出的一种目标距离走动快速补偿方法，该方法充分挖掘目标在距离...

【技术保护点】

【技术特征摘要】
1.一种目标距离走动快速补偿方法，其特征在于步骤如下：步骤1：构造原始回波矩阵雷达接收的目标回波基带信号可表示为：其中，η分别表示快时间和慢时间，A表示幅度系数，f
c
表示载波频率；假设目标相对于雷达做匀速运动，τ
n
(η)表示目标时延：其中，R0为目标初始距离，v为目标径向速度，c为光速，1＜n≤N，N为方位维脉冲个数；经过距离脉压后，其距离
‑
慢时间二维回波数据可表示为：其中，A
p
表示脉压后的幅度，ρ(
·
)表示参考信号与回波信号的互相关函数；将式(3)沿快时间维进行傅里叶变换，可得距离频域
‑
慢时间二维回波数据：其中，P(
·
)为ρ(
·
)的傅里叶变换结果；参考信号与目标回波进行互相关处理后得到所有脉冲内的距离脉冲压缩信号离散化后的数据记为D
mn
，其中1＜m≤M,1＜n≤N，M为距离维傅立叶变换后的点数，N为方位维脉冲个数；步骤2：阈值处理门限因子为T
k
，回波矩阵D
mn
的能量最大值为E
max
，在回波矩阵D
mn
中选取出幅度值大于T
k
×
E
max
的所有点，记为(d1,d2,
…
,d
K
)，其中K<M，K为预处理后选取的点数；(x
i
...

【专利技术属性】
技术研发人员：粟嘉，崔国楠，陶明亮，范一飞，李滔，王伶，张兆林，韩闯，宫延云，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：