CN111583896A - 一种多通道有源降噪头靠的降噪方法 - Google Patents

一种多通道有源降噪头靠的降噪方法 Download PDF

Info

Publication number
CN111583896A
CN111583896A CN202010401900.2A CN202010401900A CN111583896A CN 111583896 A CN111583896 A CN 111583896A CN 202010401900 A CN202010401900 A CN 202010401900A CN 111583896 A CN111583896 A CN 111583896A
Authority
CN
China
Prior art keywords
microphone
virtual
physical
signal
sound source
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
Application number
CN202010401900.2A
Other languages
English (en)
Other versions
CN111583896B (zh
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.)
Suzhou Silent Tai Technology Co ltd
Original Assignee
Suzhou Silent Tai Technology Co ltd
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 Suzhou Silent Tai Technology Co ltd filed Critical Suzhou Silent Tai Technology Co ltd
Priority to CN202010401900.2A priority Critical patent/CN111583896B/zh
Publication of CN111583896A publication Critical patent/CN111583896A/zh
Application granted granted Critical
Publication of CN111583896B publication Critical patent/CN111583896B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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
    • 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/17823—Reference signals, e.g. ambient acoustic environment
    • 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/1787—General system configurations
    • G10K11/17879—General system configurations using both a reference signal and an error signal
    • 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
    • G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30—Means
    • G10K2210/321—Physical
    • G10K2210/3221—Headrests, seats or the like, for personal ANC systems
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开了一种多通道有源降噪头靠的降噪方法，通过采集计算模块得到耳廓空间位置信息，并选定耳廓上一点坐标为虚拟误差传声器的空间坐标，通过已知空间坐标的物理传声器和次级声源，计算得到从次级声源和物理传声器到虚拟误差传声器的时延值，使用离线建模的方式获得次级声源至物理传声器的物理次级路径传递函数矩阵Gpu；自适应控制器根据参考信号、物理误差信号、传递函数矩阵Gvu、M和Gpu以及次级声源的强度u(n)计算得到虚拟误差信号ev，自适应控制器根据参考信号与虚拟误差信号进行计算，产生消声信号，并传递给次级声源；次级声源接收自适应控制器传来的消声信号并发出声音，对原有噪声进行抵消降噪，使虚拟误差传声器处的声压幅值平方最小。

Description

