申请/专利权人：北京中科海讯数字科技股份有限公司

申请日：2016-08-04

公开（公告）日：2022-12-02

公开（公告）号：CN106301288B

主分类号：H03H21/00

分类号：H03H21/00

优先权：

专利状态码：有效-授权

法律状态：2022.12.02#授权;2019.08.16#实质审查的生效;2017.01.04#公开

摘要：本发明提供一种自适应滤波方法，包括1动量项LMS迭代滤波，采用梯度下降法调节权向量；2变动量项因子LMS迭代滤波，在迭代达到收敛之前或者非平稳环境下，通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度，当逐渐达到收敛状态时，使步长相对较小从而降低稳态误差提高其稳定性；3选择不同步长变动量项因子系统；4双变动量项因子滤波；5组合因子迭代；6得到变化结果，实现滤波，得到均方误差的表示方法及动量项因子和组合因子的迭代方式，进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。

主权项：1.一种自适应滤波方法，应用于通信领域，包括以下步骤：1动量项LMS迭代滤波，采用梯度下降法调节权向量；2变动量项因子LMS迭代滤波，在迭代达到收敛之前或者非平稳环境下，通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度，当逐渐达到收敛状态时，使步长相对较小从而降低稳态误差提高其稳定性；3选择不同步长变动量项因子系统；4双变动量项因子滤波；5组合因子迭代；6得到变化结果，实现滤波，得到均方误差的表示方法及动量项因子和组合因子的迭代方式，进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果；其中所述步骤3中，选择第一动量项系统与第二动量项系统，所述第一动量项系统与第二动量项系统的步长参数分别为μ1和μ2，且μ1μ2；所述步骤4中系统步长组合在步骤3所述第一动量项系统与第二动量项系统的步长参数的基础上，组合两个不同步长参数的系统，通过动态的调整两个步长参数在组合系统中所占的比重来调节系统，进一步减小稳态误差；所述步骤4中，在滤波方法迭代过程中，根据需要实时调节总体组合步长的大小，并最终收敛后得到最优组合步长，实现滤波方法的最优化，自适应动态组合两个系统的组合因子用λn表示，具体组合公式为 上式中，wn表示权向量，λn表示组合因子，其取值范围为λn∈[0,1]；所述步骤5采用梯度下降的方法对组合因子的数值大小进行调节，自适应组合因子λn的迭代公式表示为： 上式中，ε为一很小的常数，在试验中取值0.02；Jn是LMS方法代价函数；表示Jn在λ＝λn处的梯度；组合因子的自适应迭代公式为 其中un表示输入，en是由期望响应dn与输出yn的差值所表示的估计误差，输出yn由输入un和n时刻权向量的内积得到,w1n与w2n分别表示所述第一动量项系统与第二动量项系统的权向量。

