CN106409307B - 一种具有选择进化仿射投影阶数的仿射投影方法 - Google Patents

一种具有选择进化仿射投影阶数的仿射投影方法 Download PDF

Info

Publication number
CN106409307B
CN106409307B CN201610858355.3A CN201610858355A CN106409307B CN 106409307 B CN106409307 B CN 106409307B CN 201610858355 A CN201610858355 A CN 201610858355A CN 106409307 B CN106409307 B CN 106409307B
Authority
CN
China
Prior art keywords
affine projection
signal
current time
projection order
affine
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
CN201610858355.3A
Other languages
English (en)
Other versions
CN106409307A (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN201610858355.3A priority Critical patent/CN106409307B/zh
Publication of CN106409307A publication Critical patent/CN106409307A/zh
Application granted granted Critical
Publication of CN106409307B publication Critical patent/CN106409307B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • 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
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L2021/02082—Noise filtering the noise being echo, reverberation of the speech

Landscapes

  • Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Quality & Reliability (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)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

一种具有选择进化仿射投影阶数的仿射投影方法，包括：初始化滤波器各参数以及定义阈值；更新输入信号与麦克风信号；计算先验误差信号；确定仿射投影阶数；更新滤波器系数。本发明的一种具有选择进化仿射投影阶数的仿射投影方法，用于声学回声消除器，能够使声学回声消除器在不影响仿射投影方法性能的前提下，通过自适应地选择仿射投影阶数来有效地降低计算复杂度，使仿射投影算法更易于实际应用。

Description

一种具有选择进化仿射投影阶数的仿射投影方法
技术领域
本发明涉及一种投影方法。特别是涉及一种针对于回声消除的具有选择进化仿射投影阶数的仿射投影方法。
背景技术
声学回声主要是指从扬声器发出的声音被麦克风重新拾取后再次发送到远端形成的回声。声学回声消除器(Acoustic Echo Cancellation,AEC)广泛应用于各种嵌入式设备和各种VoIP应用中，包括各种电信网络设备和终端设备、丰富的软件视频会议和VoIP软件电话等。
麦克风在时刻n接收的信号为：
d(n)＝xT(n)h(n)+v(n) 1.1
式中：d(n)代表麦克风信号,输入信号向量x(n)＝[x(n),x(n-1),...,x(n-N+1)]T，T代表矩阵的转置，N为回声路径长度，h(n)＝[h0(n),h1(n),...,hN(n)]T是系统的冲击响应，v(n)代表近端信号。回声消除的目的就是利用估计出回声路径将回声从麦克风信号中消除。
代表估计出来的回声信号，e(n)代表回声消除后得到的误差信号，其中，线性回声消除器主要依靠自适应算法来估计回声路径，由于仿射投影算法(Affine ProjectionAlgorithm,APA)能够在收敛速度和计算复杂度之间取得较好折衷而被广泛应用，对于一个最大仿射投影阶数为Kmax的APA，其更新方程是
其中μ为自适应滤波器的步长，Kn为当前时刻的仿射投影阶数，为输入信号矩阵，为单位矩阵，δ为一常数称为正则化因子，为前验误差信号向量，计算方法是
其中d(n)＝[d(n),d(n-1),...,d(n-Kn+1)]T。
发明内容
本发明所要解决的技术问题是，提供一种在不影响仿射投影方法性能的前提下，通过自适应地选择仿射投影阶数来降低计算复杂度的具有选择进化仿射投影阶数的仿射投影方法。
本发明所采用的技术方案是：一种具有选择进化仿射投影阶数的仿射投影方法，包括如下步骤：
1)初始化滤波器各参数以及定义阈值ε(K)：
对滤波器的系数步长μ、正则化因子δ和最大仿射投影阶数Kmax进行初始化，阈值其中为近端噪声功率；
2)更新输入信号与麦克风信号：
包括获取当前时刻麦克风信号d(n)与输入信号x(n)，并分别更新至麦克风信号向量d(n)与输入信号向量x(n)中，所述d(n)＝[d(n),d(n-1),...,d(n-Kmax+1)]T，所述输入信号向量x(n)＝[x(n),x(n-1),...,x(n-N+1)]T；
3)计算先验误差信号：
将n-1时刻估计的滤波器系数代入下式
得到当前时刻回声信号估计值输入信号向量x(n)＝[x(n),x(n-1),...,x(n-N+1)]T
将当前时刻回声信号估计值和当前时刻麦克风信号d(n)代入下式
得到当前时刻先验误差信号e(n)；
4)确定仿射投影阶数：
令K＝Kn-1，若当前时刻先验误差信号e(n)的平方小于阈值ε(K)，则将前一时刻仿射投影阶数Kn-1的值减1后设定为当前时刻仿射投影阶数Kn，由于当前时刻仿射投影阶数Kn不能小于1，所以当前时刻仿射投影阶数Kn＝max{Kn-1,1}；若当前时刻先验误差信号e(n)的平方大于阈值ε(K+1)，则不断地将K的值加1，直至当前时刻先验误差信号e(n)的平方小于阈值ε(K+1)或者K大于最大仿射投影阶数Kmax为止，此时当前仿射投影阶数Kn＝min{K+1,Kmax}；
5)更新滤波器系数：
根据当前投影阶数，选择对应长度的输入信号矩阵先验误差信号矩阵然后利用所述矩阵进行滤波器的更新。
步骤5)所述的输入信号矩阵表示如下：
步骤5)所述的先验误差信号矩阵表示如下：
计算方法是：
其中
本发明的一种具有选择进化仿射投影阶数的仿射投影方法，用于声学回声消除器，能够使声学回声消除器在不影响仿射投影方法性能的前提下，通过自适应地选择仿射投影阶数来有效地降低计算复杂度，使仿射投影算法更易于实际应用。
具体实施方式
下面结合实施例对本发明的一种具有选择进化仿射投影阶数的仿射投影方法做出详细说明。
本发明的一种具有选择进化仿射投影阶数的仿射投影方法，在不影响仿射投影方法性能的前提下，通过自适应地选择仿射投影阶数来降低计算复杂度(收敛速度和稳态误差)。
本发明的一种具有选择进化仿射投影阶数的仿射投影方法，包括如下步骤：
1)初始化滤波器各参数以及定义阈值ε(K)：
对滤波器的系数步长μ、正则化因子δ和最大仿射投影阶数Kmax进行初始化，阈值其中为近端噪声功率；
本发明实施例中设定，滤波器系数步长μ＝0.5、正则化因子和最大仿射投影阶数Kmax＝8。其中滤波器长度N＝1024，λ＝1-1/(6N)，近端信号中加入信噪比20dB的独立平稳高斯白噪声。
2)更新输入信号与麦克风信号：
包括获取当前时刻麦克风信号d(n)与输入信号x(n)，并分别更新至麦克风信号向量d(n)与输入信号向量x(n)中，所述d(n)＝[d(n),d(n-1),...,d(n-Kmax+1)]T，所述输入信号向量x＝[x(n),x(n-1),...,x(n-N+1)]T；
3)计算先验误差信号：
将n-1时刻估计的滤波器系数代入下式
得到当前时刻回声信号估计值输入信号向量x(n)＝[x(n),x(n-1),...,x(n-N+1)]T
将当前时刻回声信号估计值和当前时刻麦克风信号d(n)代入下式
得到当前时刻先验误差信号e(n)；
4)确定仿射投影阶数：
令K＝Kn-1，若当前时刻先验误差信号e(n)的平方小于阈值ε(K)，则将前一时刻仿射投影阶数Kn-1的值减1后设定为当前时刻仿射投影阶数Kn，由于当前时刻仿射投影阶数Kn不能小于1，所以当前时刻仿射投影阶数Kn＝max{Kn-1,1}；若当前时刻先验误差信号e(n)的平方大于阈值ε(K+1)，则不断地将K的值加1，直至当前时刻先验误差信号e(n)的平方小于阈值ε(K+1)或者K大于最大仿射投影阶数Kmax为止，此时当前仿射投影阶数Kn＝min{K+1,Kmax}；
5)更新滤波器系数：
根据当前投影阶数，选择对应长度的输入信号矩阵先验误差信号矩阵然后利用所述矩阵进行滤波器的更新。
所述的输入信号矩阵表示如下：
所述的先验误差信号矩阵表示如下：
计算方法是：其中

