[go: up one dir, main page]

CN105654959B - 一种自适应滤波的系数更新方法及装置 - Google Patents

一种自适应滤波的系数更新方法及装置 Download PDF

Info

Publication number
CN105654959B
CN105654959B CN201610046280.9A CN201610046280A CN105654959B CN 105654959 B CN105654959 B CN 105654959B CN 201610046280 A CN201610046280 A CN 201610046280A CN 105654959 B CN105654959 B CN 105654959B
Authority
CN
China
Prior art keywords
filter
coefficient
input signal
segment
coefficients
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610046280.9A
Other languages
English (en)
Other versions
CN105654959A (zh
Inventor
文昊翔
洪远泉
罗欢
周永明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaoguan University
Original Assignee
Shaoguan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shaoguan University filed Critical Shaoguan University
Priority to CN201610046280.9A priority Critical patent/CN105654959B/zh
Publication of CN105654959A publication Critical patent/CN105654959A/zh
Application granted granted Critical
Publication of CN105654959B publication Critical patent/CN105654959B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/16—Vocoder architecture
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/16—Vocoder architecture
    • G10L19/18—Vocoders using multiple modes
    • G10L19/24—Variable rate codecs, e.g. for generating different qualities using a scalable representation such as hierarchical encoding or layered encoding

Landscapes

  • Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Quality & Reliability (AREA)
  • Filters That Use Time-Delay Elements (AREA)

Abstract

本发明公开了一种自适应滤波的系数更新方法，包括步骤1：将k时刻的滤波器W(k)的长度等分为M段；步骤2：求取每段滤波器的系数的二次范数；步骤3：将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器系数的二次范数进行累加以获得最大累加值；步骤4：以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余则为非活跃系数分段Wn‑a(k)；步骤5：更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn‑a(k)。本发明更新滤波器的全部活跃系数并循环更新一小段非活跃系数，降低滤波器有效长度，提高收敛速度，降低了算法复杂度。本发明同时提供一种自适应滤波的系数更新装置。

Description

一种自适应滤波的系数更新方法及装置
技术领域
本发明涉及自适应降噪领域，尤其涉及一种自适应滤波的系数更新方法及装置。
背景技术
随着近年来技术的进步，通信技术日益发展，各种新兴通信层出不穷，但最主要的通信方式仍然是语音通信，用户也越来越关心语音通信质量。而在使用电子、网络通信设备或免提通信设备时，使用者经常会在接收端听到自己的声音，这种现象称为通信回声现象。通信回声严重影响通信质量，在极端情况下可引起啸叫，使通信无法进行。针对通信回声现象，常用的解决方案是在通信设备内置回声消除(echo cancellation，简称EC)系统以抑制回声能量，提高通信质量。
回声消除系统可简化为自适应理论中的系统辨识问题，其核心是自适应滤波器，通过自适应算法对滤波器系数进行自适应更新，以实现最佳滤波。请参阅图1，其是自适应系统辨识原理图。若将声学传播信道的冲激响应序列记为H＝[h1，h2，…，hN]，其中N为序列长度，设k时刻输入信号为X(k)＝[x(k-1)，x(k-2)，…，x(k-N)]。X(k)与H卷积加上背景噪声v(k)干扰得到回声信号d(k)，其中d(k)的计算公式如下：
d(k)＝X(k)HT+v(k)
为消除d(k)，以自适应滤波器W(k)＝[w1(k)，w2(k)，。。。，wN(k)]自适应辨识冲激响应H，以X(k)与W(k)的卷积结果y(k)作为d(k)的估计值，两者相减以抑制回声能量。
为辨识冲激响应H，需要对滤波器的系数进行自适应更新，通常采用归一化最小均方(NorMalized Least Mean Square，简称NLMS)算法对所有系数进行迭代更新，具体的系数迭代更新方程为：
e(k)＝d(k)-W(k)XT(k)；
Figure GDA0002225264410000011
式中，μ为收敛因子，取值范围为0<μ<1；
Figure GDA0002225264410000012
为输入信号序列X(k)的二次范数，计算公式为||X(k)||＝X(k)XT(k)，
Figure GDA0002225264410000013
作为NLMS算法的归一化因子，可通过自回归方式得到，自回归算法如下：
Figure GDA0002225264410000014
综上可知，NLMS算法迭代一次约需算2N次加法、2N次乘法。
因声学信道延时的不确定性，为保证充分辨识冲激响应序列，在某些极端条件下，自适应滤波器需要上千个系数，而过长的滤波器将导致自适应系统收敛速度下降、计算复杂度增加，进而将严重影响自适应系统的性能。
发明内容
本发明在于克服现有技术的缺点与不足，提供一种降低滤波器有效长度、提高系统收敛速度、降低算法复杂度的自适应滤波的系数更新方法。
本发明是通过以下技术方案实现的：一种自适应滤波的系数更新方法，包括如下步骤：
步骤1：将k时刻的滤波器W(k)的长度等分为M段；其中M为大于1的整数；
步骤2：求取每段滤波器的系数的二次范数；
步骤3：将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器系数的二次范数进行累加以获得最大累加值；
步骤4：以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余滤波器分段则为非活跃系数分段Wn-a(k)；
步骤4A：将k时刻的输入信号X(k)等分成M段，且将输入信号亦对应地分为输入信号活跃系数分段Xa(k)和输入信号非活跃系数分段Xn-a(k)，并代入滤波器系数更新公式对滤波器系数进行更新；滤波器系数更新方程为：
e(k)＝d(k)-W(k)XT(k)；
Figure GDA0002225264410000021
Figure GDA0002225264410000022
其中，μ为收敛因子，取值范围为0<μ<1；X(k)为输入信号，
Figure GDA0002225264410000023
为输入信号序列X(k)的二次范数，d(k)为回声信号，e(k)为输出信号；Wi(k)为从Wn-a(k)中循环提取的一段滤波器分段；对应地，Xi(k)为从Xn-a(k)中循环提取的一段输入信号分段；
步骤5：对滤波器系数进行迭代更新，且每次迭代时更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn-a(k)，以获得下一时刻k+1的滤波器系数。
相比于现有技术，本发明提出的自适应滤波器系数更新方法利用声学传播信道的稀疏性，在滤波器内通过移动窗积分法对信道进行延时估计以定位活跃系数位置，并通过精确辨识全部活跃系数完成声学传播信道的自适应辨识。进一步地，在每次迭代时更新滤波器的全部活跃系数，使得活跃系数更新频率得到提高，进而加快系统收敛速度；通过循环更新一段非活跃系数以降低非活跃系数更新频率，进而降低滤波器有效长度；同时，非活跃系数仍能获得一定的更新频率，从而确保了系统的跟踪能力。本发明提出的自适应滤波器系数更新方法可大幅降低滤波器有效长度，提高系统收敛速度，降低了算法复杂度。
进一步地，在步骤3中，获得每段滤波器系数的二次范数并组成新序列，将长度为L的矩形窗在该新序列上移动，每移动一个位置，都对窗口内各分段滤波器系数的二次范数进行累加，以获得最大累加值；其中L为大于1的整数。
本发明同时还提供一种自适应滤波的系数更新装置，包括滤波器，所述滤波器包括
——滤波器分段模块，用于将k时刻的滤波器W(k)的长度等分为M段；
——二次范数获取模块，用于求取每段滤波器系数的二次范数；
——最大累加值获取模块，用于将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器系数的二次范数进行累加以获得最大累加值；其中M为大于1的整数；
——滤波器活跃系数分段和非活跃系数分段获取模块，用于以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余滤波器分段则为非活跃系数分段Wn-a(k)；
——滤波器系数更新模块，用于对滤波器系数进行迭代更新，且每次迭代时用于更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn-a(k),以获得下一时刻k+1的滤波器系数。
相比于现有技术，本发明提出的自适应滤波器系数更新方法利用声学传播信道的稀疏性，在滤波器内通过移动窗积分法对信道进行延时估计以定位活跃系数位置，并通过精确辨识全部活跃系数完成声学传播信道的自适应辨识。进一步地，在每次迭代时更新滤波器的全部活跃系数，使得活跃系数更新频率得到提高，进而加快系统收敛速度；通过循环更新一段非活跃系数以降低非活跃系数更新频率，进而降低滤波器有效长度；同时，非活跃系数仍能获得一定的更新频率，从而确保了系统的跟踪能力。本发明提出的自适应滤波器系数更新方法可大幅降低滤波器有效长度，提高系统收敛速度，降低了算法复杂度。
为了更好地理解和实施，下面结合附图详细说明本发明。
附图说明
图1是自适应系统辨识原理图；
图2是典型的采样频率为8k Hz的声学传播信道冲激响应序列的坐标图；
图3是本发明自适应滤波的系数更新方法的流程图；
图4是与冲激序列位置相对应的滤波器活跃系数分段位置坐标图；
图5是本发明自适应滤波的系数更新装置的模块图。
具体实施方式
请参阅图2，其是典型的采样频率为8k Hz的声学传播信道冲激响应序列的坐标图，其中，横坐标代表时间K，纵坐标代表幅度A，符号t表示延时时间，s表示活跃系数区间。因声学传播信道主要由幅度衰减、声延时与混响构成，导致冲激响应序列具有典型的稀疏性，即序列中大部分系数幅值为零值或极小值以模拟延时；产生声强能量的大幅值系数称为活跃系数，因幅度衰减与混响效应，活跃系数在时域聚集，且数量只占系数总量极少一部分。根据稀疏性，长延时声学传播信道只有活跃系数产生声强能量；对于小幅值、零值系数非活跃系数仅用于模拟延时，并不产生能量，因此仅需辨识活跃系数即可完成对声学传播信道的精确辨识。
本发明基于声学传播信道稀疏性的性质，通过自适应滤波器对滤波器的长度和输入的输入信号进行分段，确定活跃系数和非活跃系数的位置，并利用滤波器滤除后的输出信号e(k)构造迭代更新方程，实现对滤波器系数的自适应更新。
请参阅图3，其是本发明自适应滤波的系数更新方法的流程图。该自适应滤波的系数更新方法，包括如下步骤：
步骤1：将k时刻的滤波器W(k)的长度等分为M段；其中M为大于1的整数。
获取当前时刻k总长度为N的滤波器W(k)并等分为M段，则该滤波器系数可表示为W(k)＝[W1(k)，W2(k)，…，WM(k)]，每个分段长度为N1＝N/M，则第i段滤波器系数可表示为Wi(k)＝[WiN/M(k)，W(iN/M)+1(k)，…，W((i+1)N/M)-1(k)]，其中1≤i≤M，且i为整数。
步骤2：求取每段滤波器的系数的二次范数；
利用活跃系数幅值远大于非活跃系数的特性，通过二次范数以非线性方式突出活跃系数，可有效区分活跃系数与非活跃系数，提高活跃系数定位精度。因此，本发明先对每段滤波器的系数求取二次范数，组成新序列
Figure GDA0002225264410000041
其中，二次范数定义为
Figure GDA0002225264410000042
步骤3：将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器的二次范数进行累加以获得最大累加值。
将长度为L的矩形窗在A(k)序列上移动，其中N1<L<N，即该矩形窗包含了L个滤波器分段；矩形窗每移动一个位置，均对窗口内所有滤波器分段的二次范数进行累加，其累加值可记为：
Figure GDA0002225264410000043
式中，j为矩形窗首位置在A(k)上的坐标。
步骤4：以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余滤波器分段则为非活跃系数分段Wn-a(k)。
根据稀疏性，活跃系数在时域聚集，且其幅值远大于非活跃系数，因此，以最大累加值对应的窗口位置为活跃系数位置。设活跃系数位置包含的首块滤波器分段序号为b，则
Figure GDA0002225264410000051
则活跃系数分段可确定为：
Wa(k)＝[Wb(k),Wb+1(k),…,Wb+L-1(k)]
如设N＝1024、M＝8、L＝3，可得活跃系数分段为Wa(k)＝[W5(k),W6(k),W7(k)]、非活跃系数分段为Wn-a(k)＝[W1(k),W2(k),W3(k),W4(k),W8(k)]。
步骤4A：将k时刻的输入信号X(k)等分成M段，且将输入信号亦对应地分为输入信号活跃系数分段Xa(k)和输入信号非活跃系数分段Xn-a(k)，并代入滤波器系数更新公式对滤波器系数进行更新。
步骤5：对滤波器系数进行迭代更新，且每次迭代时更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn-a(k)，以获得下一时刻k+1的滤波器系数。
当声学信道发生剧变时，特别是活跃系数在时域发生大幅度的位移时，仅对活跃系数分段Wa(k)进行自适应算法将无法跟踪声学信道的突变，为使算法具备快速跟踪能力，必须再循环更新一段非活跃系数。
具体的，迭代时从Wn-a(k)中循环提取一段滤波器分段并记为Wi(k)；对应地，从Xn-a(k)中循环提取一段输入信号分段并记为Xi(k)；则滤波器系数更新方程为：
e(k)＝d(k)-W(k)XT(k)；
Figure GDA0002225264410000052
Figure GDA0002225264410000053
其中，μ为收敛因子，取值范围为0<μ<1；X(k)为输入信号，
Figure GDA0002225264410000054
为输入信号序列X(k)的二次范数，d(k)为回声信号，e(k)为输出信号。
下面举个具体的例子来说明上述更新方程：请参阅图4，其是与冲激序列位置相对应的滤波器活跃系数分段位置坐标图。假设活跃系数分段位置已知为Wa(k)＝[W5(k)W6(k)]，则每次迭代均更新Wa(k)与循环更新一段Wi(k)，即每次迭代，依次更新以下滤波器分段：第一次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W1(k)；第二次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W2(k)；第三次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W3(k)；第四次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W4(k)；第五次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W7(k)；第六次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W8(k)。
本发明提出的系数更新算法，只需要进行6次迭代即可完成对滤波器全部系数的更新。在6次迭代中活跃系数更新6次，非活跃系数只更新1次。可见，本发明提出的自适应滤波器系数更新方法中，活跃系数的更新频率为1，非活跃系数的更新频率降为1/(M-L)，通过降低非活跃系数更新频率以降低滤波器有效长度，进而最终达到加快收敛速度、降低计算复杂度的目的
相比于现有技术，本发明提出的自适应滤波器系数更新方法利用声学传播信道的稀疏性，在滤波器内通过移动窗积分法对信道进行延时估计以定位活跃系数位置，并通过精确辨识全部活跃系数完成声学传播信道的自适应辨识。进一步地，在每次迭代时更新滤波器的全部活跃系数，使得活跃系数更新频率得到提高，进而加快系统收敛速度；通过循环更新一段非活跃系数以降低非活跃系数更新频率，进而降低滤波器有效长度；同时，非活跃系数仍能获得一定的更新频率，从而确保了系统的跟踪能力。本发明提出的自适应滤波器系数更新方法可大幅降低滤波器有效长度，提高系统收敛速度，降低了算法复杂度。进一步地，本发明提出的滤波系数更新方法用一个滤波器完成延时估计与活跃系数辨识，可有效避免双滤波器结构信息冗余的缺陷、节省计算复杂度。与传统的NLMS算法相比，新算法可大幅降低滤波器有效长度，提高系统收敛速度，降低了计算复杂度。
请参阅图5，其是本发明自适应滤波的系数更新装置的模块图。本发明同时还提供一种自适应滤波的系数更新装置，包括滤波器，所述滤波器包括：
——滤波器分段模块，用于将k时刻的滤波器W(k)的长度等分为M段；
——二次范数获取模块，用于求取每段滤波器系数的二次范数；
——最大累加值获取模块，用于将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器系数的二次范数进行累加以获得最大累加值；其中M为大于1的整数；
——滤波器活跃系数分段和非活跃系数分段获取模块，用于以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余滤波器分段则为非活跃系数分段Wn-a(k)；
——滤波器系数更新模块，用于对滤波器系数进行迭代更新，且每次迭代时用于更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn-a(k),以获得下一时刻k+1的滤波器系数。
所述滤波器分段模块获取当前时刻k总长度为N的滤波器W(k)并等分为M段，则该滤波器系数可表示为W(k)＝[W1(k)，W2(k)，…，WM(k)]，每个分段长度为N1＝N/M，则第i段滤波器系数可表示为Wi(k)＝[WiN/M(k)，W(iN/M)+1(k)，…，W((i+1)N/M)-1(k)]，其中1≤i≤M，且i为整数。
利用活跃系数幅值远大于非活跃系数的特性，通过二次范数以非线性方式突出活跃系数，可有效区分活跃系数与非活跃系数，提高活跃系数定位精度。因此，本发明通过所述二次范数获取模块对每段滤波器的系数求取二次范数，将求取的二次范数组成新序列
Figure GDA0002225264410000071
其中，二次范数定义为
Figure GDA0002225264410000072
所述最大累加值获取模块用于将长度为L的矩形窗在A(k)序列上移动，其中N1<L<N，即该矩形窗包含了L个滤波器分段；矩形窗每移动一个位置，均对窗口内所有滤波器分段的二次范数进行累加，其累加值可记为：
Figure GDA0002225264410000073
式中，j为矩形窗首位置在A(k)上的坐标。
根据稀疏性，活跃系数在时域聚集，且其幅值远大于非活跃系数。因此所述滤波器活跃系数分段和非活跃系数分段获取模块因此，以最大累加值对应的窗口位置为活跃系数位置。设活跃系数位置包含的首块滤波器分段序号为b，则
Figure GDA0002225264410000074
则活跃系数分段可确定为：
Wa(k)＝[Wb(k),Wb+1(k),…,Wb+L-1(k)]
如设N＝1024、M＝8、L＝3，可得活跃系数分段为Wa(k)＝[W5(k),W6(k),W7(k)]、非活跃系数分段为Wn-a(k)＝[W1(k),W2(k),W3(k),W4(k),W8(k)]。
本发明提供的自适应滤波的系数更新装置进一步还包括输入信号处理模块，所述输入信号处理模块包括输入信号分段模块、以及输入信号活跃系数分段和非活跃系数分段获取模块。所述输入信号分段模块用于将k时刻的输入信号X(k)等分成M段；所述输入信号活跃系数分段和非活跃系数分段获取模块用于将输入信号也分为输入信号活跃系数分段Xa(k)和输入信号非活跃系数分段Xn-a(k)，并代入滤波器系数更新公式对滤波器系数进行更新。
当声学信道发生剧变时，特别是活跃系数在时域发生大幅度的位移时，仅对活跃系数分段Wa(k)进行自适应算法将无法跟踪声学信道的突变，为使算法具备快速跟踪能力，必须再循环更新一段非活跃系数。
具体的，迭代时从Wn-a(k)中循环提取一段滤波器分段并记为Wi(k)；对应地，从Xn-a(k)中循环提取一段输入信号分段并记为Xi(k)；所述滤波器系数更新模块提供的滤波器系数更新方程为：
e(k)＝d(k)-W(k)XT(k)；
Figure GDA0002225264410000075
Figure GDA0002225264410000076
下面举个具体的例子来说明上述更新方程：假设活跃系数分段位置已知为Wa(k)＝[W5(k)W6(k)]，则每次迭代均更新Wa(k)与循环更新一段Wi(k)，即每次迭代，依次更新以下滤波器分段：第一次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W1(k)；第二次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W2(k)；第三次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W3(k)；第四次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W4(k)；第五次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W7(k)；第六次迭代时，更新Wa(k)＝[W5(k)W6(k)]与W8(k)。
本发明提出的系数更新算法，只需要进行6次迭代即可完成对滤波器全部系数的更新。在6次迭代中活跃系数更新6次，非活跃系数只更新1次。可见，本发明提出的自适应滤波器系数更新方法中，活跃系数的更新频率为1，非活跃系数的更新频率降为1/(M-L)，通过降低非活跃系数更新频率以降低滤波器有效长度，进而最终达到加快收敛速度、降低计算复杂度的目的
相比于现有技术，本发明提出的自适应滤波器系数更新方法利用声学传播信道的稀疏性，在滤波器内通过移动窗积分法对信道进行延时估计以定位活跃系数位置，并通过精确辨识全部活跃系数完成声学传播信道的自适应辨识。进一步地，在每次迭代时更新滤波器的全部活跃系数，使得活跃系数更新频率得到提高，进而加快系统收敛速度；通过循环更新一段非活跃系数以降低非活跃系数更新频率，进而降低滤波器有效长度；同时，非活跃系数仍能获得一定的更新频率，从而确保了系统的跟踪能力。本发明提出的自适应滤波器系数更新方法可大幅降低滤波器有效长度，提高系统收敛速度，降低了算法复杂度。进一步地，本发明提出的滤波系数更新方法用一个滤波器完成延时估计与活跃系数辨识，可有效避免双滤波器结构信息冗余的缺陷、节省计算复杂度。与传统的NLMS算法相比，新算法可大幅降低滤波器有效长度，提高系统收敛速度，降低了计算复杂度。
本发明并不局限于上述实施方式，如果对本发明的各种改动或变形不脱离本发明的精神和范围，倘若这些改动和变形属于本发明的权利要求和等同技术范围之内，则本发明也意图包含这些改动和变形。

Claims (4)

1.一种自适应滤波的系数更新方法，其特征在于：包括如下步骤：
步骤1：将k时刻的滤波器W(k)的长度等分为M段；其中M为大于1的整数；
步骤2：求取每段滤波器的系数的二次范数；
步骤3：将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器系数的二次范数进行累加以获得最大累加值；
步骤4：以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余滤波器分段则为非活跃系数分段Wn-a(k)；
步骤4A：将k时刻的输入信号X(k)等分成M段，且将输入信号亦对应地分为输入信号活跃系数分段Xa(k)和输入信号非活跃系数分段Xn-a(k)，并代入滤波器系数更新公式对滤波器系数进行更新；滤波器系数迭代更新的方程为：
e(k)＝d(k)-W(k)XT(k)；
Figure FDA0002225264400000011
Figure FDA0002225264400000012
其中，μ为收敛因子，取值范围为0<μ<1；X(k)为输入信号，
Figure FDA0002225264400000013
为输入信号序列X(k)的二次范数，d(k)为回声信号，e(k)为输出信号；Wi(k)为从Wn-a(k)中循环提取的一段滤波器分段；对应地，Xi(k)为从Xn-a(k)中循环提取的一段输入信号分段；
步骤5：对滤波器系数进行迭代更新，且每次迭代时更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn-a(k)，以获得下一时刻k+1的滤波器系数。
2.根据权利要求1所述的自适应滤波的系数更新方法，其特征在于：在步骤3中，获得每段滤波器系数的二次范数并组成新序列，将长度为L的矩形窗在该新序列上移动，每移动一个位置，都对窗口内各分段滤波器系数的二次范数进行累加，以获得最大累加值；其中L为大于1的整数。
3.一种自适应滤波的系数更新装置，其特征在于：包括滤波器，所述滤波器包括：
——滤波器分段模块，用于将k时刻的滤波器W(k)的长度等分为M段；
——二次范数获取模块，用于求取每段滤波器系数的二次范数；
——最大累加值获取模块，用于将一矩形窗在该M段滤波器上移动，对窗口内的各分段滤波器系数的二次范数进行累加以获得最大累加值；其中M为大于1的整数；
——滤波器活跃系数分段和非活跃系数分段获取模块，用于以最大累加值所在窗口内的各滤波器分段作为滤波器活跃系数分段Wa(k)，其余滤波器分段则为非活跃系数分段Wn-a(k)；
——滤波器系数更新模块，用于对滤波器系数进行迭代更新，且每次迭代时用于更新滤波器的全部活跃系数分段Wa(k)并循环更新一段非活跃系数分段Wn-a(k),以获得下一时刻k+1的滤波器系数；
还包括输入信号处理模块；所述输入信号处理模块包括输入信号分段模块、以及输入信号活跃系数分段和非活跃系数分段获取模块；所述输入信号分段模块用于将k时刻的输入信号X(k)等分成M段；所述输入信号活跃系数分段和非活跃系数分段获取模块用于将输入信号也分为输入信号活跃系数分段Xa(k)和输入信号非活跃系数分段Xn-a(k)，并代入滤波器系数更新公式对滤波器系数进行更新；
所述滤波器系数更新模块对滤波器系数迭代更新的方程为：
e(k)＝d(k)-W(k)XT(k)；
Figure FDA0002225264400000021
Figure FDA0002225264400000022
其中，μ为收敛因子，取值范围为0<μ<1；X(k)为输入信号，
Figure FDA0002225264400000023
为输入信号序列X(k)的二次范数，d(k)为回声信号，e(k)为输出信号；Wi(k)为从Xn-a(k)中循环提取的一段滤波器分段；对应地，Xi(k)为从Xn-a(k)中循环提取的一段输入信号分段。
4.根据权利要求3所述的自适应滤波的系数更新装置，其特征在于：所述最大累加值获取模块还进一步用于获得每段滤波器系数的二次范数并组成新序列，将长度为L的矩形窗在该新序列上移动，每移动一个位置，都对窗口内各分段滤波器系数的二次范数进行累加，以获得最大累加值；其中L为大于1的整数。
CN201610046280.9A 2016-01-22 2016-01-22 一种自适应滤波的系数更新方法及装置 Active CN105654959B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610046280.9A CN105654959B (zh) 2016-01-22 2016-01-22 一种自适应滤波的系数更新方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610046280.9A CN105654959B (zh) 2016-01-22 2016-01-22 一种自适应滤波的系数更新方法及装置

Publications (2)

Publication Number Publication Date
CN105654959A CN105654959A (zh) 2016-06-08
CN105654959B true CN105654959B (zh) 2020-03-06

Family

ID=56484294

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610046280.9A Active CN105654959B (zh) 2016-01-22 2016-01-22 一种自适应滤波的系数更新方法及装置

Country Status (1)

Country Link
CN (1) CN105654959B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107483029B (zh) * 2017-07-28 2021-12-07 广州多益网络股份有限公司 一种voip通讯中的自适应滤波器的长度调节方法及装置
CN109889187B (zh) * 2019-01-15 2023-06-02 宁波连鸿电子科技有限公司 基于自适应滤波器的信号处理方法、装置及电子设备
CN110211599B (zh) * 2019-06-03 2021-07-16 Oppo广东移动通信有限公司 应用唤醒方法、装置、存储介质及电子设备
CN110855269B (zh) * 2019-11-06 2023-05-19 韶关学院 一种自适应滤波的系数更新方法
CN115171635A (zh) * 2022-05-30 2022-10-11 华为技术有限公司 噪声控制的方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1200210A (zh) * 1995-10-18 1998-11-25 艾利森电话股份有限公司 一种自适应双滤波器回声消除方法
CN1810019A (zh) * 2003-06-27 2006-07-26 诺基亚有限公司 统计自适应滤波器控制器
CN101461257A (zh) * 2006-05-30 2009-06-17 微软公司 自适应声学回声消除
CN102132491A (zh) * 2008-08-25 2011-07-20 杜比实验室特许公司 用于通过预白化确定通过lms算法调整的自适应滤波器的更新滤波系数的方法
EP2477188A1 (en) * 2011-01-18 2012-07-18 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Encoding and decoding of slot positions of events in an audio signal frame
CN102664006A (zh) * 2012-04-14 2012-09-12 中国人民解放军国防科学技术大学 基于时频域分析的异常人声检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4509126B2 (ja) * 2007-01-24 2010-07-21 沖電気工業株式会社 エコーキャンセラ及びエコーキャンセル方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1200210A (zh) * 1995-10-18 1998-11-25 艾利森电话股份有限公司 一种自适应双滤波器回声消除方法
CN1810019A (zh) * 2003-06-27 2006-07-26 诺基亚有限公司 统计自适应滤波器控制器
CN101461257A (zh) * 2006-05-30 2009-06-17 微软公司 自适应声学回声消除
CN102132491A (zh) * 2008-08-25 2011-07-20 杜比实验室特许公司 用于通过预白化确定通过lms算法调整的自适应滤波器的更新滤波系数的方法
EP2477188A1 (en) * 2011-01-18 2012-07-18 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Encoding and decoding of slot positions of events in an audio signal frame
CN102664006A (zh) * 2012-04-14 2012-09-12 中国人民解放军国防科学技术大学 基于时频域分析的异常人声检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Averaging Analysis of Partial-Update Adaptive Filters;Kutluyl Doğançay;《TENCON 2009 - 2009 IEEE Region 10 Conference》;20100122;全文 *
用于回声消除系统的自适应延时估计算法研究;陈国志;《科学技术与工程》;20150131;全文 *
面向实时通信系统的自适应回声消除算法研究;文昊翔;《中国博士学位论文全文数据库信息科技辑》;20140715;第4-7、91-98页 *
面向实时通信系统的自适应回声消除算法研究;文昊翔;《中国博士学位论文全文数据库信息科技辑》;20140715;第83-90页 *

Also Published As

Publication number Publication date
CN105654959A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
CN105654959B (zh) 一种自适应滤波的系数更新方法及装置
JP7018130B2 (ja) 遅延時間推定に基づくエコー除去方法及び装置
US11056130B2 (en) Speech enhancement method and apparatus, device and storage medium
CN103680515B (zh) 采用系数重用的比例自适应滤波器系数向量更新方法
WO2017202292A1 (zh) 一种回声时延跟踪方法及装置
CN103561185B (zh) 一种稀疏路径的回声消除方法
CN104050971A (zh) 声学回声减轻装置和方法、音频处理装置和语音通信终端
CN112037809A (zh) 基于多特征流结构深度神经网络的残留回声抑制方法
CN110211602B (zh) 智能语音增强通信方法及装置
CN103871420B (zh) 麦克风阵列的信号处理方法及装置
CN108154885A (zh) 一种使用qr-rls算法对多通道语音信号去混响方法
CN104994249B (zh) 声回波消除方法和装置
CN110138990A (zh) 一种消除移动设备VoIP电话回声的方法
CN109688284B (zh) 一种回音延时检测方法
CN105070296A (zh) 一种活性因子集员成比例子带自适应回声消除方法
CN106875938A (zh) 一种改进的非线性自适应语音端点检测方法
CN103369162A (zh) 一种低复杂度的电话回声自适应消除方法
CN104410761B (zh) 一种仿射投影符号子带凸组合自适应回声消除方法
Yang et al. A noise reduction method based on LMS adaptive filter of audio signals
CN107333017B (zh) 一种重用权系数指数型误差自适应回声消除方法
CN110491405A (zh) 基于协同非线性自适应滤波的麦克风阵列语音增强方法
CN110708651B (zh) 一种基于分段陷波的助听器啸叫检测与抑制方法及装置
Alam et al. Performance comparison of STFT, WT, LMS and RLS adaptive algorithms in denoising of speech signal
CN105070295A (zh) 一种活性因子成比例子带的自适应方法应用于回声消除
JP7495684B2 (ja) エコー除去方法、装置および電子機器

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant