空域反转阵导向最小方差波束形成方法及装置制造方法及图纸

本发明专利技术公开了一种空域反转阵导向最小方差波束形成方法及装置，所述方法包括：将原阵列快拍与自身进行空域反转过后的阵列快拍进行卷积，再将卷积后每个阵元的功率进行均衡，得到空域反转阵列的快拍；获取空域反转阵的宽带阵列信号和方向向量，确定波束形成的协方差矩阵，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱。本发明专利技术通过将空间反转和卷积，得到孔径扩大近一倍的空域反转阵，并构造空域反转阵的导向向量，实现其导向最小方差波束形成，该方法对目标辐射信号的统计特征、频谱特征没有限制，并且对阵形没有特殊要求。并且对阵形没有特殊要求。并且对阵形没有特殊要求。

【技术实现步骤摘要】
空域反转阵导向最小方差波束形成方法及装置


[0001]本专利技术涉及水声阵列领域，具体涉及一种空域反转阵导向最小方差波束形成方法及装置。

技术介绍

[0002]对于被动水声阵列而言，阵列孔径越大则阵列性能越好。水声环境的时变空变、低信噪比等特点使得自适应波束形成成为目前水声阵列信号处理研究和应用最好的方法。但受尺寸和动力的限制，水下无人潜航器不能拖曳大孔径的水声线阵列，从而导致其阵列性能受到制约。如何提高水下无人潜航器小尺寸线阵列的性能成为目前研究的需求和热点。
[0003]自适应波束形成是利用一定的准则来设计阵列数据的权矢量来满足波束形成的要求，同时抑制干扰信号的噪声的方法。最小方差无畸变波束形成(minimum variance distortionless response，MVDR)在阵列模型和接收数据的协方差矩阵确定时被证明是最优的，但在宽带水声条件下，当目标运动速度太快时，MVDR收敛所需快拍数量不能小于阵元数量，因此收敛时间长，导致对运动目标的方位估计存在偏差。导向最小方差(steered minimum variance，STMV)波束形成通过多个频域子带的相干叠加，增大协方差矩阵的自由度，提高算法收敛速度，因此需要的时域快拍数量远少于MVDR，对于快速运动目标的方位谱估计效果更准确。
[0004]对于小尺寸阵列，阵列接收的数据相对于阵列的数量变得充沛，使得协方差矩阵获得更高的估计精确性。因此影响小尺寸阵列性能的因素就从接受数据协方差矩阵的误差转变为由于阵列孔径的减小导致的宽主瓣。主瓣宽度变宽会导致两个后果：一是造成阵列的多目标分辨性能下降；二会让强度高的目标成为强的干扰，导致附近相对弱的目标难以分辨。
[0005]为解决阵列孔径不足导致主瓣宽度变宽的问题，发展了阵列扩展技术。现有的阵列扩展主要有以下五种：1、四阶累积量法，四阶累积量在阵列扩展的基础上还能够抑制高斯噪声，但水面水下目标辐射的声信号并不符合四阶累积量法的条件，即信号的四阶累积量为零；2、内插变换法，其阵列扩展的原理是通过在特定的扫描子区域插值来拟合出实际阵元与虚拟阵元导向向量间的特定映射关系，缺点是当目标与变换区域的位置误差过大时会导致扩展阵列的波束形成出现零限漂移现象，即存在“角度敏感”问题；3、矩阵重构法，其利用某些阵列的特殊结构达到阵列扩展的目的，但不适用于水声线阵列。4、合成孔径法，其利用目标的时间相关性，只适用于目标信号存在线谱的情况，同样不适用于水声信号；5、互质阵法，需要两排阵列和特定的阵列间距。因此，现有的阵列扩展方法因为对目标特性和阵型的要求从而不适用于水下均匀线阵列的应用环境。

技术实现思路

