CN115565515A - 一种分步虚拟传感降噪方法 - Google Patents

一种分步虚拟传感降噪方法 Download PDF

Info

Publication number
CN115565515A
CN115565515A CN202211186487.8A CN202211186487A CN115565515A CN 115565515 A CN115565515 A CN 115565515A CN 202211186487 A CN202211186487 A CN 202211186487A CN 115565515 A CN115565515 A CN 115565515A
Authority
CN
China
Prior art keywords
filter
physical
virtual
microphone
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.)
Pending
Application number
CN202211186487.8A
Other languages
English (en)
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.)
Nanjing University
Original Assignee
Nanjing 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 Nanjing University filed Critical Nanjing University
Priority to CN202211186487.8A priority Critical patent/CN115565515A/zh
Publication of CN115565515A publication Critical patent/CN115565515A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods 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/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods 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/1781Methods 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
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods 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/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods 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/1781Methods 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/17813Methods 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/17817Methods 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
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods 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/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods 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/1785Methods, e.g. algorithms; Devices
    • G10K11/17853Methods, e.g. algorithms; Devices of the filter
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods 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/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods 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/1787General system configurations

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***即可在放置误差麦克风的期望位置处降低噪声。
根据使用麦克风的不同,ANC降噪***在结构上分为前馈式和反馈式。前馈式***性能稳定,适合处理宽带噪声;反馈式***只能降低单频噪声、窄带噪声,具有局限性。根据使用扬声器数量的不同,ANC降噪***还可分为单通道ANC***和多通道ANC***。
在一些应用场景,如无线蓝牙耳机,目标降噪位置处无法放置误差麦克风,虚拟传感降噪技术应运而生。虚拟传感降噪技术主要分为两个阶段:训练阶段和控制阶段。在训练阶段,需暂时在目标降噪位置处放置虚拟麦克风,得到物理位置与虚拟位置之间的传递关系;在控制阶段,目标降噪位置的虚拟麦克风被移除,根据物理麦克风采集到的信号与在训练阶段得到的物理位置与虚拟位置之间的传递关系对噪声控制滤波器更新,最终达到稳定状态,从而实现目标降噪位置的有源降噪。目前,已存在两种虚拟传感技术:远程麦克风方法(remote microphone,RM)和辅助滤波器方法(auxiliary filter,AF)。对于AF方法,训练阶段比控制阶段多使用一个麦克风,对硬件控制器的接口数和运算量均有较大要求,为此本发明旨在简化AF算法。
发明内容
本发明的发明目的在于:为了解决AF方法训练阶段比控制阶段需多使用一个麦克风,硬件资源消耗大的缺点,同时为了当参考信号音量范围较大时,目标降噪位置降噪均能达到指定阈值,本发明提出了一种分步虚拟传感降噪方法,包括如下步骤:
步骤一:建模虚拟次级通路Sv,将虚拟麦克风临时放置在目标降噪位置采集虚拟误差信号,根据LMS算法,次级扬声器发出的次级声源与虚拟麦克风接收到的虚拟误差信号共同作用更新虚拟次级通路滤波器
Figure BSA0000285341440000021
最终达到稳定状态,得到虚拟次级通路滤波器
Figure BSA0000285341440000022
步骤二:训练噪声控制滤波器W,根据FxLMS算法,虚拟麦克风接收到的虚拟误差信号为参考信号通过噪声控制滤波器W的输出和初级噪声的叠加,参考信号通过虚拟次级通路滤波器
Figure BSA0000285341440000029
的输出与虚拟麦克风接收到的虚拟误差信号共同作用更新噪声控制滤波器W,最终达到稳定状态,得到噪声控制滤波器W;
步骤三:建模物理次级通路Sm,将步骤一、二中的虚拟麦克风作为物理麦克风放置在物理位置采集物理误差信号,根据LMS算法,次级扬声器发出的次级声源与物理麦克风接收到的物理误差信号共同作用更新物理次级通路滤波器
Figure BSA0000285341440000023
最终达到稳定状态,得到物理次级通路滤波器
Figure BSA0000285341440000024
步骤四:训练辅助滤波器H,目标降噪位置不放置任何麦克风,物理位置放置物理麦克风,将步骤二训练得到的噪声控制滤波器W作为固定系数滤波器W0,物理麦克风接收到的物理误差信号为参考信号通过固定系数滤波器W0和物理次级通路滤波器
Figure BSA0000285341440000025
的输出和初级噪声的叠加,辅助滤波器产生的残余信号为物理麦克风接收到的物理误差信号与参考信号通过辅助滤波器H的输出的叠加,参考信号与辅助滤波器产生的残余信号共同作用更新辅助滤波器H,最终达到稳定状态,得到辅助滤波器H;
步骤五:进入控制阶段,目标降噪位置不放置任何麦克风,物理位置放置物理麦克风,物理麦克风接收到的物理误差信号为参考信号通过噪声控制滤波器W和物理次级通路滤波器
Figure BSA0000285341440000026
的输出和初级噪声的叠加,将步骤四训练得到的辅助滤波器H作为固定系数滤波器H0,辅助滤波器产生的残余信号为物理麦克风接收到的物理误差信号与参考信号通过固定系数滤波器H0的输出的叠加,参考信号通过物理次级通路滤波器
Figure BSA0000285341440000027
的输出与辅助滤波器产生的残余信号共同作用更新噪声控制滤波器W,最终达到稳定状态,实现目标降噪位置的有源降噪。
进一步地,步骤二中,噪声控制滤波器W更新的过程为:
通过公式
Figure BSA0000285341440000028
对权系数矢量进行更新,其中w(n)表示噪声控制滤波器更新前的权系数矢量,w(n+1)表示更新后的权系数矢量,x(n)表示参考信号,
Figure BSA0000285341440000031
表示虚拟次级通路滤波器,ev(n)表示虚拟误差信号,μw1表示噪声控制滤波器更新的收敛因子。
进一步地,步骤四中,辅助滤波器H更新的过程为:
通过公式h(n+1)=h(n)-μheh(n)x(n)对权系数矢量进行更新,其中h(n)表示辅助滤波器更新前的权系数矢量,h(n+1)表示更新后的权系数矢量,x(n)表示参考信号,eh(n)表示辅助滤波器产生的残余信号,μh表示辅助滤波器更新的收敛因子。
进一步地,步骤五中,噪声控制滤波器W更新的过程为:
通过公式w(n+1)=w(n)-μw2eh(n)xs(n)对权系数矢量进行更新,其中w(n)表示噪声控制滤波器更新前的权系数矢量,w(n+1)表示更新后的权系数矢量,xs(n)表示参考信号通过物理次级通路滤波器的输出,eh(n)表示辅助滤波器产生的残余信号,μw2表示噪声控制滤波器更新的收敛因子。
进一步地,为满足初级噪声强度范围较大时的降噪需求,收敛步长分别设置为:
Figure BSA0000285341440000032
Figure BSA0000285341440000033
Figure BSA0000285341440000034
其中,μw1表示噪声控制滤波器在步骤二中更新的收敛因子,μh表示辅助滤波器在步骤四中更新的收敛因子,μw2表示噪声控制滤波器在步骤五中更新的收敛因子,x表示参考信号。
综上所述,由于采用了上述技术方案,本发明的有益效果是:既解决了AF方法训练阶段比控制阶段多使用一个麦克风,硬件资源消耗大的问题,又使当参考信号音量范围较大时,目标降噪位置降噪均能达到指定阈值。
附图说明
图1是本发明分步虚拟传感降噪方法的过程框图a)步骤二过程框图;b)步骤四过程框图;c)步骤五过程框图。
图2是具体实施方式中,本发明在参考信号音量为中时的次级通路图。
图3是具体实施方式中,本发明在参考信号音量为中时步骤二、步骤四的误差信号图。
图4是具体实施方式中,本发明在参考信号音量为中时步骤五的误差信号图。
图5是具体实施方式中,本发明在参考信号音量分别为低、中、高时步骤五的噪声功率谱图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明做进一步的详细描述。
为了解决AF方法训练阶段比控制阶段需多使用一个麦克风,硬件资源消耗大的缺点,同时为了当参考信号音量范围较大时,目标降噪位置降噪均能达到指定阈值,本发明提出了一种分步虚拟传感降噪方法。
附图1展示了本发明分步虚拟传感降噪方法的过程框图,本具体实施方式中,步骤一将虚拟麦克风临时放置在目标降噪位置采集虚拟误差信号,建模虚拟次级通路Sv,表示为:
yv(n)=sv(n)*y(n)
其中,y(n)表示次级扬声器发出的次级声源,yv(n)表示虚拟麦克风接收到的虚拟误差信号,次级扬声器发出的次级声源与虚拟麦克风接收到的虚拟误差信号共同作用更新虚拟次级通路滤波器
Figure BSA0000285341440000041
最终达到稳定状态,得到虚拟次级通路滤波器
Figure BSA0000285341440000042
本具体实施方式中,步骤二(参考附图1(a))训练噪声控制滤波器W,根据FxLMS算法,虚拟麦克风接收到的虚拟误差信号为参考信号通过噪声控制滤波器W的输出和初级噪声的叠加,表示为:
ev(n)=dv(n)+sv(n)*y(n)
其中,dv(n)表示初级噪声,ev(n)表示虚拟麦克风接收到的虚拟误差信号。参考信号通过虚拟次级通路滤波器
Figure BSA0000285341440000043
的输出与虚拟麦克风接收到的虚拟误差信号共同作用更新噪声控制滤波器W,噪声控制滤波器W更新的过程为:
Figure BSA0000285341440000044
最终达到稳定状态,得到噪声控制滤波器W。其中w(n)表示噪声控制滤波器更新前的权系数矢量,w(n+1)表示更新后的权系数矢量,x(n)表示参考信号,
Figure BSA0000285341440000045
表示虚拟次级通路滤波器,ev(n)表示虚拟误差信号,μw1表示噪声控制滤波器更新的收敛因子。
本具体实施方式中,步骤三将步骤一、二中的虚拟麦克风作为物理麦克风放置在物理位置采集物理误差信号,建模物理次级通路Sm,表示为:
ym(n)=sm(n)*y(n)
其中,y(n)表示次级扬声器发出的次级声源,ym(n)表示物理麦克风接收到的物理误差信号,次级扬声器发出的次级声源与物理麦克风接收到的物理误差信号共同作用更新物理次级通路滤波器
Figure BSA0000285341440000051
最终达到稳定状态,得到物理次级通路滤波器
Figure BSA0000285341440000052
本具体实施方式中,步骤四(参考附图1(b))目标降噪位置不放置任何麦克风,物理位置放置物理麦克风,训练辅助滤波器H,将步骤二训练得到的噪声控制滤波器W作为固定系数滤波器W0,物理麦克风接收到的物理误差信号为参考信号通过固定系数滤波器W0和物理次级通路滤波器
Figure BSA0000285341440000053
的输出和初级噪声的叠加,表示为:
Figure BSA0000285341440000054
其中,dm(n)表示初级噪声,em(n)表示物理麦克风接收到的物理误差信号。辅助滤波器产生的残余信号为物理麦克风接收到的物理误差信号与参考信号通过辅助滤波器H的输出的叠加,参考信号与辅助滤波器产生的残余信号共同作用更新辅助滤波器H,辅助滤波器H更新的过程为:
h(n+1)=h(n)-μheh(n)x(n)
最终达到稳定状态,得到辅助滤波器H。其中h(n)表示辅助滤波器更新前的权系数矢量,h(n+1)表示更新后的权系数矢量,x(n)表示参考信号,eh(n)表示辅助滤波器产生的残余信号,μh表示辅助滤波器更新的收敛因子。
本具体实施方式中,步骤五(参考附图1(c))目标降噪位置不放置任何麦克风,物理位置放置物理麦克风,进入控制阶段,物理麦克风接收到的物理误差信号为参考信号通过噪声控制滤波器W和物理次级通路滤波器
Figure BSA0000285341440000055
的输出和初级噪声的叠加,表示为:
Figure BSA0000285341440000056
其中,dm(n)表示初级噪声,em(n)表示物理麦克风接收到的物理误差信号。将步骤四训练得到的辅助滤波器H作为固定系数滤波器H0,辅助滤波器产生的残余信号为物理麦克风接收到的物理误差信号与参考信号通过固定系数滤波器H0的输出的叠加,参考信号通过物理次级通路滤波器
Figure BSA0000285341440000061
的输出与辅助滤波器产生的残余信号共同作用更新噪声控制滤波器W,噪声控制滤波器W更新的过程为:
w(n+1)=w(n)-μw2eh(n)xs(n)
最终达到稳定状态,实现目标降噪位置的有源降噪。其中w(n)表示噪声控制滤波器更新前的权系数矢量,w(n+1)表示更新后的权系数矢量,xs(n)表示参考信号通过物理次级通路滤波器的输出,eh(n)表示辅助滤波器产生的残余信号,μw2表示噪声控制滤波器更新的收敛因子。
本发明从运算量来说(以M代表一次乘法运算,A代表一次加法运算,L为滤波器长度,LM即代表L次乘法运算,默认所以滤波器长度相同):
执行一次步骤一运算量为:(2L+1)M+2LA
执行一次步骤二运算量为:(4L+1)M+(4L-2)A
执行一次步骤三运算量为:(2L+1)M+2LA
执行一次步骤四运算量为:(4L+1)M+(4L-1)A
执行一次步骤五运算量为:(6L+1)M+(6L-2)A
具体实施方式通过先录音,后将采集得到的录音在MATLAB R2016b软件进行仿真的方法来预测目标降噪位置的降噪量。实验在一个局域半封闭空间进行。噪声源为噪声录音(音量分为低、中、高三种情况),参考麦克风放置在局域半封闭空间后方,虚拟次级通路Sv、物理次级通路Sm建模时次级扬声器发出的次级声源为白噪。将局域半封闭空间左耳耳道口设为目标降噪位置,左耳耳道口、次级扬声器及物理麦克风置于同一水平线,次级扬声器位于局域半封闭空间左部,物理麦克风距离次级扬声器1cm,左耳耳道口距离物理麦克风8cm。
本实验测得的数据,所有录音均为25s。仿真时,噪声控制滤波器、辅助滤波器的长度均设置为384阶。为使初级噪声强度范围较大时,目标降噪位置降噪均能达到指定阈值,已知收敛因子的取值与参考信号能量的反比有关,本发明收敛步长分别设置为:
Figure BSA0000285341440000071
Figure BSA0000285341440000072
Figure BSA0000285341440000073
其中,μw1表示噪声控制滤波器在步骤二中更新的收敛因子,μh表示辅助滤波器在步骤四中更新的收敛因子,μw2表示噪声控制滤波器在步骤五中更新的收敛因子,x表示参考信号。
首先具体分析参考信号音量为中时数据训练及控制阶段的情况,最后再讨论该方法在局域半封闭空间应用中当参考信号音量分别为低、中、高时的降噪效果。以下所有结果都是通过MATLAB R2016b软件仿真得到:
附图2是得到的次级通路滤波器。***的采样率设置为24k,虚拟次级通路、物理次级通路的脉冲响应长度为16ms。
附图3展示了参考信号音量为中时,步骤二及步骤四的误差信号图。观测附图3可知,步骤二在t=0.5后达到稳定状态,步骤四在t=1s后达到稳定状态。
通过附图4,可观察到步骤五在t非常小时就达到收敛状态,即***在很短的时间内就达到稳定状态。
附图5展示了参考信号音量分别为低、中、高时的降噪效果。从附图5可得知:参考信号音量为低时,目标降噪位置处的降噪效果达到了20dB左右;参考信号音量为中时,目标降噪位置处的降噪效果达到了22dB左右;参考信号音量为高时,目标降噪位置处的降噪效果达到了22dB左右。因此本发明对于参考信号在较大音量范围内时,目标降噪位置均能形成有源静区。
综上所述,使用一种分步虚拟传感降噪方法可以减少使用一个麦克风,节省硬件资源,同时对于参考信号在较大音量范围内时,目标降噪位置均能形成有源静区,说明了本发明所提方法的有效性。
需要说明的是上述实施例,并非用来限定本发明的保护范围,在上述技术方案的基础上所做出的等同变换或替换及以此为基础扩充等若干改进均落入本发明权利要求所保护的范围。

