一种用于WASN的分布式语音噪声消除系统的制作方法

本发明涉及音频处理技术领域，尤其涉及一种用于wasn的分布式语音噪声消除系统。

背景技术：

在实际应用中，音频处理设备正常接收的语音信号经常会被各种噪声干扰，从而严重破坏接收到的语音信号质量，导致工作设备的输出语音性能下降。为了避免噪声给输出语音带来的不利影响，需要从含有干扰噪声的语音信号中提取出纯净的语音信号，其中提取纯净语音信号的方法也被称为语音噪声消除技术。从麦克风数量而言，语音噪声消除技术分为基于单通道的(单麦克风)和基于多通道的(多麦克风)。其中单通道因无法凭借单个麦克风获取空间信息而限制噪声消除后的语音性能；多通道的麦克风阵列技术虽然可以通过利用空间信息克服单通道的弊端，但是它只能在阵列结构规则(阵列几何信息已知)的情况下应用。

随着无线传感器技术的飞速发展，无线声传感器网络(wasn)的应用越来越广泛。因为wasn由独立的节点(每个节点可以是一个或多个麦克风传感器)组成，所以麦克风之间无法满足空间采样定理，从而导致已有的阵列技术不能直接应用到wasn中。尽管如此，wasn能够保证在同时利用时间和空间信息的前提下，可以克服阵列的一些局限性，所以用于wasn的分布式语音噪声消除技术开始兴起。现实生活中就可以利用wifi(或蓝牙)将多个智能手机或笔记本电脑构建成一个wasn。

在现有技术中研究了用于阵列的最小方差无失真响应算法，利用加权值控制噪声功率谱密度中非对角元素的能量值，通过广义线性坐标下降(generalizedlinearcoordinate-descent)算法执行节点之间的信息传递功能，提出了最小方差无失真响应算法分布式实现的技术方案。该技术虽然实现了分布式最小方差无失真响应算法，但是该技术对语音噪声消除后的噪声残留依然很严重，语音质量的感知评估(pesq)值和短期客观清晰度(stoi)值提高不大。

另外现有技术中存在研究了gossip算法的使用，提出一种分布式延迟求和波束形成语音噪声消除技术。并且该技术在wasn任意拓扑连接的情景下提出了改进的通用分布式同步平均(improvedgeneraldistributedsynchronousaveraging)方法交换每个节点处麦克风的数据，使每个节点的输出与含有数据处理中心的效果相同。虽然该技术提出了新的分布式算法，并且可以使最终的输出结果与含有数据处理中心实现的效果相同，但是该技术的输出效果和技术二的输出效果基本一样，性能较差。

在wasn没有数据处理中心的情况下，每个节点只能和附近的节点通信(通信半径之内的节点)，并且网络节点的能源有限，这就使得须用分布式算法实现信号的语音噪声消除，以期消噪后的效果能够达到把所有传感器的数据汇集到数据处理中心统一处理的效果(含有数据处理中心的算法不能直接应用到wasn中)。现存的分布式语音噪声消除技术，有一些是无法达到含有数据处理中心的输出效果，还有一些虽然达到了含有数据处理中心的输出效果，但是每个节点麦克风的输出性能并不是很高，噪声残留依然很大。

技术实现要素：

根据现有技术存在的问题，本发明公开了一种用于wasn的分布式语音噪声消除系统，

相位对齐模块，确定每个节点到声源的距离，将距离声源最远的节点定义为参考节点，将其余节点接收的信号与参考节点接收的信号进行相位对齐获得同相位节点信号；

离散傅里叶变换模块，将相位对齐模块传送的各个节点信号分别进行分帧加窗处理、对每一帧信号进行离散傅里叶变换获得离散频谱信号；

语音活动检测模块，接收离散傅里叶变换模块传送的离散频谱信号、通过离散频谱信号进行语音活动检测、判断每一帧信号有无语音；

噪声功率谱密度估计模块，接收所述语音活动检测模块传送的检测结果、根据无语音帧信号的离散频谱信息计算噪声功率谱密度；

分布式参数多通道维纳滤波模块，接收离散傅里叶变换模块传送的离散频谱信号以及噪声功率谱密度估计模块传送的噪声功率谱密度信息、采用分布式参数多通道维纳滤波方法获得分布式参数多通道维纳滤波器系数；将分布式参数多通道维纳滤波器系数与离散频谱信号组合形成输出信号yp；

分布式算法迭代模块，接收分布式参数多通道维纳滤波模块传送的输出信号yp，将输出信号yp处理为求平均值的形式、根据metropolis权矩阵多次迭代求得初始状态值的平均值得到每个节点的输出信号yp；

离散傅里叶反变换模块，接收分布式算法迭代模块传送的输出信号yp，通过对输出信号yp进行离散傅里叶反变换得到时域当前帧输出语音信号，并将时域每帧输出信号进行重叠相加得到最终的输出信号。

作为优选的方式，所述分布式参数多通道维纳滤波器系数采用如下方式获取：

其中，h是分布式参数多通道维纳滤波器系数，[]^t表示向量或矩阵的转置，是由δi取倒数得到的，α是该算法中的参数分别为1、3、5，|xi|²表示信号功率谱密度。

由于采用了上述技术方案，本发明提供的一种用于wasn的分布式语音噪声消除系统，该系统通过对用于阵列的参数多通道维纳滤波器系数进行了修改，使得噪声消除后的语音信号性能甚至优于修改前用于阵列输出的语音信号性能。在用于阵列的参数多通道维纳滤波算法的基础上，提出一种可以用于wasn的分布式语音噪声消除技术，它是一种可以应用在任意拓扑网络连接下的分布式语音噪声消除技术，本方法通过使用metropolis权矩阵，在一定迭代次数下使每个节点处麦克风的输出信号均与含有数据处理中心的结果一致，并具有一定的鲁棒性。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明系统的结构示意图；

图2为本发明中无线声传感器网络示意图；

图3为本发明实施例中每种方法语音噪声消除后的stoi值示意图：图3(a)为无混响；图3(b)为混响时间为300ms；

图4为本发明实施例中每种方法语音噪声消除后的pesq值：图4(a)为无混响；图4(b)为混响时间为300ms

具体实施方式

为使本发明的技术方案和优点更加清楚，下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述：

如图1所示的一种用于wasn的分布式语音噪声消除系统，包括相位对齐模块、离散傅里叶变换模块、语音活动检测模块、噪声功率谱密度估计模块、分布式参数多通道维纳滤波模块、分布式算法迭代模块和离散傅里叶反变换模块。

其中相位对齐模块的作用是确定每个节点到声源的距离，将距离声源最远的节点定义为参考节点，将其余节点接收的信号与参考节点接收的信号进行相位对齐获得节点信号。

作为优选的方式，相位对齐模块的工作原理是：在wasn中，在距声源某个已知距离为d的位置放一个参考麦克风，利用该麦克风接收到的信号能量和其它节点处麦克风接收到的信号能量可以估计出wasn中每个节点到声源的距离di。其中角标i＝1,2,...,i，i表示wasn中节点的数量。距离估计公式如下：

其中，e和ei分别是参考信号和wasn中每个节点处麦克风信号的能量，ε和εi是背景噪声的能量，能量公式如下：

其中，n是节点处麦克风接收信号对应的总采样点数，fs是采样频率，即一秒信号对应的点数。该公式利用语音第一秒大部分是无语音段的特点，估计背景噪声的能量。

由上述方法确定出每个节点到声源的距离后，定义距离声源最远的那个节点为参考节点，此时该节点的输入信号定义为x″a(n)，其余节点的待对齐输入信号为x″b(n)。令x″b(n)循环通过一个单位延迟器，并且同时与x″a(n)做互相关运算，表达式如下：

rab(τ)＝e[x″a(n)x″b(n-τ)],τ＝0,1,...,t(3)

其中，t是平移最大量，可适当选取。当τ的取值使以上两路信号对齐时，互相关函数值最大。令find{}为取最大值对应的τ值运算，则待对齐信号与参考信号对齐的输出信号为：

进一步的，离散傅里叶变换模块的作用是将相位对齐模块传送的各个节点的信号分别进行分帧加窗处理、对每一帧信号进行离散傅里叶变换获得离散频谱信号。

作为优选的方式，离散傅里叶变换模块的工作原理是接收相位对齐模块传送的节点信号，对各路信号分别进行分帧加窗处理，对每一帧信号进行离散傅里叶变换dft，具体实施中在验证时语音信号的采样频率fs是16khz，用的是汉宁窗，帧移为50％，每帧数据长度为m＝320点。其中汉宁窗的表达式如下：

ω(m)＝0.5-0.5cos(2πm/m),m＝0,1,...,m-1(5)

根据汉宁窗表达式可以得到加窗后的信号为：

xi(m)＝x′i(n)ω(m)(6)

然后对各路信号加窗后的每一帧信号进行dft，变换后可得离散频谱为：

其中，k表示频点索引，l表示当前帧。

语音活动检测模块的作用是：接收离散傅里叶变换模块传送的离散频谱信号、通过离散频谱信号进行语音活动检测、判断每一帧信号有无语音。

作为优选的方式：语音活动检测模块进行语音活动检测时：同样利用语音第一秒大部分是无语音段的特点，并结合分帧加窗的处理过程，语音信号最初始的无语音帧的数量为nis帧，其中nis＝fs/(50％×m)-1＝99。那么，利用这nis帧估计的噪声平均频谱为：

式(8)表示每一帧信号的对应频点先求和再取平均。进一步，噪声帧的对数频谱估计有如下表示：

其中|·|是取模运算。然后，计算每帧信号的对数频谱：

由公式(9)和公式(10)可以得出每帧信号与噪声信号的对数频谱距离，对数频谱距离公式如下：

综上，可以得出语音活动检测的判断方法：首先，设置一个无语音段计数器，该计数器初始值可以设置为100，同时设置一个对数谱距离阈值为3。然后计算每一帧信号与噪声帧的对数谱距离dspec，判断dspec是否小于该对数谱距离阈值，如果是，则该帧是无语音帧，无语音段计数器加1，如果不是，则该帧是有语音帧，且无语音段计数器不论是多少，必须归零。最后还需要注意的是如果无语音段计数器在归零之前的值小于最小无语音长度，则认为无语音段计数器上一次归零后到这一次归零前是无语段的帧均为有语音帧。这里令最小无语音长度为10。

噪声功率谱密度估计模块的作用是接收所述语音活动检测模块传送的检测结果、根据无语音帧信号的离散频谱信息计算噪声功率谱密度。

作为优选的方式，噪声功率谱密度仅在无语音帧进行更新。每个节点处的噪声功率谱密度更新公式如下：

δi＝(1-β)|xi,noise(k,l)|²+β|xi,noise(k,l-1)|²(12)

其中，β＝0.997，δi表示第i个节点的噪声功率谱密度估计值，它对应到每个频点都会有一个估计值。如果当前帧是噪声帧，该值就会被上式更新。|xi,noise(k,l)|²表示当前第l帧是噪声帧对应的频点模值的平方。

进一步的，分布式参数多通道维纳滤波模块作用是接收离散傅里叶变换模块传送的离散频谱信号以及噪声功率谱密度估计模块传送的噪声功率谱密度信息、采用分布式参数多通道维纳滤波方法获得分布式参数多通道维纳滤波器系数；将分布式参数多通道维纳滤波器系数与离散频域信号组合形成滤波信号。

其中具体的计算方式是：

其中，h是一个向量，即分布式参数多通道维纳滤波器系数；由于δi和|xi|²的存在，它的每一个向量值都会对应到具体的频点，[]^t表示向量或矩阵的转置，是由δi取倒数得到的，α是该算法中提到的参数，在本专利中，该参数值的取值分别为1、3、5。|xi|²表示信号功率谱密度，同δi一样，每个频点都会被更新，更新公式如下：

其中，l表示当前帧。上式在每一帧都要更新，即不论是有语音帧还是无语音帧都要更新。根据式(13)可以得到第p个节点的输出信号yp(即该模块的输出)的表达式：

其中，[]^h表示向量或矩阵的共轭转置，x＝[x1(k,l),x2(k,l),...,xi(k,l)]^t。

分布式算法迭代模块作用是接收分布式参数多通道维纳滤波模块传送的滤波信号，将滤波信号处理为求平均值的形式、根据metropolis权矩阵多次迭代求得初始状态值的平均值得到每个节点的输出信号yp。作为优选的方式：实现之前需要先将yp的表达式写成求平均值的形式：

