发明名称 基于相干化处理的步进频率ISAR成像方法
摘要 本发明是一种基于相干化处理的步进频率ISAR成像方法。主要解决传统参数估计方法运算效率较低,包络补偿后存在的相位误差对频率合成有较大影响的问题。其实现步骤为:建立步进频率ISAR回波信号模型;根据回波信号模型进行运动参数估计;构造补偿函数对回波信号的子脉冲包络进行补偿;对补偿完的回波信号进行基于相干化处理的频率合成;对合成后的回波信号频谱进行IFFT变换,获取高分辨一维距离像;利用RD算法或RID算法进行方位成像,获取方位的高分辨图像。实验证明本发明具有高相位误差补偿精度和高运算效率的优点,可用于实现目标的二维高分辨ISAR成像。
申请公布号 CN102998672A 申请公布日期 2013.03.27
申请号 CN201210493356.4 申请日期 2012.11.27
申请人 西安电子科技大学 发明人 李亚超;全英汇;邢孟道;许斌;周瑞雨
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于相干化处理的步进频率ISAR成像方法,包括如下步骤:(1)建立步进频率ISAR回波信号模型;(2)根据回波信号模型进行运动参数估计:2a)将不同脉组的第一个子脉冲数据分别取出,得到回波信号相邻相关表达式为:<maths num="0001"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><msub><mi>mT</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&sigma;</mi><mi>P</mi></msub><mi>sin</mi><mi>c</mi><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>&tau;</mi><mo>+</mo><mfrac><mrow><mn>2</mn><mi>M</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>M</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>M</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,τ为相关时间,m为脉组数,0≤m≤M-1,mT<sub>a</sub>为方位时间的离散表示,σ<sub>P</sub>为目标回波后向散射系数,A为相关匹配量,c为光速,λ为信号波长,M=M<sub>1</sub>+M<sub>2</sub>,<maths num="0002"><![CDATA[<math><mrow><msub><mi>M</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>v</mi><mi>r</mi></msub><msub><mi>x</mi><mi>P</mi></msub><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>T</mi><mi>a</mi></msub><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>a</mi><mi>r</mi></msub><msubsup><mi>T</mi><mi>a</mi><mn>2</mn></msubsup><mo>,</mo></mrow></math>]]></maths>M<sub>2</sub>=T<sub>a</sub>a<sub>r</sub>mT<sub>a</sub>其中,v<sub>r</sub>是目标相对雷达运动的径向速度,x<sub>P</sub>是散射点的横向位置信息,ω为目标相对雷达运动的转动角速度,T<sub>a</sub>为脉组间的重复周期,a<sub>r</sub>是目标相对雷达运动的径向加速度。2b)对回波信号相邻相关表达式沿着方位向做FFT积累,得到方位频域回波信号相邻相关表达式为:<maths num="0003"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>&tau;</mi><mo>,</mo><msub><mi>f</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&sigma;</mi><mi>P</mi></msub><mi>sin</mi><mi>c</mi><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>&tau;</mi><mo>+</mo><mfrac><mrow><mn>2</mn><mi>M</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><mi>sin</mi><mi>c</mi><mrow><mo>(</mo><msub><mi>&pi;T</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>a</mi></msub><mo>-</mo><mfrac><msub><mrow><mn>2</mn><mi>T</mi></mrow><mi>a</mi></msub><mi>&lambda;</mi></mfrac><msub><mi>a</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>M</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,T<sub>m</sub>为子脉冲重复周期,f<sub>a</sub>为方位向频率,<img file="FDA00002477120000014.GIF" wi="324" he="115" />是回波信号相邻相关表达式的相关峰值位置;忽略M中相邻两次回波ω引起的包络变化,则相关峰值位置又可表示为<maths num="0004"><![CDATA[<math><mrow><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><msub><mi>T</mi><mi>a</mi></msub><mo>+</mo><msub><mi>a</mi><mi>r</mi></msub><msubsup><mi>T</mi><mi>a</mi><mn>2</mn></msubsup></mrow><mi>c</mi></mfrac><mo>,</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>T</mi></mrow><mi>a</mi></msub><msub><mi>a</mi><mi>r</mi></msub></mrow><mi>&lambda;</mi></mfrac><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>2c)将相关峰值位置在距离向和方位向的坐标分别设为<img file="FDA00002477120000016.GIF" wi="32" he="33" />和<img file="FDA00002477120000017.GIF" wi="46" he="35" />计算目标的径向加速度a<sub>r</sub>和速度v<sub>r</sub>:<maths num="0005"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>a</mi><mi>r</mi></msub><mo>=</mo><mfrac><mrow><mi>&lambda;prf</mi><mo>[</mo><mover><mi>n</mi><mo>~</mo></mover><mo>-</mo><mrow><mo>(</mo><mfrac><mover><mi>N</mi><mo>~</mo></mover><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><msub><mrow><mn>2</mn><mi>NT</mi></mrow><mi>a</mi></msub></mfrac></mtd><mtd><mover><mi>n</mi><mo>~</mo></mover><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mover><mi>N</mi><mo>~</mo></mover></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>r</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mi>c</mi><mo>[</mo><mover><mi>m</mi><mo>~</mo></mover><mo>-</mo><mrow><mo>(</mo><mfrac><mover><mi>M</mi><mo>~</mo></mover><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mrow><msub><mrow><mn>2</mn><mi>f</mi></mrow><mi>s</mi></msub><msub><mi>T</mi><mi>a</mi></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><mi>&lambda;</mi><mo>[</mo><mover><mi>n</mi><mo>~</mo></mover><mo>-</mo><mrow><mo>(</mo><mfrac><mover><mi>N</mi><mo>~</mo></mover><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mn>4</mn><mover><mi>N</mi><mo>~</mo></mover><msub><mi>T</mi><mi>a</mi></msub></mrow></mfrac></mtd><mtd><mover><mi>m</mi><mo>~</mo></mover><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mover><mi>M</mi><mo>~</mo></mover></mtd></mtr></mtable></mfenced></math>]]></maths>其中,<img file="FDA00002477120000022.GIF" wi="42" he="39" />和<img file="FDA00002477120000023.GIF" wi="35" he="40" />分别为距离向和方位向的离散点数,prf为方位向采样频率,f<sub>s</sub>为距离向采样频率,T<sub>r</sub>为脉冲重复周期;(3)根据径向加速度a<sub>r</sub>和速度v<sub>r</sub>,利用如下包络补偿的相位补偿函数s<sub>sref</sub>对回波信号的子脉冲包络进行补偿,<maths num="0006"><![CDATA[<math><mrow><msub><mi>s</mi><mi>sref</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>c</mi></mfrac><mrow><mo>(</mo><msub><mi>v</mi><mi>r</mi></msub><msub><mi>f</mi><mi>r</mi></msub><msub><mi>t</mi><mi>m</mi></msub><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>a</mi><mi>r</mi></msub><msub><mi>f</mi><mi>r</mi></msub><msubsup><mi>t</mi><mi>m</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>其中,f<sub>r</sub>为距离向频率,t<sub>m</sub>为方位向时间;(4)对完成子脉冲包络补偿的回波信号进行基于相干化处理的频率合成,得到合成后的回波信号频谱s(f<sub>r</sub>,T<sub>r</sub>);(5)对合成后的回波信号频谱s(f<sub>r</sub>,T<sub>r</sub>)进行逆傅里叶IFFT变换,获取高分辨一维距离像;(6)对高分辨一维距离像利用RD算法或RID算法进行方位成像,获取目标方位的高分辨图像。
地址 710071 陕西省西安市太白南路2号