一种时频记忆子带比例自适应回声消除方法技术

本发明专利技术公开了一种时频记忆子带比例自适应回声消除方法，属于电声自适应回声消除技术领域。针对回声消除问题，提供了一种时频记忆子带比例自适应回声消除方法。该发明专利技术一方面针对输入回声信号的非平稳特性，采用子带结构，在频域上对输入回声信号进行解相关，采用AP算法从时域上进行预白化，获得较好的回声消除性能；另一方面利用带有时频记忆的比例因子，抵抗回声信道的稀疏性，获得较快的收敛速率。

【技术实现步骤摘要】

本专利技术具体涉及一种时频记忆子带比例自适应回声消除方法，属于电声自适应回声消除

技术介绍
在免提电话会议或手机通话中，由于通话环境的原因，经常能从扬声器中听到自己延迟的声音，该延迟的声音称为回声。了消除此种回声，提高通话质量，通常在声学系统内部安装回声消除装置(acousticechocanceller,AEC)，产生与回声幅值相同而相位相反的声学信号，并在声音接收端将接收信号与此同幅反相位的信号相加，从而消除回声干扰。回声消除的核心问题是回声信道的跟踪和回声的估计，理想情况下当AEC产生的回声估计值与回声幅值相同，相位相反时，系统中的回声可以完全被抑制。目前，回声信道跟踪和回声估计通常采用自适应滤波技术，设计各种滤波器(FIR滤波器、IIR滤波器和子带滤波器等)对外部回声信道进行模拟，并采用自适应算法对滤波器的抽头进行实时的更新。因此，滤波器的自适应更新收敛速率和准确性是影响回声消除性能的两个主要因素。滤波器收敛速率过低，回声估计值跟踪能力下降，造成抑制性能降低；滤波器估计准确度过低，导致稳态时残留回声较大，影响正常的通话。语音信号是一种非平稳信号，它的时频域特性会随着时间变化，因此传统的自适应算法，如最小均方(leastmeansquare,LMS)算法，归一化最小均方(Normalizedleastmeansquare,NLMS)算法在语音信号背景下，其收敛速率和稳态精确性能都受到巨大的限制；另外，回声消除的另一个特点是其回声信道的稀疏特性，即回声信道的能量主要集中在少数有意义的抽头权系数上，大量的抽头权系数幅值都接近零，这种稀疏性严重的制约了传统自适应算法的收敛动态性能，针对以上两种回波消除特点，通过对现有技术文献的检索，目前采用的方法有一下几种：参考文献1:在中国专利技术专利申请号201510028006.4，名称为“一种改进的凸组合解相关成比例自适应回声消除方法”中，其针对回声信道的稀疏性，采用滤波器的凸组合结构，结合解相关自适应算法和比例自适应算法，有效的缓解了信道稀疏性对算法收敛速率的影响，获得了较快的收敛速率，低稳态误差以及较好的抗干扰性能，但是由于需要两套滤波器的构成凸组合结构，其耗费资源和计算复杂度都较大，算法的特性对于凸组合混合参数等较为敏感，增加了硬件实现的难度。参考文献2:“Animprovedmultiband-structuredsubbandadaptivefilteralgorithm”(中文名：一种改进的多自带结构子带自适应滤波算法Yang,F.,Wu,M.,Ji,P.,Yang,J.,IEEESignalProcess.Lett.,2012,19,(10),pp.647-650)针对回声信号的非平稳特性，采用子带分割的方法，将信号的频谱划分为多个子带，使得每个子带信号具有白噪声的特性，降低了其非平稳度，另外，采用仿射投影算法(Affineprojection,AP)对每个子带进行滤波器抽头权系数更新，使得整个更新结构具有时频二维特性，有效的提高了回声抑制的性能，但是其并没有考虑是回声信道的稀疏特性，而且子带个数较多的情况会带来巨大的计算复杂度。参考文献3:“Anefficientproportionateaffineprojectionalgorithmforechocancellation”(中文名：一有效的回波消除仿射投影算法Paleologu,C.,Ciochina,S.,Benesty,J.IEEESignalProcess.Lett.,2010,17,(2),pp.165–168)针对回声信号的稀疏特性和回声的非平稳性，采用滤波器抽头加权系数比例更新以及在AP算法在时域的白化作用，提高了算法的收敛速率以及回声的消除性能，但该方法由于未考虑非平稳信号的频域特性，仍存在收敛速率缓慢的问题。
技术实现思路
本专利技术针对上述回声消除问题，提供了一种时频记忆子带比例自适应回声消除方法。方法具体包括以下步骤：步骤一将远端信号x(n)和近端接收到的包含回声的信号d(n)输入至分析滤波器组F0(z),F1(z),…FN-1(z)中，其中n表示时间刻度，Fi(z)i＝1,2,…N-1表示分析滤波器组的Z域传递函数，分析滤波器组包含N个分析滤波器，分析滤波器的长度为M，将信号在频域上均匀分割为N个等带宽的子带远端信号x0(n),x1(n),…xN-1(n)和N个等带宽的子带近端信号d0(n),d1(n),…dN-1(n)；步骤二远端子带信号x0(n),x1(n),…xN-1(n)和近端信号子带d0(n),d1(n),…dN-1(n)进行降低速率的N倍抽取，得到降速后的子带信号x0(p),x1(p),…xN-1(p)、d0(p),d1(p),…dN-1(p)，其中p表示时间刻度，且p＝n/N；步骤三按公式一计算p时刻N个子带回波估值向量，其中Xi(p)＝[xi(p),xi(p-1),…xi(p-D+1)]为第i个输入子带矩阵，其维数为L×D，D是投影阶数，xi(p)＝[xi(pN),xi(pN-1),…,xi(pN-L+1)]T为滤波器第i个子带输入向量，w(p)＝[w0(p),w1(p),…wL-1(p)]T为p时刻滤波器抽头权值向量，其具体数值由p-1时刻回波消除过程中已知，且在起始0时刻，w(0)＝[0,0,…0]T，L为抽头个数，L＝256～1024，T表示共轭转置运算；步骤四按公式二计算第i个子带的误差估计向量，其中ei,D(p)＝[ei,0(pN),ei,1(pN-1),…ei,L-1(pN-L+1)]T是维数为L×1的误差子带向量，di,D(p)＝[di,0(pN),di,1(pN-1),…di,L-1(pN-L+1)]T表示近端子带信号向量；步骤五按公式三计算第i个子带比例因子向量，ci(p)＝[ci,0(p),ci,1(p),…ci,L-1(p)]T公式三ci,l(p)为ci(p)的第i个元素，其中ζ取0.5为比例尺度因子，其取值范围-1＜ζ＜1，σ＝0.001为格式化因子，θi-1,l(p)表示第i-1个子带信号对第l个权系数的增量贡献因子，为比例因子的频域记忆因子，当i＝0时，θi-1,l(p)＝0，同时按照上述方法，计算p-1时刻，p-2时刻,……,p-D+1时刻的子带比例因子向量ci(p-1)，ci(p-2)，……,ci(p-D+1)，D为比例因子的时域记忆刻度；步骤六通过子带输入信号矩阵和比例因子向量按公式四计算比例输入矩阵其中表示向量的hadmard积；步骤七按公式五计算第i个子带权向量增量θi(p)其中λ＝0.001为格式化因子，I为D×D维单位矩阵，λI为对角元素为λ的对角矩阵，其目的是保证计算的稳定性，μ＝0.1为步长因子，其取值范围0＜μ＜2；步骤八回到步骤四，根据步骤四至步骤七计算出所有子带权向量增量，即θi(p)i＝1,2,…N-1；步骤九抽头权值向量按照公式六进行更新，最终得到p+1时刻的权值向量w(p+1)，其维数为L×1；步骤十回到步骤三,计算p+1时刻N个子带回波估向量，其中Xi(p+1)为p+1时刻滤波器第i个子带输入矩阵，其维数为L×D，w(p+1)为步骤九中计算出的p+1时刻滤波器抽头权值向量，其维数为本文档来自技高网...

