[发明专利]一种宽带声学回声消除方法

权利要求书

1.一种宽带声学回声消除方法，其特征在于，该方法包含如下步骤：

步骤1：对长度为K的包含远端信号的输入移位寄存器和包含近端信号的输入移位寄存器中的数据通过M子带均匀IDFT多相分析滤波器组进行滤波，分别得到远端子带数据和近端子带数据；

步骤2：分别对上述远端子带数据和近端子带数据进行IDFT运算，并去除冗余数据，得到M/2+1个远端复子带数据和M/2+1个近端复子带数据；

步骤3：对上述远端复子带数据和近端复子带数据按子带分别进行自适应滤波及系数更新，得到复子带回声估计误差；

步骤4：对上述复子带的回声估计误差进行IDFT运算，得到M个回声估计误差；

步骤5：将上述回声估计误差分别送入M子带均匀IDFT多相合成滤波器组进行合成滤波，将滤波获得的滤波输出数据与长度为K的输出移位寄存器中的值相加，并使用相加获得的值更新输出移位寄存器；

步骤6：将上述输出移位寄存器中的前r点输出，后续K-r点依次前移r点的位置，最后r点位置补零；

步骤7：将上述包含远端信号的输入移位寄存器和包含近端信号的输入移位寄存器中的后K-r点数据前移r点的位置；并将新输入的r点远端信号和近端信号依次分别放入上述输入移位寄存器的后r点的位置；返回步骤1；

其中，M，r为整数，K为M的整数倍，M大于r。

2.如权利要求1所述的宽带声学回声消除方法，其特征在于，所述步骤3包含如下子步骤：

步骤3a：用第1～M/2子带对应的M/2个远端复子带数据更新子带自适应滤波器的状态寄存器；并用远端复子带数据与复数系数向量相乘，得到复子带的回声估计值；

步骤3b：用第1～M/2子带对应的M/2个近端复子带数据减去上述复子带的回声估计值，得到复子带回声估计误差；

步骤3c：用上述复子带回声估计误差和子带自适应滤波器的状态寄存器来更新复数系数向量。

3.如权利要求1所述的宽带声学回声消除方法，其特征在于，所述步骤2中的M/2+1个远端复子带数据和近端复子带数据分别为IDFT结果中的第0～M/2子带对应的远端复子带数据和近端复子带数据。

4.如权利要求2所述的宽带声学回声消除方法，其特征在于，所述复数系数向量的初始值为0。

5.如权利要求2所述的宽带声学回声消除方法，其特征在于，在所述步骤2c中采用如下公式对所述复数系数向量进行更新：

更新后的复数系数向量＝更新前的复数系数向量+增益×回声估计误差^*复数状态向量；

上述^*表示向量乘法；上述增益大于0。

6.如权利要求1所述的宽带声学回声消除方法，其特征在于，所述M子带均匀IDFT多相分析滤波器组E_i(z)为原型滤波器F₀(z)的多相分量；所述M子带均匀IDFT多相合成滤波器组R_i(z)为原型滤波器G₀(z)的多相分量；其中0≤i≤M-1。

7.如权利要求6所述的宽带声学回声消除方法，其特征在于，对于8KHz、16KHz和32KHz的语音数据，所述原型滤波器F₀(z)和G₀(z)的长度分别大于或等于M、2M和3M。

8.如权利要求1所述的宽带声学回声消除方法，其特征在于，所述步骤4中参加IDFT运算的是第0至第M/2子带对应的复子带回声估计误差，其中第0子带对应的复子带回声估计误差为0。

9.如权利要求1所述的宽带声学回声消除方法，其特征在于，K为r的整数倍。

10.如权利要求1所述的宽带声学回声消除方法，其特征在于，所述步骤2中的M/2+1个远端复子带数据和M/2+1个近端复子带数据为IDFT结果中的第0～M/2子带对应的复子带数据。