一种海洋环境下多潜水艇特征信号盲源分离方法

摘要

本发明公开了一种海洋环境下多潜水艇特征信号盲源分离方法，包括以下步骤：采集信号样本；在时频域中基于Hermitian角的样本向量和参考向量，采用聚类算法对混合信号样本进行聚类以获得多个聚类簇，将聚类算法中的隶属函数作为掩码；基于估计簇数目进行聚类验证；基于K‑均值聚类算法将掩码聚类成Q个掩码聚类簇；通过掩码估计，进行信号分离。本发明能够解决欠定情况下盲源分离问题，在掩码估计中不需要估计任何混合信号矩阵或源位置，简化了信号分离的复杂度，提高了运算效率。

说明书

技术领域

本发明涉及水声信号处理领域，特别涉及一种海洋环境下多潜水艇特征信号盲源分离方法及系统。

背景技术

当前的盲分离算法存在诸多因素限制，首先盲分离算法主要用在信号的瞬态混合形式，信号卷积混合形式的盲分离算法相对较少，大致可分为以下几类：1)基于消除信号分量间互相关性的信号分离网，该网络利用反馈电路来解决问题。该算法在实际应用中有较大缺陷，在信号之间的尺度相差很大或混合矩阵是病态的情况下，分离效果比较差，甚至无法分离信号。2)利用非线性传递函数对输出进行变换，如基于信息传输最大化原理的嫡最大化的方法。此算法收敛速度很慢，同时分离矩阵的求逆又带来数值不稳定。3)非线性主分量分析(PCA)算法，它是线性主分量分析方法的推广。另外，以上算法在在掩码估计中需要估计任何混合信号矩阵或源位置，使得信号分离过程相对复杂，运算效率不高。

发明内容

本发明的目的在于提供一种海洋环境下多潜水艇特征信号盲源分离方法，其能够解决欠定情况下盲源分离问题，简化了信号分离的复杂度，提高了运算效率。

为实现上述目的，本发明采用以下技术方案：

一种海洋环境下多潜水艇特征信号盲源分离方法，包括以下步骤：

S1、采集混合信号样本；

S2、在时频域中基于Hermitian角的样本向量和参考向量，采用聚类算法对混合信号样本进行聚类以获得多个聚类簇，将聚类算法中的隶属函数作为掩码；

S3、基于估计簇数目进行聚类验证；

S4、基于K-均值聚类算法将掩码聚类成Q个掩码聚类簇，记为c_q,q＝1,...Q，并满足D_q,q＝1,...,Q的总和是最小的，D_q是在第q个簇里掩码与簇质心的总距离，即：

其中，是第k频点的第i个掩码，C_q是第q簇c_q的质心，是和簇质心C_q之间的皮尔森相关系数，k_st和k_end是用于聚类的相邻频点组的开始和结束频点，即频点总的数目是k_end-k_st+1，作为距离度量，这样高度相关的掩码(较小的距离)将来自一个簇；

S5、通过掩码估计，进行信号分离。

优选地，步骤S2中的所述聚类算法采用K-均值聚类或模糊C-均值聚类算法。

采用上述技术方案后，本发明与背景技术相比，具有如下优点：

本发明基于Hermitian角的样本向量和参考向量，使用聚类算法聚集了混合信号样本，对未知信号进行分离，能够解决欠定情况下盲源分离问题，在掩码估计中不需要估计任何混合信号矩阵或源位置，简化了信号分离的复杂度，提高了运算效率。

附图说明

图1为本发明的流程示意图；

图2为源信号与接收器的排列方式示意图；

图3为信道模型示意图；

图4a给出了蓝鲸的原始声信号波形，图4b给出了座头鲸的原始声信号波形；

图5a给出了蓝鲸的原始声音频谱图，图5b给出了座头鲸的原始声音频谱图；

图6a给出了接收端R1捕获的混合声音波形，图6b给出了接收端R2捕获的混合声音波形；

图7a为分离出的蓝鲸声音波形，图7b为分离出的座头鲸声音波形；

图8a为分离出的蓝鲸声音频谱图，图8b为分离出的座头鲸声音频谱图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

实施例一

请参阅图1，本发明公开了一种海洋环境下多潜水艇特征信号盲源分离方法，包括以下步骤：

S1、采集混合信号样本。

S2、在时频域中基于Hermitian(厄米矩阵)角的样本向量和参考向量，采用聚类算法对混合信号样本进行聚类以获得多个聚类簇，将聚类算法中的隶属函数作为掩码。聚类算法采用K-均值聚类或模糊C-均值聚类算法。

S3、基于估计簇数目进行聚类验证。

S4、基于K-均值聚类算法将掩码聚类成Q个掩码聚类簇，记为c_q,q＝1,...Q，并满足D_q,q＝1,...,Q的总和是最小的，D_q是在第q个簇里掩码与簇质心的总距离，即：

其中，是第k频点的第i个掩码，C_q是第q簇c_q的质心，是和簇质心C_q之间的皮尔森相关系数，k_st和k_end是用于聚类的相邻频点组的开始和结束频点，即频点总的数目是k_end-k_st+1，作为距离度量，这样高度相关的掩码(较小的距离)将来自一个簇。

S5、通过掩码估计，进行信号分离。

实验评价

一、信道设计

在水下模拟两个信号信号接收端(R1和R2)，两个声信号发射端(S1和S2)，位置任意。其中R1与S1相距3.5km，R2与S2相距4.0km，R1与S2相距4.5km，R2与S1相距4.7km，排列方式如图2所示。

R1与S1形成信道h11、R1与S2形成信道h12、R2与S1形成信道h21、R2与S2形成信道h22，四个信道模型如图3所示。

二、仿真实验

该仿真实验是在Matlab2014b平台上进行的，电脑的处理器为AMD双核A6-4400M2.7GHz，内存为4G。

仿真实验使用的声音分别是蓝鲸、座头鲸的叫声，分别进行了不同种类鲸鱼声音分离实验和两头蓝鲸声音分离实验。

图4a给出了蓝鲸的原始声信号波形，图4b给出了座头鲸的原始声信号波形；图5a给出了蓝鲸的原始声音频谱图，图5b给出了座头鲸的原始声音频谱图；图6a给出了接收端R1捕获的混合声音波形，图6b给出了接收端R2捕获的混合声音波形。

三、分离结果分析

按本发明的盲源分离方法对混合声音信号进行分离，分离出的信号波形如图7所示，其中图7a为分离出的蓝鲸声音波形，图7b为分离出的座头鲸声音波形。分别将图7a与图4a、图7b与图4b进行比较，可以看出本发明的盲源分离方法能够较好的将混合信号进行分离。

为了更进一步检验，对分离出的信号做短时傅里叶变换(STFT)，得到的频谱图像如图8所示，其中图8a为分离出的蓝鲸声音频谱图，图8b为分离出的座头鲸声音频谱图。将图8a与图5a进行比较，二者底部的特征能量是一致的；将图8b与图5b进行比较，二者底部的特征能量是一致的。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。