【技术保护点】
一种时频记忆子带比例自适应回声消除方法，其特征在于，所述方法包括以下步骤：步骤一将远端信号x(n)和近端接收到的包含回声的信号d(n)输入至分析滤波器组F0(z),F1(z),…FN‑1(z)中，其中n表示时间刻度，Fi(z)i＝1,2,…N‑1表示分析滤波器组的Z域传递函数，分析滤波器组包含N个分析滤波器，分析滤波器的长度为M，将信号在频域上均匀分割为N个等带宽的子带远端信号x0(n),x1(n),…xN‑1(n)和N个等带宽的子带近端信号d0(n),d1(n),…dN‑1(n)；步骤二远端子带信号x0(n),x1(n),…xN‑1(n)和近端信号子带d0(n),d1(n),…dN‑1(n)进行降低速率的N倍抽取，得到降速后的子带信号x0(p),x1(p),…xN‑1(p)、d0(p),d1(p),…dN‑1(p)，其中p表示时间刻度，且p＝n/N；步骤三按公式一计算p时刻N个子带回波估值向量，其中Xi(p)＝[xi(p),xi(p‑1),…xi(p‑D+1)]为第i个输入子带矩阵，其维数为L×D，D是投影阶数，xi(p)＝[xi(pN),xi(pN‑1),…,xi(pN‑L+1)]T为滤波器第i个子带输入向量，w(p)＝[w0(p),w1(p),…wL‑1(p)]T为p时刻滤波器抽头权值向量，其具体数值由p‑1时刻回波消除过程中已知，且在起始0时刻，w(0)＝[0,0,…0]T，L为抽头个数，L＝256～1024，T表示共轭转置运算；步骤四按公式二计算第i个子带的误差估计向量，其中ei,D(p)＝[ei,0(pN),ei,1(pN‑1),…ei,L‑1(pN‑L+1)]T是维数为L×1的误差子带向量，di,D(p)＝[di,0(pN),di,1(pN‑1),…di,L‑1(pN‑L+1)]T表示近端子带信号向量；步骤五按公式三计算第i个子带比例因子向量，ci(p)＝[ci,0(p),ci,1(p),…ci,L‑1(p)]T   公式三ci,l(p)为ci(p)的第i个元素，ci,l(p)=1-ζ2L+(1+ζ)|wi(p)+θi-1,i(p)|2Σl=0L-1|wl(p)+θi-1,l(p)|+σ,]]>其中ζ取0.5为比例尺度因子，其取值范围‑1＜ζ＜1，σ＝0.001为格式化因子，θi‑1,l(p)表示第i‑1个子带信号对第l个权系数的增量贡献因子，为比例因子的频域记忆因子，当i＝0时，θi‑1,l(p)＝0，同时按照上述方法，计算p‑1时刻，p‑2时刻,……,p‑D+1时刻的子带比例因子向量ci(p‑1)，ci(p‑2)，……,ci(p‑D+1)，D为比例因子的时域记忆刻度；步骤六通过子带输入信号矩阵和比例因子向量按公式四计算比例输入矩阵其中表示向量的hadmard积；步骤七按公式五计算第i个子带权向量增量θi(p)其中λ＝0.001为格式化因子，I为D×D维单位矩阵，λI为对角元素为λ的对角矩阵，μ＝0.1为步长因子，其取值范围0＜μ＜2；步骤八回到步骤四，根据步骤四至步骤七计算出所有子带权向量增量，即θi(p)i＝1,2,…N‑1；步骤九抽头权值向量按照公式六进行更新，最终得到p+1时刻的权值向量w(p+1)，其维数为L×1；步骤十回到步骤三,计算p+1时刻N个子带回波估向量，其中Xi(p+1)为p+1时刻滤波器第i个子带输入矩阵，其维数为L×D，w(p+1)为步骤九中计算出的p+1时刻滤波器抽头权值向量，其维数为L×1；步骤十一按公式七从包含有回波信号的p+1时刻近端输入di,D(p+1)(i＝1,2,…N‑1)中减去回波干扰估计值，其中ei,D(p+1)是干扰消除后的误差信号，包含近端有用信号和残留回波，其维数为L×1，di,D(p+1)＝[di,0(pN+1),di,1(pN),…di,L‑1(pN‑L)]T表示p+1时刻近端子带信号，其维数为L×1，最终完成回波干扰消除。...