[0006]本专利技术的目的是提出一种空域反转阵导向最小方差波束形成方法及装置，至少部分地解决现有技术中的问题。
[0007]为了实现上述专利技术目的，本专利技术采用如下的技术方案：
[0008]一种空域反转阵导向最小方差空间波束形成方法，包括以下步骤：
[0009]将原阵列快拍与自身进行空域反转过后的阵列快拍进行卷积，再将卷积后每个阵元的功率进行均衡，得到空域反转阵列的快拍；
[0010]获取空域反转阵的宽带阵列信号和方向向量，确定波束形成的协方差矩阵，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱。
[0011]进一步地，所述空域反转阵列的快拍表示为：
[0012][0013]其中M为原阵列阵元数，N是空域反转阵中阵元的序号，N＝
‑
M+1,
…
M
‑
1，表示原阵列中第m+(|N|
‑
N)/2个阵元输出信号的共轭，X
(m+(N+N)/2)
(t)表示原阵列中第m+(|N|
‑
N)/2个阵元输出的信号。
[0014]进一步地，获取空域反转阵的宽带阵列信号和方向向量，确定波束形成的协方差矩阵包括：
[0015]对空域反转阵列的快拍做傅里叶变换得到空域反转阵的宽带阵列信号X
s
(f
j
)，根据阵型得到阵列的方向向量a
s
(f
j
,θ)，基于下式确定波束形成的协方差矩阵：
[0016][0017]其中：
[0018][0019]R
a
(f
j
)＝E{X
s
(f
j
)X
sH
(f
j
)}
[0020]f
j
表示宽带信号带内第j个频点，J为宽带信号频点数，θ表示方向角。
[0021]进一步地，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱包括：
[0022]基于协方差矩阵以及权向量将波束形成转换为如下问题：
[0023][0024][0025]其中表示权向量，1
s
是(2M
‑
1)
×
1的1向量，表示寻找使标量最小时的方向角θ；
[0026]求解问题得到方向角θ，则有权向量为：方位谱为：
[0027]一种空域反转阵导向最小方差空间波束形成装置，包括
[0028]空域反转阵列构建模块，用于将原阵列快拍与自身进行空域反转过后的阵列快拍进行卷积，再将卷积后每个阵元的功率进行均衡，得到空域反转阵列的快拍；
[0029]波束形成处理模块，用于获取空域反转阵的宽带阵列信号和方向向量，确定波束形成的协方差矩阵，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱。
spatial reversal array。该方法对目标辐射信号的统计特征、频谱特征没有限制，并且对阵形没有特殊要求，适用于最普遍的均匀直线阵。同时也不存在“角度敏感”问题，可以广泛适用于水下均匀线阵的应用。
附图说明
[0052]图1为空域反转阵示意图；
[0053]图2为STMV和SRVA
‑
STMV两种方法的方位谱；
[0054]图3为STMV和SRVA
‑
STMV两种方法的目标分辨概率与信噪比的关系；
[0055]图4为STMV和SRVA
‑
STMV两种方法的输出信干噪比随输入信噪比的变化曲线；
[0056]图5为STMV和SRVA
‑
STMV两种方法的方位历程图；
[0057]图6为STMV和SRVA
‑
STMV两种方法在不同时刻的方位谱对比图。
具体实施方式
[0058]下面结合附图和实施例对本专利技术的技术方案作进一步说明。
[0059]考虑一个远场条件下的M元均匀线阵列，假设信号的入射角为θ，阵元间距为d。设入射信号为窄带信号，每个阵元接收的噪声也是互不相关的高斯白噪声，信号与噪声同样相互独立，阵列信号模型为：
[0060][0061]式(1)中，X(t)是t时刻的阵列快拍；s本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种空域反转阵导向最小方差波束形成方法，其特征在于，所述方法包括：将原阵列快拍与自身进行空域反转过后的阵列快拍进行卷积，再将卷积后每个阵元的功率进行均衡，得到空域反转阵列的快拍；获取空域反转阵的宽带阵列信号和方向向量，确定波束形成的协方差矩阵，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱。2.根据权利要求1所述的方法，其特征在于，所述空域反转阵列的快拍表示为：其中M为原阵列阵元数，N是空域反转阵中阵元的序号，N＝
‑
M+1,
…
M
‑
1，表示原阵列中第m+(|N|
‑
N)/2个阵元输出信号的共轭，X
(m+(N+N)/2)
(t)表示原阵列中第m+(|N|
‑
N)/2个阵元输出的信号。3.根据权利要求1所述的方法，其特征在于，获取空域反转阵的宽带阵列信号和方向向量，确定波束形成的协方差矩阵包括：对空域反转阵列的快拍做傅里叶变换得到空域反转阵的宽带阵列信号X
s
(f
j
)，根据阵型得到阵列的方向向量a
s
(f
j
,θ)，基于下式确定波束形成的协方差矩阵：其中：R
a
(f
j
)＝E{X
s
(f
j
)X
sH
(f
j
)}f
j
表示宽带信号带内第j个频点，J为宽带信号频点数，θ表示方向角。4.根据权利要求3所述的方法，其特征在于，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱包括：基于协方差矩阵以及权向量将波束形成转换为如下问题：基于协方差矩阵以及权向量将波束形成转换为如下问题：其中表示权向量，1
s
是(2M
‑
1)
×
1的1向量，表示寻找使标量最小时的方向角θ；求解问题得到方向角θ，则有权向量为：方位谱为：5.一种空域反转阵导向最小方差空间波束形成装置，其特征在于，包括：空域反转阵列构建模块，用于将原阵列快拍与自身进行空域反转过后的阵列快拍进行卷积，再将卷积后每个阵元的功率进行均衡，得到空域反转阵列的快拍；波束形成处理模块，用于获取空域反转阵的宽带阵列信号和方向向量，确定波束形成
的协方差矩阵，基于协方差矩阵以及权向量根据STMV波束形成方法确定方位谱。6.根据权利要求5所述的装置，其特征在于，所...

【专利技术属性】
技术研发人员：陈阳，梁景山，
申请(专利权)人：南京锐声海洋科技有限公司，
类型：发明
国别省市：