一种子带自适应滤波器的权值向量更新方法
【专利摘要】本发明公开了一种子带自适应滤波器的权值向量更新方法，属于数字滤波器设计领域。该方法首先采用子带输入信号和子带误差信号来估计子带误差信号功率和子带系统噪声功率，然后使用估计得到的子带误差信号功率和子带系统噪声功率来估计全带误差信号功率和全带系统噪声功率，最后通过比较全带误差信号功率和全带系统噪声功率的大小来确定重复利用权值向量的阶数。使用该方法的子带自适应滤波器，既可获得较快的收敛速度，还能获得较低的稳态失调。
【专利说明】一种子带自适应滤波器的权值向量更新方法
【技术领域】
[0001]本发明涉及一种自适应滤波器的设计方法，具体地涉及一种子带自适应
[0002]滤波器的权值向量更新方法，属于数字滤波器设计领域。
【背景技术】
[0003]传统的数字滤波器的系数向量是固定的。传统的数字滤波器的主要任务是滤除输入信号中无用的频谱分量，而保留需要的频谱分量，因而其运行的方式是根据输入信号和滤波器的系数向量获得输出信号。与传统的系数向量固定的滤波器不同，自适应滤波器能够未知系统的输入、输出信号，来逼近该未知系统。由于解决系统辨识、回声消除、主动噪声控制、信道均衡、干扰抵消等问题的实质，是根据未知系统的输入和输出信号来求得该未知系统，因而自适应滤波器在噪声抵消器、免提电话、视频电视、电子助听器、主动噪声控制器等设备中获得了广泛应用。
[0004]衡量自适应滤波器性能的主要指标有收敛速度和稳态失调。收敛速度决定了自适应滤波器逼近未知系统需要的时间，而稳态失调决定了逼近未知系统所能达到的精度。归一化最小均方(NLMS)自适应滤波器结构简单，因而得到了广泛使用。但NLMS自适应滤波器在输入为相关信号时，其收敛速度很慢。为此，Lee等人提出了一种子带自适应滤波器[Improving convergence of the NLMS algorithm using constrained subbandupdates, IEEE Signal Processing Letters, 11 (9)，736 - 739，2004]。该滤波器称为归一化子带自适应滤波器(NSAF)。虽然NSAF加快了自适应滤波器的收敛速度，但是具有较大的稳态失调。
[0005]为了降低NSAF的稳态失调，Ni提出了 一种改进的归一化子带自适应滤 波器[Improved normalised subband adaptive filter, ElectronicsLetters, 48 (6),320 - 321，2012]，该滤波器简写为INSAF。INSAF利用权值向量重用的方法来更新自适应滤波器的权值向量，从而降低自适应滤波器的稳态失调。INSAF权值向量重用阶数会同时影响收敛速度的快慢和稳态失调的大小:权值重用阶数越小，收敛越快，但是稳态失调越高；权值重用阶数越大，稳态失调越低，但是收敛越慢。因为INSAF无法兼得快的收敛速度和低的稳态失调。
【发明内容】
[0006]本发明目的是提供一种子带自适应滤波器的权值向量更新方法，该方法解决了改进的归一化子带自适应滤波器(INSAF)不能兼得快的收敛速度和低的稳态失调的问题。
[0007]该发明的原理为:通过比较全带误差信号功率 <(幻与由全带系统噪声功率
《(幻表示的稳态均方误差的大小，来确定子带自适应滤波器所处的阶
段。当自适应滤波器处于自适应初期时，自适应滤波器进行更新时只基于当前时刻的权值向量，从而加快自适应滤波器的收敛速度；当自适应滤波器处于收敛阶段时，自适应滤波器权值向量的更新基于前Rmax个时刻的权值向量，从而降自适应滤波器的低稳态失调。
[0008]本发明公开的一种子带自适应滤波器的权值向量更新方法的方案，具体包含以下步骤:
[0009]步骤1.利用子带自适应滤波器N个子带误差信号ei，D(k)，i=0，l，"％N-l，来估计N个子带误差信号功率i=0, I, - ,N-1 ；
[0010]步骤2.利用子带自适应滤波器N个子带输入信号七(11)，1=0，1，"％^1，和N个子带误差信号ei，D(k)，i=0，I, - ,N-1,来估计N个子带系统噪声的功率? 1=0, I,…,N-1 ；
[0011]步骤3.利用步骤I中得到的子带误差信号的功率(lh i=0, I, -,N-1,来估计全带误差信号的功率 <(幻；
[0012]步骤4.利用步骤2中得到的子带系统噪声功率(幻，?=0，1，...，Ν_1，来估计全带系统噪声的功率
[0013]步骤5.根据步骤3中得到的全带误差信号功率<(衫和步骤4中得到的全带系统噪声功率来判定自适应滤波器权值向量重用的阶数R(k)；
[0014]步骤6.用步骤5中得到的滤波器权值向量重用阶数R(k)来更新自适应滤波器的权值向量。
[0015]优选的，上述步骤I中的子带误差信号功率<n(A), i=0, I,-,N-1,由公
式=+来进行估计，其中λ为平滑因子，该平滑因子在
[0.9,0.999]范围内取值。
[0016]优选的，上述步骤2中的子带系统噪声功率<n W , i=0, I,…，Ν-1，由
【权利要求】
1.一种子带自适应滤波器的权值向量更新方法，其特征在于所述方法包括以下步骤: 步骤1.利用子带自适应滤波器N个子带误差信号ei，D(k)，i=0，l，...，N-l，来估计N个子带误差信号功率，i=0, I,，，N-1 ; 步骤2.利用子带自适应滤波器N个子带输入信号ui(n)，i=0，l，...，N-l，和N个子带误差信号ei，D(k)，i=0，I,，，N-1，来估计N个子带系统噪声的功率gn(A) , i=0, I, -,N-1 ； 步骤3.利用步骤I中得到的子带误差信号的功率?=0, 1,...,Ν-1,来估计全带误差信号的功率 步骤4.利用步骤2中得到的子带系统噪声功率式^(幻，i=0，l，…，Ν-1，来估计全带系统噪声的功率<(衫； 步骤5.根据步骤3中得到的全带误差信号功率<(Α.)和步骤4中得到的全带系统噪声功率#(幻来判定自适应滤波器权值向量重用的阶数R(k)； 步骤6.用步骤5中得到的滤波器权值向量重用阶数R(k)来更新自适应滤波器的权值向量。
2.根据权利要求1所述的子带自适应滤波器的权值向量更新方法，其特征在于:所述步骤I中的子带误差信号功率<#) , i=0, 1，…，N-1，由公式
3.根据权利要求1所述的子带自适应滤波器的权值向量更新方法，其特征在于:所述步骤2中的子带系统噪声功率, i=0, I,…，N-1，由公式
4.根据权利要求1所述的子带自适应滤波器的权值向量更新方法，其特征在于: 所述步骤3中全带误差信号功率由滤波器的所有子带误差信号功率,
5.根据权利要求1所述的子带自适应滤波器的权值向量更新方法，其特征在于: 所述步骤4中全带系统噪声功率《(幻由滤波器的所有子带系统噪声功率si=0, I, - ,N-1,合成，即
6.根据权利要求1所述的子带自适应滤波器的权值向量更新方法，其特征在于:所述步骤5中自适应滤波器权值向量重用的阶数R(k)通过比较全带误差信号功率<(幻与由全带系统噪声功率#(幻表示的稳态均方误差2^)=2^/(2-/1)的大小得到， 即
7.根据权利要求1所述的子带自适应滤波器的权值向量更新方法，其特征在于:所述步骤6中的自适应滤波器权重向量的更新公式为
【文档编号】H03H21/00GK103929150SQ201410117428
【公开日】2014年7月16日 申请日期:2014年3月27日 优先权日:2014年3月27日
【发明者】倪锦根 申请人:苏州大学