一种多通道有源降噪头靠的降噪方法
技术领域
本发明涉及噪声控制领域，特别涉及一种多通道有源降噪头靠的降噪方法。
背景技术
噪声污染是一种对人体影响较大的环境污染，当人们长期在高噪声环境中时，工作效率会降低，人的听力和身心健康也会受到影响。尤其是对于汽车、高铁、轮船、和飞机等运载交通工具，舱室区域的噪声严重影响司乘人员乘坐时的舒适性，舱室噪声中的低频噪声成分很难控制，被动噪声控制对材料的尺寸、体积有较高要求，在舱室中难以满足。有源噪声控制(Active noise control,ANC)技术具有体积小、低频噪声控制效果好的特点。有源噪声控制技术基于惠更斯波形叠加原理，因此对舱室噪声进行全局控制难以满足，一般采用局部有源噪声控制的方法降低乘客头部区域的噪声。在无遮挡的条件要求下，有源头靠、头枕(Active headrest)是局部空间有源噪声控制的成功应用之一。
专利CN201620075679.5和CN201810463634.9分别公开了一种有源降噪头枕，将噪声控制系统一体化集成在头枕中，控制装置根据环境噪声和校正噪声控制扬声器产生控制声波，降低使用者头部活动区域的低频噪声，提升声学舒适性。但在一般的有源噪声控制系统中，误差传声器应尽量靠近人耳布放，次级源根据应用场景布放于人耳周围，这样才能最终在误差传声器附近产生有源静音区。通常情况下，静音区的大小与噪声波长相关，单通道有源降噪系统产生的静区直径一般小于1/10波长，为了保证降噪效果，误差传声器应尽可能靠近人耳，这就造成了误差传声器与人头移动的冲突。在实际使用过程中，司乘人员不可能一直紧贴着头靠，当人头移动时，有源降噪头靠的降噪效果就大打折扣，甚至失去效果。
为解决人头移动的影响，专利201910998134.X引入虚拟传声器技术，通过目标区域外的物理传声器拾取声信号对人耳附近虚拟传声器处(原误差传声器位置)的声压进行估计和控制，进而转移有源静区至人耳附近来提高有源降噪系统在人头移动时的鲁棒性。但该专利仅仅提出了当人头前后左右移动和转动时的虚拟传声器与物理传声器之前的传递路径计算方法，采用移动路径估计的方法计算出使人头位于不同位置时系统的噪声残余最大值最小化的最优虚拟次级路径传递函数估计值，当复杂运动时计算量急剧增大。该方法还将人头建模为一刚性球，双耳在刚性球表面，位于通过刚性球中心的直径上，与实际情况并不相符，无法取得真实的良好控制效果。
司乘人员在乘坐过程中，头部的运动是由上下、左右、前后移动和左右转动等多种方式组合而成的复杂运动，虚拟传声器与物理传声器之间的映射关系是动态变化的。因此需要提供一种能够直接获取虚拟传声器传递函数的方法，提供给有源降噪头靠系统以解决以上问题。
发明内容
本发明的目的是提供一种多通道有源降噪头靠的降噪方法，在人耳移动过程中，能够精确的得到虚拟误差传声器的位置，经过简单计算即可得到有源降噪系统所需的虚拟次级路径传递函数矩阵Gvu和初级声场传递函数矩阵M，系统响应快，降噪效果好。
为达到上述目的，本发明采用的技术方案是：
一种多通道有源降噪头靠的降噪方法，包括以下步骤：
(1)通过装载在头靠上的采集计算模块得到耳廓空间位置信息；
(2)选定耳廓上某一特征点作为抓取点，将所述抓取点的空间坐标定义为虚拟误差传声器的空间坐标；
(3)通过装载在所述头靠上的已知空间坐标的物理传声器和次级声源，计算得到从所述次级声源和所述物理传声器到所述虚拟误差传声器的时延值，所述物理传声器包括物理参考传声器和物理误差传声器；所述次级声源至所述虚拟误差传声器的虚拟次级路径传递函数矩阵为Gvu，所述物理误差传声器与所述虚拟误差传声器之间的初级声场传递函数矩阵为M，由于所述次级声源与所述物理误差传声器位置固定，使用离线建模的方式获得固定的所述次级声源至固定的所述物理传声器的物理次级路径传递函数矩阵Gpu；
(4)所述物理参考传声器采集噪声产生参考信号，所述物理误差传声器采集噪声产生物理误差信号ep；将所述参考信号和所述物理误差信号传递给自适应控制器，所述自适应控制器根据传递函数矩阵Gvu、M和Gpu以及所述次级声源的强度u(n)计算得到虚拟误差信号ev，所述自适应控制器根据所述参考信号与所述虚拟误差信号进行计算，产生消声信号，并传递给所述次级声源；所述次级声源接收所述自适应控制器传来的所述消声信号并发出声音，对原有噪声进行抵消降噪，使所述虚拟误差传声器处的声压幅值平方最小。
优选地，在步骤(1)中，所述采集计算模块包括RGB-D相机、卷积神经网络结构模型；所述卷积神经网络结构模型用于识别人耳耳廓；所述RGB-D相机包括彩色传声器和深度传声器，根据两者之间的外参构建彩色图像数据与深度数据之间的映射关系，以得到以像素坐标[u,v]形式存在的图像，利用内参进行坐标系变化以得到空间坐标P(x,y,z)。
更优选地，像素平面中P′的像素坐标为[u,v],根据所述RGB-D相机的内参K 和像素P′[u,v]对应的深度Z由公式(1)得出空间坐标P(x,y,z)；
Figure BDA0002489799820000031
其中，K矩阵为RGB-D相机的内参数矩阵，P为P′的三维空间坐标。
优选地，在步骤(2)中，所述抓取点为耳廓的三角窝角点或耳甲腔边缘的角点。
优选地，所述头靠上对应每一侧人耳均设有所述采集计算模块、所述物理传声器和所述次级声源。
优选地，所述多通道有源降噪头靠的降噪方法采用多通道自适应前馈控制系统，使用可移动的虚拟误差传声器作为误差传声器，建立模型如下：
设控制系统有L个次级声源、Mp个物理传声器、Mv个移动虚拟误差传声器；移动虚拟误差传声器的时变位置包含在大小为3×Mv的矩阵xv(n)中，定义为：
Figure BDA0002489799820000032
其中每个移动虚拟位置
Figure BDA0002489799820000033
由相对于一个参考系的三个空间坐标定义：
Figure BDA0002489799820000034
计算模型需要在每个时间步长都知道Mv个移动虚拟误差传声器的空间坐标 xv(n)，通过步骤(2)得到移动虚拟误差传声器的空间坐标；
在移动的虚拟位置xv(n)上计算虚拟误差信号的估计值
Figure BDA0002489799820000101
，
在移动虚拟传感算法中，首先需要获得
Figure BDA0002489799820000086
个空间坐标固定的虚拟误差传声器
Figure BDA0002489799820000087
的虚拟误差信号的估计值
Figure BDA0002489799820000088
设移动虚拟误差传声器的空间虚拟位置 xv(n)被限制在一个有限的三维区域中，此时空间位置固定的虚拟误差传声器
Figure BDA0002489799820000087
也就位于所述三维区域中，
Figure BDA0002489799820000086
个空间坐标固定的虚拟误差传声器位置的向量由式(4)给出；
Figure BDA0002489799820000048
则每个空间坐标固定的虚拟位置
Figure BDA0002489799820000049
由相对于一个参考系的三个空间坐标定义：
Figure BDA00024897998200000410
空间位置固定的虚拟误差传声器
Figure BDA00024897998200000411
的误差信号估计值
Figure BDA00024897998200000412
由物理误差传声器信号与次级声源强度传递函数得到：
Figure BDA00024897998200000413
其中，ep(n)为物理误差传声器的物理误差信号，
Figure BDA00024897998200000414
是阶数为Mp×L的次级声源到物理误差传声器的传递函数矩阵，
Figure BDA00024897998200000415
是阶数为
Figure BDA00024897998200000416
的次级声源到虚拟误差传声器的传递函数矩阵，u(n)为次级声源的强度矢量；
将两个次级通道传递函数矩阵建模成两个FIR或两个IIR滤波器；在系统辨识的初级阶段，将阶数为
Figure BDA00024897998200000417
的从空间坐标固定的物理传声器到虚拟误差传声器之间的次级通道传递函数矩阵
Figure BDA0002489799820000051
建模成一个FIR或IIR滤波器；Mp个物理误差传声器的初级噪声信号为：
Figure BDA0002489799820000052
空间位置固定的虚拟误差传声器
Figure BDA0002489799820000053
的初级噪声估计值为：
Figure BDA0002489799820000054
由此空间位置固定的虚拟误差传声器
Figure BDA00024897998200000510
的总的虚拟误差信号的估计值
Figure BDA0002489799820000055
为：
Figure BDA0002489799820000056
根据所述移动虚拟传感算法，通过对空间坐标固定的虚拟误差传声器
Figure BDA00024897998200000510
的虚拟误差信号
Figure BDA00024897998200000916
进行空间插值，可以获得在移动的虚拟位置上xv(n)的虚拟误差信号的估计值
Figure BDA0002489799820000101
。
更优选地，所述自适应控制器根据所述参考信号与所述虚拟误差信号，采用自适应有源控制算法进行计算，产生所述消声信号。
更进一步优选地，所述自适应有源控制算法为多通道Fx-LMS算法。
由于上述技术方案的运用，本发明与现有技术相比具有下列优点：本发明一种多通道有源降噪头靠的降噪方法，具有以下优点：
(1)采用多通道有源降噪系统，经合理设计可扩大静音区的有效范围，提升降噪性能；
(2)在人耳移动过程中，通过采集计算模块能够直接提取耳廓空间坐标信息，以精确的得到虚拟误差传声器的位置，用硬件的测量精度来保证降噪控制效果，提升系统鲁棒性；
(3)使用直接获取的虚拟误差传声器空间坐标，经过简单计算即可得到有源降噪系统所需的虚拟次级路径传递函数矩阵Gvu和初级声场传递函数矩阵M，系统响应快，降噪效果好。
附图说明
附图1为多通道有源降噪头靠的结构示意图一；
附图2为多通道有源降噪头靠的结构示意图二；
附图3为声压信号的时间历程图；
附图4为功率谱密度图。
附图5为移动虚拟传感算法框图。
其中：1、RGB-D相机；2、次级声源；3、物理参考传声器；4、物理误差传声器；5、虚拟误差传声器；6、头靠；7、数据处理终端。
具体实施方式
下面结合附图来对本发明的技术方案作进一步的阐述。
参见图1-2所示，为多通道有源降噪头靠，包括双耳廓空间位置获取系统和多通道有源噪声控制系统。
双耳廓空间位置获取系统为用于得到耳廓空间位置信息的采集计算模块。在本实施例中，该采集计算模块包括RGB-D相机1、卷积神经网络结构模型、数据处理器。该卷积神经网络结构模型用于识别人耳耳廓(为现有的成熟技术)。上述RGB-D相机1包括彩色传声器和深度传声器，该组合相机既能够获取RGB 彩色图像信息，同时也能获取深度图像信息。RGB-D相机1可以主动测量每个像素对应的距离，并能将RGB彩色图像信息和深度数据进行耦合。RGB-D相机1 模型根据RGB彩色摄像头和深度传声器两者之间的外参构建RGB彩色图像数据与深度数据两者之间的映射关系，相机最终得到的图像以像素坐标的形式存在，像素平面中像素坐标为[u,v]。RGB-D深度相机利用相机内参进行坐标系变换可以直接获得测量点的空间位置信息，获取空间坐标P(x,y,z)。
像素平面中P′的像素坐标为[u,v],根据RGB-D相机1的内参K和像素P′[u,v] 对应的深度Z由公式(1)得出空间坐标P(x,y,z)；
Figure BDA0002489799820000061
其中，K矩阵为RGB-D相机1的内参数矩阵，P为P′的三维空间坐标。
在本实施例中，头靠6上装有两只RGB-D相机1，左右对称的设于头靠6 两侧，用于分别拍摄两只人耳。
多通道有源噪声控制系统包括次级声源2(即次级扬声器声源)、物理参考传声器3、物理误差传声器4、自适应控制器。其中，次级扬声器声源和物理传声器数目可以根据控制系统尺寸大小和降噪要求增减。在本实施例中，头靠6 上装有四只次级扬声器声源、两只物理参考传声器3、两只物理误差传声器4，左右对称的设于头靠6两侧，用于分别对两只人耳作用。
其中，多通道有源噪声控制系统中的自适应控制器与双耳廓空间位置获取系统中的数据处理器共称为数据处理终端7，安装在头靠6中。
一种多通道有源降噪头靠的降噪方法，通过上述多通道有源降噪头靠实现，包括以下步骤：
(1)通过装载在头靠6上的采集计算模块得到耳廓空间位置信息；
(2)选定耳廓上某一特征点作为抓取点，将抓取点的空间坐标定义为虚拟误差传声器5的空间坐标；
具体的，耳廓空间位置获取系统在工作时，首先对耳廓部位四周进行拍照识别并进行特点提取，保存拍摄图片中检测出的各特征点的像素坐标，选定一个特征点作为抓取点，即拟测量的物体(人耳)上某个具体空间坐标的点，这个点可以是物体上任意处的点。在本实施例中，为方便计算，系统一般选择耳廓的三角窝角点或耳甲腔边缘的角点作为抓取点，借助深度相机和公式测量保存各特征点和抓取点的空间坐标。
(3)通过装载在头靠6上的已知空间坐标的物理传声器和次级声源2，计算得到从次级声源2和物理传声器到虚拟误差传声器5的时延值，物理传声器包括物理参考传声器3和物理误差传声器4；次级声源2至虚拟误差传声器5 的虚拟次级路径传递函数矩阵为Gvu，物理误差传声器4与虚拟误差传声器5 之间的初级声场传递函数矩阵为M，由于次级声源2与物理误差传声器4位置固定，使用离线建模的方式获得固定的次级声源2至固定的物理传声器的物理次级路径传递函数矩阵Gpu；
(4)物理参考传声器3采集噪声产生参考信号，物理误差传声器4采集噪声产生物理误差信号ep；将参考信号和物理误差信号传递给自适应控制器，自适应控制器根据传递函数矩阵Gvu、M和Gpu以及次级声源2的强度u(n)计算得到虚拟误差信号ev，自适应控制器根据参考信号与虚拟误差信号进行计算，产生消声信号，并传递给次级声源2；次级声源2接收自适应控制器传来的消声信号并发出声音，对原有噪声进行抵消降噪，使虚拟误差传声器5处的声压幅值平方最小。由于虚拟误差传声器5的位置在耳廓上，因此可以保证耳廓始终处于静音区，不会因为人头不规则移动等因素造成降噪失效。多通道有源降噪系统根据参考信号与虚拟误差信号产生消声信号的算法使用自适应有源控制算法，如多通道Fx-LMS算法。该算法在教科书中已有介绍，此处不再赘述。
具体的，多通道有源降噪头靠的降噪方法采用多通道自适应前馈控制系统，使用可移动的虚拟误差传声器5作为误差传声器，建立模型如下：
设控制系统有L个次级声源2、Mp个物理传声器、Mv个移动虚拟误差传声器5；移动虚拟误差传声器5的时变位置包含在大小为3×Mv的矩阵xv(n)中，定义为：
Figure BDA0002489799820000081
其中每个移动虚拟位置
Figure BDA0002489799820000082
由相对于一个参考系的三个空间坐标定义：
Figure BDA0002489799820000083
计算模型需要在每个时间步长都知道Mv个移动虚拟误差传声器5的空间坐标xv(n)，通过步骤(2)实时测量特征点空间坐标来获取移动虚拟误差传声器5 的空间坐标；
参见图5所示，在移动的虚拟位置xv(n)上计算虚拟误差信号的估计值
Figure BDA0002489799820000101
，
在移动虚拟传感算法中，首先需要获得
Figure BDA0002489799820000086
个空间坐标固定的虚拟误差传声器
Figure BDA0002489799820000087
的虚拟误差信号的估计值
Figure BDA0002489799820000088
设移动虚拟误差传声器5的空间虚拟位置 xv(n)被限制在一个有限的三维区域中，此时空间位置固定的虚拟误差传声器
Figure BDA0002489799820000089
也就位于该三维区域中，
Figure BDA00024897998200000810
个空间坐标固定的虚拟误差传声器5位置的向量由式(4)给出；
Figure BDA0002489799820000091
则每个空间坐标固定的虚拟位置
Figure BDA0002489799820000092
由相对于一个参考系的三个空间坐标定义：
Figure BDA0002489799820000093
空间位置固定的虚拟误差传声器
Figure BDA0002489799820000094
的误差信号估计值
Figure BDA0002489799820000095
由物理误差传声器4信号与次级声源2强度传递函数得到：
Figure BDA0002489799820000096
其中，ep(n)为物理误差传声器4的物理误差信号，
Figure BDA0002489799820000097
是阶数为Mp×L的次级声源2到物理误差传声器4的传递函数矩阵，
Figure BDA0002489799820000098
是阶数为
Figure BDA0002489799820000099
的次级声源2 到虚拟误差传声器5的传递函数矩阵，u(n)为次级声源2的强度矢量；
将两个次级通道传递函数矩阵建模成两个FIR或两个IIR滤波器；在系统辨识的初级阶段，将阶数为
Figure BDA00024897998200000910
的从空间坐标固定的物理传声器到虚拟误差传声器5之间的次级通道传递函数矩阵
Figure BDA00024897998200000911
建模成一个FIR或IIR滤波器；Mp个物理误差传声器4的初级噪声信号为：
Figure BDA00024897998200000912
空间位置固定的虚拟误差传声器
Figure BDA00024897998200000913
的初级噪声估计值为：
Figure BDA00024897998200000914
由此空间位置固定的虚拟误差传声器
Figure BDA00024897998200000915
的总的虚拟误差信号的估计值
Figure BDA00024897998200000916
为：
Figure BDA00024897998200000917
根据移动虚拟传感算法，通过对空间坐标固定的虚拟误差传声器
Figure BDA00024897998200000915
的虚拟误差信号
Figure BDA00024897998200000916
进行空间插值，可以获得在移动的虚拟位置上xv(n)的虚拟误差信号的估计值
Figure BDA0002489799820000101
。
移动的虚拟位置上xv(n)可以通过双耳廓空间位置获取系统实时获取，因此通过简单的插值计算即可得到虚拟误差信号的估计值
Figure BDA0002489799820000101
来进行反馈控制，极大地减少了系统运算量，提高运算效率。
使用人工头模拟司乘人员头部进行有源降噪头靠6的效能测试。
图3给出了在人工头的左耳中测得的声压信号的时间历程，横坐标为时间(S)，纵坐标为声压(Pa)。
其中前30秒采集到的是控制系统不工作时测量的噪声信号即原始环境噪声；从30秒到60秒采集到的信号是使用物理误差传声器4采集的信号作为控制系统误差信号时，对人耳区域噪声进行有源控制后的结果，可以看出一点控制效果，但不明显；从60秒以后采集到的信号是启用了头部耳廓跟踪后，使用移动虚拟传感技术进行有源控制后的结果，可以看出降噪效果极为明显，同时降噪效果不受人头移动的影响，控制效果非常稳定。在60秒后能识别正确的头部位置。
图4给出了控制前、固定物理误差传声器4控制和移动虚拟传感控制三种情况下的信号功率谱密度。
上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围内。

