主权项 |
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>‾</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>τ</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</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>·</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>τ</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>τ</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>τ</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msub></munder><mi>EC</mi><mrow><mo>(</mo><msub><mi>τ</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>‾</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>τ</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>τ</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>π</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow><msub><mover><mi>τ</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" />为矫正后相位。 |