主权项 |
一种基于干扰噪声协方差矩阵重构的波束形成方法,其特征在于,包括如下步骤: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>θ</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>Σ</mi><mrow><mi>d</mi><mo>=</mo><mn>2</mn></mrow><mi>D</mi></munderover><mi>a</mi><mrow><mo>(</mo><msub><mi>θ</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>Σ</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>。 |