CN102724152A - 基于Laguerre结构的多项式自适应有源噪声对消方法 - Google Patents

基于Laguerre结构的多项式自适应有源噪声对消方法 Download PDF

Info

Publication number
CN102724152A
CN102724152A CN2012102403347A CN201210240334A CN102724152A CN 102724152 A CN102724152 A CN 102724152A CN 2012102403347 A CN2012102403347 A CN 2012102403347A CN 201210240334 A CN201210240334 A CN 201210240334A CN 102724152 A CN102724152 A CN 102724152A
Authority
CN
China
Prior art keywords
mrow
msub
msup
math
filter
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
CN2012102403347A
Other languages
English (en)
Other versions
CN102724152B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201210240334.7A priority Critical patent/CN102724152B/zh
Publication of CN102724152A publication Critical patent/CN102724152A/zh
Application granted granted Critical
Publication of CN102724152B publication Critical patent/CN102724152B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Filters That Use Time-Delay Elements (AREA)

Abstract

本发明公开了一种基于Laguerre结构的多项式自适应有源噪声对消方法。它将具有IIR滤波器特点的Laguerre时延单元取代Volterra滤波器中的延迟单元Z-1,经过Laguerre延迟节点的输出作为基于有界约束乘积耦合Volterra近似结构的非线性滤波器的输入,能有效地改善特征收敛条件,保证系统的稳定性,计算复杂度低,硬件容易实现;并采用输入信号和瞬时误差归一化的LMS自适应算法调整其滤波系数,能够在保证非线性自适应滤波器稳定的前提下,进一步降低滤波器的输入长度,加快收敛速度;特别对于存在较强的非线性畸变情形,显示出很好的抗噪能力,有利于自适应多项式滤波器的推广与应用。

Description

