CN112562627B - 前馈滤波器设计方法、主动降噪方法、系统及电子设备 - Google Patents

前馈滤波器设计方法、主动降噪方法、系统及电子设备 Download PDF

Info

Publication number
CN112562627B
CN112562627B CN202011374437.3A CN202011374437A CN112562627B CN 112562627 B CN112562627 B CN 112562627B CN 202011374437 A CN202011374437 A CN 202011374437A CN 112562627 B CN112562627 B CN 112562627B
Authority
CN
China
Prior art keywords
secondary path
noise
noise reduction
signal
filter
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.)
Active
Application number
CN202011374437.3A
Other languages
English (en)
Other versions
CN112562627A (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.)
Shenzhen Bailing Acoustics Co ltd
Original Assignee
Shenzhen Bailing Acoustics 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 Shenzhen Bailing Acoustics Co ltd filed Critical Shenzhen Bailing Acoustics Co ltd
Priority to CN202011374437.3A priority Critical patent/CN112562627B/zh
Publication of CN112562627A publication Critical patent/CN112562627A/zh
Application granted granted Critical
Publication of CN112562627B publication Critical patent/CN112562627B/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
    • 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/17815—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 reference signals and the error signals, i.e. primary 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/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/1785—Methods, e.g. algorithms; Devices
    • G10K11/17853—Methods, e.g. algorithms; Devices of the 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/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

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

本申请适用于主动降噪技术领域，提供了一种前馈滤波器设计方法、主动降噪方法、系统及电子设备。前馈滤波器设计方法应用于主动降噪系统，前馈滤波器设计方法包括：采用内部扬声器播放全带白噪声，获取内部扬声器发送的第一次级路径输入信号、和误差麦克风接收到的第一次级路径输出信号；根据第一次级路径输入信号、和第一次级路径输出信号，对次级路径进行维纳估计，获得次级路径估计；主动降噪系统接收外部扬声器播放的带限白噪声，获取参考麦克风接收到的外部噪声信号、和误差麦克风接收到的第二次级路径输出信号；根据次级路径估计、外部噪声信号、和第二次级路径输入信号，获得前馈滤波器参数。获得的滤波器运算复杂度低、易于实现。

Description

