一种基于测量域噪声相减的压缩感知语音增强方法与流程
1.本发明涉及语音增强领域，特别涉及一种基于测量域噪声相减的压缩感知语音增强方法。
背景技术：
2.自人类文明开始以来，语音一直是日常生活和工作交流的主要形式。但是在实际情境中，来自餐厅、市场和工厂等环境的附加背景噪声会对语音信号造成污染。当噪声强度较大时，还会掩盖真实的语音信号，使信息的交流变得困难。当人与人在噪声严重的环境下进行交流时，为了对抗其他声音的干扰，不得不拉近交流的距离或者提高交流的声音，这使得交流者发音器官和听觉系统的负担加重。近年来，由于人工智能技术的兴起，智能系统和设备的出现方便了人们的生活。语音是人与智能设备进行人机交互时最便捷、最准确的方式，但是环境中噪声会大大降低智能设备中语音识别系统的性能。因此，如何降低噪声对人与人或人与系统之间的干扰，从带噪语音信号中恢复出纯净信号一直都是研究的热点和重点。
3.语音增强技术就是针对语音通信过程中的噪声干扰问题，抑制带噪语音信号中的噪声，从而获得纯净的语音信号。作为语音信号处理领域的经典研究课题，语音增强技术已经被许多学者广泛研究。语音增强算法也从最初的单通道语音增强算法扩展到麦克风阵列增强算法，由最初的时域研究扩展到变换域研究。目前，常用的传统语音增强算法有谱减法、相位补偿算法、子空间算法、听觉掩蔽效应算法等。自压缩感知(compressed sensing)理论被提出之后，针对稀疏信号的重构算法在语音信号处理领域得到了一定的发展。大多数的信号，包括语音信号，在某个特定的变换域上都是稀疏的。利用语音信号的稀疏性，压缩感知可以通过从带噪语音中重建稀疏分量来估计纯净语音的稀疏表示。基于压缩感知技术的语音增强现已经取得一些研究成果，但是，这些方法中还存在一些问题需要解决。目前的研究大多数是针对平稳的噪声，针对非平稳的噪声，重构性能不佳。此外，当输入带噪信号的信噪比较低，即噪声能量相对较大时，该方法的重构质量下降。综上所述，针对基于压缩感知的语音增强技术，如何提高重构算法的重构性能是一个亟待解决的科学问题。
技术实现要素：
4.本发明的目的在于克服压缩感知重构语音信号中的噪声干扰的问题，提出一种基于测量域噪声相减的压缩感知语音增强方法。该方法是在压缩感知算法重构之前，从测量域混合(语音和噪声)信号分量中减去噪声分量，以此来减少噪声的影响从而实现更加准确的语音信号重构。
5.为达到上述目的，本发明通过下述技术方案实现。
6.本发明提出了一种基于测量域噪声相减的压缩感知语音增强方法，该方法包括：
7.将带噪语音信号通过测量矩阵映射到测量域中得到测量域中的带噪语音信号，对测量域中不含语音信息的信号进行帧平均处理获得噪声估计信号，将测量域中的带噪语音
信号减去噪声估计信号获得纯净语音信号，再采用压缩感知算法对纯净语音信号进行重构及处理，获得增强的语音信号。
8.作为上述技术方案的改进之一，所述方法包括以下步骤：
9.s1.将带噪语音信号通过测量矩阵映射到测量域中；
10.s2.对测量域的每一帧信号逐帧判断是否含有语音信息：不含语音信息的信号为静音帧，含有语音信息的信号为语音帧；
11.s3.对于判断得到的第一个静音帧进行帧平均处理，获得噪声估计信号；
12.s4.继续判断下一帧是否含有语音信息，若下一帧为静音帧，则对该静音帧与上一帧静音帧共同进行帧平均处理并更新噪声估计信号；若下一帧为语音帧，则步骤s3获得的噪声估计信号为该语音帧的噪声估计信号，转至步骤s6；
13.s5.重复步骤s4，直到下一帧为语音帧；
14.s6.对下一帧的语音帧进行噪声相减得到纯净语音帧；
15.s7.采用压缩感知算法对语音帧进行重构；
16.s8.对重构语音帧的频谱进行逆短傅里叶变化得到增强的语音帧；
17.s9.重复步骤s3-s8，得到完整的增强的语音信号。
18.作为上述技术方案的改进之一，所述步骤s1包括以下步骤：
19.s1-1.将带噪语音信号进行语音分割；
20.s1-2.将带燥语音信号进行短时傅里叶变换映射到频域中，选择测量矩阵，将目标信号通过测量矩阵映射到测量域中。
21.作为上述技术方案的改进之一，所述步骤s1-1中，使用hamming窗函数对目标信号进行语音分割。
22.作为上述技术方案的改进之一，所述步骤s2中，采用语音活动检测技术判断测量域的每一帧信号是否含有语音信息。
23.作为上述技术方案的改进之一，所述步骤s3和s4中，对于语音信号，由于初始的几帧不包含语音信息，即判断为静音信号，故对初始静音帧的噪声估计信号进行估计，表达式为：
[0024][0025]
其中，为第i次测量的噪声估计信号，f表示帧索引，n
is
表示静音帧的长度，测量次数m为对于语音信号前几帧静音信号的测量次数，由测量矩阵φ决定，φ∈rm×n，r为实数，rm×n表示由实数构成的m行n列的矩阵，y
x
(i,f)表示第i次测量第f帧噪声语音x在测量域的投影，i＝1,2,......,m-1,m。
[0026]
作为上述技术方案的改进之一，所述步骤s4中对噪声估计信号不断更新，以及步骤s5中得到的最终噪声估计信号，均满足下式：
[0027][0028]
其中，为第j次测量的噪声估计信号，为第j-1次测量的噪声估计信
号，j＝m+1,m+2,......，y
x
(j,f)表示第j次测量第f帧噪声语音x在测量域的投影，ns表示静音帧的总数。
[0029]
作为上述技术方案的改进之一，所述步骤s6.对下一帧的语音帧进行噪声相减得到纯净语音帧，表达式为：
[0030][0031]
其中，表示测量域中第f帧的估计噪声信号，yf表示测量域与f帧的估计噪声信号对应的语音帧，且号对应的语音帧，且表示与测量域yf对应的带噪语音帧，φ表示测量矩阵，表示测量域噪声相减后的纯净语音帧。
[0032]
本发明与现有技术相比优点在于：
[0033]
1、针对基于压缩感知的语音增强技术，如何提高重构算法的重构性能是一个亟待解决的科学问题，本技术提出了从测量域混合信号中减去噪声分量，并对如何获得噪声分量进行了合理估计，以此来较少噪声的影响从而实现更加准确的语音信号重构，提高压缩感知算法的重构性能；
[0034]
2、基于测量域混合信号减去噪声分量的思想，可以有效的降低重构时噪声的干扰，实现更为准确的信号重构，因此，针对输入信号信噪比较低(噪声能量相对较大)的情况，本技术方法具有更好的重构性能。
附图说明
[0035]
图1为压缩感知的数学模型；
[0036]
图2为压缩感知的流程图；
[0037]
图3为本发明基于测量域相减的压缩感知语音增强方法流程图；
[0038]
图4为纯净语音和输入为0db白噪声时语音的波形图，其中，图4(a)为纯净语音(原始信号)的波形图，图4(b)为输入为0db白噪声时语音(带噪语音信号)的波形图；
[0039]
图5为采用两种算法(omp算法和本发明基于测量域噪声相减的omp算法)输出语音的波形图，其中，图5(a)为采用压缩感知正交匹配追踪(omp)算法输出的语音信号波形图，图5(b)为采用本发明基于测量域噪声相减的压缩感知正交匹配追踪(omp)算法输出的语音信号波形图；
[0040]
图6为添加不同信噪比白噪声的情况下两种算法输出的信噪比(snr)和感知语音质量评价(pesq)对比图，其中，图6(a)为添加不同信噪比白噪声的情况下两种算法输出的snr对比图，图6(b)为添加不同信噪比白噪声的情况下两种算法输出的pesq对比图；
[0041]
图7为添加不同信噪比粉红噪声的情况下两种算法输出的snr和pesq对比图，其中，图7(a)为添加不同信噪比粉红噪声的情况下两种算法输出的snr对比图，图7(b)为添加不同信噪比粉红噪声的情况下两种算法输出的感知语音质量评价pesq对比图。
具体实施方式
[0042]
下面结合附图给出本发明的实施流程，以详细说明本发明的技术方案。
[0043]
1、压缩感知基本原理
[0044]
压缩感知是一种根据信号的稀疏性而提出来的采样方法，它突破了奈奎斯特采样定理的限制，以远低于信号带宽两倍的速率进行采样，再使用优化算法可以将数据恢复出来。压缩感知的数学模型如图1所示，该模型可表示为：
[0045]
y＝φx
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0046]
其中，x表示长度为n的一维信号向量，x∈rn；y表示测量值，为长度为m的向量，y∈rm；φ表示测量矩阵，φ∈rm×n。压缩感知采样方法针对的信号是稀疏的，信号在稀疏域中表示为：
[0047]
x＝ψs
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0048]
其中，ψ为正交稀疏基矩阵，ψ∈rn×n；s是信号在稀疏基上的稀疏向量，s∈rn，其中有k《《n个非零值。结合公式(1)和公式(2)可得：
[0049]
y＝φψs
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0050]
压缩感知的主要目的是通过公式(3)逆向求解稀疏信号s，从而恢复目标信号x。该求解问题可以转化为求解最小化l0范数的问题，即：
[0051][0052]
然而求解l0范数最小化的问题是一个较为困难的np-hard的问题，candes和donoho等人指出，当测量矩阵满足有限等距性质(restricted isometry property,rip)条件时，求解l0范数最小化的问题可以转化为l1最小范数下的最优化问题，其求解模型为：
[0053][0054]
压缩感知的流程如图2所示，目标信号经过测量矩阵映射和稀疏变换后，采用压缩感知重构算法将目标信号恢复出来。现有的压缩感知语音增强方法是以测量域中混合(语音和噪声)信号分量来估计纯净的语音信号，当噪声能量增加时，会导致信号重构的质量下降。为了提高噪声能量较大情况下的压缩感知语音增强方法的性能，本文提出了一种基于测量域噪声相减的方法，即在信号重构之前，估计并减去测量域中的噪声信号分量，以此来实现更为准确的信号重构。
[0055]
2、基于测量域噪声相减的压缩感知语音增强方法
[0056]
在语音信号中，假设带噪信号为x，纯净语音信号为xs，加性噪声为xn，则带噪语音信号可表示为：
[0057]
x＝xs+xnꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0058]
对等式(6)同时进行短时的傅里叶(short-time fourier transform,stft)变换得：
[0059]
x＝xs+xnꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0060]
其中x、xs、xn分别表示频域上的带噪语音信号、纯净语音信号和噪声。将信号映射到压缩感知的测量域(综合公式(2)、(3)和(7))为：
[0061]
y＝φx＝φxs+φxn＝φψss+φψsn＝ys+ynꢀꢀꢀꢀꢀ
(8)
[0062]
其中ys、yn分别表示纯净语音和噪声在测量域中的向量，即xs、xn在测量域中的投影。由式(8)可知，在压缩感知算法进行信号重构时，选用的向量y是测量域混合信号分量(ys+yn)，其中噪声信号分量yn会造成重构质量下降。由于ys、yn是通过在测量矩阵φ上独立
投影得来的，所以为了提高重构质量，我们可以在测量域中减去噪声分量，选择测量域纯净语音分量进行压缩感知算法重构。
[0063]
基于测量域相减的压缩感知语音增强过程如图3所示，主要分为测量域映射和语音信号重构两个部分。首先，使用hamming窗函数对目标信号进行语音分割，对于分割的每一帧信号，进行stft变换到频域。同时，利用语音活动检测(voice activity detection,vad)技术来判断当前帧信号是否含有语音信息，并将整段信号分为语音帧和静音帧(不含语音信息)。语音活动检测(vad)又称语音端点检测，主要利用端点识别技术将有效的语音信号和静音信号或者噪声信号进行分离，使得后续的语音识别、语义识别、语音增强等工作更加高效。常见的语音活动检测的算法有双门限法、频谱熵法、自相关法、时频方差法等，本文使用双门限法进行语音端点检测。一般情况下，vad具有二进制输出。当特征值超过端点检测阈值，vad的输出为1，否则，vad的输出为0。测量域中的噪声是通过对静音语音信号进行帧平均来估计的。通常假设语音信号的前100-300ms不包含语音信息。本文假设前250ms信号为初始静音信号，则可以使用第一个n
is
静音帧来估计噪声，即：
[0064][0065]
其中，f表示帧索引，m表示测量次数，y
x
表示噪声语音x在测量域的投影。后续再检测到静音帧f(vad(f)＝0)时更新因此，测量域的帧平均噪声估计可表示为：
[0066][0067]
其中，ns表示静音帧的总数。
[0068]
在估计测量域中的噪声后，将vad掩码函数m应用于频谱x得：
[0069]
xm＝mx
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0070]
其中，掩码函数m为：
[0071][0072]
其中，b(b＝1
……
b)是频率区间索引，f(f＝1
……
f)是帧索引，乘法以元素的方法执行。对每个语音帧进行测量域噪声相减，如下所示：
[0073][0074]
其中，yf和分别表示测量域中f帧的信号y和估计噪声信号，且分别表示测量域中f帧的信号y和估计噪声信号，且表示测量域噪声相减后的f帧纯净语音信号。
[0075]
在测量域进行噪声相减后，应用压缩感知正交匹配追踪算法(orthogonal matching pursuit，omp)进行信号的重构，即：
[0076][0077]
其中，为估计稀疏信号的第f帧，ψ为稀疏基矩阵，本文选择傅里叶矩阵作作
birch canoe slid on the smooth planks”，文件名为“sp01.wav”，时长为2.67s，采样率为8000hz。
[0086]
叠加的噪声为白噪声和粉红噪声，对噪声进行下采样为8000hz。语音帧长n＝256，帧移为128。仿真软件为matlab 2021b。为了验证基于测量域噪声相减的omp算法的有效性，选用传统的omp算法(对比)和所提出的算法对带噪语音进行语音增强。
[0087]
图4为原始纯净语音和输入为0db白噪声时语音的波形图，其中，图4(a)为纯净语音(原始信号)的波形图，图4(b)为输入为0db白噪声时语音(带噪语音信号)的波形图。图5为omp算法和基于测量域相减的omp算法输出语音的波形图，其中，图5(a)为采用压缩感知正交匹配追踪(omp)算法输出的语音信号波形图，图5(b)为采用本发明基于测量域噪声相减的压缩感知正交匹配追踪(omp)算法输出的语音信号波形图。语音的波形图说明，传统的omp算法虽然能在一定程度上进行语音增强，但是对于噪声的抑制作用相对较差。本文提出的基于测量域相减的omp算法重构出的语音更加纯净和准确。对比传统的omp算法，所提出的算法在噪声抑制方面明显增强。
[0088]
对于语音增强的结果，常用的评价指标还有输出信噪比(snr)和感知语音质量评价(pesq)。为了进一步的验证基于测量域噪声相减的omp算法性能，本文分别计算了两种算法在不同信噪比下的白噪声和粉红噪声时的snr和pesq。图6表示分别添加-10db、-5db、0db、5db和10db信噪比白噪声的情况下，两种算法输出的snr和pesq，其中，图6(a)为添加不同信噪比白噪声的情况下两种算法输出的snr对比图，图6(b)为添加不同信噪比白噪声的情况下两种算法输出的pesq对比图；图7表示分别添加-10db、-5db、0db、5db和10db信噪比粉红噪声的情况下，两种算法输出的snr和pesq，其中，图7(a)为添加不同信噪比粉红噪声的情况下两种算法输出的snr对比图，图7(b)为添加不同信噪比粉红噪声的情况下两种算法输出的感知语音质量评价pesq对比图。图6、图7表明，在噪声类型不同，噪声程度不同的情况下，对比omp算法，基于测量域相减的omp语音增强算法重构出的语音质量更好。总体来说，本文所提出的方法具有更好降噪能力和重构性能。
[0089]
从上述对本发明的具体描述可以看出，本发明方法在压缩感知算法重构之前，从测量域混合(语音和噪声)信号分量中减去噪声分量，以此来减少噪声的影响从而实现了更加准确的语音信号重构。
[0090]
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
技术特征：
1.一种基于测量域噪声相减的压缩感知语音增强方法，其特征在于，该方法包括：将带噪语音信号通过测量矩阵映射到测量域中得到测量域中的带噪语音信号，对测量域中不含语音信息的信号进行帧平均处理获得噪声估计信号，将测量域中的带噪语音信号减去噪声估计信号获得纯净语音信号，再采用压缩感知算法对纯净语音信号进行重构及处理，获得增强的语音信号。2.根据权利要求1所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述方法包括以下步骤：s1.将带噪语音信号通过测量矩阵映射到测量域中；s2.对测量域的每一帧信号逐帧判断是否含有语音信息：不含语音信息的信号为静音帧，含有语音信息的信号为语音帧；s3.对于判断得到的第一个静音帧进行帧平均处理，获得噪声估计信号；s4.继续判断下一帧是否含有语音信息，若下一帧为静音帧，则对该静音帧与上一帧静音帧共同进行帧平均处理并更新噪声估计信号；若下一帧为语音帧，则步骤s3获得的噪声估计信号为该语音帧的噪声估计信号，转至步骤s6；s5.重复步骤s4，直到下一帧为语音帧；s6.对下一帧的语音帧进行噪声相减得到纯净语音帧；s7.采用压缩感知算法对语音帧进行重构；s8.对重构语音帧的频谱进行逆短傅里叶变化得到增强的语音帧；s9.重复步骤s3-s8，得到完整的增强的语音信号。3.根据权利要求2所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述步骤s1包括以下步骤：s1-1.将带噪语音信号进行语音分割；s1-2.将带燥语音信号进行短时傅里叶变换映射到频域中，选择测量矩阵，将目标信号通过测量矩阵映射到测量域中。4.根据权利要求3所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述步骤s1-1中，使用hamming窗函数对目标信号进行语音分割。5.根据权利要求2所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述步骤s2中，采用语音活动检测技术判断测量域的每一帧信号是否含有语音信息。6.根据权利要求2所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述步骤s3和s4中，对于语音信号，由于初始的几帧不包含语音信息，即判断为静音信号，故对初始静音帧的噪声估计信号进行估计，表达式为：其中，为第i次测量的噪声估计信号，f表示帧索引，n
is
表示静音帧的长度，测量次数m为对于语音信号前几帧静音信号的测量次数，由测量矩阵φ决定，φ∈r
m
×
n
，r为实数，r
m
×
n
表示由实数构成的m行n列的矩阵，y
x
(i,f)表示第i次测量第f帧噪声语音x在测量域的投影，i＝1,2,......,m-1,m。
7.根据权利要求2所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述步骤s4中对噪声估计信号不断更新，以及步骤s5中得到的最终噪声估计信号，均满足下式：其中，为第j次测量的噪声估计信号，为第j-1次测量的噪声估计信号，j＝m+1,m+2,......，y
x
(j,f)表示第j次测量第f帧噪声语音x在测量域的投影，n
s
表示静音帧的总数。8.根据权利要求2所述的基于测量域噪声相减的压缩感知语音增强方法，其特征在于，所述步骤s6.对下一帧的语音帧进行噪声相减得到纯净语音帧，表达式为：其中，表示测量域中第f帧的估计噪声信号，y
f
表示测量域与f帧的估计噪声信号对应的语音帧，且的语音帧，且表示与测量域y
f
对应的带噪语音帧，φ表示测量矩阵，表示测量域噪声相减后的纯净语音帧。
技术总结
本发明涉及一种基于测量域噪声相减的压缩感知语音增强方法，该方法包括：将带噪语音信号通过测量矩阵映射到测量域中，对测量域中不含语音信息的信号进行帧平均处理获得噪声估计信号，将测量域中的带噪语音信号减去噪声估计信号获得纯净语音信号，采用压缩感知算法对纯净语音信号进行重构及处理，获得增强的语音信号。本发明方法在压缩感知算法重构之前，从测量域混合(语音和噪声)信号分量中减去噪声分量，以此来减少噪声的影响从而实现了更加准确的语音信号重构。准确的语音信号重构。准确的语音信号重构。
