发明名称 基于单通道SAR的运动目标径向速度无模糊估计方法
摘要 本发明公开了一种基于单通道SAR的运动目标径向速度无模糊估计方法,主要解决现有技术中实现成本较高、运算量大、估计精度低的问题。其实现过程为:对雷达接收的运动目标回波信号进行距离脉压,将距离脉压结果转换到距离频率域,通过两视处理和距离向逆傅立叶变换得到运动目标距离脉压域两视信号,再将运动目标距离脉压域两视信号分别转换到多普勒频域,乘以方位向匹配函数,再经过干涉处理得到运动目标干涉结果;利用最小二乘线性拟合方法估计干涉结果的干涉相位随多普勒频率变化的斜率,根据该斜率估计运动目标径向速度。本发明具有系统实现成本低、运算量小、估计精度高、能同时估计两个运动目标径向速度的优点,可用于运动目标成像和精确定位。
申请公布号 CN103364783A 申请公布日期 2013.10.23
申请号 CN201310279714.6 申请日期 2013.07.04
申请人 西安电子科技大学 发明人 廖桂生;张学攀;朱圣棋;王兰美;王桂宝;杨东;张学敬
分类号 G01S13/90(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于单通道SAR的运动目标径向速度无模糊估计方法,包括如下步骤:(1)对利用单通道SAR接收的运动目标回波信号<img file="FDA00003465868400011.GIF" wi="148" he="89" />进行距离向脉压处理,得到运动目标距离脉压域信号<img file="FDA00003465868400012.GIF" wi="176" he="87" />其中,<img file="FDA00003465868400013.GIF" wi="35" he="64" />为快时间,t<sub>m</sub>为慢时间;(2)将运动目标距离脉压域信号<img file="FDA00003465868400014.GIF" wi="148" he="88" />经过距离向傅立叶变换,得到运动目标距离频率域信号S(f<sub>r</sub>,t<sub>m</sub>),其中,f<sub>r</sub>为距离频率;(3)对运动目标距离频率域信号S(f<sub>r</sub>,t<sub>m</sub>)通过两视处理,得到第一视数信号S<sub>1</sub>(f<sub>r</sub>,t<sub>m</sub>)和第二视数信号S<sub>2</sub>(f<sub>r</sub>,t<sub>m</sub>);对运动目标距离频率域的两视数信号再经过距离向逆傅里叶变换,得到运动目标距离脉压域的第一视数信号S<sub>1</sub>(Δf,t<sub>m</sub>)和第二视数信号S<sub>2</sub>(Δf,t<sub>m</sub>),其中,Δf为距离频率域的频谱宽度;(4)将运动目标距离脉压域的第一视数信号S<sub>1</sub>(Δf,t<sub>m</sub>)和第二视数信号S<sub>2</sub>(Δf,t<sub>m</sub>)分别通过方位向傅立叶变换转换到多普勒频域,并分别乘以第一视数信号的多普勒频域方位向匹配函数s<sub>a1</sub>(Δf,f<sub>d1</sub>)和第二视数信号的多普勒频域方位向匹配函数s<sub>a2</sub>(Δf,f<sub>d2</sub>),得到运动目标多普勒频域的第一视数信号S<sub>1</sub>(Δf,f<sub>d1</sub>)和第二视数信号S<sub>2</sub>(Δf,f<sub>d2</sub>),其中,f<sub>d1</sub>和f<sub>d2</sub>分别为运动目标多普勒频域的第一视数信号和第二视数信号的多普勒频率;(5)将运动目标多普勒频域的第一视数信号S<sub>1</sub>(Δf,f<sub>d1</sub>)和第二视数信号S<sub>2</sub>(Δf,f<sub>d2</sub>)进行干涉处理,得到运动目标干涉结果ΔS为:<maths num="0001"><![CDATA[<math><mrow><mi>&Delta;S</mi><mo>=</mo><mi>A</mi><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;&Delta;f</mi></mrow><mi>c</mi></mfrac><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>B</mi></msub><mo>-</mo><mfrac><msubsup><mi>V</mi><mi>r</mi><mn>2</mn></msubsup><msup><mi>V</mi><mn>2</mn></msup></mfrac><msub><mi>R</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>d</mi></msub><msub><mi>R</mi><mi>B</mi></msub><msub><mi>V</mi><mi>r</mi></msub></mrow><msup><mi>V</mi><mn>2</mn></msup></mfrac><mfrac><mi>&Delta;f</mi><msub><mi>f</mi><mi>c</mi></msub></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,A为运动目标干涉结果ΔS的幅度,exp(·)表示指数运算,j为虚部符号,f<sub>d</sub>为运动目标信号多普勒频率,f<sub>c</sub>为载波中心频率,c为光速,R<sub>B</sub>为运动目标到雷达平台运行轨道的最近距离,V<sub>r</sub>和V分别为运动目标径向速度和雷达平台运行速度;(6)利用最小二乘线性拟合方法估计干涉结果ΔS的干涉相位<img file="FDA00003465868400016.GIF" wi="46" he="55" />随多普勒频率f<sub>d</sub>变化的斜率<img file="FDA00003465868400017.GIF" wi="69" he="76" />由该斜率<img file="FDA00003465868400018.GIF" wi="36" he="75" />计算运动目标的径向速度估计值<img file="FDA00003465868400019.GIF" wi="46" he="84" />为:<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>^</mo></mover><mi>r</mi></msub><mo>=</mo><mover><mi>k</mi><mo>^</mo></mover><msup><mi>V</mi><mn>2</mn></msup><msub><mi>f</mi><mi>c</mi></msub><mo>/</mo><mn>2</mn><mi>&pi;</mi><msub><mi>R</mi><mi>B</mi></msub><mi>&Delta;f</mi><mo>.</mo></mrow></math>]]></maths>
地址 710071 陕西省西安市太白南路2号