[论文解读] Acoustic intensity, energy-density and diffuseness estimation in a directionally-constrained region

Archontis Politis, Ville Pulkki|arXiv (Cornell University)|Sep 12, 2016

Speech and Audio Processing参考文献 11被引用 10

一句话总结

该论文提出了一种球谐域方法，通过空间波束成形，在方向受限区域内估计声强、能量密度和扩散度。通过使用波束成形加权声场，该方法实现了局部化能量分析，使波束成形能够聚焦于特定的方位区域，实现波达方向（DOA）估计和声场表征，且对平面波、扩散场及混合场情况均提供解析解。

ABSTRACT

This work presents a method for estimation of the acoustic intensity, the energy density and the associated sound field diffuseness around the origin, when the sound field is weighted with a spatial filter. The method permits energetic DOA estimation and sound field characterization focused in a specific angular region determined by the beam pattern of the spatial filter. The formulation of the estimators is presented and their behavior is analyzed for the fundamental cases useful in parametric sound field models of a single plane wave, a uniform diffuse field and a mixture of the two.

研究动机与目标

• 通过空间滤波将估计限制在特定方位区域，实现声场的局部化能量分析。
• 建立基于球谐域表示的数学框架，用于在方向受限区域内估计声强、能量密度与扩散度。
• 分析在基本声场模型下的估计器行为：单一直角平面波、均匀扩散场及其混合场。
• 量化波束成形设计对混合声场中DOA估计偏差与扩散度估计的影响。
• 为波束成形系数提供闭式解，从而实现从空间滤波信号中准确估计质点速度与强度分量。

提出的方法

• 使用单位球面上的波束成形加权积分，在球谐（SH）域中表述声强、能量密度与扩散度。
• 将压力与速度分量表示为SH展开，波束成形表示为函数 $ w(\theta, \rho) $，从而实现对声场的方向加权。
• 推导出加权强度矢量 $ \mathbf{i}_{\text{pd,w}} $、能量密度 $ E_{\text{pd,w}} $ 与扩散度 $ \psi_{\text{pd,w}} $ 的闭式表达式，其形式基于波束成形与场分量的SH系数。
• 利用球谐变换（SHT）表示振幅函数 $ a(f, \Omega) $，从而高效计算空间加权的能量量。
• 提出一种设计框架，其中波束成形系数用于计算SH系数的线性组合，以估计质点速度分量。
• 利用余弦定理与矢量几何，推导出轴对称波束成形与混合场情况下扩散度与DOA偏差的解析表达式。

实验结果

研究问题

• RQ1 如何通过空间波束成形在方向受限区域内估计声强、能量密度与扩散度？
• RQ2 在存在扩散噪声的情况下，空间波束成形对平面波波达方向（DOA）估计有何影响？
• RQ3 波束成形如何影响由平面波与扩散成分组成的混合声场中扩散度的估计？
• RQ4 所提方法能否在标准场模型下，于球谐域中为声强与能量密度提供解析解？
• RQ5 波束成形直接性、扩散场功率与DOA估计偏差之间的关系为何？

主要发现

• 当波束成形主瓣与平面波传播方向对齐时，空间滤波下单一直角平面波的DOA保持不变；但当存在扩散场能量时，DOA估计会向波束成形方向产生偏差。
• 在均匀扩散场存在时，估计的强度矢量不会为零，而是指向波束成形的空间平均方向；若波束成形为轴对称，则方向直接指向主瓣。
• 在包含平面波与扩散成分的混合场中，估计的强度矢量会向波束成形方向产生偏差，且偏差随扩散场能量增加而增大，并依赖于平面波DOA与波束轴之间的夹角 $ \alpha $。
• 当平面波不存在或波束成形在平面波传播方向形成零点时，扩散度估计 $ \psi_{\text{pd,w}} $ 达到最大值（等于扩散场值）。
• 当平面波与波束主瓣对齐（即 $ \alpha = 0 $）时，扩散度达到最小值，表达式为 $ \psi_{\text{pd,w}}(\Gamma, 0) = \frac{\psi_{\text{df,w}}}{Q\Gamma + 1} $，其中 $ \Gamma $ 为能量比，$ Q $ 为SH系数数量。
• DOA偏差 $ \beta $ 的解析表达式为 $ \beta(\Gamma, \alpha) = \arcsin\left( \frac{K \sin\alpha}{4\pi \Delta(\Gamma, \alpha)} \right) $，表明其同时依赖于扩散场能量与相对于波束轴的偏移角。