发明名称 一种利用稀疏统计特性的机载雷达前视超分辨成像方法
摘要 本发明属于机载雷达前视超分辨成像技术领域,公开了一种利用稀疏统计特性的机载雷达前视超分辨成像方法。该利用稀疏统计特性的机载雷达前视超分辨成像方法包括以下步骤:得出距离徙动校正后数据s天线方向图矩阵H;机载雷达有M个距离单元;在距离徙动校正后数据s中,与第i距离单元对应的数据为s<sub>i</sub>,i取1至M;则将s<sub>i</sub>转化为如下形式:s<sub>i</sub>=H·σ<sub>i</sub>+n,其中,σ<sub>i</sub>表示与待求的第i距离单元对应的前视场景目标散射系数,n为设定的噪声;依据贝叶斯准则和目标场景的空域稀疏特性,并根据共轭梯度法得出σ<sub>1</sub>至σ<sub>M</sub>,完成机载雷达在前视场景下的超分辨成像过程。
申请公布号 CN103885058B 申请公布日期 2016.03.02
申请号 CN201410067314.3 申请日期 2014.02.26
申请人 西安电子科技大学 发明人 李明;陈洪猛;吴艳;张鹏;卢云龙;贾璐;张晓伟;刘高峰
分类号 G01S13/89(2006.01)I;G01S7/41(2006.01)I;G06F19/00(2011.01)I 主分类号 G01S13/89(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种利用稀疏统计特性的机载雷达前视超分辨成像方法,其特征在于,包括以下步骤:S1:机载雷达接收原始回波信号,对原始回波信号依次进行距离向脉冲压缩和距离徙动校正,得到距离徙动校正后数据s;根据实测的天线方向图构建天线方向图向量,通过对天线方向图向量进行循环移位,得到天线方向图矩阵H:H=[h<sub>‑Na/2</sub>,h<sub>‑Na/2+1</sub>...,h<sub>0</sub>,...,h<sub>Na/2</sub>]其中,h<sub>0</sub>为天线方向图向量,h<sub>k</sub>表示对h<sub>0</sub>进行循环移位,k取‑Na/2至Na/2且k≠0,Na为机载雷达方位向的锐化比;k的正负用于表示h<sub>0</sub>的两个相反的移位方向;S2:机载雷达有M个距离单元;在距离徙动校正后数据s中,与第i距离单元对应的数据为s<sub>i</sub>,i取1至M;则将s<sub>i</sub>转化为如下形式:s<sub>i</sub>=H·σ<sub>i</sub>+n其中,σ<sub>i</sub>表示与待求的第i距离单元对应的前视场景目标散射系数,n为设定的噪声;S3:依据贝叶斯准则和目标场景的空域稀疏特性,将σ<sub>i</sub>的求解问题转化为以下优化模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>=</mo><munder><mi>argmax</mi><msub><mi>&sigma;</mi><mi>i</mi></msub></munder><mo>&lsqb;</mo><mi>P</mi><mi>D</mi><mi>F</mi><mrow><mo>(</mo><msub><mi>s</mi><mi>i</mi></msub><mo>|</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>,</mo><msup><mi>&epsiv;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>P</mi><mi>D</mi><mi>F</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>|</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000806988610000013.GIF" wi="988" he="126" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>P</mi><mi>D</mi><mi>F</mi><mrow><mo>(</mo><msub><mi>s</mi><mi>i</mi></msub><mo>|</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>,</mo><msup><mi>&epsiv;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mn>2</mn><msup><mi>&pi;&epsiv;</mi><mn>2</mn></msup><mo>)</mo></mrow><mrow><mo>(</mo><mo>-</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow></msup><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msup><mi>&epsiv;</mi><mn>2</mn></msup></mrow></mfrac><mo>|</mo><mo>|</mo><msub><mi>s</mi><mi>i</mi></msub><mo>-</mo><msub><mi>H&sigma;</mi><mi>i</mi></msub><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000806988610000011.GIF" wi="1180" he="138" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>P</mi><mi>D</mi><mi>F</mi><mrow><mo>(</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>|</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><mfrac><mi>&gamma;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&gamma;</mi><mo>|</mo><mo>|</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>(</mo><mi>j</mi><mo>)</mo><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>&gamma;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mi>N</mi></msup><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>&gamma;</mi><mo>|</mo><mo>|</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000806988610000012.GIF" wi="1326" he="165" /></maths>其中,γ为拉普拉斯系数,ε<sup>2</sup>为设定的噪声的方差,σ<sub>i</sub>(j)表示第i距离单元第j个散射点的散射系数,N为散射点的个数;||·||<sub>1</sub>表示求L1范数,||·||<sub>2</sub>表示求L2范数;PDF(s<sub>i</sub>|σ<sub>i</sub>,ε<sup>2</sup>)表示回波信号的似然函数,PDF(σ<sub>i</sub>|γ)描述前视场景下目标散射系数的稀疏统计特性;S4:利用共轭梯度法对步骤S3中的优化模型进行求解,得出σ<sub>i</sub>;根据求解得出的σ<sub>1</sub>至σ<sub>M</sub>,完成机载雷达在前视场景下的超分辨成像过程。
地址 710071 陕西省西安市太白南路2号