CN112133321B - 一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 - Google Patents
一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 Download PDFInfo
- Publication number
- CN112133321B CN112133321B CN202011007852.5A CN202011007852A CN112133321B CN 112133321 B CN112133321 B CN 112133321B CN 202011007852 A CN202011007852 A CN 202011007852A CN 112133321 B CN112133321 B CN 112133321B
- Authority
- CN
- China
- Prior art keywords
- signal
- gaussian
- gaussian noise
- signals
- follows
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000926 separation method Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000001629 suppression Effects 0.000 title claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 230000005236 sound signal Effects 0.000 claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 13
- 238000004891 communication Methods 0.000 claims abstract description 7
- 238000012880 independent component analysis Methods 0.000 claims description 19
- 238000013507 mapping Methods 0.000 claims description 14
- 241000282461 Canis lupus Species 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 230000000739 chaotic effect Effects 0.000 claims description 10
- 238000010606 normalization Methods 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 5
- 238000011946 reduction process Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 241000160777 Hipparchia semele Species 0.000 abstract 1
- 230000002087 whitening effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000005534 acoustic noise Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0264—Noise filtering characterised by the type of parameter measurement, e.g. correlation techniques, zero crossing techniques or predictive techniques
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0272—Voice signal separating
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0272—Voice signal separating
- G10L21/0308—Voice signal separating characterised by the type of parameter measurement, e.g. correlation techniques, zero crossing techniques or predictive techniques
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法,用于水声通信信号去噪。首先基于水声信号稀疏特性,采用稀疏编码理论预先去除高斯噪声,包括基于拉普拉斯分布模型构造稀疏惩罚函数,求解在极大似然估计情况下的稀疏分量。然后基于水声信号的非高斯特性,将非高斯随机噪声视作源信号,基于改进灰狼优化算法求解水声信号盲源分离问题。最后基于分离后信号和接收信号之间的相关系数矩阵提取目标信号,去除非高斯噪声。本发明能有效抑制复杂水声环境中高斯和非高斯噪声,提高水声信号接收质量。
Description
技术领域
本发明属于水声信号去噪领域,具体地说,涉及一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法。
背景技术
水下无线数据传输技术是建设海洋强国的关键技术。声波是目前唯一可以在水下远距离传播的信息载体,被广泛应用在水下通信领域。然而,声波在水下传输时,受到大量复杂背景噪声影响。水声噪声的来源可以划分为两种,第一种是人类活动引起的,如船舶航运、工业活动等机器声,第二种是自然因素引起的,如雨、海洋生物和地震活动等产生的噪声。
基于自适应噪声、互信息、排列熵和小波阈值的完全集成经验模态分解的水声信号去噪新方法,这种方式能较好的抑制大的噪声分量,但是计算复杂度较高;基于集成经验模态分解、独自分量分析(Independent Component Analysis,ICA)和小波阈值去噪的信号去噪方法,可降低噪声对信号的影响,但是盲源分离问题往往伴随对分离后源信号顺序不确定等问题;针对卷积混合模型,将时域带噪信号转化为频域进行分离,达到降噪目的,并解决了分离后信号顺序不确定问题;代数优化方法已经用于求解ICA目标函数,并基于余弦距离对分离后信号进行排序,具有较高的准确性;然而,实际水声环境中,有用信号和噪声信号通常为非高斯信号,具有较小的负熵,基于负熵的ICA方法分离效果较差。且水声通信系统对能量消耗敏感,高复杂度的ICA算法难以适用实际水下工程;根据水声信道的稀疏性,采用拉普拉斯分布代替负熵运算,实现了有效信号与非高斯噪声的区分,且降低了基于负熵算法的计算复杂度;基于代数迭代优化算法优化ICA目标函数,一定程度上降低了算法计算复杂度。
总体来看,当前水声信号去噪方法仍然存在去噪不彻底、算法复杂度高、得到的目标信号顺序不确定性等问题。
发明内容
本发明的目的是提出一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法,以弥补现有技术的不足。
为实现上述发明目的,本发明采用下述技术方案予以实现:
一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法,其特征在于,包括以下步骤:
S1:获取待去噪的含高斯/非高斯噪声的水声信号;
S2:抑制高斯噪声,得到不含高斯噪声的水声信号;
S3:对S2处理后的水声信号进行盲源分离,得到分离后源信号,即非高斯噪声和目标水声信号;
S4:基于相关系数矩阵确定S3中分离后源信号中的目标水声信号,从而去除非高斯噪声,最终得到高斯/非高斯去噪后的水声信号。
进一步的,所述S1具体如下:
S1-1:信号接收模型:
假设水声通信系统模型采用一发多收模式,发送目标信号为s0(t),将信道中的非高斯噪声信号视为n-1维非高斯信号[s1(t),s1(t),...,sn-1(t)]T,接收阵元个数为n,则接收端收到的信号为:
X(t)=h(t)[s0(t),s1(t),s2(t),...sn-1(t)]T+e(t)
其中,h(t)表示信道参数,e(t)表示高斯噪声,[g]T表示转置。
进一步的,所述S2中基于稀疏编码理论抑制高斯噪声。
更进一步的,所述S2具体如下:
S2-1:信号归一化处理:
对接收信号X(t)=[x1,x2,...,xn]T归一化处理,处理过程如下所示:
S2-2:求解拉普拉斯稀疏惩罚函数的收缩函数:
拉普拉斯的概率分布如下所示:
其中,s=E{X′2},E{X′2}表示求X′2均值;由于水声信号的稀疏性和非高斯性,采用极大似然估计,最大化信号的非高斯性,去除高斯噪声;对p(X′)取对数,得:
l(X′)=lnp(X′)
对X′求导得:
S2-3:求解极大似然估计情况下信号的稀疏分量,去除高斯噪声;
极大似然估计表达式如下所示:
M(X′)=sign(X′)max(0,|X′|-σ2|l'(X′)|)
其中sign(g)为符号函数,由下式给出:
将上式代入极大似然估计表达式,得到去除高斯噪声后信号:
式中σ2表示X′中高斯噪声方差,X″是去除高斯噪后的信号;其中X″为归一化信号;为便于S3处理,需对X″进行部分还原,使其仅保持0均值特性,还原过程如下:
X″′=σ1X″
其中σ1为X(t)的标准差,X″′为基于稀疏编码去噪后信号。
进一步的,所述S3中基于改进灰狼算法优化独立分量分析方法对信号进行盲源分离。
更进一步的,所述S3具体如下:
S3-1:对信号进行白化预处理:
白化可去除信号之间的相关性,同时能够简化后续独立分量提取过程并增强算法的收敛性,白化后的向量z=(z1,L,zm)T要求满足E{zzT}=I,I为单位矩阵;因此对S2输出信号X″′乘以下式所示白化矩阵:
W0=Λ-1/2UT
其中U和Λ分别代表X″′的协方差矩阵Cx的特征向量矩阵和特征值矩阵,白化后信号为:
Xwhite=W0X″′
S3-2:基于峭度最大化分离源信号:
峭度可描述信号非高斯特性,基于独立分量分析(ICA)算法,通过最大化分离后信号峭度,分离源信号,定义随机信号x的峭度如式下所示:
kurt(x)=E{x4}-3(E{x2})2
其中,E{g}表示求期望;传统基于牛顿迭代法的求解过程计算复杂度大,因此基于群智能优化理论,采用改进灰狼算法求解ICA问题,通过最大化信号峭度,获得最优分离矩阵W;适应度函数如下所示:
fit=|kurt(WXwhite)|
其中|g|表示绝对值;
S3-3:基于改进灰狼算法优化独立分量分析方法:
为提高灰狼算法优化性能基于混沌映射初始化种群,提高初始种群个体的多样性和计算效率,然后基于逆不完全Γ函数修改狩猎模型,提高算法的探索与开发能力,具体步骤如下:
第一步:初始化种群,利用混沌系统的周期性、随机性和规律性等特点,基于混沌映射初始化种群,即初始分离矩阵Wt=[Wt1,Wt2,...,Wtn],其中t表示迭代次数,t=0时,W0表示初始种群,W01,W02,...,W0n表示初始个体,n为种群大小;本发明采用Iterative混沌序列产生灰狼的初始位置,Iteraive映射的映射方程如下所示:
其中W01是随机生成的分离矩阵,W0(k+1)是经过Iteraive映射之后生成的初始分离矩阵,k=1,2,...,n,b是控制参数,b∈(0,1);
第二步:基于改进收敛因子的狩猎模型,当迭代次数为t时,计算种群内个体适应度值,根据适应度值保留前三个最大适应度值对应的个体位置,分别记作群体最优个体为Wtα,次优个体为Wtβ,第三最优个体为Wtδ;随机挑选种群中个体Wtk,k=1,2,...,n,并根据Wtα、Wtβ、Wtδ更新产生下一代个体,更新规则如下所示:
其中A1,A2,A3、C1,C2,C3分别是Wtk对应Wtα、Wtβ、Wtδ的系数向量,A1,A2,A3和C1,C2,C3分别由下式确定:
其中r1i,r2i在[0,1]范围内随机取值,a为收敛因子;本发明基于逆不完全Γ函数改进收敛因子a:
其中amax=2,最小值amin=0,tmax是最大迭代次数,λ是随机变量,λ≥0;
第三步:收敛条件,当算法达到最大迭代次数tmax,或者满足下式时,根据当代最优个体Wtα输出最优分离矩阵Wopt=Wtα;
||Wtα-W(t-1)α||2<σ
其中||g||2表示求二范数,σ一般取值为σ=10-7;分离矩阵Wopt确定之后,分离后信号表示为下所示:
y'(t)=WoptXwhite。
进一步的,所述S4具体如下:
S4-1:计算分离信号和观测信号的相关系数矩阵:
定义第i个分离后的独立分量y'i(t)对第j个观测数据分量xj(t)的相关系数矩阵为:
其中Cij表示向量y'i(t)和xj(t)的协方差矩阵,Cii和Cjj分别为y'i(t)和xj(t)的方差矩阵;计算相关系数矩阵P得:
其主对角线元素相等且值为1,副对角线元素pi,j相等且取值范围在[-1,1]内,当且仅当y'i(t),xj(t)互相独立时为0;
S4-2:对分离后信号进行排序和相位调整:
首先分别选出矩阵P中各行元素中除主对角线之外的绝对值最大的元素Pi,j,且每列也只能取一个元素,若有同一列出现了两个行的最大绝对值,那么选取较大的绝对值元素所在的列,并排除这一列,在剩下的列中再运用相同的方法进行排序,直到将各行相应的排序依据都选出来;
依据选出的排序标记进行排序,具体准则为各行选定的绝对值最大者所在的列即为重新排序后信号所处的位置。
本发明保证了目标源信号的可靠性,最终实现短时间、高准确率的水声信号分类去噪,得到顺序正确的源信号。
本发明的优点和技术效果如下:
本发明根据水声信号的稀疏性,首先基于拉普拉斯分布的稀疏模型,利用其收缩函数,可有效抑制高斯噪声;其次基于改进灰狼优化算法提高ICA方法分离性能,同时基于分离后的信号之间的相关系数对分离后信号排序,提取水声目标信号,可有效去除非高斯噪声。本发明能够降低高斯/非高斯复杂海洋噪声对水声信号质量的影响,大幅度提升水声接收信号质量。
附图说明
图1是本发明的一种实施例的整体模型图;
图2是本发明的一种实施例的技术路线图;
图3是本发明的一种实施例中的整体流程图;
图4是本发明的一种实施例中改进灰狼优化算法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下将结合附图和实施例,对本发明作进一步详细说明。
实施例1:
声波在水下传输时,发送信号受噪声影响,水声通信质量降低,噪声的复杂性和随机性增加了对接收信号的去噪难度,使用单一的信号去噪方式不易完全去除信号的多种复杂噪声,根据噪声特点对噪声进行分类去除,能够保障噪声去除的有效性,得到准确的目标信号,具体模型如图1所示。
基于优化算法的信号盲源分离技术具有较好的信号分离性能,其主要步骤为:首先使用优化算法得到一个分离矩阵,用分离矩阵乘以接收信号得到分离信号。但水下噪声复杂,当前水声信号去噪方法存在去噪不彻底、算法复杂度高、得到的目标信号顺序不确定等问题。如何找到快速收敛的优化算法、并降低盲源分离算法计算复杂度以及对分离后信号的排序,是本实施例所要解决的技术问题。
本实施例提出了一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法,其技术路线图如图2所示,包括以下步骤:
S1:获取待去噪的含高斯/非高斯噪声的水声信号,具体步骤如下:
S1-1:信号接收模型:
假设水声通信系统模型采用一发多收模式,发送目标信号为s0(t),将信道中的非高斯噪声信号视为n-1维非高斯信号[s1(t),s1(t),...,sn-1(t)]T,接收阵元个数为n;则接收端收到的信号为:
X(t)=h(t)[s0(t),s1(t),s2(t),...sn-1(t)]T+e(t)
其中,h(t)表示信道参数,e(t)表示高斯噪声,[g]T表示转置。
S2:基于稀疏编码理论抑制高斯噪声,得到不含高斯噪声的水声信号,具体步骤如下:
S2-1:信号归一化处理:
对接收信号X(t)=[x1,x2,...,xn]T归一化处理,处理过程如下所示:
S2-2:求解拉普拉斯稀疏惩罚函数的收缩函数:
拉普拉斯的概率分布如下所示:
其中,s=E{X′2},E{X′2}表示求X′2均值;由于水声信号的稀疏性和非高斯性,采用极大似然估计,最大化信号的非高斯性,去除高斯噪声;对p(X′)取对数,得:
l(X′)=lnp(X′)
对X′求导得:
S2-3:求解极大似然估计情况下信号的稀疏分量,去除高斯噪声;
极大似然估计表达式如下所示:
M(X′)=sign(X′)max(0,|X′|-σ2|l'(X′)|)
其中sign(g)为符号函数,由下式给出:
将上式代入极大似然估计表达式,得到去除高斯噪声后信号:
式中σ2表示X′中高斯噪声方差,X″是去除高斯噪后的信号;其中X″为归一化信号;为便于S3处理,需对X″进行部分还原,使其仅保持0均值特性,还原过程如下:
X″′=σ1X″
其中σ1为X(t)的标准差,X″′为基于稀疏编码去噪后信号。
S3:基于改进灰狼算法优化独立分量分析方法,并分离S2中不含高斯噪声的水声信号,得到分离后源信号,具体步骤如下:
S3-1:对信号进行白化预处理:
白化可去除信号之间的相关性,同时能够简化后续独立分量提取过程并增强算法的收敛性,白化后的向量z=(z1,L,zm)T要求满足E{zzT}=I,I为单位矩阵;因此对S2得到信号X″′乘以下式所示白化矩阵:
W0=Λ-1/2UT
其中U和Λ分别代表X″′的协方差矩阵Cx的特征向量矩阵和特征值矩阵,白化后信号为:
Xwhite=W0X″′
S3-2:基于峭度最大化分离源信号:
峭度可描述信号非高斯特性,本实施例基于ICA算法,通过最大化分离后信号峭度,分离源信号,定义随机信号x的峭度如式下所示:
kurt(x)=E{x4}-3(E{x2})2
其中,E{g}表示求期望;传统基于牛顿迭代法的求解过程计算复杂度大,因此本实施例基于群智能优化理论,采用改进灰狼算法求解ICA问题,通过最大化信号峭度,获得最优分离矩阵W;适应度函数如下所示:
fit=|kurt(WXwhite)|
其中|g|表示绝对值;
S3-3:基于改进灰狼算法优化独立分量分析方法:
为提高灰狼算法优化性能,基于混沌映射初始化种群,提高初始种群个体的多样性和计算效率,然后基于逆不完全Γ函数修改狩猎模型,提高算法的探索与开发能力,具体步骤如下:
第一步:初始化种群,利用混沌系统的周期性、随机性和规律性等特点,基于混沌映射初始化种群,即初始分离矩阵Wt=[Wt1,Wt2,...,Wtn],其中t表示迭代次数,t=0时,W0表示初始种群,W01,W02,...,W0n表示初始个体,n为种群大小;采用Iterative混沌序列产生灰狼的初始位置,Iteraive映射的映射方程如下所示:
其中W01是随机生成的分离矩阵,W0(k+1)是经过Iteraive映射之后生成的初始分离矩阵,k=1,2,...,n,b是控制参数,b∈(0,1);
第二步:基于改进收敛因子的狩猎模型,当迭代次数为t时,计算种群内个体适应度值,根据适应度值保留前三个最大适应度值对应的个体位置,分别记作群体最优个体为Wtα,次优个体为Wtβ,第三最优个体为Wtδ;随机挑选种群中个体Wtk,k=1,2,...,n,并根据Wtα、Wtβ、Wtδ更新产生下一代个体,更新规则如下所示:
其中A1,A2,A3、C1,C2,C3分别是Wtk对应Wtα、Wtβ、Wtδ的系数向量,A1,A2,A3和C1,C2,C3分别由下式确定:
其中r1i,r2i在[0,1]范围内随机取值,a为收敛因子;基于逆不完全Γ函数改进收敛因子a:
其中amax=2,最小值amin=0,tmax是最大迭代次数,λ是随机变量,λ≥0;
第三步:收敛条件,当算法达到最大迭代次数tmax,或者满足下式时,根据当代最优个体Wtα输出最优分离矩阵Wopt=Wtα;
||Wtα-W(t-1)α||2<σ
其中||g||2表示求二范数,σ一般取值为σ=10-7,算法流程图如图4所示,伪代码如下;
分离矩阵Wopt确定之后,分离后信号表示为下所示:
y'(t)=WoptXwhite
S4:基于相关系数确定S3中分离后源信号的顺序,提取目标水声信号,去除非高斯噪声,最终得到去噪后信号,具体步骤如下:
S4-1:计算分离信号和观测信号的相关系数矩阵:
定义第i个分离后的独立分量y'i(t)对第j个观测数据分量xj(t)的相关系数矩阵为:
其中Cij表示向量y'i(t)和xj(t)的协方差矩阵,Cii和Cjj分别为y'i(t)和xj(t)的方差矩阵;计算相关系数矩阵P得:
其主对角线元素相等且值为1,副对角线元素pi,j相等且取值范围在[-1,1]内,当且仅当y'i(t),xj(t)互相独立时为0;
S4-2:对分离后信号进行排序和相位调整:
首先分别选出矩阵P中各行元素中除主对角线之外的绝对值最大的元素Pi,j,且每列也只能取一个元素,若有同一列出现了两个行的最大绝对值,那么选取较大的绝对值元素所在的列,并排除这一列,在剩下的列中再运用相同的方法进行排序,直到将各行相应的排序依据都选出来;
依据选出的排序标记进行排序,具体准则为各行选定的绝对值最大者所在的列即为重新排序后信号所处的位置。
上述方法保证了目标源信号的可靠性,最终实现短时间、高准确率的水声信号分类去噪,得到顺序正确的源信号。
例如,假设源信号为[s1,s2,s3]T,分离后信号为[y′1,y′2,y′3]T。计算分离后信号相关系数矩阵如下所示:
则除对角线元素绝对值最大的元素如下所示:
其中,第一行除对角线元素绝对值最大元素为0.98。第二行按照规则应当取0.98,而由于具体实施步骤第一步中“每列也只能取一个元素”原则,将导致第三行只剩对角线元素可取,与具体实施步骤矛盾。因此第二行取第三列数字0.87,第三行取-0.55。调整分离后分离信号如下所示:
y(t)=[y'3,y'1,y'2]
其中y'3,y'1,y'2分别对应源信号s1,s2,s3;基于此,本实施例对盲源分离后信号y'(t)进行顺序调整后,即可得到对源信号[s0(t),s1(t),s2(t),...sn-1(t)]T的准确估计[s′0(t),s1(t),s′2(t),...s′n-1(t)]T,从而选取第一路信号s′0(t),去除其他信号,作为对目标信号s0(t)的准确估计。
本实施例针对复杂海洋噪声,根据水声信号的稀疏性,基于拉普拉斯稀疏模型,求解收缩函数抑制高斯噪声,并基于ICA方法分离接收到的混合信号,同时基于相关系数确定分离后信号顺序,提取期望信号,去除非高斯噪声,最终实现基于盲源分离的水声信号分类去噪方法。
以上所述之实施例子只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (5)
1.一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法,其特征在于,包括以下步骤:
S1:获取待去噪的含高斯/非高斯噪声的水声信号;
S2:抑制高斯噪声,得到不含高斯噪声的水声信号;
S3:对S2处理后的水声信号进行盲源分离,得到分离后源信号,即非高斯噪声和目标水声信号;所述S3中基于改进灰狼算法优化独立分量分析方法对信号进行盲源分离;所述S3中改进灰狼算法优化独立分量分析方法具体为:
第一步:初始化种群,基于混沌映射初始化种群,即初始分离矩阵Wt=[Wt1,Wt2,...,Wtn],其中t表示迭代次数,t=0时,W0表示初始种群,W01,W02,...,W0n表示初始个体,n为种群大小;采用Iterative混沌序列产生灰狼的初始位置,Iteraive映射的映射方程如下所示:
其中W01是随机生成的分离矩阵,W0(k+1)是经过Iteraive映射之后生成的初始分离矩阵,k=1,2,...,n,b是控制参数,b∈(0,1);
第二步:基于改进收敛因子的狩猎模型,当迭代次数为t时,计算种群内个体适应度值,根据适应度值保留前三个最大适应度值对应的个体位置,分别记作群体最优个体为Wtα,次优个体为Wtβ,第三最优个体为Wtδ;随机挑选种群中个体Wtk,k=1,2,...,n,并根据Wtα、Wtβ、Wtδ更新产生下一代个体,更新规则如下所示:
其中A1,A2,A3、C1,C2,C3分别是Wtk对应Wtα、Wtβ、Wtδ的系数向量,A1,A2,A3和C1,C2,C3分别由下式确定:
其中r1i,r2i在[0,1]范围内随机取值,a为收敛因子;基于逆不完全Г函数改进收敛因子a:
其中amax=2,最小值amin=0,tmax是最大迭代次数,λ是随机变量,λ≥0;
第三步:收敛条件,当算法达到最大迭代次数tmax,或者满足下式时,根据当代最优个体Wtα输出最优分离矩阵Wopt=Wtα;
||Wtα-W(t-1)α||2<σ
其中||g||2表示求二范数,σ取值为σ=10-7;分离矩阵Wopt确定之后,分离后信号表示为下所示:
y′(t)=WoptXwhite;
S4:基于相关系数矩阵确定S3中分离后源信号中的目标水声信号,从而去除非高斯噪声,最终得到高斯/非高斯去噪后的水声信号。
2.如权利要求1所述的水声信号高斯/非高斯噪声抑制方法,其特征在于,所述S1具体如下:建立信号接收模型:
假设水声通信系统模型采用一发多收模式,发送目标信号为s0(t),将信道中的非高斯噪声信号视为n-1维非高斯信号[s1(t),s2(t),...,sn-1(t)]T,接收阵元个数为n,则接收端收到的信号为:
X(t)=h(t)[s0(t),s1(t),s2(t),...sn-1(t)]T+e(t)
其中,h(t)表示信道参数,e(t)表示高斯噪声,[g]T表示转置。
3.如权利要求1所述的水声信号高斯/非高斯噪声抑制方法,其特征在于,所述S2中基于稀疏编码理论抑制高斯噪声。
4.如权利要求3所述的水声信号高斯/非高斯噪声抑制方法,其特征在于,所述S2具体如下:
S2-1:信号归一化处理:
对接收信号X(t)=[x1,x2,...,xn]T归一化处理,处理过程如下所示:
S2-2:求解拉普拉斯稀疏惩罚函数的收缩函数:
拉普拉斯的概率分布如下所示:
其中,s=E{X′2},E{X′2}表示求X′2均值;由于水声信号的稀疏性和非高斯性,采用极大似然估计,最大化信号的非高斯性,去除高斯噪声;对p(X′)取对数,得:
l(X′)=lnp(X′)
对X′求导得:
S2-3:求解极大似然估计情况下信号的稀疏分量,去除高斯噪声;
极大似然估计表达式如下所示:
M(X′)=sign(X′)max(0,|X′|-σ2|l'(X′)|)
其中sign(g)为符号函数,由下式给出:
将上式代入极大似然估计表达式,得到去除高斯噪声后信号:
式中σ2表示X′中高斯噪声方差,X″是去除高斯噪后的信号;其中X″为归一化信号;为便于S3处理,需对X″进行部分还原,使其仅保持0均值特性,还原过程如下:
X″′=σ1X″
其中σ1为X(t)的标准差,X″′为基于稀疏编码去噪后信号。
5.如权利要求1所述的水声信号高斯/非高斯噪声抑制方法,其特征在于,所述S4具体如下:
S4-1:计算分离信号和观测信号的相关系数矩阵:
定义第i个分离后的独立分量y′i(t)对第j个观测数据分量xj(t)的相关系数矩阵为:
其中Cij表示y′i(t)和xj(t)的协方差矩阵,Cii和Cjj分别为y′i(t)和xj(t)的方差矩阵;计算相关系数矩阵P得:
其主对角线元素相等且值为1,副对角线元素pi,j相等且取值范围在[-1,1]内,当且仅当y′i(t),xj(t)互相独立时为0;
S4-2:对分离后信号进行排序和相位调整:
首先分别选出矩阵P中各行元素中除主对角线之外的绝对值最大的元素Pi,j,且每列也只取能一个元素,若有同一列出现了两个行的最大绝对值,那么选取较大的绝对值元素所在的列,并排除这一列,在剩下的列中再运用相同的方法进行排序,直到将各行相应的排序依据都选出来;依据选出的排序标记进行排序,具体准则为各行选定的绝对值最大者所在的列即为重新排序后信号所处的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011007852.5A CN112133321B (zh) | 2020-09-23 | 2020-09-23 | 一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011007852.5A CN112133321B (zh) | 2020-09-23 | 2020-09-23 | 一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112133321A CN112133321A (zh) | 2020-12-25 |
CN112133321B true CN112133321B (zh) | 2021-05-14 |
Family
ID=73842790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011007852.5A Active CN112133321B (zh) | 2020-09-23 | 2020-09-23 | 一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112133321B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112953607B (zh) * | 2021-02-22 | 2022-08-09 | 西安交通大学 | 一种mimo-ofdm系统量化噪声消除方法、介质及设备 |
CN113241087B (zh) * | 2021-05-10 | 2022-05-27 | 青岛科技大学 | 基于改进vmd和改进小波阈值法相结合的水声信号去噪方法 |
CN113779805B (zh) * | 2021-09-16 | 2023-11-14 | 北京中安智能信息科技有限公司 | 海洋噪声相关性仿真方法和装置、设备及存储介质 |
CN114333897B (zh) * | 2022-03-14 | 2022-05-31 | 青岛科技大学 | 基于多信道噪声方差估计的BrBCA盲源分离方法 |
CN115932732B (zh) * | 2022-12-02 | 2023-09-22 | 中国人民解放军海军航空大学 | 一种基于改进NMF和FastICA的联合水声多目标盲分离方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004109377A (ja) * | 2002-09-17 | 2004-04-08 | Nippon Hoso Kyokai <Nhk> | ブラインド信号分離方法、装置およびそのためのコンピュータプログラム |
JP4769238B2 (ja) * | 2007-08-24 | 2011-09-07 | 日本電信電話株式会社 | 信号分離装置、信号分離方法、プログラム及び記録媒体 |
CN102546494B (zh) * | 2010-12-30 | 2015-05-06 | 中国科学院声学研究所 | 一种水下盲均衡系统及方法 |
CN103679199B (zh) * | 2013-12-11 | 2017-03-01 | 中国石油大学(华东) | 基于入侵性杂草算法的有噪独立分量分析方法 |
CN107884751B (zh) * | 2017-09-29 | 2021-03-12 | 广东工业大学 | 一种利用单通道接收信号作信源数目估计的方法 |
CN108304855B (zh) * | 2017-12-05 | 2020-06-23 | 厦门大学 | 一种海洋环境下多潜水艇特征信号盲源分离方法 |
CN109243483B (zh) * | 2018-10-17 | 2022-03-08 | 西安交通大学 | 一种含噪频域卷积盲源分离方法 |
CN110580911B (zh) * | 2019-09-02 | 2020-04-21 | 青岛科技大学 | 一种可抑制多个非稳定次高斯干扰的波束形成方法 |
CN111681667B (zh) * | 2020-06-23 | 2021-05-04 | 青岛科技大学 | 基于自适应窗口滤波和小波阈值优化的水声信号去噪方法 |
-
2020
- 2020-09-23 CN CN202011007852.5A patent/CN112133321B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112133321A (zh) | 2020-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112133321B (zh) | 一种基于盲源分离的水声信号高斯/非高斯噪声抑制方法 | |
CN113241087B (zh) | 基于改进vmd和改进小波阈值法相结合的水声信号去噪方法 | |
CN109841226B (zh) | 一种基于卷积递归神经网络的单通道实时降噪方法 | |
CN111653288B (zh) | 基于条件变分自编码器的目标人语音增强方法 | |
Tan et al. | Nonlinear blind source separation using higher order statistics and a genetic algorithm | |
CN111971743A (zh) | 用于改进的实时音频处理的系统、方法和计算机可读介质 | |
CN112735460B (zh) | 基于时频掩蔽值估计的波束成形方法及系统 | |
CN109559329B (zh) | 一种基于深度去噪自动编码器的粒子滤波跟踪方法 | |
CN110223708B (zh) | 基于语音处理的语音增强方法及相关设备 | |
Sun et al. | Cluster guide particle swarm optimization (CGPSO) for underdetermined blind source separation with advanced conditions | |
CN111368710A (zh) | 一种联合深度学习的地震数据随机噪声压制方法 | |
CN110490816B (zh) | 一种水下异构信息数据降噪方法 | |
CN111816200B (zh) | 一种基于时频域二值掩膜的多通道语音增强方法 | |
CN107292855B (zh) | 一种结合自适应非局部样本和低秩的图像去噪方法 | |
Wakayama et al. | CNN-transformer with self-attention network for sound event detection | |
CN115146667A (zh) | 一种基于曲波变换与多分枝深层自编码的多尺度地震噪声压制方法 | |
CN114897002A (zh) | 基于LPINet的低截获概率雷达信号脉内调制识别方法 | |
CN108710836B (zh) | 一种基于级联特征提取的唇部检测及读取方法 | |
CN114999525A (zh) | 一种基于神经网络的轻量环境声音识别方法 | |
CN112037813B (zh) | 一种针对大功率目标信号的语音提取方法 | |
Chen et al. | Underwater target recognition method based on convolution autoencoder | |
Ansari et al. | Blind Source Separation Based on Genetic Algorithm-Optimized Multiuser Kurtosis | |
Rajput et al. | Noisy sonar image segmentation using reptile search algorithm | |
CN116756491A (zh) | 一种基于蜣螂优化算法优化小波阈值的阀门信号降噪方法 | |
CN113780301B (zh) | 防御对抗攻击的自适应去噪机器学习应用方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |