发明名称 一种对非协同目标高精度的ISAR运动补偿方法
摘要 一种对非协同目标高精度的ISAR运动补偿方法,该方法在传统的包络相关方法的基础上重新定义参考距离像,然后计算待对准距离像与参考距离像之间的互相关函数EC(τ<sub>m+1</sub>),其中m=0,1,…,M-1,M为距离像组数,τ<sub>m+1</sub>为距离偏移量。运动补偿方法的核心技术是精确的估计每一个距离像与参考距离像之间的偏移量,即:<img file="dda00001744307600011.GIF" wi="495" he="77" />从而消除距离徙动;再对每个距离像中的相位矢量进行估计和补偿,消除目标的相位误差;从而可以得到运动补偿后的高精度ISAR图像。采用本发明的方法可以显著消除目标运动引起的成像模糊问题,是一种为高速运动目标提供高分辨率的ISAR图像的有效方法。
申请公布号 CN102707282A 申请公布日期 2012.10.03
申请号 CN201210189776.3 申请日期 2012.06.08
申请人 北京航空航天大学 发明人 陈娟;袁运能
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 成金玉
主权项 1.一种对非协同目标高精度的ISAR运动补偿方法,其特征在于实现步骤如下:(1)由步进频率逆合成孔径雷达获得大小为M*N的非协同目标回波数据E[m,n],其中M表示传感器接收的总脉冲串组数,N表示每一组脉冲串中步进频率脉冲的个数;(2)对目标回波数据沿距离向进行脉冲压缩,得到M组长度为N的距离像RP<sub>m</sub>(n),m=0,1,…,M-1;n=0,1,…,N-1,其中m表示脉冲串数,n表示脉冲数;(3)定义第一组对准距离像<img file="FDA00001744307300011.GIF" wi="349" he="65" />式中RP<sub>0</sub>(n)表示第一个组未对准距离像;(4)由第m组已对准距离像<img file="FDA00001744307300012.GIF" wi="145" he="65" />和参考距离像RP<sub>refm</sub>(n),按如下公式计算出下一组参考距离像RP<sub>refm+1</sub>(n):<maths num="0001"><![CDATA[<math><mrow><msub><mi>RP</mi><mrow><mi>refm</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>m</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></mfrac><msub><mi>RP</mi><mi>refm</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>|</mo><mover><mrow><msub><mi>RP</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>(5)将RP<sub>m</sub>(n)右移τ<sub>m+1</sub>,计算RP<sub>refm+1</sub>(n)与RP<sub>m+1</sub>(n-τ<sub>m+1</sub>)之间的互相关函数:<maths num="0002"><![CDATA[<math><mrow><mi>EC</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>|</mo><msub><mi>RP</mi><mrow><mi>refm</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo><mo>&CenterDot;</mo><mo>|</mo><msub><mi>RP</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中τ<sub>m+1</sub>为[0,1,…,N-1]中的整数;(6)计算EC(τ<sub>m+1</sub>)取得最大值时,τ<sub>m+1</sub>在[0,1,…,N-1]中的值,记为τ<sub>m+1,0</sub>;(7)根据Nelder-Mead法则,将τ<sub>m+1,0</sub>作为迭代技术的初始猜测值,得到最佳距离徙动值<img file="FDA00001744307300015.GIF" wi="96" he="55" /><maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>arg</mi><munder><mi>max</mi><msub><mi>&tau;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub></munder><mi>EC</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA00001744307300017.GIF" wi="73" he="54" />不一定为整数;(8)计算<img file="FDA00001744307300018.GIF" wi="199" he="66" />如下式所示:<maths num="0004"><![CDATA[<math><mrow><mover><mrow><msub><mi>RP</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>=</mo><msub><mi>RP</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>当<img file="FDA000017443073000110.GIF" wi="72" he="54" />不是整数时,由傅里叶变换的平移性质,按如下方法实现上述过程:<maths num="0005"><![CDATA[<math><mrow><msub><mi>RP</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>FFT</mi><mo>{</mo><msup><mi>e</mi><mrow><mi>j</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mi>q</mi></mrow></msup><mi>IFFT</mi><mo>{</mo><msub><mi>RP</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>}</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>q表示向量[0,1,…,N-1]<sup>T</sup>,[·]<sup>T</sup>表示向量转置;(9)令m=m+1,如果m<M-1,重复以上(4)-(8)进行下一个距离像的对准;(10)所有距离像对准完毕后,进行相位校正:<img file="FDA00001744307300021.GIF" wi="1318" he="110" />式中c为光速,f<sub>n</sub>为第n个脉冲的频率,<img file="FDA00001744307300022.GIF" wi="78" he="55" />为矫正后相位。
地址 100191 北京市海淀区学院路37号