Claims (8)

1.一种多通道有源降噪头靠的降噪方法，其特征在于：包括以下步骤：
(1)通过装载在头靠上的采集计算模块得到耳廓空间位置信息；
(2)选定耳廓上某一特征点作为抓取点，将所述抓取点的空间坐标定义为虚拟误差传声器的空间坐标；
(3)通过装载在所述头靠上的已知空间坐标的物理传声器和次级声源，计算得到从所述次级声源和所述物理传声器到所述虚拟误差传声器的时延值，所述物理传声器包括物理参考传声器和物理误差传声器；所述次级声源至所述虚拟误差传声器的虚拟次级路径传递函数矩阵为Gvu，所述物理误差传声器与所述虚拟误差传声器之间的初级声场传递函数矩阵为M，由于所述次级声源与所述物理误差传声器位置固定，使用离线建模的方式获得固定的所述次级声源至固定的所述物理传声器的物理次级路径传递函数矩阵Gpu；
(4)所述物理参考传声器采集噪声产生参考信号，所述物理误差传声器采集噪声产生物理误差信号ep；将所述参考信号和所述物理误差信号传递给自适应控制器，所述自适应控制器根据传递函数矩阵Gvu、M和Gpu以及所述次级声源的强度u(n)计算得到虚拟误差信号ev，所述自适应控制器根据所述参考信号与所述虚拟误差信号进行计算，产生消声信号，并传递给所述次级声源；所述次级声源接收所述自适应控制器传来的所述消声信号并发出声音，对原有噪声进行抵消降噪，使所述虚拟误差传声器处的声压幅值平方最小。
2.根据权利要求1所述的一种多通道有源降噪头靠的降噪方法，其特征在于：在步骤(1)中，所述采集计算模块包括RGB-D相机、卷积神经网络结构模型；所述卷积神经网络结构模型用于识别人耳耳廓；所述RGB-D相机包括彩色传声器和深度传声器，根据两者之间的外参构建彩色图像数据与深度数据之间的映射关系，以得到以像素坐标[u,v]形式存在的图像，利用内参进行坐标系变化以得到空间坐标P(x,y,z)。
3.根据权利要求2所述的一种多通道有源降噪头靠的降噪方法，其特征在于：像素平面中P′的像素坐标为[u,v],根据所述RGB-D相机的内参K和像素P′[u,v]对应的深度Z由公式(1)得出空间坐标P(x,y,z)；
Figure FDA0002489799810000011
其中，K矩阵为RGB-D相机的内参数矩阵，P为P′的三维空间坐标。
4.根据权利要求1所述的一种多通道有源降噪头靠的降噪方法，其特征在于：在步骤(2)中，所述抓取点为耳廓的三角窝角点或耳甲腔边缘的角点。
5.根据权利要求1所述的一种多通道有源降噪头靠的降噪方法，其特征在于：所述头靠上对应每一侧人耳均设有所述采集计算模块、所述物理传声器和所述次级声源。
6.根据权利要求1所述的一种多通道有源降噪头靠的降噪方法，其特征在于：所述多通道有源降噪头靠的降噪方法采用多通道自适应前馈控制系统，使用可移动的虚拟误差传声器作为误差传声器，建立模型如下：
设控制系统有L个次级声源、Mp个物理传声器、Mv个移动虚拟误差传声器；移动虚拟误差传声器的时变位置包含在大小为3×Mv的矩阵xv(n)中，定义为：
Figure FDA0002489799810000021
其中每个移动虚拟位置
Figure FDA0002489799810000024
由相对于一个参考系的三个空间坐标定义：
Figure FDA0002489799810000022
计算模型需要在每个时间步长都知道Mv个移动虚拟误差传声器的空间坐标xv(n)，通过步骤(2)得到移动虚拟误差传声器的空间坐标；
在移动的虚拟位置xv(n)上计算虚拟误差信号的估计值
Figure FDA0002489799810000025
在移动虚拟传感算法中，首先需要获得
Figure FDA0002489799810000026
个空间坐标固定的虚拟误差传声器
Figure FDA0002489799810000038
的虚拟误差信号的估计值
Figure FDA0002489799810000037
设移动虚拟误差传声器的空间虚拟位置xv(n)被限制在一个有限的三维区域中，此时空间位置固定的虚拟误差传声器
Figure FDA0002489799810000039
也就位于所述三维区域中，
Figure FDA00024897998100000310
个空间坐标固定的虚拟误差传声器位置的向量由式(4)给出；
Figure FDA0002489799810000031
则每个空间坐标固定的虚拟位置
Figure FDA00024897998100000311
由相对于一个参考系的三个空间坐标定义：
Figure FDA0002489799810000032
空间位置固定的虚拟误差传声器
Figure FDA00024897998100000312
的误差信号估计值
Figure FDA00024897998100000313
由物理误差传声器信号与次级声源强度传递函数得到：
Figure FDA0002489799810000033
其中，ep(n)为物理误差传声器的物理误差信号，
Figure FDA00024897998100000314
是阶数为Mp×L的次级声源到物理误差传声器的传递函数矩阵，
Figure FDA00024897998100000315
是阶数为
Figure FDA00024897998100000316
的次级声源到虚拟误差传声器的传递函数矩阵，u(n)为次级声源的强度矢量；
将两个次级通道传递函数矩阵建模成两个FIR或两个IIR滤波器；在系统辨识的初级阶段，将阶数为
Figure FDA00024897998100000317
的从空间坐标固定的物理传声器到虚拟误差传声器之间的次级通道传递函数矩阵
Figure FDA00024897998100000318
建模成一个FIR或IIR滤波器；Mp个物理误差传声器的初级噪声信号为：
Figure FDA0002489799810000034
空间位置固定的虚拟误差传声器
Figure FDA00024897998100000319
的初级噪声估计值为：
Figure FDA0002489799810000035
由此空间位置固定的虚拟误差传声器
Figure FDA00024897998100000320
的总的虚拟误差信号的估计值
Figure FDA00024897998100000321
为：
Figure FDA0002489799810000036
根据所述移动虚拟传感算法，通过对空间坐标固定的虚拟误差传声器
Figure FDA0002489799810000041
的虚拟误差信号
Figure FDA0002489799810000042
进行空间插值，可以获得在移动的虚拟位置上xv(n)的虚拟误差信号的估计值
Figure FDA0002489799810000043
7.根据权利要求6所述的一种多通道有源降噪头靠的降噪方法，其特征在于：所述自适应控制器根据所述参考信号与所述虚拟误差信号，采用自适应有源控制算法进行计算，产生所述消声信号。
8.根据权利要求7所述的一种多通道有源降噪头靠的降噪方法，其特征在于：所述自适应有源控制算法为多通道Fx-LMS算法。
CN202010401900.2A 2020-05-13 2020-05-13 一种多通道有源降噪头靠的降噪方法 Active CN111583896B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010401900.2A CN111583896B (zh) 2020-05-13 2020-05-13 一种多通道有源降噪头靠的降噪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010401900.2A CN111583896B (zh) 2020-05-13 2020-05-13 一种多通道有源降噪头靠的降噪方法

