CN103956992A - 一种基于多步梯度下降的自适应信号处理方法 - Google Patents

一种基于多步梯度下降的自适应信号处理方法 Download PDF

Info

Publication number
CN103956992A
CN103956992A CN201410116777.4A CN201410116777A CN103956992A CN 103956992 A CN103956992 A CN 103956992A CN 201410116777 A CN201410116777 A CN 201410116777A CN 103956992 A CN103956992 A CN 103956992A
Authority
CN
China
Prior art keywords
signal processing
gradient descent
mrow
adaptive signal
size
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.)
Granted
Application number
CN201410116777.4A
Other languages
English (en)
Other versions
CN103956992B (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.)
SHANGHAI FUKONG HUALONG MICROSYSTEM TECHNOLOGY Co Ltd
Fudan University
Original Assignee
SHANGHAI FUKONG HUALONG MICROSYSTEM TECHNOLOGY Co Ltd
Fudan 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 SHANGHAI FUKONG HUALONG MICROSYSTEM TECHNOLOGY Co Ltd, Fudan University filed Critical SHANGHAI FUKONG HUALONG MICROSYSTEM TECHNOLOGY Co Ltd
Priority to CN201410116777.4A priority Critical patent/CN103956992B/zh
Publication of CN103956992A publication Critical patent/CN103956992A/zh
Application granted granted Critical
Publication of CN103956992B publication Critical patent/CN103956992B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Filters That Use Time-Delay Elements (AREA)

Abstract

本发明公开了一种基于多步梯度下降的自适应信号处理方法,属于自适应信号处理技术领域,该方法将多步梯度下降方法和自适应信号处理的NLMS技术相结合,提出了多步梯度下降的NLMS技术。具体利用多步梯度下降方法来加快基于标准梯度下降方法的NLMS技术的收敛速度,以此来提高信号处理的跟踪能力。实验结果表明,该技术的性能十分理想,能够有效解决现有技术所存在的问题。

Description