基于Laguerre结构的多项式自适应有源噪声对消方法
技术领域
本发明涉及一种基于Laguerre结构的多项式自适应有源噪声对消方法,用于噪声、码间串扰以及非线性畸变的消除,具有良好的抗噪声性能。
背景技术
近年来,随着现代信号处理理论的发展,各种非线性技术已经成功的广泛应用于存在着非平稳信号、非高斯噪声及非线性干扰的扩频通信、系统辨识、噪声控制、图像处理和混沌预测等非线性信号处理领域,因此,非线性信号处理技术受到国内外信号处理研究人员的越来越多的关注,同时也是现代信号处理理论中一个重要研究方向。研究表明:在非线性、非高斯或高斯背景下的最优滤波是非线性的。非线性自适应滤波器常具有比线性自适应滤波器更好的性能,并且由于其本质上的非线性特性,能够克服线性自适应滤波器在探究非线性信号的高阶统计冗于性方面能力和逼近非线性函数能力的限制,显示出其良好的非线性处理能力。现有的非线性滤波技术中,常见的如下两种方法:
(1)基于神经网络的非线性自适应滤波器。
由于神经网络能够以任意精度逼近非线性函数,并且具有自适应能力,因此,各种神经网络已经成功地应用在自适应有源噪声对消系统中。参考文献[1]“Recurrent radial basis function networks for adaptive noise cancellation”(S.A.Billings,C.F.Fung,Neural Networks[J],vol.8,No.2,pp.273-290,1995)利用径向基神经网络的反馈结构构造了递归径向基神经网络自适应对消器,通过基于混合K-均值聚类的梯度算法对其进行训练。参考文献[2]“神经网络自适应有源噪声对消器”(莫玮,蒋洪睿,谢维信,信号处理[J],Vol.16,No.3,pp.286-289,2000)研究了RNN自适应滤波器组成的有源噪声对消器,采用具有类IIR属性的递归神经网络(RNN),并通过RTRL算法对滤波系数进行调整,因此,在对消系统中只需较少的输入长度即可达到较好的抗噪效果。但是,由于神经网络本身故有的结构复杂、参数空间过大以及参数的初始化需谨慎等缺点,限制了其在有源噪声对消系统中实用化;而且由于IIR结构本身固有不稳定的缺陷,致使系统的稳定性难以保证。
(2)基于多项式的非线性自适应滤波器。
Volterra滤波器是非线性自适应滤波理论中一个重要的研究方向,由于其输出线性依赖于滤波器本身系数,从而使Volterra滤波器与其核矢量具有线性关系,便于分析线性系统;同时,Volterra滤波器存在大量的非线性耦合项,因此使基于截断Volterra模型的非线性技术能够有效的处理非线性噪声问题。然而,随着Volterra滤波器的阶数或记忆单元m增大,致使其权系数按幂次快速增加,相应所需的计算次数亦成幂次快速增加,这将限制其在有源噪声对消系统中实用化应用。参考文献[3]“基于Sigmoid函数的Volterra自适应有源噪声对消器”(张家树,肖先赐,电子与信息学报,Vol.24,No.4,pp.461-466,2002)采用Volterra的乘积耦合近似结构,构造了一种低复杂度的非线性自适应对消器。参考文献[4]“一种改进的Volterra自适应噪声对消器”(张秀梅,赵知劲,尚俊娜,电讯技术,Vol.34,No.2,pp.77-81,2010)利用Sigmoid函数对背景噪声进行预处理,改进了上述Volterra的近似结构,并用它构造了非线性自适应噪声对消器,具有较好的抗噪性能,但是其系统的稳定性不能够保证,而且为了得到较好的性能,增加其输入信号的长度,然而这将导致计算复杂度的增加。
综上所述,由于存在计算复杂及系统稳定性不能够保证等缺点,将限制上述的非线性自适应滤波器在有源噪声对消系统中实用化应用。
发明内容
本发明的目的就是提供一种基于Laguerre结构的多项式白适应有源噪声对消方法,该方法稳定性好、收敛速度快,计算复杂度低,硬件实现容易。
本发明实现其发明目的,所采用的技术方案是,一种基于Laguerre结构的多项式白适应有源噪声对消方法,其步骤如下:
(1):参考噪声源的产生
将取值在[-1,1]区间且方差σ2=3的白噪声,作为对含有用信号s(n)和噪声v(n)的接收信号d(n)进行噪声对消用的参考噪声源r(n),其中n为时刻;
(2):基于Laguerre延迟节点的输出信号矢量的计算
将参考噪声源r(n)通过Laguerre延迟节点得到m个延迟输出信号分量uk(n),0≤k≤m-1;
当k=0时, u 0 ( n ) = α u 0 ( n - 1 ) + 1 - α 2 r ( n - 1 ) ;
当k=1,2,3,...,k-1时,uk(n)=uk-1(n-1)+α[uk(n-1)-uk-1(n)]
其中,k为Laguerre延迟节点的序号,α为Laguerre延迟节点的极点;
由m个延迟输出信号分量uk(n)及常数1构成延迟输出信号向量U(n),
U(n)=[1,u0(n),u1(n),…,um(n)]T
(3):噪声的估计
将第(2)步得到的延迟输出信号向量U(n)输入有界约束乘积耦合Volterra近似结构的滤波器,滤波得到接收信号d(n)中的噪声v(n)的估计y(n),
y ( n ) = a 0 ( n ) + Σ k = 0 m - 1 a k + 1 ( n ) u k ( n ) + [ b 0 ( n ) + Σ k = 0 m - 1 b k + 1 ( n ) u k ( n ) ] × Sigm ( c 0 ( n ) + Σ k = 0 m - 1 c k + 1 ( n ) u k ( n ) )
= [ a 0 ( n ) , a 1 ( n ) , · · · , a m ( n ) , b 0 ( n ) , b 1 ( n ) , · · · , b m ( n ) ] [ U ( n ) , U ( n ) Sigm ( C T ( n ) U ( n ) ) ] T
= AB T ( n ) X ( n )
其中,Sigm(·)为有界函数,其表达式为Sigm(·)=(1-exp(·))/(1+exp(·))且|Sigm(·)|<1,X(n)=[U(n),U(n)Sigm(CT(n)U(n))]T是滤波器的输入信号向量,其对应的滤波器权系数为AB(n)=[a0(n),a1(n),…,am(n),b0(n),b1(n),…,bm(n)]T,而C(n)=[c0(n),c1(n),…,cm(n)]T为n时刻滤波器有界约束分解子项的权系数向量;
(4):有用信号的估计
将接收信号d(n)减去第(3)步得到的噪声v(n)的估计y(n),得到n时刻接收信号d(n)中的有用信号s(n)的估计e(n):
e(n)=d(n)-y(n)=s(n)+v(n)-y(n)
(5):权系数更新
将第(4)步得到的有用信号估计e(n)的平方,求期望得到代价函数J(n),J(n)=E{e2(n)}=E{|d(n)-y(n)|2},E{·}为求期望,再用代价函数J(n)对非线性滤波器的权系数AB(n)和C(n)求瞬时梯度,并利用输入信号和瞬时误差归一化的LMS算法准则,得到权系数AB(n)和C(n)的如下迭代公式:
σ X 2 ( n ) = ρ σ X 2 ( n - 1 ) + ( 1 - ρ ) [ | | X ( n ) | | 2 + ( 2 m + 2 ) e 2 ( n ) ]
AB ( n + 1 ) = AB ( n ) + η 1 e ( n ) η 1 σ X 2 ( n ) X ( n )
σ U 2 ( n ) = ρ σ U 2 ( n - 1 ) + ( 1 - ρ ) ( | | U ( n ) | | 2 + ( m + 1 ) e 2 ( n ) [ b ( n ) T U ( n ) ( 1 - Sigm 2 ( C T ( n ) U ( n ) ) ) ] 2 )
C ( n + 1 ) = C ( n ) + η 2 e ( n ) η 2 σ U 2 ( n ) b ( n ) T U ( n ) ( 1 - Sigm 2 ( C T ( n ) U ( n ) ) )
其中
Figure BDA00001879139100048
Figure BDA00001879139100049
分别为滤波器输入信号向量X(n)和U(n)的功率;b(n)=[b0(n),b1(n),…,bm(n)]T;η1为控制滤波器权系数向量AB(n)收敛的步长参数,取值区间为(0,2);η2为控制滤波器有界约束分解子项的权系数向量C(n)收敛的步长参数,取值区间为(0,2);ρ为平滑因子其取值区间为(0 0.1];
(6):令n=n+1,重复2~5的步骤,即可实时消除接受信号d(n)中的污染噪声v(n),得到实时的有用信号s(n)。
与现有技术相比,本发明的有益效果是:
一、改善特征收敛条件,滤波器的输入长度少,计算复杂度低,硬件实现容易。
将参考噪声源r(n)通过Laguerre延迟节点得到m+1个延迟输出信号分量uk(n),0≤k≤m;
当k=0时:Laguerre延时节点的输出可以通过以下公式计算得到:
u 0 ( n ) = α u 0 ( n - 1 ) + 1 - α 2 r ( n - 1 ) - - - ( 1 )
该延时节点为单极点α的低通滤波,且转移函数记为L0(z,a),即
L 0 ( z , a ) 1 - a 2 1 - a z - 1 , ( 0 < &alpha; < 1 ) - - - ( 2 )
当k=1,2,3,...,m时:Laguerre延时节点的输出可以通过以下公式计算得到:
uk(n)=uk-1(n-1)+α[uk(n-1)-uk-1(n)]       (3)
该延时节点为全通滤波L(z,a),其转移函数记为:
L ( z , a ) z - 1 - a 1 - a z - 1 , ( 0 < &alpha; < 1 ) - - - ( 4 )
则第k个Laguerre延迟节点的抽头输出处的转移函数为
L k ( z , &alpha; ) = 1 - &alpha; 2 ( z - 1 - &alpha; ) k ( 1 - &alpha; z - 1 ) k + 1 , k = 0,1 , &CenterDot; &CenterDot; &CenterDot; , m - - - ( 5 )
Laguerre序列lk(n,α)(即Lk(z,α)的逆Z变换)的表达式为:
l k ( n , &alpha; ) = 1 - &alpha; 2 &Sigma; j = 0 k ( - 1 ) k + j &zeta; j k &zeta; k n + k - j &alpha; n + k - 2 j - - - ( 6 )
其中,
Figure BDA00001879139100062
因该Laguerre序列lk(n,α)满足
&Sigma; n = 0 &infin; l k ( n , &alpha; ) l j ( n , &alpha; ) = &delta; k , j = 1 , k = j 0 , k &NotEqual; j - - - ( 7 )
并且(5)式又可以写为:Lk(z)=L0(z)Lk(z),则该离散Laguerre函数构成了空间的一个完备的标准正交集。因此,将参考噪声源r(n)通过Laguerre延迟节点得到m+1个延迟输出信号分量uk(n)作为有界约束乘积耦合Volterra近似结构的滤波器的输入信号,可以有效的改善特征收敛条件,从而加快滤波器的收敛速度,而且由于Laguerre延迟节点具有IIR特征,从而可以有效地减少滤波器的输入长度,降低滤波器的计算复杂度。
二、系统的稳定性好
本发明滤波器的输出y(n)可以进一步写成如下表达式:
y ( n ) = AB T ( n ) X ( n )
= AB T ( n ) [ U ( n ) , U ( n ) Sigm ( C T U ( n ) ) ] T
= a 0 ( n ) + &Sigma; k = 0 m - 1 a k + 1 ( n ) u k ( n ) + [ b 0 ( n ) + &Sigma; k = 0 m - 1 b k + 1 ( n ) u k ( n ) ] &times; Sigm ( c 0 ( n ) + &Sigma; k = 0 m - 1 c k + 1 ( n ) u k ( n ) ) - - - ( 8 )
= a 0 ( n ) + b 0 ( n ) + &Sigma; k = 1 m ( a k ( n ) + Sigm ( C T ( n ) U ( n ) ) b k ( n ) ) U ( n - m + k )
= a 0 ( n ) + b 0 ( n ) + g ( n )
其中, g ( n ) = &Sigma; k = 1 m ( a k ( n ) + Sigm ( C T ( n ) U ( n ) ) b k ( n ) ) U ( n - m + k ) .
要使整个系统稳定,则只需y(n)中的g(n)项稳定。对g(n)两边取Z变换得到:
G ( Z ) = &Sigma; k = 1 m ( a k ( n ) + Sigm ( C T ( n ) U ( n ) ) b k ( n ) ) L k ( Z ) X ( Z ) - - - ( 9 )
因此,传递函数H(Z)为:
H ( Z ) = G ( Z ) X ( Z )
= &Sigma; k = 1 m ( a k ( n ) + Sigm ( C T ( n ) U ( n ) ) b k ( n ) ) L k ( Z ) - - - ( 10 )
= &Sigma; k = 1 m ( a k ( n ) + Sigm ( C T ( n ) U ( n ) ) b k ( n ) ) L 0 ( Z ) L k ( Z )
由于当0<α<1时,(8)式满足Laguerre滤波稳定理论,因此传递函数H(Z)稳定,从而保证整个滤波器系统稳定。
下面结合附图和具体实施方式对本发明作进一步说明。
附图说明
图1a是本发明方法与文献[4]方法的仿真实验一的收敛曲线比较图。
图1b本发明方法的仿真实验一输出的有用信号估计。
图2a是本发明方法与文献[4]方法的仿真实验二的收敛曲线比较图。
图2b本发明方法的仿真实验二输出的有用信号估计。
具体实施方式
实施例
本发明的一种具体实施方式是,一种基于Laguerre结构的多项式自适应有源噪声对消方法,其步骤如下:
(1):参考噪声源的产生
将取值在[-1,1]区间且方差σ2=3的白噪声,作为对含有用信号s(n)和噪声v(n)的接收信号d(n)进行噪声对消用的参考噪声源r(n),其中n为时刻;
(2):基于Laguerre延迟节点的输出信号矢量的计算
将参考噪声源r(n)通过Laguerre延迟节点得到m个延迟输出信号分量uk(n),0≤k≤m-1;
当k=0时, u 0 ( n ) = &alpha; u 0 ( n - 1 ) + 1 - &alpha; 2 r ( n - 1 ) ;
当k=1,2,3,...,m-1时,uk(n)=uk-1(n-1)+α[uk(n-1)-uk-1(n)]
其中,k为Laguerre延迟节点的序号,α为Laguerre延迟节点的极点;
由m个延迟输出信号分量uk(n)及常数1构成延迟输出信号向量U(n),
U(n)=[1,u0(n),u1(n),…,um(n)]T
(3):噪声的估计
将第(2)步得到的延迟输出信号向量U(n)输入有界约束乘积耦合Volterra近似结构的滤波器,滤波得到接收信号d(n)中的噪声v(n)的估计y(n),
y ( n ) = a 0 ( n ) + &Sigma; k = 0 m - 1 a k + 1 ( n ) u k ( n ) + [ b 0 ( n ) + &Sigma; k = 0 m - 1 b k + 1 ( n ) u k ( n ) ] &times; Sigm ( c 0 ( n ) + &Sigma; k = 0 m - 1 c k + 1 ( n ) u k ( n ) )
= [ a 0 ( n ) , a 1 ( n ) , &CenterDot; &CenterDot; &CenterDot; , a m ( n ) , b 0 ( n ) , b 1 ( n ) , &CenterDot; &CenterDot; &CenterDot; , b m ( n ) ] [ U ( n ) , U ( n ) Sigm ( C T ( n ) U ( n ) ) ] T
= AB T ( n ) X ( n )
其中,Sigm(·)为有界函数,其表达式为Sigm(·)=(1-exp(·))/(1+exp(·))且|Sigm(·)|<1,X(n)=[U(n),U(n)Sigm(CT(n)U(n))]T是滤波器的输入信号向量,其对应的滤波器权系数为AB(n)=[a0(n),a1(n),…,am(n),b0(n),b1(n),…,bm(n)]T,而C(n)=[c0(n),c1(n),…,cm(n)]T为n时刻滤波器有界约束分解子项的权系数向量;
(4):有用信号的估计
将接收信号d(n)减去第(3)步得到的噪声v(n)的估计y(n),得到n时刻接收信号d(n)中的有用信号s(n)的估计e(n):
e(n)=d(n)-y(n)=s(n)+v(n)-y(n)
(5):权系数更新
将第(4)步得到的有用信号估计e(n)的平方,求期望得到代价函数J(n),J(n)=E{e2(n)}=E{|d(n)-y(n)|2},E{·}为求期望,再用代价函数J(n)对非线性滤波器的权系数AB(n)和C(n)求瞬时梯度,并利用输入信号和瞬时误差归一化的LMS算法准则,得到权系数AB(n)和C(n)的如下迭代公式:
&sigma; X 2 ( n ) = &rho; &sigma; X 2 ( n - 1 ) + ( 1 - &rho; ) [ | | X ( n ) | | 2 + ( 2 m + 2 ) e 2 ( n ) ]
AB ( n + 1 ) = AB ( n ) + &eta; 1 e ( n ) &eta; 1 &sigma; X 2 ( n ) X ( n )
&sigma; U 2 ( n ) = &rho; &sigma; U 2 ( n - 1 ) + ( 1 - &rho; ) ( | | U ( n ) | | 2 + ( m + 1 ) e 2 ( n ) [ b ( n ) T U ( n ) ( 1 - Sigm 2 ( C T ( n ) U ( n ) ) ) ] 2 )
C ( n + 1 ) = C ( n ) + &eta; 2 e ( n ) &eta; 2 &sigma; U 2 ( n ) b ( n ) T U ( n ) ( 1 - Sigm 2 ( C T ( n ) U ( n ) ) )
其中
Figure BDA00001879139100095
Figure BDA00001879139100096
分别为滤波器输入信号向量X(n)和U(n)的功率;b(n)=[b0(n),b1(n),…,bm(n)]T;η1为控制滤波器权系数向量AB(n)收敛的步长参数,取值区间为(0,2);η2为控制滤波器有界约束分解子项的权系数向量C(n)收敛的步长参数,取值区间为(0,2);ρ为平滑因子其取值区间为(0 0.1];
(6):令n=n+1,重复2~5的步骤,即可实时消除接受信号d(n)中的污染噪声v(n),得到实时的有用信号s(n)。
仿真实验
为了验证本发明的基于Laguerre结构的多项式自适应有源噪声对消方法(LSV)的有效性,进行了本发明方法与现有的基于Sigmoid函数的Volterra自适应有源对消方法(SV)的两个仿真实验。
两个仿真实验中,SV滤波器的阶数取6(参数η1=η2=0.8,权系数初始值为小于10-3的随机数),本发明方法滤波器的阶数为3(参数α=0.1,平滑因子ρ=0.1,权系数AB(n)和C(n)均取初始值为小于10-3的随机数)。
实验一
接收信号d(n)中无有用信号,即s(n)=0。
接收信号d(n)中的基本噪声v’(n)由以下方法产生;
v(n)=r(n)+fn(v(n-1))
其中,非线性滤波器fn(·)定义为:
f n ( v ( n - 1 ) ) = 1 2 exp ( - ( v ( n - 1 ) - 1 ) 2 2 &sigma; 2 ) - 1 2 exp ( - ( v ( n - 1 ) + 1 ) 2 2 &sigma; 2 ) ;
对基本噪声v’(n)附加较强的码间窜扰,即将基本噪声v’(n)畸变为v(n),v(n)=v’(n)+0.2v’(n-1),畸变后的噪声v(n)即作为接收信号d(n)中含有(混入)的噪声v(n)。
图1a是本发明方法与SV方法收敛曲线比较图,图中纵坐标为噪声估计均方误差(100次重复实验得到的噪声v(n)的估计y(n)与输入的噪声v(n)的差值的平方的平均值),横坐标表示迭代次数;其中由黑色的点构成的曲线为本发明方法的收敛曲线,由灰色的点构成的曲线为SV方法的收敛曲线。由该图可知,本发明的收敛速度明显优于SV方法,而其降噪比稍好于SV方法,且降噪比为-62dB左右,在极短的时间内就可以收敛且大约在100次迭代时达到较低的稳态误差。
图1b是本发明一次实验所得的有用信号的估计e(n)。由图1b可见,本发明得出的有用信号的估计e(n)与输入的有用信号s(n)=0相比,失真极小。
实验二:
接收信号d(n)中无有用信号,即s(n)=0。
接收信号d(n)中的基本噪声v’(n)产生方法与实验一相同。
对基本噪声v’(n)附加较强的非线性畸变,即将基本噪声v’(n)畸变为v(n),v(n)=v’(n)+0.2v’2(n),畸变后的噪声v(n)即作为接收信号d(n)中含有(混入)的噪声v(n)。
图2a是本发明方法与SV方法收敛曲线比较图,图中纵坐标为噪声估计均方误差(100次重复实验得到的噪声v(n)的估计y(n)与输入的噪声v(n)的差值的平方的平均值),横坐标表示迭代次数;其中由黑色的点构成的曲线为本发明方法的收敛曲线,由灰色的点构成的曲线为SV方法的收敛曲线。由该图可知,由图可知本发明的降噪比、收敛速度性能均明显优于SV方法,并且其降噪比为-57dB左右,相对于SV方法减少了4dB左右,大约在时刻200处就可以收敛到非常小的稳态误差。
图2b是本发明一次实验所得的有用信号估计e(n)。由图2b可见,本发明得出的有用信号的估计e(n)与输入的有用信号s(n)=0相比,失真极小。
上述仿真实验一、二的结果表明,本发明能够有效地消除噪声。特别地,仿真实验二表明,对于噪声中有附加的强的非线性畸变的情况下,它在降噪比、收敛速度等方面表现出更好的性能。而且,由于其信号输入长度为3,而现有的SV输入长度为6,因此,本发明的计算量也得到明显降低。