前馈滤波器设计方法、主动降噪方法、系统及电子设备
技术领域
本申请属于主动降噪技术领域，尤其涉及一种前馈滤波器设计方法、主动降噪方法、系统及电子设备。
背景技术
噪声是发声体做无规则震动时产生的声音。噪声不仅干扰人们日常的工作生活，严重时甚至会损伤听力影响身体健康。
有源噪声控制是一种通过引入可控次级声源来抑制噪声的方法，这种方法也称为主动噪声控制(Active noise control，ANC)。
由于外界声源与环境具有时变性，ANC系统必须不断调整滤波器系数实现误差信号的最小化才能达到良好的降噪效果，但由于需要不断更新和迭代系数对数字信号处理(Digital Signal Processing，DSP)硬件平台要求较高，并且运算复杂度高、缺乏鲁棒性，想要在实际系统中应用存在困难。
发明内容
本申请实施例提供了一种前馈滤波器设计方法、主动降噪方法、系统及电子设备，通过采集数据实现对内部扬声器到误差麦克风的次级路径的维纳估计，利用次级路径求解固定前馈滤波器系数，获得的滤波器运算复杂度低、易于实现。
第一方面，本申请实施例提供了一种前馈滤波器设计方法，包括：
应用于主动降噪系统，所述主动降噪系统包括参考麦克风、内部扬声器、误差麦克风；
所述前馈滤波器设计方法包括：
采用所述内部扬声器播放全带白噪声，获取所述内部扬声器发送的第一次级路径输入信号ys(n)、和所述误差麦克风接收到的第一次级路径输出信号y′s(n)；根据所述第一次级路径输入信号ys(n)、和所述第一次级路径输出信号y′s(n)，对所述次级路径进行维纳估计，获得次级路径估计s′(n)；
所述主动降噪系统接收外部扬声器播放的带限白噪声，获取所述参考麦克风接收到的外部噪声信号x(n)、和所述误差麦克风接收到的第二次级路径输出信号y′p(n)；根据所述次级路径估计s′(n)、所述外部噪声信号x(n)、和所述第二次级路径输入信号y′p(n)，获得前馈滤波器参数w(n)。
应理解，本申请实施例通过采集数据实现对内部扬声器到误差麦克风的次级路径的维纳估计，利用次级路径求解固定前馈滤波器系数，获得的滤波器运算复杂度低、易于实现。
第二方面，本申请实施例提供了一种主动降噪方法，包括：
应用于主动降噪系统，所述主动降噪系统包括：参考麦克风、内部扬声器、和前馈滤波器；
所述前馈滤波器通过所述参考麦克风接收外部噪声信号后，对所述外部噪声信号进行滤波，并将滤波后的信号进行相位反转，获得噪声抵消信号；
通过所述内部扬声器播放所述噪声抵消信号；
所述前馈滤波器为采用上述第一方面所述前馈滤波器设计方法设计的固定参数滤波器。
第三方面，本申请实施例提供了一种主动降噪系统，包括：
参考麦克风、内部扬声器、和前馈滤波器；所述前馈滤波器为采用上述第一方面所述前馈滤波器设计方法设计的固定参数滤波器。
第四方面，本申请实施例提供了一种电子设备，包括：上述第三方面所述的主动降噪系统。
第五方面，本申请实施例提供了一种计算机可读存储介质，包括：所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述第二方面所述的方法步骤。
第六方面，本申请实施例提供了一种计算机程序产品，当计算机程序产品在电子设备上运行时，使得电子设备执行上述第二方面所述的方法步骤。
可以理解的是，上述第二方面至第六方面的有益效果可以参见上述第一方面中的相关描述，在此不再赘述。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一实施例提供的主动降噪系统示意图；
图2是本申请一实施例提供的次级路径辨示意图
图3是本申请一实施例提供的初级路径辨识示意图
图4是本申请一实施例提供的一种前馈滤波器设计方法的流程示意图；
图5是本申请另一实施例提供的一种前馈滤波器设计方法的流程示意图；
图6是本申请另一实施例提供的一种前馈滤波器设计方法的流程示意图；
图7是本申请一实施例提供的一种主动降噪方法流程示意图；
图8是本申请一实施例提供的次级路径的幅频响应
图9是本申请一实施例提供的次级路径的相频响应
图10是本申请一实施例提供的前馈噪声控制滤波器幅频响应
图11是本申请一实施例提供的前馈噪声控制滤波器相频响应
图12是本申请一实施例提供的滤波器理论仿真降噪效果；
图13是本申请实施例提供的电子设备的结构示意图。
具体实施方式
以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本申请实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本申请的描述。
应当理解，当在本申请说明书和所附权利要求书中使用时，术语“包括”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
还应当理解，在本申请说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。
如在本申请说明书和所附权利要求书中所使用的那样，术语“如果”可以依据上下文被解释为“当...时”或“一旦”或“响应于确定”或“响应于检测到”。类似地，短语“如果确定”或“如果检测到[所描述条件或事件]”可以依据上下文被解释为意指“一旦确定”或“响应于确定”或“一旦检测到[所描述条件或事件]”或“响应于检测到[所描述条件或事件]”。
另外，在本申请说明书和所附权利要求书的描述中，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
在本申请说明书中描述的参考“一个实施例”或“一些实施例”等意味着在本申请的一个或多个实施例中包括结合该实施例描述的特定特征、结构或特点。由此，在本说明书中的不同之处出现的语句“在一个实施例中”、“在一些实施例中”、“在其他一些实施例中”、“在另外一些实施例中”等不是必然都参考相同的实施例，而是意味着“一个或多个但不是所有的实施例”，除非是以其他方式另外特别强调。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”，除非是以其他方式另外特别强调。
噪声是发声体做无规则震动时产生的声音。由于工业化进程的飞速发展，人类创造出大量机械设备，给社会带来了繁荣和进步，但各类设备的频繁使用也带来噪声污染日趋严重的问题。其中噪声污染的具体来源包括交通噪声(机动车辆、火车、飞机等产生的噪声)、工业噪声(工厂中各类设备产生的噪声)和建筑噪声(建筑机械产生的噪声)等。这些噪声不仅干扰人们日常的工作生活，严重时甚至会损伤听力影响身体健康。
目前常用的噪声控制方法主要分为两类：无源噪声控制和有源噪声控制。无源噪声控制又称被动噪声控制，主要通过声学材料或者结构实现对声波的吸收或控制。这种噪声控制技术对高频噪声的消除有着良好效果，但是由于低频噪声波长较长、穿透力强，被动噪声控制的方法在低频段内降噪增益有限。与无源噪声控制相比，有源噪声控制是一种通过引入可控次级声源来抑制噪声的方法。引入可控次级声源是指生成与原始噪声幅度相等相位相反的噪声，实现噪声的抵消。因此这种方法也称为主动噪声控制(Active noisecontrol，ANC)。主动噪声控制用于效消除低频噪声，与被动噪声控制相结合可在较宽频带范围内可达到良好降噪效果。
由于外界声源与环境具有时变性，ANC系统必须不断调整滤波器系数实现误差信号的最小化才能达到良好的降噪效果。近些年来将自适应滤波算法，例如最小均方算法(Least Mean Square,LMS)和FxLMS(Filter-x LMS)自适应算法，引入到ANC系统中成为人们的主要研究方向。这些自适应算法应用在ANC系统中具有良好降噪特性，但由于需要不断更新和迭代系数对DSP硬件平台要求较高，并且运算复杂度高、缺乏鲁棒性，想要在实际系统中应用存在困难。
针对以上问题本申请提出一种用于主动降噪系统的前馈滤波器设计方法，应用于主动降噪系统，该主动降噪系统系统主要组成部分包含参考麦克风、内部扬声器和误差麦克风。通过采集数据实现对次级路径(内置扬声器到误差麦克风的信道)的维纳估计，最后利用次级路径求解前馈滤波器系数。本申请提供的前馈滤波器算法具有运算复杂度低、易于实现等优点。本申请既可用于耳机，也可用于智能家居车载系统等需要实现主动噪声控制的系统。
在说明本申请实施例提供的前馈滤波器设计方法之前，为方便对本申请实施例的理解，下面结合图1至图3对本申请实施例提供的前馈滤波器设计方法的原理以及本申请实施例中涉及的相关概念进行说明。
图1示出的是本申请实施例提供的一种主动降噪系统100的示意图。该主动降噪系统也可以称为前馈主动噪声控制系统。
如图1所示主动降噪系统100包括：参考麦克风101、内部扬声器102、误差麦克风103、和前馈滤波器104。
前馈滤波器104通过所述参考麦克风101接收外部噪声信号后，对所述外部噪声信号进行滤波，并将滤波后的信号进行相位反转，获得噪声抵消信号；通过所述内部扬声器102播放所述噪声抵消信号；从达到现降低外部噪声的目的。所述前馈滤波器104可以为固定参数滤波器。该前馈滤波器可以为专用滤波器硬件实体，也可以为采用通用处理器或数字信号处理器运行计算机程序实现的滤波器。
在Z变换域上，假设从参考麦克风到达误差麦克风的信道是初级路径，初级路径表示为P(z)，误差麦克风接收到的空间噪声用D(z)表示，因此有
D(z)＝P(z)X(z)； (1)
其中，X(z)为参考麦克风接收的信号。前馈主动噪声控制通过内部扬声器播放与空间噪声D(z)反向的信号进行抵消，得到误差信号E(z)，即，
E(z)＝D(z)Y′(z)； (2)
其中，Y′(z)为次级路径S(z)的输入信号Y(z)通过次级路径到达误差麦克风的信号。
Y′(z)＝S(z)Y(z)； (3)
Y(z)为次级路径输入信号，在主动降噪系统工作时Y(z)由参考麦克风的接收信号得到，在设计前馈滤波器时，Y(z)由内部扬声器产生。
Y(z)＝W(z)X(z)； (4)
这里W(z)为估计的前馈滤波器。我们期望通过算法尽量抑制到达参考麦克风的空间噪声。
P(z)X(z)＝S(z)W(z)X(z)； (5)
所以理想的前馈滤波器为
Figure BDA0002807794930000071
对于固定滤波器参数的主动噪声控制，在本申请实施例中，首先辨识次级路径。在时域，参考麦克风接收信号x(n)、误差麦克风接收的刺激路径的输出信号y′(n)和扬声器输出的次级路径输入信号y(n)的信号模型分别表示为，
x(n)＝[x(n)x(n-1)…x(n-L+1)]T；
y′(n)＝[y′(n)y′(n-1)…y′(n-L+1)]T；
y(n)＝[y(n)y(n-1)…y(n-L+1)]T；
次级路径s和初级路径p分别表示为，
s＝[s0s1…sL-1]T；
p＝[p0p1…pL-1]T；
滤波器w在时域表示形式为
w＝[w0w1…wL-1]T；
其中L表示滤波器长度，在一些实施例中L＝1024。
其中，初级路径(Primary path)又称为主路径、主信道、初级通道。次级路径(Secondary path)又称为次级信道、次级通道。
图2示出的是次级路径辨示意图。如图2所示，可以采用估计的线性滤波器
Figure BDA0002807794930000081
来逼近次级路径S(z)。误差麦克风接收的次级路径输出信号y′(n)可以表示为：
y′(n)＝sTy(n)； (7)
令
Figure BDA0002807794930000082
为
Figure BDA0002807794930000083
的时域表示，则：
Figure BDA0002807794930000084
本申请实施例中，通过维纳估计得到次级路径估计
Figure BDA0002807794930000085
Figure BDA0002807794930000086
其中，
Ryy＝E[y(n)yT(n)]； (9)
ryy′＝E[y(n)y′(n)]； (10)
图3示出的是初级路径辨识示意图。参考图3，在得到次级路径估计
Figure BDA0002807794930000087
后，接下来根据以上得到的
Figure BDA0002807794930000088
求解滤波器w。图3中的e(n)为误差信号E(z)的时域表示。
参考图3，空间噪声D(z)的Z变换域表示为：
D(z)＝W(z)X′(z)； (11)
Z变换域的外部噪声估计X′(z)可以表示为：
X′(z)＝S(z)X(z)； (12)
所以对参考麦克风接收到外部噪声信号x(n)进行滤波，得到外部噪声估计x′(n)：
Figure BDA0002807794930000089
x′(n)的向量表示为：x′(n)＝[x′(n)x′(n-1)…x′(n-L+1)]T。
参考图3，可以得到空间噪声d(n)，
d(n)＝wTx′(n)； (14)
由此，可以计算得出前馈控制滤波器w，
w＝Rx′x′ -1rx′d； (15)
其中，
Rx′x′＝E[x′(n)x′T(n)]； (16)
rx′d＝E[x′(n)d(n)]； (17)
本申请实施例中的E[]表示求取数学期望。
主动降噪效果采用误差麦开启算法前后功率比值进行评价。记降噪开启和关闭时的误差麦克风信号分别为e(n),c(n)分别采用平均周期图法估计两信号的功率谱，
Pe(f)＝jFFT[e(n)]j2； (18)
Pc(f)＝jFFT[c(n)]j2； (19)
计算出功率谱值，并换算成dB值，
Figure BDA0002807794930000091
由于主动降噪算法主要针对低频噪声，因此取频率f的范围为20Hz<f<1kHz的SNR(f)，求平均计算平均降噪量，
NR＝E[JNR(F)j]；(21)
其中，FFT[]为快速傅里叶变换。J为虚数单位。
图4是本申请实施例提供的一种前馈滤波器设计方法。如图4所示，基于以上的理论推导过程，本申请实施例提供的前馈滤波器设计方法，包括步骤S1和步骤S2，各个步骤的具体实现原理如下：
S1，采用所述内部扬声器播放全带白噪声，获取所述内部扬声器发送的第一次级路径输入信号ys(n)、和所述误差麦克风接收到的第一次级路径输出信号y′s(n)；根据所述第一次级路径输入信号ys(n)、和所述第一次级路径输出信号y′s(n)，对所述次级路径进行维纳估计，获得次级路径估计s′(n)。
可选的，如图5所示，本申请实施例提供的一种实施方式中，步骤S1，采用所述内部扬声器播放全带白噪声，获取第一次级路径输入信号ys(n)、和第一次级路径输出信号y′s(n)；根据所述第一次级路径输入信号ys(n)、和所述第一次级路径输出信号y′s(n)，对所述次级路径进行维纳估计，获得次级路径估计s′(n)，包括步骤S11至步骤S14：
步骤S11，通过所述内部扬声器播放全带白噪声，测得所述内部扬声器产生的第一次级路径输入信号ys(n)、和所述误差麦克风接收到的第一次级路径输出信号y′s(n)；
步骤S12，采用Ryy＝E[y(n)yT(n)]获得所述第一次级路径输入信号ys(n)的自相关矩阵Ryy；其中E[]表示求取数学期望；
步骤S13，采用ryy′＝E[y(n)y′(n)]获得所述第一次级路径输入信号ys(n)和所述第一次级路径输出信号y′s(n)的互相关向量ryy′；
步骤S14，采用
Figure BDA0002807794930000101
获得所述次级路径估计s′(n)。
可选的，在图5所示的实施例的基础上，将步骤S11至步骤S14重复K次，获得K个次级路径估计s′(n)；其中K为大于或等于2的整数；
去除错误的次级路径估计后，将剩余的次级路径的估计取均值，作为次级路径估计s′(n)。
S2，所述主动降噪系统接收外部扬声器播放的带限白噪声，获取所述参考麦克风接收到的外部噪声信号x(n)、和所述误差麦克风接收到的第二次级路径输出信号y′p(n)；根据所述次级路径估计s′(n)、所述外部噪声信号x(n)、和所述第二次级路径输入信号y′p(n)，获得前馈滤波器参数w(n)。
可选的，如图6所示，本申请实施例提供的一种实施方式中，步骤S2，采用外部扬声器播放带限白噪声，获取外部噪声信号x(n)、和第二次级路径输出信号y′p(n)；根据所述次级路径估计s′(n)、所述外部噪声信号x(n)、和所述第二次级路径输入信号y′p(n)，获得所述前馈滤波器参数w(n)，包括步骤S21至步骤S26：
步骤S21，所述主动降噪系统接收外部扬声器播放的带限白噪声，；测得所述参考麦克风的接收到的外部噪声信号x(n)，和所述误差麦克风接收到的第二次级路径输出信号y′p(n)；
步骤S22，根据外部噪声信号x(n)和次级路径估计s′(n)，采用x′(n)＝s′(n)Tx(n)得到所述误差麦克风接收到的空间噪声x′(n)；
步骤S23，令空间噪声d(n)＝y′p(n)；
步骤S24，采用Rx′x′＝E[x′(n)x′T(n)]获得所述外部噪声估计x′(n)的自相关矩阵Rx′x′；
步骤S25，采用rx′d＝E[x′(n)d(n)]获得所述外部噪声估计x′(n)和所述空间噪声d(n)的互相关向量rx′d；
步骤S26，采用w(n)＝Rx′x′ -1rx′d，获得所述前馈滤波器参数w(n)。
可选的，在图6所示的实施例的基础上，针对所述外部扬声器的N个摆放位置中的每个摆放位置，和每个所述外部扬声器的摆放位置对应的M个类型的次级路径，将步骤S21至步骤S26重复M×N次，获得的M×N个滤波器参数；其中M和N为大于或等于2的整数；在一些实施例中，外部扬声器的外部摆放位置包括相对于主动降噪系统的距离和朝向。
去除M×N个第一滤波器参数中的错误数据后，将剩余的滤波器参数取均值，获得前馈噪声控制滤波器参数。本申请实施例中所称错误数据可以为明显偏离其他数据的数据，可以通过设置方差阈值，将方差明显大于方差阈值的数据看做错误数据。通过多组测试数据求平均以提高鲁棒性，并且对不同位置噪声和不同佩戴人员均有良好降噪效果
可选的，在以上各个实施例的基础上，所述带限白噪声为0至2kHz的带限白噪声。
图7示出了本申请实施例提供的主动降噪方法，应用于上述图1所示的主动降噪系统，可由所述主动降噪系统的软件和/或硬件实现。如图7所示，该方法包括步骤S710至S720。各个步骤的具体实现原理如下：
所述主动降噪系统包括：参考麦克风、内部扬声器、和前馈滤波器；
S710，所述前馈滤波器通过所述参考麦克风接收外部噪声信号后，对所述外部噪声信号进行滤波，并将滤波后的信号进行相位反转，获得噪声抵消信号。
S720，通过所述内部扬声器播放所述噪声抵消信号。
所述前馈滤波器为采用以上各个实施例提供的所述前馈滤波器设计方法设计的的固定参数滤波器。该前馈滤波器可以为专用滤波器硬件实体，也可以为采用通用处理器或数字信号处理器运行计算机程序实现的滤波器。
为了更好地展示本申请的效果，给出如下具体的实例来验证本申请实施例提供的前馈滤波器设计方法设计的滤波器的效果。
实验条件如下：测试信号采用0-2kHz带限白噪声，四名测试人员，每人采三组数据，采样率768kHz，扬声器有三种不同摆放位置，分别是测试人员正前方1m，正前方2m，左侧1m。
由四名测试人员右耳佩戴耳机在消声室中进行测试，每人采三组数据。在本示例中，四名测试人员右耳佩戴耳机在消声室中进行测试提供了四种不同的次级路径。在将本申请实施例应用在汽车或建筑环境内，可以参考次级路径的定义和本申请的实施例进行设置。
根据测试数据得到如图8所示的次级路径的幅频响应，和如图9所示的次级路径的相频响应。图8至图12中横坐标为频率(HZ)，纵坐标为幅度(dB)。
接下来，按照三种不同外部扬声器摆放位置，由四名测试人员佩戴耳机在消声室中进行测试，每人采三组数据。计算得到平均后的前馈噪声控制滤波器，图10示出的是前馈噪声控制滤波器幅频响应，图11示出的是前馈噪声控制滤波器相频响应。图12示出的是相应的滤波器理论仿真降噪效果。如表1所示，采用本申请实施例设计的滤波器的算法取得了非常好的降噪效果，仿真最大降噪量达到了48.71dB,在较低的频点如100Hz处，降噪量也能达到30.98dB的降噪量。如果我们同时得到次级路径与初级路径(参考麦到误差麦的信道)的维纳估计，根据公式(6)可计算得到前馈滤波器系数，我们把这种方法记为方法一。本申请与方法一相比，只进行一次信道估计引入误差更小，因此具有更好的降噪效果，具体降噪量对比如表2所示。并且由于本申请实施例设计实施固定系数滤波器，获得的滤波器运算复杂度低、易于实现。
表1：
0-2kHz平均降噪量 0-1kHz平均降噪量 最大降噪量 10kHz平均降噪量
33.97dB 40.09dB 48.71dB 30.98dB
表2：
方法一0-2kHz平均降噪量 本申请实施例0-2kHz平均降噪量
24.37dB 33.97dB
本申请的一实施例提供了一种电子设备，该电子设备可以为主动降噪耳机、车载降噪设备或室内降噪设备。该电子设备包括上述各个实施例提供的主动降噪系统。
图13为本申请一实施例提供的电子设备的结构示意图。如图13所示，该实施例的电子设备D10包括：至少一个处理器D100(图13中仅示出一个)处理器、存储器D101以及存储在所述存储器D101中并可在所述至少一个处理器D100上运行的计算机程序D102，所述处理器D100执行所述计算机程序D102时实现上述任意各个方法实施例中的步骤。该电子设备中还包含上述主动降噪系统的参考麦克风D103、内部扬声器D104和误差麦克风D105。该电子设备还包括前馈滤波器，该前馈滤波器可以为专用滤波器硬件实体，也可以为采用通用处理器或数字信号处理器运行计算机程序实现的滤波器。如图13所示，该电子设备中的各个组成部件/单元可以通过总线D110通信，总线D110的类型不做具体限定。
所称处理器D100可以是中央处理单元(Central Processing Unit，CPU)，该处理器D100还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器D101在一些实施例中可以是所述电子设备D10的内部存储单元，例如电子设备D10的硬盘或内存。所述存储器D101在另一些实施例中也可以是所述电子设备D10的外部存储设备，例如所述电子设备D10上配备的插接式硬盘，智能存储卡(Smart MediaCard,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。进一步地，所述存储器D101还可以既包括所述电子设备D10的内部存储单元也包括外部存储设备。所述存储器D101用于存储操作系统、应用程序、引导装载程序(BootLoader)、数据以及其他程序等，例如所述计算机程序的程序代码等。所述存储器D101还可以用于暂时地存储已经输出或者将要输出的数据。
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。
需要说明的是，上述装置/单元之间的信息交互、执行过程等内容，由于与本申请方法实施例基于同一构思，其具体功能及带来的技术效果，具体可参见方法实施例部分，此处不再赘述。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
本申请实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现可实现上述各个方法实施例中的步骤。
本申请实施例提供了一种计算机程序产品，当计算机程序产品在电子设备上运行时，使得电子设备执行时实现可实现上述各个方法实施例中的步骤。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实现上述实施例方法中的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质至少可以包括：能够将计算机程序代码携带到拍照装置/终端设备的任何实体或装置、记录介质、计算机存储器、只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random AccessMemory，RAM)、电载波信号、电信信号以及软件分发介质。例如U盘、移动硬盘、磁碟或者光盘等。在某些司法管辖区，根据立法和专利实践，计算机可读介质不可以是电载波信号和电信信号。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
在本申请所提供的实施例中，应该理解到，所揭露的装置/网络设备和方法，可以通过其它的方式实现。例如，以上所描述的装置/网络设备实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。

