一种语音信号处理方法及装置，所述方法包括：对输入的语音信号做T采样的延迟，得到延迟信号；对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号；将所述输入的语音信号与所述滤波信号时域相减，得到消除信号；对所述消除信号进行维纳滤波，得到处理后的语音信号，并将所述维纳滤波的增益以及所述消除信号作为控制信号控制所述加窗频域自适应滤波的滤波系数更新。所述方法及装置可以在消除鸣号声干扰的同时，保证语音信号质量。
对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号；
将所述输入的语音信号与所述滤波信号时域相减，得到消除信号；
对所述消除信号进行维纳滤波，得到处理后的语音信号，并将所述维纳滤波的增益以及所述消除信号作为控制信号控制所述加窗频域自适应滤波的滤波系数更新；
所述将所述维纳滤波的增益以及所述消除信号作为控制信号控制所述加窗频域自适应滤波的滤波系数更新，通过如下公式进行：其中：Ht为所述加窗频域自适应滤波
器的频域系数，Ht为长度为N的一维向量且Ht＝F[ht 0N-M]；F为离散傅里叶矩阵；ht为所述加窗频域自适应滤波器的时域系数，ht为长度为M的一维向量且ht＝[ht(1),ht(2),…,ht(M-
1),ht(M)]T，M满足：M>fs*5ms；F-为逆离散傅里叶矩阵，RM为约束矩阵，M个上对角元素为1，其余元素为0；λ为非线性控制因子且λ＞0；μu为更新步长且0＜μu＜2；μL为泄露系数且0＜μL＜
1；ε为防止除零错误系数且ε＞0；上标*为共轭操作，Zw,t为z(t)的N点加窗傅里叶变换，Ew,t为e(t)的N点加窗傅里叶变换，E[|Zw,t|2]是|Zw,t|2在t时刻的短期数学期望值，e(t)为所述消除信号，Gt为维纳滤波器增益。
2.如权利要求1所述的语音信号处理方法，其特征在于，所述将所述输入的语音信号与所述滤波信号时域相减，得到消除信号，通过如下公式进行：e(t)＝x(t)-y(t)；
其中：x(t)为所述输入的语音信号，e(t)为所述消除信号，y(t)为所述滤波信号。
3.如权利要求2所述的语音信号处理方法，其特征在于，所述对所述延迟信号进行加窗频域自适应滤波处理包括以下任一种：在时域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理；
在频域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理。
4.如权利要求3所述的语音信号处理方法，其特征在于，所述在时域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理，通过如下公式进行：其中，z(t)为所述延迟信号， 为卷积操作符号。
5.如权利要求3所述的语音信号处理方法，其特征在于，所述在频域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理，通过如下公式进行：[0M,y(t-N+M+1),…,y(t-1),y(t)]T＝RN-M{F-[Zt·Ht]}；
其中：·为点乘操作符号；RN-M为约束矩阵，N-M个下对角元素为1，其余元素为0；Zt＝F[z(t-N+1),…,z(t-1),z(t)]T，N
延迟单元，用于对输入的语音信号做T采样的延迟，得到延迟信号；
加窗频域自适应滤波器，用于对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号；
减法器，用于将所述输入的语音信号与所述滤波信号时域相减，得到消除信号；
维纳滤波器，用于对所述消除信号进行维纳滤波，得到处理后的语音信号；
反馈单元，用于将所述维纳滤波器的滤波增益以及所述消除信号反馈至所述加窗频域自适应滤波器，以控制所述加窗频域自适应滤波器的滤波系数；
所述反馈单元用于采用如下公式将所述维纳滤波器的滤波增益以及所述消除信号反馈至所述加窗频域自适应滤波器： 其中：Ht为所述加窗频域自适应滤波器的频域系数，Ht为长度为N的一维向量且Ht＝F[ht
0N-M]；F为离散傅里叶矩阵；ht为所述加窗频域自适应滤波器的时域系数，ht为长度为M的一维向量且ht＝[ht(1),ht(2),…,ht(M-1),ht(M)]T，M满足：M>fs*5ms；F-为逆离散傅里叶矩阵，RM为约束矩阵，M个上对角元素为1，其余元素为0；λ为非线性控制因子且λ＞0；μu为更新步长且0＜μu＜2；μL为泄露系数且0＜μL＜1；ε为防止除零错误系数且ε＞0；上标*为共轭操作，Zw,t为z(t)的N点加窗傅里叶变换，Ew,t为e(t)的N点加窗傅里叶变换，E[|Zw,t|2]是|Zw,t|2在t时刻的短期数学期望值，e(t)为所述消除信号，Gt为维纳滤波器增益。
7.如权利要求6所述的语音信号处理装置，其特征在于，所述减法器用于采用公式e(t)＝x(t)-y(t)将输入的语音信号与所述滤波信号时域相减，得到消除信号，其中：x(t)为所述输入的语音信号，e(t)为所述消除信号，y(t)为所述滤波信号。
8.如权利要求7所述的语音信号处理装置，其特征在于，所述加窗频域自适应滤波器用于在时域上或频域上对所述延迟信号进行加窗频域自适应滤波处理。
9.如权利要求8所述的语音信号处理装置，其特征在于，所述加窗频域自适应滤波器用于采用如下公式对所述延迟信号在时域上进行加窗频域自适应滤波处理：其中，z(t)为所述延迟信号， 为卷积操作符号。
10.如权利要求8所述的语音信号处理装置，其特征在于，所述加窗频域自适应滤波器用于采用如下公式对所述延迟信号在频域上进行加窗频域自适应滤波处理：[0M,y(t-N+M+1),…,y(t-1),y(t)]T＝RN-M{F-[Zt·Ht]}；
其中：·为点乘操作符号；RN-M为约束矩阵，N-M个下对角元素为1，其余元素为0；Zt＝F[z(t-N+1),…,z(t-1),z(t)]T，N
语音信号处理方法及装置
技术领域
[0001] 本发明涉及通信领域，尤其涉及一种语音信号处理方法及装置。
背景技术
[0002] 在日常生活中，鸣号声(Horn)较为常见，通常源自交通工具，例如汽车的鸣笛声。当交通状况较为恶劣时，鸣号声出现的几率大大增加，会严重干扰语音通话的质量。
[0003] 现有的智能电话通常配备有双麦克风，通过双麦克风能够有效消除鸣号声。然而，针对一些低端的智能电话，出于成本考虑，可能只配置一个麦克风，导致单麦克风的智能电话容易受到鸣号声的影响。
[0004] 现有技术中，为减少鸣号声的影响，根据鸣号声的可预测性，使用预测器将之前的语音信号进行处理，生成一个预测信号，用于抵消鸣号声。然而，语音信号同样具有一定的可预测性，传统的预测器会消除语音信号中的一些语音成分，导致语音信号质量较差。
发明内容
[0005] 本发明解决的技术问题是单麦克风终端，如何在消除鸣号声干扰的同时保证语音信号质量。
[0006] 为解决上述技术问题，本发明实施例提供一种语音信号处理方法，包括：对输入的语音信号做T采样的延迟，得到延迟信号；对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号；将所述输入的语音信号与所述滤波信号时域相减，得到消除信号；对所述消除信号进行维纳滤波，得到处理后的语音信号，并将所述维纳滤波的增益以及所述消除信号作为控制信号控制所述加窗频域自适应滤波的滤波系数更新。
[0007] 可选的，所述对输入的语音信号做T采样的延迟，得到延迟信号，采用如下公式进行：z(t)＝x(t-T)；其中：z(t)为所述延迟信号，x(t)为所述输入的语音信号，T为延迟时长且满足0.1ms≤T/fs≤1s，fs为x(t)的采样频率。
[0008] 可选的，所述将所述输入的语音信号与所述滤波信号时域相减，得到消除信号，通过如下公式进行：e(t)＝x(t)-y(t)；其中：e(t)为所述消除信号，y(t)为所述滤波信号。
[0009] 可选的，所述对所述延迟信号进行加窗频域自适应滤波处理包括以下任一种：在时域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理；在频域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理。
[0010] 可选的，所述在时域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理，通过如下公式进行：
[0011]
[0012] 其中，ht为所述加窗频域自适应滤波器的时域系数，ht为长度为M的一维向量且ht＝[ht(1),ht(2),…,ht(M-1),ht(M)]T， 为卷积操作符号，M满足：M>fs*5ms。
[0013] 可选的，所述在频域上采用加窗频域自适应滤波器对所述延迟信号进行加窗频域自适应滤波处理，通过如下公式进行：
[0014] [0M,y(t-N+M+1),…,y(t-1),y(t)]T＝RN-M{F-[Zt·Ht]}；
[0015] 其中：Ht为所述加窗频域自适应滤波器的频域系数；Ht为长度为N的一维向量且Ht＝F[ht 0N-M]；F为离散傅里叶矩阵；·为点乘操作符号；RN-M为约束矩阵，N-M个下对角元素为1，其余元素为0；F-为逆离散傅里叶矩阵，Zt＝F[z(t-N+1),…,z(t-1),z(t)]T，N
[0016] 可选的，所述将所述维纳滤波的增益以及所述消除信号作为控制信号控制所述加窗频域自适应滤波的滤波系数更新，通过如下公式进行：
[0017]
[0018] 其中：RM为约束矩阵，M个上对角元素为1，其余元素为0；λ为非线性控制因子且λ＞0；μu为更新步长且0＜μu＜2；μL为泄露系数且0＜μL＜1；ε为防止除零错误系数且ε＞0；上标*为共轭操作，Zw,t为z(t)的N点加窗傅里叶变换，Ew,t为e(t)的N点加窗傅里叶变换，E[|Zw,t|2]是|Zw,t|2在t时刻的短期数学期望值。
[0019] 本发明实施例还提供了一种语音信号处理装置，包括：延迟单元，用于对输入的语音信号做T采样的延迟，得到延迟信号；加窗频域自适应滤波器，用于对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号；减法器，用于将所述输入的语音信号与所述滤波信号时域相减，得到消除信号；维纳滤波器，用于对所述消除信号进行维纳滤波，得到处理后的语音信号；反馈单元，用于将所述维纳滤波器的滤波增益以及所述消除信号反馈至所述加窗频域自适应滤波器，以控制所述加窗频域自适应滤波器的滤波系数。
[0020] 可选的，所述延迟单元用于采用公式z(t)＝x(t-T)对输入的语音做T采样的延迟，得到延迟信号，其中：z(t)为所述延迟信号，x(t)为所述输入的语音信号，T为延迟时长且满足0.1ms≤T/fs≤1s，fs为x(t)的采样频率。
[0021] 可选的，所述减法器用于采用公式e(t)＝x(t)-y(t)将输入的语音信号与所述滤波信号时域相减，得到消除信号，其中：e(t)为所述消除信号，y(t)为所述滤波信号。
[0022] 可选的，所述加窗频域自适应滤波器用于在时域上或频域上对所述延迟信号进行加窗频域自适应滤波处理。
[0023] 可选的，所述加窗频域自适应滤波器用于采用如下公式对所述延迟信号在时域上进行加窗频域自适应滤波处理：
[0024] 其中，ht为所述加窗频域自适应滤波器的时域系数，ht为长度为M的一维向量且ht＝[ht(1),ht(2),…,ht(M-1),ht(M)]T， 为卷积操作符号，M满足：M>fs*5ms。
[0025] 可选的，所述加窗频域自适应滤波器用于采用如下公式对所述延迟信号在频域上进行加窗频域自适应滤波处理：
[0026] [0M,y(t-N+M+1),…,y(t-1),y(t)]T＝RN-M{F-[Zt·Ht]}；
[0027] 其中：Ht为所述加窗频域自适应滤波器的频域系数；Ht为长度为N的一维向量且Ht＝F[ht 0N-M]；F为离散傅里叶矩阵；·为点乘操作符号；RN-M为约束矩阵，N-M个下对角元素为1，其余元素为0；F-为逆离散傅里叶矩阵，Zt＝F[z(t-N+1),…,z(t-1),z(t)]T，N
[0028] 可选的，所述反馈单元用于采用如下公式将所述维纳滤波器的滤波增益以及所述消除信号反馈至所述加窗频域自适应滤波器：
[0029]
[0030] 其中：RM为约束矩阵，M个上对角元素为1，其余元素为0；λ为非线性控制因子且λ＞0；μu为更新步长且0＜μu＜2；μL为泄露系数且0＜μL＜1；ε为防止除零错误系数且ε＞0；上标*为共轭操作，Zw,t为z(t)的N点加窗傅里叶变换，Ew,t为e(t)的N点加窗傅里叶变换，E[|Zw,t|2]是|Zw,t|2在t时刻的短期数学期望值。
[0031] 与现有技术相比，本发明实施例的技术方案具有以下有益效果：
[0032] 对延迟信号进行加窗频域自适应滤波处理，得到滤波信号，将输入的语音信号与滤波信号时域相减，即可将输入的语音信号中的鸣号声消除，将维纳滤波的增益以及消除信号进行反馈，从而控制加窗频域自适应滤波的滤波器系数更新，使得整个鸣号声的消除流程可控，避免消除输入的语音信号中的语音部分。因此，在有效消除鸣号声的同时，保证语音信号质量。
附图说明
[0033] 图1是本发明实施例中的一种语音信号处理方法的流程图；
[0034] 图2(a)是存在鸣号声干扰的原始语音信号波形图；
[0035] 图2(b)是消除鸣号声干扰的语音信号波形图；
[0036] 图3是本发明实施例中的一种语音信号处理装置的结构示意图。
具体实施方式
[0037] 现有技术中，鸣号声具有较强的可预测性。为减少鸣号声的干扰，可以使用预测器将之前的语音信号进行处理，生成一个预测信号，用于抵消鸣号声。然而，在实际应用中，语音信号同样具有一定的可预测性，因此，传统的预测器在消除鸣号声的同时，会消除语音信号中的一些语音成分，导致语音信号质量较差。
[0038] 在本发明实施例中，对延迟信号进行加窗频域自适应滤波处理，得到滤波信号。将输入的语音信号与滤波信号时域相减，即可将输入的语音信号中的鸣号声消除。将维纳滤波的增益以及消除信号进行反馈，从而控制加窗频域自适应滤波的滤波器系数更新，使得整个鸣号声的消除流程是可控的，从而避免消除输入的语音信号中的语音部分。因此，在消除鸣号声的同时，保证语音信号质量。
[0039] 为使本发明的上述目的、特征和有益效果能够更为明显易懂，下面结合附图对本发明的具体实施例做详细的说明。
[0040] 本发明实施例提供了一种语音信号处理方法，参照图1，以下通过具体步骤进行详细说明。
[0041] 步骤S101，对输入的语音信号做T采样的延迟，得到延迟信号。
[0042] 在本发明实施例中，可以采用公式z(t)＝x(t-T)对输入的语音信号做T采样的延迟，其中：z(t)为经过T采样的延迟之后的延迟信号，x(t)为输入的语音信号，T为延迟时长，且T满足如下条件：0.1ms≤T/fs≤1s，其中，fs为x(t)的采样频率。
[0043] 在实际应用中，T的取值不宜过大也不宜过小。当T过大时，鸣号声的相关性会降低；当T过小时，T易落在语音信号的强相关时间内，导致语音信号存在较大的损失。在本发明一实施例中，T/fs的取值为2ms。
[0044] 在得到延迟信号后，执行步骤S102。
[0045] 步骤S102，对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号。
[0046] 在具体实施中，可以采用预设的加窗频域自适应滤波器对延迟信号z(t)进行加窗频域自适应滤波。可以在时域上对延迟信号z(t)进行加窗频域自适应滤波，也可以在频域上对延迟信号z(t)进行加窗频域自适应滤波。
[0047] 在本发明一实施例中，加窗频域自适应滤波器在时域上对延迟信号z(t)进行加窗频域自适应滤波。
[0048] 加窗频域自适应滤波器的时域系数为ht，ht为长度为M的一维向量，且ht＝[ht(1),ht(2),…,ht(M-1),ht(M)]T。
[0049] 在时域上，加窗频域自适应滤波器的滤波过程可以通过如下公式实现：
[0050]
[0051] 其中， 为卷积操作符号，M的取值满足如下关系：M>fs*5ms，也即5ms＜M/fs。在本发明一实施例中，M/fs＝20ms。
[0052] 在本发明另一实施例中，加窗频域自适应滤波器在频域上对延迟信号z(t)进行加窗频域自适应滤波。加窗频域自适应滤波器的频域系数为Ht，Ht为长度为N的一维向量，通过如下方式生成：在ht的尾部补N-M个0，得到长度为N的一维向量，对得到的向量进行快速傅里叶变换，得到的结果即为Ht，也即：Ht＝F[ht 0N-M]，F为离散傅里叶矩阵。
[0053] 在频域上，加窗频域自适应滤波器的滤波过程可以通过如下公式实现：
[0054] [0M,y(t-N+M+1),…,y(t-1),y(t)]T＝RN-M{F-[Zt·Ht]}； (2)
[0055] 其中，i为点乘操作符号；RN-M为约束矩阵，只有N-M个下对角元素为1，其余元素为0；F-为逆离散傅里叶矩阵，Zt＝F[z(t-N+1),…,z(t-1),z(t)]T；N的取值满足如下关系：N
[0056] 在本发明一实施例中，N/fs＜100＝32ms。
[0057] 在获取到滤波信号之后，可以执行步骤S103。
[0058] 步骤S103，将所述输入的语音信号与所述滤波信号时域相减，得到消除信号。
[0059] 在具体实施中，将输入的语音信号x(t)与滤波信号y(t)时域相减之后，可以将输入的语音信号x(t)中的鸣号声消除，也即得到的消除信号中，不存在鸣号声的干扰。
[0060] 在本发明实施例中，可以采用如下公式获取消除信号：e(t)＝x(t)-y(t)，其中：e(t)为消除信号。
[0061] 在得到消除信号e(t)后，可以执行步骤S104。
[0062] 步骤S104，对所述消除信号进行维纳滤波，得到处理后的语音信号，并将所述维纳滤波的增益以及所述消除信号作为控制信号控制所述加窗频域自适应滤波的滤波系数更新。
[0063] 在具体实施中，可以采用预设的维纳滤波器对消除信号e(t)进行维纳滤波。
[0064] 在实际应用中，维纳滤波器是一种经典的最优滤波器，通常的实现形式是通过估计噪声信号(信号中的无用成分)来求取维纳滤波器增益G[k]t。维纳滤波器增益G[k]t应用到被滤波信号后，被滤波信号中的有用成分保留，无用成分会被衰减乃至消除，即：Dt＝Et·Gt，其中Dt为维纳滤波器的输出信号，也即处理后的频域上的语音信号，Et为被滤波信号的短时傅里叶变换。
[0065] 在本发明实施例中，被滤波信号为消除信号e(t)，Et可以采用如下公式计算：
[0066]
[0067] 其中，wana(k)是N点的分析窗，k＝1,2,3,……,N；分析窗一般用正弦窗或者Hanning窗；Ps是消除信号e(t)的自功率谱。
[0068] 在实际应用中常用的一种频域维纳滤波器的迭代计算形式如下，t时刻的维纳滤波器增益为：
[0069]
[0070] 信噪比SNR为：
[0071] SNR[k]t＝αG[k]t-12SNRP[k]t-1+(1-α)max(SNRP[k]t-1,0)； (5)
[0072] SNRp为先验信噪比，采用如下公式计算：
[0073]
[0074] 其中，α取值范围为0<α<1，在本发明一实施例中α＝0.5。
[0075] 可以采用跟踪信号频谱短时最小值方法计算噪声自功率谱Pn，也可以采用时间递归平均方法计算噪声自功率谱Pn，还可以采用其他的方法计算噪声自功率谱Pn，此处不做赘述。
[0076] 本发明一实施例中，采用G.Doblinger噪声估计法计算噪声自功率谱Ps，G.Doblinger噪声估计法是频谱短时最小值方法与时间递归平均方法的结合，具体计算公式为：
[0077]
[0078] 其中η1，η2，η3是平滑因子，它们的取值范围为0<η1，η2，η3<1，典型值η1＝0.99，η2＝0.99，η3＝0.8。
[0079] 在本发明实施例中，在获取到维纳滤波器的增益后，可以将维纳滤波器的滤波增益以及消除信号e(t)作为控制信号，反馈至加窗频域自适应滤波器，以对加窗频域自适应滤波器的滤波系数进行更新。
[0080] 在本发明一实施例中，采用如下公式对加窗频域自适应滤波器的滤波系数进行更新：
[0081]
[0082] 公式(8)中：RM为约束矩阵，M个上对角元素为1，其余元素为0。λ为大于0的非线性控制因子，在本发明实施例中，λ＝1。μu为更新步长，取值范围过大时稳定性较差，取值过小则更新时长较长，在本发明实施例中，μu的取值范围为0＜μu＜2。μL为泄露系数，且μL的取值范围为0＜μL＜1，在本发明实施例中，μL比1略小，用于消除远期信号的影响。ε为防止除零错误系数的微小正数；上标*为共轭操作，Zw,t为z(t)的N点加窗傅里叶变换，Ew,t为e(t)的N点加窗傅里叶变换，E[|Zw,t|2]是|Zw,t|2在t时刻的短期数学期望值。
[0083] 在本发明实施例中，Zw,t与Ew,t分别采用如下公式计算：
[0084]
[0085] 公式(9)中，上标T为转置操作，wadf(k)是N点的用于自适应滤波器的窗函数，k＝1,2,……,N，窗函数可以为矩形窗、Sine窗、Hanning窗、Hamming窗以及Tukey窗等中的任一种，也可以为其他类型的窗函数。在本发明一实施例中，wadf(k)采用Tukey窗。
[0086] 在本发明一实施例中，采用自回归方法计算获取E[|Zw,t|2]，计算公式如下：
[0087] E[|Zw,t|2]＝γE[|Zw,t-1|2]+(1-γ)|Zw,t|2,0<γ<1； (10)
[0088] 在本发明其他实施例中，还可以存在其他的方法来计算获取E[|Zw,t|2]，此处不做赘述。
[0089] 参照图2，图2(a)给出了原始语音信号波形图，图2(b)给出了采用本发明实施例中提供的语音信号处理方法处理后的语音信号波形图。
[0090] 参照图2(a)，受鸣号声干扰的语音信号波形段201中的黑色斑点状2021为鸣号声干扰。从201中可以获知，当前的语音信号存在较大的鸣号声干扰。
[0091] 而在图2(b)中，201中的黑色斑点状已经被完全滤除，也即当前的语音信号不存在鸣号声干扰。此外，与图2(a)相比，在消除鸣号声的同时，语音信号的质量并没有受到影响。
[0092] 由此可见，对延迟信号进行加窗频域自适应滤波处理，得到滤波信号。将输入的语音信号与滤波信号时域相减，即可将输入的语音信号中的鸣号声消除。将维纳滤波的增益以及消除信号进行反馈，从而控制加窗频域自适应滤波的滤波器系数更新，使得整个鸣号声的消除流程是可控的，从而避免消除输入的语音信号中的语音部分。因此，在消除鸣号声的同时，保证语音信号质量。
[0093] 参照图3，给出了本发明实施例中的一种语音信号处理装置，包括：延迟单元301、加窗频域自适应滤波器302、减法器303、维纳滤波器304以及反馈单元305，其中：
[0094] 延迟单元301，用于对输入的语音信号做T采样的延迟，得到延迟信号；
[0095] 加窗频域自适应滤波器302，用于对所述延迟信号进行加窗频域自适应滤波处理，得到滤波信号；
[0096] 减法器303，用于将所述输入的语音信号与所述滤波信号时域相减，得到消除信号；
[0097] 维纳滤波器304，用于对所述消除信号进行维纳滤波，得到处理后的语音信号Dt；
[0098] 反馈单元305，用于将所述维纳滤波器的滤波增益以及所述消除信号反馈至所述加窗频域自适应滤波器，以控制所述加窗频域自适应滤波器的滤波系数。
[0099] 在具体实施中，所述延迟单元301可以用于采用公式z(t)＝x(t-T)对输入的语音做T采样的延迟，得到延迟信号，其中：z(t)为所述延迟信号，x(t)为所述输入的语音信号，T为延迟时长且满足0.1ms≤T/fs≤1s，fs为x(t)的采样频率。
[0100] 在具体实施中，所述减法器303可以用于采用公式e(t)＝x(t)-y(t)将输入的语音信号与所述滤波信号时域相减，得到消除信号，其中：e(t)为所述消除信号，y(t)为所述滤波信号。
[0101] 在具体实施中，所述加窗频域自适应滤波器302可以用于在时域上对所述延迟信号进行加窗频域自适应滤波处理，也可以用于在频域上对所述延迟信号进行加窗频域自适应滤波处理。
[0102] 在具体实施中，所述加窗频域自适应滤波器302可以用于采用如下公式对所述延迟信号在时域上进行加窗频域自适应滤波处理：
[0103]
[0104] 其中，ht为所述加窗频域自适应滤波器的时域系数，ht为长度为M的一维向量且ht＝[ht(1),ht(2),…,ht(M-1),ht(M)]T， 为卷积操作符号，M满足：M>fs*5ms。
[0105] 在具体实施中，所述加窗频域自适应滤波器用于采用如下公式对所述延迟信号在频域上进行加窗频域自适应滤波处理：
[0106] [0M,y(t-N+M+1),…,y(t-1),y(t)]T＝RN-M{F-[Zt·Ht]}；
[0107] 其中：Ht为所述加窗频域自适应滤波器的频域系数；Ht为长度为N的一维向量且Ht＝F[ht 0N-M]；F为离散傅里叶矩阵；·为点乘操作符号；RN-M为约束矩阵，N-M个下对角元素为1，其余元素为0；F-为逆离散傅里叶矩阵，Zt＝F[z(t-N+1),…,z(t-1),z(t)]T，N
[0108] 在具体实施中，所述反馈单元305可以用于采用如下公式将所述维纳滤波器的滤波增益以及所述消除信号反馈至所述加窗频域自适应滤波器：
[0109]
[0110] 其中：RM为约束矩阵，M个上对角元素为1，其余元素为0；λ为非线性控制因子且λ＞0；μu为更新步长且0＜μu＜2；μL为泄露系数且0＜μL＜1；ε为防止除零错误系数且ε＞0；上标*为共轭操作，Zw,t为z(t)的N点加窗傅里叶变换，Ew,t为e(t)的N点加窗傅里叶变换，E「|Zw,t|2]是|Zw,t|2在t时刻的短期数学期望值。
[0111] 本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：ROM、RAM、磁盘或光盘等。
[0112] 虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。
