一种双端通话鲁棒的变阶变步长仿射投影方法技术

一种双端通话鲁棒的变阶变步长仿射投影方法，包括：初始化滤波器的参数；更新输入信号与麦克风信号；计算先验信息；计算当前时刻的麦克风信号功率估计、回声信号功率估计和误差信号功率估计；计算步长；长自信息及其相关计算；确定仿射投影阶数；更新滤波器系数。本发明专利技术用于声学回声消除器，能够使声学回声消除器在双端讲话时选择较小的仿射投影阶数，而在收敛时选择较大的仿射投影阶数。因此本方法提出的方法，可以有效解决收敛速度与双端讲话鲁棒性之间的矛盾，从而获得更好的性能。本方法的鲁棒性更强且没有降低其它性能。

【技术实现步骤摘要】

本专利技术涉及一种仿射投影方法。特别是涉及一种针对于声学回声消除中双端通话问题自适应滤波的双端通话鲁棒的变阶变步长仿射投影方法。
技术介绍
声学回声主要是指从扬声器发出的声音被麦克风重新拾取后再次发送到远端形成的回声。声学回声消除器(AcousticEchoCancellation,AEC)广泛应用于各种嵌入式设备和各种VoIP应用中，包括各种电信网络设备和终端设备、丰富的软件视频会议和VoIP软件电话等。一般的声学回声消除器至少包含双端通话检测模块(DoubleTalkDetector,DTD)、线性回声消除两个部分。图1是一个典型的声学回声消除器。麦克风在时刻n接收的信号为：d(n)＝y(n)+v(n)+w(n)1.1式中：d(n)代表麦克风信号，y(n)代表输入信号x(n)经过扬声器播放后产生的回声信号，v(n)代表近端语音信号，w(n)代表近端噪声信号。输入信号x(n)经过系统传递函数h滤波后形成回声y(n)＝x(n)T*h1.2其中，h＝[h0(n),h1(n),...,hN(n)]T1.3x(n)＝[x(n),x(n-1),...,x(n-N+1)]T1.4x(n)是输入信号向量，N为回声路径长度，T代表矩阵的转置。回声消除的目标就是设计一个自适应有限冲击响应滤波器估计出介于麦克风与扬声器之间的回声路径然后根据这一估计路径得到回声估计值将其从d(n)中消除，从而保留v(n)。e(n)代表线性回声滤波器消除后得到的误差信号，其中，L为自适应滤波器长度，实际中一般L<N。一般认为，声学回声消除器处理的场景分为三种情况：远端情况，只存在回声信号而不存在近端语音信号；近端情况，不存在回声，只存在近端语音信号；双端通话情况，回声信号与近端语音信号同时存在。线性回声消除器主要依靠自适应算法来估计回声路径，由于仿射投影算法(AffineProjectionAlgorithm,APA)能够在收敛速度和计算复杂度之间取得较好折衷而被广泛应用，其更新方程是其中μ为自适应滤波器的步长，为单位矩阵，Kmax为投影阶数最大值，为输入信号矩阵，δ为以常数称为正则化因子，为前验误差信号向量，计算方法是其中d(n)＝[d(n),...,d(n-Kmax+1)]T。近几年，有人考虑了近端信号的变化，提出将步长μ替换为对角矩阵μ(n)，表示为：其更新方程是其中l＝0,1,...,Kmax-1，λ＝1-1/(6*N)。代表麦克风信号功率估计，代表回声信号功率估计，代表误差信号功率估计。该算法虽然对于近端信号有一定的鲁棒性，但是随着仿射投影阶数的增大，该方法的收敛速度将提升，而鲁棒性会变差。
技术实现思路
本专利技术所要解决的技术问题是，提供一种在声学回声消除中自适应地调整仿射投影阶数，使之既能取得较快的收敛速度，又能取得较好的双端通话鲁棒性的双端通话鲁棒的变阶变步长仿射投影方法。本专利技术所采用的技术方案是：一种双端通话鲁棒的变阶变步长仿射投影方法，包括如下步骤：1)初始化滤波器的参数：将滤波器系数麦克风信号功率估计回声信号功率估计误差信号功率估计均初始化为0，对步长μ、正则化因子δ、双端讲话状态仿射投影阶数Ks和最大仿射投影阶数Kmax进行初始化；2)更新输入信号与麦克风信号：获取当前时刻麦克风信号d(n)和输入信号x(n)，并更新至麦克风信号向量d(n)和输入信号向量x(n)中，其中，所述的麦克风信号向量d(n)表示为d(n)＝[d(n),...,d(n-Kmax+1)]T，所述的输入信号向量x(n)表示为x(n)＝[x(n),x(n-1),...,x(n-N+1)]T，N为回声路径长度，T代表矩阵的转置；3)计算先验信息：根据前一时刻估计的滤波器系数计算出先验误差信号矩阵和回声信号所述的先验误差信号矩阵其中，是输入信号矩阵，所述的回声信号计算当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功率估计4)计算步长：根据当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功率估计获得当前时刻的步长矩阵μ(n)如下：其中，式中，l＝0,1,...,Kmax-1；5)长自信息及其相关计算：利用步骤4)得到的当前时刻的步长矩阵μ(n)，计算步长自信息量I(μl(n))，l＝0,1,...,Kmax-1，然后求出第一个步长自信息量在整个步长矩阵中所占比例prop0(n)，6)确定仿射投影阶数：若第一个步长自信息量所占比例prop0(n)小于1/Kmax-ε，则认为自适应滤波器处于双端讲话状态，令仿射投影阶数Kn＝Ks，其中ε为常数称为约束因子；7)更新滤波器系数：根据仿射投影阶数Kn，选择对应长度的输入信号矩阵先验误差信号矩阵和步长矩阵代入下式进行滤波器的更新：步骤3)所述的计算当前时刻的麦克风信号功率估计公式如下：所述的计算当前时刻的回声信号功率估计公式如下：所述的计算当前时刻的误差信号功率估计公式如下：其中，l＝0,1,...,Kmax-1，λ＝1-1/(6*N)。本专利技术的一种双端通话鲁棒的变阶变步长仿射投影方法，用于声学回声消除器，能够使声学回声消除器在双端讲话时选择较小的仿射投影阶数，而在收敛时选择较大的仿射投影阶数。因此本方法提出的方法，可以有效解决收敛速度与双端讲话鲁棒性之间的矛盾，从而获得更好的性能。本方法的鲁棒性更强且没有降低其它性能。附图说明图1是典型的声学回声消除系统；图2是本发吸的方法与其它方法的性能对比；图3是本发吸的方法与其它方法的性能对比。具体实施方式下面结合实施例和附图对本专利技术的一种双端通话鲁棒的变阶变步长仿射投影方法做出详细说明。本专利技术的一种双端通话鲁棒的变阶变步长仿射投影方法，主要解决在声学回声消除中滤波器受双端讲话影响而发散问题。本专利技术通过定义的步长自信息量自适应地调整仿射投影阶数，使之能够在双端讲话时维持较小的仿射投影阶数，获得较好的鲁棒性。而在收敛时使用较大的仿射投影阶数，获得较快的收敛速度。本专利技术的一种双端通话鲁棒的变阶变步长仿射投影方法，包括如下步骤：1)初始化滤波器的参数：将滤波器系数麦克风信号功率估计回声信号功率估计误差信号功率估计均初始化为0，对步长μ、正则化因子δ、双端讲话状态仿射投影阶数Ks和最大仿射投影阶数Kmax进行初始化；本专利技术实施例设定对数：滤波器系数近端信号功率估计回声信号功率估计误差信号功率估计双端讲话仿射投影阶数Ks＝2，约束因子ε＝1/(5Kmax)，最大仿射投影阶数Kmax＝8，滤波器长度N＝1024，正则化因子其中λ＝1-1/(6N)。近端信号中加入信噪比20dB的独立平稳高斯白噪声。2)更新输入信号与麦克风信号：获取当前时刻麦克风信号d(n)和输入信号x(n)，并更新至麦克风信号向量d(n)和输入信号向量x(n)中，其中，所述的麦克风信号向量d(n)表示为d(n)＝[d(n),...,d(n-Kmax+1)]T，所述的输入信号向量x(n)表示为x(n)＝[x(n),x(n-1),...,x(n-N+1)]T，N为回声路径长度，T代表矩阵的转置；3)计算先验信息：根据前一时刻估计的滤波器系数计算出先验误差信号矩阵和回声信号所述的先验误差信号矩阵其中，是输入信号矩阵，所述的回声信号计算当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功本文档来自技高网...

