发明名称 一种Capon波束形成器的加权稀疏约束稳健设计方法
摘要 本发明提供了一种Capon波束形成器的加权稀疏约束稳健设计方法,根据所采样的数据构建协方差矩阵,使用特征分解得到信号子空间和噪声子空间,利用信号子空间与噪声子空间的正交性构造加权矩阵;利用所构造的加权矩阵对SCCB中的稀疏约束项进行加权,得到WSCCB,获得比SCB和SCCB更高的输出SINR,以及更优的抗角度失配能力。
申请公布号 CN104931941A 申请公布日期 2015.09.23
申请号 CN201510324985.8 申请日期 2015.06.15
申请人 西北工业大学 发明人 孙超;刘振;刘雄厚;卓颉;杨益新
分类号 G01S7/36(2006.01)I 主分类号 G01S7/36(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 一种Capon波束形成器的加权稀疏约束稳健设计方法,其特征在于包括下述步骤:1)根据所采样的数据构建协方差矩阵,使用特征分解得到信号子空间和噪声子空间,利用信号子空间与噪声子空间的正交性构造加权矩阵;具体内容如下:利用阵列采集信号,得到数据采样矩阵X,利用X求得阵列采样数据的协方差矩阵<img file="FDA0000737696990000011.GIF" wi="339" he="127" />其中,L代表采样的快拍数;对R<sub>x</sub>进行特征分解,R<sub>x</sub>=UΓU<sup>H</sup>,其中,U和Γ分别为R<sub>x</sub>的特征向量和特征值矩阵;在入射信号和干扰互不相关的情形下,R<sub>x</sub>的特征值满足λ<sub>1</sub>≥λ<sub>2</sub>≥…≥λ<sub>J+1</sub>>λ<sub>J+2</sub>≥…≥λ<sub>M</sub>,其中,λ<sub>m</sub>为第m个特征值,m=1,2,…,M,J为干扰个数,M为阵元个数;M个特征值中,较大的J+1个特征值所对应的特征向量张成的子空间为信号子空间,记为U<sub>S</sub>;剩余特征值对应的特征向量张成的子空间为噪声子空间,记为U<sub>N</sub>;A<sup>H</sup>U<sub>N</sub>=0,其中,A为信号和干扰所在方位上的阵列流型向量所组成的矩阵,每一列对应着信号或某个干扰方位上的阵列流型向量;设干扰信号的阵列流型向量矩阵为A<sub>I</sub>,<img file="FDA0000737696990000016.GIF" wi="273" he="85" />设旁瓣区域阵列流型向量矩阵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>A</mi><mi>&Theta;</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>A</mi><mi>I</mi></msub></mtd><mtd><msub><mi>A</mi><msup><mi>I</mi><mi>C</mi></msup></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000737696990000017.GIF" wi="343" he="73" /></maths>其中<img file="FDA0000737696990000018.GIF" wi="83" he="73" />为A<sub>Θ</sub>中除A<sub>Ι</sub>之外的阵列流型向量组成的矩阵,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>A</mi><mi>&Theta;</mi><mi>H</mi></msubsup><msub><mi>U</mi><mi>N</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>A</mi><mi>I</mi><mi>H</mi></msubsup><msub><mi>U</mi><mi>N</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>A</mi><msup><mi>I</mi><mi>C</mi></msup><mi>H</mi></msubsup><msub><mi>U</mi><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>D</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>D</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mi>D</mi><mo>,</mo></mrow>]]></math><img file="FDA0000737696990000012.GIF" wi="735" he="197" /></maths>其中,D<sub>1</sub>为干扰阵列流型向量与噪声子空间的内积,其值均为0,D<sub>2</sub>为旁瓣区域中除干扰外的阵列流型向量与噪声子空间向量的内积,D为整个旁瓣区域阵列流型向量与噪声子空间的内积;对矩阵D的每一行取l<sub>2</sub>范数,得到列向量<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>d</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>d</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>d</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000737696990000013.GIF" wi="242" he="179" /></maths>其中,d<sub>1</sub>对应于干扰方位,其值均为0,d<sub>2</sub>对应于旁瓣区域中除干扰外的其它方位;最终,取加权矩阵<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><mi>diag</mi><mrow><mo>(</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn><mo>/</mo><msub><mi>d</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>1</mn><mo>/</mo><msub><mi>d</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>)</mo></mrow><mo>=</mo><mi>diag</mi><mrow><mo>(</mo><mn>1</mn><mo>/</mo><mi>d</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000737696990000014.GIF" wi="742" he="193" /></maths>其中,diag()表示利用括号内的向量生成对角阵,该对角阵主对角线元素为括号内的向量,其它元素均为0;2)利用所构造的加权矩阵对SCCB中的稀疏约束项进行加权,得到WSCCB的优化模型<img file="FDA0000737696990000015.GIF" wi="717" he="86" />s.t.w<sup>H</sup>a(θ<sub>0</sub>)=1,其中,||||<sub>1</sub>代表l<sub>1</sub>范数约束,γ≥0为稀疏约束项w<sup>H</sup>A<sub>Θ</sub>的加权系数,θ<sub>0</sub>代表波束图主瓣指向;求解获得WSCCB的加权向量w。
地址 710072 陕西省西安市友谊西路127号