一种基于Radon变换的运动目标径向速度估计方法技术

一种基于Radon变换的运动目标径向速度估计方法，该方法避免了Radon变换的角度搜索过程，利用少数几次Radon变换即可实现运动目标径向速度估计，将现有基于Radon变换的运动目标径向速度估计方法由二维搜索降低为数次一维搜索，提高了参数估计的实时性；由于不存在Radon变换角度搜索过程，参数估计的精度不受角度搜索步长的影响，可以使得估计精度得到进一步提升；该方法通过最小二乘方法估计参数，使得本发明专利技术具有较好的稳健性。采用本发明专利技术处理方式运算量较小，估计精度更高，而且在较低SNR条件下即可实现较高的估计精度。因此，本发明专利技术可以满足高实时性和稳健性要求。

【技术实现步骤摘要】

本专利技术提供了一种基于Radon变换的运动目标径向速度估计方法，属于雷达信号处理
，用于雷达目标检测、参数估计、运动目标成像和定位等。
技术介绍
在合成孔径雷达(syntheticapertureradar,SAR)-运动目标检测系统中，运动目标径向速度无模糊估计是实现运动目标成像和精确定位的必要前提。运动目标径向速度无模糊估计分为多通道和单通道方法。多通道方法在沿航迹干涉系统构型下，通过多载频或多基线构造多组包含运动目标径向速度的干涉相位，利用不同载频或基线的参差关系解径向速度模糊。但多通道方法面临以下问题：系统实现成本较高，通道间误差和配准误差等非理想因素对估计性能影响较大。传统的单通道方法基于包络信息，利用距离走动与运动目标径向速度线性相关的思想，估计运动目标的径向速度。由于包络信息不受相位缠绕问题影响，所以估计的径向速度是无模糊的。在单通道SAR系统下，利用Radon变换通过二维搜索估计运动目标信号包络的斜率，可以实现运动目标径向速度无模糊估计。但基于包络信息的单通道估计方法存在以下三点问题：(1)由于目标包络信息并不敏感于运动目标参数，所以估计精度较低；(2)通过二维搜索估计参数，运算量较大；(3)二维搜索步长和参数估计精度需要折衷考虑。
技术实现思路
本专利技术的目的在于：克服现有技术的不足，提出一种基于Radon变换的运动目标径向速度快速且稳健的估计方法，以提高系统实现的实时性和稳健性。本专利技术采用的技术方案为：一种基于Radon变换的运动目标径向速度估计方法，步骤如下：(1)对利用SAR接收的运动目标回波信号进行距离向脉压处理，得到运动目标距离脉压域信号其中，为快时间，tm为慢时间；(2)计算运动目标距离走动倾角θ的最大范围[-θm,θm]，该最大范围即为Radon变换角度范围；其中，θm＝atan(2fsrvrm/c·fsa)，fsr为SAR系统采样频率，fsa为SAR系统脉冲重复频率，c为光速，vrm为所述运动目标径向速度vr所属范围[-vrm,vrm]的边界值，其中，vr＝c·tanθ·fsa/2fsr；(3)根据步骤(2)得到的Radon变换角度范围，得到N次Radon变换的变换角度分别为：[α1,α2,α3,…,αN]，其中，α1＝-θm，αN＝θm，[α1,α2,α3,…,αN]为等间隔排列或者为不等间隔排列；(4)对运动目标距离脉压域信号进行变换角度分别为[α1,α2,α3,…,αN]的N次Radon变换，得到N次Radon变换结果Rn，n＝1,2,…,N，Rn对应的归一化Radon变换结果分别为且其中，max(·)表示取最大值操作；(5)将步骤(4)中得到的归一化Radon变换结果分别与预设门限值进行比较，中大于门限值的数值的个数分别记为L1,L2,L3,…,LN，即为N次Radon变换结果的归一化长度；所述预设门限值为0.4-0.6。(6)由N次Radon变换的变换角度[α1,α2,α3,…,αN]构造观测矩阵A由N次Radon变换的归一化长度L1,L2,L3,…,LN构造观测值L(7)利用最小二乘法，估计包含运动目标距离走动倾角θ信息的矩阵其中，[·]-1表示矩阵求逆运算，[·]T表示矩阵转置操作，Lθ为运动目标距离走动轨迹在二维平面中的长度；(8)通过公式计算运动目标径向速度估计值从而完成所述基于Radon变换的运动目标径向速度估计。本专利技术与现有技术相比具有以下优点：1)与多通道估计方法相比，本专利技术利用单通道实现径向速度估计，降低了系统实现成本，并且不受通道间误差和配准误差的影响；2)与基于相位的估计方法相比，本专利技术不存在模糊估计的问题；3)与基于包络的单通道估计方法相比，在运算量方面，避免了Radon变换的角度搜索过程，利用少数几次Radon变换即可实现运动目标径向速度的快速估计，将现有Radon变换方法的二维搜索降低为数次一维搜索，提高了参数估计的实时性；4)由于不存在Radon变换角度搜索过程，运动目标径向速度的估计精度不受角度搜索步长的限制，可以使得估计精度得到进一步提升；5)通过最小二乘方法，使得本专利技术具有较好的稳健性。附图说明图1是本专利技术方法的流程图；图2是本专利技术与现有Radon变换方法在不同Radon变换次数下的对比仿真图；图3为本专利技术与现有Radon变换方法在不同信噪比条件下的对比仿真图；具体实施方式合成孔径雷达-运动目标检测(SAR-GMTI)系统广泛应用于战场侦察、交通管制等应用中。在SAR-GMTI系统中，运动目标的径向速度导致SAR成像处理后运动目标的方位位置发生偏移；运动目标径向速度导致的运动目标距离走动使得运动目标成像结果散焦，因此，估计运动目标径向速度是实现运动目标精确定位和高分辨率成像的必要前提。沿航迹干涉(ATI)系统利用了径向速度与干涉相位的关系，常用于估计运动目标径向速度。在常规沿航迹干涉系统中，由于干涉相位存在相位缠绕问题，当运动目标径向速度超过系统最大不模糊速度时，径向速度估计就会出现模糊。为了在ATI系统中实现径向速度无模糊估计，研究人员通过利用多波长、多基线、多速度SAR等方法实现解模糊。但这些方法需要在多通道SAR系统上利用多波长、多基线，或者利用不同速度的SAR运动平台，构成多组径向速度与干涉相位关系，实现相位解缠绕，从而无模糊的估计运动目标径向速度。但这些方法不仅增加系统成本、灵活性较差，而且估计性能精度受通道间误差影响很大。为避免以上问题，研究人员利用Radon变换估计运动目标径向速度。由于运动目标导致距离走动，该距离走动轨迹在二维平面上呈现为一条斜线，而运动目标距离走动轨迹的倾斜角度θ(简称距离走动倾角)与运动目标径向速度vr具有如下关系vr＝c·tanθ·fsa/2fsr，其中，fsr为SAR系统采样频率，fsa为SAR系统脉冲重复频率，c为光速。由于Radon变换能够估计二维平面上直线的倾角，研究人员利用Radon变换通过搜索所有可能的角度：当Radon变换角度与运动目标轨迹距离走动倾角最为接近时，在Radon变换域的结果最大，此时的Radon变换角度即可用于估计运动目标径向速度，这就是现有Radon变换估计方法的核心思想。但在实际应用中，由于Radon变换涉及到距离和角度的二维搜索问题，其运算量较大；而且搜索步长与径向速度估计精度和运算量都有关系：搜索步长越小，估计径向速度所需的运算量越大，实现的估计精度越高；而搜索步长越大，运算量越小，估计精度越差该方法实现简单，但本文档来自技高网...

