一种基于机器学习的水声信号处理算法制造技术

本发明专利技术涉及一种基于机器学习的水声信号处理算法，该算法采用EM准则估计信号和噪声的功率，使水声信号处理算法可得到超分辨的谱估计。与传统方法相比，该算法无需已知信号的数量，可用于低信噪比、快拍数较少情况下的相干信号的处理。信号的处理。信号的处理。

【技术实现步骤摘要】
一种基于机器学习的水声信号处理算法


[0001]本专利技术涉及雷达、声纳
，尤其涉及一种基于机器学习的水声信号处理算法。

技术介绍

[0002]信号的波达方向(DOA)估计是水声信号处理中一个十分重要的方向，其主要目的是估计信号的空间位置，这也是雷达、声纳等许多领域的重要任务之一。其实现方法为得到信号在空间各个方向上的能量分布，进而估计出信号的波达方向。常用的方法为常规波束形成法、MVDR、MUSIC等。
[0003]现有的DOA估计方法存在的问题主要有以下几个：1)常规波束形成法的分辨力受阵列孔径的限制，不具备很高的分辨能力，无法分辨两个相距较近的目标；2)一些高分辨算法在利用信号协方差矩阵时需要对其进行矩阵分解等操作，必须要求其为满秩，而包含相干信号的接收数据协方差矩阵不再为满秩矩阵，导致算法性能下降甚至失效；3)需要较大的快拍数，而实际工程应用中，快拍数一般十分有限，这将导致算法性能明显下降。

技术实现思路

[0004]本专利技术的目的在于提供一种基于机器学习的水声信号处理算法，以解决上述
技术介绍
中遇到的问题。
[0005]为实现上述目的，本专利技术的技术方案如下：
[0006]一种基于机器学习的水声信号处理算法，包括以下步骤：
[0007]步骤一、设定本方案的约束条件；阵列为均匀线列阵，阵元数为N，阵元间距为d；信源信号为远场窄带信号，阵元数要大于信号个数；传播介质是均匀且各向同性的；观测噪声为加性噪声，符合均值为0，方差为σ2的高斯分布，噪声与入射信号相互独立。
[0008]步骤二、设定每个波束上的信号方差为γ＝[γ1,γ2,...γ
M
]T
∈R
M
，令γ
m
＝1，m＝1,2,...M，设定加性噪声方差初始值σ2＝0.1，设定收敛速度阈值ε
min
＝0.001以及迭代次数最大值j
max
＝500；
[0009]步骤三、计算观测矩阵Y的协方差矩阵Σ
y
，
[0010]Σ
y
＝σ2I
N
+AΓA
H
[0011]其中，I
N
为N阶单位矩阵，Γ＝diag(γ)，A为阵列导向矢量矩阵，且A＝[a1,...,a
M
]∈C
N
×
M
，c为水中声速，ω为角频率，M为波束数，
[0012]步骤四、采用EM准则更新γ的第m个元素和噪声方差σ2；
[0013][0014][0015]其中，L为快拍数，Y为阵列观测数据矩阵，μ
X
为信源X的均值且
[0016]步骤五、更新收敛速度ε＝||γ
new
‑
γ
old
||1/||γ
old
||1；
[0017]步骤六、迭代步骤三至步骤五，直到收敛速度大于阈值或者达到最大迭代次数；
[0018]步骤七、输出γ
new
和(σ2)
new
，根据γ
new
确定来波方向。
[0019]与现有技术相比，本专利技术的有益效果是：采用EM准则估计信号和噪声的功率，使水声信号处理算法可得到超分辨的谱估计。与传统方法相比，该算法无需已知信号的数量，可用于低信噪比、快拍数较少情况下的相干信号的处理。
附图说明
[0020]参照附图来说明本专利技术的公开内容。应当了解，附图仅仅用于说明目的，而并非意在对本专利技术的保护范围构成限制。在附图中，相同的附图标记用于指代相同的部件。其中：
[0021]图1为本专利技术基于机器学习的水声信号处理算法的流程图。
具体实施方式
[0022]为了使本专利技术实现的技术手段、创作特征、达成目的与功效易于明白了解，现在结合附图对本专利技术作进一步详细的说明。这些附图均为简化的示意图，仅以示意方式说明本专利技术的基本结构，因此其仅显示本专利技术有关的构成。
[0023]根据本专利技术的技术方案，在不变更本专利技术实质精神下，本领域的一般技术人员可以提出可相互替换的多种结构方式以及实现方式。因此，以下具体实施方式以及附图仅是对本专利技术的技术方案的示例性说明，而不应当视为本专利技术的全部或者视为对本专利技术技术方案的限定或限制。
[0024]下面结合附图和实施例对本专利技术的技术方案做进一步的详细说明。
[0025]如图1所示，一种基于机器学习的水声信号处理算法，包括以下步骤：
[0026]步骤一、设定本方案的约束条件：阵列为均匀线列阵，阵元数为N，阵元间距为d；信源信号为远场窄带信号，阵元数要大于信号个数；传播介质是均匀且各向同性的；观测噪声为加性噪声，符合均值为0，方差为σ2的高斯分布，噪声与入射信号相互独立。
[0027]设定波束数为M，快拍数为L，信源信号矩阵可表示为X＝[x1,...,x
L
]∈C
M
×
L
。线性回归模型，阵元数为N的阵列观测矩阵Y可表示为:Y＝AX+N
[0028]其中，A为导向矢量矩阵A＝[a1,a2,...,a
M
]∈C
N
×
M
，，N为加性高斯噪声N∈C
N
×
L
。
[0029]步骤二、设定每个波束上的信号方差功率为γ＝[γ1,γ2,...γ
M
]T
∈R
M
，令γ
m
＝1，m＝1,2,...M，设定加性噪声方差初始值σ2＝0.1，设定收敛速度阈值ε
min
＝0.001以及迭代次数最大值j
max
＝500；
[0030]步骤三、计算观测矩阵Y的协方差矩阵Σ
y
，
[0031]Σ
y
＝σ2I
N
+AΓA
H
[0032]其中，I
N
为N阶单位矩阵，Γ＝diag(γ)，c为水中声速，ω为角频率，M为波束数，
[0033]观测矩阵Y服从高斯分布：
[0034][0035]根据先验信息，信源信号符合零均值的高斯分布，每个波束上的信号的方差(能量)为γ
m
∈γ＝[γ1,...γ
M
]T
。信源信号的均值和方差为：
[0036][0037][0038]其中Γ＝diag(γ)，Σ
Y
为观测矩阵Y的协方差矩阵：Σ
y
＝σ2I
N
+AΓA
H
，I
N
为N阶单位矩阵。
[0039]在本步骤中，本算法中并未对协方差矩阵做矩阵分解等操作，所以协方差矩阵是否满秩并无影响。
[0040]步骤四、采用EM准则更新γ的第m个元素和噪声方差σ2；
[0041][0042][0本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于机器学习的水声信号处理算法，其特征在于，包括以下步骤：步骤一、设定本方案的约束条件；步骤二、设定每个波束上的信号方差为γ＝[γ1,γ2,...γ
M
]
T
∈R
M
，令γ
m
＝1，m＝1,2,...M，设定加性噪声方差初始值σ2＝0.1，设定收敛速度阈值ε
min
＝0.001以及迭代次数最大值j
max
＝500；步骤三、计算观测矩阵Y的协方差矩阵Σ
y
，Σ
y
＝σ2I
N
+AΓA
H
其中，I
N
为N阶单位矩阵，Γ＝diag(γ)，A为阵列导向矢量矩阵，且A＝[a1,...,a
M
]∈C
N
×
M
，c为水中声速，ω为角频率，M为波束数，步骤四、采用EM准则...

【专利技术属性】
技术研发人员：陆桦，
申请(专利权)人：北京中安智能信息科技有限公司，
类型：发明
国别省市：