CN113096629A - 用于单通道反馈有源噪声控制***的相对路径虚拟传感方法 - Google Patents
用于单通道反馈有源噪声控制***的相对路径虚拟传感方法 Download PDFInfo
- Publication number
- CN113096629A CN113096629A CN202110234745.4A CN202110234745A CN113096629A CN 113096629 A CN113096629 A CN 113096629A CN 202110234745 A CN202110234745 A CN 202110234745A CN 113096629 A CN113096629 A CN 113096629A
- Authority
- CN
- China
- Prior art keywords
- virtual
- relative
- model
- microphone
- physical
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000012549 training Methods 0.000 claims abstract description 38
- 238000005316 response function Methods 0.000 claims abstract description 27
- 230000009467 reduction Effects 0.000 claims abstract description 24
- 230000003044 adaptive effect Effects 0.000 claims abstract description 18
- 238000012546 transfer Methods 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000012544 monitoring process Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 208000006735 Periostitis Diseases 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 210000003460 periosteum Anatomy 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 206010042602 Supraventricular extrasystoles Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000012814 acoustic material Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003584 silencer Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
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/1785—Methods, e.g. algorithms; Devices
- G10K11/17853—Methods, e.g. algorithms; Devices of the filter
- G10K11/17854—Methods, e.g. algorithms; Devices of the filter the filter being an adaptive filter
-
- 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/17813—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 acoustic paths, e.g. estimating, calibrating or testing of transfer functions or cross-terms
- G10K11/17817—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 acoustic paths, e.g. estimating, calibrating or testing of transfer functions or cross-terms between the output signals and the error signals, i.e. secondary path
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明公开了一种用于单通道反馈有源噪声控制***的相对路径虚拟传感方法,属于有源噪声控制技术领域。本发明包括:训练步骤:将虚拟麦克风临时放置在期望的降噪位置,物理麦克风放置在远离有源静区的位置,建立并训练物理麦克风与虚拟麦克风之间的传递模型,包括相对初级路径模型和相对次级路径模型,用以获取相对初级路径模型和相对次级路径模型的脉冲响应函数;控制步骤:在目标有源静区不放置在误差麦克风的情况下,自适应控制滤波器基于物理麦克风的实时输出、相对初级路径模型和相对次级路径模型的脉冲响应函数进行权系数的迭代更新,直至目标虚拟位置处形成有源静区。本发明可以避免噪声源无法定位的问题,又能够减少模型误差。
Description
技术领域
本发明属于有源噪声控制技术领域,具体涉及一种用于单通道反馈有源噪声控制***的 相对路径虚拟传感技术。
背景技术
传统的噪声控制属于无源噪声控制(passive noise control,PNC),主要的技术手段包 括吸声处理、隔声处理、***等,它的噪声控制机理在于使噪声声波与声学材料或结构相 互作用从而达到降低噪声的目的。当噪声信号为低频信号时,PNC在成本、部署难度和有效 性上均有局限。因此,有源噪声控制(Active noise control,ANC)的提出与发展填补上了 这一部分的空白,近年来已发展成为噪声控制中的主要研究方向。
ANC是利用声波的叠加特性来降低噪声的技术。当控制源辐射出与噪声源声波幅度相同、 相位相反的“抗噪声波”时,ANC***就可以在放置误差麦克风的期望位置处降低噪声。在 该***中,误差麦克风的作用是监测误差信号并将该信号不断地反馈给采用了自适应控制滤 波器的ANC***,此处采用自适应控制滤波器的目的是为了解决噪声信号的时变问题,不断 跟踪噪声声波的幅度、相位变化,直至放置误差麦克风的位置处形成有源静区(超过10dB 的降噪区域)。采用了自适应控制滤波器的ANC***具有实现简单、稳定性好的优点,因而在 有源噪声控制中被广泛使用。训练自适应控制滤波器常采用最小均方算法(Least mean square, LMS)。
根据使用麦克风的不同,ANC***在结构上分为前馈控制和反馈控制。前馈ANC*** 包含参考麦克风和误差麦克风,其中参考麦克风靠近噪声源并采集噪声声波的先验信息,以 便为ANC控制器提供与误差麦克风获得的误差信号高度相关的参考信号。因此,前馈ANC ***更适合处理宽带噪声,***性能更加稳定。反馈ANC***仅由误差麦克风组成,通常使 用内部模型控制噪声,该内部模型利用误差信号预测原始噪声以得到估计的参考信号。因此, 使用内部模型的反馈ANC***只能降低可预测噪声(包括单频噪声和窄带噪声)。
ANC***也可以根据使用扬声器的数量分类。单通道ANC***由一个扬声器组成,通常 与一个误差麦克风、至多一个参考麦克风一起使用。多通道ANC***包含多个扬声器和多个 麦克风,它可以形成一个相对较大的有源静区。本发明只考虑在单通道ANC***上的实现。
在一些应用场景中,误差麦克风无法放置在所需的降噪位置处(虚拟位置),比如说无线 蓝牙耳机的应用,由于无法在希望得到有源静区的人耳骨膜处放置麦克风,因而就需要利用 虚拟传感技术在人耳的骨膜处形成有源静区。该技术主要包括训练阶段和控制阶段。在训练 阶段,需要事先在虚拟位置处放置误差麦克风,得到表示监测位置与虚拟位置之间传递关系 的先验信息。在控制阶段,将原本放置在虚拟位置处的误差麦克风移除,并利用在训练阶段 得到的先验信息使虚拟位置处形成有源静区。目前,存在两种已有的虚拟传感技术,它们分 别是辅助滤波器方法(the auxiliary filter,AF)和远程麦克风方法(remote microphone, RM)。AF方法的训练阶段较为复杂,主要由两部分组成。首先在虚拟位置处放置误差麦克风, 利用FxLMS(filtered-x least mean square)算法获得此时ANC***稳定后的滤波器;然 后利用该滤波器训练包含传递关系先验信息的辅助滤波器。RM方法可以在训练阶段估计物理 初级通路和虚拟初级通路之间的关系,从而能够在控制阶段根据物理麦克风处的干扰信号估 计期望位置处的干扰信号,继而得到虚拟位置处的有源静区。Edamoto S,Shi C,Kajikawa Y于2016年发表的论文《Virtual sensing techniquefor feedforward active noise control》 和Jung W,Elliott S J,Cheer J于2018年发表的论文《Estimation of the pressure at a listener's ears in an activeheadrest system using the remote microphone technique》分别在前馈ANC***中评估了AF方法、RM法在虚拟位置处的降噪效果,从论文 的实验结果得知这两种方法均可以在不能放置误差麦克风的虚拟位置处获得超过10dB的降 噪量。
目前虚拟传感技术较多的应用于前馈ANC***中,但在噪声源无法定位的情况下,不能 确定参考麦克风的位置,因此无法使用前馈ANC***。相对AF方法,RM方法的训练阶段较 为简单,但其需要在控制阶段使用已离线建模好的虚拟次级通路模型。当声场环境变化造成 虚拟次级通路改变时,由于路径无法在不能放置误差麦克风的位置处进行在线建模,因此虚 拟次级通路模型的脉冲响应函数会和真实路径之间产生偏差,从而在***中引入模型误差, 影响***的降噪效果。
发明内容
本发明的发明目的在于:为了解决RM方法的缺点的同时降低前馈ANC***的规模,本 发明提出了一种用于单通道反馈ANC***的相对路径虚拟传感(Relative Path,RP)方法。
本发明的用于单通道反馈有源噪声控制***的相对路径虚拟传感方法,包括下列步骤:
训练步骤:将虚拟麦克风临时放置在期望的降噪位置,物理麦克风(即监测麦克风)放 置在远离有源静区的位置,建立并训练物理麦克风与虚拟麦克风之间的传递模型,包括相对 初级路径模型和相对次级路径模型,用以获取相对初级路径模型的脉冲响应函数和相对次级 路径模型的脉冲响应函数;
其中,相对初级路径模型指:物理初级通路与虚拟初级通路之间的模型,相对次级路径 模型指:物理次级通路与虚拟次级通路之间的模型;
控制步骤:在目标有源静区不放置在误差麦克风的情况下,自适应控制滤波器基于物理 麦克风的实时输出和训练步骤获得的相对初级路径模型的脉冲响应函数和相对次级路径模型 的脉冲响应函数进行权系数的迭代更新,直至目标虚拟位置处形成有源静区。
在一种可能的实现方式中,自适应控制滤波器为噪声控制滤波器。
在一种可能的实现方式中,训练步骤中,对所述传递模型的训练过程为:
对所建立的相对初级路径模型和相对次级路径模型进行迭代更新,当相对初级路径模型 中的虚拟位置处的误差信号和相对次级路径模型中的虚拟位置处的误差信号满足预置的收敛 条件时,停止迭代更新,并基于当前迭代更新后的初级路径模型和相对次级路径模型得到对 应的脉冲响应函数。
在一种可能的实现方式中,控制步骤中,自适应控制滤波器进行更新迭代的处理过程具 体为:
基于相对初级路径模型的脉冲响应函数和相对次级路径模型的脉冲响应函数获取虚拟位 置估计处的干扰信号和抗噪信号,两者的叠加得到虚拟位置处估计的误差信号,自适应控制 滤波器通过该误差信号和滤波参考信号进行迭代更新。
综上所述,由于采用了上述技术方案,本发明的有益效果是:既可以避免噪声源无法定 位的问题,又能够减少由于不能对次级通路进行在线建模所导致的模型误差。
附图说明
图1是具体实施方式中,本发明提出的用于单通道反馈ANC***的相对路径虚拟传感方 法的过程框图,其中图1(a)是训练阶段,物理初级通道域虚拟初级通路之间的关系图;图 1(b)是训练阶段,物理次级通路与虚拟次级通路之间的关系图;图1(c)是控制阶段的过 程框图;
图2是具体实施方式中,本发明在单通道反馈ANC***上的结构部署示意图;
图3是具体实施方式中,对应图2所示的部署结构的声场路径设置图;
图4是具体实施方式中,本发明在训练阶段的误差信号图;
图5是具体实施方式中,本发明在控制阶段的误差信号图;
图6是具体实施方式中,本发明在控制阶段的降噪效果图;
图7是具体实施方式中,本发明在控制阶段收敛后的噪声功率谱图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作 进一步地详细描述。
当虚拟次级通路发生变化时,由于控制阶段不能在虚拟位置处放置临时麦克风,因此无 法用在线通路建模技术获得RM方法所需要的精确次级通路模型。为了解决RM方法的缺点的 同时降低前馈ANC***的规模,本发明提出了用于单通道反馈ANC***的相对路径虚拟传感 方法。
RP方法的第一阶段是训练阶段。在这一阶段,为了建立物理麦克风与误差麦克风之间的 传递模型,即物理初级通路与虚拟初级通路之间的模型、物理次级通路与虚拟次级通路之间 的模型,虚拟麦克风临时放置在期望的降噪位置处。物理麦克风,也称为监测麦克风,放置 在远离有源静区的位置。即监测麦克风与有源静区的距离大于某个指定的距离门限。
RP方法的第二阶段是控制阶段,在没有任何误差麦克风放置在目标有源静区的情况下, 自适应控制滤波器利用物理麦克风的实时输出和训练阶段获得的先验信息进行更新迭代,直 至虚拟位置处形成有源静区。
图1展示了本发明的相对路径虚拟传感方法的过程框图,RP方法的训练阶段包括两个部 分。在训练阶段(a)中(参考图1(a))),利用物理麦克风测量得到的干扰信号去估计虚拟 麦克风处的干扰信号,这种估计依赖于相对初级通路模型Cp(z),可表示为:
其中,下标m、v分别表示物理位置和虚拟位置,cp(n)是相对初级路径模型的脉冲响应 函数,x(n)是噪声源信号的矢量形式,*表示卷积运算。在公式(1)中,dm(n)、dv(n)分别是当前时刻物理麦克风和虚拟麦克风测得的干扰信号,它们的矢量形式分别表示为:
dv(n)=[dv(n),dv(n-1),...,dv(n-N+1)]T (2)
dm(n)=[dm(n),dm(n-1),...,dm(n-N+1)]T (3)
其中,N是采样长度,取决于相对初级路径模型Cp(z)、相对次级路径模型Cs(z)的长度。 n表示采样点,物理初级通路pm(n)、虚拟初级通路pv(n)分别是噪声源到物理麦克风和虚拟 麦克风之间的脉冲响应函数,上标T表示矩阵转置。
在训练阶段(b)中(参考图1(b)),相对次级通路模型Cs(z)能够根据物理麦克风测量 得到的抗噪信号估计出虚拟麦克风处的抗噪信号,因此有:
其中,y(n)是控制源信号的矢量表示,cs(n)是相对次级路径模型的脉冲响应函数。ym(n)、 yv(n)分别表示物理麦克风和虚拟麦克风测得的抗噪信号,物理次级通路sm(n)、虚拟次级通 路sv(n)分别是次级扬声器到物理麦克风和虚拟麦克风之间的脉冲响应函数。
相对次级路径模型Cs(z)只与物理次级通路及虚拟次级通路有关。因此,在控制阶段可利 用物理次级通路模型及相对次级路径模型间接得到虚拟次级通路模型,避免环境发生变化时, 虚拟次级通路无法在线建模的弊端。
本具体实施方式中,利用虚拟位置处测得的误差信号ecp(n)、ecs(n)来判断训练阶段是否 达到稳定状态,它们分别表示为:
ecp(n)=dv(n)-dm(n)*cp(n) (5)
ecs(n)=yv(n)-ym(n)*cs(n) (6)
当误差信号ecp(n)、ecs(n)随着相对路径模型的更新迭代逐渐收敛并保持收敛状态,那么 表示训练阶段逐渐达到稳定状态并保持稳定状态。
在控制阶段(参考图1(c)),图1中所示的带背景色的方框(灰色方框)分别表示相对 路径模型,其中,是物理次级通路模型,假设即利用训练 阶段得到的cp(n)和cs(n),可以估计得到虚拟位置处的误差信号:
其中,是当前时刻估计的虚拟位置处的干扰信号和抗噪信号。在公式(7) 中,是虚拟次级通路模型的脉冲响应函数,假设通过该 式得知虚拟次级通路模型可以通过cs(n)间接获得,由于可以在物理位置处进行在线建模 来获取,因此,当声场环境改变虚拟次级通路时,可以通过对sm(n)进行在线建模得到最优的 虚拟次级通路模型
在公式(8)中,em(n)是物理麦克风检测得到的误差信号,可表示为:
em(n)=dm(n)+sm(n)*y(n) (9)
同样地,可以得到虚拟位置处的误差信号ev(n),可表示为:
ev(n)=dv(n)+sv(n)*y(n) (10)
当误差信号ev(n)随着噪声控制滤波器的更新迭代逐渐收敛并保持收敛状态,那么表示控 制阶段逐渐达到稳定状态并保持稳定。
在控制阶段,wRP(z)是长度为Lw的噪声控制滤波器,它输出的控制源信号y(n)可表示为:
y(n)=wRP(n)Tx(n) (11)
其中,wRP(n)=[wRP(n),wRP(n-1),...,wRP(n-Lw+1)]T是控制滤波器的权系数矢量,T表示转 置运算。y(n)的矢量表示为:
y(n)=[y(n),y(n-1),...,y(n-Ls+1)]T (12)
其中,Ls是取决于物理次级通路模型的长度。
因此,在该反馈***中,可利用em(n)估计得到噪声控制滤波器wRP(z)的输入参考信号 x(n):
它的矢量形式表示为:
x(n)=[x(n),x(n-1),...,x(n-Lw+1)]T (14)
因此,噪声控制滤波器可以通过下式进行更新迭代:
r(n)=[r(n),r(n-1),...,r(n-Lw+1)]T (16)
综上所述,RP方法的控制阶段通过利用训练阶段得到的cp(n)和cs(n),可以得到虚拟位 置处估计的干扰信号和抗噪信号进而通过声波的叠加特性得到虚拟位置处估计 的误差信号,自适应控制滤波器WRP(z)通过和滤波参考信号r(n)进行迭代更新,直至虚 拟位置处形成有源静区。可通过下式观察虚拟位置处的降噪效果:
即,通过式(17)计算降噪效果Reduction,当其大于或等于指定阈值(例如0、10等,可以基于实际应用场景进行设置,本发明不做具体限定)时,则确认当前虚拟位置形成有源静区。
设置实验所需参数。图2中展示了本发明的RP方法在单通道反馈ANC***上的结构部署, 图3是其声场路径设置,由实时DSP(Digital Signal Processing)平台测量所需全部路径 的脉冲响应。其中,***的采样率是16kHz;物理初级通路、虚拟初级通路的脉冲响应长度 为75ms;物理次级通路、虚拟次级通路的脉冲响应长度为25ms。将噪声控制滤波器、相对路 径模型的长度设置为400阶。本发明的RP方法在两个阶段所用的自适应算法均为归一化FxLMS 算法,收敛步长设置可参考表1。
表1 RP方法在反馈ANC***上的收敛部长设置
由于反馈ANC***能够降低窄带噪声或单频噪声。因此,本发明将噪声源作如下设置:
首先,设置噪声源是带宽为600-800Hz的窄带高斯白噪声,引入宽带干扰,设置信噪比 为64。在训练阶段,把虚拟麦克风临时放置在有源静区处,通过观察误差信号ecp(n)、ecs(n)的 收敛情况,判断训练阶段是否达到稳定状态。
图4展示了当噪声源是窄带信号时,训练阶段的误差信号图。观测图3可知,训练阶段 在t=200s后达到稳定状态。因此,可得到相对路径模型cp(n)和cs(n)的滤波器系数矩阵。随 后将虚拟位置处的误差麦克风撤离,利用训练阶段得到的cp(n)和cs(n)估计虚拟位置处的干扰 信号和抗噪信号,从而能够在控制阶段得到虚拟位置处的误差信号,通过式(15)对滤波器 进行更新迭代,直至虚拟位置处形成有源静区。
通过图5,可观察到控制阶段在t=20s后达到收敛状态,即***在t=20s后达到稳定状态。 根据式(17),可以计算出在单通道反馈ANC***中使用RP方法的降噪效果,如图6所示。 从图6可得知:虚拟位置处的降噪效果达到了15dB左右,因此RP算法能够在虚拟位置处形 成有源静区。而在物理位置处,噪声增加了约2dB左右,因此物理位置不适合作为降噪区域。 同样的结论可通过图7所示的噪声功率谱图观察到,通过图7,可得知,在相应频段内虚拟位 置处的噪声得到了处理,但是降噪后的噪声功率并没有低于底噪,符合真实场景(任何环境 中都存在底噪,降噪后的噪声功率不低于底噪)。
综上所述,在单通道反馈ANC***中使用RP方法可以在期望的降噪位置处形成有源静区, 同时该方法除了可以避免噪声源无法定位的问题外还能够减少虚拟次级通路无法在线建模所 带来的模型误差。因此,可考虑该方法可在具体应用中的实现。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述, 均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过 程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (7)
1.用于单通道反馈有源噪声控制***的相对路径虚拟传感方法,其特征在于,包括下列步骤:
训练步骤:将虚拟麦克风临时放置在期望的降噪位置,物理麦克风放置在远离有源静区的位置,建立并训练物理麦克风与虚拟麦克风之间的传递模型,所述传递模包括相对初级路径模型和相对次级路径模型,用以获取相对初级路径模型的脉冲响应函数和相对次级路径模型的脉冲响应函数;
其中,相对初级路径模型指:物理初级通路与虚拟初级通路之间的模型,相对次级路径模型指:物理次级通路与虚拟次级通路之间的模型;
控制步骤:在目标有源静区不放置误差麦克风的情况下,自适应控制滤波器基于物理麦克风的实时输出和训练步骤获得的相对初级路径模型的脉冲响应函数和相对次级路径模型的脉冲响应函数进行权系数的迭代更新,直至目标虚拟位置处形成有源静区。
2.如权利要求1所述的方法,其特征在于,所述自适应控制滤波器为噪声控制滤波器。
3.如权利要求1所述的方法,其特征在于,训练步骤中,对所述传递模型的训练过程为:
对所建立的相对初级路径模型和相对次级路径模型进行迭代更新,当相对初级路径模型中的虚拟位置处的误差信号和相对次级路径模型中的虚拟位置处的误差信号满足预置的收敛条件时,停止迭代更新,并基于当前迭代更新后的初级路径模型和相对次级路径模型得到对应的脉冲响应函数。
4.如权利要求1所述的方法,其特征在于,控制步骤中,自适应控制滤波器进行更新迭代的处理过程具体为:
基于相对初级路径模型的脉冲响应函数和相对次级路径模型的脉冲响应函数获取虚拟位置估计处的干扰信号和抗噪信号,两者的叠加得到虚拟位置处估计的误差信号,自适应控制滤波器通过该误差信号和滤波参考信号进行迭代更新。
5.如权利要求1至4任一项所述的方法,其特征在于,所述相对初级通路模型为:
其中,下标m、v分别表示物理位置和虚拟位置,cp(n)表示相对初级路径模型的脉冲响应函数,x(n)表示噪声源信号的矢量形式,物理初级通路pm(n)、虚拟初级通路pv(n)分别是噪声源到物理麦克风和虚拟麦克风之间的脉冲响应函数,*表示卷积运算,dm(n)、dv(n)分别表示当前时刻物理麦克风和虚拟麦克风测得的干扰信号,其矢量形式分别表示为:
dv(n)=[dv(n),dv(n-1),...,dv(n-N+1)]T
dm(n)=[dm(n),dm(n-1),...,dm(n-N+1)]T
其中,N表示采样长度;
所述相对次级通路模型为:
其中,y(n)表示控制源信号的矢量表示,cs(n)表示相对次级路径模型的脉冲响应函数。ym(n)、yv(n)分别表示物理麦克风和虚拟麦克风测得的抗噪信号,物理次级通路sm(n)、虚拟次级通路sv(n)分别是次级扬声器到物理麦克风和虚拟麦克风之间的脉冲响应函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110234745.4A CN113096629B (zh) | 2021-03-03 | 2021-03-03 | 用于单通道反馈有源噪声控制***的相对路径虚拟传感方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110234745.4A CN113096629B (zh) | 2021-03-03 | 2021-03-03 | 用于单通道反馈有源噪声控制***的相对路径虚拟传感方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113096629A true CN113096629A (zh) | 2021-07-09 |
CN113096629B CN113096629B (zh) | 2022-11-04 |
Family
ID=76668055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110234745.4A Expired - Fee Related CN113096629B (zh) | 2021-03-03 | 2021-03-03 | 用于单通道反馈有源噪声控制***的相对路径虚拟传感方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113096629B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113763917A (zh) * | 2021-09-30 | 2021-12-07 | 展讯通信(上海)有限公司 | 滤波器处理方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190035380A1 (en) * | 2017-07-28 | 2019-01-31 | Harman Becker Automotive Systems Gmbh | Silent zone generation |
CN109448686A (zh) * | 2018-12-13 | 2019-03-08 | 重庆邮电大学 | 基于次级通道在线辨识新算法交叉更新有源噪声控制*** |
CN110599996A (zh) * | 2019-09-09 | 2019-12-20 | 电子科技大学 | 一种利用上位机实时离线更新的有源噪声控制*** |
CN110719550A (zh) * | 2019-10-21 | 2020-01-21 | 南京南大电子智慧型服务机器人研究院有限公司 | 一种双通道有源降噪头靠的虚拟传声器优化设计方法 |
CN111031442A (zh) * | 2019-12-31 | 2020-04-17 | 湖南景程电子科技有限公司 | 一种基于dsp虚拟传感的耳机主动降噪自适应*** |
CN111583896A (zh) * | 2020-05-13 | 2020-08-25 | 苏州静声泰科技有限公司 | 一种多通道有源降噪头靠的降噪方法 |
CN112102806A (zh) * | 2020-09-06 | 2020-12-18 | 西安艾科特声学科技有限公司 | 列车驾驶室基于虚拟传感技术的主动噪声控制***及方法 |
-
2021
- 2021-03-03 CN CN202110234745.4A patent/CN113096629B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190035380A1 (en) * | 2017-07-28 | 2019-01-31 | Harman Becker Automotive Systems Gmbh | Silent zone generation |
CN109448686A (zh) * | 2018-12-13 | 2019-03-08 | 重庆邮电大学 | 基于次级通道在线辨识新算法交叉更新有源噪声控制*** |
CN110599996A (zh) * | 2019-09-09 | 2019-12-20 | 电子科技大学 | 一种利用上位机实时离线更新的有源噪声控制*** |
CN110719550A (zh) * | 2019-10-21 | 2020-01-21 | 南京南大电子智慧型服务机器人研究院有限公司 | 一种双通道有源降噪头靠的虚拟传声器优化设计方法 |
CN111031442A (zh) * | 2019-12-31 | 2020-04-17 | 湖南景程电子科技有限公司 | 一种基于dsp虚拟传感的耳机主动降噪自适应*** |
CN111583896A (zh) * | 2020-05-13 | 2020-08-25 | 苏州静声泰科技有限公司 | 一种多通道有源降噪头靠的降噪方法 |
CN112102806A (zh) * | 2020-09-06 | 2020-12-18 | 西安艾科特声学科技有限公司 | 列车驾驶室基于虚拟传感技术的主动噪声控制***及方法 |
Non-Patent Citations (3)
Title |
---|
HASEEB ABDUL: "" A robust approach for online feedback path modeling in single-channel narrow-band active noise control systems using two distinct variable step size methods"", 《APPLIED ACOUSTICS》 * |
LIU L: ""Virtual sensing techniques and their applications"", 《2009 INTERNATIONAL CONFERENCE ON NETWORKING, SENSING AND CONTROL》 * |
SHI C: ""An active noise control casing using the multi-channel feedforward control system and the relative path based virtual sensing method"", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113763917A (zh) * | 2021-09-30 | 2021-12-07 | 展讯通信(上海)有限公司 | 滤波器处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113096629B (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108428445B (zh) | 一种无误差传声器的自适应主动降噪方法 | |
EP0836736B1 (en) | Digital feed-forward active noise control system | |
US9633645B2 (en) | Adaptive noise control system with improved robustness | |
US10373600B2 (en) | Active noise control system | |
JP5629372B2 (ja) | 聴取者に対する環境雑音の効果を低減させる方法および装置 | |
US7885417B2 (en) | Active noise tuning system | |
US20160300562A1 (en) | Adaptive feedback control for earbuds, headphones, and handsets | |
EP1003154A2 (en) | Acoustic system identification using acoustic masking | |
JPH09501779A (ja) | 適応フィードフォワード及びフィードバック制御装置 | |
Shen et al. | A wireless reference active noise control headphone using coherence based selection technique | |
CN109346052B (zh) | 一种利用主动降噪优化车内声品质的装置及方法 | |
CN113257214A (zh) | 一种用于风扇管道***的主动降噪方法 | |
CN113096629B (zh) | 用于单通道反馈有源噪声控制***的相对路径虚拟传感方法 | |
CN114127845A (zh) | 用于消除麦克风信号中的道路噪声的***和方法 | |
Kajikawa et al. | Comparison of virtual sensing techniques for broadband feedforward active noise control | |
EP1074971B1 (en) | Digital feed-forward active noise control system | |
Huang et al. | Implementation of feedforward active noise control techniques for headphones | |
Fan et al. | Robust adaptive algorithms for active noise and vibration control | |
Sachau et al. | Real-time implementation of the frequency-domain FxLMS algorithm without block delay for an adaptive noise blocker | |
US11948546B2 (en) | Feed-forward adaptive noise-canceling with dynamic filter selection based on classifying acoustic environment | |
US20240038211A1 (en) | Adaptive noise-canceling with dynamic filter selection based on multiple noise sensor signal phase differences | |
Fabry et al. | Online secondary path estimation with masked auxiliary noise for active noise control | |
JP4906787B2 (ja) | 能動型振動騒音制御装置 | |
Nagi et al. | Active noise cancellation with TMS320C5402 DSP starter kit | |
Narahari | Noise cancellation in headphones |
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 | ||
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: 20221104 |