Publications (2)

Publication Number Publication Date
CN111583896A true CN111583896A (zh) 2020-08-25
CN111583896B CN111583896B (zh) 2023-09-08

Family

ID=72112685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010401900.2A Active CN111583896B (zh) 2020-05-13 2020-05-13 一种多通道有源降噪头靠的降噪方法

Country Status (1)

Country Link
CN (1) CN111583896B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112102806A (zh) * 2020-09-06 2020-12-18 西安艾科特声学科技有限公司 列车驾驶室基于虚拟传感技术的主动噪声控制系统及方法
CN112331173A (zh) * 2020-10-26 2021-02-05 Tcl通力电子（惠州）有限公司 车内降噪方法、控制器和车载枕头和计算机可读存储介质
CN112669804A (zh) * 2020-12-11 2021-04-16 西北工业大学 有源噪声控制系统降噪效果估计方法
CN112669806A (zh) * 2020-12-24 2021-04-16 浙江大学 一种集成灶主动降噪装置
CN112771606A (zh) * 2020-12-31 2021-05-07 华为技术有限公司 一种汽车空调噪声控制系统、方法及相关车载设备
CN113080928A (zh) * 2021-03-30 2021-07-09 上海联影医疗科技股份有限公司 主动降噪系统和具有其的医疗设备
CN113096629A (zh) * 2021-03-03 2021-07-09 电子科技大学 用于单通道反馈有源噪声控制系统的相对路径虚拟传感方法
CN113436599A (zh) * 2021-05-31 2021-09-24 新东方教育科技集团有限公司 信号处理方法、装置、存储介质及电子设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017000814A1 (zh) * 2015-06-30 2017-01-05 芋头科技（杭州）有限公司 一种有源噪声控制系统
CN108428445A (zh) * 2018-03-15 2018-08-21 中国科学院声学研究所 一种无误差传声器的自适应主动降噪方法
CN108806664A (zh) * 2018-05-03 2018-11-13 清华大学苏州汽车研究院（相城） 一种汽车车内噪声控制方法
CN110719550A (zh) * 2019-10-21 2020-01-21 南京南大电子智慧型服务机器人研究院有限公司 一种双通道有源降噪头靠的虚拟传声器优化设计方法
CN111031442A (zh) * 2019-12-31 2020-04-17 湖南景程电子科技有限公司 一种基于dsp虚拟传感的耳机主动降噪自适应系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017000814A1 (zh) * 2015-06-30 2017-01-05 芋头科技（杭州）有限公司 一种有源噪声控制系统
CN108428445A (zh) * 2018-03-15 2018-08-21 中国科学院声学研究所 一种无误差传声器的自适应主动降噪方法
CN108806664A (zh) * 2018-05-03 2018-11-13 清华大学苏州汽车研究院（相城） 一种汽车车内噪声控制方法
CN110719550A (zh) * 2019-10-21 2020-01-21 南京南大电子智慧型服务机器人研究院有限公司 一种双通道有源降噪头靠的虚拟传声器优化设计方法
CN111031442A (zh) * 2019-12-31 2020-04-17 湖南景程电子科技有限公司 一种基于dsp虚拟传感的耳机主动降噪自适应系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘波: "单通道有源噪声控制技术研究" *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112102806B (zh) * 2020-09-06 2024-04-26 西安艾科特声学科技有限公司 列车驾驶室基于虚拟传感技术的主动噪声控制系统及方法
CN112102806A (zh) * 2020-09-06 2020-12-18 西安艾科特声学科技有限公司 列车驾驶室基于虚拟传感技术的主动噪声控制系统及方法
CN112331173A (zh) * 2020-10-26 2021-02-05 Tcl通力电子（惠州）有限公司 车内降噪方法、控制器和车载枕头和计算机可读存储介质
CN112331173B (zh) * 2020-10-26 2024-02-23 通力科技股份有限公司 车内降噪方法、控制器和车载枕头和计算机可读存储介质
CN112669804B (zh) * 2020-12-11 2024-02-09 西北工业大学 有源噪声控制系统降噪效果估计方法
CN112669804A (zh) * 2020-12-11 2021-04-16 西北工业大学 有源噪声控制系统降噪效果估计方法
CN112669806A (zh) * 2020-12-24 2021-04-16 浙江大学 一种集成灶主动降噪装置
CN112771606A (zh) * 2020-12-31 2021-05-07 华为技术有限公司 一种汽车空调噪声控制系统、方法及相关车载设备
WO2022141562A1 (zh) * 2020-12-31 2022-07-07 华为技术有限公司 一种汽车空调噪声控制系统、方法及相关车载设备
CN113096629A (zh) * 2021-03-03 2021-07-09 电子科技大学 用于单通道反馈有源噪声控制系统的相对路径虚拟传感方法
CN113096629B (zh) * 2021-03-03 2022-11-04 电子科技大学 用于单通道反馈有源噪声控制系统的相对路径虚拟传感方法
CN113080928A (zh) * 2021-03-30 2021-07-09 上海联影医疗科技股份有限公司 主动降噪系统和具有其的医疗设备
CN113436599A (zh) * 2021-05-31 2021-09-24 新东方教育科技集团有限公司 信号处理方法、装置、存储介质及电子设备