Claims (1)

1.一种基于Laguerre结构的多项式自适应有源噪声对消方法,其步骤如下:
(1)参考噪声源的产生
将取值在[-1,1]区间且方差σ2=3的白噪声,作为对含有用信号s(n)和噪声v(n)的接收信号d(n)进行噪声对消用的参考噪声源r(n),其中n为时刻;
(2)基于Laguerre延迟节点的输出信号矢量的计算
将参考噪声源r(n)通过Laguerre延迟节点得到m个延迟输出信号分量uk(n),0≤k≤m-1;
当k=0时, u 0 ( n ) = &alpha; u 0 ( n - 1 ) + 1 - &alpha; 2 r ( n - 1 ) ;
当k=1,2,3,...,m-1时,uk(n)=uk-1(n-1)+α[uk(n-1)-uk-1(n)]
其中,k为Laguerre延迟节点的序号,α为Laguerre延迟节点的极点;
由m个延迟输出信号分量uk(n)及常数1构成延迟输出信号向量U(n),
U(n)=[1,u0(n),u1(n),…,um(n)]T
(3)噪声的估计
将第(2)步得到的延迟输出信号向量U(n)输入有界约束乘积耦合Volterra近似结构的滤波器,滤波得到接收信号d(n)中的噪声v(n)的估计y(n),
y ( n ) = a 0 ( n ) + &Sigma; k = 0 m - 1 a k + 1 ( n ) u k ( n ) + [ b 0 ( n ) + &Sigma; k = 0 m - 1 b k + 1 ( n ) u k ( n ) ] &times; Sigm ( c 0 ( n ) + &Sigma; k = 0 m - 1 c k + 1 ( n ) u k ( n ) )
= [ a 0 ( n ) , a 1 ( n ) , &CenterDot; &CenterDot; &CenterDot; , a m ( n ) , b 0 ( n ) , b 1 ( n ) , &CenterDot; &CenterDot; &CenterDot; , b m ( n ) ] [ U ( n ) , U ( n ) Sigm ( C T ( n ) U ( n ) ) ] T
= AB T ( n ) X ( n )
其中,Sigm(·)为有界函数,其表达式为Sigm(·)=(1-exp(·))/(1+exp(·))且|Sigm(·)|<1,X(n)=[U(n),U(n)Sigm(CT(n)U(n))]T是滤波器的输入信号向量,其对应的滤波器权系数为AB(n)=[a0(n),a1(n),…,am(n),b0(n),b1(n),…,bm(n)]T,而C(n)=[c0(n),c1(n),…,cm(n)]T为n时刻滤波器有界约束分解子项的权系数向量;
(4)有用信号的估计
将接收信号d(n)减去第(3)步得到的噪声v(n)的估计y(n),得到n时刻接收信号d(n)中的有用信号s(n)的估计e(n):
e(n)=d(n)-y(n)=s(n)+v(n)-y(n)
(5)权系数更新
将第(4)步得到的有用信号估计e(n)的平方,求期望得到代价函数J(n),J(n)=E{e2(n)}=E{|d(n)-y(n)|2},E{·}为求期望,再用代价函数J(n)对非线性滤波器的权系数AB(n)和C(n)求瞬时梯度,并利用输入信号和瞬时误差归一化的LMS算法准则,得到权系数AB(n)和C(n)的如下迭代公式:
&sigma; X 2 ( n ) = &rho; &sigma; X 2 ( n - 1 ) + ( 1 - &rho; ) [ | | X ( n ) | | 2 + ( 2 m + 2 ) e 2 ( n ) ]
AB ( n + 1 ) = AB ( n ) + &eta; 1 e ( n ) &eta; 1 &sigma; X 2 ( n ) X ( n )
&sigma; U 2 ( n ) = &rho; &sigma; U 2 ( n - 1 ) + ( 1 - &rho; ) ( | | U ( n ) | | 2 + ( m + 1 ) e 2 ( n ) [ b ( n ) T U ( n ) ( 1 - Sigm 2 ( C T ( n ) U ( n ) ) ) ] 2 )
C ( n + 1 ) = C ( n ) + &eta; 2 e ( n ) &eta; 2 &sigma; U 2 ( n ) b ( n ) T U ( n ) ( 1 - Sigm 2 ( C T ( n ) U ( n ) ) )
其中
Figure FDA00001879139000025
Figure FDA00001879139000026
分别为滤波器输入信号向量X(n)和U(n)的功率;b(n)=[b0(n),b1(n),…,bm(n)]T;η1为控制滤波器权系数向量AB(n)收敛的步长参数,取值区间为(0,2);η2为控制滤波器有界约束分解子项的权系数向量C(n)收敛的步长参数,取值区间为(0,2);ρ为平滑因子其取值区间为(0 0.1];
(6)令n=n+1,重复2~5的步骤,即可实时消除接受信号d(n)中的污染噪声v(n),得到实时的有用信号s(n)的估计e(n)。
CN201210240334.7A 2012-07-12 2012-07-12 基于Laguerre结构的多项式自适应有源噪声对消方法 Expired - Fee Related CN102724152B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210240334.7A CN102724152B (zh) 2012-07-12 2012-07-12 基于Laguerre结构的多项式自适应有源噪声对消方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210240334.7A CN102724152B (zh) 2012-07-12 2012-07-12 基于Laguerre结构的多项式自适应有源噪声对消方法

