一种双层圆环阵稳健超指向性波束形成方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN107238829A(43)申请公布日2017.10.10(21)申请号201710479784.4(22)申请日2017.06.22(71)申请人西北工业大学地址710072陕西省西安市友谊西路127号(72)发明人杨益新汪勇马远良雷波(74)专利代理机构西北工业大学专利中心61204代理人王鲜凯(51)Int.Cl.G01S7/52(2006.01)权利要求书4页说明书9页附图3页(54)发明名称一种双层圆环阵稳健超指向性波束形成方法(57)摘要本发明涉及一种双层圆环阵稳健超指向性波束形成方法，首先将阵列流形和空间相关矩阵表示为两个圆环阵各自阵列流形和空间相关矩阵的组合形式，然后利用分块矩阵求逆公式得到空间相关矩阵逆矩阵，并结合循环矩阵性质将最优权值向量分解为有限阶子分量叠加的形式，由此进一步将最优超指向性波束、最大指向性因子和总误差敏感度函数分别分解为各阶特征波束及其指向性因子和误差敏感度函数叠加的形式，其中随阶数升高对应的特征波束的误差敏感度函数变大，稳健性变差。依据实际情况确定适用的最高阶数，通过舍去对误差敏感的高阶项并保留稳健的低阶项，合成得到最终的超指向性结果。CN107238829ACN107238829A权利要求书1/4页1.一种双层圆环阵稳健超指向性波束形成方法，其特征在于步骤如下：步骤1：以双层圆环阵的阵列流形向量P表示两层圆环阵各自阵列流形：T所述P1＝[p1,0,p1,1,...,p1,M-1]为第1层圆环阵阵列流形；T所述P2＝[p2,0,p2,1,...,p2,M-1]为第2层圆环阵阵列流形；其中：k＝-k[sinθcosφ,sinθsinφ,cosθ]T，Tr1,m＝[a1sinθmcosφm,a1sinθmsinφm,a1cosθm]，r2,m＝[a2sinθmcos(φm+δ),a2sinθmsin(φmT+δ),a2cosθm+h]，a1和a2分别表示第1和第2层圆环阵的半径，θm和φm分别表示第m个阵元所在位置的垂直俯仰角和水平方位角，θ和φ分别表示平面波入射的垂直俯仰角和水平方位角，δ为第1和第2层圆环阵相同序号阵元所在位置方位角的差值，h为第2层圆环阵在z轴上的坐标值；波数k＝2π/λ，λ为波长，φm＝mβ，β＝2π/M，下标“1”和“2”分别表示第1和第2层圆环阵，下标“m”表示第m个阵元，两层圆环阵的阵元数均为M，双层圆环阵总阵元数为2M，(·)T表示转置；以空间相关矩阵ρn表示空间相关矩阵的组合形式为：所述矩阵ρ11和ρ22分别是第1和第2层圆环阵各自的空间相关矩阵，其维度均为M×M，元素依次为：ρ11,m,n＝ρ11,s＝sinc(k·d1s)ρ22,m,n＝ρ22,s＝sinc(k·d2s)其中：d1s＝2a1sin(sβ/2)，d2s＝2a2sin(sβ/2)，s＝|m-n|；所述矩阵ρ12和ρ21是由第1层和第2层圆环阵相互之间的空间相关矩阵，其维度均为M×M，元素为：ρ12,m,n＝ρ12,s＝ρ21,m,n＝ρ21,s＝sinc(k·d3s)其中：所述矩阵ρ11、ρ22、ρ12和ρ21均是循环矩阵，其特征值分别为：以上特征值都是实数且满足对称关系λm＝λM-m；步骤2：利用分块矩阵求逆公式将空间相关矩阵ρn的逆矩阵表示为如下形式：2CN107238829A权利要求书2/4页其中：步骤3：计算最优权值向量其中：-1/2imβi(M-1)mβTvm＝M[1，e,...,e]H归一化参数(θ0,φ0)为设定的波束指向方向，(·)表示共轭转置，(·)*表示求共轭；步骤4：将最优权值向量代入波束形成公式得到超指向性波束，并进一步将其分解为：其中：3CN107238829A权利要求书3/4页最优指向性因子为：误差敏感度函数为：其中：和分别为第m阶模态波束及其指向性因子和误差敏感度函数；步骤5：将模态波束及其指向性因子和误差敏感度函数进行组合叠加：其中阵元数M为奇数；当阵元数M为偶数时，组合叠加形式为：4CN107238829A权利要求书4/4页其中Bm、Dm和Tm分别为第m阶特征波束及其指向性因子和误差敏感度函数；步骤6：舍去大于最大阶数N的特征波束，由公式和分别得到最终超指向性波束及其指向性因子和误差敏感度函数，其中α＝1/DF。5CN107238829A说明书1/9页一种双层圆环阵稳健超指向性波束形成方法技术领域[0001]本发明属于声学阵列信号处理和声呐技术等领域，涉及波束形成方法，特别是涉及一种双层圆环阵稳健超指向性波束形成方法，特别是涉及一种稳健超指向性波束形成方法，适用于双层圆环形阵列的低信噪比目标检测以及目标方位的高分辨率估计，。背景技术[0002]超指向性波束形成方法可以在不改变阵列尺