【技术特征摘要】
1.一种时频记忆子带比例自适应回声消除方法，其特征在于，所述方法包括以下步骤：步骤一将远端信号x(n)和近端接收到的包含回声的信号d(n)输入至分析滤波器组F0(z),F1(z),…FN-1(z)中，其中n表示时间刻度，Fi(z)i＝1,2,…N-1表示分析滤波器组的Z域传递函数，分析滤波器组包含N个分析滤波器，分析滤波器的长度为M，将信号在频域上均匀分割为N个等带宽的子带远端信号x0(n),x1(n),…xN-1(n)和N个等带宽的子带近端信号d0(n),d1(n),…dN-1(n)；步骤二远端子带信号x0(n),x1(n),…xN-1(n)和近端信号子带d0(n),d1(n),…dN-1(n)进行降低速率的N倍抽取，得到降速后的子带信号x0(p),x1(p),…xN-1(p)、d0(p),d1(p),…dN-1(p)，其中p表示时间刻度，且p＝n/N；步骤三按公式一计算p时刻N个子带回波估值向量，其中Xi(p)＝[xi(p),xi(p-1),…xi(p-D+1)]为第i个输入子带矩阵，其维数为L×D，D是投影阶数，xi(p)＝[xi(pN),xi(pN-1),…,xi(pN-L+1)]T为滤波器第i个子带输入向量，w(p)＝[w0(p),w1(p),…wL-1(p)]T为p时刻滤波器抽头权值向量，其具体数值由p-1时刻回波消除过程中已知，且在起始0时刻，w(0)＝[0,0,…0]T，L为抽头个数，L＝256～1024，T表示共轭转置运算；步骤四按公式二计算第i个子带的误差估计向量，其中ei,D(p)＝[ei,0(pN),ei,1(pN-1),…ei,L-1(pN-L+1)]T是维数为L×1的误差子带向量，di,D(p)＝[di,0(pN),di,1(pN-1),…di,L-1(pN-L+1)]T表示近端子带信号向量；步骤五按公式三计算第i个子带比例因子向量，ci(p)＝[ci,0(p),ci,1(p),…ci,L-1(p)]T公式三ci,l(p)为ci(p)的第i个元素，ci,l(p)=1-ζ2L+(1+ζ)|wi(p)+θi-1,i(p)|2Σl=0L-1|wl(p)+θi-1,l(p)|+σ,]]>其中...

【专利技术属性】
技术研发人员：刘畅，张志，唐校，王彩申，
申请(专利权)人：东莞理工学院，
类型：发明
国别省市：广东;44