Publications (2)

Publication Number Publication Date
CN102724152A true CN102724152A (zh) 2012-10-10
CN102724152B CN102724152B (zh) 2014-10-08

Family

ID=46949817

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210240334.7A Expired - Fee Related CN102724152B (zh) 2012-07-12 2012-07-12 基于Laguerre结构的多项式自适应有源噪声对消方法

Country Status (1)

Country Link
CN (1) CN102724152B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105005197A (zh) * 2015-05-29 2015-10-28 北京航空航天大学 基于切比雪夫多项式展开的时变神经动力学系统辨识方法
CN106170049A (zh) * 2016-05-12 2016-11-30 西南交通大学 一种有偏补偿的归一化子带自适应回声消除方法
CN109977350A (zh) * 2017-12-27 2019-07-05 核动力运行研究所 一种超声涡流信号自适应nnt对消方法
CN112466271A (zh) * 2020-11-30 2021-03-09 声耕智能科技(西安)研究院有限公司 一种分布式主动噪声控制方法、系统、设备及存储介质
CN114186435A (zh) * 2021-12-20 2022-03-15 大连理工大学 一种基于Volterra级数和极点-留数操作的非线性动力响应预测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1773853A (zh) * 2004-11-08 2006-05-17 西南交通大学 一种基于Laguerre结构的自适应非线性滤波器
CN102187387A (zh) * 2008-10-20 2011-09-14 伯斯有限公司 调节有源降噪滤波器的泄漏因子

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1773853A (zh) * 2004-11-08 2006-05-17 西南交通大学 一种基于Laguerre结构的自适应非线性滤波器
CN102187387A (zh) * 2008-10-20 2011-09-14 伯斯有限公司 调节有源降噪滤波器的泄漏因子

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵海全,张家树: "非线性通信信道的神经FIR自适应幅值Laguerre均衡器", 《中国科学》, no. 10, 31 October 2009 (2009-10-31) *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105005197A (zh) * 2015-05-29 2015-10-28 北京航空航天大学 基于切比雪夫多项式展开的时变神经动力学系统辨识方法
CN106170049A (zh) * 2016-05-12 2016-11-30 西南交通大学 一种有偏补偿的归一化子带自适应回声消除方法
CN106170049B (zh) * 2016-05-12 2019-05-17 西南交通大学 一种有偏补偿的归一化子带自适应回声消除方法
CN109977350A (zh) * 2017-12-27 2019-07-05 核动力运行研究所 一种超声涡流信号自适应nnt对消方法
CN109977350B (zh) * 2017-12-27 2023-01-24 核动力运行研究所 一种超声涡流信号自适应径向基函数神经网络对消方法
CN112466271A (zh) * 2020-11-30 2021-03-09 声耕智能科技(西安)研究院有限公司 一种分布式主动噪声控制方法、系统、设备及存储介质
CN114186435A (zh) * 2021-12-20 2022-03-15 大连理工大学 一种基于Volterra级数和极点-留数操作的非线性动力响应预测方法
CN114186435B (zh) * 2021-12-20 2024-10-08 大连理工大学 一种基于Volterra级数和极点-留数操作的非线性动力响应预测方法

