自适应反馈有源控制系统水床效应抑制的方法

摘要

本发明公开了一种自适应反馈有源控制系统水床效应抑制的方法，方法包括以下步骤：(1)次级路径建模；(2)设计频谱整形滤波器幅频响应，并计算其单位冲激响应；(3)将白噪声经过频谱整形滤波器滤波得到信号a(n)；(4)将信号a(n)经过控制滤波器滤波，得到信号b(n)；(5)使用误差传感器所拾取的误差信号e(n)减去馈给次级声源的控制信号y(n)经次级路径模型滤波后的信号合成参考信号x(n)；(6)根据迭代公式对控制滤波器系数迭代；(7)不断迭代控制滤波器系数，使得代价函数最小，本发明能够在自适应系统中抑制反馈有源控制系统因水床效应引起的噪声放大，并且计算量小。

说明书

技术领域

本发明属于有源噪声控制的技术领域，具体涉及一种自适应反馈有源控制系统水床效应抑制的方法。

背景技术

有源噪声控制(Active Noise Control)已广泛应用于有源降噪耳机、汽车有源降噪、降低机器噪声、有源声屏障和降低室内噪声等场合，在新兴领域也有应用前景，如降低机器人使用环境的噪声，提高其识别语音的能力。

有源噪声控制结构可分为前馈有源噪声控制和反馈有源噪声控制，其中前者性能依赖于参考信号与初级噪声信号之间的相干性。在许多ANC系统的应用场景中，无法获取高质量的参考信号，降噪性能弱或无法降噪，可使用反馈有源控制系统。然而反馈系统存在“水床”效应，即在某特定频段降噪的同时会引起其他频段噪声放大。

针对该问题，人们通常使用离线算法设计控制滤波器，并针对离线算法设计改善水床效应的方法，如Rafaely等提出用H

基于此，有必要针对反馈有源噪声控制系统，提供一种可以通过对特定频段的控制滤波器幅值抑制来抑制水床效应引起的噪声放大，且计算量较小的算法。

发明内容

发明目的：为了克服现有技术中存在的不足，本发明提供一种自适应反馈有源控制系统水床效应抑制的方法，与泄漏滤波-x最小均方算法相比，可在自适应反馈有源控制系统中实现对一个或多个特定频段的水床抑制，与现有可对特定频段水床效应进行抑制的算法即广义泄漏滤波-x最小均方算法相比，运算量大幅降低，更利于实用。

技术方案：为实现上述目的，本发明采用的技术方案为：

一种自适应反馈有源控制系统水床效应抑制的方法，包括如下步骤：

步骤1，测量有源控制系统次级路径传递函数单位冲激响应，得到次级路径传递函数的估计向量

步骤2，根据需要抑制水床效应的频段设计频谱整形滤波器A(ω)的幅频响应，在需要抑制的频段频谱整形滤波器A(ω)的幅值大于其他频段的幅值，并根据需要抑制的频段频谱整形滤波器A(ω)的幅频响应计算其单位冲激响应。

步骤3，将白噪声σ(n)经过频谱整形滤波器A(ω)滤波得到频谱整形滤波信号a(n)。

步骤4，将频谱整形滤波信号a(n)经过控制滤波器W(ω)滤波，得到控制滤波信号b(n)。

步骤5，使用误差传感器所拾取的误差信号e(n)减去馈给次级声源的控制信号y(n)经次级路径建模滤波器

步骤6，通过公式