全文数据：_种自适应滤波方法技术领域[0001]本发明属于自适应滤波领域，尤其涉及提高自适应滤波器的收敛速度及降低滤波方法稳态误差方面，一种可以平衡并优化这两方面性能的滤波方法。背景技术[0002]在消除回波、自动均衡和雷达声呐的波束形成等通信领域方面，及其他的识别参数、抑制噪声、谱估计等信号处理领域方面，自适应滤波器都发挥着重要的作用。而在这些领域，在实际的应用问题中，接收设备所获取的接收信号当中往往伴随着由于环境所引起的干扰和噪声，由此导致的信号误码率上升会明显影响获取信号的准确性。干扰和噪声几乎存在于现实当中的所有应用领域。而自适应滤波器所要解决的核心问题就是在信号处理过程中从充满着干扰和噪声的混合信号中估计恢复出原始的纯信号。随着时间的推移，自适应滤波器自动的调节其自身的参数以此适应外界环境的变化。过去的几十年里，数字信号处理器得到了极大的发展，速度的加快、复杂性的提高、功耗的降低等都对我们的通信领域的自适应滤波方法提出了更高的要求。在自适应滤波器领域，随着人们研究的深入，其技术理论与实际操作也越来越成熟。[0003]在之前提出的变动量项因子自适应滤波方法当中，通过变化的动量项因子来控制LMS方法的收敛速度和稳态误差，以期望在保持其稳态误差较小的情况下提高其收敛速度，但由于固定步长参数的限制，其最终的稳态误差最理想的情况下也是与传统LMS方法的稳态误差大约相等，因此我们说它的效果还是有一定的瓶颈的，滤波方法无法更大限度的降低稳态误差。发明内容[0004]1发明目的[0005]本发明所提出的滤波方法在加快收敛速度的同时，更大程度上降低稳态误差，最终的效果是相对于之前提出的变动量项因子自适应滤波方法，收敛速度相当，稳态误差减小，对比传统LMS方法，收敛速度加快，稳态误差减小。[0006]2技术方案[0007]本发明将两个具有不同大小步长的变动量因子系统组合起来，同时引入组合因子，通过组合因子的自适应调节来自适应的选择两个系统的步长参数在整个过程中所占的比重，具体包括以下步骤:1动量项LMS迭代滤波，采用梯度下降法调节权向量;2变动量项因子LMS迭代滤波，在迭代达到收敛之前或者非平稳环境下，通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度，当逐渐达到收敛状态时，使步长相对较小从而降低稳态误差提高其稳定性;3选择不同步长变动量项因子系统;4双变动量项因子滤波;5组合因子迭代;6得到变化结果，实现滤波，得到均方误差的表示方法及动量项因子和组合因子的迭代方式，进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。[0008]3优点及积极效果[0009]通过仿真，经过700次的蒙特卡洛实验，得到结果如图3,4,5,6所示。从仿真结果中可以看出，新的滤波方法在迭代次数大约为200的位置上实现收敛，与滤波方法1相当，相对而言维持了快的收敛速度，与此同时，稳态误差也比滤波方法1和滤波方法2都小，达到了作者所期待的效果。[0010]该滤波方法通过自适应的调节两个不同步长的变动量项因子LMS系统的组合因子的方法调节步长所占比重，将两个动量项系统快速收敛与小稳态误差的优点合二为一，解决了两个系统所各自具有的缺陷，有效提高了滤波方法的整体滤波性能。附图说明[0011]图1变动量项因子自适应滤波的双系统步长组合滤波方法原理框图[0012]图2滤波方法步骤流程图[0013]图3μι=0·012，μ2=0·003时滤波方法1，2与步长组合滤波方法的收敛性能比较[0014]图4yi=0.012时，滤波方法1的动量项因子变化情况[0015]图5yi=0.003时，滤波方法2的动量项因子变化情况[0016]图6新步长组合滤波方法的组合因子变化情况具体实施方式具体实施方式[0017]包括以下步骤：[0018]1动量项LMS迭代滤波:为使均方误差达到最小，采用梯度下降法调节权向量，因此权向量的迭代公式为：[0020]wn表示权向量，μ表示步长，Jw表示代价函数。[0021]可经过推导表示为-2p+2R·wn，p=E[un·d*n]为自相关向量，R=E[un·uHn]为互相关矩阵，un表示输入，dn表示期望响应，p，R的瞬态估计表示为un·cfn，un·uHn，代入权向量迭代公式中，得到：[0023]为提高收敛速度，在以上公式的基础上加入一个随着不同时刻变化的量，称为动量项。公式表示为[0025]其中en是由期望响应dn与输出yn的差值所表示的估计误差，代表其动量项。[0026]通过以上的权向量迭代公式，将各个时刻的输入与权值进行内积得到输出，输出与期望响应的差值得到估计误差，从而得到均方误差。[0027]2变动量项因子LMS迭代滤波:在迭代达到收敛之前或者非平稳环境下，通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度，当逐渐达到收敛状态时，使步长相对较小从而降低稳态误差提高其稳定性。[0028]变动量项因子LMS方法的权向量迭代公式为：[0030]其中un表示输入，eη是由期望响应dn与输出yη的差值所表示的估计误差，输出yη由输入uη和η时刻权向量的内积得到，|V"-i]代表其动量项，αη即代表变化的动量项因子。[0031]3选择不同步长变动量项因子系统:选择两个动量项系统，它们的步长参数取不同的固定值，其中μιμ2,也即第一个系统的速度性能高于第二个系统，而第二个系统的稳态性能要高于第一个系统，在后面的实验当中分别记为系统1和系统2,它们的权向量迭代公式分别表不为：[0034]优选的，步长参数分别取μι=0·012，μ2=0·003，为保证初始收敛速度和合适的均方误差使不至于发散，动量项因子aKnhc^n的初值均取0.85。[0035]4双变动量项因子滤波系统步长组合:在这两个步长参数的基础上，组合两个不同步长参数的系统，通过动态的调整两个步长参数在组合系统中所占的比重来调节系统，在原有变动量项因子滤波方法的基础上突破常规LMS方法的步长参数对它的限制，进一步减小稳态误差。在滤波方法迭代过程中，其可以自适应的根据需要实时调节总体组合步长的大小，并最终收敛后得到最优组合步长，实现滤波方法的最优化。变动量项因子组合滤波方法的原理框图如图1所示。[0036]自适应动态组合两个系统的组合因子用λη表示，具体组合公式如下：[0038]上式中，λη表示组合因子，其取值范围为ληε[0，1]。[0039]5组合因子迭代实现:为使得组合因子最终达到新的步长组合滤波方法所具有的收敛速度和稳态性能，采用梯度下降的方法对组合因子的数值大小进行调节。因此自适应组合因子λη的迭代公式可表示为：[0041]上式中，ε为一很小的常数，在试验中取值0.02;Jη是LMS方法代价函数;▽,J?表示Jn在λ=λη处的梯度。梯度可用公式表示为：[0045]而根据上面的组合公式可以得到：[0047]将上面两式带入到梯度的表达式当中，进而得到组合因子的自适应迭代公式：[0049]6得到变化结果，实现滤波:得到均方误差的表示方法及动量项因子和组合因子的迭代方式，进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。