【技术保护点】
一种基于Radon变换的运动目标径向速度估计方法，其特征在于步骤如下：(1)对利用SAR接收的运动目标回波信号进行距离向脉压处理，得到运动目标距离脉压域信号其中，为快时间，tm为慢时间；(2)计算运动目标距离走动倾角θ的最大范围[‑θm,θm]，该最大范围即为Radon变换角度范围；(3)根据步骤(2)得到的Radon变换角度范围，得到N次Radon变换的变换角度分别为：[α1,α2,α3,…,αN]，其中，α1＝‑θm，αN＝θm，[α1,α2,α3,…,αN]为等间隔排列或者为不等间隔排列；(4)对运动目标距离脉压域信号进行变换角度分别为[α1,α2,α3,…,αN]的N次Radon变换，得到N次Radon变换结果Rn，n＝1,2,…,N，Rn对应的归一化Radon变换结果分别为(5)将步骤(4)中得到的归一化Radon变换结果分别与预设门限值进行比较，中大于门限值的数值的个数分别记为L1,L2,L3,…,LN，即为N次Radon变换结果的归一化长度；(6)由N次Radon变换的变换角度[α1,α2,α3,…,αN]构造观测矩阵A，由N次Radon变换的归一化长度L1,L2,L3,…,LN构造观测值L；(7)对步骤(6)中得到的观测矩阵A和观测值L，利用最小二乘法，估计包含运动目标距离走动倾角θ信息的矩阵Lθsinθ^Lθcosθ^;]]>(8)根据步骤(7)得到的包含运动目标距离走动倾角θ信息的矩阵Lθsinθ^Lθcosθ^,]]>计算运动目标径向速度估计值从而完成所述基于Radon变换的运动目标径向速度估计。...

