CN104575512B - 基于Volterra滤波器的非线性窄带主动噪声控制方法 - Google Patents
基于Volterra滤波器的非线性窄带主动噪声控制方法 Download PDFInfo
- Publication number
- CN104575512B CN104575512B CN201410834232.7A CN201410834232A CN104575512B CN 104575512 B CN104575512 B CN 104575512B CN 201410834232 A CN201410834232 A CN 201410834232A CN 104575512 B CN104575512 B CN 104575512B
- Authority
- CN
- China
- Prior art keywords
- linear
- noise
- module
- signal
- volterra
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000009021 linear effect Effects 0.000 claims abstract description 98
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 34
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 34
- 230000009022 nonlinear effect Effects 0.000 claims abstract description 15
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims abstract description 7
- 230000000644 propagated effect Effects 0.000 claims description 5
- 230000002401 inhibitory effect Effects 0.000 claims description 3
- 230000000717 retained effect Effects 0.000 claims description 2
- 238000011478 gradient descent method Methods 0.000 claims 1
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 3
- 208000036829 Device dislocation Diseases 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Noise Elimination (AREA)
Abstract
本发明提出了一种基于Volterra滤波器的非线性窄带主动噪声控制方法，针对窄带噪声源受到非线性扭曲产生高次谐波噪声的情形；其结构包括线性信号合成子系统（1）、一个非线性信号合成子系统（2）以及参考信号生成模块以及最小均方算法模块；该方法的参考信号由同步信号发生器根据非声学传感器获得的源噪声频率产生；通过一个线性组合器合成一部分次级源，通过另一个线性组合器和非线性Volterra滤波器合成另一部分次级源，在初级源噪声受到非线性扭曲时，线性和非线性滤波器分别独立抑制基频分量和高次谐波分量；该方法能有效抑制受非线性扭曲的窄带噪声，且独立工作的结构易于分析非线性对系统性能的影响。
Description
技术领域
本发明涉及一种基于Volterra滤波器的非线性窄带主动噪声控制方法，属于主动噪声控制领域。
背景技术
主动噪声控制相对传统被动噪声控制，具有低频性好，控制空间小，已经广泛的应用到噪声控制领域。其中，窄带主动噪声控制应用较多，窄带噪声具有规律性的低频信号特征，往往是由旋转设备或具有往复运动的装置产生。在窄带主动噪声控制中，窄带噪声频率可以通过非声学传感器直接或间接测得，再经信号发生器可产生同频率的正弦或余弦参考信号。
在线性或者近似线性声学环境，线性窄带主动噪声控制可以有效应对目标噪声。然而，当实际声学环境具有较强非线性时，源噪声经非线性通道传播在相消点将产生严重非线性扭曲或扰动。此时，线性窄带主动噪声控制系统将难以应对。
考虑到窄带噪声经由非线性通道传播后产生的非线性扭曲实际上激发出了高次谐波噪声，通过非线性初级通道的窄带噪声在相消点形成的目标噪声可分为原噪声源(即基波分量噪声)和受到非线性影响产生的高次谐波分量噪声。因此，消除两种分量噪声可分别使用线性滤波器和非线性滤波器。基于此，可以构造一种线性信号合成子系统和非线性信号合成子系统分别独立工作的非线性窄带主动噪声控制系统。
发明内容
所要解决的技术问题：
本发明的目的是针对上述问题，设计一种有效降低经过非线性初级通道扭曲的窄带噪声的非线性窄带主动噪声控制方法。
技术方案：
为了实现以上功能，本发明提供了一种基于Volterra滤波器的非线性窄带主动噪声控制方法，其特征在于：该方法提供一个线性信号合成子系统1、一个非线性信号合成子系统2以及参考信号生成模块；
所述线性信号合成子系统1包括第一线性组合器模块21和次级通道，用于抑制基频噪声成份；所述非线性信号合成子系统2包括第二线性组合器模块22和非线性Volterra滤波器模块3，非线性Volterra滤波器模块3仅保留2阶以上核及其系数用于抑制高次谐波噪声成份；
该控制方法依照以下步骤：
步骤1，参考信号生成模块，根据由非声学传感器获得的同步信号产生正弦型参考信号，提供给线性信号合成子系统1和非线性信号合成子系统2；
步骤2，在所述线性信号合成子系统1中，接收到的参考信号由第一线性组合器模块21进行信号合成，输出信号y0(n)；
在所述非线性信号合成子系统2中，接收到的参考信号由第二线性组合器模块22进行信号合成，其输出信号为非线性Volterra滤波器模块3的参考信号；
步骤3，所述非线性Volterra滤波器模块3对其参考信号进行合成，输出信号y1(n)；该输出信号y1(n)叠加至第一线性组合器模块21的输出信号y0(n)上，并与目标噪声p(n)进行相消处理，得到相消后的误差e0(n)；
步骤4，对相消后的误差e0(n)经过次级通道s(z)，并引入不可消除的加性背景噪声，得到残余噪声e(n)。
该方法还提供一个最小均方算法模块4，将次级通道s(z)采用离线辨识方法进行估计，并得到估计模型对所述正弦型参考信号经估计模型处理，并且对所述非线性Volterra滤波器模块3的参考信号经估计模型处理，得到滤波-X信号，连同残余噪声e(n)一起共同作为最小均方算法模块4的输入信号；
所述第一线性组合器模块21、第二线性组合器模块22为线性滤波器，最小均方算法模块4是第一线性组合器模块21、第二线性组合器模块22中的滤波器权值的更新控制器，并且最小均方算法模块4为非线性Volterra滤波器模块3提供系数更新控制；使用梯度下降法对权值和系数进行更新，通过调节更新滤波器的权值使合成的次级噪声与目标噪声相消后所得的残余噪声能量最小。
所述目标噪声p(n)由窄带噪声源经过非线性初级通道传播产生非线性扭曲后形成的。
有益效果：
本发明提出的非线性窄带主动噪声控制方法中，线性信号合成子系统和非线性信号合成子系统分别独立工作的结构，能有效降低经过非线性初级通道扭曲的窄带噪声。
本方法中，所使用的非线性滤波器，即Volterra滤波器，能够以较低的阶数应对较强的初级通道非线性；在传统线性窄带主动噪声控制方法基础上增添了非线性Volterra滤波器，该滤波器在实际的主动噪声控制器中易于实现；本方法通过简单的并行扩展，易于实现多频率通道的主动噪声控制。
附图说明
下面结合附图和实施例对本发明进一步说明：
图1本发明提供的基于Volterra滤波器的非线性窄带主动噪声控制方法的流程图；
其中，1-线性信号合成子系统，2-非线性信号合成子系统，3-非线性Volterra滤波器模块，4-最小均方算法模块，21-第一线性组合器模块，22-第二线性组合器模块。
具体实施方式
本发明提供一种基于Volterra滤波器的非线性窄带主动噪声控制方法，为使本发明的目的，技术方案及效果更加清楚，明确，以及参照附图并举实例对本发明进一步详细说明。应当理解，此处所描述的具体实施仅用以解释本发明，并不用于限定本发明。
本发明提供的基于Volterra滤波器的非线性窄带主动噪声控制方法具体流程如图1所示，该方法提供一个线性信号合成子系统1、一个非线性信号合成子系统2以及参考信号生成模块；
所述线性信号合成子系统1包括第一线性组合器模块21和次级通道，用于抑制基频噪声成份；所述非线性信号合成子系统2包括第二线性组合器模块(22)和非线性Volterra滤波器模块3，非线性Volterra滤波器模块3仅保留2阶以上核及其系数用于抑制高次谐波噪声成份；
该控制方法依照以下步骤：
步骤1，参考信号生成模块，根据由非声学传感器获得的同步信号产生正弦型参考信号，提供给线性信号合成子系统1和非线性信号合成子系统2。
目标噪声p(n)是初级窄带源噪声经由非线性初级通道传播，产生非线性扭曲后在相消点形成的噪声。实际上，对于窄带噪声，这种非线性扭曲激发出了高次谐波分量噪声。通常，p(n)可表示为参考信号x(n)的多项式，即
p(n)＝a1x(n)+a2x2(n)+a3x3(n)
其中x(n)＝axa(n)+bxb(n)，xa(n)＝cos(ωn)，xb(n)＝sin(ωn)
图中参考信号生成模块(1)中的cos(ωn)和sin(ωn)由函数信号发生器依据噪声源的频率ω产生，频率ω通过非声学传感器测得。
步骤2，在所述线性信号合成子系统1中，接收到的参考信号由第一线性组合器模块21进行信号合成，输出信号y0(n)；
在所述非线性信号合成子系统2中，接收到的参考信号由第二线性组合器模块22进行信号合成，其输出信号为非线性Volterra滤波器模块3的参考信号。
第一线性组合器模块21作为一个线性滤波器用于控制受非线性扭曲的目标噪声的基波分量。输入参考信号xa(n)、xb(n)通过第一线性组合器模块21产生的次级噪声源信号为
第二线性组合器模块22产生的输出信号作为非线性Volterra滤波器的参考输入信号。Volterra滤波器保留2阶以上核及其系数即能能有效抑制高次谐波分量噪声。另一方面，考虑到对称特性，Volterra滤波器权值删减一半并不影响滤波器的性能。因此，Volterra滤波器输出的另一次级噪声源信号为
该方法还提供一个最小均方算法模块4，将次级通道s(z)采用离线辨识方法进行估计，得到估计模型对所述正弦参考信号经估计模型处理，并且对所述非线性Volterra滤波器模块3的参考信号经估计模型处理，得到滤波-X信号，连同残余噪声e(n)一起共同作为最小均方算法模块4的输入信号；
所述第一线性组合器模块21、第二线性组合器模块22为线性滤波器，最小均方算法模块4是第一线性组合器模块21、第二线性组合器模块22中的滤波器权值的更新控制器，并且最小均方算法模块4为非线性Volterra滤波器模块3提供系数更新控制；使用梯度下降法对权值和系数进行更新，通过调节更新滤波器的权值使合成的次级噪声与目标噪声相消后所得的残余噪声能量最小。
通过滤波-X最小均方(FXLMS)算法更新滤波器权值，包括使产生的次级噪声源能有效抑制目标噪声p(n)。
图中的第一线性组合器模块21和非线性Volterra滤波器模块3的权值按如下方程更新(Volterra滤波器权值更新方程仅给出前3阶)：
其中，滤波-X信号是参考信号xa(n)、xb(n)、xr(n)分别经过次级通道估计所得。真实次级通道s(z)及其估计通常是由FIR滤波器构成，即
经过次级通道估计滤波后的滤波-X信号为
步骤3，所述非线性Volterra滤波器模块3对其参考信号进行合成，输出信号y1(n)；该输出信号y1(n)叠加至第一线性组合器模块21的输出信号y0(n)上，并与实际噪声p(n)进行相消处理，得到相消后的噪声e0(n)。
使J(n)分别对求偏导
根据FXLMS算法，系数的更新方程为
其中
步骤4，对相消后的误差e0(n)经过次级通道s(z)，并引入不可消除的加性背景噪声，得到残余噪声e(n)。
附图中目标噪声抵消后的残余噪声信号e(n)可计算为：
其中
e0(n)＝p(n)-y0(n)-y1(n)；
vp(n)是加性背景噪声，通常是由均值为零的高斯白噪声组成，其方差为
可以理解的是，对本领域普通技术人员来说，可以根据本发明的技术方案及其发明构思加以等同替换或改变，而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。
Claims (3)
1.基于Volterra滤波器的非线性窄带主动噪声控制方法，其特征在于：该方法提供一个线性信号合成子系统(1)、一个非线性信号合成子系统(2)以及参考信号生成模块；
所述线性信号合成子系统(1)包括第一线性组合器模块(21)和次级通道，用于抑制基频噪声成份；所述非线性信号合成子系统(2)包括第二线性组合器模块(22)和非线性Volterra滤波器模块(3)，非线性Volterra滤波器模块(3)仅保留2阶以上核及其系数用于抑制高次谐波噪声成份；
该控制方法依照以下步骤：
步骤1，参考信号生成模块，根据由非声学传感器获得的同步信号产生正弦型参考信号，提供给线性信号合成子系统(1)和非线性信号合成子系统(2)；
步骤2，在所述线性信号合成子系统(1)中，接收到的参考信号由第一线性组合器模块(21)进行信号合成，输出信号y0(n)，为滤波器权值，xa(n)、xb(n)为输入参考信号；
在所述非线性信号合成子系统(2)中，接收到的参考信号由第二线性组合器模块(22)进行信号合成，其输出信号为非线性Volterra滤波器模块(3)的参考信号；所述第二线性组合器模块(22)产生的输出信号为
步骤3，所述非线性Volterra滤波器模块(3)对其参考信号进行合成，输出信号y1(n)；该输出信号y1(n)叠加至第一线性组合器模块(21)的输出信号y0(n)上，并与目标噪声p(n)进行相消处理，得到相消后的误差e0(n)；
步骤4，对相消后的误差e0(n)经过次级通道s(z)，并引入不可消除的加性背景噪声，得到残余噪声e(n)。
2.根据权利要求1所述的基于Volterra滤波器的非线性窄带主动噪声控制方法，其特征在于：该方法还提供一个最小均方算法模块(4)，将次级通道s(z)采用离线辨识方法进行估计，并得到估计模型对所述正弦型参考信号经估计模型处理，并且对所述非线性Volterra滤波器模块(3)的参考信号经估计模型处理，得到滤波-X信号，连同残余噪声e(n)一起共同作为最小均方算法模块(4)的输入信号；
所述第一线性组合器模块(21)、第二线性组合器模块(22)为线性滤波器，最小均方算法模块(4)是第一线性组合器模块(21)、第二线性组合器模块(22)中的滤波器权值的更新控制器，并且最小均方算法模块(4)为非线性Volterra滤波器模块(3)提供系数更新控制；使用梯度下降法对权值和系数进行更新，通过调节更新滤波器的权值使合成的次级噪声与目标噪声相消后所得的残余噪声能量最小。
3.根据权利要求1或2所述的基于Volterra滤波器的非线性窄带主动噪声控制方法，其特征在于：所述目标噪声p(n)由窄带噪声源经过非线性初级通道传播产生非线性扭曲后形成的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410834232.7A CN104575512B (zh) | 2014-12-29 | 2014-12-29 | 基于Volterra滤波器的非线性窄带主动噪声控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410834232.7A CN104575512B (zh) | 2014-12-29 | 2014-12-29 | 基于Volterra滤波器的非线性窄带主动噪声控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104575512A CN104575512A (zh) | 2015-04-29 |
| CN104575512B true CN104575512B (zh) | 2018-06-26 |
Family
ID=53091412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410834232.7A Expired - Fee Related CN104575512B (zh) | 2014-12-29 | 2014-12-29 | 基于Volterra滤波器的非线性窄带主动噪声控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104575512B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105575382B (zh) * | 2015-12-29 | 2019-09-13 | 哈尔滨工业大学 | 一种快速稳定收敛的完全并行窄带主动噪声控制方法 |
| CN107315343B (zh) * | 2017-02-21 | 2020-02-07 | 中国人民解放军海军工程大学 | 一种机械主动隔振的多振源多参考窄带自适应方法 |
| CN107102658B (zh) * | 2017-02-24 | 2022-10-04 | 苏州东菱智能减振降噪技术有限公司 | 一种离线次级通道的辨识方法及装置 |
| CN108417220B (zh) * | 2018-02-11 | 2019-06-25 | 陕西师范大学 | 基于代理模型Volterra建模的语音信号编码和解码方法 |
| CN110794681B (zh) * | 2019-11-22 | 2022-05-10 | 西北工业大学 | 一种针对窄带主动振动系统的FxLMS滤波方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101819766A (zh) * | 2010-01-15 | 2010-09-01 | 浙江万里学院 | 一种用于消减噪声的多通道有源噪声控制方法 |
| CN101976560A (zh) * | 2010-09-29 | 2011-02-16 | 哈尔滨工业大学 | 前馈型窄带主动噪声控制系统性能提高的方法 |
| CN102280102A (zh) * | 2010-06-14 | 2011-12-14 | 哈曼贝克自动系统股份有限公司 | 自适应噪声控制 |
| CN103872988A (zh) * | 2014-03-14 | 2014-06-18 | 中国人民解放军海军工程大学 | 一种频率失调下的前馈式窄带主动控制系统 |
-
2014
- 2014-12-29 CN CN201410834232.7A patent/CN104575512B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101819766A (zh) * | 2010-01-15 | 2010-09-01 | 浙江万里学院 | 一种用于消减噪声的多通道有源噪声控制方法 |
| CN102280102A (zh) * | 2010-06-14 | 2011-12-14 | 哈曼贝克自动系统股份有限公司 | 自适应噪声控制 |
| CN101976560A (zh) * | 2010-09-29 | 2011-02-16 | 哈尔滨工业大学 | 前馈型窄带主动噪声控制系统性能提高的方法 |
| CN103872988A (zh) * | 2014-03-14 | 2014-06-18 | 中国人民解放军海军工程大学 | 一种频率失调下的前馈式窄带主动控制系统 |
Non-Patent Citations (1)
| Title |
|---|
| NONLINEAR NARROWBAND ACTIVE NOISE CONTROL USING VOLTERRA FILTER;Jian Liu等;《The 21st International Congress on Sound and Vibration》;20140717;第3-4页第2.2节、图2 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104575512A (zh) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104575512B (zh) | 基于Volterra滤波器的非线性窄带主动噪声控制方法 | |
| CN108428445B (zh) | 一种无误差传声器的自适应主动降噪方法 | |
| CN102843633B (zh) | 扬声器输出的控制 | |
| CN105916079B (zh) | 一种扬声器非线性补偿方法及装置 | |
| US9704470B2 (en) | Method and apparatus for nonlinear compensation in an active noise control system | |
| CN104980877A (zh) | 识别和补偿机电换能器中的非线性振动的装置和方法 | |
| Wu et al. | An M-estimator based algorithm for active impulse-like noise control | |
| Zhu et al. | A novel FELMS-based narrowband active noise control system and its convergence analysis | |
| Edamoto et al. | Virtual sensing technique for feedforward active noise control | |
| Ahmed et al. | A novel approach for improved noise reduction performance in feed-forward active noise control systems with (loudspeaker) saturation non-linearity in the secondary path | |
| Boz et al. | IIR filtering based adaptive active vibration control methodology with online secondary path modeling using PZT actuators | |
| Mazur et al. | Hammerstein nonlinear active noise control with the Filtered-Error LMS algorithm | |
| JP2008250131A (ja) | 能動型騒音制御装置 | |
| CN107315343B (zh) | 一种机械主动隔振的多振源多参考窄带自适应方法 | |
| CN106782594A (zh) | 一种抑制窄带主动噪声控制中频率不匹配的方法 | |
| Abdelrhman et al. | A generalized robust logarithmic family-based adaptive filtering algorithms | |
| Muradore et al. | On the rejection of vibrations in Adaptive Optics Systems | |
| Kukde et al. | Development of a novel narrowband active noise controller in presence of sensor error | |
| Noll et al. | Adaptive equalization of vibrotactile actuators | |
| Sun et al. | A narrowband active noise control system with frequency mismatch compensation | |
| Nakao et al. | An estimation method of parameters for closed-box loudspeaker system | |
| CN103872988A (zh) | 一种频率失调下的前馈式窄带主动控制系统 | |
| Bambang et al. | Recent progress in adaptive nonlinear active noise control | |
| Mohapatra et al. | A sigmoid function based feedback filtered-X-LMS algorithm with improved offline modelling | |
| VC | Novel algorithms based on Legendre neural network for nonlinear active noise control with nonlinear secondary path |
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 | ||
| 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: 20180626 Termination date: 20201229 |