发明名称 基于双重稀疏约束的逆合成孔径雷达自聚焦成像方法
摘要 本发明公开了一种基于双重稀疏约束的逆合成孔径雷达自聚焦成像的方法。主要解决传统逆合成孔径雷达成像方法中分辨率受发射带宽的限制以及回波受到相位误差干扰时导致成像结果出现散焦的问题。其技术方案是:1.雷达照射稀疏目标场景H,获得该场景的回波信号S<sub>r</sub>(τ,t);2.通过对回波信号S<sub>r</sub>(τ,t)二维离散采样得到离散回波矩阵S<sub>r</sub>,并构造含相位误差的回波矩阵S<sub>re</sub>;3.根据相位误差的回波矩阵S<sub>re</sub>,获得双重稀疏约束函数,并构造拉格朗日函数f(U,Z);4.用交替迭代法求解拉格朗日函数f(U,Z),得到场景H的图像。本发明能实现逆合成孔径雷达的高分辨成像,且结合相位误差校正实现自聚焦成像。
申请公布号 CN105093225A 申请公布日期 2015.11.25
申请号 CN201510527352.7 申请日期 2015.08.25
申请人 西安电子科技大学 发明人 赵光辉;刘永飞;沈方芳;石光明;左功玉
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;黎汉华
主权项 基于双重稀疏约束的逆合成孔径雷达自聚焦成像方法,包括如下步骤:(1)雷达不断向稀疏目标场景H发射线性调频脉冲信号S,并在发射脉冲的同时接收该场景H的回波脉冲,得到回波信号为S<sub>r</sub>(τ,t);(2)对回波信号S<sub>r</sub>(τ,t)进行二维离散采样,并构造含相位误差的回波矩阵:S<sub>re</sub>=D<sub>2</sub>AFBD<sub>1</sub>其中,A为距离向基矩阵,B为方位向基矩阵,F为散射系数矩阵,D<sub>1</sub>为方位向对角阵,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>D</mi><mn>1</mn></msub><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mo>&lsqb;</mo><mtable><mtr><mtd><msup><mi>e</mi><mrow><msub><mi>j&Phi;</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mtd><mtd><mn>...</mn></mtd><mtd><msup><mi>e</mi><mrow><msub><mi>j&Phi;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup></mtd><mtd><mn>...</mn></mtd><mtd><msup><mi>e</mi><mrow><msub><mi>j&Phi;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></msup></mtd></mtr></mtable><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0000788880160000011.GIF" wi="899" he="91" /></maths>Φ<sub>1</sub>为待求的方位向相位误差,Φ<sub>1</sub>=[Φ<sub>1</sub>(1) … Φ<sub>1</sub>(n) … Φ<sub>1</sub>(N)]<sup>T</sup>,D<sub>2</sub>为距离向对角阵,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>D</mi><mn>2</mn></msub><mo>=</mo><mi>diag</mi><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>e</mi><mrow><msub><mi>j&Phi;</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msup><mi>e</mi><mrow><msub><mi>j&Phi;</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></msup></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&Phi;</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>M</mi><mo>)</mo></mrow></mrow></msup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000788880160000012.GIF" wi="926" he="92" /></maths>Φ<sub>2</sub>为待求的距离向相位误差,Φ<sub>2</sub>=[Φ<sub>2</sub>(1) … Φ<sub>2</sub>(n) … Φ<sub>2</sub>(N)]<sup>T</sup>,T表示矩阵转置;(3)令方位向中间变量矩阵ψ=BD<sub>1</sub>、距离向中间变量矩阵<img file="FDA0000788880160000016.GIF" wi="214" he="93" />分离变量矩阵U=Fψ和中间变量矩阵Z=F<sup>T</sup>,得到如下双重稀疏约束函数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>U</mi><mo>,</mo><mi>Z</mi></mrow></munder><mrow><mo>{</mo><mrow><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>U</mi><mo>|</mo><msubsup><mo>|</mo><mrow><mn>2</mn><mo>,</mo><mi>l</mi></mrow><mi>l</mi></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><mi>Z</mi><mo>|</mo><msubsup><mo>|</mo><mrow><mn>2</mn><mo>,</mo><mi>l</mi></mrow><mi>l</mi></msubsup></mrow><mo>}</mo></mrow></mrow>]]></math><img file="FDA0000788880160000013.GIF" wi="478" he="115" /></maths><img file="FDA0000788880160000014.GIF" wi="344" he="162" />其中,l为范数参数,λ为正则参数,||·||<sub>2,l</sub>表示先对矩阵的各行向量求2范数,再将所得的2范数结果构成列向量求l范数;(4)根据步骤(3)中的双重稀疏约束函数,构造拉格朗日函数f(U,Z):<img file="FDA0000788880160000015.GIF" wi="1366" he="155" />其中β<sub>1</sub>为距离向惩罚因子,β<sub>2</sub>方位向为惩罚因子,||·||<sub>F</sub>表示求矩阵的Frobenius范数,λ为正则参数;(5)用交替迭代法对步骤(4)中的拉格朗日函数f(U,Z)进行求解,得到中间变量矩阵Z、方位向相位误差Φ<sub>1</sub>和距离向相位误差Φ<sub>2</sub>;(6)对步骤(5)得到中间变量矩阵Z的最终迭代结果取转置,得到稀疏目标场景的散射系数矩阵F=Z<sup>T</sup>。
地址 710071 陕西省西安市太白南路2号