一种基于多步梯度下降的自适应信号处理方法
技术领域
本发明属于自适应信号处理技术领域,具体为一种基于多步梯度下降的变步长的归一化最小均方(Normalized Least Mean Squares,NLMS)技术。
背景技术
在过去的30年中,信号处理技术获得了巨大的发展。由于数字电路设计的进步,促使人们对数字信号处理领域越来越感兴趣。数字信号处理系统之所以具有吸引力,是因为它具有可靠性高、精度高、体积小及灵活性强等特点。随着数字超大规模集成(VLSI,very large scale integration)技术的发展,自适应信号处理技术在许多领域得到了广泛的应用,例如自适应信号处理在噪声和回声抵消、谱线增强、信道均衡、系统辨识和延迟时间估计等应用中起着非常重要的作用。针对这些应用,在过去的几十年中,人们提出了许多自适应滤波的技术。其中,基于标准梯度下降的最小均方(least mean squares,LMS)技术是经典的方法。由它发展而来的归一化最小均方(normalized least mean squares,NLMS)自适应滤波技术,由于其具有较好的鲁棒性和易于使用,从而使其成为最常用的自适应滤波技术。
NLMS技术成功应用的关键在于其步长的选取,这是因为人们需要它来折中信号处理过程中的收敛速度和稳态失调。当步长取的越大,处理过程中的收敛性就越快,但最终的稳态失调也越大;步长取的越小,最终的稳态失调越小,但处理的收敛速度比较慢。可是,在应用中如何选取可靠且最优的步长一直都是一个未能很好解决。在众多不同的步长选取策略中,大部分都以变步长为特点。通过分析不难发现:变步长技术虽然解决了信号处理过程中收敛速度的问题,但由于时变的步长不能最终收敛到零,在大多数情况下其是一个未知的数值,因此人们一般并不能知道最终稳态失调的大小。
发明内容
针对现有自适应信号处理技术所存在的问题,本发明的目的在于提供一种基于多步梯度下降的自适应信号处理方法,以提高自适应信号处理的性能。
为了达到上述目的,本发明采用如下的技术方案:
一种基于多步梯度下降的自适应信号处理方法,所述方法利用多步梯度下降方法来加快基于标准梯度下降方法的NLMS技术的收敛速度,以此来提高信号处理的跟踪能力。
在处理方法的优选实例中,所述方法中具体用多步梯度中的动量项对基于标准梯度下降方法的NLMS技术进行加速;当基于标准梯度下降方法的NLMS技术中步长μ大小给定时,通过添加该动量项,以此控制信号处理过程中的收敛速度和最终失调。
进一步的,当梯度步长μ给定时,根据动量项步长β对信号处理过程中的收敛速度和最终失调的影响,利用滤波器的输出误差实时的调整步长β的大小。
再进一步的,在利用滤波器的输出误差实时调整步长β大小时,使得当误差较大的时候,β的值也比较大,由此可以加速滤波器的收敛速度;当误差较小的时候,则滤波器接近收敛,此时β的值比较小,使得稳态的失调变小;最终β的值趋近于零,稳态失调最终只有由步长μ控制。
再进一步的,所述梯度步长μ小于0.1。
本发明的优点不止加快了标准梯度下降的速度,滤波器最终的稳态失调也是可以控制的。这为设计一定规范的滤波器提供了指导。
附图说明
以下结合附图和具体实施方式来进一步说明本发明。
图1a为SNR=20dB时,系统未发生突变时的识别效果图;
图1b为SNR=5dB时,系统未发生突变时的识别效果图;
图2a为SNR=20dB时,在第5秒系统发生突变时的识别结果图;
图2b为SNR=5dB时,在第5秒系统发生突变时的识别结果图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
众所周知,传统自适应信号处理中的NLMS技术性能取决于其算法的步长μ:当步长取的越大,NLMS技术处理的收敛性就越快,但最终的稳态失调也越大;步长取的越小,最终的稳态失调越小,但处理的收敛速度比较慢。
由此,本发明将多步梯度下降方法和自适应信号处理的NLMS技术相结合,提出了多步梯度下降的NLMS方案。
该方案的原理为在小步长μ给定的情况下,利用多步梯度下降方法来加快基于标准梯度下降方法的NLMS技术的收敛速度,以此来提高信号处理的跟踪能力。
在基础上,本发明通过设定基于标准梯度下降方法的NLMS技术中步长μ大小,通过添加多步梯度中的动量项,以此控制信号处理过程中的收敛速度和最终失调。
具体的,当梯度步长μ设定很小时(小于0.1),根据动量项步长β对信号处理过程中的收敛速度和最终失调的影响,利用滤波器的输出误差实时的调整步长β的大小,以此来实现变步长的特点。
在具体调整步长β大小时,使得当滤波器的输出误差较大的时候,β的值也比较大,由此可以加速滤波器的收敛速度;当滤波器的输出误差较小的时候,则滤波器接近收敛,此时步长β的值比较小,使得稳态的失调变小;最终当步长β的值趋近于零,稳态失调最终只有由步长μ控制。
针对上述方案,以下通过一具体实现过程对其进行具体说明。
首先,将多步梯度下降的方法和NLMS技术相结合,更新相应的滤波器,具体的更新方程如式(1):
w ( n + 1 ) = w ( n ) + μe ( n ) x ( n ) x T ( n ) x ( n ) + δ + β [ w ( n ) - w ( n - 1 ) ] , - - - ( 1 )
其中,w(n)为滤波器系数,x(n)是输入信号,e(n)是滤波器的输出误差,δ为正则参数,一般取很小的值,μ和β分别为控制梯度的步长和控制动量项的步长。
针对式(1)可知,当步长μ(小于0.1)给定时,控制动量项的步长β对式(1)所示滤波器的收敛性、收敛速度以及稳态失调之间的影响。具体为:滤波器收敛时步长β的取值范围;当在收敛范围内取值时,步长β越大,收敛速度越快,稳态失调也越大;反之,步长β越小,收敛速度越慢,但稳态失调越小。
由此,本方案通过设定步长μ(小于0.1)的大小,并利用滤波器的输出误差实时的调整控制步长β的大小。而调整控制可采用如下的时变关系:使得当误差较大的时候,β的值也比较大,这样一来可以加速滤波器的收敛速度;当误差较小的时候,说明滤波器接近收敛,此时β的值比较小,使得稳态的失调变小。最终β的值趋近于零,稳态失调最终只有由步长μ控制。
针对式(1)所示的滤波器更新方程,其具体调整时,首先初始化滤波器的系数w(0)=0,设定正则参数δ的大小为输入信号功率,设定步长μ的大小,取μ<0.1;
接着,计算滤波器输出误差的瞬时功率,并对此功率进行平滑,估计误差的功率如下:
&sigma; ^ e 2 ( n ) = &lambda; &sigma; ^ e 2 ( n - 1 ) + ( 1 - &lambda; ) e 2 ( n ) , - - - ( 2 )
其中 &lambda; = 1 - 1 LK , ( K &GreaterEqual; 2 ) , L为滤波器的长度;
再接着,将式(2)估计出的误差功率和实际的噪声进行比较,当按照式(2)估计出的误差功率大于实际噪声功率时,将估计得出的误差功率带入式(3)的输出函数计算β值的大小;
&beta; = 1.96 [ 1 1 + e ( p - 1 ) / ( m&mu; ) - 1 2 ] , - - - ( 3 )
其中,σv为噪声实际功率,取0.001,m为输出函数的特性参数且满足条件0.05≤mμ≤0.5。
其它情况时,β的值取为零。
最后,利用第(3)步计算出β值的大小并带入式(1)中进行滤波器系数的更新。
以下再通过一些实验进一步说明本方案的实施效果。
具体将上述方案应用在系统识别的场景下,识别系统将基于多步梯度下降进行自适应信号的自动处理,并据此进行了大量的实验。其中性能指标采用自适应滤波领域中一个通用的标准,称之为归一化偏差(dB):
20 log 10 ( | | w - w ^ ( n ) | | 2 | | w | | 2 ) . - - - ( 4 )
具体使得实验结果参见图1和图2,图中MSGVSS-NLMS技术为本发明提出的技术,其它作为对比技术。
图1是在系统未发生突变时的识别效果图,其中图1a为SNR=20dB时,系统未发生突变时的识别效果图;图1b为SNR=5dB时,系统未发生突变时的识别效果图。
图2是在第5秒系统发生突变时的识别结果图,其中图2a为SNR=20dB时,在第5秒系统发生突变时的识别结果图;图2b为SNR=5dB时,在第5秒系统发生突变时的识别结果图。
通过对比可知,在不同场景中,本发明所提出的这种自适应信号技术都是十分有效的。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (5)

