一种适用于圆形阵列的高增益波束形成方法技术

利用大规模水下阵列进行目标观测是水下警戒的有效手段。一般而言，我们希望通过增加阵列孔径提升阵列的观测能力,但是，阵列的孔径受到声场相干性的限制。在水平方向上，声场相关性主要分为水平纵向相关性和水平横向相关性。水平横向相关性的影响因素主要是介质的非均匀特性，水平纵向相关性的主要影响因素是简正波的干涉。其中，水平纵向相关性的下降会对阵列的有效孔径造成影响较大。本发明专利技术提出了一种适用于水下均匀圆形水平阵列的波束形成处理方法，该方法不需要测量额外的环境参数，对信号的频谱也不需要做过多假设，仅利用均匀圆阵的对称特性，将阵列分成沿水平纵向对称的两个子阵，通过分裂阵互相关波束形成算法，便可以有效消除声场纵向相关性下降带来的不利影响，提升阵增益，从而提高阵列的观测能力。从而提高阵列的观测能力。从而提高阵列的观测能力。

【技术实现步骤摘要】
一种适用于圆形阵列的高增益波束形成方法


[0001]本专利技术属于阵列信号处理领域，具体涉及一种适用于圆形阵列的高增益波束形成方法。

技术介绍

[0002]利用大规模水下阵列进行目标观测是水下警戒的有效手段。一般而言，我们希望通过增加阵列孔径提升阵列的观测能力,但是，阵列的孔径受到声场相干性的限制。在水平方向上，声场相关性主要分为水平纵向相关性和水平横向相关性。水平横向相关性的影响因素主要是介质的非均匀特性，水平纵向相关性的主要影响因素是简正波的干涉。其中，水平纵向相关性的下降会对阵列的有效孔径造成影响较大。针对该问题，苏小星等人利用声场的波导不变量提出了频移补偿方法(苏晓星,张仁和,李风华.利用波导不变性提高声场的水平纵向相关[J].声学学报,2006(4):305
‑
309)。论文作者根据声场的波导不变性，指出了声场频域响应与水平纵向距离之间的线性相位关系，然后在频谱缓变的假设下，提出了通过频移补偿提高水平纵向相关性的方法。该方法需要有相对精确的波导不变量的先验值或者测量值，以确定频移补偿量。但是实际应用中，水文环境参数随时间、季节的变化剧烈，精确测量和预测水文参数的难度较大；另外，水下目标辐射噪声一般具有线谱特征，很难满足频谱缓变的假设，这些问题限制了该算法的广泛应用。

技术实现思路

[0003]本专利技术提出了一种适用于水下均匀圆形水平阵列的波束形成处理方法，该方法不需要测量额外的环境参数，对信号的频谱也不需要做过多假设，仅利用均匀圆阵的对称特性，将阵列分成沿水平纵向对称的两个子阵，通过分裂阵互相关波束形成算法，便可以有效消除声场纵向相关性下降带来的不利影响，提升阵增益，从而提高阵列的观测能力。具体包括以下步骤：
[0004]步骤1：从圆形阵列中选取最优的阵段的中心阵元
[0005]将N元均匀圆形阵列的阵形用极坐标表示为：
[0006][0007]其中，R为圆形阵列的极径，γ
n
为第n号阵元的极角，n＝1...N其中，R为圆形阵列的极径，γ
n
为第n号阵元的极角，n＝1...N，i表示虚数单位；
[0008]设点声源信号的频谱为S(ω)，通过N元均匀圆形阵列接收信号，则接收信号表示为：
[0009][0010]其中H(ω,r
n
)是由声源到第n号阵元的信道传输函数，r
n
表示收发距离，对于点源声场，由简正波理论可得：
[0011][0012]其中：
[0013][0014]在式(3)、式(4)中，ω为声源角频率，z
s
为声源深度，z
r
为接收深度，r为传播距离，μ
l
、β
l
和ψ
l
分别为第l号简正波的本征值、衰减系数和本征函数，μ
m
是m号主要简正波的本征值；
[0015]进行最有波束形成阵列段选取的策略是使阵元在声传播的水平纵向上的分布距离差异最小，当扫描角度为θ时，按照上述原则，最优阵列的中心阵元可以通过下式确定：
[0016][0017]步骤2：将阵段分成两个子阵，并构造两个子阵的导向矢量
[0018]以n0阵元为中心选取M个阵元，M为奇数，并分为两个等长的子阵，两个子阵的序号分别为：
[0019]N
sub1
＝[n0,n0‑
1,
…
,n0‑
(M
‑
1)/2][0020]N
sub2
＝[n0,n0+1,
…
,n0+(M
‑
1)/2]ꢀꢀ
(6)
[0021]相应的接收信号分别为S
array1
＝{S
n
(ω)|n∈N
sub1
}，S
array2
＝{S
n
(ω)|n∈N
sub2
}；
[0022]两个子阵波束形成的导向矢量分别表示为：
[0023][0024]导向矢量中的元素a
n
(θ,ω)为：
[0025][0026]其中p
n
＝(x
n
,y
n
)为阵元n的直角坐标c
m
为参考声速，一般取主要简正波对应的群速度。
[0027]步骤3：分别对两个子阵进行波束形成
[0028]根据步骤2得出的两个子阵导向矢量进行波束形成，第一个子阵的波束输出为：
[0029][0030]同理，第二个子阵的波束输出为
[0031][0032]步骤4：计算两个子阵输出信号的互相关函数，输出方位谱
[0033]计算两波束输出信号的互相关函数，可直接在频域共轭相乘，然后反变换到时域：
[0034][0035]表示傅里叶反变换；
[0036]将输出波束(9)、(10)看作以中心阵元n0进行的时延补偿，因此直接输出时延为0时刻的互相关值作为当前扫描方位θ的方位谱强度：
[0037]bf
out
(θ)＝|c
b1,b2
(0)| (12)
[0038]最后，改变θ的值并重复步骤1～步骤4得到方位谱。
[0039]本专利技术的有益效果在于：
[0040]依据本专利技术方法的声场水平纵向相关半径比常规处理方法提升了2.5倍，对声场导致的相干损失有5dB左右的抑制能力，输出信噪比增益显著提高。
附图说明
[0041]图1为通道间传输函数的相关系数图。
[0042]图2为阵元选择示意图。
[0043]图3为常规处理方法与分裂互相关波束形成方法波束图。
[0044]图4(a)为无声场时波束输出增益与阵元数量关系图。
[0045]图4(b)为有声场时波束输出增益与阵元数量关系图。
[0046]图5(a)为有声场时125阵元输出时方位谱。
[0047]图5(b)为无声场时125阵元输出时方位谱。
[0048]图6(a)为常规处理方法的方位历程图。
[0049]图6(b)为本专利技术所提方法的方位历程图。
具体实施方式
[0050]下面结合符合和具体实施方式，对本专利技术做进一步详细说明。
[0051]仿真中，海水深度100m，声速梯度剖面为典型负跃层。接收阵为均匀圆阵布设在海底，阵元间距为3.75m，共有阵元250个。噪声源深度设置为35m，噪声源距离接收阵10km，位于阵列的90
°
方向上。仿真频段为100～200Hz。利用Kraken声场计算程序生成浅海声场数据。
[0052]首先分析通道间传输函数的相关性，如图1所示。从仿真结果中可以看出，沿入射方向，随着阵元之间的水平纵向距离增加，通道间信号的相关性也随之下降，导致波束形成增益损失。
[0053]因此在工程应用中会选取部分相关性较高的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种适用于圆形阵列的高增益波束形成方法，其特征在于：具体包括以下步骤：步骤1：从圆形阵列中选取最优的阵段的中心阵元将N元均匀圆形阵列的阵形用极坐标表示为：其中，R为圆形阵列的极径，γ
n
为第n号阵元的极角，n＝1...N，i表示虚数单位；设点声源信号的频谱为S(ω)，通过N元均匀圆形阵列接收信号，则接收信号表示为：其中H(ω,r
n
)是由声源到第n号阵元的信道传输函数，r
n
表示收发距离，对于点源声场，由简正波理论可得：其中在式(3)、式(4)中，ω为声源角频率，z
s
为声源深度，z
r
为接收深度，r为传播距离，μ
l
、β
l
和ψ
l
分别为第l号简正波的本征值、衰减系数和本征函数，μ
m
是m号主要简正波的本征值；进行最有波束形成阵列段选取的策略是使阵元在声传播的水平纵向上的分布距离差异最小，当扫描角度为θ时，按照上述原则，最优阵列的中心阵元可以通过下式确定：步骤2：将阵段分成两个子阵，并构造两个子阵的导向矢量以n0阵元为中心选取M个阵元，M为奇数，并分为两个等长的子阵，两个子阵的序号分别为：N
sub1
＝[n0,n0‑
1,
…
,n0‑
(M
‑
1)/2]N
sub2
＝[n0,n0+1,
…
,n0+(M
‑
1)/2]
ꢀꢀꢀꢀ

【专利技术属性】
技术研发人员：张蒙，陈亚伟，
申请(专利权)人：中国电子科技集团公司第十四研究所，
类型：发明
国别省市：