[发明专利]一种应用于回声消除的自适应子带滤波方法、装置

权利要求书

1.一种应用于回声消除的自适应子带滤波方法，其特征在于，包括：

将输入信号x(t)和近端信号d(t)进行子带分割和抽样,并进行离散傅立叶变换得到频域分析数据X_k(nR)和频域近端数据D_k(nR)；

将所述频域分析数据X_k(nR)进行分块组合，得到B个分析数据块；

将所述B个分析数据块输入自适应滤波器进行滤波，得到预测数据结果Y(nR)；

计算所述预测数据结果Y(nR)与频域近端子带数据D(nR)之间的误差数据E(nR)，其中，频域近端子带数据D(nR)＝[D₁(nR),D₂(nR),...,D_M-1(nR),D_M(nR)]^T；

根据所述误差数据E(nR)更新所述自适应滤波器的系数；

根据所述误差数据E(nR)合成时域输出信号y(t)。

2.根据权利要求1所述的应用于回声消除的自适应子带滤波方法，其特征在于，所述将输入信号x(t)和近端信号d(t)进行子带分割和抽样,并进行离散傅立叶变换得到频域分析数据X_k(nR)和频域近端数据D_k(nR)的步骤，包括：

将输入信号x(t)和近端信号d(t)通过分析窗函数h(t)的反转函数进行加权，并进行大小为M的块分割，分块数为K＝P/M，对分块数据叠加求和，并对求和数据做离散傅立叶变换(DFT)运算，得到频域分析数据X_k(nR)和频域近端数据D_k(nR)，具体地，

其中n为抽样时间系数，k为频点,

h(t)为分析窗函数，长度为P,

M为子带数，R为抽样因子，

K＝P/M。

3.根据权利要求2所述的应用于回声消除的自适应子带滤波方法，其特征在于，所述将所述频域分析数据X_k(nR)进行分块组合，得到B个分析数据块的步骤中，第b个分析数据块为：

S_b(nR)＝S_b-1(nR)，其中，

S₁(nR)＝[X₁(nR),X₂(nR),...,X_M-1(nR),X_M(nR)]^T。

4.根据权利要求2所述的应用于回声消除的自适应子带滤波方法，其特征在于，所述将所述分析数据输入到自适应滤波器进行滤波，得到预测数据结果Y(nR)的步骤中，自适应滤波器通过下述滤波算法进行滤波：

其中

H_b(nR)＝diag{H_b(nR)}

H_b(nR)＝DFT([h_nR(bN),...,h_nR(bN+N-1),0_1×(M-N)]^T),H_b(nR)为对所述自适应滤波器系数做分块处理后的滤波器系数子块，N为所述滤波器系数子块的大小，N＝L/B，L为所述自适应滤波器系数的长度。

5.根据权利要求4所述的应用于回声消除的自适应子带滤波方法，其特征在于，所述根据所述误差数据E(nR)更新所述自适应滤波器的系数，根据下述算法实现：

其中，

M_b(nR)＝diag{M_b(nR)}

M_b(nR)＝α/[R_b(nR)+ρ]

S_b(nR)＝diag{S_b(nR)}

β为平滑因子,0β1，

α为收敛因子,2α0，ρ为正则化因子,0ρ1。