CN106301288A - 一种自适应滤波方法 - Google Patents
一种自适应滤波方法 Download PDFInfo
- Publication number
- CN106301288A CN106301288A CN201610633176.XA CN201610633176A CN106301288A CN 106301288 A CN106301288 A CN 106301288A CN 201610633176 A CN201610633176 A CN 201610633176A CN 106301288 A CN106301288 A CN 106301288A
- Authority
- CN
- China
- Prior art keywords
- connector
- momentum term
- filtering
- convergence
- factor
- 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
Links
Classifications
-
- 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
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H2021/007—Computation saving measures; Accelerating measures
Landscapes
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明提供一种自适应滤波方法,包括1)动量项LMS迭代滤波,采用梯度下降法调节权向量;2)变动量项因子LMS迭代滤波,在迭代达到收敛之前或者非平稳环境下,通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度,当逐渐达到收敛状态时,使步长相对较小从而降低稳态误差提高其稳定性;3)选择不同步长变动量项因子系统;4)双变动量项因子滤波;5)组合因子迭代;6)得到变化结果,实现滤波,得到均方误差的表示方法及动量项因子和组合因子的迭代方式,进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。
Description
技术领域
本发明属于自适应滤波领域,尤其涉及提高自适应滤波器的收敛速度及降低滤波方法稳态误差方面,一种可以平衡并优化这两方面性能的滤波方法。
背景技术
在消除回波、自动均衡和雷达声呐的波束形成等通信领域方面,及其他的识别参数、抑制噪声、谱估计等信号处理领域方面,自适应滤波器都发挥着重要的作用。而在这些领域,在实际的应用问题中,接收设备所获取的接收信号当中往往伴随着由于环境所引起的干扰和噪声,由此导致的信号误码率上升会明显影响获取信号的准确性。干扰和噪声几乎存在于现实当中的所有应用领域。而自适应滤波器所要解决的核心问题就是在信号处理过程中从充满着干扰和噪声的混合信号中估计恢复出原始的纯信号。随着时间的推移,自适应滤波器自动的调节其自身的参数以此适应外界环境的变化。过去的几十年里,数字信号处理器得到了极大的发展,速度的加快、复杂性的提高、功耗的降低等都对我们的通信领域的自适应滤波方法提出了更高的要求。在自适应滤波器领域,随着人们研究的深入,其技术理论与实际操作也越来越成熟。
在之前提出的变动量项因子自适应滤波方法当中,通过变化的动量项因子来控制LMS方法的收敛速度和稳态误差,以期望在保持其稳态误差较小的情况下提高其收敛速度,但由于固定步长参数的限制,其最终的稳态误差最理想的情况下也是与传统LMS方法的稳态误差大约相等,因此我们说它的效果还是有一定的瓶颈的,滤波方法无法更大限度的降低稳态误差。
发明内容
(1)发明目的
本发明所提出的滤波方法在加快收敛速度的同时,更大程度上降低稳态误差,最终的效果是相对于之前提出的变动量项因子自适应滤波方法,收敛速度相当,稳态误差减小,对比传统LMS方法,收敛速度加快,稳态误差减小。
(2)技术方案
本发明将两个具有不同大小步长的变动量因子系统组合起来,同时引入组合因子,通过组合因子的自适应调节来自适应的选择两个系统的步长参数在整个过程中所占的比重,具体包括以下步骤:1)动量项LMS迭代滤波,采用梯度下降法调节权向量;2)变动量项因子LMS迭代滤波,在迭代达到收敛之前或者非平稳环境下,通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度,当逐渐达到收敛状态时,使步长相对较小从而降低稳态误差提高其稳定性;3)选择不同步长变动量项因子系统;4)双变动量项因子滤波;5)组合因子迭代;6)得到变化结果,实现滤波,得到均方误差的表示方法及动量项因子和组合因子的迭代方式,进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。
(3)优点及积极效果
通过仿真,经过700次的蒙特卡洛实验,得到结果如图3,4,5,6所示。从仿真结果中可以看出,新的滤波方法在迭代次数大约为200的位置上实现收敛,与滤波方法1相当,相对而言维持了快的收敛速度,与此同时,稳态误差也比滤波方法1和滤波方法2都小,达到了作者所期待的效果。
该滤波方法通过自适应的调节两个不同步长的变动量项因子LMS系统的组合因子的方法调节步长所占比重,将两个动量项系统快速收敛与小稳态误差的优点合二为一,解决了两个系统所各自具有的缺陷,有效提高了滤波方法的整体滤波性能。
附图说明
图1变动量项因子自适应滤波的双系统步长组合滤波方法原理框图
图2滤波方法步骤流程图
图3μ1=0.012,μ2=0.003时滤波方法1,2与步长组合滤波方法的收敛性能比较
图4μ1=0.012时,滤波方法1的动量项因子变化情况
图5μ1=0.003时,滤波方法2的动量项因子变化情况
图6新步长组合滤波方法的组合因子变化情况
具体实施方式
具体实施方式包括以下步骤:
1)动量项LMS迭代滤波:为使均方误差达到最小,采用梯度下降法调节权向量,因此权向量的迭代公式为:
w(n)表示权向量,μ表示步长,J(w)表示代价函数。
其中,梯度可经过推导表示为-2p+2R·w(n),p=E[u(n)·d*(n)]为自相关向量,R=E[u(n)·uH(n)]为互相关矩阵,u(n)表示输入,d(n)表示期望响应,p,R的瞬态估计表示为u(n)·d*(n),u(n)·uH(n),代入权向量迭代公式中,得到:
为提高收敛速度,在以上公式的基础上加入一个随着不同时刻变化的量,称为动量项。公式表示为
其中e(n)是由期望响应d(n)与输出y(n)的差值所表示的估计误差,代表其动量项。
通过以上的权向量迭代公式,将各个时刻的输入与权值进行内积得到输出,输出与期望 响应的差值得到估计误差,从而得到均方误差。
2)变动量项因子LMS迭代滤波:在迭代达到收敛之前或者非平稳环境下,通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度,当逐渐达到收敛状态时,使步长相对较小从而降低稳态误差提高其稳定性。
变动量项因子LMS方法的权向量迭代公式为:
其中u(n)表示输入,e(n)是由期望响应d(n)与输出y(n)的差值所表示的估计误差,输出y(n)由输入u(n)和n时刻权向量的内积得到,代表其动量项,α(n)即代表变化的动量项因子。
3)选择不同步长变动量项因子系统:选择两个动量项系统,它们的步长参数取不同的固定值,其中μ1>μ2,也即第一个系统的速度性能高于第二个系统,而第二个系统的稳态性能要高于第一个系统,在后面的实验当中分别记为系统1和系统2,它们的权向量迭代公式分别表示为:
优选的,步长参数分别取μ1=0.012,μ2=0.003,为保证初始收敛速度和合适的均方误差使不至于发散,动量项因子α1(n),α2(n)的初值均取0.85。
4)双变动量项因子滤波系统步长组合:在这两个步长参数的基础上,组合两个不同步长参数的系统,通过动态的调整两个步长参数在组合系统中所占的比重来调节系统,在原有变动量项因子滤波方法的基础上突破常规LMS方法的步长参数对它的限制,进一步减小稳态误差。在滤波方法迭代过程中,其可以自适应的根据需要实时调节总体组合步长的大小,并最终收敛后得到最优组合步长,实现滤波方法的最优化。变动量项因子组合滤波方法的原理框图如图1所示。
自适应动态组合两个系统的组合因子用λ(n)表示,具体组合公式如下:
上式中,λ(n)表示组合因子,其取值范围为λ(n)∈[0,1]。
5)组合因子迭代实现:为使得组合因子最终达到新的步长组合滤波方法所具有的收敛速度和稳态性能,采用梯度下降的方法对组合因子的数值大小进行调节。因此自适应组合因子λ(n)的迭代公式可表示为:
上式中,ε为一很小的常数,在试验中取值0.02;J(n)是LMS方法代价函数;表示J(n) 在λ=λ(n)处的梯度。梯度可用公式表示为:
其中根据LMS方法已知
而根据上面的组合公式可以得到:
将上面两式带入到梯度的表达式当中,进而得到组合因子的自适应迭代公式:
6)得到变化结果,实现滤波:得到均方误差的表示方法及动量项因子和组合因子的迭代方式,进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。
Claims (8)
1.一种自适应滤波方法,包括以下步骤:1)动量项LMS迭代滤波,采用梯度下降法调节权向量;2)变动量项因子LMS迭代滤波,在迭代达到收敛之前或者非平稳环境下,通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度,当逐渐达到收敛状态时,使步长相对较小从而降低稳态误差提高其稳定性;3)选择不同步长变动量项因子系统;4)双变动量项因子滤波;5)组合因子迭代;6)得到变化结果,实现滤波,得到均方误差的表示方法及动量项因子和组合因子的迭代方式,进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。
2.一种如权利要求1所述的自适应滤波方法,其特征在于,所述步骤1)中权向量的迭代公式为:
w(n)表示权向量,μ表示步长,J(w)表示代价函数。
3.一种如权利要求2所述的自适应滤波方法,其特征在于,所述步骤2)中变动量项因子LMS方法的权向量迭代公式为:
其中u(n)表示输入,e(n)是由期望响应d(n)与输出y(n)的差值所表示的估计误差,输出y(n)由输入u(n)和n时刻权向量的内积得到,代表其动量项,α(n)即代表变化的动量项因子。
4.一种如权利要求3所述的自适应滤波方法,其特征在于,所述步骤3)中所述第一个系统和第二和系统权向量迭代公式分别表示为:
5.一种如权利要求4所述的自适应滤波方法,其特征在于,所述步骤4)中系统步长组合在步骤3所述两个步长参数的基础上,组合两个不同步长参数的系统,通过动态的调整两个步长参数在组合系统中所占的比重来调节系统,进一步减小稳态误差。
6.一种如权利要求5所述的自适应滤波方法,其特征在于,所述步骤5)采用梯度下降的方法对组合因子的数值大小进行调节,自适应组合因子λ(n)的迭代公式表示为:
上式中,ε为一很小的常数,在试验中取值0.02;J(n)是LMS方法代价函数;表示J(n)在λ=λ(n)处的梯度;组合因子的自适应迭代公式为
7.如权利要求6所述一种自适应滤波方法,其特征在于所述步骤(4)中,在滤波方法迭代过程中,根据需要实时调节总体组合步长的大小,并最终收敛后得到最优组合步长,实现滤波方法的最优化,自适应动态组合两个系统的组合因子用λ(n)表示,具体组合公式为
上式中,λ(n)表示组合因子,其取值范围为λ(n)∈[0,1]。
8.如权利要求7所述一种自适应滤波方法,其特征在于所述步骤(3)中,所述步长参数分别取μ1=0.012,μ2=0.003,为保证初始收敛速度和合适的均方误差使不至于发散,动量项因子α1(n),α2(n)的初值均取0.85。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610633176.XA CN106301288B (zh) | 2016-08-04 | 2016-08-04 | 一种自适应滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610633176.XA CN106301288B (zh) | 2016-08-04 | 2016-08-04 | 一种自适应滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106301288A true CN106301288A (zh) | 2017-01-04 |
CN106301288B CN106301288B (zh) | 2022-12-02 |
Family
ID=57665388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610633176.XA Active CN106301288B (zh) | 2016-08-04 | 2016-08-04 | 一种自适应滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106301288B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107248412A (zh) * | 2017-05-18 | 2017-10-13 | 歌尔科技有限公司 | 音频降噪方法和装置 |
CN107947761A (zh) * | 2017-12-18 | 2018-04-20 | 西安理工大学 | 基于最小均方四阶的变阈值比例更新自适应滤波算法 |
CN110010116A (zh) * | 2018-11-23 | 2019-07-12 | 重庆邮电大学 | 一种基于动量FxLMS算法的主动噪声控制系统 |
CN110048694A (zh) * | 2019-04-22 | 2019-07-23 | 哈尔滨理工大学 | 基于变元步长的随机傅里叶特征核最小均方算法 |
CN113556072A (zh) * | 2021-07-29 | 2021-10-26 | 西南大学 | 一种超高速电机参数在线辨识方法 |
CN113595528A (zh) * | 2021-07-30 | 2021-11-02 | 重庆长安汽车股份有限公司 | 一种基于幂函数的自适应变步长lms滤波器及其实现方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178846A (zh) * | 2013-03-29 | 2013-06-26 | 华南理工大学 | 一种用于adc校准的lms算法 |
CN104009734A (zh) * | 2014-05-07 | 2014-08-27 | 西安理工大学 | 梯度变步长lms自适应滤波方法 |
CN104283528A (zh) * | 2014-09-18 | 2015-01-14 | 河海大学 | 一种变步长lms自适应滤波方法 |
US20150263701A1 (en) * | 2013-04-10 | 2015-09-17 | King Fahd University Of Petroleum And Minerals | Adaptive filter for system identification |
-
2016
- 2016-08-04 CN CN201610633176.XA patent/CN106301288B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178846A (zh) * | 2013-03-29 | 2013-06-26 | 华南理工大学 | 一种用于adc校准的lms算法 |
US20150263701A1 (en) * | 2013-04-10 | 2015-09-17 | King Fahd University Of Petroleum And Minerals | Adaptive filter for system identification |
CN104009734A (zh) * | 2014-05-07 | 2014-08-27 | 西安理工大学 | 梯度变步长lms自适应滤波方法 |
CN104283528A (zh) * | 2014-09-18 | 2015-01-14 | 河海大学 | 一种变步长lms自适应滤波方法 |
Non-Patent Citations (2)
Title |
---|
曾召华: "一种具有双瞬变因子的LMS自适应滤波算法", 《现代雷达》 * |
李晓波 等: "双变因子LMS自适应滤波算法", 《火力与指挥控制》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107248412A (zh) * | 2017-05-18 | 2017-10-13 | 歌尔科技有限公司 | 音频降噪方法和装置 |
CN107947761A (zh) * | 2017-12-18 | 2018-04-20 | 西安理工大学 | 基于最小均方四阶的变阈值比例更新自适应滤波算法 |
CN107947761B (zh) * | 2017-12-18 | 2021-09-10 | 西安理工大学 | 基于最小均方四阶的变阈值比例更新自适应滤波方法 |
CN110010116A (zh) * | 2018-11-23 | 2019-07-12 | 重庆邮电大学 | 一种基于动量FxLMS算法的主动噪声控制系统 |
CN110048694A (zh) * | 2019-04-22 | 2019-07-23 | 哈尔滨理工大学 | 基于变元步长的随机傅里叶特征核最小均方算法 |
CN110048694B (zh) * | 2019-04-22 | 2022-09-23 | 哈尔滨理工大学 | 基于变元步长的随机傅里叶特征核最小均方算法 |
CN113556072A (zh) * | 2021-07-29 | 2021-10-26 | 西南大学 | 一种超高速电机参数在线辨识方法 |
CN113595528A (zh) * | 2021-07-30 | 2021-11-02 | 重庆长安汽车股份有限公司 | 一种基于幂函数的自适应变步长lms滤波器及其实现方法 |
CN113595528B (zh) * | 2021-07-30 | 2023-08-04 | 重庆长安汽车股份有限公司 | 一种基于幂函数的自适应变步长lms滤波器及其实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106301288B (zh) | 2022-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106301288A (zh) | 一种自适应滤波方法 | |
Heerikhuisen et al. | The effect of new interstellar medium parameters on the heliosphere and energetic neutral atoms from the interstellar boundary | |
CN103888145B (zh) | 一种重构信号的方法 | |
CN102139766B (zh) | 一种亚轨道飞行器再入飞行的设计攻角获取方法 | |
CN103413557B (zh) | 语音信号带宽扩展的方法和装置 | |
Fisk et al. | On theories for stochastic acceleration in the solar wind | |
CN108364659A (zh) | 基于多目标优化的频域卷积盲信号分离方法 | |
CN106330219A (zh) | 一种基于递归网络的无源互调干扰对消系统及方法 | |
CN109375253A (zh) | 基于全部发震构造最大可信地震的地震动参数评价方法 | |
CN104036331A (zh) | 一种基于改进粒子群算法的电力系统动态经济调度方法 | |
CN107947761A (zh) | 基于最小均方四阶的变阈值比例更新自适应滤波算法 | |
CN103888104A (zh) | Fir数字滤波器设计方法和系统 | |
CN108169794B (zh) | 近断层人工地震波与规范反应谱匹配的调整方法 | |
CN104866866A (zh) | 改进的自然梯度变步长盲源分离算法 | |
CN106297778A (zh) | 数据驱动的基于奇异值分解的神经网络声学模型裁剪方法 | |
CN110059286A (zh) | 一种基于fft的结构非平稳响应高效分析方法 | |
CN105895104A (zh) | 说话人自适应识别方法及系统 | |
CN105812116B (zh) | 全双工时延估计方法和全双工时延估计装置 | |
CN109274352A (zh) | 基于最大相关熵的多凸组合自适应滤波方法 | |
CN106503456B (zh) | 基于超球体变换的集合卡尔曼滤波油藏动态历史拟合方法 | |
CN108613737A (zh) | 基于小波包与stft的飞行器多频振动信号的辨识方法 | |
Ali et al. | An improved gain vector to enhance convergence characteristics of recursive least squares algorithm | |
CN106878211A (zh) | 一种卫星信道多模盲均衡算法 | |
CN102006105B (zh) | 深空接收天线组阵相关加权方法及系统 | |
Merritt | Gravitational Encounters and the Evolution of Galactic Nuclei. III. Anomalous Relaxation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |