一种声学回声消除器启动阶段稳健的步长调整方法
一种声学回声消除器启动阶段稳健的步长调整方法
【技术领域】
[0001] 本发明设及一种声学回声消除器。特别是设及一种声学回声消除器启动阶段稳健 的步长调整方法。
【背景技术】
[0002] 随着通信技术的不断发展，人们对便捷交流方式的要求越来越高，各种手持电信 设备、视频会议和VoIP软件电话应用越来越广。在运类通信终端中，近端语音由扬声器传到 远端，由于远端传声器和扬声器之间的禪合，使得近端语音传回本地，形成声学回声。声学 回声严重影响语音传输质量，因此声学回声抵消系统是不可或缺的。声学回声消除器 (Acoustic Echo Cancellation,AEC)是解决声学回声消除的最佳解决方案之一，一般的声 学回声消除器至少包含双端通话检测模块化ouble Talk Detector，DTD)、线性回声消除两 个部分。图1是一个典型的声学回声消除系统。
[0003] 麦克风在时刻η接收的信号为：
[0004] d(n) =y(n)+v(n)+w(n) 1.1
[0005] 式中：d(n)代表近端信号，y(n)代表远端信号x(n)经过扬声器播放后产生的回声， v(n)代表近端语音信号，w(n)代表近端噪声信号。远端信号x(n)经过系统传递函数h滤波后 形成y(n)
[0010] X是远端信号向量，N为房间冲击响应长度，T代表矩阵的转置。
[0011]回声消除的目标就是设计一个自适应有限冲击响应滤波器估计出介于麦克风与 扬声器之间的回声路径&(n)，然后根据运一估计路径得到回声估计信号r(?)，将其从d(n) 中消除，而保留v(n)。
[0014] e(n)代表线性回声滤波器消除后得到的误差信号，其中，
[0015]
1,5
[0016] L为自适应滤波器长度，实际中一般L<N。一般认为，声学回声消除器处理的场景分 为Ξ种情况:远端情况，只存在回声信号而不存在近端语音信号;近端情况，不存在回声，只 存在近端语音信号;双端通话情况，回声信号与近端语音信号同时存在。
[0017] 欠定模型双端通话鲁棒的变步长归一化最小均方误差算法(DOUBLE-TALK R0脚ST Variable Step Size Normalized Least-mean-square For Under-modeling,VSS-NLMS-UMDT)是一种新型较实用的双端通话鲁棒的回声消除算法，与其它归一化最小均方误差算 法(VSS-NLMS)算法相比不需要双端通话检测器(DTD),就能够在欠定和双端通话情况下稳 定工作，对近端信号干扰不敏感，依然保持较小且稳定的稳态失调，并且在实际应用中易于 实施与控制，不需要声学环境的任何参数，鲁棒性很强。其提出的控制步长和滤波器更新方 程是
[0020]式中，μ(η)为自适应滤波器的步长，丫 ed是e(n)与d(n)之间的互相关估计，ri(n)是 滤波器的收敛性统计参数，巧、if和晏分别代表d(n)、K")和e(n)的能量期望估计，δ、ξ 均为一常数。W上参数可由式1.8和1.9得到
[002引 Ε{ ·}表示数学期望，λ是一极小正常数，标记《表示序列ρ(η)的能量估计，它可W 通过指数递归公式来计算，计算方法是
[0024]
1.10
[0025] 此算法虽然优势明显，但仍存在一些弊端。此算法在系统启动阶段，需要用较大的 步长
来更新滤波器，使其快速收敛到某一稳定状态。若系统启动阶 段处于双端通话情况，此时滤波器会将近端语音过滤掉，运是不被系统允许的，大大降低系 统性能;若系统启动阶段采用小控制步长
调整的话，滤波器收敛速度将会变慢。由于此算法不需要DTD模块，所W系统自身无法判定 是处于近端抑或是远端情况，所W不能选择合适得步长初始化滤波器。
【发明内容】
[0026] 本发明所要解决的技术问题是，提供一种在双端通话情况不仅能保证滤波器快速 收敛到某一稳定状态，还能在启动阶段完整保留近端语音，提高了系统性能，实用性更强的 声学回声消除器启动阶段稳健的步长调整方法。
[0027] 本发明所采用的技术方案是：一种声学回声消除器启动阶段稳健的步长调整方 法，包括如下步骤：
[002引1)确定算法参数；
[0029] 2)先验滤波，所述的先验滤波，包括:去除远端信号向量X中距离η时刻最远的数据 x(n-L)，将当前远端信号χ(η)更新至远端信号向量X;利用滤波器的前一状态估计当 前回声信号，并将估计出的当前回声信号从当前近端信号d(n)中滤除，得到先验误差信号ε (η)；
[0030] 3)确定滤波器初始化步长；
[0031] 4)进行回声消除，根据步骤3)所选取的步长，利用欠定模型双端通话鲁棒的变步 长归一化最小均方误差算法进行声学回声消除。
[0032] 步骤1)所述的算法参数包括有:包括语音采样频率fs，滤波器长度L，滤波器状态 h^; <，滤波器步长μ (η)，滤波器步长极大值ymax，系统启动时间ini tTime，远端信号向量X， 先验误差信号ε(η)，后验误差信号e(n)，近端信号的能量期望估计巧、远端信号的能量期望 估计驾、估计回声信号的能量期望估计^和误差信号的能量期望估计彩，先验误差信号与 近端信号之间的互相关估计丫 ed，收敛性统计参数η(η)，收敛性统计参数期望值e邱ri(n)，收 敛性统计参数期望阔值thres。
[0033] 步骤3)所述的确定滤波器初始化步长包括如下步骤：
[0034] (1)计算出先验误差信号与近端信号之间的互相关估计丫 ed，近端信号的能量期望 估计巧，估计回声信号的能量期望估计驾:和误差信号的能量期望估计冬
[0035] (2)将步骤（1)中的参数，代入收敛性统计参数ri(n)计算公式，得到收敛性统计参 数η(η)，
[0036]
1.9
[0037] (3)定义收敛性统计参数期望值expri(n)，计算公式是
[0038] θχρη(η)=λ*θχρη(η-1)+(1-λ)*η(η) 1.11
[0039] 在系统启动时间initTime内，若e邱η(η)小于收敛性统计参数期望阔值thres，则 认为系统当前处于远端情况，采取滤波器步长极大值ymax更新滤波器;否则对步长进行更 新。
[0040] 步骤(3)所述的对步长进行更新是采用如下公式：
[0041]
1.6
[0042] 式中，δ、ξ均为一常数。
[0043] 本发明的一种声学回声消除器启动阶段稳健的步长调整方法，不需要DTD模块，在 滤波器收敛速度和系统性能之间取得折中。通过收敛性统计参数期望值在声学回声消除器 启动阶段选择合适的步长调整方式，实现在远端情况下大部分时间里选择大步长更新滤波 器，令其快速收敛，而在其它情况下保持小步长更新，既能保留近端语音，同时也能维持滤 波器快速更新。本发明的方法在双端通话情况不仅能保证滤波器快速收敛到某一稳定状 态，还能在启动阶段完整保留近端语音，提高了系统性能，实用性更强。本发明解决了 VSS-NLMS-UMDT算法在启动阶段的系统容易失调的弊端，在滤波器收敛速度和系统性能之间取 得折中，增强了系统的稳定性与实用性。
【附图说明】
[0044] 图1是一个典型的声学回声消除系统结构；
[0045] 图2a是远端语音信号时域图；
[0046] 图化是近端信号时域图；
[0047] 图2c是近端语音信号时域图；
[004引图3a是VSS-NLMS处理结果；
[0049]图 3b 是 VSS-NLMS-UMDT 处理结果；
[(K)加]图3c是VSS-NLMS-UMDT-C沈处理结果；
[0051 ]图4是本发明实施例的算法流程图。
【具体实施方式】
[0052] 下面结合实施例和附图对本发明的一种声学回声消除器启动阶段稳健的步长调 整方法做出详细说明。
[0053] 本发明的一种声学回声消除器启 动阶段稳健的步长调整方法，所要解决的关键是 声学回声消除器启动阶段如何选择合适的步长调整方法。启动阶段若在远端情况下，选取 大步长更新滤波器，使滤波器快速收敛;否则选取小步长，既能保留近端语音信号同时也会 维持滤波器的更新。
[0054] 本发明的一种声学回声消除器启动阶段稳健的步长调整方法，包括如下步骤：
[0055] 1)确定算法参数；
[0化6]所述的算法参数包括有:包括语音采样频率fs，滤波器长度L，滤波器状态斬的，滤 波器步长μ(η)，滤波器步长极大值ymax，系统启动时间initTime，远端信号向量X，先验误差 信号ε(η)，后验误差信号e(n)，近端信号的能量期望估计驾、远端信号的能量期望估计癸、 估计回声信号的能量期望估计巧和误差信号的能量期望估计勞，先验误差信号与近端信 号之间的互相关估计丫 ed，收敛性统计参数η(η)，收敛性统计参数期望值expri(n)，收敛性统 计参数期望阔值thres。
[0057] 2)先验滤波，所述的先验滤波，包括:去除远端信号向量X中距离η时刻最远的数据 x(n-L)，将当前远端信号χ(η)更新至远端信号向量X;利用滤波器的前一状态?ι(/;-1)估计当 前回声信号，并将估计出的当前回声信号从当前近端信号d(n)中滤除，得到先验误差信号ε (η)；
[005引3)确定滤波器初始化步长；
[0059] 所述的确定滤波器初始化步长包括如下步骤：
[0060] (1)计算出先验误差信号与近端信号之间的互相关估计丫 ed，近端信号的能量期望 估计遊，估计回声信号的能量期望估计和误差信号的能量期望估计勞；
[0061] (2)将步骤（1)中的参数，代入收敛性统计参数ri(n)计算公式，得到收敛性统计参 数η(η)，
[0062]
IJ
[0063] (3)收敛性统计参数ri(n)仅能够代表滤波器的当前状态，跳跃性较大，不能很好地 跟踪描述滤波器。为了更好地统计滤波器的收敛性，因此定义收敛性统计参数期望值expri (η)，计算公式是
[0064] θχρη(η)=λ*θχρη(η-1)+(1-λ)*η(η) 1.11
[0065] 实验研究发现，若滤波器在远端情况下稳定更新，则收敛性统计参数期望e邱ri(n) 大部分时间将处在某一较小值thresW下。因此在系统启动时间initTime内，若e邱ri(n)小 于收敛性统计参数期望阔值thres，则认为系统当前处于远端情况，采取滤波器步长极大值 ymax更新滤波器;否则对步长进行更新。所述的对步长进行更新是采用如下公式：
[0066]
1.6
[0067] 式中，δ、ξ均为一常数。
[0068] 4)进行回声消除，根据步骤3)所选取的步长，利用欠定模型双端通话鲁棒的变步 长归一化最小均方误差算法进行声学回声消除。
[0069] 下面，W语音采样率为16Κ、滤波器长度为1000阶的系统为例实施本专利所提方 案，处理步骤参考图4算法流程图。
[0070] 算法参数设置：语音采样频率f S = 16Κ，滤波器长度L= 1000，滤波器初始状态 h(内)=〇，滤波器初始步长μ ( η ) = 0，常数泣=30*式2，滤波器步长极大值
和滤波器初始时间initTime = 6s,远端信号向量χ = 0,先验误差信 号ε(η)=0,后验误差信号e(n)=0,近端信号、估计回声信号和误差信号的能量期望估计 巧:=巧:=安=0,先验误差信号与近端信号之间的互相关估计丫 ed = 0,收敛性统计参数η(η) =0，收敛性统计参数期望值expri(n) = 0，收敛性统计参数期望阔值thres = 0.05。
[0071 ]本算法具体实施步骤如下：
[0072] 1、读取当前远端信号x(n)和当前近端信号d(n)，将x(n)更新至远端信号向量XL (η)中；
[0073] 2、通过公式1.3和1.4,计算先验误差信号ε(η);
[0074] 3、通过公式l.s-l.io,计算出特征值丫 ed,巧，聲，繁,η(η)和ex阿(n);
[0075] 4、若当前时刻n<initTime*L，则滤波器处在启动阶段，此时若收敛性统计参数期 望值expn ( η )小于收敛性统计参数期望阔值thr e S，则令μ ( η ) = ymax ;否则，
[0076] 1、若11〉= ;[]1；[1:111116礼，贝11系统已启动结束，令
[0077]
[007引 5、将得到的步长代入式1.6得到滤波器状态}!(/?)，代入式1.4和1.3后得到后验误 差信号e(n)。若读取语音已全部处理完毕，算法终止;否则跳转回第1步。
[0079]后验误差信号e(n)即为系统的最终输出。
【主权项】
1. 一种声学回声消除器启动阶段稳健的步长调整方法，其特征在于，包括如下步骤： 1) 确定算法参数； 2) 先验滤波，所述的先验滤波，包括:去除远端信号向量X中距离η时刻最远的数据x(n-L)，将当前远端信号x(n)更新至远端信号向量X;利用滤波器的前一状态估计当前回 声信号，并将估计出的当前回声信号从当前近端信号d(n)中滤除，得到先验误差信号ε(η); 3) 确定滤波器初始化步长； 4) 进行回声消除，根据步骤3)所选取的步长，利用欠定模型双端通话鲁棒的变步长归 一化最小均方误差算法进行声学回声消除。2. 根据权利要求1所述的一种声学回声消除器启动阶段稳健的步长调整方法，其特征 在于，步骤1)所述的算法参数包括有：包括语音采样频率fs，滤波器长度L，滤波器状态 ，滤波器步长Ii(n)，滤波器步长极大值Iimax，系统启动时间ini tTime，远端信号向量X， 先验误差信号ε(η)，后验误差信号e(n)，近端信号的能量期望估计H、远端信号的能量期望 估计式、估计回声信号的能量期望估计f和误差信号的能量期望估计乾，先验误差信号与 近端信号之间的互相关估计γ ed，收敛性统计参数n(n)，收敛性统计参数期望值expn(n)，收 敛性统计参数期望阈值thres。3. 根据权利要求1所述的一种声学回声消除器启动阶段稳健的步长调整方法，其特征 在于，步骤3)所述的确定滤波器初始化步长包括如下步骤： (1) 计算出先验误差信号与近端信号之间的互相关估计yed，近端信号的能量期望估计 句，估计回声信号的能量期望估计有和误差信号的能量期望估计式2; (2) 将步骤（1)中的参数，代入收敛性统计参数η( η)计算公式，得到收敛性统计参数η (η)，(3) 定义收敛性统计参数期望值expn(n)，计算公式是 θχρη(η)=λ*θχρη(η_1)+(1_λ)*η(η) 1.11 在系统启动时间initTime内，若expn(n)小于收敛性统计参数期望阈值thres，则认为 系统当前处于远端情况，采取滤波器步长极大值ymax更新滤波器;否则对步长进行更新。4. 根据权利要求3所述的一种声学回声消除器启动阶段稳健的步长调整方法，其特征 在于，步骤(3)所述的对步长进行更新是采用如下公式：式中，均为一常数。
【专利摘要】一种声学回声消除器启动阶段稳健的步长调整方法，包括：确定算法参数；先验滤波：去除远端信号向量中距离n时刻最远的数据，将当前远端信号更新至远端信号向量；利用滤波器的前一状态估计当前回声信号，并将估计出的当前回声信号从当前近端信号中滤除，得到先验误差信号；确定滤波器初始化步长；进行回声消除，根据所选取的步长，利用欠定模型双端通话鲁棒的变步长归一化最小均方误差算法进行声学回声消除。本发明解决了VSS-NLMS-UMDT算法在启动阶段的系统容易失调的弊端，在滤波器收敛速度和系统性能之间取得折中，增强了系统的稳定性与实用性。
【IPC分类】H04M9/08
【公开号】CN105491256
【申请号】CN201510915783
【发明人】张涛, 焦海泉, 唐伟, 赵鑫
【申请人】天津大学
【公开日】2016年4月13日
【申请日】2015年12月9日
【技术领域】
[0001] 本发明设及一种声学回声消除器。特别是设及一种声学回声消除器启动阶段稳健 的步长调整方法。
【背景技术】
[0002] 随着通信技术的不断发展，人们对便捷交流方式的要求越来越高，各种手持电信 设备、视频会议和VoIP软件电话应用越来越广。在运类通信终端中，近端语音由扬声器传到 远端，由于远端传声器和扬声器之间的禪合，使得近端语音传回本地，形成声学回声。声学 回声严重影响语音传输质量，因此声学回声抵消系统是不可或缺的。声学回声消除器 (Acoustic Echo Cancellation,AEC)是解决声学回声消除的最佳解决方案之一，一般的声 学回声消除器至少包含双端通话检测模块化ouble Talk Detector，DTD)、线性回声消除两 个部分。图1是一个典型的声学回声消除系统。
[0003] 麦克风在时刻η接收的信号为：
[0004] d(n) =y(n)+v(n)+w(n) 1.1
[0005] 式中：d(n)代表近端信号，y(n)代表远端信号x(n)经过扬声器播放后产生的回声， v(n)代表近端语音信号，w(n)代表近端噪声信号。远端信号x(n)经过系统传递函数h滤波后 形成y(n)
[0010] X是远端信号向量，N为房间冲击响应长度，T代表矩阵的转置。
[0011]回声消除的目标就是设计一个自适应有限冲击响应滤波器估计出介于麦克风与 扬声器之间的回声路径&(n)，然后根据运一估计路径得到回声估计信号r(?)，将其从d(n) 中消除，而保留v(n)。
[0014] e(n)代表线性回声滤波器消除后得到的误差信号，其中，
[0015]
1,5
[0016] L为自适应滤波器长度，实际中一般L<N。一般认为，声学回声消除器处理的场景分 为Ξ种情况:远端情况，只存在回声信号而不存在近端语音信号;近端情况，不存在回声，只 存在近端语音信号;双端通话情况，回声信号与近端语音信号同时存在。
[0017] 欠定模型双端通话鲁棒的变步长归一化最小均方误差算法(DOUBLE-TALK R0脚ST Variable Step Size Normalized Least-mean-square For Under-modeling,VSS-NLMS-UMDT)是一种新型较实用的双端通话鲁棒的回声消除算法，与其它归一化最小均方误差算 法(VSS-NLMS)算法相比不需要双端通话检测器(DTD),就能够在欠定和双端通话情况下稳 定工作，对近端信号干扰不敏感，依然保持较小且稳定的稳态失调，并且在实际应用中易于 实施与控制，不需要声学环境的任何参数，鲁棒性很强。其提出的控制步长和滤波器更新方 程是
[0020]式中，μ(η)为自适应滤波器的步长，丫 ed是e(n)与d(n)之间的互相关估计，ri(n)是 滤波器的收敛性统计参数，巧、if和晏分别代表d(n)、K")和e(n)的能量期望估计，δ、ξ 均为一常数。W上参数可由式1.8和1.9得到
[002引 Ε{ ·}表示数学期望，λ是一极小正常数，标记《表示序列ρ(η)的能量估计，它可W 通过指数递归公式来计算，计算方法是
[0024]
1.10
[0025] 此算法虽然优势明显，但仍存在一些弊端。此算法在系统启动阶段，需要用较大的 步长
来更新滤波器，使其快速收敛到某一稳定状态。若系统启动阶 段处于双端通话情况，此时滤波器会将近端语音过滤掉，运是不被系统允许的，大大降低系 统性能;若系统启动阶段采用小控制步长
调整的话，滤波器收敛速度将会变慢。由于此算法不需要DTD模块，所W系统自身无法判定 是处于近端抑或是远端情况，所W不能选择合适得步长初始化滤波器。
【发明内容】
[0026] 本发明所要解决的技术问题是，提供一种在双端通话情况不仅能保证滤波器快速 收敛到某一稳定状态，还能在启动阶段完整保留近端语音，提高了系统性能，实用性更强的 声学回声消除器启动阶段稳健的步长调整方法。
[0027] 本发明所采用的技术方案是：一种声学回声消除器启动阶段稳健的步长调整方 法，包括如下步骤：
[002引1)确定算法参数；
[0029] 2)先验滤波，所述的先验滤波，包括:去除远端信号向量X中距离η时刻最远的数据 x(n-L)，将当前远端信号χ(η)更新至远端信号向量X;利用滤波器的前一状态估计当 前回声信号，并将估计出的当前回声信号从当前近端信号d(n)中滤除，得到先验误差信号ε (η)；
[0030] 3)确定滤波器初始化步长；
[0031] 4)进行回声消除，根据步骤3)所选取的步长，利用欠定模型双端通话鲁棒的变步 长归一化最小均方误差算法进行声学回声消除。
[0032] 步骤1)所述的算法参数包括有:包括语音采样频率fs，滤波器长度L，滤波器状态 h^; <，滤波器步长μ (η)，滤波器步长极大值ymax，系统启动时间ini tTime，远端信号向量X， 先验误差信号ε(η)，后验误差信号e(n)，近端信号的能量期望估计巧、远端信号的能量期望 估计驾、估计回声信号的能量期望估计^和误差信号的能量期望估计彩，先验误差信号与 近端信号之间的互相关估计丫 ed，收敛性统计参数η(η)，收敛性统计参数期望值e邱ri(n)，收 敛性统计参数期望阔值thres。
[0033] 步骤3)所述的确定滤波器初始化步长包括如下步骤：
[0034] (1)计算出先验误差信号与近端信号之间的互相关估计丫 ed，近端信号的能量期望 估计巧，估计回声信号的能量期望估计驾:和误差信号的能量期望估计冬
[0035] (2)将步骤（1)中的参数，代入收敛性统计参数ri(n)计算公式，得到收敛性统计参 数η(η)，
[0036]
1.9
[0037] (3)定义收敛性统计参数期望值expri(n)，计算公式是
[0038] θχρη(η)=λ*θχρη(η-1)+(1-λ)*η(η) 1.11
[0039] 在系统启动时间initTime内，若e邱η(η)小于收敛性统计参数期望阔值thres，则 认为系统当前处于远端情况，采取滤波器步长极大值ymax更新滤波器;否则对步长进行更 新。
[0040] 步骤(3)所述的对步长进行更新是采用如下公式：
[0041]
1.6
[0042] 式中，δ、ξ均为一常数。
[0043] 本发明的一种声学回声消除器启动阶段稳健的步长调整方法，不需要DTD模块，在 滤波器收敛速度和系统性能之间取得折中。通过收敛性统计参数期望值在声学回声消除器 启动阶段选择合适的步长调整方式，实现在远端情况下大部分时间里选择大步长更新滤波 器，令其快速收敛，而在其它情况下保持小步长更新，既能保留近端语音，同时也能维持滤 波器快速更新。本发明的方法在双端通话情况不仅能保证滤波器快速收敛到某一稳定状 态，还能在启动阶段完整保留近端语音，提高了系统性能，实用性更强。本发明解决了 VSS-NLMS-UMDT算法在启动阶段的系统容易失调的弊端，在滤波器收敛速度和系统性能之间取 得折中，增强了系统的稳定性与实用性。
【附图说明】
[0044] 图1是一个典型的声学回声消除系统结构；
[0045] 图2a是远端语音信号时域图；
[0046] 图化是近端信号时域图；
[0047] 图2c是近端语音信号时域图；
[004引图3a是VSS-NLMS处理结果；
[0049]图 3b 是 VSS-NLMS-UMDT 处理结果；
[(K)加]图3c是VSS-NLMS-UMDT-C沈处理结果；
[0051 ]图4是本发明实施例的算法流程图。
【具体实施方式】
[0052] 下面结合实施例和附图对本发明的一种声学回声消除器启动阶段稳健的步长调 整方法做出详细说明。
[0053] 本发明的一种声学回声消除器启 动阶段稳健的步长调整方法，所要解决的关键是 声学回声消除器启动阶段如何选择合适的步长调整方法。启动阶段若在远端情况下，选取 大步长更新滤波器，使滤波器快速收敛;否则选取小步长，既能保留近端语音信号同时也会 维持滤波器的更新。
[0054] 本发明的一种声学回声消除器启动阶段稳健的步长调整方法，包括如下步骤：
[0055] 1)确定算法参数；
[0化6]所述的算法参数包括有:包括语音采样频率fs，滤波器长度L，滤波器状态斬的，滤 波器步长μ(η)，滤波器步长极大值ymax，系统启动时间initTime，远端信号向量X，先验误差 信号ε(η)，后验误差信号e(n)，近端信号的能量期望估计驾、远端信号的能量期望估计癸、 估计回声信号的能量期望估计巧和误差信号的能量期望估计勞，先验误差信号与近端信 号之间的互相关估计丫 ed，收敛性统计参数η(η)，收敛性统计参数期望值expri(n)，收敛性统 计参数期望阔值thres。
[0057] 2)先验滤波，所述的先验滤波，包括:去除远端信号向量X中距离η时刻最远的数据 x(n-L)，将当前远端信号χ(η)更新至远端信号向量X;利用滤波器的前一状态?ι(/;-1)估计当 前回声信号，并将估计出的当前回声信号从当前近端信号d(n)中滤除，得到先验误差信号ε (η)；
[005引3)确定滤波器初始化步长；
[0059] 所述的确定滤波器初始化步长包括如下步骤：
[0060] (1)计算出先验误差信号与近端信号之间的互相关估计丫 ed，近端信号的能量期望 估计遊，估计回声信号的能量期望估计和误差信号的能量期望估计勞；
[0061] (2)将步骤（1)中的参数，代入收敛性统计参数ri(n)计算公式，得到收敛性统计参 数η(η)，
[0062]
IJ
[0063] (3)收敛性统计参数ri(n)仅能够代表滤波器的当前状态，跳跃性较大，不能很好地 跟踪描述滤波器。为了更好地统计滤波器的收敛性，因此定义收敛性统计参数期望值expri (η)，计算公式是
[0064] θχρη(η)=λ*θχρη(η-1)+(1-λ)*η(η) 1.11
[0065] 实验研究发现，若滤波器在远端情况下稳定更新，则收敛性统计参数期望e邱ri(n) 大部分时间将处在某一较小值thresW下。因此在系统启动时间initTime内，若e邱ri(n)小 于收敛性统计参数期望阔值thres，则认为系统当前处于远端情况，采取滤波器步长极大值 ymax更新滤波器;否则对步长进行更新。所述的对步长进行更新是采用如下公式：
[0066]
1.6
[0067] 式中，δ、ξ均为一常数。
[0068] 4)进行回声消除，根据步骤3)所选取的步长，利用欠定模型双端通话鲁棒的变步 长归一化最小均方误差算法进行声学回声消除。
[0069] 下面，W语音采样率为16Κ、滤波器长度为1000阶的系统为例实施本专利所提方 案，处理步骤参考图4算法流程图。
[0070] 算法参数设置：语音采样频率f S = 16Κ，滤波器长度L= 1000，滤波器初始状态 h(内)=〇，滤波器初始步长μ ( η ) = 0，常数泣=30*式2，滤波器步长极大值
和滤波器初始时间initTime = 6s,远端信号向量χ = 0,先验误差信 号ε(η)=0,后验误差信号e(n)=0,近端信号、估计回声信号和误差信号的能量期望估计 巧:=巧:=安=0,先验误差信号与近端信号之间的互相关估计丫 ed = 0,收敛性统计参数η(η) =0，收敛性统计参数期望值expri(n) = 0，收敛性统计参数期望阔值thres = 0.05。
[0071 ]本算法具体实施步骤如下：
[0072] 1、读取当前远端信号x(n)和当前近端信号d(n)，将x(n)更新至远端信号向量XL (η)中；
[0073] 2、通过公式1.3和1.4,计算先验误差信号ε(η);
[0074] 3、通过公式l.s-l.io,计算出特征值丫 ed,巧，聲，繁,η(η)和ex阿(n);
[0075] 4、若当前时刻n<initTime*L，则滤波器处在启动阶段，此时若收敛性统计参数期 望值expn ( η )小于收敛性统计参数期望阔值thr e S，则令μ ( η ) = ymax ;否则，
[0076] 1、若11〉= ;[]1；[1:111116礼，贝11系统已启动结束，令
[0077]
[007引 5、将得到的步长代入式1.6得到滤波器状态}!(/?)，代入式1.4和1.3后得到后验误 差信号e(n)。若读取语音已全部处理完毕，算法终止;否则跳转回第1步。
[0079]后验误差信号e(n)即为系统的最终输出。
【主权项】
1. 一种声学回声消除器启动阶段稳健的步长调整方法，其特征在于，包括如下步骤： 1) 确定算法参数； 2) 先验滤波，所述的先验滤波，包括:去除远端信号向量X中距离η时刻最远的数据x(n-L)，将当前远端信号x(n)更新至远端信号向量X;利用滤波器的前一状态估计当前回 声信号，并将估计出的当前回声信号从当前近端信号d(n)中滤除，得到先验误差信号ε(η); 3) 确定滤波器初始化步长； 4) 进行回声消除，根据步骤3)所选取的步长，利用欠定模型双端通话鲁棒的变步长归 一化最小均方误差算法进行声学回声消除。2. 根据权利要求1所述的一种声学回声消除器启动阶段稳健的步长调整方法，其特征 在于，步骤1)所述的算法参数包括有：包括语音采样频率fs，滤波器长度L，滤波器状态 ，滤波器步长Ii(n)，滤波器步长极大值Iimax，系统启动时间ini tTime，远端信号向量X， 先验误差信号ε(η)，后验误差信号e(n)，近端信号的能量期望估计H、远端信号的能量期望 估计式、估计回声信号的能量期望估计f和误差信号的能量期望估计乾，先验误差信号与 近端信号之间的互相关估计γ ed，收敛性统计参数n(n)，收敛性统计参数期望值expn(n)，收 敛性统计参数期望阈值thres。3. 根据权利要求1所述的一种声学回声消除器启动阶段稳健的步长调整方法，其特征 在于，步骤3)所述的确定滤波器初始化步长包括如下步骤： (1) 计算出先验误差信号与近端信号之间的互相关估计yed，近端信号的能量期望估计 句，估计回声信号的能量期望估计有和误差信号的能量期望估计式2; (2) 将步骤（1)中的参数，代入收敛性统计参数η( η)计算公式，得到收敛性统计参数η (η)，(3) 定义收敛性统计参数期望值expn(n)，计算公式是 θχρη(η)=λ*θχρη(η_1)+(1_λ)*η(η) 1.11 在系统启动时间initTime内，若expn(n)小于收敛性统计参数期望阈值thres，则认为 系统当前处于远端情况，采取滤波器步长极大值ymax更新滤波器;否则对步长进行更新。4. 根据权利要求3所述的一种声学回声消除器启动阶段稳健的步长调整方法，其特征 在于，步骤(3)所述的对步长进行更新是采用如下公式：式中，均为一常数。
【专利摘要】一种声学回声消除器启动阶段稳健的步长调整方法，包括：确定算法参数；先验滤波：去除远端信号向量中距离n时刻最远的数据，将当前远端信号更新至远端信号向量；利用滤波器的前一状态估计当前回声信号，并将估计出的当前回声信号从当前近端信号中滤除，得到先验误差信号；确定滤波器初始化步长；进行回声消除，根据所选取的步长，利用欠定模型双端通话鲁棒的变步长归一化最小均方误差算法进行声学回声消除。本发明解决了VSS-NLMS-UMDT算法在启动阶段的系统容易失调的弊端，在滤波器收敛速度和系统性能之间取得折中，增强了系统的稳定性与实用性。
【IPC分类】H04M9/08
【公开号】CN105491256
【申请号】CN201510915783
【发明人】张涛, 焦海泉, 唐伟, 赵鑫
【申请人】天津大学
【公开日】2016年4月13日
【申请日】2015年12月9日
最新回复(0)