CN111161700A - 一种空间降噪系统中主通道时延估计方法 - Google Patents
一种空间降噪系统中主通道时延估计方法 Download PDFInfo
- Publication number
- CN111161700A CN111161700A CN201911347173.XA CN201911347173A CN111161700A CN 111161700 A CN111161700 A CN 111161700A CN 201911347173 A CN201911347173 A CN 201911347173A CN 111161700 A CN111161700 A CN 111161700A
- Authority
- CN
- China
- Prior art keywords
- delay
- signal
- time delay
- time
- noise
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1781—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions
- G10K11/17821—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions characterised by the analysis of the input signals only
- G10K11/17825—Error signals
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1785—Methods, e.g. algorithms; Devices
- G10K11/17853—Methods, e.g. algorithms; Devices of the filter
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了一种空间降噪系统中主通道时延估计方法,包括以FxLMS算法的频域通道拟合模块以及采用相关法时延估计的时域通道拟合模块,相关法时延估计部分为将参考传声器处的原始噪声与误差传感器处的噪声做时延估计,并将得到的主通道时延返回给系统输出链路,保证输出波形的同步。本发明沿用传统FxLMS算法作为主框架,并在其基础上加入相关法时延估计,通过估计主通道时延并反馈给系统,使得系统具备对主通道的追踪能力。
Description
技术领域
本发明涉及空间降噪系统技术领域,特别涉及到一种空间降噪系统中主通道时延估计方法。
背景技术
噪声控制技术近年来受到广泛研究,对于一定体积且信道复杂的噪音设备,可在一定空间内采用主动噪声控制。采用若干次级声源、参考传声器以及误差传感器包裹目标区域,利用次级声源输出与原有噪声频率相同、相位相反的声波与原有的噪声抵消,达到空间内噪声降低的目的。以上便是现有的空间降噪技术的基本原理。
目前空间降噪技术已在空调柜机、汽车发动机、电力变压器等中大型噪音设备上有了一定的研究与应用,在实际测试中可以很大程度上拟合1000Hz以下的低频噪声,取得了一定的降噪效果。
现有的空间降噪技术需要多个次级声源、参考传声器以及误差传感器,各个单元间的协调问题是目前研究的重点。参考传声器与误差传感器间所构成的物理通道成为主通道,系统传递函数对主通道的拟合程度决定了系统输出与原始噪声抵消后的剩余能量的大小。
对于系统传递函数的拟合,目前的研究一般采用FxLMS(滤波x最小均方误差)算法以获取主通道的频率响应,但是均未考虑实际应用场景下在主通道与次级通道的时延差对系统的影响。
在理想状态下,主通道与次级通道的传输时延差为零,使得次级声源处的系统输出与原始噪声在时域上保持一致,降噪量将完全取决于系统对主通道的估计精度。实际场景下时延差往往达到若干扫描周期,使得抵消波形与原始噪声无法保持时域上的同步。同时,降噪目标区域或伴随一定的机械振动,或在使用过程中受到人为影响,参考传声器与次级声源间的相对位置可能发生变化,致使该时延差无法视作常量或先验获取。
发明内容
本发明的目的在于针对现有技术中的不足,提供一种空间降噪系统中主通道时延估计方法,以解决上述问题。
本发明所解决的技术问题可以采用以下技术方案来实现:
一种空间降噪系统中主通道时延估计方法,包括以FxLMS算法的频域通道拟合模块以及采用相关法时延估计的时域通道拟合模块,相关法时延估计部分为将参考传声器处的原始噪声与误差传感器处的噪声做时延估计,并将得到的主通道时延返回给系统输出链路,保证输出波形的同步。
进一步的,系统输入参考传声器采集到的噪声信号x(n),经过主通道传递函数P(z)后形成误差传感器处输入的原始噪声信号d(n-D),其中D表示该信号相较于x(n)的时延所对应的采样点数;
x(n)输入序列延时n个采样间隔形成Z-n,而后t时刻输入的两路序列做相关运算Rxd(t);
Tc表示系统估计时延的最大范围,根据实际场景决定,对于同一时刻不同估计时延所对应的相关矩阵Rxd(t),选取最大值反馈给次级通路的延时器Z-n,W(n)为噪声抵消部分的自适应滤波器,采用最小均方误差准则,通过误差传感器接收到的系统误差e(n)以及x(n)经过次级通道估计得到的滤波器输入信号xf(n)进行参数调节,滤波器输出u(n)经过次级通道S(z)得到次级输出信号y(n-Δ),其中Δ表示次级通道延时所对应的采样点数,再经延时器D(z)得到系统输出 为时延估计部分对延时D的估计。
进一步的,加入时延估计部分以及延时器单元拓展了系统对主通道时延的追踪性,并须具备以下条件:
D>Δ
即主通道时延D与主次通道时延差d均为正值;
对于时延估计部分,在当前时刻t分别计算时延从1到Tc之间的相关矩阵Rxd(t):
其中,Rpp表示信号p的自相关函数,Rpq表示信号p与信号q的互相关函数,假设信号x(n-t)与信号d(n-D)中的噪声互不相关,则
Rxd(t)=Rxx(t-D)
由自相关函数性质:|Rss(τ-D)|≤Rss(0)
可知,当τ-D=0时,Rss(·)取得最大值,此时τ为当前时刻系统时延的估计
式中,arg为取自变量运算;
相关运算选取的序列长度L越长,估计结果越精确,但运算量会相应增大,考虑系统时延估计的刷新速率有着一定的要求,系统设计时应根据实际情况选取;
Δ为次级通道延迟所对应的采样点数;
对于噪声抵消部分,滤波器W(n)的输出u(n):
u(n)=ωT(n)x(n)
其中,ω(n)=[ω0(n) ω1(n) ω2(n) … ωN-1(n)]T为在时刻n的长度为N的滤波器权系数向量,x(n)=[x(n) x(n-1) x(n-2) … x(n-N+1)]T为时刻n的长度为N的输入信号向量;
u(n)经过经过次级通道S(z)得到系统输出y(n):
y(n-Δ)=s(n)*u(n)
其中,s(n)为S(z)的脉冲响应,*表示卷积运算,y(n-Δ)为由次级声源输出为抵消波形;
计算滤波x信号xf(n):
更新滤波器W(n)的权系数ω(n):
ω(n+1)=ω(n)-μe(n)xf(n)
其中,μ为步长因子,其收敛范围如下:
其中,N为滤波器长度,Px为滤波器输入信号功率;
与现有技术相比,本发明的有益效果如下:
本发明沿用传统FxLMS算法作为主框架,并在其基础上加入相关法时延估计,通过估计主通道时延并反馈给系统,使得系统具备对主通道的追踪能力,从而抑制空间降噪系统中主通道时延无法确定或时变所导致的系统抵消性能下降问题。
相关法时延估计旨在将空间降噪系统中的主通道时延由不可测变量变为可测变量,且该方法具有计算量较小、算法结构相对简单等优点,与FxLMS算法框架形成了很好的衔接。
加入时延估计使系统拥有对主通道时延的估计能力以及对主通道时延变化的追踪能力,实现了原始噪声与抵消波形在时域上更高的同步精度,从而提高噪声抵消效果。
附图说明
图1为本发明所述的空间降噪系统中主通道时延估计方法的系统框图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
参见图1,本发明所述的一种空间降噪系统中主通道时延估计方法,包括以FxLMS算法的频域通道拟合模块以及采用相关法时延估计的时域通道拟合模块。其中FxLMS算法部分沿用传统噪声抵消系统中的算法,原理过程与以往的应用一致,相关法时延估计部分为将参考传声器处的原始噪声与误差传感器处的噪声做时延估计,并将得到的主通道时延返回给系统输出链路,保证输出波形的同步。
图1中x(n)表示系统输入,即参考传声器采集到的噪声信号,d(n-D)表示经过主通道传递函数P(z)后的噪声信号,即误差传感器处输入的原始噪声信号。D表示该信号相较于x(n)的时延所对应的采样点数,Z-n表示输入序列延时n个采样间隔,Rxd(t)表示t时刻输入的两路序列做相关运算,Tc表示系统估计时延的最大范围,根据实际场景决定,对于同一时刻不同估计时延所对应的相关矩阵Rxd(t),选取最大值反馈给次级通路的延时器Z-n,W(n)为噪声抵消部分的自适应滤波器,采用最小均方误差准则,通过误差传感器接收到的系统误差e(n)以及x(n)经过次级通道估计得到的滤波器输入信号xf(n)进行参数调节,滤波器输出u(n)经过次级通道S(z)得到次级输出信号y(n-Δ),其中Δ表示次级通道延时所对应的采样点数,再经延时器D(z)得到系统输出 为时延估计部分对延时D的估计。
加入时延估计部分以及延时器单元拓展了系统对主通道时延的追踪性,并须具备以下条件:
D>Δ
即主通道时延D与主次通道时延差d均为正值。
对于时延估计部分,在当前时刻t分别计算时延从1到Tc之间的相关矩阵Rxd(t):
其中,Rpp表示信号p的自相关函数,Rpq表示信号p与信号q的互相关函数,假设信号x(n-t)与信号d(n-D)中的噪声互不相关,则
Rxd(t)=Rxx(t-D)
由自相关函数性质:|Rss(τ-D)|≤Rss(0)
可知,当τ-D=0时,Rss(·)取得最大值,此时τ为当前时刻系统时延的估计
式中,arg为取自变量运算。
相关运算选取的序列长度L越长,估计结果越精确,但运算量会相应增大,考虑系统时延估计的刷新速率有着一定的要求,系统设计时应根据实际情况选取。
Δ为次级通道延迟所对应的采样点数,一般由ADC、DAC以及系统处理时间所带来的延时确定。
对于噪声抵消部分,滤波器W(n)的输出u(n):
u(n)=ωT(n)x(n)
其中,ω(n)=[ω0(n) ω1(n) ω2(n) … ωN-1(n)]T为在时刻n的长度为N的滤波器权系数向量,x(n)=[x(n) x(n-1) x(n-2) … x(n-N+1)]T为时刻n的长度为N的输入信号向量。
u(n)经过经过次级通道S(z)得到系统输出y(n):
y(n-Δ)=s(n)*u(n)
其中,s(n)为S(z)的脉冲响应,*表示卷积运算,y(n-Δ)为由次级声源输出为抵消波形。
计算滤波x信号xf(n):
更新滤波器W(n)的权系数ω(n):
ω(n+1)=ω(n)-μe(n)xf(n)
其中,μ为步长因子,其收敛范围如下:
其中,N为滤波器长度,Px为滤波器输入信号功率。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (3)
1.一种空间降噪系统中主通道时延估计方法,其特征在于:包括以FxLMS算法的频域通道拟合模块以及采用相关法时延估计的时域通道拟合模块,相关法时延估计部分为将参考传声器处的原始噪声与误差传感器处的噪声做时延估计,并将得到的主通道时延返回给系统输出链路,保证输出波形的同步。
2.根据权利要求1所述的空间降噪系统中主通道时延估计方法,其特征在于:系统输入参考传声器采集到的噪声信号x(n),经过主通道传递函数P(z)后形成误差传感器处输入的原始噪声信号d(n-D),其中D表示该信号相较于x(n)的时延所对应的采样点数;
x(n)输入序列延时n个采样间隔形成Z-n,而后t时刻输入的两路序列做相关运算Rxd(t);
3.根据权利要求2所述的空间降噪系统中主通道时延估计方法,其特征在于:加入时延估计部分以及延时器单元拓展了系统对主通道时延的追踪性,并须具备以下条件:
D>Δ
即主通道时延D与主次通道时延差d均为正值;
对于时延估计部分,在当前时刻t分别计算时延从1到Tc之间的相关矩阵Rxd(t):
其中,Rpp表示信号p的自相关函数,Rpq表示信号p与信号q的互相关函数,假设信号x(n-t)与信号d(n-D)中的噪声互不相关,则
Rxd(t)=Rxx(t-D)
由自相关函数性质:|Rss(τ-D)|≤Rss(0)
可知,当τ-D=0时,Rss(·)取得最大值,此时τ为当前时刻系统时延的估计
式中,arg为取自变量运算;
相关运算选取的序列长度L越长,估计结果越精确,但运算量会相应增大,考虑系统时延估计的刷新速率有着一定的要求,系统设计时应根据实际情况选取;
Δ为次级通道延迟所对应的采样点数;
对于噪声抵消部分,滤波器W(n)的输出u(n):
u(n)=ωT(n)x(n)
其中,ω(n)=[ω0(n) ω1(n) ω2(n) … ωN-1(n)]T为在时刻n的长度为N的滤波器权系数向量,x(n)=[x(n) x(n-1) x(n-2) … x(n-N+1)]T为时刻n的长度为N的输入信号向量;
u(n)经过经过次级通道S(z)得到系统输出y(n):
y(n-Δ)=s(n)*u(n)
其中,s(n)为S(z)的脉冲响应,*表示卷积运算,y(n-Δ)为由次级声源输出为抵消波形;
误差信号e(n)表示误差传感器接收到的抵消剩余信号:
计算滤波x信号xf(n):
更新滤波器W(n)的权系数ω(n):
ω(n+1)=ω(n)-μe(n)xf(n)
其中,μ为步长因子,其收敛范围如下:
其中,N为滤波器长度,Px为滤波器输入信号功率;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911347173.XA CN111161700B (zh) | 2019-12-24 | 2019-12-24 | 一种空间降噪系统中主通道时延估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911347173.XA CN111161700B (zh) | 2019-12-24 | 2019-12-24 | 一种空间降噪系统中主通道时延估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111161700A true CN111161700A (zh) | 2020-05-15 |
CN111161700B CN111161700B (zh) | 2023-03-31 |
Family
ID=70558374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911347173.XA Active CN111161700B (zh) | 2019-12-24 | 2019-12-24 | 一种空间降噪系统中主通道时延估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111161700B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115296784A (zh) * | 2022-08-10 | 2022-11-04 | 横川机器人(深圳)有限公司 | 伺服驱动器低延时同步方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080144853A1 (en) * | 2006-12-06 | 2008-06-19 | Sommerfeldt Scott D | Secondary Path Modeling for Active Noise Control |
CN105489212A (zh) * | 2014-09-17 | 2016-04-13 | 中兴通讯股份有限公司 | 下行主动降噪装置、方法及移动终端 |
CN106782594A (zh) * | 2016-12-14 | 2017-05-31 | 南京航空航天大学 | 一种抑制窄带主动噪声控制中频率不匹配的方法 |
-
2019
- 2019-12-24 CN CN201911347173.XA patent/CN111161700B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080144853A1 (en) * | 2006-12-06 | 2008-06-19 | Sommerfeldt Scott D | Secondary Path Modeling for Active Noise Control |
CN105489212A (zh) * | 2014-09-17 | 2016-04-13 | 中兴通讯股份有限公司 | 下行主动降噪装置、方法及移动终端 |
CN106782594A (zh) * | 2016-12-14 | 2017-05-31 | 南京航空航天大学 | 一种抑制窄带主动噪声控制中频率不匹配的方法 |
Non-Patent Citations (2)
Title |
---|
JOHANNES N. STROHM ET AL.: "A Fast Convergence FxLMS Algorithm for Vibration Damping of a Quarter Car", 《2018 IEEE CONFERENCE ON DECISION AND CONTROL (CDC)》 * |
张昊天: "基于主通道时延估计的变压器有源降噪研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115296784A (zh) * | 2022-08-10 | 2022-11-04 | 横川机器人(深圳)有限公司 | 伺服驱动器低延时同步方法 |
CN115296784B (zh) * | 2022-08-10 | 2023-10-20 | 横川机器人(深圳)有限公司 | 伺服驱动器低延时同步方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111161700B (zh) | 2023-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH1152988A (ja) | 適応アレイの制御方法および適応アレイ装置 | |
CN110402540B (zh) | 主动降噪方法、装置、芯片、主动控制系统和存储介质 | |
CN110265054B (zh) | 语音信号处理方法、装置、计算机可读存储介质和计算机设备 | |
CN106301288B (zh) | 一种自适应滤波方法 | |
CN105721729B (zh) | 基于块稀疏成比例重用权系数仿射投影的回声消除方法 | |
JP2018531555A (ja) | ビーム形成用途のための適応的位相歪曲のない振幅応答等化 | |
CN109379652B (zh) | 一种耳机有源噪声控制的次级通道离线辨识方法 | |
CN111161700B (zh) | 一种空间降噪系统中主通道时延估计方法 | |
CN113241053A (zh) | 简化窄带无次级路径建模有源控制方法 | |
WO2007124117A1 (en) | Methods and systems for reducing acoustic echoes in communication systems | |
CN108983144B (zh) | 改进维纳滤波器及基于该滤波器进行目标方位的估计方法 | |
EP1425853A1 (en) | Method and apparatus for generating a set of filter coefficients for a time updated adaptive filter | |
CN109147753B (zh) | 基于误差平方与误差平方对数之差最小的凸组合降噪方法 | |
CN113676157B (zh) | 一种抵消直流偏置的自适应滤波算法 | |
CN102227096A (zh) | 一种非高斯环境下的变步长最小p-范数系统辨识方法 | |
CN112881447B (zh) | 一种扫描电子显微镜真空腔室内电磁干扰消除系统及方法 | |
Chen et al. | The stability and convergence characteristics of the delayed-x LMS algorithm in ANC systems | |
Lampl | Implementation of adaptive filtering algorithms for noise cancellation | |
TWI579833B (zh) | 信號處理裝置與信號處理方法 | |
Chen et al. | Evaluation of the convergence characteristics of the filtered-x LMS algorithm in the frequency domain | |
CN112385247A (zh) | 声学延迟估计 | |
TWI385941B (zh) | 干擾消除裝置及其方法 | |
Wu et al. | Dynamic discrimination of convergence of the LMS time delay estimation in complicated noisy environments | |
CN108877829A (zh) | 一种信号处理的方法和装置 | |
CN113381730B (zh) | 一种鲁棒性自适应滤波系统 |
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 |