步骤7，不断迭代控制滤波器系数w(n)，使得代价函数J＝E[e

优选的：步骤6控制滤波器系数迭代公式的得到方法包括以下步骤：

步骤61，代价函数的频域形式定义为：

式中，P

步骤62，将误差信号e(n)＝p(n)+w

J＝E[p

式中，p(n)为初始噪声信号，R为滤波-x信号向量r(n)自相关矩阵，A为频谱整形滤波信号向量a(n)的自相关矩阵，P为初始噪声信号p(n)和滤波-x信号向量r(n)的互相关向量。

步骤63，代价函数J得瞬时梯度

步骤64，使用随机梯度下降法，得控制滤波器系数迭代公式。

优选的：初始噪声信号p(n)和滤波-x信号向量r(n)的互相关向量P：P＝E[p(n)r(n)]。

优选的：滤波-x信号向量r(n)自相关矩阵R：R＝E[r(n)r

优选的：频谱整形滤波信号向量a(n)的自相关矩阵A：A＝E[a(n)a

本发明相比现有技术，具有以下有益效果：

(1)本发明提出的方法与非自适应滤波器设计方法相比，可用于实时自适应系统。

(2)本发明与LFxLMS算法相比，可在自适应反馈有源控制系统中实现对一个或多个特定频段的水床抑制，降噪效果更好。

(3)本发明与GLFxLMS算法相比，避免了迭代过程中矩阵与向量相乘的复杂操作，运算量大幅降低，更具有实用性。

附图说明

图1是本发明的算法框图。

图2是不同阶数的频谱整形滤波器幅频响应。

图3是本发明方法与传统FxLMS算法，LFxLMS算法和GLFxLMS算法降噪前后噪声信号功率谱密度对比图。

图4是本发明算法与传统FxLMS算法，LFxLMS算法和GLFxLMS算法稳态滤波器幅频响应对比图。

具体实施方式

下面结合附图和具体实施例，进一步阐明本发明，应理解这些实例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

一种自适应反馈有源控制系统水床效应抑制的方法，如图1所示，通过以下技术方案来实现：

(1)测量有源控制系统次级路径传递函数单位冲激响应，得到次级路径函数的估计

(2)根据需要抑制水床效应的频段设计频谱整形滤波器A(ω)的幅频响应，在需要抑制的频段A(ω)幅值设计为较大，其他频段A(ω)幅值设计为相对较小，根据幅频响应计算其单位冲击响应。

(3)将白噪声σ(n)经过频谱整形滤波器A(ω)滤波得到信号a(n)。

(4)将信号a(n)经过控制滤波器W(ω)滤波，得到信号b(n)。

(5)使用误差传感器所拾取的误差信号e(n)减去馈给次级声源的控制信号y(n)经次级路径建模滤波器

(6)通过公式

(7)不断迭代控制滤波器系数w(n)，使得代价函数J＝E[e

根据图1所示的本发明算法框图，其原理简介如下：

设定控制滤波器长度L，频谱整形滤波器长度L

本发明方法的代价函数的时域形式为：

J＝E[e

式中b(n)＝w

根据维纳-辛钦定理，代价函数的频域形式定义为：

式中P

将误差信号e(n)＝p(n)+w

J＝E[p

式中R和A分别为滤波-x信号向量r(n)和a(n)的自相关矩阵，P为初始噪声信号p(n)和r(n)的互相关向量，分别定义如下：

P＝E[p(n)r(n)] (4)

R＝E[r(n)r

A＝E[a(n)a

式中r(n)＝[r(n),r(n-1),…,r(n-L+1)]

使用随机梯度下降法，得控制滤波器系数迭代公式为：

式中μ为迭代步长。因该算法可实现对特定频段控制滤波器约束，故称其为频段约束滤波-x最小均方(Frequency band constrained Filtered-x Least Mean Square)算法，简称FcFxLMS算法。

下面以有源降噪头靠系统为例说明本发明的效果，有源降噪头靠系统一般布置在人耳附近，使用扬声器作为次级声源，通过降低人耳附近误差传声器的噪声，在人耳附近产生静区。每只耳朵使用一个控制通道，包含一个次级声源和一个误差传声器，误差传声器距离次级声源振膜中心8cm。2只耳朵共2个通道，一般而言，有源降噪头靠系统的2个通道耦合较弱，每一通道可视为单通道系统。由于头靠系统使用场合限制，扬声器尺寸不能太大，如集成于汽车、飞机座位的扬声器，其低频响应有限，例如本实施例中作为次级声源的扬声器在小于100Hz的低频无法产生足够大的声压，若次级源在小于100Hz的低频输出较大，有可能导致次级源输出信号失真，额外产生噪声，甚至损坏扬声器，因此有必要对该频段的滤波器幅值进行约束。另一方面，人耳对某些中高频段的信号相对更为敏感，不希望该频段的噪声被放大。因此故将频谱整形滤波器A(ω)设置为一个低通滤波器与一个带通滤波器组合的形式。其中低通滤波器的通带截止频率为150Hz，阻带截止频率为200Hz，通带幅值和阻带幅值分别为0.01和0.0018。带通滤波器的通带上、下截止频率分别为1000Hz和2000Hz，过渡带宽为100Hz，通带幅值和阻带幅值分别为0.01和0.0018。根据幅频响应计算其单位冲激响应。设计了频谱整形滤波器A(ω)，可使用上述本发明技术方案(1)～(7)来实现本发明，其中方案中的第(2)项即将频谱整形滤波器A(ω设置为此处描述的低通滤波器与带通滤波器组合。

为了说明本发明方法的优点，以下利用仿真对传统的LFxLMS算法、GLFxLMS算法和本发明提出的FcFxLMS算法的稳态性能进行比较验证。使用上述有源降噪头靠系统实测的传递函数进行仿真，使用一个音箱作为初级声源，初级噪声为方差为10

表1不同算法降噪性能和计算量对比

图3和表1显示：使用FxLMS算法可获得最大降噪量(18.6dB)，但其在1000-2000Hz噪声放大量较大，同时控制滤波器在100Hz以下幅值较大，其峰值为-1.9dB。使用LFxLMS可有效降低由水床效应引起的噪声放大，同时较好地抑制100Hz以下控制滤波器幅值，但会牺牲降噪量(仅有10.6dB)。使用GLFxLMS和本发明提出的FcFxLMS获得效果类似，降噪量为13.1dB，控制滤波器100Hz以下幅度峰值为-16.0dB，相对于LFxLMS，以牺牲部分降噪量降低了1000-2000Hz的噪声放大。但FcFxLMS的计算量远小于GLFxLMS，当滤波器阶数为512频谱整型滤波器的数L

使用以上几种算法，得到的收敛后的控制滤波器w(n)的幅频响应如图4所示。本发明所提的FcFxLMS算法和传统的GLFxLMS算法可实现对控制滤波器不同频段的加权约束，权重项可自由设计。本例中，在100Hz以下和1000-2000Hz频段设置了较大权重以抑制该频段内控制滤波器的幅值，在其他频段使用相对较小权重以减小降噪量的下降。FcFxLMS算法和GLFxLMS算法收敛后的控制滤波器基本完全重合，如图4中实线和方形标记线所示。

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。