发明名称 一种基于干扰噪声协方差矩阵重构的波束形成方法
摘要 本发明属于阵列信号处理领域,主要涉及基于最差性能最优化准则的Capon自适应波束形成算法的稳健性。本发明提供一种基于干扰噪声协方差矩阵重构的波束形成方法,基于最大化最差情况下的波束输出SINR的目标,首先采用阵列接收数据对干扰噪声协方差矩阵进行重新估计<img file="DDA0000641240970000011.GIF" wi="111" he="64" />然后用<img file="DDA0000641240970000012.GIF" wi="87" he="64" />代替原始最差性能最佳化算法中的样本协方差矩阵<img file="DDA0000641240970000014.GIF" wi="77" he="60" />按照WCPO准则使得最差情况下的波束输出OSINR最大化。本发明用重构的干扰噪声协方差矩阵重构<img file="DDA0000641240970000013.GIF" wi="83" he="64" />代替原始WCPO优化问题中的样本协方差矩阵<img file="DDA0000641240970000015.GIF" wi="76" he="60" />可以大大减弱或避免期望信号自零陷现象,最大化最差情况的输出SINR,提高波束形成算法的稳健性。
申请公布号 CN104535969A 申请公布日期 2015.04.22
申请号 CN201410809621.4 申请日期 2014.12.23
申请人 电子科技大学 发明人 袁晓垒;甘露;杨金金;廖红舒;张花国
分类号 G01S7/285(2006.01)I;G01S7/35(2006.01)I;G01S7/38(2006.01)I 主分类号 G01S7/285(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李玉兴
主权项 一种基于干扰噪声协方差矩阵重构的波束形成方法,其特征在于,包括如下步骤:S1、由M个阵元构成的均匀线阵接收到D个来自远场信源的信号,所述D个信号的来波方向为θ<sub>d</sub>,设第1个信号为期望信号,其余D‑1个均为干扰信号,所述D个信号之间互不相关,且信号与噪声之间互不相关,则第n个快拍下阵列接收数据记为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>d</mi><mo>=</mo><mn>2</mn></mrow><mi>D</mi></munderover><mi>a</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mo>)</mo></mrow><msub><mi>s</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>As</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000641240940000011.GIF" wi="1143" he="141" /></maths>其中,A=[a(θ<sub>1</sub>),…,a(θ<sub>D</sub>)]为阵列流型矩阵,s(n)为阵列接收到的信号源矢量,v(n)表示阵列接收到的噪声矢量,v(n)为零均值高斯白噪声,则阵列接收到的N个快拍数据可表示为如下的形式:X=[x(1),…,x(N)]=AS+VS=[s(1),…,s(N)]V=[v(1),…,v(N)]由阵列接收数据矩阵X可以得到阵列接收数据的样本协方差矩阵<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><msup><mi>XX</mi><mi>H</mi></msup><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>x</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000641240940000012.GIF" wi="681" he="141" /></maths>设期望信号的预估计导向矢量为<img file="FDA0000641240940000013.GIF" wi="152" he="90" />真实的信号导向矢量a(θ<sub>1</sub>)位于椭圆不确定集合<img file="FDA0000641240940000014.GIF" wi="801" he="101" />中,其中,ε表示期望信号预估计导向矢量<img file="FDA0000641240940000015.GIF" wi="122" he="84" />与真实导向矢量a(θ<sub>1</sub>)之间估计误差δ的范数上界,其中,M为整数,D为整数,N为整数,d=1,2,3,…,D,‑90°≤θ<sub>d</sub>≤90°,n=1,2,3,...,N;S2、将信号来波方向的整个区间Θ=[‑90°,90°]划分为两个区间Θ<sub>1</sub>,Θ<sub>2</sub>,所述两个区间Θ<sub>1</sub>,Θ<sub>2</sub>满足<img file="FDA0000641240940000016.GIF" wi="564" he="73" />其中,期望信号来波方向在区间Θ<sub>1</sub>中,即θ<sub>1</sub>∈Θ<sub>1</sub>,而干扰信号来波方向均在区间Θ<sub>2</sub>中,即θ<sub>d</sub>∈Θ<sub>2</sub>,d=2,…,D,对区间Θ<sub>2</sub>进行离散化得到离散角度区间<img file="FDA0000641240940000017.GIF" wi="110" he="76" />在<img file="FDA0000641240940000018.GIF" wi="66" he="77" />采用Capon空间谱估计算法进行干扰噪声协方差矩阵的重构,得到<img file="FDA0000641240940000019.GIF" wi="124" he="92" />S3、利用S2所述<img file="FDA00006412409400000110.GIF" wi="120" he="85" />S1所述<img file="FDA00006412409400000111.GIF" wi="128" he="83" />和a(θ<sub>1</sub>)所属的椭圆不确定集合S,采用最差性能最佳化准则进行WCPO波束形成设计,得到稳健的阵列加权w<sub>INCMR</sub>。
地址 611731 四川省成都市高新区(西区)西源大道2006号