主权项 |
基于双重稀疏约束的逆合成孔径雷达自聚焦成像方法,包括如下步骤:(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>[</mo><mtable><mtr><mtd><msup><mi>e</mi><mrow><msub><mi>jΦ</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Φ</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Φ</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></msup></mtd></mtr></mtable><mo>]</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Φ</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Φ</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>Φ</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>λ</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>。 |