CN106788337B - 稳健的仿射投影符号自适应滤波算法 - Google Patents

稳健的仿射投影符号自适应滤波算法 Download PDF

Info

Publication number
CN106788337B
CN106788337B CN201611059496.5A CN201611059496A CN106788337B CN 106788337 B CN106788337 B CN 106788337B CN 201611059496 A CN201611059496 A CN 201611059496A CN 106788337 B CN106788337 B CN 106788337B
Authority
CN
China
Prior art keywords
algorithm
noise
step size
proportional
vector
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
CN201611059496.5A
Other languages
English (en)
Other versions
CN106788337A (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.)
Shenyang University of Technology
Original Assignee
Shenyang University of Technology
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 Shenyang University of Technology filed Critical Shenyang University of Technology
Publication of CN106788337A publication Critical patent/CN106788337A/zh
Application granted granted Critical
Publication of CN106788337B publication Critical patent/CN106788337B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

Landscapes

  • Filters That Use Time-Delay Elements (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

稳健的仿射投影符号自适应滤波算法，该方法引入步长函数，并将变步长的方法和比例矩阵的思想融合到一起，即变步长的改进比例仿射投影符号算法,本申请将变步长的方法和比例矩阵的思想融合到一起，引入步长函数，提出了一种稳健的仿射投影符号自适应滤波算法——变步长的改进比例仿射投影符号算法（VSS‑IPAPSA，Variable Step‑Size Improved Proportionate Affine Projection Sign Algorithm）。该算法不仅可以缓解收敛速度与稳态失调之间的矛盾，同时也可以增加其对系统的不同稀疏特性和噪声特性的适应性。理论分析和仿真结果验证了其稳健性和有效性。

Description

稳健的仿射投影符号自适应滤波算法
技术领域：本发明涉及一种稳健的仿射投影符号自适应滤波算法。
背景技术：非高斯的冲击噪声在现实世界广泛存在，严重影响了基于l2范数优化准则的自适应滤波算法的性能。在各类自适应滤波算法中仿射投影符号算法(APSA，AffineProjection Sign Algorithm)结合了仿射投影算法(APA，Affine Projection Algorithm)良好的收敛特性和符号算法对非高斯冲击噪声干扰的抑制能力，因而其在非高斯冲击噪声条件下具有良好的性能。但是该算法的步长选择是固定的，且未考虑系统稀疏特性，因而在参数选择和收敛速度方面有一定的局限性。具体的说，自适应滤波算法在很多领域都有广泛的应用，如声学和网络回声消除、噪声抑制，信道估计等[1-4]。文献[5]提出的仿射投影算法(APA，Affine Projection Algorithm)及其改进算法是一类重要的自适应滤波算法，该类算法在输入信号相关性较高的情况下仍具有良好的收敛性能。目前大多数算法是基于高斯噪声假设和l2范数优化准则的。然而，现实世界中广泛存在着各类非高斯的冲击噪声，如图像中的椒盐噪声、电力开关中的脉动噪声等，这些噪声的存在破坏了基于l2范数优化准则的自适应滤波算法的性能。为了提高算法对非高斯噪声的抑制性能，文献[6-8]相继提出了一些代替基于l2范数优化准则的自适应滤波算法，如基于l1范数优化准则的仿射投影符号算法(APSA，Affine Projection Sign Algorithm)，该算法在存在非高斯噪声干扰时具有良好的稳健性，但是该算法的步长选择是固定的，因而在参数选择和收敛速度方面有一定的局限性，因此，固定步长参数的选择势必要在收敛性能、跟踪性能、稳态误差等之间做出妥协。针对上述问题，文献[9]引入了变步长的概念，提出了基于最小均方偏差(MSD，MeanSquare Deviation)准则的变步长仿射投影符号算法(VSS-APSA，Variable Step-SizeAPSA)，该算法通过引入步长函数，采用随机逼近和移动平均的方法得到算法收敛的最优步长，该算法很大程度上减少了收敛速度与稳态失调之间的矛盾。
然而，上述算法均并未考虑系统的稀疏特性[10](车载电话、电话会议等系统中声环境的温度和压力的变化或者电话持有者的位置变化会导致系统的稀疏性在一个较大的范围内变化)。
上述参考文献为：
[1]HAYKIN S.Adaptive filter theory[M].Upper Saddle River，NJ，USA：Prentice-Hall，2002.
[2]蔡宇，洪缨，原建平，等.语音系统中的子带自适应回声消除技术[J].仪器仪表学报，2013，34(7)：1448-1453.
[3]张红梅，韩万刚.一种新的变步长LMS自适应滤波算法研究及其应用[J].仪器仪表学报，2015，36(8)：1822-1830.
[4]倪锦根，马兰申.分布式子带自适应滤波算法[J].电子学报，2015，43(11)：2225-2231.
[5]OZEKI K，UMEDA T.An adaptive filtering algorithm using anorthogonal projection to an affine subspace and its properties[J].ElectronicsCommunication，1984，67(A)：19-27.
[6]XIAO L S，WU M，YANG J.A new efficient filtered-x affine projectionsign algorithm for active control of impulsive noise[J].Signal Processing，2016，120(3)：456-461.
[7]NI J G，LI F.Efficient implementation of the affine projection signalgorithm[J].IEEE Signal Processing Letters，2012，19(1)：24-26.
[8]SHAO TIANGE，ZHENG Y R，BENESTY J.An affine projection signalgorithm robust against impulsive interferences[J].IEEE Signal ProcessingLetters，2010，17(4)：327-330.
[9]SHIN J，YOO J，PARK P.Variable step-size affine projection signalgorithm[J].Electronics Letters，2012,48(9)：483-485.
[10]郭莹，侯明云.基于指数梯度和凸组合的稀疏自适应滤波算法[J].仪器仪表学报，2014，35(4)：894-900。
发明内容：
发明目的：发明提供一种稳健的仿射投影符号自适应滤波算法，其目的是解决以往所存在的问题。
技术方案：发明是通过以下技术方案实现的：
一种稳健的仿射投影符号自适应滤波算法，其特征在于：该方法引入步长函数，并将变步长的方法和比例矩阵的思想融合到一起，即变步长的改进比例仿射投影符号算法(VSS-IPAPSA,Variable Step-Size Improved Proportionate Affine Projection SignAlgorithm)。具体的，本申请主要进行了三方面的工作：1)引入步长函数的概念，采用随机逼近的方法推导得最优步长；2)在实际情况分析过程中，采用移动平均法处理期望项，得出准确的可变步长更新公式；3)将比例矩阵的思想引入传统的仿射投影算法中，提高算法在不同稀疏系统中的适用性，同时利用符号算法对非高斯噪声干扰的抑制能力。
该方法中将比例归一化最小均方算法(PNLMS，Proportionate Normalized LeastMean Square)中成比例的方法引入仿射投影符号算法(APSA，Affine Projection SignAlgorithm)算法中，得到比例仿射投影符号算法(PAPSA,Proportionate AffineProjection Sign Algorithm)算法的更新公式：
Figure GDA0002578613180000041
其中，sgn(·)表示符号运算，μ是步长因子，L是滤波器长度，n是时间系数，
Figure GDA0002578613180000042
表示所估计的自适应滤波器权系数向量，用于辨识未知系统w0，输入信号向量为x(n)＝[x(n)x(n-1)…x(n-L+1)]T，未知系统的期望响应信号为d(n)＝xT(n)w0+v(n)(v(n)是噪声)，则先验输出误差向量和后验输出误差向量分别为
Figure GDA0002578613180000043
和
Figure GDA0002578613180000044
d(n)＝[d(n)d(n-1)…d(n-M+1)]T是期望输出信号向量，M为投影阶数，X(n)是APA滤波器结构中包含M维输入向量的输入矩阵，即X(n)＝[x(n)x(n-1)…x(n-M+1)]，L×L对角矩阵G(n)表示成比例矩阵，gi(n)是其对角矩阵G(n)的对角元素，表示滤波器第i个抽头参数在n时刻的成比例比重，G(n)和gi(n)的计算方法如下:
G(n)＝diag{g0(n),…gL-1(n)} (2)
Figure GDA0002578613180000045
Figure GDA0002578613180000046
Figure GDA0002578613180000047
其中，max(·)表示求最大值，||·||∞表示求无穷范数，||·||2表示求2范数，参数δp是一个很小的正数，它在所有滤波器抽头权值为零的时候启动更新；ρ和f(n)能够防止远小于最大抽头的权值停止更新；
同理，将IPNLMS算法的思想应用于APSA算法中，得到本申请的IPAPSA算法，其与PAPSA算法的权系数更新公式相同，不同之处在于对角元素上,
Figure GDA0002578613180000048
式中:||·||1表示求1范数，θ是与系统稀疏度有关的参数，稀疏度越大表示回声路径越稀疏，ε是一个比较小的正数，以防止分母为0。
(1)最优步长推导
本申请的算法推导是基于最小均方偏差(MSD，Mean Square Deviation)准则的。定义滤波器系数误差向量
Figure GDA0002578613180000051
并将固定步长μ替代为可变步长μ(n)。为了得到最优步长，与传统方法不同，本申请不作噪声向量与误差向量相关性的假设，即不删除噪声项，采用随机逼近的方法推导出步长，具体推导过程如下：
对式
Figure GDA0002578613180000052
两边取2范数平方的期望得：
Figure GDA0002578613180000053
其中，f(μ(n))是关于步长的一个函数，且
Figure GDA0002578613180000054
v(n)是噪声向量，符合
Figure GDA0002578613180000055
表示定义为。
为了使MSD的值达到最小化，本申请的改进算法从n到n+1的迭代过程中通过选择最优步长来使步长函数最小化，即步长函数写为：
Figure GDA0002578613180000056
当自适应滤波器收敛到最佳状态时，
Figure GDA0002578613180000057
此时，sgn(eT(n))≈sgn(v(n))，且在稳定状态时：
Figure GDA0002578613180000061
由于准确地计算出sgn(eT(n))v(n)的值比较困难，所以不能直接获得f(μ(n))；因此，本申请根据sgn(eT(n))v(n)的上限，运用随机逼近方法得到下式：
Figure GDA0002578613180000062
但是，v(n)的绝对值不是精确的测量值，本申请将v(n)的绝对值近似为其期望值，v(n)有半正态分布的属性，且
Figure GDA0002578613180000063
即将||v(n)||1改为：
Figure GDA0002578613180000064
由式(12)得，
Figure GDA0002578613180000065
其中，M为投影阶数，σv为噪声信号v(n)的方差。
运用随机逼近的方法，步长函数的上限重写为：
Figure GDA0002578613180000066
从n到n+1的迭代过程中，最小化步长函数f(μ(n))大幅度降低了MSD的值；因此，对式(14)关于μ(n)求偏导，得，
Figure GDA0002578613180000067
令上式求导结果为0，得到的派生步长为：
Figure GDA0002578613180000071
(2)最优步长分析
通过式(16)推导出了步长最小化步长函数的上限，然而由于式(16)中期望项的存在，要得到准确的步长是非常困难的；需要考虑以下两种情况，首先定义：
Figure GDA0002578613180000072
情况1：
Figure GDA0002578613180000073
即β(n)>0时，步长因子按本申请中的迭代公式进行更新；
情况2:当
Figure GDA0002578613180000074
即β(n)≤0时，此时μ(n)停止更新，即μ(n)＝μ(n-1)。
本申请采用移动平均法来处理上式中的期望项，得到如下：
Figure GDA0002578613180000075
其中，α(0≤α＜1)为平滑因子，本申请平滑因子的值设为
Figure GDA0002578613180000076
k为一个恒定的值(本申请中k取4)；
因此，本申请所提出的VSS-IPAPSA算法的滤波器系数更新公式为：
Figure GDA0002578613180000077
G(n)＝diag{g0(n),…gL-1(n)}
Figure GDA0002578613180000078
Figure GDA0002578613180000081
Figure GDA0002578613180000082
本申请提出的算法总结如表1所示：
表1 VSS-IPAPSA算法总结
Figure GDA0002578613180000083
Figure GDA0002578613180000091
由于上述迭代过程中需要估算噪声信号v(n)的方差σv，考虑到冲击噪声干扰对噪声方差估算方法的影响，本申请通过使用中值滤波器来提高输出噪声方差对冲击噪声的稳健性，过程如下：
if mod(n,Nw)＝0
Figure GDA0002578613180000092
Figure GDA0002578613180000093
Figure GDA0002578613180000094
Figure GDA0002578613180000101
end
其中，Ae(n)＝[e(n)e(n-1)…e(n-Nw+1)]，Be(n)＝[e2(n)e2(n-1)…E2(N-Nw+1)]。
Figure GDA0002578613180000102
表示n时刻估计的误差方差的平方，
Figure GDA0002578613180000103
表示n时刻估计的输入信号方差的平方，
Figure GDA0002578613180000104
表示n时刻估计的噪声方差的平方，
Figure GDA0002578613180000105
表示n时刻所估计的中间向量，模mod(n,Nw)表示整数n和Nw之间的余数，Nw为估计窗的长度，
Figure GDA0002578613180000106
为有限样值校正因子，目的是减少输出噪声方差估计算法在每n次迭代时计算复杂性。
优点效果：
本申请将变步长的方法和比例矩阵的思想融合到一起，引入步长函数，提出了一种稳健的仿射投影符号自适应滤波算法——变步长的改进比例仿射投影符号算法(VSS-IPAPSA，Variable Step-Size Improved Proportionate Affine Projection SignAlgorithm)。该算法不仅可以缓解收敛速度与稳态失调之间的矛盾，同时也可以增加其对系统的不同稀疏特性和噪声特性的适应性。理论分析和仿真结果验证了其稳健性和有效性。
附图说明：
图1为不同特征指数条件下的概率密度函数。
图2(a)、(b)、(c)、(d)分别为α＝0.5,1.0,1.5,2.0时的α-稳定分布；
图3为自适应滤波器系统辨识原理图；
图4为仿真中用的信道，图4中(a)和(b)分别为仿真中所用的稀疏信道和非稀疏信道，(a)稀疏度为0.85349的稀疏信道，(b)稀疏度为0.32883的非稀疏信道；
图5为高斯噪声条件下各类算法在非稀疏系统下性能比较；
图6为非高斯噪声条件下各个算法在非稀疏系统下性能比较；
图7为非高斯噪声条件下各个算法在稀疏系统下性能比较；
图8为非高斯噪声条件下各个算法在不同系统下的跟踪性能。
具体实施方式：
如图1所示，发明提供一种稳健的仿射投影符号自适应滤波算法，
传统算法分析
α-稳定分布
α-稳定分布作为非高斯冲击噪声模型既可以满足随机噪声产生过程的合理假设又计算方便，从而可以当作信号处理中噪声的理想模型来使用。α-稳定分布有很多种定义方式[13-14]，本申请主要介绍基于特征函数的定义[15]。
α-稳定分布的概率密度函数并没有统一的闭合形式，但它的特征函数具有统一的表达式，这是表示α-稳定分布最简便的方法。若一个随机变量X服从α-稳定分布，那么其特征函数可以这样描述：
Figure GDA0002578613180000111
式中，-∞＜δ＜∞，γ>0，0＜α≤2，-1≤β≤1，四个参数δ，γ，α，β分别代表不同的物理意义，δ表示位置参数，γ表示尺度参数，α表示特征值数，用来表征α-稳定分布概率密度函数拖尾的厚重程度(当α＝2时为高斯分布)，β为对称参数(当β＝0时，稳定分布是关于δ的对称分布，称为对称α-稳定分布(SαS))。
本申请所用的非高斯冲击噪声模型是用α-稳定分布来描述，图1所示为标准SαS分布在不同特征指数α下的概率密度函数曲线图，图2分别给出了不同α值的SαS稳定噪声。
由图1，图2可以看到，当α＝2时，SαS概率密度函数曲线实质上为均值为零且方差为2的高斯分布，由此可见，SαS概率密度函数曲线拥有许多的特征与高斯分布相同，例如曲线光滑、单峰分布、关于中值对称、钟形等特征。而与高斯分布不同之处在于SαS概率密度函数的拖尾要比高斯分布的更厚重，且特征指数的值越小，其拖尾越厚重，即大幅度样本发生的概率越大。因此，SαS分布非常适合来描述那些类似于高斯分布，却具有很强冲击性的非高斯分布。由于非高斯冲击噪声条件下的稳定分布不存在二阶距，因此服从SαS分布的随机变量的方差也不存在，即该分布破坏了基于l2范数优化准则的自适应滤波算法的性能。
APSA
本申请以系统辨识[4]为例来分析各类自适应滤波算法的性能，图3为自适应滤波器系统辨识原理图。
图3中，
Figure GDA0002578613180000121
表示所估计的自适应滤波器权系数向量，用于辨识未知系统w0，N是滤波器长度，n是时间系数，输入信号向量为x(n)＝[x(n)x(n-1)…x(n-N+1)]T，则
未知系统的期望响应信号为：
d(n)＝xT(n)w0+v(n) (21)
设先验输出误差向量和后验输出误差向量分别为：
Figure GDA0002578613180000122
Figure GDA0002578613180000123
式(22)和(23)中，d(n)＝[d(n)d(n-1)…d(n-M+1)]T是期望输出信号向量，M为投影阶数，X(n)是APA滤波器结构中包含M维输入向量的输入矩阵，即X(n)＝[x(n)x(n-1)…x(n-M+1)]。
传统APSA的滤波器系数更新公式可根据式(24)和(25)推导：
Figure GDA0002578613180000131
Figure GDA0002578613180000132
式中，||.||1表示某矢量的l1范数，τ2是保证滤波器系数在一次迭代中不会发生太大的变化的参数，一般τ2取较小值。利用拉格朗日乘数法，APSA算法的权系数更新公式为：
Figure GDA0002578613180000133
其中，μ为步长参数，式(26)中用先验误差向量e(n)代替后验误差向量ep(n)。这是由于ep(n)在实践中无法获得，而先验误差向量e(n)是对其较好的近似。
由式(26)可见，APSA利用输出误差向量e(n的符号进行滤波器系数更新，因此，它可有效地抑制非高斯冲激噪声干扰。此外，APSA通过重复利用过去的输入向量提高了算法在输入信号严重相关时的收敛性能。但是，式(26)中的固定步长使得APSA算法不能同时满足高速收敛与低稳态误差的性能需求，VSS-APSA算法[9]可有效改善这一缺陷。
VSS-APSA
该算法的推导是基于最小MSD准则的。定义滤波器系数误差向量
Figure GDA0002578613180000134
并将固定步长μ替代为可变步长μ(n)。据此，式(26)可改写为
Figure GDA0002578613180000135
对式(27)两边取2范数平方的期望，并对所得等式关于μ(n)求导并令其为零，得变步长表达式：
Figure GDA0002578613180000136
其中，
Figure GDA0002578613180000141
式(28)可进一步变换为
Figure GDA0002578613180000142
该方法首先假定噪声向量与误差向量不相关即
Figure GDA0002578613180000143
其次用μ(n)的时间平均替代统计平均，从而进一步得到最优步长：
其中，α(0≤α＜1)为平滑因子，最小函数min(·)是为保证算法在冲激噪声干扰下的稳定性，即保证步长在迭代过程中一直减小。
由此得出变步长的仿射投影符号算法权系数的更新表达式为：
Figure GDA0002578613180000145
新算法的提出及分析
上述VSS-APSA算法虽然提高了APSA的收敛性能与稳态误差性能，但其对噪声项
Figure GDA0002578613180000146
进行了删除，导致其步长的计算并非最优的。同时，式(30)中最小函数min(·)的应用使得算法在非稳定环境下(如待辨识系统冲激响应在某时刻从w0变为-w0)性能严重衰减。
本申请引入步长函数，并将变步长的方法和比例矩阵的思想融合到一起，提出了一种新的算法，即变步长的改进比例仿射投影符号算法(VSS-IPAPSA)。该算法不仅降低了IPAPSA算法在非稀疏系统下稳态误差，而且在存在冲击噪声环境下稳健性和收敛速度都要优于传统的仿射投影算法，之后的理论分析和仿真结果验证了其有效性和稳健性。
本申请算法中成比例思想的引入
将PNLMS算法[11]中的成比例的方法引入APSA算法中，可得到PAPSA算法[16]的更新公式
Figure GDA0002578613180000151
其中，L×L对角矩阵G(n)表示成比例矩阵，gi(n)表示滤波器第i个抽头参数在n时刻的成比例比重。其对角矩阵G(n)的对角元素gi(n)计算方法如下:
G(n)＝diag{g0(n),…gL-1(n)}
Figure GDA0002578613180000152
Figure GDA0002578613180000153
Figure GDA0002578613180000154
其中,参数δp是一个很小的正数,它在所有滤波器抽头权值为零的时候启动更新；ρ和f(n)能够防止远小于最大抽头的权值停止更新。
同理，将IPNLMS算法[12]的思想应用于APSA算法中，得到本申请的IPAPSA算法，其与PAPSA算法的权系数更新公式相同，不同之处在于对角元素上：
Figure GDA0002578613180000155
式中:θ是与系统稀疏度(稀疏度越大表示回声路径越稀疏)有关的参数，ε是一个比较小的正数，以防止分母为0。
3.2本申请算法的最优步长推导及分析
(1)最优步长推导
由于式(29)中噪声项的删除，使得传统的滤波算法中所得步长不是最优的，因此，为了得到最优步长，本申请不作噪声向量与误差向量相关性的假设，即不删除噪声项，采用随机逼近的方法推导出步长，定义滤波器系数误差向量
Figure GDA0002578613180000161
并将固定步长μ替代为可变步长μ(n)。具体推导过程如下：
对式(26)两边取2范数平方的期望得：
Figure GDA0002578613180000162
其中，f(μ(n))是关于步长的一个函数，且
Figure GDA0002578613180000163
v(n)是噪声向量，符号
Figure GDA0002578613180000164
表示定义为。
为了使MSD的值达到最小化，本申请的改进算法从n到n+1的迭代过程中通过选择最优步长来使步长函数最小化，即步长函数写为：：
Figure GDA0002578613180000165
当自适应滤波器收敛到最佳状态时，
Figure GDA0002578613180000166
此时，sign(eT(n))≈sign(v(n))，且在稳定状态时：
Figure GDA0002578613180000167
由于准确地计算出sign(eT(n))v(n)的值比较困难，所以不能直接获得f(μ(n))；因此，本申请根据sign(eT(n))v(n)的上限，运用随机逼近方法得到：
Figure GDA0002578613180000168
但是，v(n)的绝对值不是精确的测量值，本申请将v(n)的绝对值近似为其期望值，v(n)有半正态分布的属性，且
Figure GDA0002578613180000169
即将||v(n)||1改为：
Figure GDA0002578613180000171
由式(39)得，
Figure GDA0002578613180000172
其中，M为投影阶数，τv为噪声信号v(n)的方差。
运用随机逼近的方法，步长函数的上限重写为：
Figure GDA0002578613180000173
从n到n+1的迭代过程中，最小化步长函数f(μ(n))大幅度降低了MSD的值。因此，对式(41)关于μ(n)求偏导，得，
Figure GDA0002578613180000174
令上式求导结果为0，得到的派生步长为：
Figure GDA0002578613180000175
(2)最优步长分析
由式(43)推导出的是步长最小化步长函数的上限。然而，由于期望项的存在，要得到准确的步长是非常困难的。需要考虑以下两种情况：
Figure GDA0002578613180000176
情况1：
Figure GDA0002578613180000177
即β(n)>0时，步长因子按本申请中的迭代公式进行更新；
情况2:当
Figure GDA0002578613180000181
即β(n)≤0时，此时μ(n)停止更新，即
μ(n)＝μ(n-1)。本申请采用移动平均法来处理上式中的期望项，得到如下：
Figure GDA0002578613180000182
其中，α(0≤α＜1)为平滑因子，本申请平滑因子的值设为
Figure GDA0002578613180000183
k为一个恒定的值(本申请中k取4)；
因此，本申请所提出的VSS-IPAPSA算法的滤波器系数更新公式为：
Figure GDA0002578613180000184
G(n)＝diag{g0(n),…gL-1(n)}
Figure GDA0002578613180000185
Figure GDA0002578613180000186
Figure GDA0002578613180000187
3.3算法总结
综上所述，本申请提出的算法总结如表2所示。
表2 VSS-IPAPSA算法总结
Figure GDA0002578613180000188
Figure GDA0002578613180000191
由于上述迭代过程中需要估算噪声信号v(n)的方差σv，考虑到冲击噪声干扰对噪声方差估算方法的影响，本申请通过使用中值滤波器来提高输出噪声方差对冲击噪声的稳健性。推导如下：
if mod(n,Nw)＝0
Figure GDA0002578613180000192
Figure GDA0002578613180000201
Figure GDA0002578613180000202
Figure GDA00025786131800002010
end (28)
其中，Ae(n)＝[e(n)e(n-1)…e(n-Nw+1)],Be(n)＝[e2(n)e2(n-1)…e2(n-Nw+1)]。
Figure GDA0002578613180000203
表示n时刻估计的误差方差的平方，
Figure GDA0002578613180000204
表示n时刻估计的输入信号方差的平方，
Figure GDA0002578613180000205
表示n时刻估计的噪声方差的平方，
Figure GDA0002578613180000206
表示n时刻所估计的中间向量，模mod(n,Nw)表示整数n和Nw之间的余数，Nw为估计窗的长度，
Figure GDA0002578613180000207
为有限样值校正因子，目的是减少输出噪声方差估计算法在每n次迭代时计算复杂性。
仿真实验
仿真条件
本申请在MATLAB环境下进行仿真实验。输入的有色信号均由零均值高斯白噪声通过一阶AR系统
Figure GDA0002578613180000208
产生，这里假设自适应滤波器的长度和未知系统的长度相等，均为120。实验中变步长算法的初始步长统一设置为μ(0)＝0.1，各个算法的投影阶数均为4，且算法的比较都在公平原则下进行。每个仿真均是20次实验的平均结果。信道h由式(47)产生:
Figure GDA0002578613180000209
其中，k表示自适应滤波器的权系数序号，本申请中k取1至120。
性能指标及参数设置
实验是使算法在等效步长的条件下，以权误差向量(WEVN，Weight Error Vector)和归一化均方偏差(NMSD，Normalized MSD)的收敛曲线来评价算法的收敛性能。
Figure GDA0002578613180000211
WEVN(权误差向量)的值愈小意味着自适应滤波器愈逼近未知系统。
Figure GDA0002578613180000212
NMSD(归一化的均方偏差)用于度量自适应滤波器与目标系统的逼近程度。
仿真中的参数是在等效步长[17]的条件下设定的，如表3所示。
表3实验中的参数设置
Figure GDA0002578613180000213
仿真结果与分析
(1)高斯噪声条件下各类算法在非稀疏系统中的性能比较
图5为APA、APSA、VSS-APSA和本申请改进的新算法(VSS-IPAPSA)在高斯噪声条件下的NMSD收敛曲线，各算法参数设置如表2所示。
由图5可以看到，步长参数对APA算法和APSA算法收敛速度有明显的影响，VSS-APSA算法利用变步长的方法使得其在保证收敛速度加快的同时，稳态误差降低。本申请提出的VSS-IPAPSA算法结合了变步长方法和比例矩阵的优势，既具有快的收敛速度又具有小的稳态误差。
(2)非高斯噪声条件下各类算法在非稀疏系统中的性能比较
图6为APA、APSA、VSS-APSA、IPAPSA和本申请改进的新算法(VSS-IPAPSA)在非高斯噪声条件下的NMSD收敛曲线，各算法参数设置如表3所示。
由图6可以看到，当干扰噪声为非高斯噪声时，APA算法失效，APSA算法的收敛速度变慢，VSS-APSA算法的收敛速度要优于APSA算法，IPAPSA算法的稳健性和收敛速度要优于APSA算法和其VSS-APSA算法，但是IPAPSA算法的稳态性能比APSA和VSS-APSA差，而本申请提出的新的算法，即变步长IPAPSA算法(VSS-IPAPSA)，在保证其他性能优于传统算法的同时稳态性能得到了大大的提高。
(3)非高斯噪声条件下各类算法在稀疏系统中的性能比较
图7为APSA、VSS-APSA、IPAPSA和本申请改进的新算法(VSS-IPAPSA)在稀疏系统且存在非高斯噪声条件下的NMSD收敛曲线。
由图7可以看到，在稀疏系统下IPAPSA算法的性能要优于APSA算法和其VSS-APSA算法，本申请提出的VSS-IPAPSA算法结合了变步长算法和比例矩阵对收敛速度的优势，使其收敛速度、稳健性及稳态误差都要优于传统的仿射投影算法。
(4)非高斯噪声条件下各类算法的跟踪性能比较
图8中(a)和(b)为非高斯噪声条件下的APSA、VSS-APSA、IPAPSA和新算法(VSS-IPAPSA)分别在稀疏系统与非稀疏系统下的跟踪性能比较图。当迭代次数n≤3000时，所使用的信道是稀疏度为0.85349的稀疏信道(见图4(a))，当迭代次数n≥3000时，所使用的信道是稀疏度为0.32883的非稀疏信道(见图4(b))。
由图8可以看到，本申请改进的新算法(VSS-IPAPSA)在系统突变时的跟踪性能稳定，且优于其他算法。
(5)不同信噪比下，各类算法的权误差向量
在以上仿真实验中，测量噪声v(n)加到输出信号(XT(n)w0)上，其信噪比为(SNR)25dB。SNR的定义如下：
Figure GDA0002578613180000231
然而，不同的信噪比对自适应算法的影响也是不一样的，表4为各类算法在不同信噪比下的权误差向量对比表，其中，其他实验条件同实验3一致。
表4各类算法在不同信噪比下的权误差向量(WEVN)
Figure GDA0002578613180000232
由表4可得，通过对比，发现本申请改进的新算法(VSS-IPAPSA)在不同信噪比条件下的稳态性能要优于同一条件下的其他算法。
综上：
本申请提出的基于l1范数的改进的变步长比例仿射投影符号自适应滤波算法(VSS-IPAPSA)是针对非高斯噪声干扰对自适应滤波算法性能的影响和系统的不同稀疏性所提出的。具体地，本申请主要进行了三方面的工作：1)引入步长函数的概念，采用随机逼近的方法推导得最优步长；2)在实际情况分析过程中，采用移动平均法处理期望项，得出准确的步长；3)将比例矩阵的思想引入传统的仿射投影算法中，提高算法在不同稀疏系统中的适用性，同时利用符号算法对非高斯噪声干扰的抑制能力。实验结果和理论分析表明，本申请的算法对不同稀疏特性的系统具有更好的适应性，且对非高斯噪声的抑制能力、收敛速度、稳态性能、系统发生突变时的跟踪性能均优于其他的自适应滤波算法。
本申请将IPNLMS中比例系数的思想应用于APSA，得到了改进的比例仿射投影符号算法(IPAPSA，Improved Proportionate APSA)，同时结合变步长的方法，提出了一种变步长的改进比例仿射投影符号算法(VSS-IPAPSA，Variable Step-Size IPAPSA)，本算法一方面引入步长函数，采用随机逼近和移动平均的方法得到算法收敛的最优步长，另一方面结合比例矩阵的优势，使得算法对不同系统稀疏性具有广泛的适用性，仿真结果和理论分析表明本申请提出的算法在非高斯噪声和高斯噪声干扰下均具有良好的收敛性能、稳态性能和跟踪性能。
上述的描述中：
[11]单志明.α-稳定分布参数估计及自适应算法研究[D].哈尔滨：哈尔滨工程大学，2012.
[12]张凯，单甘霖.基于高斯和SCKF的非线性非高斯滤波算法[J].仪器仪表学报，2015,35(11)：2524-2530.
[13]SAMORODNITSKY G，TAQQU M S.Stable non-Gaussian random processes：stochastic models with infinite variance[D].New York：Chapman and Hall，1994.
[14]DUTTWEILER D L.Proportionate normalized least-mean-squaresadaptation in echo cancellers[J].IEEE Transactions on Speech and AudioProcessing，2000，8(5)：508-518.
[15]ALBU FELIX，KWAN H.K.New proportionate affine projection signalgorithms[J].IEEE International Symposium on Circuits and Systems，2013：521-524.
[16]BENESTY J，GAY S L.An improved PNLMS algorithm[C].IEEEInternational Conference On Acoustics，Speech，and Signal Processing，Orlando，FL，USA，2002：1881-1884.
[17]郭莹，高媛.几种稀疏自适应滤波算法的对比研究[J].微型机与应用，2014，33(8)：1-3。