【技术保护点】
一种双端通话鲁棒的变阶变步长仿射投影方法，其特征在于，包括如下步骤：1)初始化滤波器的参数：将滤波器系数麦克风信号功率估计回声信号功率估计误差信号功率估计均初始化为0，对步长μ、正则化因子δ、双端讲话状态仿射投影阶数Ks和最大仿射投影阶数Kmax进行初始化；2)更新输入信号与麦克风信号：获取当前时刻麦克风信号d(n)和输入信号x(n)，并更新至麦克风信号向量d(n)和输入信号向量x(n)中，其中，所述的麦克风信号向量d(n)表示为d(n)＝[d(n),...,d(n‑Kmax+1)]T，所述的输入信号向量x(n)表示为x(n)＝[x(n),x(n‑1),...,x(n‑N+1)]T，N为回声路径长度，T代表矩阵的转置；3)计算先验信息：根据前一时刻估计的滤波器系数计算出先验误差信号矩阵和回声信号所述的先验误差信号矩阵其中，是输入信号矩阵，所述的回声信号计算当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功率估计4)计算步长：根据当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功率估计获得当前时刻的步长矩阵μ(n)如下：μ(n)=diag(μ0(n),μ1(n),...,μKmax-1(n))]]>其中，式中，l＝0,1,...,Kmax‑1；5)长自信息及其相关计算：利用步骤4)得到的当前时刻的步长矩阵μ(n)，计算步长自信息量I(μl(n))，I(μl(n))=log2μl(n)Σi=0Kmax-1μi(n)]]>l＝0,1,...,Kmax‑1，然后求出第一个步长自信息量在整个步长矩阵中所占比例prop0(n)，prop0(n)=I(μ0(n))Σi=0Kmax-1I(μi(n));]]>6)确定仿射投影阶数：若第一个步长自信息量所占比例prop0(n)小于1/Kmax‑ε，则认为自适应滤波器处于双端讲话状态，令仿射投影阶数Kn＝Ks，其中ε为常数称为约束因子；7)更新滤波器系数：根据仿射投影阶数Kn，选择对应长度的输入信号矩阵先验误差信号矩阵和步长矩阵代入下式进行滤波器的更新：h^(n)=h^(n-1)+μ(n)XKn(n)[XKnT(n)XKn(n)+δIKn]-1eKn(n).]]>...