Claims (3)

1.一种具有选择进化仿射投影阶数的仿射投影方法，其特征在于，包括如下步骤：
1)初始化滤波器各参数以及定义阈值ε(K)：
对滤波器的系数步长μ、正则化因子δ和最大仿射投影阶数Kmax进行初始化，阈值其中为近端噪声功率，K为仿射投影阶数；
2)更新输入信号与麦克风信号：
包括获取当前时刻麦克风信号d(n)与输入信号x(n)，并分别更新至麦克风信号向量d(n)与输入信号向量x(n)中，所述d(n)＝[d(n),d(n-1),...,d(n-Kmax+1)]T，所述输入信号向量x(n)＝[x(n),x(n-1),...,x(n-N+1)]T；
3)计算先验误差信号：
将n-1时刻估计的滤波器系数代入下式
得到当前时刻回声信号估计值输入信号向量x(n)＝[x(n),x(n-1),...,x(n-N+1)]T
将当前时刻回声信号估计值和当前时刻麦克风信号d(n)代入下式
得到当前时刻先验误差信号e(n)；
4)确定仿射投影阶数：
令K＝Kn-1，若当前时刻先验误差信号e(n)的平方小于阈值ε(K)，则将前一时刻仿射投影阶数Kn-1的值减1后设定为当前时刻仿射投影阶数Kn，由于当前时刻仿射投影阶数Kn不能小于1，所以当前时刻仿射投影阶数Kn＝max{Kn-1,1}；若当前时刻先验误差信号e(n)的平方大于阈值ε(K+1)，则不断地将K的值加1，直至当前时刻先验误差信号e(n)的平方小于阈值ε(K+1)或者K大于最大仿射投影阶数Kmax为止，此时当前仿射投影阶数Kn＝min{K+1,Kmax}；
5)更新滤波器系数：
根据当前投影阶数，选择对应长度的输入信号矩阵先验误差信号矩阵然后利用所述矩阵进行滤波器的更新。
2.根据权利要求1所述的一种具有选择进化仿射投影阶数的仿射投影方法，其特征在于，步骤5)所述的输入信号矩阵表示如下：
3.根据权利要求2所述的一种具有选择进化仿射投影阶数的仿射投影方法，其特征在于，步骤5)所述的先验误差信号矩阵表示如下：
计算方法是：
其中
CN201610858355.3A 2016-09-28 2016-09-28 一种具有选择进化仿射投影阶数的仿射投影方法 Active CN106409307B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610858355.3A CN106409307B (zh) 2016-09-28 2016-09-28 一种具有选择进化仿射投影阶数的仿射投影方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610858355.3A CN106409307B (zh) 2016-09-28 2016-09-28 一种具有选择进化仿射投影阶数的仿射投影方法

