基于参数化超增益波束形成器的音频信号提取方法及装置的制造方法
【专利摘要】本发明提供了一种基于参数化超增益波束形成器的音频信号提取方法及装置，其中方法包括：获取传感器阵列的结构信息，基于传感器阵列的结构信息，构造各向同性噪声的归一化自相关矩阵；在归一化自相关矩阵中引入参数，根据引入参数后的归一化自相关矩阵构建参数化超增益波束形成器；获取传感器阵列中各个传感器采集到的音频信号，该音频信号为时域信号；利用参数化超增益波束形成器对音频信号进行提取，得到来自期望方向的音频信号。通过本发明中的基于参数化超增益波束形成器的音频信号提取方法及装置，只需要在一个很小的范围内选取一个合适的参数，就能够在白噪声增益和指向性因子之间获得很好地平衡，实施简单，实用性强。
【专利说明】
基于参数化超増益波束形成器的音频信号提取方法及装置
技术领域
[0001] 本发明涉及声学信号处理领域，具体而言，涉及一种基于参数化超增益波束形成 器的音频信号提取方法及装置。
【背景技术】
[0002] 波束形成技术在信号处理领域有着广泛的应用，它通常是利用传感器阵列在空间 形成一个滤波器，该滤波器响应最大的方向对准期望方向，该滤波器对干扰的抑制取决于 干扰方向上滤波器响应的幅度。图1示出了相关技术中的传感器阵列的响应示意图，如图1 所示，该传感器阵列在0°方向的响应最大。
[0003] 波束形成大致可以分为两大类，一类是固定波束形成，另一类是自适应波束形成。 固定波束形成的滤波器系数不随应用环境的变化而变化，如延迟相加波束形成和超增益波 束形成，而自适应波束形成的滤波器系数随信号统计特性的变化而自适应的改变，如线性 约束最小方差波束形成。
[0004] 超增益波束形成器属于固定波束形成器，和其它的固定波束形成器相比，在同样 传感器数目的情况下，超增益波束形成器可以获得最大的指向性，因此能更为有效地抑制 来自主波束以外其它方向上的噪声。但在实际应用中，超增益波束形成面临最大的问题就 是白噪声放大问题，尤其是对于阵列间距很小的超增益波束形成，其低频处的白噪声增益 非常低，因而限制了其在实际系统中的应用。如何开发出稳健的超增益波束形成技术一直 是一个很具挑战性的问题。常见的解决办法是对噪声的归一化自相关矩阵进行对角加载， 通过在噪声的归一化自相关矩阵的对角线元素上加载一个小的常数(又称加载因子)来改 变该矩阵的条件数，从而提高白噪声增益。
[0005] 发明人在研究中发现，相关技术中基于对角加载的超增益波束形成虽然能够提高 白噪声增益，但是理论上加载因子的选择范围是从零到无穷大，在实际应用中，将很难去选 择一个合适的加载因子，导致基于对角加载的超增益波束形成技术实施困难，实用性不足。
【发明内容】
[0006] 有鉴于此，本发明的目的在于提供一种基于参数化超增益波束形成器的音频信号 提取方法及装置，只需要在一个很小的范围内选取一个合适的参数，就能够在白噪声增益 和指向性因子之间获得很好地平衡，实施简单，实用性强。
[0007] 第一方面，本发明实施例提供了一种基于参数化超增益波束形成器的音频信号提 取方法，所述方法包括:获取传感器阵列的结构信息，基于所述传感器阵列的结构信息，构 造各向同性噪声的归一化自相关矩阵;在所述归一化自相关矩阵中引入参数，根据引入参 数后的所述归一化自相关矩阵构建参数化超增益波束形成器;获取所述传感器阵列中各个 传感器采集到的音频信号，所述音频信号为时域信号；利用所述参数化超增益波束形成器 对所述音频信号进行提取，得到来自期望方向的音频信号。
[0008] 结合第一方面，本发明实施例提供了第一方面第一种可能的实施方式，其中，基于 所述传感器阵列的结构信息，构造各向同性噪声的归一化自相关矩阵，包括:根据所述传感器 阵列中传感器的数量Μ和相邻传感器之间的间距δ构造 Μ XM的各向同性噪声的归一化自相关矩 阵Γ ( ω )，矩阵Γ ( ω )的第（i j)个元素表示为
其中，ω为角频率，τ〇为相邻传感器之间的最大声传播时延，τ〇 = δ/c，c为所述音频信号在空 气中的传播速度。
[0009] 结合第一方面，本发明实施例提供了第一方面第二种可能的实施方式，其中，在所 述归一化自相关矩阵中引入参数，根据引入参数后的所述归一化自相关矩阵构建参数化超 增益波束形成器，包括:在所述归一化自相关矩阵中引入参数，根据引入参数后的所述归一 化自相关矩阵定义参数化信噪比增益;在期望方向信号不失真约束条件下最大化所述参数 化信噪比增益，得到所述参数化超增益波束形成器。
[0010] 结合第一方面第二种可能的实施方式，本发明实施例提供了第一方面第三种可能 的实施方式，其中，在所述归一化自相关矩阵中引入参数，根据引入参数后的所述归一化自 相关矩阵定义参数化信噪比增益，包括:对所述归一化自相关矩阵Γ ( ω )进行特征分解，得 到Γ ( ω )=υ( ω ) Λ ( ω )υΤ( ω )，其中，ω为角频率，υ( ω )为正交矩阵，Λ ( ω )为对角矩阵， UT( ω )为矩阵U( ω )的转置;基于所述特征分解的结果，在所述归一化自相关矩阵Γ ( ω )中 引入参数Ρ，定义Ι/p阶归一化自相关矩阵为
3所述 Ι/p阶归一化自相关矩阵，P e [ 1，3];根据所述Ι/p阶归一化自相关矩阵定义所述参数化信 噪比增益为
其中，GP[h(co)]为所述参数化信噪比增益，h (ω )为长度为Μ的线性滤波器，Μ为所述传感器阵列中传感器的数量，hH( ω )为h( ω )的共辄 转置，d( ω )为所述传感器阵列的导向矢量。
[0011] 结合第一方面第二种可能的实施方式，本发明实施例提供了第一方面第四种可能 的实施方式，其中，在期望方向信号不失真约束条件下最大化所述参数化信噪比增益，得到 所述参数化超增益波束形成器，包括:求解优化问I
subject to hH (ω )d( ω ) = 1的解，得到所述参数化超增益波束形成器
；其中， hP,P( ω )为所述参数化超增益波束形成器，ω为角频率，h( ω )表示长度为Μ的线性滤波器，Μ 为所述传感器阵列中传感器的数量，hH( ω )为h( ω )的共辄转置，阶归一化自相 关矩阵，d( ω )为所述传感器阵列的导向矢量，dH( ω )为d( ω )的共辄转置
归一化自相关矩阵，pe [1，3]。
[0012] 结合第一方面，本发明实施例提供了第一方面第五种可能的实施方式，其中，在利 用所述参数化超增益波束形成器对所述音频信号进行提取之前，所述方法还包括:应用短 时傅里叶变换将所述音频信号由时域信号转换为频域信号;利用所述参数化超增益波束形 成器对所述音频信号进行提取，包括:利用所述参数化超增益波束形成器对所述频域信号 进行提取。
[0013] 第二方面，本发明实施例提供了一种基于参数化超增益波束形成器的音频信号提 取装置，矩阵构造模块，用于获取传感器阵列的结构信息，基于所述传感器阵列的结构信 息，构造各向同性噪声的归一化自相关矩阵;波束形成器构建模块，用于在所述归一化自相 关矩阵中引入参数，根据引入参数后的所述归一化自相关矩阵构建参数化超增益波束形成 器;信号获取模块，用于获取所述传感器阵列中各个传感器采集到的音频信号，所述音频信 号为时域信号;信号提取模块，用于利用所述参数化超增益波束形成器对所述音频信号进 行提取，得到来自期望方向的音频信号。
[0014] 结合第二方面，本发明实施例提供了第二方面第一种可能的实施方式，其中，所述 波束形成器构建模块包括:信噪比增益定义单元，用于在所述归一化自相关矩阵中引入参 数，根据引入参数后的所述归一化自相关矩阵定义参数化信噪比增益;波束形成器构建单 元，用于在期望方向信号不失真约束条件下最大化所述参数化信噪比增益，得到所述参数 化超增益波束形成器。
[0015] 结合第二方面第一种可能的实施方式，本发明实施例提供了第二方面第二种可能 的实施方式，其中，所述信噪比增益定义单元包括:特征分解子单元，用于对所述归一化自 相关矩阵r ( ω )进行特征分解，得到Γ ( ω ) =υ( ω ) Λ ( ω )υτ( ω )，其中，ω为角频率，u (ω )为正交矩阵，Λ ( ω )为对角矩阵，UT( ω )为矩阵U( ω )的转置;矩阵定义子单元，用于基 于所述特征分解的结果，在所述归一化自相关矩阵Ρ(ω)中引入参数p，定义Ι/p阶归一化 自相关矩阵为
为所述Ι/p阶归一化自相关矩阵，P e
[1，3];信噪比增益定义子单元，用于根据所述l/ρ阶归一化自相关矩阵定义所述参数化信 噪比增益为
其中，GP[h( ω )]为所述参数化信噪比增益，h( ω ) 为长度为Μ的线性滤波器，Μ为所述传感器阵列中传感器的数量，hH( ω )为h( ω )的共辄转 置，d( ω )为所述传感器阵列的导向矢量。
[0016] 结合第二方面第一种可能的实施方式，本发明实施例提供了第二方面第三种可能的 实施方式，其中，所述波束形成器构建单元用于:求解优化问题min h" 〇>ΓΛsub ject h(./M to hH( ω )d( ω ) = 1的解，得到所述参数化超增益波束形成署
其中，hP,P( ω )为所述参数化超增益波束形成器，ω为角频率，h( ω )表示长度为Μ的线性滤 波器，Μ为所述传感器阵列中传感器的数量，hH( ω )为h( ω )的共辄转置，为l/ρ阶归一 化自相关矩阵，d( ω )为所述传感器阵列的导向矢量，dH( ω )为d( ω )的共辄转置， 阶归一化自相关矩阵，pe [1，3]。 P
[0017] 本发明实施例中，首先构造各向同性噪声的归一化自相关矩阵，然后在归一化自 相关矩阵中引入参数构建参数化超增益波束形成器，最后利用参数化超增益波束形成器对 传感器阵列中各个传感器采集到的音频信号进行提取，得到来自期望方向的音频信号。通 过本发明实施例中的基于参数化超增益波束形成器的音频信号提取方法及装置，只需要在 一个很小的范围内，如1至3的范围内，选取一个合适的参数构建参数化超增益波束形成器， 就能够在白噪声增益和指向性因子之间获得很好地平衡，实现复杂环境下对期望方向信号 的获取，实施简单，实用性强，从而有效缓解相关技术中加载因子选择困难，基于对角加载 的超增益波束形成技术实施困难，实用性不足的问题。
[0018] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合 所附附图，作详细说明如下。
【附图说明】
[0019] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附 图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对 范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这 些附图获得其他相关的附图。
[0020] 图1示出了相关技术中的传感器阵列的响应示意图；
[0021] 图2示出了本发明实施例1提供的一种基于参数化超增益波束形成器的音频信号 提取方法的流程示意图；
[0022] 图3示出了本发明实施例1提供的另一种基于参数化超增益波束形成器的音频信 号提取方法的流程示意图；
[0023] 图4示出了本发明实施例2提供的一种基于参数化超增益波束形成器的音频信号 提取装置的模块组成示意图；
[0024] 图5示出了本发明实施例2提供的另一种基于参数化超增益波束形成器的音频信 号提取装置的模块组成示意图。
【具体实施方式】
[0025]为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅 是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实 施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的 实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实 施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所 有其他实施例，都属于本发明保护的范围。
[0026]考虑到相关技术中基于对角加载的超增益波束形成技术实施困难，实用性不足的 问题，本发明提供了一种基于参数化超增益波束形成器的音频信号提取方法及装置，下面 结合实施例进行具体描述。
[0027] 实施例1
[0028] 本实施例提供了一种基于参数化超增益波束形成器的音频信号提取方法，本发明 实施例的执行主体为音频信号提取服务器，音频信号提取服务器首先构造各向同性噪声的 归一化自相关矩阵，其次在归一化自相关矩阵中引入参数构建参数化超增益波束形成器， 然后获取传感器阵列中各个传感器采集到的音频信号，最后利用参数化超增益波束形成器 对音频信号进行提取，得到来自期望方向的音频信号。其中，来自期望方向的音频信号简称 为期望方向信号，音频信号提取服务器可以使用现有的任何型号的计算终端，这里不再一 一赘述。
[0029] 图2示出了本发明实施例1提供的一种基于参数化超增益波束形成器的音频信号 提取方法的流程示意图，如图2所示，该方法包括以下步骤：
[0030] 步骤S202，获取传感器阵列的结构信息，基于该传感器阵列的结构信息，构造各向 同性噪声的归一化自相关矩阵。
[0031] 本实施例中，传感器阵列可以为任意形状，包括但不限于方形阵列、圆形阵列等。 本实施例中的方法既可以用来处理窄带信号亦可以处理宽带信号如语音信号。
[0032] 本步骤中，根据传感器阵列中传感器的数量Μ和相邻传感器之间的间距δ构造 MXM 的各向同性噪声的归一化自相关矩阵Γ ( ω )，矩阵Γ ( ω )的第(ij)个元素表示为：
[0033]
⑴
[0034]公式（1)，ω为角频率，ω =2Jif，f>〇为时间频率为相邻传感器之间的最大声传 播时延，τ〇 = δ/c，c为音频信号在空气中的传播速度。
[0035]具体地，步骤S202中，假设一均匀线性传感器阵列含有Μ个全指向性传感器，声源 以平面波的形式传播，其入射方位和传感器阵列呈角度Θ。理想情况下各个传感器上接收到 的信号刚好相差一个时延，此时阵列的导向矢量可以表不为：
[0036]
(2)
[0037] 上式(2)中上标Τ表示对向量的转置操作，/ = #是虚数单位。
[0038] 对于超增益波束形成，假设相邻传感器的间距δ很小，且声源信号从端射方向入 射，g卩θ = 〇,因此，麦克风阵列接收信号可以表示为：
[0039] γ(ω) = [Υ!(ω) Υ2(ω) ... ΥΜ(ω)]τ
[0040] =χ( ω )+ν( ω )
[0041] =d( ω )Χ( ω )+ν( ω )， (3)
[0042] 这里：α?) = 是第m个传感器接收到的信号，Χ( ω )是期望方 向信号，Vm( ω )是第m个传感器上的加性噪声，χ( ω ) = d( ω )Χ( ω )，这里d( ω ) = d( ω，〇)，ν (ω )类似y( ω )定义。
[0043] 通常，波束形成是对各传感器接收到的信号乘以一个复数权值再相加得到输出， 即
[0044]
(4)
[0045] 这里Ζ( ω )是对期望方向信号Χ( ω )的估计，上标$表示取共辄，h( ω )是长度为Μ的 线性滤波器，上标η是共辄转置操作。同时需要期望方向信号不失真，该不失真约束定义为hH (ω)(1(ω) = 1〇
[0046] 为了推导出最优波束形成器，本实施例定义几个基本的性能评价参数，它们分别 是波束图、白噪声增益、指向性因子。
[0047] (1)波束图反应了阵列对平面波从Θ方向入射时的响应，它定义为 [00481
(5)
[0049] (2)白噪声增益反应了阵列对阵列和传感器对各种不匹配的鲁棒性，它定义为如 下形式
[0050]
(6)
[0051] (3)指向性因子反应了波束形成器的方向性增益，它被定义为
[0052]
. (7)
[0053] 这里上式中ΜΧΜ的各向同性噪声的归一化自相关矩阵Γ(ω)的第（ij)个元素如 公式(1)所示，表示为
[0054]步骤S204,在归一化自相关矩阵中引入参数，根据引入参数后的归一化自相关矩 阵构建参数化超增益波束形成器。
[0055] 本步骤能够通过以下两个子步骤(1)和(2)实现：
[0056] (1)在归一化自相关矩阵中引入参数，根据引入参数后的归一化自相关矩阵定义 参数化信噪比增益。
[0057]本子步骤(1)具体包括：
[0058] (11)对归一化自相关矩阵Γ ( ω )进行特征分解，得到
[0059] Γ ( ω )=U( ω ) Λ ( ω )UT( ω )， (8)
[0060] 其中，ω为角频率，U( ω )为正交矩阵，UT( ω )为矩阵U( ω )的转置，UT( ω )U( ω )=U (ω )UT( ω ) = !"，Λ ( ω )为对角矩阵，Λ ( ω )的所有对角元素都是正值。
[0061] (12)基于特征分解的结果，在归一化自相关矩阵Γ(ω)中引入参数p，定义l/ρ阶 归一化自相关矩阵为
[0062]
(9)
[0063] 其中，为l/ρ阶归一化自相关矩阵，Pe [ 1，00]，优选Pe [ 1，3]。
[0064] (13)根据l/ρ阶归一化自相关矩阵定义参数化信噪比增益为
[0065]
(.10.)
[0066] 其中，GP[h( ω )]为参数化信噪比增益，h( ω )为长度为Μ的线性滤波器，Μ为传感器 阵列中传感器的数量，hH( ω )为h( ω )的共辄转置，d( ω )为传感器阵列的导向矢量。这里很 显然有G~[h( ω ) ] =W[h( ω ) ]andGi[h( ω ) ] =D[h( ω ) ] 0
[0067] (2)在期望方向信号不失真约束条件下最大化参数化信噪比增益，得到参数化超 增益波束形成器。
[0068] 本子步骤具体执行为：求解优化问题
. h(.iy) (ω ) = 1的解，得到参数化超增益波束形成器
[0069]
(11)
[0070] 其中，hP,P( ω )为参数化超增益波束形成器，ω为角频率，h( ω )表示长度为Μ的线 性滤波器，Μ为传感器阵列中传感器的数量，hH( ω )为h( ω )的共辄转置，为l/ρ阶归一 化自相关矩阵，d( ω )为传感器阵列的导向矢量，dH( ω )为d( ω )的共辄转置
阶归一化自相关矩阵，p e [ 1，00 ]，优选p e [ 1，3]。
[0071] 这里l/ρ是参数超增益波束形成的参数阶数，参数矩阵!··^#可如下计算：
[0072]
(12)
[0073] 至此，能够得到参数化超增益波束形成器的白噪声增益为
[0074]
(13)
[0075] 参数化超增益波束形成器的指向性系数为
[0076]
(M)
[0077] 对于任意的Pl$p2,总有
[0078] (615)
[0079] (16)
[0080] 很显然，通过调节参数P，可以得到不同形式的波束形成器：当P=1时，能够得到传 统的超增益波束形成器，即h P>1( ω )=hs( ω )，它在给定传感器数目的前提下能得到最大的 指向性因子;当P = m时，能够得到延时相加波束形成器，即，它在给定传感器数目的前提下 能得到最大的白噪声增益；当l〈P〈m时，能够得到折衷的波束形成器，它的白噪声增益随着 参数P的增加而增加，指向性因子随着参数P增加而减小，如果选择合适的参数P，参数化超 增益波束形成器h P,P( ω )能够在很好地改善白噪声放大问题的前提下同时保持合理的指向 性。
[0081] 步骤S206,获取上述传感器阵列中各个传感器采集到的音频信号，该音频信号为 时域信号。
[0082] 步骤S208,利用参数化超增益波束形成器对上述音频信号进行提取，得到来自期 望方向的音频信号。
[0083] 本步骤具体可以为，利用公对上述音频信号进行提 取，得到来自期望方向的音频信号。
[0084] 本步骤中，能够利用参数化超增益波束形成器对时域信号进行提取，得到来自期 望方向的音频信号，还能够利用参数化超增益波束形成器对频域信号进行提取，得到来自 期望方向的音频信号。当对频域信号提取时，在步骤S208之前，本实施例中的方法还包括： 应用短时傅里叶变换将上述音频信号由时域信号转换为频域信号;此时，利用参数化超增 益波束形成器对音频信号进行提取，具体执行为:利用参数化超增益波束形成器对频域信 号进行提取。
[0085]需要说明的是，本实施例中，对传感器接收到的信号进行傅里叶变换(对于非平稳 的信号如语音信号，采取分帧处理进行短时傅里叶变换），对变换后的频域信号在各个频段 分开处理，如果是窄带信号就只需要在有信号能量的频段处理，如果是宽带信号则在全频 段或者信号能量集中的频段进行处理。
[0086]本发明实施例中，首先构造各向同性噪声的归一化自相关矩阵，然后在归一化自 相关矩阵中引入参数构建参数化超增益波束形成器，最后利用参数化超增益波束形成器对 传感器阵列中各个传感器采集到的音频信号进行提取，得到来自期望方向的音频信号。通 过本发明实施例中的基于参数化超增益波束形成器的音频信号提取方法，只需要在一个很 小的范围内，如1至3的范围内，选取一个合适的参数构建参数化超增益波束形成器，就能够 在白噪声增益和指向性因子之间获得很好地平衡，实现复杂环境下对期望方向信号的获 取，实施简单，实用性强，从而有效缓解相关技术中加载因子选择困难，基于对角加载的超 增益波束形成技术实施困难，实用性不足的问题。
[0087] 需要说明的是，本实施例中步骤S202至步骤S204所示的构建参数化超增益波束形 成器的过程与步骤S206所示的获取音频信号的过程可以同时进行，也可以先获取音频信 号，后构建参数化超增益波束形成器，也可以如图2所示先构建参数化超增益波束形成器， 后获取音频信号。因此，本领域技术人员应当理解本实施例中的方案中，步骤S202至步骤 S204所示的构建参数化超增益波束形成器的过程与步骤S206所示的获取音频信号的过程 的先后顺序可以不进行具体限定。
[0088] 为进一步说明构建参数化超增益波束形成器的过程与获取音频信号的过程之间 的逻辑关系，本实施例还提供了如图3所示的方法流程，该流程中，构建参数化超增益波束 形成器与获取音频信号同时进行，图3所示的方法流程，具备图2所示的方法流程相同的技 术效果。
[0089] 为了更好的展示本发明实施例的效果，这里利用一个含有8个全指向性麦克风的 均匀线性阵列，麦克风的间距为1.5cm，依照如图2所示的方法设计参数化超增益波束形成 器，并用波束图、白噪声增益和指向性因子对其性能进行评价。
[0090] 通过实验能够发现，本实施例中的参数化超增益波束形成器当参数p = 1时，是传 统的超增益波束形成器，波束图在θ = 〇°度上的响应为1，在0°到180°的范围上有三个零点， 此时的波束图对应为三阶的超心形，这和理论上Μ个麦克风设计的超增益波束形成器对应 Μ-1超心形相符合。随着参数ρ的增加，波束图由超心形向延时相加波束图过度，当参数ρ的 取值较大(如P = l〇)时，此时的波束图已经非常接近延时相加波束图。
[0091] 通过实验还能够发现，与传统的超增益波束形成器相比，本实施例中的参数化超 增益波束形成器能够在高的指向性因子和低的白噪声增益之间获得很好的平衡。随着参数 Ρ的值从1.1增加到10,指向性因子的值相应减小而白噪声增益的值显著增加。
[0092] 通过实验还能够发现，指向性因子随着参数ρ的增加而减小，白噪声增益随着参数 Ρ的增加而增加，尤其值得注意的是，其随着参数Ρ的显著变化集中在1到3这个很小的范围 内，当Ρ的值大于3时，指向性因子和白噪声增益都不再显著变化，这说明本实施例可以通过 在一个很小的范围内选取合适的参数Ρ去设计最优的参数化超增益波束形成器。
[0093] 实施例2
[0094] 对应实施例1中的基于参数化超增益波束形成器的音频信号提取方法，本实施例 提供了 一种基于参数化超增益波束形成器的音频信号提取装置，如图4所示，该装置包括：
[0095] 矩阵构造模块31，用于获取传感器阵列的结构信息，基于传感器阵列的结构信息， 构造各向同性噪声的归一化自相关矩阵；
[0096] 波束形成器构建模块32,用于在归一化自相关矩阵中引入参数，根据引入参数后 的归一化自相关矩阵构建参数化超增益波束形成器；
[0097]信号获取模块33,用于获取上述传感器阵列中各个传感器采集到的音频信号，该 音频信号为时域信号；
[0098] 信号提取模块34,用于利用参数化超增益波束形成器对上述音频信号进行提取， 得到来自期望方向的音频信号。
[0099] 上述矩阵构造模块31具体用于:根据传感器阵列中传感器的数量Μ和相邻传感器 之间的间距S构造 MXM的各向同性噪声的归一化自相关矩阵Γ ( ω )，矩阵Γ ( ω )的第（ij) 个元素表示为
1，其中，ω为角频率，为相邻传感 器之间的最大声传播时延，t〇 = S/c，c为音频信号在空气中的传播速度。
[0100]上述波束形成器构建模块32包括:信噪比增益定义单元，用于在归一化自相关矩 阵中引入参数，根据引入参数后的归一化自相关矩阵定义参数化信噪比增益;波束形成器 构建单元，用于在期望方向信号不失真约束条件下最大化参数化信噪比增益，得到参数化 超增益波束形成器。
[0101 ]其中，信噪比增益定义单元包括:特征分解子单元，用于对归一化自相关矩阵Γ ( ω ) 进行特征分解，得到Γ ( ω ) =U( ω ) Λ ( ω )UT( ω )，其中，ω为角频率，U( ω )为正交矩阵，Λ ( ω ) 为对角矩阵，UT( ω )为矩阵U( ω )的转置;矩阵定义子单元，用于基于特征分解的结果，在归一化 自相关矩阵F ( ω )中引入参数p，定义Ι/p阶归一化自相关矩阵为
阶归一化自相关矩阵定义参数化信噪比增益为
其中，GP[h (ω )]为参数化信噪比增益，h( ω )为长度为Μ的线性滤波器，Μ为传感器阵列中传感器的数 量，hH( ω )为h( ω )的共辄转置，d( ω )为传感器阵列的导向矢量。
[0102] 其中，波束形成器构建单元用于:求解优化问是
ubject to hH( ω )d( ω ) = 1的解，得到参数化超增益波束形成署
其中， hP,P( ω )为参数化超增益波束形成器，ω为角频率，h( ω )表示长度为Μ的线性滤波器，Μ为传 感器阵列中传感器的数量，hH( ω )为h( ω )的共辄转置，j4(iy)Sl/p阶归一化自相关矩阵，d (ω )为传感器阵列的导向矢量，dH( ω )为d( ω )的共辄转置，丨^^为-^阶归一化自相关 矩阵，Pe [1，3]。
[0103] 当对频域信号进行处理时，本实施例中的装置还包括:傅里叶变换模块，用于应用 短时傅里叶变换将音频信号由时域信号转换为频域信号。该种情况下，上述信号提取模块 34具体用于:利用参数化超增益波束形成器对频域信号进行提取。
[0104] 本发明实施例中，首先构造各向同性噪声的归一化自相关矩阵，然后在归一化自 相关矩阵中引入参数构建参数化超增益波束形成器，最后利用参数化超增益波束形成器对 传感器阵列中各个传感器采集到的音频信号进行提取，得到来自期望方向的音频信号。通 过本发明实施例中的基于参数化超增益波束形成器的音频信号提取装置，只需要在一个很 小的范围内，如1至3的范围内，选取一个合适的参数构建参数化超增益波束形成器，就能够 在白噪声增益和指向性因子之间获得很好地平衡，实现复杂环境下对期望方向信号的获 取，实施简单，实用性强，从而有效缓解相关技术中加载因子选择困难，基于对角加载的超 增益波束形成技术实施困难，实用性不足的问题。
[0105] 需要说明的是，本实施例中矩阵构造模块31和波束形成器构建模块32所示的构建 参数化超增益波束形成器的过程与信号获取模块33所示的获取音频信号的过程可以同时 进行，也可以先获取音频信号，后构建参数化超增益波束形成器，也可以如图4所示先构建 参数化超增益波束形成器，后获取音频信号。因此，本领域技术人员应当理解本实施例中的 方案中，波束形成器构建模块32可以通过信号获取模块33与信号提取模块34相连，也可以 直接与信号提取模块34相连。
[0106] 为进一步说明构建参数化超增益波束形成器的过程与获取音频信号的过程之间 的逻辑关系，本实施例还提供了如图5所示的装置模块组成，该装置中，构建参数化超增益 波束形成器与获取音频信号同时进行，波束形成器构建模块32直接与信号提取模块34相 连，图5所示的装置组成，具备图4所示的装置组成相同的技术效果。
[0107] 本发明实施例所提供的基于参数化超增益波束形成器的音频信号提取装置可以 为设备上的特定硬件或者安装于设备上的软件或固件等。本发明实施例所提供的装置，其 实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及 之处，可参考前述方法实施例中相应内容。所属领域的技术人员可以清楚地了解到，为描述 的方便和简洁，前述描述的系统、装置和单元的具体工作过程，均可以参考上述方法实施例 中的对应过程，在此不再赘述。
[0108] 在本发明所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方 式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻 辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可 以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间 的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连 接，可以是电性，机械或其它的形式。
[0109] 所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显 示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个 网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目 的。
[0110] 另外，在本发明提供的实施例中的各功能单元可以集成在一个处理单元中，也可 以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
[0111] 所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以 存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说 对相关技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计 算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个 人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。 而前述的存储介质包括:U盘、移动硬盘、只读存储器(R0M，Read-0nly Memory)、随机存取存 储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0112] 应注意到:相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一 个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释，此外，术语"第 一"、"第二"、"第三"等仅用于区分描述，而不能理解为指示或暗示相对重要性。
[0113]最后应说明的是：以上所述实施例，仅为本发明的【具体实施方式】，用以说明本发明 的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发 明进行了详细的说明，本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员 在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻 易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使 相应技术方案的本质脱离本发明实施例技术方案的精神和范围。都应涵盖在本发明的保护 范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
【主权项】
1. 一种基于参数化超增益波束形成器的音频信号提取方法，其特征在于，所述方法包 括： 获取传感器阵列的结构信息，基于所述传感器阵列的结构信息，构造各向同性噪声的 归一化自相关矩阵； 在所述归一化自相关矩阵中引入参数，根据引入参数后的所述归一化自相关矩阵构建 参数化超增益波束形成器； 获取所述传感器阵列中各个传感器采集到的音频信号，所述音频信号为时域信号； 利用所述参数化超增益波束形成器对所述音频信号进行提取，得到来自期望方向的音 频信号。2. 根据权利要求1所述的方法，其特征在于，基于所述传感器阵列的结构信息，构造各 向同性噪声的归一化自相关矩阵，包括： 根据所述传感器阵列中传感器的数量Μ和相邻传感器之间的间距δ构造 MXM的各向同 性噪声的归一化自相关矩阵Γ ( ω )，矩阵Γ ( ω )的第α j)个元素表示为：其中，ω为角频率，το为相邻传感器之间的最大声传播时延，T〇 = S/c，c为所述音频信号 在空气中的传播速度。3. 根据权利要求1所述的方法，其特征在于，在所述归一化自相关矩阵中引入参数，根 据引入参数后的所述归一化自相关矩阵构建参数化超增益波束形成器，包括： 在所述归一化自相关矩阵中引入参数，根据引入参数后的所述归一化自相关矩阵定义 参数化信噪比增益； 在期望方向信号不失真约束条件下最大化所述参数化信噪比增益，得到所述参数化超 增益波束形成器。4. 根据权利要求3所述的方法，其特征在于，在所述归一化自相关矩阵中引入参数，根 据引入参数后的所述归一化自相关矩阵定义参数化信噪比增益，包括： 对所述归一化自相关矩阵Γ ( ω )进行特征分解，得到Γ ( ω ) =U( ω ) Λ ( ω )IJT( ω )，其 中，ω为角频率，U( ω )为正交矩阵，Λ ( ω )为对角矩阵，ljT( ω )为矩阵U( ω )的转置； 基于所述特征分解的结果，在所述归一化自相关矩阵Γ ( ω )中引入参数P，定义1/p阶 归一化自相关矩阵夫^所述1/P阶归一化自相关矩 阵，pe [1,3]; 根据所述1/p阶归一化自相关矩阵定义所述参数化信噪比增益为其中，Gp比（ω )]为所述参数化信噪比增益，K ω )为长度为Μ的线性滤波器，Μ为所述传感器 阵列中传感器的数量，hH( ω )为Κ ω )的共辆转置，d( ω )为所述传感器阵列的导向矢量。5. 根据权利要求3所述的方法，其特征在于，在期望方向信号不失真约束条件下最大化 所述参数化信噪比增益，得到所述参数化超增益波束形成器，包括： 求解优化问题的解，得到所述参数 化超增益波束形成署其中，hp,p( ω )为所述参数化超增益波束形成器，ω为角频率，Μ ω )表示长度为Μ的线 性滤波器，Μ为所述传感器阵列中传感器的数量，hH( ω )为Μ ω )的共辆转置，阶归一化自相关矩阵，d( ω )为所述传感器阵列的导向矢量，dH( ω )为d( ω )的共辆转置，济归一化自相关矩阵，Ρ e [ 1,3 ]。6. 根据权利要求1所述的方法，其特征在于，在利用所述参数化超增益波束形成器对所 述音频信号进行提取之前，所述方法还包括： 应用短时傅里叶变换将所述音频信号由时域信号转换为频域信号； 利用所述参数化超增益波束形成器对所述音频信号进行提取，包括： 利用所述参数化超增益波束形成器对所述频域信号进行提取。7. -种基于参数化超增益波束形成器的音频信号提取装置，其特征在于，所述装置包 括： 矩阵构造模块，用于获取传感器阵列的结构信息，基于所述传感器阵列的结构信息，构 造各向同性噪声的归一化自相关矩阵； 波束形成器构建模块，用于在所述归一化自相关矩阵中引入参数，根据引入参数后的 所述归一化自相关矩阵构建参数化超增益波束形成器； 信号获取模块，用于获取所述传感器阵列中各个传感器采集到的音频信号，所述音频 信号为时域信号； 信号提取模块，用于利用所述参数化超增益波束形成器对所述音频信号进行提取，得 到来自期望方向的音频信号。8. 根据权利要求7所述的装置，其特征在于，所述波束形成器构建模块包括： 信噪比增益定义单元，用于在所述归一化自相关矩阵中引入参数，根据引入参数后的 所述归一化自相关矩阵定义参数化信噪比增益； 波束形成器构建单元，用于在期望方向信号不失真约束条件下最大化所述参数化信噪 比增益，得到所述参数化超增益波束形成器。9. 根据权利要求8所述的装置，其特征在于，所述信噪比增益定义单元包括： 特征分解子单元，用于对所述归一化自相关矩阵Γ ( ω )进行特征分解，得到Γ ( ω )=U (ω ) Λ ( ω )IJT( ω )，其中，ω为角频率，U( ω )为正交矩阵，Λ ( ω )为对角矩阵，IJT( ω )为矩阵 υ(ω)的转置； 矩阵定义子单元，用于基于所述特征分解的结果，在所述归一化自相关矩阵Γ ( ω )中 引入参数Ρ，定义1/Ρ阶归一化自相关矩阵天弓所述 1/ρ阶归一化自相关矩阵，pe [1,3]; 信噪比增益定义子单元，用于根据所述1/p阶归一化自相关矩阵定义所述参数化信噪 比增益为庚中，Gp[h( ω )]为所述参数化信噪比增益，Μ ω ) 为长度为Μ的线性滤波器，Μ为所述传感器阵列中传感器的数量，hH( ω )为h( ω )的共辆转 置，d( ω )为所述传感器阵列的导向矢量。10.根据权利要求8所述的装置，其特征在于，所述波束形成器构建单元用于： 求解优化问题的解，得到所述参数化 超增益波束形成I其中，hp,p( ω )为所述参数化超增益波束形成器，ω为角频率，Μ ω )表示长度为Μ的线 性滤波器，Μ为所述传感器阵列中传感器的数量，hH( ω )为Κ ω )的共辆转置，为1/Ρ阶 归一化自相关矩阵，d( ω )为所述传感器阵列的导向矢量，dH( ω )为d( ω )的共辆转置，新归一化自相关矩阵，pe[l，3]。
【文档编号】G10K11/34GK105976822SQ201610545565
【公开日】2016年9月28日
【申请日】2016年7月12日
【发明人】陈景东, 黄公平, 雅各布·贝内斯蒂
【申请人】西北工业大学