Also Published As

Publication number Publication date
CN102724152B (zh) 2014-10-08

Similar Documents

Publication Publication Date Title
CN109741727B (zh) 基于有源噪声控制算法的主动降噪耳机、降噪方法及存储介质
CN102724152B (zh) 基于Laguerre结构的多项式自适应有源噪声对消方法
CN110492868B (zh) 一种多参数变步长lms自适应滤波方法
CN107342751B (zh) 一种基于互相关熵的变步长sa自适应滤波算法
CN105955025B (zh) 一种广义不确定系统的干扰补偿与抑制方法
Liu et al. Robust Geman-McClure based nonlinear spline adaptive filter against impulsive noise
Zhao et al. A novel adaptive nonlinear filter-based pipelined feedforward second-order Volterra architecture
Yang et al. Efficient μ-law improved proportionate affine projection algorithm for echo cancellation
Talebi et al. Complex-Valued Nonlinear Adaptive Filters With Applications in $\alpha $-Stable Environments
CN113111505A (zh) 基于非线性Hammerstein系统的变遗忘因子递推最小二乘方法及系统
CN112953463B (zh) 受约束的带遗忘因子递归最大相关熵自适应滤波器
Fang et al. On the convergence of Volterra filter equalizers using a pth-order inverse approach
Prongnuch et al. Stability and steady-state performance of hammerstein spline adaptive filter based on stochastic gradient algorithm
Qian et al. Fixed-point maximum total complex correntropy algorithm for adaptive filter
Feng et al. Fast approximate inverse power iteration algorithm for adaptive total least-squares FIR filtering
CN112886947A (zh) 一种变步长鲁棒仿射投影自适应滤波器
Sankaran et al. Fast generalized affine projection algorithm
Waheed et al. Blind source recovery: A framework in the state space
Aref et al. Variable step size modified clipped LMS algorithm
CN114094987A (zh) 基于改进tanh函数的变步长LMP自适应滤波算法
Zang et al. Continuous-time envelope-constrained filter design via Laguerre filters and/spl Hscr//sub/spl infin//optimization methods
Moir A control theoretical approach to the polynomial spectral-factorization problem
CN117040489B (zh) 一种稀疏约束的样条自适应滤波器
CN109639258A (zh) 一种基于Hopfield神经网络二维FIR陷波滤波器的设计方法
Barros et al. An algorithm based on the even moments of the error

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141008

Termination date: 20170712