主权项 |
1.一种基于相干化处理的步进频率ISAR成像方法,包括如下步骤:(1)建立步进频率ISAR回波信号模型;(2)根据回波信号模型进行运动参数估计:2a)将不同脉组的第一个子脉冲数据分别取出,得到回波信号相邻相关表达式为:<maths num="0001"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>τ</mi><mo>,</mo><msub><mi>mT</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>σ</mi><mi>P</mi></msub><mi>sin</mi><mi>c</mi><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>τ</mi><mo>+</mo><mfrac><mrow><mn>2</mn><mi>M</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>×</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>π</mi></mrow><mi>λ</mi></mfrac><msub><mi>M</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>×</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>π</mi></mrow><mi>λ</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>ω</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>τ</mi><mo>,</mo><msub><mi>f</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>σ</mi><mi>P</mi></msub><mi>sin</mi><mi>c</mi><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>τ</mi><mo>+</mo><mfrac><mrow><mn>2</mn><mi>M</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>×</mo><mi>sin</mi><mi>c</mi><mrow><mo>(</mo><msub><mi>π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>λ</mi></mfrac><msub><mi>a</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>×</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>π</mi></mrow><mi>λ</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>λ</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>λ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>λ</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>π</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算法进行方位成像,获取目标方位的高分辨图像。 |