Claims (8)

1.一种前馈滤波器设计方法，其特征在于，应用于主动降噪系统，所述主动降噪系统包括参考麦克风、内部扬声器、误差麦克风；
所述前馈滤波器设计方法包括：
采用所述内部扬声器播放全带白噪声，获取所述内部扬声器发送的第一次级路径输入信号ys(n)、和所述误差麦克风接收到的第一次级路径输出信号y′s(n)；根据所述第一次级路径输入信号ys(n)、和所述第一次级路径输出信号y′s(n)，对所述次级路径进行维纳估计，获得次级路径估计s′(n)；
所述主动降噪系统接收外部扬声器播放的带限白噪声，获取所述参考麦克风接收到的外部噪声信号x(n)、和所述误差麦克风接收到的第二次级路径输出信号y′p(n)；根据所述次级路径估计s′(n)、所述外部噪声信号x(n)、和所述第二次级路径输出 信号y′p(n)，获得固定参数滤波器的前馈滤波器参数w(n)；
所述主动降噪系统接收外部扬声器播放的带限白噪声，获取所述参考麦克风接收到的外部噪声信号x(n)、和所述误差麦克风接收到的第二次级路径输出信号y′p(n)；根据所述次级路径估计s′(n)、所述外部噪声信号x(n)、和所述第二次级路径输出 信号y′p(n)，获得前馈滤波器参数w(n)，包括：
步骤S21，所述主动降噪系统接收外部扬声器播放的带限白噪声；测得所述参考麦克风的接收到的外部噪声信号x(n)，和所述误差麦克风接收到的第二次级路径输出信号y′p(n)；
步骤S22，根据外部噪声信号x(n)和次级路径估计s′(n)，采用x′(n)＝s′(n)Tx(n)得到所述误差麦克风接收到的外部噪声估计x′(n)；
步骤S23，令空间噪声d(n)＝y′p(n)；
步骤S24，采用Rx′x′＝E[x′(n)x′T(n)]，获得所述外部噪声估计x′(n)的自相关矩阵Rx′x′；
步骤S25，采用rx′d＝E[x′(n)d(n)]获得所述外部噪声估计x′(n)和所述空间噪声的互相关向量rx′d；
步骤S26，采用w(n)＝Rx′x′ -1rx′d，获得所述前馈滤波器参数w(n)；
所述前馈滤波器设计方法还包括：
针对所述外部扬声器的N个摆放位置中的每个摆放位置，和每个所述外部扬声器的摆放位置对应的M个类型的次级路径，将步骤S21至步骤S26重复M×N次，获得的M×N个滤波器参数；其中M和N为大于或等于2的整数；去除M×N个第一滤波器参数中的错误数据后，将剩余的滤波器参数取均值，获得前馈噪声控制滤波器参数。
2.如权利要求1所述的前馈滤波器设计方法，其特征在于，采用所述内部扬声器播放全带白噪声，获取所述内部扬声器发送的第一次级路径输入信号ys(n)、和所述误差麦克风接收到的第一次级路径输出信号y′s(n)；根据所述第一次级路径输入信号ys(n)、和所述第一次级路径输出信号y′s(n)，对所述次级路径进行维纳估计，获得次级路径估计s′(n)，包括：
步骤S11，通过所述内部扬声器播放全带白噪声，测得所述内部扬声器产生的第一次级路径输入信号ys(n)、和所述误差麦克风接收到的第一次级路径输出信号y′s(n)；
步骤S12，采用Ryy＝E[y(n)yT(n)]获得所述第一次级路径输入信号ys(n)的自相关矩阵Ryy；其中E[]表示求取数学期望；
步骤S13，采用ryy′＝E[y(n)y′(n)]获得所述第一次级路径输入信号ys(n)和所述第一次级路径输出信号y′s(n)的互相关向量ryy′；
步骤S14，采用
Figure FDA0003115256930000021
获得所述次级路径估计s′(n)。
3.如权利要求2所述的前馈滤波器设计方法，其特征在于，还包括：
将步骤S11至步骤S14重复K次，获得K个次级路径估计s′(n)；其中K 为大于或等于2的整数；
去除错误的次级路径估计后，将剩余的次级路径的估计取均值，作为次级路径估计s′(n)。
4.如权利要求所1至3任一项所述的馈滤波器设计方法，其特征在于，所述带限白噪声为0至2kHz的带限白噪声。
5.一种主动降噪方法，其特征在于，应用于主动降噪系统，所述主动降噪系统包括：参考麦克风、内部扬声器、和前馈滤波器；
所述前馈滤波器通过所述参考麦克风接收外部噪声信号后，对所述外部噪声信号进行滤波，并将滤波后的信号进行相位反转，获得噪声抵消信号；
通过所述内部扬声器播放所述噪声抵消信号；
所述前馈滤波器为采用权利要求1至3任一项所述前馈滤波器设计方法设计的固定参数滤波器。
6.一种主动降噪系统，其特征在于，包括：参考麦克风、内部扬声器、和前馈滤波器；所述前馈滤波器为采用权利要求1至3任一项所述前馈滤波器设计方法设计的固定参数滤波器。
7.一种电子设备，其特征在于，包括权利要求6所述的主动降噪系统。
8.如权利要求7所述的电子设备，其特征在于，所述电子设备包括主动降噪耳机、或车载降噪设备或室内降噪设备。
CN202011374437.3A 2020-11-30 2020-11-30 前馈滤波器设计方法、主动降噪方法、系统及电子设备 Active CN112562627B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011374437.3A CN112562627B (zh) 2020-11-30 2020-11-30 前馈滤波器设计方法、主动降噪方法、系统及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011374437.3A CN112562627B (zh) 2020-11-30 2020-11-30 前馈滤波器设计方法、主动降噪方法、系统及电子设备