Publications (2)

Publication Number Publication Date
CN106409307A CN106409307A (zh) 2017-02-15
CN106409307B true CN106409307B (zh) 2019-09-27

Family

ID=57997853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610858355.3A Active CN106409307B (zh) 2016-09-28 2016-09-28 一种具有选择进化仿射投影阶数的仿射投影方法

Country Status (1)

Country Link
CN (1) CN106409307B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109151237B (zh) * 2018-08-23 2020-10-09 西南交通大学 基于零吸引的仿射投影自适应回声消除方法
CN109767779A (zh) * 2018-11-17 2019-05-17 沈阳工业大学 基于最小误差熵的成比例仿射投影方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002369222A (ja) * 2001-06-06 2002-12-20 Kddi Corp 勾配法を利用した三次元動画像の動き補償予測における動きベクトル検出方法
CN1859519A (zh) * 2005-11-19 2006-11-08 华为技术有限公司 一种自适应滤波器及回波抵消器
CN105407243A (zh) * 2015-10-26 2016-03-16 南京邮电大学 一种Android平台上使用改进仿射投影算法的回声消除VOIP系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6909782B2 (en) * 2000-09-08 2005-06-21 Intel Corporation Fast converging affine projection based echo canceller for sparse multi-path channels

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002369222A (ja) * 2001-06-06 2002-12-20 Kddi Corp 勾配法を利用した三次元動画像の動き補償予測における動きベクトル検出方法
CN1859519A (zh) * 2005-11-19 2006-11-08 华为技术有限公司 一种自适应滤波器及回波抵消器
CN105407243A (zh) * 2015-10-26 2016-03-16 南京邮电大学 一种Android平台上使用改进仿射投影算法的回声消除VOIP系统

Also Published As

Publication number Publication date
CN106409307A (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN105577961B (zh) 增益控制器的自动调谐
CN105825864B (zh) 基于过零率指标的双端说话检测与回声消除方法
CN106657507B (zh) 一种声学回声消除方法及装置
CN109754813A (zh) 基于快速收敛特性的变步长回声消除方法
CN105308942B (zh) 回波消除
CN106448695B (zh) 一种双端通话鲁棒的变阶变步长仿射投影方法
US10115411B1 (en) Methods for suppressing residual echo
CN111742541B (zh) 声学回波抵消方法、装置、存储介质
CN103534942B (zh) 处理音频信号
CN111145771A (zh) 语音信号处理方法、处理装置、终端及其存储介质
WO2021110175A1 (zh) 回声消除方法及装置
CN107026950B (zh) 一种频域自适应回声消除方法
CN102739886A (zh) 基于回声频谱估计和语音存在概率的立体声回声抵消方法
CN108877824B (zh) 一种跟踪性能高的组合步长回声消除方法
CN106409307B (zh) 一种具有选择进化仿射投影阶数的仿射投影方法
EP2101480A3 (en) Echo canceller and echo cancelling method
CN110191245A (zh) 一种基于时变参数的自适应回声消除方法
CN112929506B (zh) 音频信号的处理方法及装置，计算机存储介质及电子设备
US9729968B2 (en) Method and system for acoustic echo cancellation using cascaded kalman filtering
CN107071196B (zh) 一种自适应回声消除方法
CN105721729A (zh) 基于块稀疏成比例重用权系数仿射投影的回声消除方法
CN105788605B (zh) 一种基于正交变换的仿射投影符号回声消除方法
CN109151237B (zh) 基于零吸引的仿射投影自适应回声消除方法
CN106157965B (zh) 一种基于权向量重用的零范数集员仿射投影自适应回声消除方法
Diniz et al. Set-membership affine projection algorithm for echo cancellation

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