Claims (3)

1.一种稳健的仿射投影符号自适应滤波方法，其特征在于：引入步长函数，并将变步长的方法和比例矩阵的思想融合到一起，即变步长的改进比例仿射投影符号算法VSS-IPAPSA；包括三方面的工作：1)引入步长函数的概念，采用随机逼近的方法推导得最优步长；2)在实际情况分析过程中，采用移动平均法处理期望项，得出准确的可变步长更新公式；3)将比例矩阵的思想引入传统的仿射投影算法中，提高算法在不同稀疏系统中的适用性，同时利用符号算法对非高斯噪声干扰的抑制能力；
该VSS-IPAPSA算法中将比例归一化最小均方算法PNLMS中成比例的方法引入仿射投影符号算法APSA算法中，得到比例仿射投影符号算法PAPSA的更新公式：
Figure FDA0002578613170000011
其中，sgn(·)表示符号运算，μ是步长因子，L是滤波器长度，n是时间系数，
Figure FDA0002578613170000012
表示所估计的自适应滤波器权系数向量，用于辨识未知系统w0，输入信号向量为x(n)＝[x(n) x(n-1) ... x(n-L+1)]T，未知系统的期望响应信号为d(n)＝xT(n)w0+v(n)，v(n)是噪声，则先验输出误差向量和后验输出误差向量分别为
Figure FDA0002578613170000013
和
Figure FDA0002578613170000014
d(n)＝[d(n) d(n-1) ... d(n-M+1)]T是期望输出信号向量，M为投影阶数，X(n)是APA滤波器结构中包含M维输入向量的输入矩阵，即X(n)＝[x(n) x(n-1) ... x(n-M+1)]，L×L对角矩阵G(n)表示成比例矩阵，gi(n)是其对角矩阵G(n)的对角元素，表示滤波器第i个抽头参数在n时刻的成比例比重，G(n)和gi(n)的计算方法如下：
G(n)＝diag{g0(n)，...gL-1(n)} (2)
Figure FDA0002578613170000021
Figure FDA0002578613170000022
Figure FDA0002578613170000023
其中，max(·)表示求最大值，||·||∞表示求无穷范数，参数δp是一个很小的正数，它在所有滤波器抽头权值为零的时候启动更新；ρ和f(n)能够防止远小于最大抽头的权值停止更新；
同理，将IPNLMS算法的思想应用于APSA算法中，得到IPAPSA算法，其与PAPSA算法的权系数更新公式相同，不同之处在于对角元素上，
Figure FDA0002578613170000024
式中：||·||1表示求1范数，θ是与系统稀疏度有关的参数，稀疏度越大表示回声路径越稀疏，ε是一个比较小的正数，以防止分母为0。
2.根据权利要求1所述的稳健的仿射投影符号自适应滤波方法，其特征在于：
(1)最优步长推导
基于最小均方偏差MSD；定义滤波器系数误差向量
Figure FDA0002578613170000025
并将固定步长μ替代为可变步长μ(n)；为了得到最优步长，不删除噪声项，采用随机逼近的方法推导出步长，具体推导过程如下：
对式
Figure FDA0002578613170000026
两边取2范数平方的期望得：
Figure FDA0002578613170000027
其中，||·||2表示求2范数；f(μ(n))是关于步长的一个函数，且
Figure FDA0002578613170000031
Figure FDA0002578613170000032
v(n)是噪声向量，符号
Figure FDA0002578613170000033
表示定义为；
为了使MSD的值达到最小化，从n到n+1的迭代过程中通过选择最优步长来使步长函数最小化，即步长函数写为：
Figure FDA0002578613170000034
当自适应滤波器收敛到最佳状态时，
Figure FDA0002578613170000035
此时，sgn(eT(n))≈sgn(v(n))，且在稳定状态时：
Figure FDA0002578613170000036
根据sgn(eT(n))v(n)的上限，运用随机逼近方法得到下式：
Figure FDA0002578613170000037
将v(n)的绝对值近似为其期望值，v(n)有半正态分布的属性，且
Figure FDA0002578613170000038
即将||v(n)||1改为：
Figure FDA0002578613170000039
由式(12)得，
Figure FDA00025786131700000310
其中，M为投影阶数，σv为噪声信号v(n)的方差；
运用随机逼近的方法，步长函数的上限重写为：
Figure FDA0002578613170000041
从n到n+1的迭代过程中，最小化步长函数f(μ(n))大幅度降低了MSD的值；因此，对式(14)关于μ(n)求偏导，得，
Figure FDA0002578613170000042
令上式求导结果为0，得到的派生步长为：
Figure FDA0002578613170000043
(2)最优步长分析
通过式(16)推导出了步长最小化步长函数的上限，考虑以下两种情况，首先定义：
Figure FDA0002578613170000044
情况1：
Figure FDA0002578613170000045
即β(n)＞0时，步长因子按迭代公式进行更新；
情况2：当
Figure FDA0002578613170000046
即β(n)≤0时，此时μ(n)停止更新，即μ(n)＝μ(n-1)；
采用移动平均法来处理上式中的期望项，得到如下：
Figure FDA0002578613170000047
，
k为一个恒定的值；
因此，VSS-IPAPSA算法的滤波器系数更新公式为：
Figure FDA0002578613170000051
G(n)＝diag{g0(n)，...gL-1(n)}
Figure FDA0002578613170000052
Figure FDA0002578613170000053
Figure FDA0002578613170000054
3.根据权利要求2所述的稳健的仿射投影符号自适应滤波方法，其特征在于：
VSS-IPAPSA算法为：
初始化：
Figure FDA0002578613170000055
μ(0)＝0.1
参数设置：θ＝0 or-0.5，ε＝0.01，
Figure FDA0002578613170000056
Nw＝8，
Figure FDA0002578613170000057
迭代过程：
Figure FDA0002578613170000058
G(n)＝diag{g0(n)，...gL-1(n)}
Figure FDA0002578613170000059
Ae(n)＝[e(n) e(n-1) ... e(n-Nw+1)]
Be(n)＝[e2(n) e2(n-1) ... e2(n-Nw+1)]
if mod(n，Nw)＝0
Figure FDA00025786131700000510
Figure FDA0002578613170000061
Figure FDA0002578613170000062
Figure FDA0002578613170000063
end
Figure FDA0002578613170000064
if β(n)＞0
μ(n)＝αμ(n-1)+(1-α)min(β(n)，μ(n-1))
else
μ(n)＝μ(n-1)
end
Figure FDA0002578613170000065
end
由于上述迭代过程中需要估算噪声信号v(n)的方差σv，通过使用中值滤波器来提高输出噪声方差对冲击噪声的稳健性，过程如下：
if mod(n，Nw)＝0
Figure FDA0002578613170000066
Figure FDA0002578613170000067
Figure FDA0002578613170000071
Figure FDA0002578613170000072
end
其中，Ae(n)＝[e(n) e(n-1) ... e(n-Nw+1)]，Be(n)＝[e2(n) e2(n-1) ... e2(n-Nw+1)]；
Figure FDA0002578613170000073
表示n时刻估计的误差方差的平方，
Figure FDA0002578613170000074
表示n时刻估计的输入信号方差的平方，
Figure FDA0002578613170000075
表示n时刻估计的噪声方差的平方，
Figure FDA0002578613170000076
表示n时刻所估计的中间向量，模mod(n，Nw)表示整数n和Nw之间的余数，Nw为估计窗的长度，
Figure FDA0002578613170000077
为有限样值校正因子，为了减少所提出的输出噪声方差估计算法在每n次迭代时计算复杂性。
CN201611059496.5A 2016-10-11 2016-11-28 稳健的仿射投影符号自适应滤波算法 Active CN106788337B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610886005 2016-10-11
CN2016108860058 2016-10-11

