CN107483029B - 一种voip通讯中的自适应滤波器的长度调节方法及装置 - Google Patents
一种voip通讯中的自适应滤波器的长度调节方法及装置 Download PDFInfo
- Publication number
- CN107483029B CN107483029B CN201710627615.0A CN201710627615A CN107483029B CN 107483029 B CN107483029 B CN 107483029B CN 201710627615 A CN201710627615 A CN 201710627615A CN 107483029 B CN107483029 B CN 107483029B
- Authority
- CN
- China
- Prior art keywords
- filter
- length
- coefficient
- module
- max
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/0009—Time-delay networks
- H03H17/0018—Realizing a fractional delay
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal 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
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0043—Adaptive algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M9/00—Arrangements for interconnection not involving centralised switching
- H04M9/08—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
- H04M9/085—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic using digital techniques
Abstract
本发明提供一种自适应滤波器的长度调节方法,所述滤波器设有参考信号延时模块和可变长度模块;所述滤波器的长度调节方法包括以下步骤:S1:获取滤波器系数W(k);S2:判断||W(k)||是否小于En,若是,则结束;若否,则执行步骤S3;其中,所述En为能量下限值;S3:对滤波器系数进行归一化处理,获得归一化系数WN(k);S4:计算加窗的滤波器归一化系数WNW(k);S5:计算最大系数WNW,i max(k),并标记为imax;S6:计算点p和点b。S7:如果1≤b≤imax‑kb,增加远端滤波器延时,并且滤波器长度变为L(k)=L(k‑1)‑b;如果imax+kp≤p<L,滤波器长度变为L(k)=p;kb和kp分别表示表示不能抛弃长度的最小值和最大值,是预设值,本发明还提供了一种用于实现上述方法的长度调节装置。
Description
技术领域
本发明涉及一种voip通讯领域,特别是一种voip通讯中的自适应滤波器的长度调节方法及装置。
背景技术
在voip通信中,回声消除是一个很重要的部分,回声消除模块一般至少包括一个线性滤波器,用于估计回声并消除回声。线性滤波器的长度是一个重要的参数,如果滤波器长了,能够处理的延时相对比较长,但是要付出更多的运算量,更长的收敛时间;如果滤波器过短,只能处理短延时的回声,可能没有办法处理回声的拖尾效应。如果系统是非因果的,长的滤波器反而会导致误差增大,导致多重回声等反效果。
请同时参阅图1,其为现有技术的自适应滤波器基本流程示意图。
这里用NLMS算法为例说明滤波器的更新,其中NLMS算法的滤波器更新公式如下:
其中k=0,1,...表示时刻;W(k)=[w1(k) w2(k) … wL(k)]表示滤波器系数,L是滤波器长度;X(k)=[x1(k) x2(k) … xL(k)]是参考信号,e(k)=d(k)-W*(k)X(k)是误差信号d(k),W*(k)分别是期望信号和最优滤波器系数;μ是更新步长。
一般的算法使用误差的能量e2(k)值及其函数变换来调节滤波器的长度。
然而,现有技术仍然存在以下的缺点和不足:
通过收敛后的误差来判断滤波器长度的方法特别容易受到信号非线性畸变和信噪比降低的影响。因为这些情况滤波器无法收敛到合适的值,同时因此无法降低滤波器的长度。而另外如果系统因为延时不准确,造成了非因果的情况,滤波器系数越长,误差反而会越大。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供了一种voip通讯中的自适应滤波器长度调节的方法。
本发明通过以下的方案实现:一种voip通讯中的自适应滤波器的长度调节方法,所述滤波器设有参考信号延时模块和可变长度模块;所述滤波器的长度调节方法包括以下步骤:
S1:滤波器获取滤波器系数W(k);
S2:滤波器判断||W(k)||是否小于En,若是,则结束长度调节;若否,则滤波器执行步骤S3;其中,所述En为能量下限值;
S3:滤波器对滤波器系数进行归一化处理,获得归一化系数WN(k);
S4:滤波器计算加窗的滤波器归一化系数WNW(k),具体通过使用一个长度为2Lw的窗口KW来计算加窗后的归一化系数WNW(k),KW是设定的窗系数;
S5:滤波器计算最大系数WNW,i max(k),并标记为i max;
λb,λp表示能量聚集的度,是预设值;
S7:如果1≤b≤i max-kb,通过所述参考信号延时模块增加远端滤波器延时,并且通过所述可变长度模块将滤波器长度变为L(k)=L(k-1)-b;
如果i max+kp≤p<L,通过所述可变长度模块将滤波器长度变为L(k)=p;
kb和kp分别表示不能抛弃长度的最小值和最大值,是预设值。
作为本发明的进一步改进,所述步骤S5中的最大系数WNW,i max(k)=max(WNW,n(k))。
本发明还提供了一种voip通讯中的自适应滤波器的长度调节装置,所述滤波器中包括参考信号延时模块和可变长度模块,该滤波器的长度调节装置包括:
滤波器系数获取模块,用于获取滤波器系数W(k);
判断模块,用于判断||W(k)||是否小于En,若是,这结束操作,若否,则继续执行;所述En为能量下限值;
归一化模块,用于对滤波器系数进行归一化处理,获得归一化系数WN(k);
加窗模块,用于计算加窗的滤波器归一化系数WNW(k),具体通过使用一个长度为2Lw的窗口KW来计算加窗后的归一化系数WNW(k),KW是设定的窗系数;
最大系数计算模块,用于计算最大系数WNW,i max(k),并标记为i max;
滤波器长度调整模块,用于当1≤b≤i max-kb,通过所述参考信号延时模块增加远端滤波器延时,并且通过所述参考信号延时模块将滤波器长度变为L(k)=L(k-1)-b;当imax+kp≤p<L,通过所述参考信号延时模块将滤波器长度变为L(k)=p;所述kb和kp分别表示不能抛弃长度的最小值和最大值,是预设值。
作为本发明的进一步改进,所述最大系数计算模块中的最大系数WNW,i max(k)=max(WNW,n(k))。
综上,相比于现有技术,本发明在传统的滤波器的结构上进行了调整,增加了延时模块和可变长度。然后,通过本发明的长度调节方法,可以实现动态调节滤波器长度,降低运算量,从而消除了现有技术的缺陷。
为了更好地理解和实施,下面结合附图详细说明本发明。
附图说明
图1是现有技术的自适应滤波器的基本流程框图。
图2是本发明的自适应滤波器的基本流程框图。
图3是本发明的自适应滤波器的长度调节方法的步骤流程图。
图4是本发明的自适应滤波器的长度调节装置的模块框图。
具体实施方式
以下结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
为了解决现有技术中的缺陷,本发明通过分析滤波器系数的分布关系,来判断滤波器的收敛性。本发明的方法只能用在最优滤波器W*(k)符合以下特征的情况使用:
a)系数具有明显单峰值;
b)峰值前的系数比较小;
c)峰值后的系数呈现出衰减状态。
这三个条件符合声学回声消除应用的场景。第一条,表示有明显的主回声回路;第二条,表示主回声回路之前的回声小;第三条,表示回声回路有多条路径,并且为衰减状况。
请参阅图2,其为本发明的自适应滤波器的基本流程框图。本发明在传统的滤波器中加入参考信号延时模块和可变长度模块,提出的滤波器的结构如图2所示,在参考信号x(k)之后加入了延时模块,同时增加了可变长度函数L(k)。
具体的请同时参阅图3,其为本发明的自适应滤波器的长度调节方法的步骤流程图。基于上述的滤波器的结构,本发明提供了一种自适应滤波器的长度调节方法,具体包括以下步骤:
S1:获取滤波器系数W(k)。
S2:判断||W(k)||是否小于En,若是,则结束;若否,则执行步骤S3;其中,所述En为能量下限值。
通过WN(k)计算滤波器分布。使用WN(k)的时候需要保证||W(k)||≥En,其中En表示能量下限,小于该下限表示能量比较小,不进入后面的计算流程。
S5:计算最大系数WNW,i max(k)=max(WNW,n(k)),并标记为imax。
λb,λp表示能量聚集的度,是预设值。
例如当前滤波器长度L=256,计算i max=64,λb=0.1,λp=0.8,通过上面计算方法获得b=50,p=128。那么表示前面50个点只有占权重的10%,第64个点到第128个点占有权重的80%。
S7:如果1≤b≤i max-kb,增加远端滤波器延时,并且滤波器长度变为L(k)=L(k-1)-b。
如果i max+kp≤p<L,滤波器长度变为L(k)=p。
kb和kp分别表示表示不能抛弃长度的最小值和最大值,是预设值,保护i max前后的值不被截断。
另外,请同时参阅图4,其为本发明的自适应滤波器的长度调节装置的模块框图。为了实现上述的方法,本发明还提供了一种自适应滤波器的长度调节装置该滤波器的长度调节装置包括:滤波器系数获取模块1、判断模块2、归一化模块3、加窗模块4、最大系数计算模块5、p点计算模块6、b点计算模块7和滤波器长度调整模块8。
滤波器系数获取模块,用于获取滤波器系数W(k);
判断模块,用于判断||W(k)||是否小于En,若是,这结束操作,若否,则继续执行;所述En为能量下限值;
最大系数计算模块,用于计算最大系数WNW,i max(k)=max(WNW,n(k)),并标记为imax。
滤波器长度调整模块,用于当1≤b≤i max-kb,增加远端滤波器延时,并且滤波器长度变为L(k)=L(k-1)-b;当i max+kp≤p<L,滤波器长度变为L(k)=p;所述kb和kp分别表示表示不能抛弃长度的最小值和最大值,是预设值。
本发明的长度调节装置的原理与上述长度调节方法的原理相同,故不赘述。
综上,相比于现有技术,本发明在传统的滤波器的结构上进行了调整,增加了延时模块和可变长度。然后,通过本发明的长度调节方法,可以实现动态调节滤波器长度,降低运算量,从而消除了现有技术的缺陷。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (8)
1.一种voip通讯中的自适应滤波器的长度调节方法,其特征在于:所述滤波器设有参考信号延时模块和可变长度模块;该滤波器的长度调节方法包括以下步骤:
S1:滤波器获取滤波器系数W(k);
S2:滤波器判断||W(k)||是否小于En,若是,则结束长度调节;若否,则滤波器执行步骤S3;其中,所述En为能量下限值;
S3:滤波器对滤波器系数进行归一化处理,获得归一化系数WN(k);
S4:滤波器计算加窗的滤波器归一化系数WNW(k),具体通过使用一个长度为2Lw的窗口KW来计算加窗后的归一化系数WNW(k),KW是设定的窗系数;
S5:滤波器计算最大系数WNW,imax(k),并标记为i max;
λb,λp表示能量聚集的度,是预设值;
S7:如果1≤b≤i max-kb,通过所述参考信号延时模块增加远端滤波器延时,并且通过所述可变长度模块将滤波器长度变为L(k)=L(k-1)-b;
如果i max+kp≤p<L,通过所述可变长度模块将滤波器长度变为L(k)=p;
kb和kp分别表示不能抛弃长度的最小值和最大值,是预设值。
4.根据权利要求1所述的voip通讯中的自适应滤波器的长度调节方法,其特征在于:所述步骤S5中的最大系数WNW,imax(k)=max(WNW,n(k))。
5.一种voip通讯中的自适应滤波器的长度调节装置,其特征在于:所述滤波器中包括参考信号延时模块和可变长度模块,该滤波器的长度调节装置包括:
滤波器系数获取模块,用于获取滤波器系数W(k);
判断模块,用于判断||W(k)||是否小于En,若是,这结束操作,若否,则继续执行;所述En为能量下限值;
归一化模块,用于对滤波器系数进行归一化处理,获得归一化系数WN(k);
加窗模块,用于计算加窗的滤波器归一化系数WNW(k),具体通过使用一个长度为2Lw的窗口KW来计算加窗后的归一化系数WNW(k),KW是设定的窗系数;
最大系数计算模块,用于计算最大系数WNW,imax(k),并标记为i max;
滤波器长度调整模块,用于当1≤b≤i max-kb,通过所述参考信号延时模块增加远端滤波器延时,并且通过所述可变长度模块将滤波器长度变为L(k)=L(k-1)-b;当i max+kp≤p<L,通过所述可变长度模块将滤波器长度变为L(k)=p;所述kb和kp分别表示不能抛弃长度的最小值和最大值,是预设值。
8.根据权利要求5所述的voip通讯中的自适应滤波器的长度调节装置,其特征在于:所述最大系数计算模块中的最大系数WNW,imax(k)=max(WNW,n(k))。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710627615.0A CN107483029B (zh) | 2017-07-28 | 2017-07-28 | 一种voip通讯中的自适应滤波器的长度调节方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710627615.0A CN107483029B (zh) | 2017-07-28 | 2017-07-28 | 一种voip通讯中的自适应滤波器的长度调节方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107483029A CN107483029A (zh) | 2017-12-15 |
CN107483029B true CN107483029B (zh) | 2021-12-07 |
Family
ID=60597202
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710627615.0A Active CN107483029B (zh) | 2017-07-28 | 2017-07-28 | 一种voip通讯中的自适应滤波器的长度调节方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107483029B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111277718B (zh) * | 2020-01-21 | 2021-10-08 | 上海推乐信息技术服务有限公司 | 一种回声消除系统及其方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1691716A (zh) * | 2004-04-23 | 2005-11-02 | 北京三星通信技术研究有限公司 | 回声消除装置 |
CN102195616A (zh) * | 2010-03-02 | 2011-09-21 | 哈曼贝克自动系统股份有限公司 | 高效子带自适应有限脉冲响应滤波 |
CN102377453A (zh) * | 2010-08-06 | 2012-03-14 | 联芯科技有限公司 | 一种控制自适应滤波器更新的方法、装置及回声消除器 |
CN102571650A (zh) * | 2011-12-20 | 2012-07-11 | 东南大学 | 一种应用于3gpp lte系统的自适应信道估计方法 |
CN102938254A (zh) * | 2012-10-24 | 2013-02-20 | 中国科学技术大学 | 一种语音信号增强系统和方法 |
CN103714810A (zh) * | 2013-12-09 | 2014-04-09 | 西北核技术研究所 | 基于Gammatone滤波器组的车型特征提取方法 |
KR101390551B1 (ko) * | 2012-09-24 | 2014-04-30 | 충북대학교 산학협력단 | 저 지연 변형된 이산 코사인 변환 방법 |
EP2860730A1 (en) * | 2013-10-10 | 2015-04-15 | Nokia Corporation | Speech processing |
CN104520925A (zh) * | 2012-08-01 | 2015-04-15 | 杜比实验室特许公司 | 噪声降低增益的百分位滤波 |
CN105654959A (zh) * | 2016-01-22 | 2016-06-08 | 韶关学院 | 一种自适应滤波的系数更新方法及装置 |
CN105825864A (zh) * | 2016-05-19 | 2016-08-03 | 南京奇音石信息技术有限公司 | 基于过零率指标的双端说话检测与回声消除方法 |
US9691378B1 (en) * | 2015-11-05 | 2017-06-27 | Amazon Technologies, Inc. | Methods and devices for selectively ignoring captured audio data |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10306389B2 (en) * | 2013-03-13 | 2019-05-28 | Kopin Corporation | Head wearable acoustic system with noise canceling microphone geometry apparatuses and methods |
CN105849804A (zh) * | 2013-12-23 | 2016-08-10 | 美国亚德诺半导体公司 | 过滤噪声的计算高效方法 |
-
2017
- 2017-07-28 CN CN201710627615.0A patent/CN107483029B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1691716A (zh) * | 2004-04-23 | 2005-11-02 | 北京三星通信技术研究有限公司 | 回声消除装置 |
CN102195616A (zh) * | 2010-03-02 | 2011-09-21 | 哈曼贝克自动系统股份有限公司 | 高效子带自适应有限脉冲响应滤波 |
CN102377453A (zh) * | 2010-08-06 | 2012-03-14 | 联芯科技有限公司 | 一种控制自适应滤波器更新的方法、装置及回声消除器 |
CN102571650A (zh) * | 2011-12-20 | 2012-07-11 | 东南大学 | 一种应用于3gpp lte系统的自适应信道估计方法 |
CN104520925A (zh) * | 2012-08-01 | 2015-04-15 | 杜比实验室特许公司 | 噪声降低增益的百分位滤波 |
KR101390551B1 (ko) * | 2012-09-24 | 2014-04-30 | 충북대학교 산학협력단 | 저 지연 변형된 이산 코사인 변환 방법 |
CN102938254A (zh) * | 2012-10-24 | 2013-02-20 | 中国科学技术大学 | 一种语音信号增强系统和方法 |
EP2860730A1 (en) * | 2013-10-10 | 2015-04-15 | Nokia Corporation | Speech processing |
CN103714810A (zh) * | 2013-12-09 | 2014-04-09 | 西北核技术研究所 | 基于Gammatone滤波器组的车型特征提取方法 |
US9691378B1 (en) * | 2015-11-05 | 2017-06-27 | Amazon Technologies, Inc. | Methods and devices for selectively ignoring captured audio data |
CN105654959A (zh) * | 2016-01-22 | 2016-06-08 | 韶关学院 | 一种自适应滤波的系数更新方法及装置 |
CN105825864A (zh) * | 2016-05-19 | 2016-08-03 | 南京奇音石信息技术有限公司 | 基于过零率指标的双端说话检测与回声消除方法 |
Non-Patent Citations (3)
Title |
---|
Energy behavior in time-varying fractional delay filters for physical modeling synthesis of musical instruments;J. Pakarinen 等;《Proceedings. (ICASSP "05). IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005.》;20050323;III-1--III-4 * |
单滤波器延时估计-局部迭代算法;文昊翔 等;《信号处理》;20160325;第32卷(第3期);321-326 * |
用于回声消除系统的自适应延时估计算法研究;陈国志 等;《科学技术与工程》;20150128;第15卷(第3期);244-249 * |
Also Published As
Publication number | Publication date |
---|---|
CN107483029A (zh) | 2017-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU752841B2 (en) | Noise cancelling method and noise cancelling unit | |
US8184828B2 (en) | Background noise estimation utilizing time domain and spectral domain smoothing filtering | |
KR20190085927A (ko) | 적응성 빔포밍 | |
CN109273019B (zh) | 用于回声抑制的双重通话检测的方法及回声抑制 | |
WO2012049986A1 (ja) | 信号処理装置、信号処理方法、並びに信号処理プログラム | |
US9361874B2 (en) | Varying adaptive filter step size in acoustic echo cancellation | |
US9443530B2 (en) | Acoustic echo cancellation | |
KR20050074503A (ko) | 잔여 반향 필터를 사용하여 음향 반향 소거 시스템을향상시키기 위한 방법 | |
JPH08241086A (ja) | 雑音消去装置 | |
US9036816B1 (en) | Frequency domain acoustic echo cancellation using filters and variable step-size updates | |
CN110191245B (zh) | 一种基于时变参数的自适应回声消除方法 | |
US9172791B1 (en) | Noise estimation algorithm for non-stationary environments | |
CN107483029B (zh) | 一种voip通讯中的自适应滤波器的长度调节方法及装置 | |
US7177416B1 (en) | Channel control and post filter for acoustic echo cancellation | |
CN112151060B (zh) | 单通道语音增强方法及装置、存储介质、终端 | |
CN109326301A (zh) | 自适应后滤波 | |
CN106782595B (zh) | 一种降低语音泄露的鲁棒阻塞矩阵方法 | |
WO2021016000A2 (en) | Frequency domain adaptation with dynamic step size adjustment based on analysis of statistic of adaptive filter coefficient movement | |
CN108462481B (zh) | 一种μ律函数下基于参数调整的比例LMP滤波方法 | |
US10582315B2 (en) | Feedback canceller and hearing aid | |
CN115834778A (zh) | 一种回声消除方法、装置、电子设备及存储介质 | |
JP5105912B2 (ja) | 音声明瞭度改善装置およびその騒音レベル推定方法 | |
JP4068182B2 (ja) | 適応フィルタ | |
CN109308907B (zh) | 单信道降噪 | |
CN112039498B (zh) | 基于多态可变步长最小均方的自适应信号处理方法及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |