本发明提供的是一种基于贪婪深度字典学习的单信道语音去混响方法。步骤一、将输入的混沌、纯净语音进行信号预处理，进行特征的提取；步骤二、结合稀疏和贪婪深度字典学习对单信道语音去混响问题进行模型建立；步骤三、将建立好的模型与提取的特征进行模型训练，得到基本系数；步骤四、再次输入混沌、纯净语音信号数据进行测试，经过傅里叶逆变换后，最终得到纯净语音。本发明将传统的单信道语音去混响方法与深度字典学习相结合，提高单信道去混响的效果，使其具有良好的去混响效果，提高深度网络结构的可推导性。
1.一种基于贪婪深度字典学习的单信道语音去混响方法，包括如下步骤：步骤一、将输入的混沌、纯净语音进行信号预处理，进行特征的提取；
步骤二、结合稀疏和贪婪深度字典学习对单信道语音去混响问题进行模型建立；
步骤三、将建立好的模型与提取的特征进行模型训练，得到基本系数；
步骤四、再次输入混沌、纯净语音信号数据进行测试，经过傅里叶逆变换后，最终得到纯净语音；
其特征是所述的结合稀疏和贪婪深度字典学习对单信道语音去混响问题进行模型建立具体包括：
第1、对单信道语音去混响进行建模，具体描述为：(1)将n时刻的语音信号利用下式进行描述：其中D是原始语音，nk是k次反射后的延时单元数，ρk是第k次反射时的发射系数，s(n-knk)是第k次反射后的语音信号，u(n)表示周围环境的噪声，x(n)表示混响语音；
(2)将环境噪声u(n)去掉简化得：(3)根据单位冲激函数的卷积性质，改为卷积的形式：其中，h(n)是冲击响应,*是卷积操作符号，δ(n)表示单位冲激函数，h(n)反应混响的物理特征、是位于knk并且都衰减ρk倍后的许多冲击函数δ(n)的叠加，房间脉冲响应h(n)和纯净语音s(n)卷积后获得混响语音x(n)；
第2、将单信道语音去混响模型转化为基于稀疏表示的单信道语音去混响模型，其转化过程如下：
(1)通过稀疏表示方法对单信道语音去混响进行建模分析，在忽略环境噪声的情况下，利用下式简单描述语音去混响问题：Y＝HX
(2)根据X＝AS，假定纯净语音含有正交基稀疏矩阵表示矩阵A，则单信道语音去混响的稀疏表示为：
第3、引入深度语音字典，将字典学习问题转化为以下两种表示方法，任选其一：方法1，将求解目标函数定义为：X＝DZ
其中X是纯净语音，D是语音字典，Z是稀疏表示矩阵；
字典学习问题表示为求解目标函数：方法2，采用K-SVD字典学习算法，将求解目标函数定义为：其中zi是矩阵Z中的第i列，||.||0表示zi中非零元素的个数，τ表示稀疏程度，用于控制Z的每一列中所含非零元素的最大个数；
第4、引入深度语音字典后，进行求解，具体描述为：(1)根据单层字典学习进行求解，在求解过程中利用稀疏编码和字典更新，具体描述如下：
1)稀疏编码，运用LARC方法，已知字典D(t-1)，在第t次字典学习迭代中稀疏矩阵的每个列向量更新如下表示：
其中X(:,n)表示给定的观测矩阵，D(t-1)表示字典矩阵， 表示满足稀疏约束条件的稀疏表示矩阵；
2)字典更新，根据所获得稀疏表示矩阵Z，从训练信号集X中寻找当次迭代中更适合训练信号的字典D，字典更新如下式表示：
3)针对每一个原子 残差表示为：式中将DZ分解为L个秩为1的矩阵(d(:,l)z(l:,))之和，Rl代表除第l个原子之外的原子所带来的残差。通过SVD分解找到与Rl最为接近的秩为1的矩阵；
4)进行预处理，定义Λl为下式：Λl＝{i|1≤i≤K,c(l,:)(i)≠0}
5)定义M为D×|Λl|的矩阵，其中(Λl(i),i)设置为1，其他设置为0，则：上式 代表R(l)保留z(l,:)中非零元素所对应的原子所构成的矩阵，代表z(l,:)只保留非零元素所构成的向量，字典更新如下式表示：
7)使用U中的第一列来跟新d(:,l)，V的第一列乘以Δ(l,l)来更新相应的系数，完成一个原子的更新；
2)采用逐层训练进行逐层分解，首先训练除第一特征层Z1，字典D1：X＝D1*Z1
其次对Z1进行分解，求解第二层字典D2、特征Z2：Z1＝D2*Z2
3)基于贪婪深度字典学习的单信道语音去混响模型由此建立。
一种基于贪婪深度字典学习的单信道语音去混响方法
技术领域
[0001] 本发明涉及的是一种单信道去混响方法。
背景技术
[0002] 混响是指在声学上，除延迟时间达到约50ms以上的反应波称为回声外，其余的反射波产生的效应统称为混响。随着互联网发展及智能设备的普及，语音或声音信号的可理
解性与可感知性对于人机交互的用户体验变得尤为重要。然而混响严重影响了语音识别的
准确率，降低了用户智能交互的体验，同时增加了有听力障碍者对于语音理解的困难度。
[0003] 语音去混响的目的是抑制或移除混响语音中的混响成分，尽可能的实现源信号的恢复。根据系统麦克风的数量分为：单信道语音去混响和多信道语音去混响。其中单信道语
音去混响利用接收位置一点的语音信号时频域特性作为先验信息对混响信号进行处理，具
有系统相对简单，计算复杂度低以及成本较为低廉等优势。因此，单信道语音去混响问题具
有重要的研究价值。
[0004] 传统的单信道语音去混响的方法大多通过先建立语音去混响模型，通过加入不同的约束条件在非适定问题中求出可确定的解，从而获得符合条件的估计信号。尽管在求解
估计信号时具有可解释性，但这些方法都存在适用范围不足，实际效果不佳等缺点。基于深
度学习的语音去混响方法通过训练大量样本，直接学习混响语音和纯净语音之间的非线性
映射关系。虽然能取得很好的去混响结果，但具有数据获得困难和参数冗余的缺点，导致其
在实际应用中性能下降。
发明内容
[0005] 本发明的目的在于提供一种能够提高单信道去混响效果的基于贪婪深度字典学习的单信道语音去混响方法。
[0006] 本发明的目的是这样实现的：
[0007] 步骤一、将输入的混沌、纯净语音进行信号预处理，进行特征的提取；
[0008] 步骤二、结合稀疏和贪婪深度字典学习对单信道语音去混响问题进行模型建立；
[0009] 步骤三、将建立好的模型与提取的特征进行模型训练，得到基本系数；
[0010] 步骤四、再次输入混沌、纯净语音信号数据进行测试，经过傅里叶逆变换后，最终得到纯净语音。
[0011] 所述的结合稀疏和贪婪深度字典学习对单信道语音去混响问题进行模型建立具体包括：
[0012] 1、对单信道语音去混响进行建模，具体描述为：
[0013] (1)将n时刻的语音信号利用下式进行描述：
[0014]
[0015] 其中D是原始语音，nk是k次反射后的延时单元数，ρk是第k次反射时的发射系数，s(n-knk)是第k次反射后的语音信号，u(n)表示周围环境的噪声，x(n)表示混响语音；
[0016] (2)将环境噪声u(n)去掉简化得：
[0017]
[0018] (3)根据单位冲激函数的卷积性质，改为卷积的形式：
[0019]
[0020] 其中，h(n)是冲击响应,*是卷积操作符号，δ(n)表示单位冲激函数，h(n)反应混响的物理特征、是位于knk并且都衰减ρk倍后的许多冲击函数δ(n)的叠加，房间脉冲响应h(n)
和纯净语音s(n)卷积后获得混响语音x(n)；
[0021] 2、将单信道语音去混响模型转化为基于稀疏表示的单信道语音去混响模型，其转化过程如下：
[0022] (1)通过稀疏表示方法对单信道语音去混响进行建模分析，在忽略环境噪声的情况下，利用下式简单描述语音去混响问题：
[0023] Y＝HX
[0024] 其中Y是混响语音，H是RIR，且有X＝AS；
[0025] (2)根据X＝AS，假定纯净语音含有正交基稀疏矩阵表示矩阵A，则单信道语音去混响的稀疏表示为：
[0026]
[0027] 其中λ＞0是用于控制稀疏性与测量空间的参数；
[0028] (3)使用ISTA算法进行求解；
[0029] 3、引入深度语音字典，将字典学习问题转化为以下两种表示方法，任选其一：
[0030] 方法1，将求解目标函数定义为：
[0031] X＝DZ
[0032] 其中X是纯净语音，D是语音字典，Z是稀疏表示矩阵；
[0033] 字典学习问题表示为求解目标函数：
[0034]
[0035] 方法2，采用K-SVD字典学习算法，将求解目标函数定义为：
[0036]
[0037] 其中zi是矩阵Z中的第i列，||.||0表示zi中非零元素的个数，τ表示稀疏程度，用于控制Z的每一列中所含非零元素的最大个数；
[0038] 4、引入深度语音字典后，进行求解，具体描述为：
[0039] (1)根据单层字典学习进行求解，在求解过程中利用稀疏编码和字典更新，具体描述如下：
[0040] 1)稀疏编码，运用LARC方法，已知字典D(t-1)，在第t次字典学习迭代中稀疏矩阵的每个列向量更新如下表示：
[0041]
[0042] 其中X(:,n)表示给定的观测矩阵，D(t-1)表示字典矩阵， 表示满足稀疏约束条件的稀疏表示矩阵；
[0043] 2)字典更新，根据所获得稀疏表示矩阵Z，从训练信号集X中寻找当次迭代中更适合训练信号的字典D，字典更新如下式表示：
[0044]
[0045]
[0046] 3)针对每一个原子||D(:,l)||2＝1, 残差表示为：
[0047]
[0048] 式中将DZ分解为L个秩为1的矩阵(d(:,l)z(l:,))之和，Rl代表除第l个原子之外的原子所带来的残差。通过SVD分解找到与Rl最为接近的秩为1的矩阵；
[0049] 4)进行预处理，定义Λl为下式：
[0050] Λl＝{i|1≤i≤K,c(l,:)(i)≠0}
[0051] 5)定义M为D×|Λl|的矩阵，其中(Λl(i),i)设置为1，其他设置为0，则：
[0052]
[0053] 上式 代表R(l)保留z(l,:)中非零元素所对应的原子所构成的矩阵，代表z(l,:)只保留非零元素所构成的向量，字典更新如下式表示：
[0054]
[0055] 6)对 进行SVD分解得到下式：
[0056]
[0057] 7)使用U中的第一列来跟新d(:,l)，V的第一列乘以Δ(l,l)来更新相应的系数，完成一个原子的更新；
[0058] (2)根据深度字典学习进行求解，具体过程如下：
[0059] 1)由下式描述多层字典学习，
[0060] X＝D1*D2*Z
[0061] 其中D2，D2分别是双层神经网络中每层的字典；
[0062] 2)采用逐层训练进行逐层分解，首先训练除第一特征层Z1，字典D1：
[0063] X＝D1*Z1
[0064] 其次对Z1进行分解，求解第二层字典D2、特征Z2：
[0065] Z1＝D2*Z2
[0066] 之后以此类推，得到更深层的语音字典学习；
[0067] 3)基于贪婪深度字典学习的单信道语音去混响模型由此建立。
[0068] 本发明将传统的单信道语音去混响方法与深度字典学习相结合，提高单信道去混响的效果，使其具有良好的去混响效果，提高深度网络结构的可推导性。
附图说明
[0069] 图1本发明的方法流程图。
[0070] 图2单层字典学习的示意图。
[0071] 图3双层学习示意图。
[0072] 图4贪婪的分层学习示意图。
具体实施方式
[0073] 下面举例对本发明做更详细的描述。
[0074] 如图1所示，本发明基于贪婪深度字典学习的单信道语音去混响技术，首先将输入的混沌、纯净语音进行信号预处理，进行特征的提取；其次，对单信道语音去混响问题结合
贪婪深度字典学习进行模型的建立，具体描述如下：
[0075] 1、对单信道语音去混响问题进行数学建模。
[0076] 通过混响产生的具体的物理原理，对单信道语音去混响问题进行尽可能的数学建模。混响是指在声学上，除延迟时间达到约50ms以上的反应波称为回声外，其余的反射波产
生的效应统称为混响。早期的混响会产生一种声染色现象，它会导致语音信号在频域失真。
后期混响的存在会导致语音模糊不清，降低语音可懂度。因此，提出该发明能够有效的解决
上述问题。具体描述方法如下所示。
[0077] (1)将n时刻的语音信号利用下式进行描述：
[0078]
[0079] 其中D是原始语音，nk是k次反射后的延时单元数，ρk是第k次反射时的发射系数，s(n-knk)是第k次反射后的语音信号，u(n)表示周围环境的噪声，x(n)表示混响语音。
[0080] (2)为了对混响问题进行建模分析，将步骤1中(1)过程的环境噪声u(n)去掉简化得：
[0081]
[0082] (3)根据单位冲激函数de卷积性质，将步骤1中(2)的公式改写为卷积的形式：
[0083]
[0084] 其中，h(n)是冲击响应(Room Impulse Response,RIR),*是卷积操作符号，δ(n)表示单位冲激函数，h(n)可以表示混响的物理特征，它是位于knk并且都衰减ρk倍后的许多冲
击函数δ(n)的叠加。
[0085] h(n)是冲击响应(Room Impulse Response,RIR)也可称为房间脉冲响应。房间脉冲响应是一个由生源位置发出到麦克风位置的脉冲响应撩号，它反映了声学环境所有声学
特性。由步骤1中(3)式可知，房间脉冲响应h(n)和纯净语音s(n)卷积后获得混响语音x(n)。
因此，去混响的过程就是在混响语音x(n)中去除房间的脉冲响应h(n)的影响，最终获得纯
净的语音s(n),因此这一过程的关键点就是求出h(n)的逆。
[0086] 2、将单信道语音去混响模型转化为基于稀疏表示的单信道语音去混响模型，转化过程如下：
[0087] (1)通过稀疏表示方法对单信道语音去混响问题进行建模分析。在忽略环境噪声的情况下，利用下式描述语音去混响问题：
[0088] Y＝HX
[0089] 其中Y是混响语音，H是RIR，且有X＝AS。
[0090] (2)根据步骤2中(1)的X＝AS，假定纯净语音含有正交基稀疏矩阵表示矩阵A，则单信道语音去混响的稀疏表示可以写为：
[0091]
[0092] 其中λ＞0是用于控制稀疏性与测量空间的参数。
[0093] (3)步骤2中的(2)式，使用ISTA算法进行求解。
[0094] 3、在步骤2的求解基础上，将深度语音字典引入混响的问题中。为了获得深度语音字典，将字典学习问题可以转化为以下两种表示方法，任选其一：
[0095] 方法1，将求解目标函数定义为：
[0096] X＝DZ
[0097] 其中X是纯净语音，D是语音字典，Z是稀疏表示矩阵。从矩阵分解角度看字典学习过程：给定样本数据集X，X的每一列表示一个样本。如图2所示，字典的学习目标就是把X矩
阵分解为D、Z矩阵。
[0098] 为了获得字典，字典学习问题通常可以表示为求解目标函数：
[0099]
[0100] 方法2，根据K-SVD字典学习算法，求解目标函数定义为：
[0101]
[0102] 其中zi是矩阵Z中的第i列，||·||0表示zi中非零元素的个数，τ表示稀疏程度，用于控制Z的每一列中所含非零元素的最大个数。
[0103] K-SVD是一种经典的字典训练算法，依据误差最小原则，对误差进行奇异值分解(singular value decomposition,SVD)分解，选择使误差最小的分解项作为更新的字典原
子和对应的原子字数，经过不断的迭代从而得到优化的解。
[0104] 4、引入深度语音字典后，进行问题的求解。具体描述为：
[0105] (1)在求解过程中利用稀疏编码和字典更新进行单层字典学习的求解。稀疏编码是指给定观测矩阵X和字典矩阵Z，估计满足稀疏约束条件的稀疏表示矩阵Z。利用LARC的方
法，具体描述如下：
[0106] 1)稀疏编码，已知字典D(t-1)，在第t次字典学习迭代中稀疏矩阵的每个列向量更新如下表示：
[0107]
[0108] 其中X(:,n)表示给定的观测矩阵，D(t-1)表示字典矩阵， 表示满足稀疏约束条件的稀疏表示矩阵。
[0109] 2)字典更新，根据步骤4(1)中1)所获得稀疏表示矩阵Z，从训练信号集X中寻找档次迭代中更适合训练信号的字典D，字典更新如下式表示：
[0110]
[0111]
[0112] 该问题可以通过K-SVD字典更新来解决，K-SVD字典更新是对字典中的每一个原子逐个更新。
[0113] 3)针对每一个原子||D(:,l)||2＝1, 残差可以表示为：
[0114]
[0115] 式中将DZ分解为L个秩为1的矩阵(d(:,l)z(l,:))之和。Rl代表除第l个原子之外的原子所带来的残差。通过SVD分解找到与Rl最为接近的秩为1的矩阵
[0116] 4)为了保证编码结果的有效性，在SVD分解前，对步骤4(1)中3)式进行预处理，定义Λl为下式：
[0117] Λl＝{i|1≤i≤K,c(l,:)(i)≠0}
[0118] 5)定义M为D×|Λl|的矩阵，其中(Λl(i),i)设置为1，其他设置为0，则步骤4(1)中2)式等价于下式：
[0119]
[0120] 上式 代表式3)中R(l)保留z(l,:)中非零元素所对应的原子所构成的矩阵。同样 代表z(l,:)只保留非零元素所构成的向量，字典更新如下式表示：
[0121]
[0122] 6)对上式 进行SVD分解得到下式：
[0123]
[0124] 7)使用U中的第一列来更新d(:,l)，V的第一列乘以Δ(l,l)来更新相应的系数，至此完成了一个原子的更新。
[0125] (2)根据深度字典学习进行求解过程，具体过程如下：
[0126] 1)由下式描述多层字典学习，其网络结构如图3所示。
[0127] X＝D1*D2*Z
[0128] 其中D2，D2分别是双层神经网络中每层的字典。
[0129] 2)采用逐层训练的思想对步骤4(2)中1)进行逐层分解。以双层分解为例进行分解如图4所示，首先训练除第一特征层Z1，字典D1：
[0130] X＝D1*Z1
[0131] 其次对Z1进行分解，求解第二层字典D2、特征Z2：
[0132] Z1＝D2*Z2
[0133] 之后以此类推，得到更深层的字典学习。
[0134] 3)先进的基于贪婪深度字典学习的单信道语音去混响模型由此建立。
[0135] 最后，将建立好的模型与提取的特征进行模型训练，得到基本系数后，再次输入混沌、纯净语音信号数据进行测试，然后经过傅里叶逆变换，最终得到纯净语音。
[0136] 一种基于贪婪深度字典学习的单信道语音去混响技术将传统的单信道语音去混响方法与深度语音字典学习相结合，能够提高模型的可解释性和可推导性，有效的提高语
音去混响的效果。