Publications (2)

Publication Number Publication Date
CN106788337A CN106788337A (zh) 2017-05-31
CN106788337B true CN106788337B (zh) 2020-08-28

Family

ID=58911544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611059496.5A Active CN106788337B (zh) 2016-10-11 2016-11-28 稳健的仿射投影符号自适应滤波算法

Country Status (1)

Country Link
CN (1) CN106788337B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107452393B (zh) * 2017-07-27 2020-09-11 桂林电子科技大学 一种回声消除系统及方法
CN107800403B (zh) * 2017-09-14 2021-04-23 苏州大学 一种鲁棒样条自适应滤波器
CN108039179B (zh) * 2017-12-13 2021-08-13 重庆邮电大学 一种用于麦克风阵列广义旁瓣消除器的高效自适应算法
CN107947761B (zh) * 2017-12-18 2021-09-10 西安理工大学 基于最小均方四阶的变阈值比例更新自适应滤波方法
CN108279564B (zh) * 2018-01-30 2020-11-03 苏州大学 一种鲁棒的稀疏多任务自适应系统及迭代方法
CN108512529B (zh) * 2018-03-05 2022-10-25 成都优艾维智能科技有限责任公司 一种μ律函数下的比例控制和归一化LMP滤波方法
CN109151237B (zh) * 2018-08-23 2020-10-09 西南交通大学 基于零吸引的仿射投影自适应回声消除方法
CN109767779A (zh) * 2018-11-17 2019-05-17 沈阳工业大学 基于最小误差熵的成比例仿射投影方法
CN110211071A (zh) * 2019-06-05 2019-09-06 陕西师范大学 一种基于L1/2范数的高分三号SAR图像的Frost滤波算法
CN110133425A (zh) * 2019-06-10 2019-08-16 集美大学 一种海底电缆故障信号滤波方法、终端设备及存储介质
CN110890875B (zh) * 2019-09-29 2021-08-10 电子科技大学 基于广义相关诱导度量改进的比例仿射投影滤波方法
CN113381731B (zh) * 2020-06-22 2023-10-20 沈阳工业大学 一种针对非高斯噪声的扩散式变步长自适应参数估计方法
CN112003588B (zh) * 2020-08-27 2023-11-17 重庆邮电大学 基于多态可变步长归一化均方的自适应信号滤波方法
CN113225045B (zh) * 2021-03-25 2023-06-23 苏州大学 低计算复杂度的稀疏促进仿射投影自适应滤波器
CN113315490A (zh) * 2021-06-04 2021-08-27 哈尔滨工程大学 一种基于后滤波器的滤波参考仿射投影符号主动控制系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101432426B1 (ko) * 2013-02-14 2014-08-22 한양대학교 산학협력단 적응필터의 필터계수 결정 방법 및 장치
CN104683614B (zh) * 2015-03-24 2016-03-02 西南交通大学 基于m估计的记忆成比例仿射投影自适应回声消除方法
CN105721729B (zh) * 2016-04-26 2019-01-15 西南交通大学 基于块稀疏成比例重用权系数仿射投影的回声消除方法