Publications (2)

Publication Number Publication Date
CN112562627A CN112562627A (zh) 2021-03-26
CN112562627B true CN112562627B (zh) 2021-08-31

Family

ID=75045380

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011374437.3A Active CN112562627B (zh) 2020-11-30 2020-11-30 前馈滤波器设计方法、主动降噪方法、系统及电子设备

Country Status (1)

Country Link
CN (1) CN112562627B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113299263B (zh) * 2021-05-21 2024-05-24 北京安声浩朗科技有限公司 声学路径确定方法、装置、可读存储介质及主动降噪耳机
CN113345400B (zh) * 2021-05-31 2024-06-21 锐迪科微电子（上海）有限公司 可穿戴设备的主动降噪系统的校准方法及装置、存储介质、终端
CN113315490A (zh) * 2021-06-04 2021-08-27 哈尔滨工程大学 一种基于后滤波器的滤波参考仿射投影符号主动控制系统
CN113299265B (zh) * 2021-07-26 2022-09-27 北京安声浩朗科技有限公司 主动降噪方法、装置及主动降噪耳机
CN113421540B (zh) * 2021-07-26 2023-10-31 北京安声浩朗科技有限公司 主动降噪方法、主动降噪装置及半入耳式主动降噪耳机
CN113707120A (zh) * 2021-08-02 2021-11-26 杭州萤石软件有限公司 声音信号的处理系统、方法及装置
CN113766387B (zh) * 2021-09-30 2024-04-09 展讯通信（上海）有限公司 反馈滤波器的设计方法、主动降噪方法、系统及电子设备
CN113920973B (zh) * 2021-10-11 2024-09-10 聊城大学 建筑空间降噪方法、装置和系统
CN116208879B (zh) * 2021-11-30 2024-08-09 华为技术有限公司 具有主动降噪功能的耳机及主动降噪方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101203905A (zh) * 2005-04-22 2008-06-18 阿诺克系统股份公司 在主动降噪系统中模仿次级路径的方法
CN105284126A (zh) * 2013-04-15 2016-01-27 美国思睿逻辑有限公司 包括自适应消噪系统的系数的动态偏置的自适应消噪系统及方法
CN105788604A (zh) * 2016-04-07 2016-07-20 虞安波 一种基于fxlms的优化主动降噪方法
CN106839386A (zh) * 2017-01-22 2017-06-13 会听声学科技（北京）有限公司 一种用于塔式空气净化器的主动降噪系统及方法
CN110933554A (zh) * 2019-12-13 2020-03-27 恒玄科技(上海)股份有限公司 主动降噪方法、系统以及耳机

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8718291B2 (en) * 2011-01-05 2014-05-06 Cambridge Silicon Radio Limited ANC for BT headphones
US8693700B2 (en) * 2011-03-31 2014-04-08 Bose Corporation Adaptive feed-forward noise reduction
US9462376B2 (en) * 2013-04-16 2016-10-04 Cirrus Logic, Inc. Systems and methods for hybrid adaptive noise cancellation
FR3039311B1 (fr) * 2015-07-24 2017-08-18 Orosound Dispositif de controle actif de bruit
US9578415B1 (en) * 2015-08-21 2017-02-21 Cirrus Logic, Inc. Hybrid adaptive noise cancellation system with filtered error microphone signal
CN109545181A (zh) * 2018-12-13 2019-03-29 四川长虹电器股份有限公司 一种自适应数字主动降噪架构
CN110996209B (zh) * 2019-12-13 2022-02-15 恒玄科技(上海)股份有限公司 主动降噪方法、系统以及耳机

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101203905A (zh) * 2005-04-22 2008-06-18 阿诺克系统股份公司 在主动降噪系统中模仿次级路径的方法
CN105284126A (zh) * 2013-04-15 2016-01-27 美国思睿逻辑有限公司 包括自适应消噪系统的系数的动态偏置的自适应消噪系统及方法
CN105788604A (zh) * 2016-04-07 2016-07-20 虞安波 一种基于fxlms的优化主动降噪方法
CN106839386A (zh) * 2017-01-22 2017-06-13 会听声学科技（北京）有限公司 一种用于塔式空气净化器的主动降噪系统及方法
CN110933554A (zh) * 2019-12-13 2020-03-27 恒玄科技(上海)股份有限公司 主动降噪方法、系统以及耳机

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
车内ANC系统次级通道辨识及仿真研究;陆森林等;《重庆理工大学学报(自然科学)》;20160331;第22-27页 *

Also Published As

Publication number Publication date
CN112562627A (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
CN112562627B (zh) 前馈滤波器设计方法、主动降噪方法、系统及电子设备
CN112562626B (zh) 混合降噪滤波器的设计方法、降噪方法、系统及电子设备
KR101798120B1 (ko) 능동 잡음 제거와 지각적 잡음 보상의 결합에 의해 음향 재생의 지각 품질을 향상시키기 위한 장치 및 방법
RU2698639C2 (ru) Субполосный алгоритм с пороговой величиной для устойчивой широкополосной системы активного шумоподавления
CN102461206B (zh) 便携式通信设备和其中处理信号的方法
CN112562625B (zh) 主动降噪滤波器设计方法、降噪方法、系统及电子设备
US20150281840A1 (en) Single-channel, binaural and multi-channel dereverberation
CN111971975B (zh) 主动降噪的方法、系统、电子设备和芯片
CN108174320A (zh) 一种自适应耳道主动降噪耳机及自适应耳道主动降噪方法
CN112562624B (zh) 主动降噪滤波器设计方法、降噪方法、系统及电子设备
CN110021305B (zh) 一种音频滤波方法、音频滤波装置及可穿戴设备
CN113593612B (zh) 语音信号处理方法、设备、介质及计算机程序产品
CN106303783A (zh) 降噪方法和装置
CN110956973A (zh) 一种回声消除方法、装置及智能终端
CN111951818B (zh) 一种基于改进功率差噪声估计算法的双麦克风语音增强方法
CN111756906B (zh) 一种语音信号的回声抑制方法、装置和计算机可读介质
CN103916810B (zh) 一种时域声能量对比控制方法和系统
US20170148453A1 (en) Bass enhancement and separation of an audio signal into a harmonic and transient signal component
CN112201273A (zh) 一种噪声功率谱密度计算方法、系统、设备及介质
CN116405823A (zh) 一种骨传导耳机音频智能去噪增强方法
CN110708651B (zh) 一种基于分段陷波的助听器啸叫检测与抑制方法及装置
CN114584908A (zh) 助听器的声学测试方法、装置以及设备
CN113242491B (zh) 降噪处理方法、装置、设备、存储介质及程序
US11355096B1 (en) Adaptive feedback processing for consistent headphone acoustic noise cancellation
Miyazaki et al. Theoretical analysis of parametric blind spatial subtraction array and its application to speech recognition performance prediction

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