Also Published As

Publication number Publication date
CN111583896B (zh) 2023-09-08

Similar Documents

Publication Publication Date Title
CN111583896B (zh) 一种多通道有源降噪头靠的降噪方法
JP6216096B2 (ja) 雑音減衰のためのマイクロホン配置のシステムおよび方法
CN109141620B (zh) 声源分离信息检测装置、机器人、声源分离信息检测方法和存储介质
US9848273B1 (en) Head related transfer function individualization for hearing device
EP3435372B1 (en) Silent zone generation
CN110719550B (zh) 一种双通道有源降噪头靠的虚拟传声器优化设计方法
CN113692750A (zh) 使用声音场景分析和波束形成的声传递函数个性化
US10070224B1 (en) Crosstalk cancellation for bone conduction transducers
WO2020069275A2 (en) Binaural sound source localization
CN108538304A (zh) 车内噪声主动控制系统
CN111599337B (zh) 一种外放式有源降噪头靠的降噪方法
US11495205B2 (en) Silent zone generation
JP6977448B2 (ja) 機器制御装置、機器制御プログラム、機器制御方法、対話装置、及びコミュニケーションシステム
JP2020170938A (ja) 頭部伝達関数学習装置および頭部伝達関数推論装置
CN114582312B (zh) 一种车内抗干扰自适应路噪主动控制方法及控制系统
US20220355063A1 (en) Hearing assistance devices with motion sickness prevention and mitigation features
CN114080820A (zh) 用于选择传感器阵列的声学传感器子集的方法及其系统
JP6879144B2 (ja) 機器制御装置、機器制御プログラム、機器制御方法、対話装置、及びコミュニケーションシステム
WO2020170789A1 (ja) ノイズキャンセル信号生成装置および方法、並びにプログラム
CN114360482A (zh) 一种降噪方法、装置、设备及介质
TW202315427A (zh) 用於空間化虛擬聲音來源之音訊系統
JP2008312113A (ja) 頭部伝達関数補間装置
US20230027663A1 (en) Audio method and system for a seat headrest
WO2023050253A1 (zh) 一种噪声控制的方法和装置
CN113168842B (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