CN105165026B - 使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法 - Google Patents

使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法 Download PDF

Info

Publication number
CN105165026B
CN105165026B CN201380073406.6A CN201380073406A CN105165026B CN 105165026 B CN105165026 B CN 105165026B CN 201380073406 A CN201380073406 A CN 201380073406A CN 105165026 B CN105165026 B CN 105165026B
Authority
CN
China
Prior art keywords
information
time
noise
frequency band
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
CN201380073406.6A
Other languages
English (en)
Other versions
CN105165026A (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.)
Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV
Original Assignee
Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV
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 Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV filed Critical Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV
Publication of CN105165026A publication Critical patent/CN105165026A/zh
Application granted granted Critical
Publication of CN105165026B publication Critical patent/CN105165026B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/16—Sound input; Sound output
    • 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/18—Methods or devices for transmitting, conducting or directing sound
    • G10K11/26—Sound-focusing or directing, e.g. scanning
    • G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
    • G10K11/341—Circuits therefor
    • G10K11/346—Circuits therefor using phase variation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/02—Circuits for transducers, loudspeakers or microphones for preventing acoustic reaction, i.e. acoustic oscillatory feedback

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Otolaryngology (AREA)
  • Signal Processing (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Human Computer Interaction (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Circuit For Audible Band Transducer (AREA)
  • Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

提供了一种滤波器(100)，用于基于两个或更多个输入麦克风信号生成音频输出信号，音频输出信号包括多个音频输出信号样本。在时频域中表示音频输出信号以及两个或更多个输入麦克风信号，其中，多个音频输出信号样本中的每个被分配给多个时频段((k,n))中的时频段((k,n))。滤波器(100)包括权重生成器(110)，其被适配成为多个时频段((k,n))中的每个，接收一个或多个声源的一个或多个声音分量的到达方向信息或一个或多个声源的位置信息，并且被适配为根据所述时频段((k,n))的一个或多个声源的一个或多个声音分量的到达方向信息或者根据所述时频段((k,n))的一个或多个声源的位置信息，为所述多个时频段((k,n))中的每个，生成加权信息。而且，滤波器包括输出信号生成器(120)，其用于通过根据所述时频段((k,n))的加权信息并且根据两个或更多个输入麦克风信号中的每个的被分配给所述时频段((k,n))的音频输入样本，为多个时频段((k,n))中的每个生成分配给所述时频段((k,n))的多个音频输出信号样本中的一个，来生成音频输出信号。

Description

使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法
技术领域
本发明涉及音频信号处理，并且尤其涉及使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法。
背景技术
在现代通信系统中通常发现在噪声和混响条件下提取声源。在过去的四十年，提出了大量空间滤波技术，以完成这个任务。在观察的信号符合信号模型时并且在计算滤波器所需要的信息精确时，现有空间滤波器是最佳的。然而，实际上，通常违背信号模型，并且估计所需要的信息是一项重大的挑战。
现有空间滤波器可广义地分成线性空间滤波器(例如，参照[1,2,3,4])和参数空间滤波器(parametric spatial filter)(例如，参照[5,6,7,8])。通常，线性空间滤波器需要估计一个或多个传播向量或者期望的一个或多个源的二阶统计量(SOS)加上干扰的SOS。一些空间滤波器被设计为提取混响或去混响的单个源信号(例如，参照[9,10,11,12,13,14,15,16])，而其他空间滤波器被设计为提取两个或更多个混响源信号的和(例如，参照[17,18])。上述方法需要一个或多个期望源的方向的先验知识或仅仅期望源单独地或者同时激活的周期。
这些方法的缺点在于，不能足够快速地适应新情况，例如，在期望源激活时变激活的源移动或竞争扬声器。参数空间滤波器通常基于比较简单的信号模型，例如，在时频域内接收的信号由单个平面波加上扩散声构成，并且根据模型参数的瞬时估计，计算参数空间滤波器。参数空间滤波器的优点是高度灵活的方向响应、扩散声音和干涉源的较强抑制、以及快速适合于新情况的能力。然而，如在[19]中所示，实际上可以容易地违背基本单平面波信号模型，这显著降低参数空间滤波器的性能。应注意的是，现有技术参数空间滤波器使用所有可用的麦克风信号，来估计模型参数，而仅单个麦克风信号和实值增益用于计算最终的输出信号。用于使多个可用麦克风信号相结合以便找出增强的输出信号的扩展不直接(straight forward)。
因此，如果提供用于获得对声源的期望空间响应的改进概念，则被高度支持。
发明内容
因此，本发明的目标在于，提供用于提取声源的改进概念。由根据权利要求1所述的滤波器、根据权利要求17所述的方法以及由根据权利要求18所述的计算机程序，解决本发明的目标。
提供了一种滤波器，该滤波器用于基于两个或更多个输入麦克风信号生成音频输出信号(其包括多个音频输出信号样本)。在时频域中表示音频输出信号和两个或更多个输入麦克风信号，其中，多个音频输出信号样本中的每个被分配给多个时频段(time-frequency bin)中的一个时频段。
该滤波器包括权重生成器，权重生成器被适配成为多个时频段中的每个，接收一个或多个声源的一个或多个声音分量的到达方向信息或一个或多个声源的位置信息，并且被适配为根据所述时频段的一个或多个声源的一个或多个声音分量的到达方向信息或者根据所述时频段的一个或多个声源的位置信息，为多个时频段中的每个生成加权信息。
而且，该滤波器包括输出信号生成器，输出信号生成器用于通过为多个时频段中的每个，生成多个音频输出信号样本中的分配给所述时频段的一个音频输出信号样本，来根据所述时频段的加权信息并且根据分配给所述时频段的音频输入样本，生成两个或更多个输入麦克风信号中的每个的音频输出信号。
实施方式提供了一种空间滤波器，用于获得至多L个同时激活的声源的期望的响应。通过最小化经受L个线性约束的在滤波器的输出处的扩散加噪声功率，获得所提供的空间滤波器。与现有技术概念相反，L个约束基于瞬时窄带到达方向估计。此外，提供用于扩散噪声比/扩散功率的新型估计器，其呈现足够高的时间和光谱分辨率，以实现去混响和噪声降低。
根据一些实施方式，提供了用于获得在每个时间-频率瞬间同时激活的至多L个声源的期望的任意空间响应的概念。为此，关于声音场景的瞬时参数信息(IPI)结合到空间滤波器的设计内，产生“知情空间滤波器(informed spatial filter)”。
在一些实施方式中，例如，这种知情空间滤波器基于复杂权重(complex weight)使所有可用的麦克风信号组合，以提供增强的输出信号。
根据实施方式，例如，知情空间滤波器可实现为线性约束最小方差(LCMV)空间滤波器或者参数多通道维纳滤波器。
在一些实施方式中，例如，通过最小化经受L个线性约束的扩散加自噪声功率，获得所提供的知情空间滤波器。
在一些实施方式中，相对于现有技术，L个约束基于瞬时到达方向(DOA)估计，并且所产生的对L个DOA的响应与特定的期望方向性对应。
而且，提供用于所需要的信号和噪声统计(例如，扩散噪声比(DNR))的新型估计器，其呈现足够高的时间和光谱分辨率，例如，以减少混响和噪声。
而且，提供了一种用于基于两个或更多个输入麦克风信号生成音频输出信号的方法，所述音频输出信号包括多个音频输出信号样本。在时频域中表示音频输出信号和两个或更多个输入麦克风信号，其中，多个音频输出信号样本中的每个被分配给多个时频段中的一个时频段。该方法包括：
-为多个时频段((k,n))中的每个，接收一个或多个声源的一个或多个声音分量的到达方向信息或一个或多个声源的位置信息，
-根据所述时频段的一个或多个声源的一个或多个声音分量的到达方向信息或者根据所述时频段的一个或多个声源的位置信息，为多个时频段中的每个生成加权信息，并且
-通过为多个时频段((k,n))中的每个，生成分配给所述时频段((k,n))的多个音频输出信号样本中的一个，来根据所述时频段((k,n))的加权信息并且根据分配给所述时频段((k,n))的音频输入样本，生成两个或更多个输入麦克风信号中的每个的音频输出信号。
而且，提供了一种计算机程序，用于在计算机或信号处理器上执行时，实施上述方法。
附图说明
在下文中，参照附图，更详细地描述本发明的实施方式，其中：
图1a示出了根据实施方式的滤波器；
图1b示出了根据实施方式的滤波器的一个可能的应用场景；
图2示出了根据实施方式的滤波器以及多个麦克风；
图3示出了根据实施方式的权重生成器；
图4示出了根据实施方式的两个实例响应的幅度；
图5示出了根据实施线性约束最小方差方法的另一个实施方式的权重生成器；
图6示出了根据实施参数多通道维纳滤波器方法的进一步实施方式的权重生成器；
图7示出了根据时间和频率的真实的和估计的扩散噪声比；
图8示出了比较空间滤波器的方向性指数和白噪声增益；
图9示出了估计的到达方向以及所产生的增益；以及
图10示出了立体扬声器再现的情况的实例。
具体实施方式
图1a示出了滤波器100，其用于基于两个或更多个输入麦克风信号生成音频输出信号(其包括多个音频输出信号样本)。在时频域中表示音频输出信号以及两个或更多个输入麦克风信号，其中，多个音频输出信号样本中的每个被分配给多个时频段(k,n)中的一个时频段(k,n)。
滤波器100包括权重生成器110，其被适配成为多个时频段(k,n)中的每个，接收一个或多个声源的一个或多个声音分量的到达方向信息或一个或多个声源的位置信息，并且被适配为根据所述时频段(k,n)的一个或多个声源的一个或多个声音分量的到达方向信息或者根据所述时频段(k,n)的一个或多个声源的位置信息，为多个时频段(k,n)中的每个，生成加权信息。
而且，滤波器包括输出信号生成器120，其用于通过根据所述时频段(k,n)的加权信息并且根据分配给所述时频段(k,n)的音频输入样本，为多个时频段(k,n)中的每个生成分配给所述时频段(k,n)的所述多个音频输出信号样本中的一个，来生成两个或更多个输入麦克风信号中的每个的音频输出信号。
例如，两个或更多个输入麦克风信号中的每个包括多个音频输入样本，其中，每个音频输入样本被分配给一个时频段(k,n)，并且音频信号生成器120可被适配为根据所述时频段(k,n)的加权信息并且根据两个或更多个输入麦克风信号中的每个的音频输入样本中的一个，即，根据两个或更多个输入麦克风信号中的每个的音频输入样本中的被分配给所述时频段(k,n)的一个音频输入样本，生成多个音频输出信号样本中的被分配给所述时频段(k,n)的一个音频输出信号样本。
对于每个时频段(k,n)的要生成的每个音频输出信号样本，权重生成器110重新生成单独的加权信息。然后，音频信号生成器120基于为该时频段生成的加权信息，生成考虑的时频段(k,n)的音频输出信号样本。换言之，通过权重生成器110为要生成音频输出信号样本的每个时频段计算新加权信息。
在生成加权信息时，权重生成器110被适配为考虑一个或多个声源的信息。
例如，权重生成器110可考虑第一声源的位置。在实施方式中，权重生成器还可考虑第二声源的位置。
例如，或者，第一声源可发射具有第一声音分量的第一声波。具有第一声音分量的第一声波到达麦克风，并且权重生成器110可考虑第一声音分量的/声波的到达方向。借此，权重生成器110将关于第一声源的信息考虑在内。而且，第二声源可发射具有第二声音分量的第二声波。具有第二声音分量的第二声波到达麦克风，并且权重生成器110可考虑第二声音分量的/第二声波的到达方向。借此，权重生成器110还将关于第二声源的信息考虑在内。
图1b示出了根据实施方式的滤波器100的可能的应用场景。具有第一声音分量的第一声波由第一扬声器121(第一声源)发射并且到达第一麦克风111。考虑在第一麦克风111处的第一声音分量的到达方向(＝第一声波的到达方向)。而且，具有第二声音分量的第二声波由第二扬声器122(第二声源)发射并且到达第一麦克风111。权重生成器110还能够考虑在第一麦克风111处的第二声音分量的到达方向，以确定加权信息。而且，权重生成器还可考虑在其他麦克风(例如，麦克风112)处的声音分量的到达方向(＝声波的到达方向)，以确定加权信息。
应注意的是，声源可以(例如)是在环境中物理上存在的物理声源，例如，扬声器、乐器或个人说话。
然而，应注意的是，镜像源(mirror image source)也是声源。例如，由扬声器122发射的声波可被壁125反射，并且然后，声波似乎从与实际上发射声波的扬声器的位置不同的位置123发射。这种镜像源123也被视为声源。权重生成器110可被适配为根据与镜像源相关的到达方向信息或者根据关于一个、两个或更多个镜像源的位置信息来生成加权信息。
图2示出了根据实施方式的滤波器100以及多个麦克风111、112、113、…、11n。在图2的实施方式中，滤波器100进一步包括滤波器组101。而且，在图2的实施方式中，权重生成器110包括信息计算模块102、权重计算模块103以及传递函数选择模块104。
在时频域中进行处理，分别地，k表示频率指数，并且n表示时间指数。将M个时域麦克风信号x1...M(t)从麦克风111、112、13、…、11n中输入设备(滤波器100)中，通过滤波器组101将这些时域麦克风信号转换成时频域。由以下向量给出转换的麦克风信号：
x(k,n)＝[X1(k,n)X2(k,n)...XM(k,n)]T。
滤波器100输出期望信号Y(k,n)(音频输出信号)。音频输出信号(期望信号)Y(k,n)可例如表示用于单声道再现的增强信号、用于双耳声音再现的耳机信号或者用于具有任意的扬声器设置的空间声音再现的扬声器信号。
由输出信号生成器120生成期望信号Y(k,n)，例如，通过例如利用以下公式基于瞬时复杂权重w(k,n)＝[W1(k,n)W2(k,n)…WM(k,n)]T进行M个麦克风信号x(k,n)的线性组合：
Y(k，n)＝wH(k，n)x(k，n) (1)
通过权重计算模块103确定权重w(k,n)。对于每个k和每个n，重新确定确定w(k,n)。换言之，对于每个时频段(k,n)，进行权重w(k,n)的确定。更具体而言，例如，基于瞬时参数信息(IPI)并且基于相应的期望传递函数G(k,n)，计算权重w(k,n)。
信息计算模块102被配置为从麦克风信号x(k,n)计算IPIIPI描述在给定的时频瞬间(k,n)的信号以及包含在麦克风信号x(k,n)内的噪声分量的特定特征。
图3示出了根据实施方式的权重生成器110。权重生成器110包括信息计算模块102、权重计算模块103以及传递函数选择模块104。
如在图3中的实例中所示，IPI主要包括例如由DOA估计模块201计算的一个或方向性声音分量(例如，平面波)的瞬时到达方向(DOA)。
如下面所解释的，DOA信息可通过空间频率(例如，通过 )、通过相移(例如，通过)、通过麦克风之间的时间延迟、通过传播向量(例如，由)、通过双耳间强度差(ILD)或者由两耳时差(ITD)表示为角度(例如，通过[方位角仰角θ(k,n)])。
而且，IPI可例如包括其他信息，例如，信号或噪声分量的二阶统计量(SOS)。
在实施方式中，权重生成器110被适配为根据关于两个或更多个输入麦克风信号的信号或噪声分量的统计信息，为多个时频段(k,n)中的每个生成加权信息。例如，这种统计信息是在此处提及的二阶统计量。例如，统计信息可以是噪声分量的功率、信号扩散(signal-to-diffuse)信息、信号噪声(signal-to-noise)信息、扩散噪声(diffuse-to-noise)信息、信号分量的功率、扩散分量的功率、或者两个或更多个输入麦克风信号的信号分量或噪声分量的功率谱密度矩阵。
二阶统计量可以由统计计算模块205计算。该二阶统计量信息可例如包括固定噪声分量(例如，自噪声)的功率、非固定噪声分量(例如，扩散噪声)的功率、信号扩散比(SDR)、信噪比(SNR)或扩散噪声比(DNR)。该信息允许根据特定的优化标准计算最佳权重w(k,n)。
“固定噪声分量”/“缓慢变化的噪声分量”例如是具有随着时间不变化或缓慢变化的统计特性的噪声分量。
“非固定噪声分量”例如是具有随着时间快速变化的统计特性的噪声分量。
在实施方式中，权重生成器110被适配为根据表示关于两个或更多个输入麦克风信号的第一噪声分量的信息的第一噪声信息并且根据表示关于两个或更多个输入麦克风信号的第二噪声分量的信息的第二噪声信息，为多个时频段(k,n)中的每个生成加权信息。
例如，第一噪声分量可以是非固定噪声分量，并且第一噪声信息可以是关于非固定噪声分量的信息。
例如，第二噪声分量可以是固定噪声分量/缓慢变化的噪声分量，并且第二噪声信息可以是关于固定的/缓慢变化的噪声分量的信息。
在实施方式中，权重生成器110被配置为通过利用例如预定义的统计信息(例如，由非固定噪声分量引起的关于两个或更多个输入麦克风信号之间的空间相干性的信息)生成第一噪声信息(例如，关于非固定的/非缓慢变化的噪声分量的信息)，并且其中，权重生成器110被配置为在不利用统计信息的情况下生成第二噪声信息(例如，关于固定的/缓慢变化的噪声分量的信息)。
关于快速变化的噪声分量，输入麦克风信号不能单独提供足够的信息来确定关于这种噪声分量的信息。例如，此外，需要统计信息来确定关于快速变化的噪声分量的信息。
然而，关于不变化或不快速变化的噪声分量，不需要统计信息来确定关于这些噪声分量的信息。相反，评估麦克风信号是足够的。
应注意的是，可利用估计的DOA信息计算统计信息，如在图3中所示。应进一步注意的是，还可在外部提供IPI。例如，声音(分别是声源的位置)的DOA可由摄影机以及脸部识别算法确定，假设人类谈话者形成声音场景。
传递函数选择模块104被配置为提供传递函数G(k,n)。图2和图3的(可能复杂的)传递函数G(k,n)描述了系统的期望响应，给定(例如，当前参数)IPI例如，G(k,n)可以描述用于在单声道再现中的信号增强的期望空间麦克风的任意拾取图案(pick-uppattern)、用于扬声器再现的DOA相关的扬声器增益、或者双耳再现的头相关传递函数(HRTF)。
应注意的是，通常，记录的声音场景的统计信息随着时间和频率快速地变化。结果，IPI以及相应的最佳加权w(k,n)仅仅对于特定的时频指数有效，因此，给每个k和n重新计算。因此，系统可立即适配于当前的记录情况。
应进一步注意的是，M个输入麦克风可形成单个麦克风阵列，或者可分布为在不同的位置形成多个阵列。而且，IPI可包括位置信息，而非DOA信息，例如，声源在三维空间中的位置。借此，可定义空间滤波器，这些滤波器不仅根据需要对特定的方向滤波，而且对记录场景的三维空间区域进行滤波。
在声源的位置信息是可用的时，对DOA提供的所有解释同样适用。例如，位置信息可由DOA(角度)和距离表示。在使用这种位置表示时，可从位置信息中立即获得DOA。或者，位置信息可例如由x、y、z坐标描述。然后，基于声源的位置信息并且基于记录各自的输入麦克风信号的麦克风的位置，可容易计算DOA。
在下文中，描述进一步实施方式。
一些实施方式允许通过去混响和噪声降低进行在空间上的选择性声音记录。在该背景下，提供实施方式，应用空间滤波，用于在源提取、去混响以及噪声降低方面的信号增强。这种实施方式的目的在于，计算与具有任意拾取图案的方向麦克风的输出对应的信号Y(k,n)。这意味着方向性声音(例如，单个平面波)根据其DOA根据需要衰减或保存，同时抑制扩散声音或麦克风自噪声。根据实施方式，所提供的空间滤波器尤其结合使现有技术空间滤波器的优点，提供在具有高DNR的情况下的高方向性指数(DI)以及高白噪声增益(WNG)。根据一些实施方式，空间滤波器可仅仅受到线性约束，这允许快速计算权重。例如，图2和图3的传递函数G(k,n)可例如表示方向性麦克风的期望拾取图案。
在下文中，提供该问题的公式表示。然后，提供权重计算模块103和IPI计算模块102的实施方式，以利用去混响和噪声降低来进行在空间上的选择性声音记录。而且，描述了相应TF选择模块104的实施方式。
首先，提供该问题的公式表示。考虑位于d1...M处的M个全方向性麦克风的阵列。对于每个(k,n)，假设声场由在各向同性以及在空间上均匀的扩散声场中传播的L<M个平面波(方向声音)构成。麦克风信号x(k,n)可写为：
其中，xl(k,n)＝[Xl(k,n,d1)...Xl(k,n,dM)]T包括与第l个平面波的声压成比例的麦克风信号，xd(k,n)是所测量的非固定噪声(例如，扩散噪声)，并且xn(k,n)是固定噪声/缓慢变化的噪声(例如，麦克风自噪声)。
假设在公式(2)中的三个分量彼此不相关，则麦克风信号的功率谱密度(PSD)矩阵可由以下公式描述：
其中，Φd(k，n)＝φd(k，n)Γd(k) (4)
其中，Φn(k,n)是固定噪声/缓慢变化的噪声的PSD矩阵，并且φd(k,n)是非固定噪声的预期功率，该功率可随着时间和频率快速地变化。由γij(k)表示的相干矩阵Γd(k)的第ij个元素是由非固定噪声引起的在麦克风i和j之间的相干性。例如，对于球面各向同性扩散场，γij(k)＝sinc(κrij)[20]，其中，波数k和rij＝||dj-di||。相干矩阵Γd(k)的第ij个元素是由固定噪声/缓慢变化的噪声引起的在麦克风i和j之间的相干性。对于麦克风自噪声，Φn(k,n)＝φn(k,n)I，其中，I是单位矩阵并且φn(k,n)是自噪声的期望功率。
在(2)中的方向性声音xl(k,n)可写为：
其中，是第l个平面波的DOA的方位角(表示阵列垂射)并且是传播向量。 的第i个元素
描述了从第一到第i个麦克风的第l个平面波的相移。应注意的是，ri＝||di-d1||等于第一与第i个麦克风之间的距离。
角度通常称为空间频率。第l个波的DOA可由或者由表示。
如上所述，实施方式的目的在于，对麦克风信号x(k,n)进行滤波，使得从特定空间区域到达的方向性声音根据需要衰减或放大，同时抑制固定和非固定噪声。因此，期望的信号可表示为：
其中，是实值或复值任意的(例如，预定义的)方向性函数，该函数可以与频率相关。
图4涉及根据实施方式的具有两个任意的方向性函数和源位置的场景。尤其地，图4显示了两个实例方向性和的幅度。在使用(参照在图4中的实线)时，从中达到的方向声音衰减21dB，而其他方向的方向声音不衰减。原则上，任意方向性可被设计为偶函数，例如，(参照在图4中的虚线)。而且，可以设计为随着时间变化，例如，一旦定位，就提取移动的或出现的声源。
通过麦克风信号的线性组合，获得信号Y(k,n)的估计，例如，通过
其中，w(k,n)是长度M的复杂权重向量。在下面获得相应的最佳权重向量w(k,n)。在下文中，为了简单起见，省略权重w(k,n)对k和n的依赖性。
现在，描述在图2和图3中的权重计算模块103的两个实施方式。
从(5)和(7)中，遵循w(k,n)应满足线性约束：
而且，在滤波器的输出处的非固定和固定/缓慢变化的噪声功率应尽可能减小。
图5描述了应用空间滤波的本发明的实施方式。尤其地，图5示出了根据另一个实施方式的权重生成器110。再次，权重生成器110包括信息计算模块102、权重计算模块103以及传递函数选择模块104。
更尤其地，图5示出了线性约束的最小方差(LCMV)方法。在该实施方式中(见图5)，根据包括L个平面波的DOA的IPI I(k,n)以及固定和非固定的噪声的统计信息，计算权重w(k,n)。后面的信息可包括DNR、两个噪声分量的单独功率φn(k,n)以及φd(k,n)、或者两个噪声分量的PSD矩阵Φn以及Φd。
例如，Φd可被视为关于两个噪声分量的第一噪声分量的第一噪声信息，并且Φn可被视为关于两个噪声分量的第二噪声分量的第二噪声信息。
例如，权重生成器110可被配置为根据在一个或多个麦克风输入信号的至少一些第一噪声分量之间的一个或多个相干性，来确定第一噪声信息Φd。例如，权重生成器110可被配置为根据表示由两个或更多个输入麦克风信号的第一噪声分量引起的相干性的相干矩阵Γd(k)，来确定第一噪声信息，例如，通过应用公式Φd(k，n)＝φd(k，n)Γd(k)。
通过在滤波器的输出处最小化自噪声(固定噪声/缓慢变化的噪声)和扩散声音功率(非固定噪声)的总和，找出用于解决在(8)中的问题的权重w(k,n)，即，
使用(4)并且假设Φn(k,n)＝φn(k,n)I，优化问题可表示为：
其中，
是在麦克风处的时变输入DNR。给定约束条件(9)，对于(10)和(12)的解决方法是[21]。
＝C-1A[AHC-1A]-1g， (15)
其中，包括根据传播向量的L个平面波的DOA信息。相应的期望增益由以下公式给出：
下面描述Ψ(k,n)和其他需要的IPI的估计的实施方式。
其他实施方式基于参数多通道维纳滤波器。在这种实施方式中，如图6所示，IPI进一步包括关于信号统计的信息，例如，包括L个平面波(方向声音)的功率的信号PSD矩阵Φs(k,n)。而且，可选控制参数λ1...L(k,n)被视为控制L个平面波中的每个的信号失真的程度。
图6示出了实施利用参数多通道维纳滤波器的权重生成器110的空间滤波器的应用的实施方式。在此，权重生成器110包括信息计算模块102、权重计算模块103以及传递函数选择模块104。
通过参数多通道维纳滤波器方法，计算权重w(k,n)。维纳滤波器最小化在输出处的残余信号的功率，即，
被最小化的价值函数(cost function)C(k,n)可写为：
其中，Φs(k,n)＝E{xs(k,n)xs(k,n)H}包括方向声音PSD，并且xs(k,n)＝[X1(k,n,d1)X2(k,n,d1)….XL(k,n,d1)]包括与在参考麦克风处的L个平面波的声压成比例的信号。要注意的是，Φs(k,n)是对角矩阵，其中，对角元素diag{Φs(k,n)}＝[φ1(k,n)...φL(k,n)]T是到达的平面波的功率。为了控制引入的信号失真，可以包括对角矩阵Λ(k,n)，其包括时间和频率相关的控制参数diag{Λ}＝[λ1(k,n)λ2(k,n)...λL(k,n)]T，即，
CPW(k，n)＝[g-AH(k，n)w]HΛ(k，n)Φs(k，n)[g-AH(k，n)w]
wHΦu(k，n)w. (20)
考虑CPW(k,n)的在(17)中的最小化问题的解决方案是：
w＝[AHΛ(k，n)Φs(k，n)A+Φu]-1AΛ(k，n)Φs(k，n)g. (21)
这等同于
应注意的是，对于，Λ-1＝0，获得在(14)中的LCMV解决方案。对于Λ-1＝I，获得多通道维纳滤波器。对于其他值λ1...L(k,n)，可分别控制相应源信号的失真程度以及残余噪声抑制的程度。因此，通常根据可用的参数信息，限定即，
其中，f(·)是任意用户定义的函数。例如，可根据以下公式选择λ1...L(k,n)：
其中，φl(k,n)是第l个信号(第l个平面波)的功率，并且φu(k,n)＝φn(k,n)+φd(k,n)是不期望信号(固定噪声/缓慢变化的噪声加上非固定的噪声)的功率。借此，参数维纳滤波器取决于关于两个或更多个输入麦克风信号的信号分量的统计信息，因此，参数维纳滤波器进一步取决于关于两个或更多个输入麦克风信号的噪声分量的统计信息。
如果与噪声相比，源1强烈，则获得接近0的意味着获得LCMV解决方案(没有源信号的失真)。如果与源功率相比，噪声强烈，则获得接近1的意味着获得参数维纳滤波器(强烈抑制噪声)。
下面描述Φs(k,n)和Φu(k,n)的估计。
在下文中，描述瞬时参数估计模块102的实施方式。
在能够计算权重之前，需要估计不同的IPI。可通过众所周知的窄带DOA估计器(例如，ESPRIT[22]或根MUSIC[23])或其他现有技术估计器获得在模块201中计算的L个平面波。例如，这些算法可为到达阵列的一个或多个波提供(例如)方位角空间频率相移或传播向量将不进一步讨论DOA估计，因为DOA估计本身在本领域中是熟知的。
在下文中，描述扩散噪声比(DNR)估计。尤其地，描述输入DNRΨ(k,n)的估计，即，在图5中的模块202的实现。DNR估计利用在模块201中获得的DOA信息。为了估计Ψ(k,n)，可使用额外的空间滤波器，该滤波器去除L个平面波，以便仅仅捕捉扩散声音。例如，通过最大化阵列的WNG，找出这个空间滤波器的权重，即，
服从
约束条件(27)确保非零加权wΨ。传播向量与特定的方向对应，方向与L个平面波的DOA不同。在下文中，对于选择与所有相距最大距离的方向，即，
其中，给定权重wΨ，由以下公式给出额外空间滤波器的输出功率：
现在，可通过(13)和(29)计算输入DNR，即，
假设功率随着时间是恒定的或者缓慢地变化，在静音(silence)时，例如，可以估计所需要的麦克风自噪声φn(k,n)的预期功率。要注意的是，由于所选择的优化标准(45)，所以所提出的DNR估计器实际上不必提供最低估计方差，但是提供公正的结果。
在下文中，讨论非固定PSDφd(k,n)的估计，即，在图5中的模块(202)的另一个实现。可以使用以下公式，估计非固定噪声的功率(PSD)：
其中，在前一个段落中，定义了wΨ。应注意的是，在静音时(即，在没有信号和非固定噪声时)，可估计固定/缓慢变化的噪声PSD矩阵Φn(k,n)，即，
Φn(k，n)＝E{x(k，n)xH(k，n)}， (32)
其中，通过对静音帧(silent frame)n求平均来逼近期望值。可以通过现有技术方法检测静音帧。
在下文中，讨论不期望信号PSD矩阵(见模块203)的估计。
可通过以下公式：
Φu(k，n)＝φn(k，n)(Ψ(k，n)Γd(k)+Γn(k))， (33)
或者更通常通过以下公式：
Φu(k，n)＝φd(k，n)Γd(k)+Φn(k，n)， (34)，
获得不期望信号(固定/缓慢变化的噪声加上非固定的噪声)Φu(k,n)的PSD矩阵。
其中，Γd(k)和Γn(k)可用作先验信息(见上文)。如上所解释的，可计算DNRΨ(k,n)、固定/缓慢变化的噪声功率φn(k,n)以及其他需要的量。因此，Φu(k,n)估计利用由模块201获得的DOA信息。
在下文中，描述信号PSD矩阵(见模块204)的估计。
可通过以下公式，计算用于计算Φs(k,n)所需要的到达平面波的功率φ1...L(k,n)：
其中，权重wl抑制所有到达的平面波，除了第l个波，即，
例如，
经受(36)。Φs(k,n)估计利用在模块(201)中获得的DOA信息。如在前一段落中所解释的，可计算不期望信号Φu(k,n)的所需要的PSD矩阵。
现在，描述根据实施方式的传递函数选择模块104。
在本申请中，可以根据DOA信息找出相应的平面波1的增益不同的传递函数可用于该系统，例如，作为用户定义的先验信息。还可基于图像的分析，例如，使用所检测的脸部的位置，计算增益。在图4中描述两个实例。这些传递函数与方向麦克风的期望的拾取图案对应。可提供传递函数例如，作为查找表，即，对于估计的从查找表中选择相应的增益 要注意的是，还可根据空间频率而非方位角即，根据G(k,μ)，而非定义传递函数。还可基于源位置信息而非DOA信息，计算增益
现在，提供实验结果。以下模拟结果展示了上述实施方式的可适用性。比较所提出的系统和现有技术系统，下面将进行解释。然后，讨论实验设置并且提供结果。
首先，考虑现有空间滤波器。
虽然在静音期间，可估计PSDφn(k,n)，但是通常假设φd(k,n)未知并且难以察觉。因此，考虑可在无这种知识的情况下计算的两个现有空间滤波器。
第一空间滤波器已知为延时求和波束成形器，并且最小化在滤波器的输出处的自噪声功率[即，最大化WNG][1]。然后，通过以下公式，获得使经受(9)的在(7)与(8)之间的均方误差(MSE)最小化的优化权重向量：
存在针对(38)的封闭式解决方案[1]，该解决方案允许快速计算wn。应注意的是，该滤波器不必提供最大的DI。
第二空间滤波器已知为稳健的超方向性(SD)波束成形器，并且通过在WNG上的下限，最小化在滤波器的输出处的扩散声音功率[即，最大化DI][24]。在WNG上的下限增大对在传播向量内的误差的稳健性，并且限制自噪声的放大[24]。然后，通过以下公式，获得最小化经受(9)的在(7)与(8)之间的MSE并且满足在WNG上的下限的优化加权向量：
并且经受二次约束wH w<β。参数β-1定义最小WNG并且确定滤波器的可实现的DI。实际上，通常难以在低SNR情况下的充分WNG与在高SNR情况下的足够高的DI之间找出最佳权衡。而且，由于二次约束，所以解决(39)，造成非凸优化问题，这需要时间来解决。由于时变约束(9)造成需要给每个k和n重新计算复杂加权向量，所以这尤其是个问题。
现在，考虑实验设置。假设在(2)中的模型中的L＝2个平面波以及具有麦克风间的间距3cm的M＝4个麦克风的均匀线性阵列(ULA)，使用源图像方法[25、26]，模拟鞋盒空间(7.0×5.4×2.4m3、RT60≈380ms)，在以及时分别具有两个声源(距离1.75m，比较图4)。信号包括0.6s静音，其后跟随模糊词(double talk)。将白高斯噪声加到麦克风信号中，产生26dB的分段信噪比(SSNR)。声音以16kHz取样并且使用512点STFT转换成时频域，具有50％重叠。
考虑图4的方向性函数即，应在不失真的情况下提取源A，同时源B的功率衰减21dB。考虑上面这两个空间滤波器以及所提供的空间滤波器。对于稳健的SD波束成形器(39)，最小的WNG设为-12dB。对于提供的空间滤波器(12)，如上面所解释的，估计DNRΨ(k,n)。在开始时，从静音信号部分计算自噪声功率φn(k,n)。在(3)中的预期由超过τ＝50ms的递归平均时间逼近。
在下文中，考虑非时变方向约束。
为了这个模拟，假设关于这两个源位置和的先验知识。在所有处理步骤中，使用和因此，在(9)和(26)中的方向约束不随着时间变化。
图7示出了真实的和估计的DNRΨ(k,n)。这两个标记的区域分别表示信号的静音和活动部分。尤其地，图7描述了根据时间和频率的真实的和估计的DNRΨ(k,n)。由于混响环境，所以在语音活动期间，获得较高的DNR。由于结合的时间平均化处理，所以在图7(b)中的估计的DNR处理有限的时间分辨率。然而，Ψ(k,n)估计足够精确，如以下结果所示。
图8(a)描述了wn和wd(均与信号无关)的以及所提出的空间滤波器wnd(与信号相关)的平均DI。对于所提出的空间滤波器，显示了信号的静音部分以及在语音活动期间[在图7(b)中标记这两个信号部分]的DI。在静音期间，所提出的空间滤波器(虚线wnd)提供与wn相同的低DI。在语音活动期间(实线wnd)，所获得的DI与稳健的SD波束成形器(wd)一样高。图8(b)显示了相应的WNG。在静音期间，所提出的空间滤波器(虚线wnd)实现了高WNG，而在信号活动期间，WNG较低。
图8：比较空间滤波器的DI和WNG。对于wd，最小WNG设为-12dB，以使空间滤波器稳健，抵抗麦克风自噪声。
通常，图8示出了所提出的空间滤波器使这两个现有空间滤波器的优点相结合：在静音部分期间，提供最大WNG，导致最小自噪声放大，即，高稳健性。
在通常掩蔽自噪声的信号活动和高混响期间，提供高DI(以低WNG为代价)，导致扩散声音的最佳降低。在这种情况下，甚至相当小的WNG可容忍。
要注意的是，对于更高的频率(f>5kHz)，由于在(39)和(12)中的相关矩阵Γd(k)几乎等于单位矩阵，所以所有空间滤波器几乎相同地运行。
在下文中，考虑瞬时方向约束。
对于该模拟，假设没有关于和的可用先验信息。通过ESPRIT估计和因此，约束条件(9)随着时间变化。仅仅对于稳健的SD波束成形器(wd)，使用与固定观测方向对应的单个和非时变约束条件(9)。这个波束成形器用作参考。
图9描述了估计的和所产生的增益尤其地，图9示出了估计的和所产生的增益|如果DOA在图4中位于空间窗口内部(实线)，则到达的平面波不衰减。否则，波的功率衰减21dB。
表1示出了所有空间滤波器(*未处理的)的性能。在括号内的值表示非时变方向约束条件，不在括号内的值表示瞬时方向约束条件。在计算SIR、SRR以及SSNR之前，信号进行A加权。
表1
尤其地，表1总结了空间滤波器在滤波器的输出处的在信号干扰比(SIR)、信号混响比(SRR)以及SSNR方面的总体性能。在SIR和SRR(源分离、去混响)方面，所提出的方法(wnd)和稳健的SD波束成形器(wd)提供最高性能。然而，所提出的wnd的SSNR比wd的SSNR高6dB，这表示可以清楚听见的优点。使用wn，获得在SSNR方面的最佳性能。在PESQ方面，wnd和wd胜过wn。使用瞬时方向约束条件，而非时变约束条件(在括号内的值)主要降低可实现的SIR，但是在改变源位置的情况下，提供快速适应。应注意的是，每个时间帧的所有需要的复杂加权的计算时间对于wd大于80s(CVX工具箱[27,28])并且对于所提出的方法小于0.08s(MATLAB R2012b、MacBook Pro 2008)。
在下文中，描述空间声音再现的实施方式。实施方式的目的在于，(例如)通过麦克风阵列捕捉声音场景，并且通过任意的声音再现系统(例如，5.1扬声器设置、耳机再现)再现空间声音，以便再创建原始空间印象。假设声音再现系统包括N个通道，即，计算N个输出信号Y(k,n)。
首先，提供问题公式表示。考虑信号模型(见以上公式(2))，并且公式表示相似的问题。固定/缓慢变化的噪声与不期望的麦克风自噪声对应，而非固定的噪声与期望的扩散声音对应。由于再现记录场景的原始空间印象至关重要，所以在本申请中期望扩散声音。
在下文中，应在没有失真的情况下从相应的实现再现方向声音Xl(k,n,d1)。而且，应利用来自所有方向的正确能量再现扩散声音，同时抑制麦克风自噪声。因此，在(7)中的期望的信号Y(k,n)现在表示为：
其中，Yi(k,n)是声音再现系统的第i个信道的信号(i＝{1,...,N})，Xd,i(k,n,d)是在从扬声器i中再现的任意点(例如，在第一麦克风d1)处测量的扩散声音，并且Gd(k,n)是扩散声音的增益函数，以在再现期间确保扩散声音的正确功率(通常)。理想地，信号Xd,i(k,n)具有正确的扩散声音功率，并且在信道i上相互不相关，即，
方向声音分量的传递函数与和DOA相关的扬声器增益函数对应。在图10中描述用于立体扬声器再现的情况的实例。如果波1从中到达，那么G1＝1并且G2＝0。这表示仅从再现系统的信道i＝1(左信道)中再现这个方向声音。对于具有 即，通过这两个扬声器的相等功率，再现方向声音。或者，如果期望双耳再现，那么可以与HRTF对应。
如上所述，基于复杂权重w(k,n)，通过麦克风信号的线性组合，估计信号Yi(k,n)，即，
经受特定的约束条件。在接下来的小部分中解释加权w(k,n)的约束条件和计算。
在下文中，在本背景下，考虑根据相应的实施方式的权重计算模块103，提供图2的权重计算模块103的两个实施方式。从公式(5)和公式(40)中得出以下结论：wi(k,n)应满足线性约束：
而且，应保持扩散声音功率。因此，wi(k,n)可以满足二次约束：
而且，在滤波器的输出处的自噪声功率应最小化。因此，最佳权重可计算为：
经受公式(43)和公式(44)。这导致(例如)通过熟知的数值方法[29]可以解决的凸优化问题。
对于瞬时参数估计模块102，根据相应的实施方式，通过熟知的窄带DOA估计器(例如，ESPRIT[22]或根MUSIC[23])或其他现有技术估计器，可以获得L个平面波的
现在，考虑根据相应的实施方式的传递函数选择模块104。在本申请中，根据DOA信息为相应的方向声音1找出信道i的增益 用于不同的和信道i的传递函数可用于系统，例如，作为用户定义的先验信息。还可基于图像的分析，例如，使用所检测的脸部的位置，计算增益。
通常提供传递函数作为查找表，即，为估计的从查找表中选择相应的增益要注意的是，还可以根据空间频率而非方位角即，根据Gi(k,μ)，而非 定义传递函数。要注意的是，传递函数还可与能够进行双耳声再现的HRTF对应。在这种情况下，通常是复值。要注意的是，还可以根据源位置信息，而非DOA信息，计算增益或传递函数。
在图10中描述了立体扬声器再现的一个实例。尤其地，图10示出了立体再现的增益函数。
虽然在设备的背景下描述了一些方面，但是显然，这些方面还表示相应方法的描述，其中，方框或装置与方法步骤或方法步骤的特征对应。相似地，在方法步骤的背景下描述的方面还表示相应设备的相应方框或项目或特征的描述。
本发明的分解信号可存储在数字储存介质上或者可在传输介质上传输，例如，无线传输介质或有线传输介质，诸如，互联网。
根据某些实施要求，本发明的实施方式可以在硬件中或者在软件中实施。可使用数字储存介质(例如，软盘、DVD、CD、ROM、PROM、EPROM、EEPROM或FLASH执行实施)，在所述数字储存介质上存储了电子可读的控制信号，其与可编程计算机系统协作(或者能够协作)，以便执行相应的方法。
根据本发明的一些实施方式包括非瞬时性数据载体，其具有电子可读控制信号，这些信号能够与可编程计算机系统协和，以便执行在本文中描述的一种方法。
通常，本发明的实施方式可实施为具有程序代码的计算机程序产品，当计算机程序产品在计算机上运行时，该程序代码可操作为执行一种方法。例如，该程序代码可存储在机器可读载体上。
其他实施方式包括存储在机器可读载体上的用于执行在本文中描述的一种方法的计算机程序。
换言之，因此，本发明的方法的实施方式是具有程序代码的计算机程序产品，当计算机程序在计算机上运行时，所述程序代码用于执行在本文中描述的一种方法。
因此，本发明的方法的进一步实施方式是数据载体(或数字存储介质或计算机可读介质)，其包括存记录在其上的用于执行在本文中描述的一种方法的计算机程序
因此，本发明的方法的进一步实施方式是数据流或信号序列，其表示用于执行在本文中描述的一种方法的计算机程序。例如，数据流或信号序列可被配置为通过数据通信连接(例如，通过互联网)传输。
进一步实施方式包括处理器件，例如，计算机或可编程逻辑装置，其被配置为或者适配为执行在本文中描述的一种方法。
进一步实施方式包括具有安装在其上的计算机程序的计算机，所述计算机程序用于执行在本文中描述的一种方法。
在一些实施方式中，可编程逻辑装置(例如，现场可编程门阵列)可用于执行在本文中描述的方法的一些或所有功能。在一些实施方式中，现场可编程门阵列可与微处理器协作，以便执行在本文中描述的一种方法。通常，这些方法优选地由任何硬件设备执行。
上述实施方式仅仅说明本发明的原理。要理解的是，在本文中描述的设置和细节的修改和变化对于本领域的技术人员将是显而易见的。因此，其目的在于，仅受到待审专利权利要求的范围的限制，而不受到通过在本文中的实施方式的描述和解释呈现的具体细节的限制。
参考文献
[1]J.Benesty,J.Chen,and Y.Huang,Microphone Array SignalProcessing.Berlin,Germany:Springer-Verlag,2008.
[2]S.Doclo,S.Gannot,M.Moonen,and A.Spriet,“Acoustic beamforming forhearing aid applications,”in Handbook on Array Processing and SensorNetworks,S.Haykin and K.Ray Liu,Eds.Wiley,2008,ch.9.
[3]S.Gannot and I.Cohen,“Adaptive beamforming and postfiltering,”inSpringer Handbook of Speech Processing,J.Benesty,M.M.Sondhi,and Y.Huang,Eds.Springer-Verlag,2008,ch.47.
[4]J.Benesty,J.Chen,and E.A.P.Habets,Speech Enhancement in the STFTDomain,ser.SpringerBriefs in Electrical and Computer Engineering.Springer-Verlag,2011.
[5]I.Tashev,M.Seltzer,and A.Acero,“Microphone array for headset withspatial noise suppressor,”in Proc.Ninth International Workshop on Acoustic,Echo and Noise Control(IWAENC),Eindhoven,The Netherlands,2005.
[6]M.Kallinger,G.Del Galdo,F.Kuech,D.Mahne,and R.Schultz-Amling,“Spatial filtering using directional audio coding parameters,”in Proc.IEEEIntl.Conf.on Acoustics,Speech and Signal Processing(ICASSP),Apr.2009,pp.217–220.
[7]M.Kallinger,G.D.Galdo,F.Kuech,and O.Thiergart,“Dereverberation inthe spatial audio coding domain,”in Audio Engineering Society Convention 130,London UK,May 2011.
[8]G.Del Galdo,O.Thiergart,T.Weller,and E.A.P.Habets,“Generatingvirtual microphone signals using geometrical information gathered bydistributed arrays,”in Proc.Hands-Free Speech Communication and MicrophoneArrays(HSCMA),Edinburgh,United Kingdom,May 2011.
[9]S.Nordholm,I.Claesson,and B.Bengtsson,“Adaptive array noisesuppression of handsfree speaker input in cars,”IEEE Trans.Veh.Technol.,vol.42,no.4,pp.514–518,Nov.1993.
[10]O.Hoshuyama,A.Sugiyama,and A.Hirano,“A robust adaptive beamformerfor microphone arrays with a blocking matrix using constrained adaptivefilters,”IEEE Trans.Signal Process.,vol.47,no.10,pp.2677–2684,Oct.1999.
[11]S.Gannot,D.Burshtein,and E.Weinstein,“Signal enhancement usingbeamforming and nonstationarity with applications to speech,”IEEETrans.Signal Process.,vol.49,no.8,pp.1614–1626,Aug.2001.
[12]W.Herbordt and W.Kellermann,“Adaptive beamforming for audiosignal acquisition,”in Adaptive Signal Processing:Applications to real-worldproblems,ser.Signals and Communication Technology,J.Benesty and Y.Huang,Eds.Berlin,Germany:Springer-Verlag,2003,ch.6,pp.155–194.
[13]R.Talmon,I.Cohen,and S.Gannot,“Convolutive transfer functiongeneralized sidelobe canceler,”IEEE Trans.Audio,Speech,Lang.Process.,vol.17,no.7,pp.1420–1434,Sep.2009.
[14]A.Krueger,E.Warsitz,and R.Haeb-Umbach,“Speech enhancement with aGSC-like structure employing eigenvector-based transfer function ratiosestimation,”IEEE Trans.Audio,Speech,Lang.Process.,vol.19,no.1,pp.206–219,Jan.2011.
[15]E.A.P.Habets and J.Benesty,“Joint dereverberation and noisereduction using a two-stage beamforming approach,”in Proc.Hands-Free SpeechCommunication and Microphone Arrays(HSCMA),2011,pp.191–195.
[16]M.Taseska and E.A.P.Habets,“MMSE-based blind source extraction indiffuse noise fields using a complex coherence-based a priori SAP estimator,”in Proc.Intl.Workshop Acoust.Signal Enhancement(IWAENC),Sep.2012.
[17]G.Reuven,S.Gannot,and I.Cohen,“Dual source transfer-functiongeneralized sidelobe canceller,”IEEE Trans.Speech Audio Process.,vol.16,no.4,pp.711–727,May 2008.
[18]S.Markovich,S.Gannot,and I.Cohen,“Multichannel eigenspacebeamforming in a reverberant noisy environment with multiple interferingspeech signals,”IEEE Trans.Audio,Speech,Lang.Process.,vol.17,no.6,pp.1071–1086,Aug.2009.
[19]O.Thiergart and E.A.P.Habets,“Sound field model violations inparametric spatial sound processing,”in Proc.Intl.Workshop Acoust.SignalEnhancement(IWAENC),Sep.2012.
[20]R.K.Cook,R.V.Waterhouse,R.D.Berendt,S.Edelman,and M.C.ThompsonJr.,“Measurement of correlation coefficients in reverberant sound fields,”TheJournal of the Acoustical Society of America,vol.27,no.6,pp.1072–1077,1955.
[21]O.L.Frost,III,“An algorithm for linearly constrained adaptivearray processing,”Proc.IEEE,vol.60,no.8,pp.926–935,Aug.1972.
[22]R.Roy and T.Kailath,“ESPRIT-estimation of signal parameters viarotational invariance techniques,”Acoustics,Speech and Signal Processing,IEEETransactions on,vol.37,no.7,pp.984–995,July 1989.
[23]B.Rao and K.Hari,“Performance analysis of root-music*,”inSignals,Systems and Computers,1988.Twenty-Second Asilomar Conference on,vol.2,1988,pp.578–582.
[24]H.Cox,R.M.Zeskind,and M.M.Owen,“Robust adaptive beamforming,”IEEETrans.Acoust.,Speech,Signal Process.,vol.35,no.10,pp.1365–1376,Oct.1987.
[25]J.B.Allen and D.A.Berkley,“Image method for efficientlysimulating small-room acoustics,”J.Acoust.Soc.Am.,vol.65,no.4,pp.943–950,Apr.1979.
[26]E.A.P.Habets.(2008,May)Room impulse response(RIR)generator.[Online].Available:http://home.tiscali.nl/ehabets/rirgenerator.html；还见:http://web.archive.org/web/20120730003147/http://home.tiscali.nl/ehabets/rir_generator.html
[27]I.CVX Research,“CVX:Matlab software for disciplined convexprogramming,version 2.0beta,”http://cvxr.com/cvx,September 2012.
[28]M.Grant and S.Boyd,“Graph implementations for nonsmooth convexprograms,”in Recent Advances in Learning and Control,ser.Lecture Notes inControl and Information Sciences,V.Blondel,S.Boyd,and H.Kimura,Eds.Springer-Verlag Limited,2008,pp.95–110.
[29]H.L.Van Trees,Detection,Estimation,and Modulation Theory:Part IV:Optimum Array Processing.John Wiley&Sons,April 2002,vol.1.

Claims (15)

1.一种滤波器(100)，用于基于两个或更多个输入麦克风信号生成音频输出信号，所述音频输出信号包括多个音频输出信号样本，其中，在时频域中表示所述音频输出信号以及所述两个或更多个输入麦克风信号，其中，所述多个音频输出信号样本中的每个被分配给多个时频段中的一个时频段，并且其中，所述滤波器(100)包括：
权重生成器(110)，被适配成：对于所述多个时频段中的每个，接收一个或多个声源的一个或多个声音分量的到达方向信息或所述一个或多个声源的位置信息，并且所述权重生成器(110)被适配为根据所述时频段的所述一个或多个声源的所述一个或多个声音分量的所述到达方向信息或者根据所述时频段的所述一个或多个声源的所述位置信息，为所述多个时频段中的每个生成加权信息；其中，所述权重生成器(110)被适配为根据第一噪声信息并根据第二噪声信息为所述多个时频段中的每个生成所述加权信息，所述第一噪声信息指示关于所述两个或更多个输入麦克风信号的第一噪声分量的第一相干矩阵的信息，所述第二噪声信息指示关于所述两个或更多个输入麦克风信号的第二噪声分量的第二相干矩阵的信息；以及
输出信号生成器(120)，用于通过根据所述时频段的所述加权信息并且根据所述两个或更多个输入麦克风信号中的每个的被分配给所述时频段的音频输入样本，为所述多个时频段中的每个生成分配给所述时频段的所述多个音频输出信号样本中的一个，来生成所述音频输出信号。
2.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被配置为通过利用统计信息来生成所述第一噪声信息，并且其中，所述权重生成器(110)被配置为在不利用所述统计信息的情况下生成所述第二噪声信息，其中，预定义所述统计信息。
3.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被适配为根据以下公式为所述多个时频段中的每个生成所述加权信息：
其中，Φu＝Φd+Φn，
其中，Φd是所述两个或更多个输入麦克风信号的所述第一噪声分量的第一功率谱密度矩阵，
其中，Φn是所述两个或更多个输入麦克风信号的所述第二噪声分量的第二功率谱密度矩阵，
其中，A表示所述到达方向信息，
其中，wnd是表示所述加权信息的向量，
其中，g表示增益，
其中，是根据所述到达方向信息的第一实值或复值的预定义的方向性函数，并且
其中，是根据所述到达方向信息的其他实值或复值的预定义的方向性函数，
其中，k表示频率，并且其中，n表示时间。
4.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被配置为根据所述两个或更多个输入麦克风信号的所述第一噪声分量的至少一些之间的一个或多个相干性，来确定所述第一噪声信息，其中，预定义所述一个或多个相干性。
5.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被配置为根据相干矩阵Γd(k)来确定所述第一噪声信息，所述相干矩阵Γd(k)指示由于所述两个或更多个输入麦克风信号的所述第一噪声分量引起的相干性，其中，预定义所述相干矩阵Γd(k)，其中，k表示频率。
6.根据权利要求5所述的滤波器(100)，其中，所述权重生成器(110)被配置为根据以下公式来确定所述第一噪声信息：
Φd(k，n)＝φd(k，n)Γd(k)，
其中，Γd(k)是所述相干矩阵，其中，预定义所述相干矩阵，
其中，Φd(k,n)是所述第一噪声信息，并且
其中，φd(k,n)是所述两个或更多个输入麦克风信号的所述第一噪声分量的预期功率，
其中，k表示频率，并且其中，n表示时间。
7.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被配置为根据所述第二噪声信息并且根据所述到达方向信息，来确定所述第一噪声信息。
8.根据权利要求1所述的滤波器(100)，
其中，所述权重生成器(110)被配置为生成所述加权信息作为第一加权信息wΨ，并且
其中，所述权重生成器(110)被配置为通过确定第二加权信息来生成所述第一加权信息，
其中，所述权重生成器(110)被配置为通过应用以下公式来生成所述第一加权信息wΨ
以便满足以下公式：
其中，表示所述到达方向信息，
其中，表示传播向量，并且
其中，w表示所述第二加权信息，
其中，k表示频率，并且其中，n表示时间。
9.根据权利要求8所述的滤波器(100)，其中，所述权重生成器(110)被配置为根据所述第二加权信息并且根据所述两个或更多个输入麦克风信号来生成扩散噪声信息或扩散分量的功率，以确定所述第一加权信息。
10.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被配置为通过应用参数维纳滤波器来确定所述加权信息，其中，所述参数维纳滤波器取决于关于所述两个或更多个输入麦克风信号的信号分量的统计信息，并且其中，所述参数维纳滤波器取决于关于所述两个或更多个输入麦克风信号的噪声分量的统计信息。
11.根据权利要求1所述的滤波器(100)，其中，所述权重生成器(110)被配置为根据指示一个或多个平面波的到达方向的所述到达方向信息，来确定所述加权信息。
12.根据权利要求1所述的滤波器(100)，
其中，所述权重生成器(110)包括用于提供预定义的传递函数的传递函数选择模块(104)，并且
其中，所述权重生成器(110)被配置为根据所述到达方向信息并且根据所述预定义的传递函数，来生成所述加权信息。
13.根据权利要求12所述的滤波器(100)，其中，所述传递函数选择模块(104)被配置为提供所述预定义的传递函数，以便所述预定义的传递函数根据所述到达方向信息指示任意拾取图案，以便所述预定义的传递函数根据所述到达方向信息指示扬声器增益，或者以便所述预定义的传递函数根据所述到达方向信息指示头相关传递函数。
14.一种用于基于两个或更多个输入麦克风信号生成音频输出信号的方法，所述音频输出信号包括多个音频输出信号样本，其中，在时频域中表示所述音频输出信号以及所述两个或更多个输入麦克风信号，其中，所述多个音频输出信号样本中的每个被分配给多个时频段中的一个时频段，并且其中，所述方法包括：
为所述多个时频段中的每个，接收一个或多个声源的一个或多个声音分量的到达方向信息或所述一个或多个声源的位置信息，
根据所述时频段的所述一个或多个声源的所述一个或多个声音分量的所述到达方向信息或者根据所述时频段的所述一个或多个声源的所述位置信息，为所述多个时频段中的每个生成加权信息，其中，根据第一噪声信息并根据第二噪声信息为所述多个时频段中的每个生成所述加权信息，所述第一噪声信息指示关于所述两个或更多个输入麦克风信号的第一噪声分量的第一相干矩阵的信息，所述第二噪声信息指示关于所述两个或更多个输入麦克风信号的第二噪声分量的第二相干矩阵的信息；并且
通过根据所述时频段的所述加权信息并且根据所述两个或更多个输入麦克风信号中的每个的被分配给所述时频段的音频输入样本，为所述多个时频段中的每个生成分配给所述时频段的所述多个音频输出信号样本中的一个，来生成所述音频输出信号。
15.一种计算机可读存储介质，存储有计算机可执行代码，用于当在计算机或信号处理器上执行时实施根据权利要求14所述的方法。
CN201380073406.6A 2012-12-21 2013-11-25 使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法 Active CN105165026B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201261740866P 2012-12-21 2012-12-21
US61/740,866 2012-12-21
EP13169163.6A EP2747451A1 (en) 2012-12-21 2013-05-24 Filter and method for informed spatial filtering using multiple instantaneous direction-of-arrivial estimates
EP13169163.6 2013-05-24
PCT/EP2013/074650 WO2014095250A1 (en) 2012-12-21 2013-11-25 Filter and method for informed spatial filtering using multiple instantaneous direction-of-arrivial estimates

Publications (2)

Publication Number Publication Date
CN105165026A CN105165026A (zh) 2015-12-16
CN105165026B true CN105165026B (zh) 2019-08-13

Family

ID=48607016

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380073406.6A Active CN105165026B (zh) 2012-12-21 2013-11-25 使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法

Country Status (8)

Country Link
US (1) US10331396B2 (zh)
EP (2) EP2747451A1 (zh)
JP (1) JP6196320B2 (zh)
CN (1) CN105165026B (zh)
BR (1) BR112015014380B1 (zh)
ES (1) ES2612528T3 (zh)
RU (1) RU2641319C2 (zh)
WO (1) WO2014095250A1 (zh)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9788128B2 (en) * 2013-06-14 2017-10-10 Gn Hearing A/S Hearing instrument with off-line speech messages
GB2521649B (en) * 2013-12-27 2018-12-12 Nokia Technologies Oy Method, apparatus, computer program code and storage medium for processing audio signals
EP2975609A1 (en) 2014-07-15 2016-01-20 Ecole Polytechnique Federale De Lausanne (Epfl) Optimal acoustic rake receiver
US9949041B2 (en) † 2014-08-12 2018-04-17 Starkey Laboratories, Inc. Hearing assistance device with beamformer optimized using a priori spatial information
US10602266B2 (en) * 2014-10-10 2020-03-24 Sony Corporation Audio processing apparatus and method, and program
US10242690B2 (en) 2014-12-12 2019-03-26 Nuance Communications, Inc. System and method for speech enhancement using a coherent to diffuse sound ratio
EP3220386A1 (en) * 2016-03-18 2017-09-20 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus and method for harmonic-percussive-residual sound separation using a structure tensor on spectrograms
EP3430620B1 (en) 2016-03-18 2020-03-25 Fraunhofer Gesellschaft zur Förderung der Angewand Encoding by reconstructing phase information using a structure tensor on audio spectrograms
CN106028227B (zh) * 2016-07-08 2019-05-24 乐鑫信息科技（上海）股份有限公司 分布式麦克风阵列及其适用的声源定位系统
CN106060743A (zh) * 2016-08-03 2016-10-26 上海山景集成电路股份有限公司 麦克风、麦克风组合及麦克风信号处理方法
CN106569773A (zh) * 2016-10-31 2017-04-19 努比亚技术有限公司 一种终端和语音交互的处理方法
CN106782590B (zh) * 2016-12-14 2020-10-09 南京信息工程大学 基于混响环境下麦克风阵列波束形成方法
GB2563606A (en) 2017-06-20 2018-12-26 Nokia Technologies Oy Spatial audio processing
GB2563952A (en) * 2017-06-29 2019-01-02 Cirrus Logic Int Semiconductor Ltd Speaker identification
EP3460795A1 (en) * 2017-09-21 2019-03-27 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Signal processor and method for providing a processed audio signal reducing noise and reverberation
CN111201784B (zh) 2017-10-17 2021-09-07 惠普发展公司，有限责任合伙企业 通信系统、用于通信的方法和视频会议系统
US10679617B2 (en) * 2017-12-06 2020-06-09 Synaptics Incorporated Voice enhancement in audio signals through modified generalized eigenvalue beamformer
TWI690218B (zh) 2018-06-15 2020-04-01 瑞昱半導體股份有限公司 耳機
CN110636400B (zh) * 2018-06-25 2021-03-16 瑞昱半导体股份有限公司 耳机
US11335357B2 (en) * 2018-08-14 2022-05-17 Bose Corporation Playback enhancement in audio systems
CN109839612B (zh) * 2018-08-31 2022-03-01 大象声科（深圳）科技有限公司 基于时频掩蔽和深度神经网络的声源方向估计方法及装置
CN109286875B (zh) * 2018-09-29 2021-01-01 百度在线网络技术（北京）有限公司 用于定向拾音的方法、装置、电子设备和存储介质
JP7407580B2 (ja) 2018-12-06 2024-01-04 シナプティクス インコーポレイテッド システム、及び、方法
US20200184994A1 (en) * 2018-12-07 2020-06-11 Nuance Communications, Inc. System and method for acoustic localization of multiple sources using spatial pre-filtering
CN111025233B (zh) * 2019-11-13 2023-09-15 阿里巴巴集团控股有限公司 一种声源方向定位方法和装置、语音设备和系统
US11064294B1 (en) 2020-01-10 2021-07-13 Synaptics Incorporated Multiple-source tracking and voice activity detections for planar microphone arrays
CN111175693A (zh) * 2020-01-19 2020-05-19 河北科技大学 一种波达方向估计方法及波达方向估计装置
CN112116920B (zh) * 2020-08-10 2022-08-05 北京大学 一种说话人数未知的多通道语音分离方法
CN113203987A (zh) * 2021-07-05 2021-08-03 成都启英泰伦科技有限公司 一种基于k均值聚类的多声源方位估计方法
US11823707B2 (en) 2022-01-10 2023-11-21 Synaptics Incorporated Sensitivity mode for an audio spotting system
DE202022105574U1 (de) 2022-10-01 2022-10-20 Veerendra Dakulagi Ein System zur Klassifizierung mehrerer Signale für die Schätzung der Ankunftsrichtung

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1565144A (zh) * 2001-08-08 2005-01-12 数字信号处理工厂有限公司 使用过采样滤波器组的定向音频信号处理
WO2005004532A1 (en) * 2003-06-30 2005-01-13 Harman Becker Automotive Systems Gmbh Handsfree system for use in a vehicle
CN101288334A (zh) * 2005-08-26 2008-10-15 思德普通信公司 用于使用衰减系数改进噪声识别的方法和设备

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2793629B1 (fr) * 1999-05-12 2001-08-03 Matra Nortel Communications Procede et dispositif d'annulation d'echo stereophonique a filtrage dans le domaine frequentiel
US8098844B2 (en) * 2002-02-05 2012-01-17 Mh Acoustics, Llc Dual-microphone spatial noise suppression
CA2399159A1 (en) * 2002-08-16 2004-02-16 Dspfactory Ltd. Convergence improvement for oversampled subband adaptive filters
JP4247037B2 (ja) * 2003-01-29 2009-04-02 株式会社東芝 音声信号処理方法と装置及びプログラム
EP1538867B1 (en) * 2003-06-30 2012-07-18 Nuance Communications, Inc. Handsfree system for use in a vehicle
US20080004729A1 (en) * 2006-06-30 2008-01-03 Nokia Corporation Direct encoding into a directional audio coding format
DE602006006664D1 (de) * 2006-07-10 2009-06-18 Harman Becker Automotive Sys Reduzierung von Hintergrundrauschen in Freisprechsystemen
RS49875B (sr) * 2006-10-04 2008-08-07 Micronasnit, Sistem i postupak za slobodnu govornu komunikaciju pomoću mikrofonskog niza
US8180062B2 (en) * 2007-05-30 2012-05-15 Nokia Corporation Spatial sound zooming
US20080298610A1 (en) * 2007-05-30 2008-12-04 Nokia Corporation Parameter Space Re-Panning for Spatial Audio
JP4455614B2 (ja) * 2007-06-13 2010-04-21 株式会社東芝 音響信号処理方法及び装置
ES2461601T3 (es) * 2007-10-09 2014-05-20 Koninklijke Philips N.V. Procedimiento y aparato para generar una señal de audio binaural
US8457328B2 (en) * 2008-04-22 2013-06-04 Nokia Corporation Method, apparatus and computer program product for utilizing spatial information for audio signal enhancement in a distributed network environment
US8724829B2 (en) * 2008-10-24 2014-05-13 Qualcomm Incorporated Systems, methods, apparatus, and computer-readable media for coherence detection
EP2380033B1 (en) * 2008-12-16 2017-05-17 Koninklijke Philips N.V. Estimating a sound source location using particle filtering
EP2375779A3 (en) * 2010-03-31 2012-01-18 Fraunhofer-Gesellschaft zur Förderung der Angewandten Forschung e.V. Apparatus and method for measuring a plurality of loudspeakers and microphone array
US9082391B2 (en) * 2010-04-12 2015-07-14 Telefonaktiebolaget L M Ericsson (Publ) Method and arrangement for noise cancellation in a speech encoder
US8639499B2 (en) * 2010-07-28 2014-01-28 Motorola Solutions, Inc. Formant aided noise cancellation using multiple microphones
RU2570359C2 (ru) * 2010-12-03 2015-12-10 Фраунхофер-Гезелльшафт Цур Фердерунг Дер Ангевандтен Форшунг Е.Ф. Прием звука посредством выделения геометрической информации из оценок направления его поступления
CN103688522B (zh) * 2011-05-18 2015-11-25 谷歌公司 时钟漂移补偿方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1565144A (zh) * 2001-08-08 2005-01-12 数字信号处理工厂有限公司 使用过采样滤波器组的定向音频信号处理
WO2005004532A1 (en) * 2003-06-30 2005-01-13 Harman Becker Automotive Systems Gmbh Handsfree system for use in a vehicle
CN101288334A (zh) * 2005-08-26 2008-10-15 思德普通信公司 用于使用衰减系数改进噪声识别的方法和设备

Also Published As

Publication number Publication date
RU2641319C2 (ru) 2018-01-17
EP2936830B8 (en) 2017-01-25
CN105165026A (zh) 2015-12-16
US10331396B2 (en) 2019-06-25
US20150286459A1 (en) 2015-10-08
EP2936830B1 (en) 2016-10-05
JP2016506664A (ja) 2016-03-03
ES2612528T3 (es) 2017-05-17
BR112015014380B1 (pt) 2022-10-11
RU2015129784A (ru) 2017-01-27
JP6196320B2 (ja) 2017-09-13
WO2014095250A1 (en) 2014-06-26
BR112015014380A2 (pt) 2020-01-28
EP2936830A1 (en) 2015-10-28
EP2747451A1 (en) 2014-06-25

Similar Documents

Publication Publication Date Title
CN105165026B (zh) 使用多个瞬时到达方向估计的知情空间滤波的滤波器及方法
TWI713844B (zh) 用於語音處理的方法及積體電路
Thiergart et al. An informed parametric spatial filter based on instantaneous direction-of-arrival estimates
Brandstein et al. Microphone arrays: signal processing techniques and applications
Habets et al. A two-stage beamforming approach for noise reduction and dereverberation
EP2868117B1 (en) Systems and methods for surround sound echo reduction
CA2819394C (en) Sound acquisition via the extraction of geometrical information from direction of arrival estimates
JP6703525B2 (ja) 音源を強調するための方法及び機器
US8958572B1 (en) Adaptive noise cancellation for multi-microphone systems
Thiergart et al. An informed LCMV filter based on multiple instantaneous direction-of-arrival estimates
US10250975B1 (en) Adaptive directional audio enhancement and selection
CN104902418A (zh) 用于估计目标和噪声谱方差的多传声器方法
Enzner Bayesian inference model for applications of time-varying acoustic system identification
CN111078185A (zh) 录制声音的方法及设备
JP2001309483A (ja) 収音方法および収音装置
Thiergart et al. An informed MMSE filter based on multiple instantaneous direction-of-arrival estimates
CN111354368B (zh) 补偿处理后的音频信号的方法
Thiergart et al. An acoustical zoom based on informed spatial filtering
As’ad et al. Beamforming designs robust to propagation model estimation errors for binaural hearing aids
Reindl et al. An acoustic front-end for interactive TV incorporating multichannel acoustic echo cancellation and blind signal extraction
Thiergart Flexible Multi-Microphone Acquisition and Processing of Spatial Sound Using Parametric Sound Field Representations
Herzog et al. Direction preserving wind noise reduction of b-format signals
Garre et al. An Acoustic Echo Cancellation System based on Adaptive Algorithm
Habets Towards multi-microphone speech dereverberation using spectral enhancement and statistical reverberation models
KR102573148B1 (ko) 사운드 교정을 위한 2-채널 공간 전달 함수의 지각적으로-투명한 추정

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant