语音增强原理之增益因子

版权声明：本文来自于icoolmedia，欢迎转载，转载请注明出处，相关的音视频算法讨论请加群（374737122）一起交流！

　　上次关于语音增强的原理讲说了噪声估计问题，这次打算说下增益因子如何确定，也就是当噪声已知后，如何进行去噪的问题（把增益因子与带噪语音相乘即可）。这里主要说下MMSE滤波，顺带说下谱减法、维纳滤波。当然也有其它方式来实现语音增强的，比如基于矩阵分解原理的子空间法、基于自适应滤波器的降噪，有的方法icoolmedia比较清楚，有的也还在学习之中，同时也欢迎各位朋友就不足之处批评指正。

　　先说下经典的谱减法。我们还是假设带噪语音y(n)由纯净语音x(n)和加性噪声d(n)组成，其时域表示与频域表示为

　　在语音增强领域中，最常使用的是频域功率谱形式，那么，带噪语音的功率谱可以表示为Y与其共轭相乘，展开可得

　　展开式的第三项被称为交叉项，当纯净语音与加性噪声不相关时，交叉项为0，那么，当我们已经估计出噪声信号的功率谱时，纯净语音信号的估值就可以表示为

　　根据线性滤波理论，可以将这个滤波过程建模为

　　综合以上两式，H可以表示为

　　这里的H，就是线性滤波系统的传递函数，在语音增强领域，通常也称为增益函数/抑制函数、或者增益因子/抑制因子，都是是同一个意思。另外，要注意，上面的这个增益因子表示只是一个理想的过程。因为当从带噪语音中减去估计噪声后，总会遗留一些或长或短的小谱峰，这些谱峰比较影响听感。这种现象就是出现了音乐噪声。因此，如果谱减法要实际使用时，必须做如下改变。

　　经过这两方面的改变，谱减法具有如下形式：

 　　这里当存在语音时，用第一个式子，当没有语音存在时，就用下面的式子，其中，alpha就是过减因子，取为一个大于1的值，beta为一个远小于1的值，具体取值范围请参考Berouti等人的论文：Enhancement of speech corrupted by acoustic noise，这里不再详加叙述。

　　这时，当存在语音时谱减法的增益因子就变为

　　这里的gamma是后验信噪比，为带噪语音与噪声的功率之比。谱减法增益因子就说完了。频域维纳滤波的增益因子可以参考我以前写的博客一个频域语音降噪算法实现及改进方法中的内容，里面有详细的推导过程，这里就不详加叙述了。下面重点说下MMSE降噪算法是如何确定增益因子的。

　　MMSE估计器用在语音增强之中，就是在贝叶斯准则下估计出来的纯净语音频幅度与实际幅度的均方误差最小，因此，也可以称为贝叶斯MSE。而要做到这一点，我们可以充分利用带噪语音的先验信息来提高估计的准确性。即，假定我们已知信号的噪声DFT系数的概率密度，在此基本上，充分利用这种已知的先验信息，提高估计的准确性。因此，贝叶斯MSE用公式表示如下：

　　我们来推导下使贝叶斯MSE最小的估计量，首先应用贝叶斯原理，联合概率密度可以写为：

所以

对中括号中的积分求导

令等式等于0，得

　　在此我们假设傅里叶变换系数之音是统计独立的。因此上式可以表示为

　　可以看到，要想得到MMSE估计器，我们首先需要计算纯净语音第k个分量的后验概率密度函数，它可以通过贝叶斯准则得到：

　　这里x_k是随机变量X_k的实际值。把上面这个后验概率密度函数表达式代入我们推导出来的MMSE估计器中

 其中

代入MMSE估计器中，我们最终得到MMSE幅度谱估计器（推导过程请参考：语音增强-理论与实践中的附录B）

其中，Γ(.)为伽马函数，Φ（a,b;c)为合流超几何函数，ξ为先验信噪比、最后一个式子为后验信噪比。

最后，把合流超几何函数写成贝塞尔函数的形式，我们就得到了最终的MMSE估计器的表达式：

　　如果我们定义：

的话，这里G就是我们要求的MMSE幅度估计器的增益。

　　另外想说一下，MMSE估计的推导思路我弄明白了，主要是通过参考《语音增强-理论与实践》、《统计信号处理基础-估计与检测理论》这两本书做到的，但关于合流超几何函数与贝塞尔函数的推导内容还没完全搞明白，如果不是对理论推导过程非常感兴趣的话，这里也没有必要深究，只要会使用这个结果就行了。

　　使用MMSE做语音增强，经典的出处应该是Speech enhancement using minimum mean-square error这篇论文，但里面讲的并不详细，这里尽可能的给出能让大家理解流程的推导。当然，如果感兴趣的话，icoolmedia还是推荐大家最好都认真看一遍上面提到的资料。