Also Published As

Publication number Publication date
CN106788337A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106788337B (zh) 稳健的仿射投影符号自适应滤波算法
CN109040499B (zh) 一种抗冲击干扰的自适应回声消除方法
CN108200522B (zh) 一种变正则化比例归一化子带自适应滤波方法
Albu Improved variable forgetting factor recursive least square algorithm
Resende et al. A variable step-size NLMS algorithm with adaptive coefficient vector reusing
Yazdanpanah et al. A simple set-membership affine projection algorithm for sparse system modeling
CN109102794A (zh) 基于凸组合的m估计成比例类仿射投影的回声消除方法
Yang et al. Efficient μ-law improved proportionate affine projection algorithm for echo cancellation
Albu et al. New variable step size affine projection algorithms
Shi et al. A frequency domain step-size control method for LMS algorithms
Gil-Cacho et al. Nonlinear acoustic echo cancellation based on a parallel-cascade kernel affine projection algorithm
CN109767779A (zh) 基于最小误差熵的成比例仿射投影方法
Yu et al. A band-independent variable step size proportionate normalized subband adaptive filter algorithm
Mohanaprasad et al. Wavelet-based ICA using maximum likelihood estimation and information-theoretic measure for acoustic echo cancellation during double talk situation
Arezki et al. A new algorithm with low complexity for adaptive filtering
Gao et al. Exponentially weighted kernel recursive least p-power algorithm
Loganathan et al. A sparseness controlled proportionate algorithm for acoustic echo cancellation
CN112803919B (zh) 改进nlms算法的稀疏系统辨识方法和滤波器和系统
Werner et al. A set-membership approach to normalized proportionate adaptation algorithms
Hamidia et al. Double-talk detection using the singular value decomposition for acoustic echo cancellation
Mohammed et al. An efficient RLS algorithm for output-error adaptive IIR filtering and its application to acoustic echo cancellation
GIRIKA et al. ADAPTIVE SPEECH ENHANCEMENT TECHNIQUES FOR COMPUTER BASED SPEAKER RECOGNITION.
Turan et al. A block LMS-type algorithm with a function controlled variable step-size for sparse system identification
Albu et al. Intermittently-updated affine projection algorithm
Yu et al. General Robust Subband Adaptive Filtering for Echo Cancellation

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