[欢迎访问我的博客原文](单通道语音增强之统计信号模型)
1. 信号估计理论简述
信号估计理论是现代统计处理的基础课题[@ZhangXianDa2002ModernSP]，在通信、语音、图像领域均有广泛应用。语音增强，就是从带噪的语音测量信号中估计原始的无噪语音，这是典型的信号估计问题。 《语音增强–理论与实践》[@loizouspeech]一书中列举了用于语音增强的一系列统计模型。
假设麦克风采集到的带噪语音序列为
，并且噪声都是加性噪声。则带噪语音序列为无噪语音序列与噪声序列的和。原始语音信号与噪声均可视为随机信号。
在时域对
进行估计是非常困难的，通过傅立叶变换，我们可以将信号分解为频域上互相独立的系数。信号估计模型转变为对每一个频点的系数进行估计的模型，不同频点之间的参数是相互独立的。
这个方法就叫做统计信号谱分析（Statistical Spectral Analysis)。显然地，纯净信号谱
带有幅度与相位两参数，我们实际上是对幅度 和相位 进行参数估计。
重组谱幅度和谱相位估计值即可恢复纯净语音谱，估计值用上标来表示：
。
实际信号的幅度和相位是不方便直接用在运算过程中的，因为信号取值范围不定，且瞬时变化。在噪声抑制领域，更常用的语音谱估计方法是对抑制增益(Suppression Gain)进行估计，不同的估计准则称为抑制准则(Supression Rule)。
通常会根据先验信噪比、后验信噪比来估计抑制增益
。 并且可以在只有噪声出现的时刻更新 ， 在语音存在的时刻进行抑制，无须每帧去调用噪声抑制算法，计算过程比直接估计信号谱灵活。
综上，语音增强的典型流程就是：
对带噪语音 分帧， 每一帧进行DFT得到 。 估计或者沿用上一帧的抑制增益 ，得到纯净语音谱 。 对 进行IDFT,得到纯净语音序列的估计 。
为了估计模型的建模，对测量信号、估计信号、噪声信号都需要作一些数学上的假设和简化。其中对噪声一般会作以下假设：
噪声是与语音独立的加性噪声； 每一帧噪声的统计分布是稳态的； 噪声的傅立叶级数是零均值复高斯分布。
2. 最大似然估计ML
基如果不考虑信号的先验分布，即认为信号值是确定信号，而不是随机信号，我们只需要分析含有信号
为参数的带噪信号 的概率分布 ，并使之最大。这种估计方法叫最大似然估计器（Maximum Likelihood Estimation）。将 的带参概率分布 称为似然函数（Likelihood Function）。对纯净信号 的估计，表达为求解合适的 值，使得似然函数 最大。
文献[@mcaulay1980speech]最早将最大似然估计法用在语音增强领域。对于纯净语音，可以假设纯净语音幅度
和相位 是未知但确定，无需考虑其先验概率分布。最大似然语音增强模型表达为：
把
代入噪声零均值复高斯分布公式中，得到：
3.贝叶斯估计
如果比最大似然估计更进一步，考虑待估计量
也是随机变量，且的先验分布为 ，这种假设下的估计方法叫做贝叶斯估计[@ZhangXianDa2002ModernSP]。定义估计值 与实际值 之间的误差函数为 ，贝叶斯估计器的目标即为找出是平均误差 最小的估计值 。
对于待估计的纯净语音谱，贝叶斯估计器可以表达为：
误差
的期望值取决于待测信号与测量信号的联合概率分布。
对联合概率密度进行分解，提取
与条件概率密度 ：
因为估计值
与 相互独立，对 的求解可以等价于对 求极大值，也就是：
误差函数
的计算模型，会引出不同种类的估计器。典型的误差函数有几种类型[@ZhangXianDa2002ModernSP]： 平方误差函数，对应最小均方估计。 绝对值误差函数，对应条件中位数估计。 均匀误差函数，对应最大后验估计。
3.1 最小均方估计（MMSE）
最小均方估计(Minimum Mean Square Error Estimation)[@ephraim1985speech]使用平方误差函数，平均误差为：
为求平均误差的极大值，可对估计量
求导，并求极值点。
显然极值点的
为被估计量 的条件均值：
亦可以使用贝叶斯公式展开，得到：
在语音增强的模型里，纯净语音谱的估计为其在带噪语音谱下的条件均值。
上式中，为得到纯净语音谱需要分别估计谱幅度
和谱相位 。
同时估计谱幅度和谱相位是很难的，研究者提出了许多分别估计谱幅度和谱相位的方法，估计完成后再用两者重组复语音信号。
3.1.1 MMSE谱幅度估计
最小均方根估计器MMSE short-time spectral amplitude
def mmse_stsa_gain(parameters=None):gamma = parameters['gamma']ksi = parameters['ksi']vk = ksi * gamma / (1 + ksi)j0 = np.i0(vk / 2)j1 = np.i1(vk / 2)A = sqrt(pi * vk) / 2 / gammaB = (1 + vk) * j0 + vk * j1C = np.exp(-0.5 * vk)gain = A * B * Creturn gain
3.1.2 MMSE对数谱幅度估计
对数最小均方根估计器The MMSE log spectral amplitude (MMSE-LSA)， 或者缩写为LogMMSE估计器。
def logmmse_gain(parameters=None):gamma = parameters['gamma']ksi = parameters['ksi']A = ksi / (1 + ksi)vk = A * gammaei_vk = 0.5 * expn(1, vk)gain = A * np.exp(ei_vk)return gain
3.1.3 MMSE平方谱幅度估计
频谱幅度平方估计器MMSE magnitude squared[@wolfeefficient]
def mmse_sqr_gain(parameters=None):gamma = parameters['gamma']ksi = parameters['ksi']vk = ksi * gamma / (1 + ksi)j0 = np.i0(vk / 2)j1 = np.i1(vk / 2)A = ksi / (1 + ksi)B = (1 + vk) / gammagain = sqrt(A * B)return gain
3.2 最大后验估计 MAP
当贝叶斯估计器采用均匀误差函数时，平均误差为：
显然要使得平均误差最小，就是要求目标估计
，使得 最大。这种估计模型称作最大后验估计(Maximum A Posteriori Estimation)。这个模型的意思是只有估计值 与原始值 相等，误差才为0，其他时候误差均匀为1。估计值可以表达为：
在语音增强的模型里，纯净语音谱的估计为其在带噪语音谱下的条件均值。
文献[@wolfeefficient]提出了两种基于最大后验估计(Maximum A Posteriori Estimation)的语音增强算法。一种是同时求解幅度和相位的混合最大后验估计： 另外一种是单纯估计幅度的方法，两种估计器最后的噪声抑制增益略有不同。
3.2.1 幅度和相位混合最大后验估计
[@wolfeefficient]
def map_joint_gain(parameters=None):gamma = parameters['gamma']ksi = parameters['ksi']eps = 1e-6gain = (ksi + sqrt(ksi^ 2 + 2 * (1.0 + ksi)* ksi/ (gamma + eps))) / 2.0/ (1.0 + ksi)return gain
3.2.2 纯幅度最大后验估计
[@wolfeefficient]
def map_sa_gain(parameters=None):gamma = parameters['gamma']ksi = parameters['ksi']eps = 1e-6gain = (ksi + sqrt(ksi^ 2 + (1.0 + ksi)* ksi/ (gamma + eps))) / 2.0/ (1.0 + ksi)return gain
参考文献
[1] 张贤达, 现代信号处理. 清华大学出版社有限公司, 2002.
[2] P. C. Loizou, Speech enhancement: theory and practice. CRC press, .
[3] R. McAulay and M. Malpass, Speech enhancement using a soft-decision noise suppression filter," IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 28, no. 2, pp. 137-145,1980.
[4] Y. Ephraim and D. Malah, Speech enhancement using a minimum mean-square error log-spectral amplitude estimator," IEEE transactions on acoustics, speech, and signal processing,vol. 33, no. 2, pp. 443-445, 1985.
[5] P. J. Wolfe and S. J. Godsill, Efficient alternatives to the ephraim and malah suppression rule for audio signal enhancement," EURASIP Journal on Applied Signal Processing, vol. , pp.1043-1051, .
如果觉得《功率谱 魏凤英统计程序_单通道语音增强之统计信号模型》对你有帮助，请点赞、收藏，并留下你的观点哦！
