一种基于Radon变换的运动目标径向速度估计方法
【技术领域】
[0001]本发明提供了一种基于Radon变换的运动目标径向速度估计方法，属于雷达信号 处理技术领域，用于雷达目标检测、参数估计、运动目标成像和定位等。
【背景技术】
[0002] 在合成孔径雷达（synthetic aperture radar,SAR)-运动目标检测系统中，运动 目标径向速度无模糊估计是实现运动目标成像和精确定位的必要前提。
[0003]运动目标径向速度无模糊估计分为多通道和单通道方法。多通道方法在沿航迹干 涉系统构型下，通过多载频或多基线构造多组包含运动目标径向速度的干涉相位，利用不 同载频或基线的参差关系解径向速度模糊。但多通道方法面临以下问题：系统实现成本较 高，通道间误差和配准误差等非理想因素对估计性能影响较大。传统的单通道方法基于包 络信息，利用距离走动与运动目标径向速度线性相关的思想，估计运动目标的径向速度。由 于包络信息不受相位缠绕问题影响，所以估计的径向速度是无模糊的。在单通道SAR系统 下，利用Radon变换通过二维搜索估计运动目标信号包络的斜率，可以实现运动目标径向速 度无模糊估计。但基于包络信息的单通道估计方法存在以下三点问题：（1)由于目标包络信 息并不敏感于运动目标参数，所以估计精度较低；（2)通过二维搜索估计参数，运算量较大；
[3] 二维搜索步长和参数估计精度需要折衷考虑。
【发明内容】
[0004] 本发明的目的在于:克服现有技术的不足，提出一种基于Radon变换的运动目标径 向速度快速且稳健的估计方法，以提高系统实现的实时性和稳健性。
[0005] 本发明采用的技术方案为：
[0006] 一种基于Radon变换的运动目标径向速度估计方法，步骤如下：
[0007] (1)对利用SAR接收的运动目标回波信号进行距离向脉压处理，得到运动目 标距离脉压域信号，其中，f为快时间，tm为慢时间；
[0008] (2)计算运动目标距离走动倾角Θ的最大范围[-0m，0m]，该最大范围即为Radon变换 角度范围；
[0009] 其中，9m=atan(2fsrvrm/c · fsa)，fsr为SAR系统采样频率，fsASAR系统脉冲重复频 率，C为光速，v rm为所述运动目标径向速度vj/f属范围[_vrm，vrm]的边界值，其中，v r = c · tan θ · fsa/2fsr；
[0010] (3)根据步骤(2)得到的Radon变换角度范围，得到N次Radon变换的变换角度分别 为：…，αΝ]，其中，= ΘΜ，[αι，α2，α3,…，αΝ]为等间隔排列或者为不等间 隔排列；
[0011] (4)对运动目标距离脉压域信号.以/乂,；)进行变换角度分别为[αι，α2，α3,…，α Ν]的Ν 次Radon变换，得到Ν次Radon变换结果Rn，η = 1，2，…，Ν，Rn对应的归一化Radon变换结果分别 为m · · _ J、；且瓦=兄/max (/?")，其中，max( ·)表示取最大值操作；
[0012] (5)将步骤(4)中得到的归一化Radon变换结果H疙,…,A分别与预设门限值 进行比较，6,/U,,·.·,瓦中大于门限值的数值的个数分别记为。，^，。，…，。，即为啦欠 Radon变换结災H…Λ的归一化长度;所述预设门限值为0.4-0.6。
[0013] (6)由Ν次Radon变换的变换角度[cn，α2，α3，…，α Ν]构造观测矩阵A
[0015]由N次Radon变换的归一化长度U，L2，L3，…，Ln构造观测值L
[0017] (7)利用最小二乘法，估计包含运动目标距离走动倾角Θ信息的矩阵
[0019]其中，[· Γ1表示矩阵求逆运算，[· ]τ表示矩阵转置操作，Le为运动目标距离走动 轨迹在二维平面中的长度；
[0020]
计算运动目标径向速度估计值 ?ν，从而完成所述基于Radon变换的运动目标径向速度估计。
[0021] 本发明与现有技术相比具有以下优点：
[0022] 1)与多通道估计方法相比，本发明利用单通道实现径向速度估计，降低了系统实 现成本，并且不受通道间误差和配准误差的影响；
[0023] 2)与基于相位的估计方法相比，本发明不存在模糊估计的问题；
[0024] 3)与基于包络的单通道估计方法相比，在运算量方面，避免了Radon变换的角度搜 索过程，利用少数几次Radon变换即可实现运动目标径向速度的快速估计，将现有Radon变 换方法的二维搜索降低为数次一维搜索，提高了参数估计的实时性；
[0025] 4)由于不存在Radon变换角度搜索过程，运动目标径向速度的估计精度不受角度 搜索步长的限制，可以使得估计精度得到进一步提升；
[0026] 5)通过最小二乘方法，使得本发明具有较好的稳健性。
【附图说明】
[0027] 图1是本发明方法的流程图；
[0028]图2是本发明与现有Radon变换方法在不同Radon变换次数下的对比仿真图；
[0029]图3为本发明与现有Radon变换方法在不同信噪比条件下的对比仿真图；
【具体实施方式】
[0030]合成孔径雷达-运动目标检测（SAR-GMTI)系统广泛应用于战场侦察、交通管制等 应用中。在SAR-GMTI系统中，运动目标的径向速度导致SAR成像处理后运动目标的方位位置 发生偏移;运动目标径向速度导致的运动目标距离走动使得运动目标成像结果散焦，因此， 估计运动目标径向速度是实现运动目标精确定位和高分辨率成像的必要前提。
[0031]沿航迹干涉(ATI)系统利用了径向速度与干涉相位的关系，常用于估计运动目标 径向速度。在常规沿航迹干涉系统中，由于干涉相位存在相位缠绕问题，当运动目标径向速 度超过系统最大不模糊速度时，径向速度估计就会出现模糊。为了在ATI系统中实现径向速 度无模糊估计，研究人员通过利用多波长、多基线、多速度SAR等方法实现解模糊。但这些方 法需要在多通道SAR系统上利用多波长、多基线，或者利用不同速度的SAR运动平台，构成多 组径向速度与干涉相位关系，实现相位解缠绕，从而无模糊的估计运动目标径向速度。但这 些方法不仅增加系统成本、灵活性较差，而且估计性能精度受通道间误差影响很大。为避免 以上问题，研究人员利用Radon变换估计运动目标径向速度。由于运动目标导致距离走动， 该距离走动轨迹在二维平面上呈现为一条斜线，而运动目标距离走动轨迹的倾斜角度Θ(简 称距离走动倾角）与运动目标径向速度Vr具有如下关系Vr = C · tan9 · fsa/2fsr，其中，fsr为 SAR系统采样频率，fsa*SAR系统脉冲重复频率，C为光速。由于Radon变换能够估计二维平面 上直线的倾角，研究人员利用Radon变换通过搜索所有可能的角度：当Radon变换角度与运 动目标轨迹距离走动倾角最为接近时，在Radon变换域的结果最大，此时的Radon变换角度 即可用于估计运动目标径向速度，这就是现有Radon变换估计方法的核心思想。但在实际应 用中，由于Radon变换涉及到距离和角度的二维搜索问题，其运算量较大;而且搜索步长与 径向速度估计精度和运算量都有关系:搜索步长越小，估计径向速度所需的运算量越大，实 现的估计精度越高;而搜索步长越大，运算量越小，估计精度越差该方法实现简单，但面临 由于二维搜索而带来的运算量大的问题，而且二维搜索步长与参数估计精度有关，在快速 估计参数的应用场合需要折衷选择二维搜索步长，不便于工程实现。
[0032]针对以上问题，本发明公开了一种基于Radon变换的运动目标径向速度估计方法， 该方法避免了 Radon变换的角度搜索过程，利用少数几次Radon变换即可实现运动目标径向 速度估计，将现有基于Radon变换的的运动目标径向速度估计方法由二维搜索降低为数次 一维搜索，提高了参数估计的实时性；由于不存在Radon变换角度搜索过程，参数估计的精 度不受角度搜索步长的影响，可以使得估计精度得到进一步提升;该方法通过最小二乘方 法估计参数，使得本发明具有较好的稳健性。采用本发明处理方式运算量较小，估计精度更 高，而且在较低SNR条件下即可实现较高的估计精度。因此，本发明可以满足高实时性和稳 健性要求。
[0033] 参照图1，本发明的具体实现步骤如下：
[0034] 步骤1，获取运动目标距离脉压域信号5@/,"):。
[0035] la)利用单通道合成孔径雷达(SAR)接收运动目标回波信号xf/j,,,；);
[0037] 其中，〇为运动目标的散射系数，,为距离窗函数，W(tm)为方位窗函数，exp (·)表示指数运算，j为虚部符号，f为快时间，tm为慢时间，γ为发射线性调频信号的调频 率，fc为载波中心频率，c为光速，Rb为运动目标到雷达平台运行轨道的最近距离，vr和V分别 为运动目标径向速度和雷达平台运行速度：
为运动目标到雷达平台的 瞬时斜距；
[0038] lb)对接收到的运动目标回波信号进行距离向傅立叶变换，得到运动目标 距离向频率域信号X( f r，tm)为：
[0039]
[0040] 其中，A(fr，tm)为运动目标距离向频率域信号X(fr，t m)的幅度，fr为距离频率；
[0041] 1C)将运动目标距离向频率域信号X(fr，U)乘以距离向匹配函数Sr(fr)，并通过距 离向逆傅里叶变换，得到运动目标距离脉压域信号5'(?，/,,,；)为：
[0043]其中，人;)为运动目标距离脉压域信号的幅度，距离向匹配函数Sr(fr) 的表不公式为：
[0045] 步骤2，计算运动目标距离走动倾角Θ的最大范围[-0m，0m]。
[0