Claims (5)

1.一种分步虚拟传感降噪方法,其特征在于,包括如下步骤:
步骤一:建模虚拟次级通路Sv,将虚拟麦克风临时放置在目标降噪位置采集虚拟误差信号,根据LMS算法,次级扬声器发出的次级声源与虚拟麦克风接收到的虚拟误差信号共同作用更新虚拟次级通路滤波器
Figure FSA0000285341430000011
最终达到稳定状态,得到虚拟次级通路滤波器
Figure FSA0000285341430000012
步骤二:训练噪声控制滤波器W,根据FxLMS算法,虚拟麦克风接收到的虚拟误差信号为参考信号通过噪声控制滤波器W的输出和初级噪声的叠加,参考信号通过虚拟次级通路滤波器
Figure FSA0000285341430000013
的输出与虚拟麦克风接收到的虚拟误差信号共同作用更新噪声控制滤波器W,最终达到稳定状态,得到噪声控制滤波器W;
步骤三:建模物理次级通路Sm,将步骤一、二中的虚拟麦克风作为物理麦克风放置在物理位置采集物理误差信号,根据LMS算法,次级扬声器发出的次级声源与物理麦克风接收到的物理误差信号共同作用更新物理次级通路滤波器
Figure FSA0000285341430000014
最终达到稳定状态,得到物理次级通路滤波器
Figure FSA0000285341430000015
步骤四:训练辅助滤波器H,目标降噪位置不放置任何麦克风,物理位置放置物理麦克风,将步骤二训练得到的噪声控制滤波器W作为固定系数滤波器W0,物理麦克风接收到的物理误差信号为参考信号通过固定系数滤波器W0和物理次级通路滤波器
Figure FSA0000285341430000016
的输出和初级噪声的叠加,辅助滤波器产生的残余信号为物理麦克风接收到的物理误差信号与参考信号通过辅助滤波器H的输出的叠加,参考信号与辅助滤波器产生的残余信号共同作用更新辅助滤波器H,最终达到稳定状态,得到辅助滤波器H;
步骤五:进入控制阶段,目标降噪位置不放置任何麦克风,物理位置放置物理麦克风,物理麦克风接收到的物理误差信号为参考信号通过噪声控制滤波器W和物理次级通路滤波器
Figure FSA0000285341430000017
的输出和初级噪声的叠加,将步骤四训练得到的辅助滤波器H作为固定系数滤波器H0,辅助滤波器产生的残余信号为物理麦克风接收到的物理误差信号与参考信号通过固定系数滤波器H0的输出的叠加,参考信号通过物理次级通路滤波器
Figure FSA0000285341430000018
的输出与辅助滤波器产生的残余信号共同作用更新噪声控制滤波器W,最终达到稳定状态,实现目标降噪位置的有源降噪。
2.根据权利要求1所述的一种分步虚拟传感降噪方法,其特征在于,步骤二中,噪声控制滤波器W更新的过程为:
通过公式
Figure FSA0000285341430000021
对权系数矢量进行更新,其中w(n)表示噪声控制滤波器更新前的权系数矢量,w(n+1)表示更新后的权系数矢量,x(n)表示参考信号,
Figure FSA0000285341430000022
表示虚拟次级通路滤波器,ev(n)表示虚拟误差信号,μw1表示噪声控制滤波器更新的收敛因子。
3.根据权利要求1所述的一种分步虚拟传感降噪方法,其特征在于,步骤四中,辅助滤波器H更新的过程为:
通过公式h(n+1)=h(n)-μheh(n)x(n)对权系数矢量进行更新,其中h(n)表示辅助滤波器更新前的权系数矢量,h(n+1)表示更新后的权系数矢量,x(n)表示参考信号,eh(n)表示辅助滤波器产生的残余信号,μh表示辅助滤波器更新的收敛因子。
4.根据权利要求1所述的一种分步虚拟传感降噪方法,其特征在于,步骤五中,噪声控制滤波器W更新的过程为:
通过公式w(n+1)=w(n)-μw2eh(n)xs(n)对权系数矢量进行更新,其中w(n)表示噪声控制滤波器更新前的权系数矢量,w(n+1)表示更新后的权系数矢量,xs(n)表示参考信号通过物理次级通路滤波器的输出,eh(n)表示辅助滤波器产生的残余信号,μw2表示噪声控制滤波器更新的收敛因子。
5.根据权利要求2至4所述的一种分步虚拟传感降噪方法,其特征在于,收敛因子设置为:
Figure FSA0000285341430000023
Figure FSA0000285341430000024
Figure FSA0000285341430000025
其中,μw1表示噪声控制滤波器在步骤二中更新的收敛因子,μh表示辅助滤波器在步骤四中更新的收敛因子,μw2表示噪声控制滤波器在步骤五中更新的收敛因子,x表示参考信号。
CN202211186487.8A 2022-09-27 2022-09-27 一种分步虚拟传感降噪方法 Pending CN115565515A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211186487.8A CN115565515A (zh) 2022-09-27 2022-09-27 一种分步虚拟传感降噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211186487.8A CN115565515A (zh) 2022-09-27 2022-09-27 一种分步虚拟传感降噪方法

Publications (1)

Publication Number Publication Date
CN115565515A true CN115565515A (zh) 2023-01-03

Family

ID=84742881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211186487.8A Pending CN115565515A (zh) 2022-09-27 2022-09-27 一种分步虚拟传感降噪方法

Country Status (1)

Country Link
CN (1) CN115565515A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116959472A (zh) * 2023-09-20 2023-10-27 武汉海微科技有限公司 主动降噪方法、装置、设备及存储介质

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116959472A (zh) * 2023-09-20 2023-10-27 武汉海微科技有限公司 主动降噪方法、装置、设备及存储介质

Similar Documents

Publication Publication Date Title
Kuo et al. Active noise control: Open problems and challenges
Chang et al. Active noise control in headsets by using a low-cost microcontroller
Shi et al. Block coordinate descent based algorithm for computational complexity reduction in multichannel active noise control system
CN111627415A (zh) 一种基于自适应MFxLMS算法的主动降噪装置及FPGA实现
Shen et al. Alternative switching hybrid ANC
Zhang et al. Deep MCANC: A deep learning approach to multi-channel active noise control
Mostafavi et al. Deep learning-based active noise control on construction sites
CN101345783B (zh) 一种噪声消除的***和方法
Ma et al. Active control of outgoing broadband noise fields in rooms
Shi et al. Computation-efficient solution for fully-connected active noise control window: Analysis and implementation of multichannel adjoint least mean square algorithm
CN105913836B (zh) 一种基于dsp的定点化实时降噪方法
WO2023124629A1 (zh) 一种车辆的主动降噪方法及设备、存储介质
CN115565515A (zh) 一种分步虚拟传感降噪方法
Ren et al. An improved artificial bee colony algorithm for model-free active noise control: algorithm and implementation
Hu et al. Invisible speakers in home windows for simultaneous auxiliary audio playback and active noise cancellation
Tang et al. Stability guaranteed active noise control: Algorithms and applications
Shen et al. The principle underlying the wireless reference microphone enhancing noise reduction performance in multi-channel active noise control windows
Xie et al. A feedforward and feedback composite active noise reduction headset based on inverse filter frequency equalization and its DSP system implementation
CN113096629B (zh) 用于单通道反馈有源噪声控制***的相对路径虚拟传感方法
Luo et al. EC-ANC: Edge Case-Enhanced Active Noise Cancellation for True Wireless Stereo Earbuds
Kuo et al. Active noise control systems with optimized secondary path
Yuan et al. Design of Active Noise Control System Based on FPGA
Kalaivani et al. Analysis of Adaptive Filter Algorithms in Real Time Signals
Huynh et al. A novel adaptive neural controller for narrowband active noise control systems
Zhaohan A Low-cost and Portable Active Noise Control Unit

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