1.一种基于多步梯度下降的自适应信号处理方法,其特征在于,所述方法利用多步梯度下降方法来加快基于标准梯度下降方法的NLMS技术的收敛速度,以此来提高信号处理的跟踪能力。
2.根据权利要求1所述的一种基于多步梯度下降的自适应信号处理方法,其特征在于,所述方法中具体用多步梯度中的动量项对基于标准梯度下降方法的NLMS技术进行加速;当基于标准梯度下降方法的NLMS技术中步长μ大小给定时,通过添加该动量项,以此控制信号处理过程中的收敛速度和最终失调。
3.根据权利要求2所述的一种基于多步梯度下降的自适应信号处理方法,其特征在于,当梯度步长μ给定时,根据动量项步长β对信号处理过程中的收敛速度和最终失调的影响,利用滤波器的输出误差实时的调整步长β的大小。
4.根据权利要求3所述的一种基于多步梯度下降的自适应信号处理方法,其特征在于,在利用滤波器的输出误差实时调整步长β大小时,使得当误差较大的时候,β的值也比较大,由此可以加速滤波器的收敛速度;当误差较小的时候,则滤波器接近收敛,此时β的值比较小,使得稳态的失调变小;最终β的值趋近于零,稳态失调最终只有由步长μ控制。
5.根据权利要求2或3或4所述的一种基于多步梯度下降的自适应信号处理方法,其特征在于,所述梯度步长μ小于0.1。
CN201410116777.4A 2014-03-26 2014-03-26 一种基于多步梯度下降的自适应信号处理方法 Active CN103956992B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410116777.4A CN103956992B (zh) 2014-03-26 2014-03-26 一种基于多步梯度下降的自适应信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410116777.4A CN103956992B (zh) 2014-03-26 2014-03-26 一种基于多步梯度下降的自适应信号处理方法

Publications (2)

Publication Number Publication Date
CN103956992A true CN103956992A (zh) 2014-07-30
CN103956992B CN103956992B (zh) 2017-09-12

Family

ID=51334234

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410116777.4A Active CN103956992B (zh) 2014-03-26 2014-03-26 一种基于多步梯度下降的自适应信号处理方法