其中，和观察式(16)可知dpmwf-α的结果只需要每个节点的麦克风获得所有节点麦克风初始状态值的平均值，就可以获得和上式同样的输出结果。在分布式算法下，是通过令每个节点之间交换特定数据的方式，将初始状态值不断迭代更新求得该初始状态值的平均值，迭代公式如下：

其中，ξ(t)＝[ξ1(t),ξ2(t),...,ξi(t)]^t，ζ(t)＝[ζ1(t),ζ2(t),...,ζi(t)]^t，t表示迭代次数。w是metropolis权矩阵，定义如下：

式(18)中e代表两个不同节点处麦克风可以互相通信的连接集合，即(i,j)∈e(i,j＝1,2,…,i,i≠j)。ηi表示第i个节点可以与附近节点通信的数量。在上述迭代计算使得每个节点处麦克风的输出信号收敛时，输出信号的结果可以达到含有数据处理中心的解。本专利在验证时设置迭代次数上限为100次，当迭代到上限时默认收敛。

离散傅里叶反变换模块的作用是接收分布式算法迭代模块传送的输出信号yp，通过对输出信号yp进行离散傅里叶反变换得到时域当前帧输出语音信号，并将时域每帧输出信号进行重叠相加得到最终的输出信号。作为优选的方式：进行idft得到时域当前帧输出语音信号yp(m,l)。idft公式如下：

因为本专利在模块2对各路信号进行了分帧加窗处理，且帧移是50％，所以从得到第一帧输出语音信号yp(m,1)时，就要与第二帧输出语音信号yp(m,2)进行重叠相加运算，重叠部分占50％，具体公式如下：

其中[·]是取整运算，[a]表示不超过数a的最大整数。

本发明公开的一种用于wasn的分布式语音噪声消除系统，为验证本发明方法的有效性，通过imgae模型模拟了一个5×5×3的封闭式房间，该房间分为无混响和混响时间为300ms两种情况。在该wasn中随机分布了10个节点，每个节点是1个麦克风，声源分别设置了5个不同的位置，节点与声源的高度均设置为1米。仿真模拟的二维wasn如图2所示，节点之间的通信距离上限设为2.2米。

声源是从timit数据库[https://download.csdn.net/download/sdhyfxh/4086482]随机选取的一个6秒纯净语音信号，采样频率是16khz。在每个节点处麦克风接收到的语音信号同时加上了一个互不相关的高斯白噪声作为输入噪声信号，该噪声可导致节点接收信号的信噪比在5db左右。

此时，采用本系统所提出的dpmwf-α(其中α值分别为1、3、5)语音噪声消除技术对每个节点处麦克风接收信号进行降噪，同时采用文献[1]和文献[3]中的方法对该实验也进行了语音噪声消除。实验结果表明，不论哪种方法都能使wasn中每个节点的输出结果达到一致。图3、图4给出了这三种方法在声源位置分别位于ⅰ、ⅱ、ⅲ、ⅳ、ⅴ时的性能比较。其中图3表示的是分别在无混响和有混响条件下，语音噪声消除后stoi值的性能比较，图4表示的是分别在无混响和有混响条件下，语音噪声消除后pesq值的性能比较。可以看出不论是在无混响条件下还是有混响条件下，亦或是声源位于哪个位置，在stoi值和pesq值性能方面，本专利所提方法和文献[1]、文献[3]中的方法相比，要占据绝大优势。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。

附件：

[1]a.bertrand,j.callebautandm.moonen,"adaptivedistributednoisereductionforspeechenhancementinwirelessacousticsensornetworks,"inproc.oftheinternationalworkshoponacousticechoandnoisecontrol(iwaenc),telaviv,israel,aug.2010.

[2]r.heusdens,g.zhang,r.c.hendriks,y.zengandw.b.kleijn,``distributedmvdrbeamformingfor(wireless)microphonenetworksusingmessagepassing,”presentedattheemiwaenc2012；internationalworkshoponacousticsignalenhancement,aachen,germany,2012,pp.1-4.

[3]y.zengandr.c.hendriks,"distributeddelayandsumbeamformerforspeechenhancementviarandomizedgossip,"inieee/acmtransactionsonaudio,speech,andlanguageprocessing,vol.22,no.1,pp.260-273,jan.2014.