【技术特征摘要】
1.一种双端通话鲁棒的变阶变步长仿射投影方法，其特征在于，包括如下步骤：1)初始化滤波器的参数：将滤波器系数麦克风信号功率估计回声信号功率估计误差信号功率估计均初始化为0，对步长μ、正则化因子δ、双端讲话状态仿射投影阶数Ks和最大仿射投影阶数Kmax进行初始化；2)更新输入信号与麦克风信号：获取当前时刻麦克风信号d(n)和输入信号x(n)，并更新至麦克风信号向量d(n)和输入信号向量x(n)中，其中，所述的麦克风信号向量d(n)表示为d(n)＝[d(n),...,d(n-Kmax+1)]T，所述的输入信号向量x(n)表示为x(n)＝[x(n),x(n-1),...,x(n-N+1)]T，N为回声路径长度，T代表矩阵的转置；3)计算先验信息：根据前一时刻估计的滤波器系数计算出先验误差信号矩阵和回声信号所述的先验误差信号矩阵其中，是输入信号矩阵，所述的回声信号计算当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功率估计4)计算步长：根据当前时刻的麦克风信号功率估计回声信号功率估计和误差信号功率估计获得当前时刻的步长矩阵μ(n)如下：μ(n)=diag(μ0(n),μ1(n),...,μKmax-1(n))]]>其中，式中，l＝0,1,...,Kmax-1；5)长自信息及其相关计算：利用步骤4)得到的当前时刻的步长矩阵μ(n)，计算步长自信息量I(μl(n))，I(μl(n))=log2μl(n)Σi=0Kmax-1μi(n)]]...

【专利技术属性】
技术研发人员：张涛，焦海泉，任相赢，安鑫祺，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12