发明名称 基于二次拟合距离方程的地面加速运动目标成像方法
摘要 本发明属于雷达技术领域,公开了一种基于二次拟合距离方程的地面加速运动目标成像方法,用于合成孔径雷达对地面加速运动目标成像,其步骤为:1)对目标距离方程进行二次拟合;2)将目标原始回波信号变到二维频域;3)根据拟合出的二次距离方程,在二维频域构造距离匹配滤波器,并进行距离压缩和距离徙动校正;4)将距离压缩和距离徙动校正后的信号变到距离多普勒域;5)在距离多普勒域构造方位匹配滤波器,并进行方位压缩;6)对方位压缩后的信号进行方位向逆傅里叶变换,完成对目标的成像。本发明不需知道目标的速度参数、位置参数和加速度,就能对地面加速运动目标进行精确成像,而且计算效率高。
申请公布号 CN103983969B 申请公布日期 2017.02.15
申请号 CN201410143271.2 申请日期 2014.04.10
申请人 西安电子科技大学 发明人 王彤;李永康;鲁缘政;任倩倩
分类号 G01S13/90(2006.01)I;G06F19/00(2011.01)I 主分类号 G01S13/90(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于二次拟合距离方程的地面加速运动目标成像方法,其特征在于,包括以下步骤:步骤1,对雷达的目标距离方程进行二次拟合,得到二次拟合后的距离方程:<maths num="0001"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mn>0</mn></msub><mo>+</mo><msub><mi>a</mi><mn>1</mn></msub><msub><mi>t</mi><mi>a</mi></msub><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><msubsup><mi>t</mi><mi>a</mi><mn>2</mn></msubsup></mrow>]]></math><img file="FDA0001120268950000011.GIF" wi="438" he="63" /></maths>其中,雷达工作在正侧视模式下,雷达平台速度为v<sub>a</sub>,t<sub>a</sub>为慢时间;t<sub>a</sub>=0时,目标方位向速度、方位向加速度、距离向速度、距离向加速度分别为v<sub>x</sub>、a<sub>x</sub>、v<sub>y</sub>和a<sub>y</sub>,且此时,雷达位于坐标原点,目标位于(0,y<sub>0</sub>);R(t<sub>a</sub>)为t<sub>a</sub>时刻目标到雷达的瞬时距离,t<sub>a</sub>为慢时间,a<sub>0</sub>、a<sub>1</sub>和a<sub>2</sub>分别为常数项、一次项和二次项系数;目标的原始回波信号表示为:<maths num="0002"><math><![CDATA[<mrow><mi>s</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>,</mo><msub><mi>t</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow><msub><mi>w</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>r</mi></msub><mo>-</mo><mn>2</mn><mi>R</mi><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo><mo>/</mo><mi>c</mi><mo>)</mo></mrow><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mn>4</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mo>+</mo><msub><mi>j&pi;K</mi><mi>r</mi></msub><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>r</mi></msub><mo>-</mo><mfrac><mrow><mn>2</mn><mi>R</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>a</mi></msub><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>}</mo></mrow>]]></math><img file="FDA0001120268950000012.GIF" wi="1478" he="180" /></maths>其中,t<sub>r</sub>为快时间,c为光速,w<sub>a</sub>(t<sub>a</sub>)、w<sub>r</sub>(t<sub>r</sub>)分别为目标原始回波信号的方位包络和距离包络,f<sub>c</sub>为雷达发射信号的载频,K<sub>r</sub>为雷达发射的线性调频信号的调频率;步骤2,对地面加速度运动目标的原始回波信号依次进行距离向傅里叶变换和方位向傅里叶变换,得到二维频域目标回波信号;其中,在步骤2中,地面加速度运动目标的原始回波信号经距离向傅里叶变换和方位向傅里叶变换后,得到二维频域目标回波信号为<maths num="0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>S</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>a</mi></msub><mo>,</mo><msub><mi>f</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>W</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>a</mi></msub><mo>-</mo><msub><mi>f</mi><mrow><mi>a</mi><mi>c</mi></mrow></msub><mo>)</mo></mrow><msub><mi>W</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>r</mi></msub><mo>)</mo></mrow><mi>exp</mi><mo>{</mo><mrow><mo>-</mo><mi>j</mi><mi>&pi;</mi><mrow><mo>(</mo><mrow><mfrac><msubsup><mi>f</mi><mi>r</mi><mn>2</mn></msubsup><msub><mi>K</mi><mi>r</mi></msub></mfrac><mo>-</mo><mfrac><mrow><msubsup><mi>cf</mi><mi>a</mi><mn>2</mn></msubsup></mrow><mrow><mn>4</mn><msub><mi>a</mi><mn>2</mn></msub><msubsup><mi>f</mi><mi>c</mi><mn>2</mn></msubsup></mrow></mfrac><msub><mi>f</mi><mi>r</mi></msub><mo>+</mo><mfrac><mrow><mn>4</mn><msub><mi>a</mi><mn>0</mn></msub></mrow><mi>c</mi></mfrac><msub><mi>f</mi><mi>r</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&times;</mo><mi>exp</mi><mrow><mo>{</mo><mrow><mi>j</mi><mi>&pi;</mi><mrow><mo>&lsqb;</mo><mrow><mfrac><mi>c</mi><mrow><mn>4</mn><msub><mi>f</mi><mi>c</mi></msub><msub><mi>a</mi><mn>2</mn></msub></mrow></mfrac><msubsup><mi>f</mi><mi>a</mi><mn>2</mn></msubsup><mo>+</mo><mfrac><msub><mi>a</mi><mn>1</mn></msub><msub><mi>a</mi><mn>2</mn></msub></mfrac><msub><mi>f</mi><mi>a</mi></msub><mo>-</mo><mfrac><mrow><mn>4</mn><msub><mi>a</mi><mn>0</mn></msub><msub><mi>a</mi><mn>2</mn></msub><mo>-</mo><msubsup><mi>a</mi><mn>1</mn><mn>2</mn></msubsup></mrow><mrow><msub><mi>ca</mi><mn>2</mn></msub></mrow></mfrac><msub><mi>f</mi><mi>r</mi></msub></mrow><mo>&rsqb;</mo></mrow></mrow><mo>}</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001120268950000013.GIF" wi="1356" he="328" /></maths>其中,W<sub>r</sub>(f<sub>r</sub>)是目标信号距离频谱的包络,W<sub>a</sub>(f<sub>a</sub>)是目标信号方位频谱的包络,f<sub>r</sub>为距离频率,f<sub>a</sub>是方位频率,f<sub>ac</sub>为目标多普勒中心频率;步骤3,根据拟合出的二次距离方程,在二维频域构造距离匹配滤波器,对二维频域目标回波信号进行距离向匹配滤波,实现距离压缩和距离徙动校正;步骤4,对二维频域距离匹配滤波后的目标回波信号进行距离向逆傅里叶变换,得到距离多普勒域目标回波信号;步骤5,在距离多普勒域构造方位匹配滤波器,对距离多普勒域目标回波信号进行方位向匹配滤波,实现方位压缩;步骤6,对方位匹配滤波后的距离多普勒域目标回波信号进行方位向逆傅里叶变换,完成对地面运动目标的成像。
地址 710071 陕西省西安市太白南路2号