一种双基前视SAR线性时间同步误差估计方法技术

本发明专利技术提供了一种双基前视SAR线性时间同步误差估计方法，属于雷达技术领域。本发明专利技术首先基于方位向变量代换去除徙动引起的距离史斜率，然后根据最小熵理论，乘以校正因子校正残余斜率，当取得最小熵时，距离徙动校正效果达到最好，同时实现线性时间同步误差精确估计，本发明专利技术步骤简单，可嵌套在成像过程中进行，在得到线性时间同步误差的同时，也能够实现距离徙动的精确校正。

【技术实现步骤摘要】
一种双基前视SAR线性时间同步误差估计方法
本专利技术属于雷达
，特别涉及一种双基前视SAR线性时间同步误差估计方法。
技术介绍
合成孔径雷达利用微波遥感技术，能够克服云雾雨雪和黑暗条件的限制对地面目标成像，具有全天时、全天候工作特点，在军事侦察、地形测绘、地物分类识别、海洋应用以及资源勘探等领域被广泛应用。双基SAR收发分置，根据不同空间几何关系，可以灵活配置，具有作用距离远、抗干扰能力强和隐蔽性好等优点，而且能够实现前视成像；但同时由于收发系统分置，也带来空间同步、时间同步、相位同步等问题，其中时间同步误差，主要由固定误差、线性误差和随机误差构成，其中随机误差的影响通常可以忽略，固定误差主要造成距离测量误差，而线性误差不仅会造成距离测量的误差，还会造成回波跨距离单元，严重影响成像聚焦质量。因而线性时间同步误差是影响双基前视SAR成像质量好坏的关键因素，获得精确的线性时间同步误差对于双基前视SAR高精度成像具有重要意义。现有技术中，有一种利用直达波信号脉压峰值位置和相位信息提取线性时间同步误差的方法，该方法由于涉及到时间参考点、该时间参考点对应的收发站距离以及发射系统脉冲重复周期三个参数的估计，误差来源多，估计精度受限；另外，还有一种依靠直达波和基于幅度的多普勒中心频率估计来实现时间同步误差估计的方法，该方法需要通过求解三元方程组来实现时间同步误差估计，由于该方程组中包含取余运算，相当于增加了一个未知量，导致自由度不足，求解过程复杂且解不唯一，虽然可以利用有限个可列解逐一对回波信号进行成像，然后将成像结果对应最好的解作为最终估计结果，但实际中需要消耗更多时间和资源，对于实时性要求较高的系统，往往难以达到要求。
技术实现思路
本专利技术的目的是针对双基前视SAR收、发平台时钟源不同步导致的时间同步误差问题，研究设计一种线性时间同步误差估计方法，实现对线性时间同步误差估计的同时，也能够实现双基前视SAR距离徙动的精校正，对于双基前视SAR高精度成像具有重要意义。一种双基前视SAR线性时间同步误差估计方法，包括以下步骤：步骤1，发射线性调频信号，获取被照射区域内目标的回波信号，对所述回波信号进行距离向脉冲压缩；步骤2，对脉冲压缩后的回波信号进行方位向变量代换；步骤3，设置初始校正参数和迭代步长；步骤4，进行残余斜率校正；步骤5，沿方位向投影并获取投影序列熵；步骤6，当熵值为最小熵值时，输出线性时间同步误差系数。进一步地，所述步骤1包括以下流程：发射线性调频信号，获取被照射区域内目标的二维回波基带信号其中，τ表示快时间，t表示慢时间，Tp为脉冲宽度，fc为发射信号的载波频率，c为光速，Kr为线性调频信号的调频斜率，R(t)为双基距离和，ε(t)＝β·t为线性时间同步误差，β为线性时间同步误差系数；利用匹配滤波方法对回波信号S(t,τ)进行脉冲压缩，脉冲压缩后的信号表达式为其中，fτ为距离向频率，B为发射信号带宽。进一步地，所述步骤2包括以下流程：获取双基前视SAR接收机、发射机的运动速度、角度信息，得到多普勒中心频率fdc其中，为发射天线斜视角，为接收天线的下视角，vT为发射机速度，vR为接收机速度，λ为发射信号波长；计算多普勒模糊估计值MM＝round(fdc/PRF)其中，PRF为脉冲重复频率，round()为向上取整函数；根据所述多普勒模糊估计值M，构造多普勒模糊相位因子PP＝exp(j2π·PRF·M·t)所述多普勒模糊相位因子P与脉冲压缩后信号S(t,fτ)相乘，去除多普勒模糊，其表达式为将S'(t,fτ)中R(t)在波束中心穿越目标时刻进行泰勒展开，保留至一阶项可得其中，RR0、RT0为收发平台波束中心对准目标时的距离，v为收发平台速度；收发平台速度相等；将R(t)表达式代入S'(t,fτ)表达式，按照进行变量变换，以t'代替S'(t,fτ)中t的得到S'(t',fτ)，由于fτ＜＜fc，通过变量代换和近似处理后，信号表达式为对S'(t',fτ)进行距离向IFFT，可得进一步地，所述步骤3包括以下流程：设置初始校正系数，确定校正系数迭代步长进一步地，所述步骤4包括以下流程：将方位向变量代换后的信号，在方位时域、距离时域乘以残余斜率校正因子，进行残余斜率校正，残余斜率因子表达式为H(t',fτ)＝exp(j2π·fτ·β'·t')。进一步地，所述步骤5包括以下流程：将残余斜率校正后的信号沿方位向投影得到能量序列x(n)，n＝0,1,…,N-1，其中，N为距离向采样点数，计算残余斜率校正后熵值H(x)进一步地，所述步骤6包括以下流程：将校正系数β'分别相加Δβ和相减Δβ得到β'+和β'-，分别用β'+和β'-替换校正系数β'重复步骤4和5，分别得到对应熵值H+和H-，比较H和H+、H-值大小，若H为最小值，则当前值H即为最小熵；若H不是最小值，对校正系数进行更新。进一步地，所述对校正系数进行更新的流程包括：若H-为最小值，则将H更新为H-，β'更新为β'-，更新后的校正系数β'相减Δβ得到β'-，重复步骤4和5得到熵值H-；若H-≤H，逐次迭代，直到所得熵H-大于上一次迭代所得熵H为止，则H即为最小熵；若H+为最小值，则将H更新为H+，β'更新为β'+，更新后的校正系数β'相加Δβ得到β'+，重复步骤4和5得到熵值H+；若H+≤H，逐次迭代，直到所得熵H+大于上一次迭代所得熵H为止，则H即为最小熵；最小熵所对应的校正系数β'即为线性时间同步误差系数β，输出线性时间同步误差系数β。本专利技术的有益效果：本专利技术提供了一种双基前视SAR线性时间同步误差估计方法，本专利技术首先基于方位向变量代换去除徙动引起的距离史斜率，然后根据最小熵理论，乘以校正因子校正残余斜率，当取得最小熵时，距离徙动校正效果达到最好，同时实现线性时间同步误差精确估计，本专利技术步骤简单，可嵌套在成像过程中进行，在得到线性时间同步误差的同时，也能够实现距离徙动的精确校正。附图说明图1为本专利技术实施例的流程图。图2为本专利技术实施例中采用的双基前视SAR系统几何模型。图3为本专利技术实施例中回波脉冲压缩结果图。图4为本专利技术实施例中变量变换后脉冲压缩结果图。图5为本专利技术实施例中残余斜率校正后序列熵趋势图。图6为本专利技术实施例中最小熵对应的距离徙动校正结果图。具体实施方式本专利技术采用仿真实验验证所提出线性时间同步误差估计方法的正确性，所有步骤、结论均在Matlab2016平台验证正确。下面结合附图对本专利技术的实施例做进一步的说明。请参阅图1，本专利技术提供的一种双基前视SAR线性时间同步误差估计方法，具体通过以下步骤实现：步骤1，发射线性调频信号，获取被照射区域内目标的回波信号，对所述回波信号进行距离向脉冲压缩。本实施例中，双基前视SAR的几何模型如图2所示，其中P为场景中心，VT、VR为收发平台速度，且VT＝VR＝v，RR0、RT0为收发平台波束中心对准目标时的距离，为接收天线的下视角，为发射天线斜视角，RT(t)、RR(t)为发射机与接收机的距离历史。发射线性调频信号的调频信号为其中，τ表示快时间，Tp为脉冲宽度，fc为发射信号的载波频率，c为光速，Kr为线性调频信号的调频斜率。由于双基前视SAR收发分置，导致发射雷达脉冲重复周期与接收雷达脉冲重复周期之间存在时本文档来自技高网...

【技术保护点】
1.一种双基前视SAR线性时间同步误差估计方法，其特征在于，包括以下步骤：步骤1，发射线性调频信号，获取被照射区域内目标的回波信号，对所述回波信号进行距离向脉冲压缩；步骤2，对脉冲压缩后的回波信号进行方位向变量代换；步骤3，设置初始校正参数和迭代步长；步骤4，进行残余斜率校正；步骤5，沿方位向投影并获取投影序列熵；步骤6，当熵值为最小熵值时，输出线性时间同步误差系数。

【技术特征摘要】
1.一种双基前视SAR线性时间同步误差估计方法，其特征在于，包括以下步骤：步骤1，发射线性调频信号，获取被照射区域内目标的回波信号，对所述回波信号进行距离向脉冲压缩；步骤2，对脉冲压缩后的回波信号进行方位向变量代换；步骤3，设置初始校正参数和迭代步长；步骤4，进行残余斜率校正；步骤5，沿方位向投影并获取投影序列熵；步骤6，当熵值为最小熵值时，输出线性时间同步误差系数。2.如权利要求1所述的双基前视SAR线性时间同步误差估计方法，其特征在于，所述步骤1包括以下流程：发射线性调频信号，获取被照射区域内目标的二维回波基带信号exp(jπKr(τ-R(t)/c-ε(t))2)其中，τ表示快时间，t表示慢时间，Tp为脉冲宽度，fc为发射信号的载波频率，c为光速，Kr为线性调频信号的调频斜率，R(t)为双基距离和，ε(t)＝β·t为线性时间同步误差，β为线性时间同步误差系数；利用匹配滤波方法对回波信号S(t,τ)进行脉冲压缩，脉冲压缩后的信号表达式为其中，fτ为距离向频率，B为发射信号带宽。3.如权利要求2所述的双基前视SAR线性时间同步误差估计方法，其特征在于，所述步骤2包括以下流程：获取双基前视SAR接收机、发射机的运动速度、角度信息，得到多普勒中心频率fdc其中，为发射天线斜视角，为接收天线的下视角，vT为发射机速度，vR为接收机速度，λ为发射信号波长；计算多普勒模糊估计值MM＝round(fdc/PRF)其中，PRF为脉冲重复频率，round()为向上取整函数；根据所述多普勒模糊估计值M，构造多普勒模糊相位因子PP＝exp(j2π·PRF·M·t)所述多普勒模糊相位因子P与脉冲压缩后信号S(t,fτ)相乘，去除多普勒模糊，其表达式为exp(j2π·PRF·M·t)将S'(t,fτ)中R(t)在波束中心穿越目标时刻进行泰勒展开，保留至一阶项可得其中，RR0、RT0为收发平台波束中心对准目标时的距离，v为收发平台速度；将R(t)表达式代入S'(t,fτ)表达式，按照进行变量变换，以...

【专利技术属性】
技术研发人员：李文超，邹德锁，李易，缪羽宣，李中余，武俊杰，黄钰林，杨建宇，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51