【技术特征摘要】
1.一种基于Radon变换的运动目标径向速度估计方法，其特征在于步骤如下：
(1)对利用SAR接收的运动目标回波信号进行距离向脉压处理，得到运动目标距
离脉压域信号其中，为快时间，tm为慢时间；
(2)计算运动目标距离走动倾角θ的最大范围[-θm,θm]，该最大范围即为Radon变换角度
范围；
(3)根据步骤(2)得到的Radon变换角度范围，得到N次Radon变换的变换角度分别为：
[α1,α2,α3,…,αN]，其中，α1＝-θm，αN＝θm，[α1,α2,α3,…,αN]为等间隔排列或者为不等间隔排
列；
(4)对运动目标距离脉压域信号进行变换角度分别为[α1,α2,α3,…,αN]的N次
Radon变换，得到N次Radon变换结果Rn，n＝1,2,…,N，Rn对应的归一化Radon变换结果分别为
(5)将步骤(4)中得到的归一化Radon变换结果分别与预设门限值进行
比较，中大于门限值的数值的个数分别记为L1,L2,L3,…,LN，即为N次Radon
变换结果的归一化长度；
(6)由N次Radon变换的变换角度[α1,α2,α3,…,αN]构造观测矩阵A，由N次Radon变换的归
一化长度L1,L2,L3,…,LN构造观测值L；
(7)对步骤(6)中得到的观测矩阵A和观测值L，利用最小二乘法，估计包含运动目标距
离走动倾角θ信息的矩阵Lθsinθ^Lθcosθ^;]]>(8)根据步骤(7)得到的包含运动目标距离走动倾角θ信息的矩阵Lθsinθ^Lθcosθ^,]]>计算运动
目标径向速度估计值从而完成所述基于Radon变换的运动目标径向速度估计。
2.根据权利要求1所述的一种基于Radon变换的运动目标径向速度估计方法，其特征在
于：步骤(5)中所述预设门限值为0.4-0.6。
3.根据权利要求1所述的一种基于R...

【专利技术属性】
技术研发人员：张学攀，刘波，刘露，刘敏，
申请(专利权)人：中国空间技术研究院，
类型：发明
国别省市：北京;11