权利要求：1.一种自适应滤波方法，包括以下步骤：1动量项LMS迭代滤波，采用梯度下降法调节权向量；2变动量项因子LMS迭代滤波，在迭代达到收敛之前或者非平稳环境下，通过自适应调节组合因子使步长相对较大从而提高滤波方法收敛速度，当逐渐达到收敛状态时，使步长相对较小从而降低稳态误差提高其稳定性;3选择不同步长变动量项因子系统;4双变动量项因子滤波;5组合因子迭代;6得到变化结果，实现滤波，得到均方误差的表示方法及动量项因子和组合因子的迭代方式，进而得到均方误差的收敛特性及动量项因子和组合因子的自适应变化结果。2.-种如权利要求1所述的自适应滤波方法，其特征在于，所述步骤1中权向量的迭代公式为：wn表示权向量，μ表示步长，Jw表示代价函数。3.-种如权利要求2所述的自适应滤波方法，其特征在于，所述步骤2中变动量项因子LMS方法的权向量迭代公式为其中un表示输入，en是由期望响应dn与输出yn的差值所表示的估计误差，输出yη由输入uη和η时刻权向量的内积得到，]代表其动量项，an即代表变化的动量项因子。4.一种如权利要求3所述的自适应滤波方法，其特征在于，所述步骤3中所述第一个系统和第二和系统权向量迭代公式分别表不为：5.-种如权利要求4所述的自适应滤波方法，其特征在于，所述步骤4中系统步长组合在步骤3所述两个步长参数的基础上，组合两个不同步长参数的系统，通过动态的调整两个步长参数在组合系统中所占的比重来调节系统，进一步减小稳态误差。6.-种如权利要求5所述的自适应滤波方法，其特征在于，所述步骤5采用梯度下降的方法对组合因子的数值大小进行调节，自适应组合因子λη的迭代公式表示为：上式中，ε为一很小的常数，在试验中取值〇.02;Jη是LMS方法代价函数;▽,J«'表示Jη在λ=λη处的梯度;组合因子的自适应迭代公式为7.如权利要求6所述一种自适应滤波方法，其特征在于所述步骤4中，在滤波方法迭代过程中，根据需要实时调节总体组合步长的大小，并最终收敛后得到最优组合步长，实现滤波方法的最优化，自适应动态组合两个系统的组合因子用λη表示，具体组合公式为上式中，λη表示组合因子，其取值范围为ληe[〇，1]。8.如权利要求7所述一种自适应滤波方法，其特征在于所述步骤3中，所述步长参数分别取^=0.012，此=0.003，为保证初始收敛速度和合适的均方误差使不至于发散，动量项因子〇11〇，€1211的初值均取〇.85。

百度查询： 北京中科海讯数字科技股份有限公司 一种自适应滤波方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD