发明名称 基于稀疏表示的雷达稳健自适应波束形成方法
摘要 本发明公开了一种基于稀疏表示的雷达稳健自适应波束形成方法,其主要思路为:分别设定雷达的发射信号波长λ和雷达阵元间距d,进而获取L次快拍的雷达回波信号,并计算得到雷达回波信号的协方差矩阵<img file="DDA0001028091220000011.GIF" wi="67" he="62" />设定目标所在角度空域范围为[θ<sub>d</sub>,θ<sub>u</sub>],并将该角度空域范围[θ<sub>d</sub>,θ<sub>u</sub>]划分为K个角度值,分别计算目标所在角度空域范围内K个角度值的常规波束形成器W<sub>c</sub>和目标所在角度空域范围内K个角度值的直接矩阵求逆波束形成器W<sub>s</sub>,进而依次计算得到雷达回波信号中的噪声功率<img file="DDA0001028091220000012.GIF" wi="57" he="71" />和雷达回波信号中的干扰信号子空间P<sub>I</sub>,并根据W<sub>c</sub>和W<sub>s</sub>,依次计算得到雷达回波信号的基矩阵B和雷达回波信号的期望波束形成器w<sub>d</sub>,进而计算雷达回波信号中的稀疏表示矢量<img file="DDA0001028091220000013.GIF" wi="54" he="53" />并根据B计算雷达回波信号的最终波束形成器<img file="DDA0001028091220000014.GIF" wi="61" he="63" />。
申请公布号 CN106093887A 申请公布日期 2016.11.09
申请号 CN201610473214.X 申请日期 2016.06.24
申请人 西安电子科技大学 发明人 赵永波;张盼盼;程增飞;水鹏朗;李慧;何学辉
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于稀疏表示的雷达稳健自适应波束形成方法,其特征在于,包括以下步骤:步骤1,分别设定雷达的发射信号波长λ和雷达阵元间距d,进而获取L次快拍的雷达回波信号,其中第l次快拍的雷达回波信号为x<sub>l</sub>,l∈{1,2,…,L},N为雷达阵元个数,L为雷达回波信号包含的快拍次数;步骤2,利用L次快拍的雷达回波信号,计算得到雷达回波信号的协方差矩阵<img file="FDA0001028091190000016.GIF" wi="67" he="64" />步骤3,设定目标所在角度空域范围为[θ<sub>d</sub>,θ<sub>u</sub>],并将该角度空域范围[θ<sub>d</sub>,θ<sub>u</sub>]划分为K个角度值,其中第k个角度值为θ<sub>k</sub>,然后根据雷达的发射信号波长λ和雷达阵元间距d,分别计算得到目标所在角度空域范围内K个角度值的常规波束形成器W<sub>c</sub>和目标所在角度空域范围内K个角度值的直接矩阵求逆波束形成器W<sub>s</sub>,其中第k个角度值θ<sub>k</sub>对应的常规波束形成器为w<sub>ck</sub>,第k个角度值θ<sub>k</sub>对应的直接矩阵求逆波束形成器为w<sub>sk</sub>;k∈{1,2,…,K},K为目标所在角度空域范围包含的角度值个数,θ<sub>d</sub>为设定的目标所在角度空域范围的下限,θ<sub>u</sub>为设定的目标所在角度空域范围的上限;步骤4,根据目标所在角度空域范围为[θ<sub>d</sub>,θ<sub>u</sub>]和雷达回波信号的协方差矩阵<img file="FDA0001028091190000011.GIF" wi="67" he="63" />计算得到雷达回波信号中干扰信号的入射方向,其中第m个干扰信号入射方向为<img file="FDA0001028091190000012.GIF" wi="90" he="70" />并计算得到雷达回波信号中的噪声功率<img file="FDA0001028091190000013.GIF" wi="83" he="62" />步骤5,根据雷达回波信号中干扰信号的入射方向,计算得到雷达回波信号中的干扰信号子空间P<sub>I</sub>;步骤6,根据目标所在角度空域范围内K个角度值的常规波束形成器W<sub>c</sub>和目标所在角度空域范围内K个角度值的直接矩阵求逆波束形成器W<sub>s</sub>,依次计算得到雷达回波信号的基矩阵B和雷达回波信号的期望波束形成器w<sub>d</sub>;步骤7,根据雷达回波信号中的干扰信号子空间P<sub>I</sub>、雷达回波信号中的噪声功率<img file="FDA0001028091190000014.GIF" wi="59" he="71" />和雷达回波信号中待求的稀疏表示矢量u,计算得到雷达回波信号的稀疏模型,其表达式为:<maths num="0001"><math><![CDATA[<mrow><munder><mrow><mi>min</mi><mi>i</mi><mi>m</mi><mi>i</mi><mi>z</mi><mi>e</mi></mrow><mi>u</mi></munder><mo>|</mo><mo>|</mo><msup><mi>u</mi><mi>H</mi></msup><msub><mover><mi>R</mi><mo>^</mo></mover><mi>y</mi></msub><mi>u</mi><mo>-</mo><msubsup><mi>&sigma;</mi><mi>u</mi><mn>2</mn></msubsup><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow>]]></math><img file="FDA0001028091190000015.GIF" wi="494" he="95" /></maths>subject to||P<sub>I</sub>Bu||<sub>2</sub>≤ξ,||u||<sub>1</sub>≤δ满足所述稀疏模型条件的雷达回波信号中待求的稀疏表示矢量u为雷达回波信号中的稀疏表示矢量<img file="FDA0001028091190000021.GIF" wi="53" he="53" />其中,<img file="FDA0001028091190000022.GIF" wi="270" he="77" />B=[W<sub>c</sub>,W<sub>s</sub>],<img file="FDA0001028091190000023.GIF" wi="218" he="62" /><img file="FDA0001028091190000024.GIF" wi="61" he="63" />为雷达回波信号中的噪声功率,||·||<sub>2</sub>为二范数,||·||<sub>1</sub>为一范数,α为人为确定的大于零的系数;ξ为约束雷达回波信号中的干扰信号子空间P<sub>I</sub>和雷达回波信号中待求的稀疏表示矢量u之间正交性的参数,δ为约束雷达回波信号中待求的稀疏表示矢量u稀疏性的参数,且α、ξ和δ均大于零,P<sub>I</sub>为雷达回波信号中的干扰信号子空间,B为雷达回波信号的基矩阵;步骤8,根据雷达回波信号中的稀疏表示矢量<img file="FDA0001028091190000025.GIF" wi="35" he="46" />和雷达回波信号的基矩阵B,计算得到雷达回波信号的最终波束形成器<img file="FDA0001028091190000026.GIF" wi="91" he="62" />
地址 710071 陕西省西安市太白南路2号