Country Status (1)

Country Link
CN (1) CN103956992B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017185257A1 (zh) * 2016-04-27 2017-11-02 北京中科寒武纪科技有限公司 一种用于执行Adam梯度下降训练算法的装置及方法
CN107315570A (zh) * 2016-04-27 2017-11-03 北京中科寒武纪科技有限公司 一种用于执行Adam梯度下降训练算法的装置及方法
CN111795791A (zh) * 2020-07-06 2020-10-20 浙江大学 一种变步长液压振动台自适应幅相控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1859519A (zh) * 2005-11-19 2006-11-08 华为技术有限公司 一种自适应滤波器及回波抵消器
CN101043560A (zh) * 2006-03-22 2007-09-26 北京大学深圳研究生院 回声消除器及回声消除方法
CN102903368A (zh) * 2011-07-29 2013-01-30 杜比实验室特许公司 用于卷积盲源分离的方法和设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1859519A (zh) * 2005-11-19 2006-11-08 华为技术有限公司 一种自适应滤波器及回波抵消器
CN101043560A (zh) * 2006-03-22 2007-09-26 北京大学深圳研究生院 回声消除器及回声消除方法
CN102903368A (zh) * 2011-07-29 2013-01-30 杜比实验室特许公司 用于卷积盲源分离的方法和设备

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017185257A1 (zh) * 2016-04-27 2017-11-02 北京中科寒武纪科技有限公司 一种用于执行Adam梯度下降训练算法的装置及方法
CN107315570A (zh) * 2016-04-27 2017-11-03 北京中科寒武纪科技有限公司 一种用于执行Adam梯度下降训练算法的装置及方法
CN107315570B (zh) * 2016-04-27 2021-06-18 中科寒武纪科技股份有限公司 一种用于执行Adam梯度下降训练算法的装置及方法
CN111795791A (zh) * 2020-07-06 2020-10-20 浙江大学 一种变步长液压振动台自适应幅相控制方法
CN111795791B (zh) * 2020-07-06 2021-06-04 浙江大学 一种变步长液压振动台自适应幅相控制方法

Also Published As

Publication number Publication date
CN103956992B (zh) 2017-09-12

Similar Documents

Publication Publication Date Title
CN109257030B (zh) 一种变步长lp范数LMS算法的稀疏系统辨识方法
CN110687784B (zh) 一种针对非线性不确定系统的自适应事件触发控制方法
CN108646570B (zh) 一种改进极点配置的混沌轨迹跟踪方法
CN105891810B (zh) 一种快速自适应联合时延估计方法
CN109298636B (zh) 一种改进的积分滑模控制方法
CN105306010A (zh) 一种基于最小误差熵的凸组合自适应滤波器设计方法
CN103716013B (zh) 变参数比例自适应滤波器
CN106059531B (zh) 一种非负自适应滤波器
CN103956992B (zh) 一种基于多步梯度下降的自适应信号处理方法
CN107342751B (zh) 一种基于互相关熵的变步长sa自适应滤波算法
CN109617538A (zh) 鲁棒的变参数稀疏自适应滤波器
Benallal et al. A fast convergence normalized least‐mean‐square type algorithm for adaptive filtering
Xu et al. Variable tap-length LMS algorithm based on adaptive parameters for TDL structure adaption
CN110784423A (zh) 一种基于稀疏约束的水声信道估计方法
CN109347457A (zh) 一种变参数零吸引子自适应滤波器
MENGÜÇ et al. A novel adaptive filter design using Lyapunov stability theory
Gannot et al. On the application of the unscented Kalman filter to speech processing
CN111181531B (zh) 一种变正则化偏差补偿符号子带自适应滤波器
Aref et al. Variable step size modified clipped LMS algorithm
CN103929150A (zh) 一种子带自适应滤波器的权值向量更新方法
Tang et al. A variable step-size mixed modulated Lagrange explicit time delay estimation algorithm
Shafiq et al. Direct adaptive inverse control
CN116127266B (zh) 一种基于最大类间方差的次级系数稀疏方法
Bensouda et al. Improved Fast NLMS algorithm based on variable forgetting factor technique for acoustic echo cancellation
Zhang et al. Use estimation of performance index to realize adaptive blind source separation

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