用于麦克风阵列球面谐波分解.docVIP

用于麦克风阵列球面谐波分解 　　摘要：随着麦克风阵列在室内声场分析、声场再现、语音通信中的广泛应用，球面谐波分解在阵列信号处理中的地位变得更加重要。针对球面阵列的球面谐波分解是近年发展起来的波束形成新技术，它还可以扩展到线阵和环形阵等其它阵形中。从球面傅里叶变换着手，在球面谐波域内对不同条件下波束形成算法的原理、性能及其在阵列信号处理中可能的应用进行了全面的综述。 关键词：麦克风阵列；球面谐波；波束形成； 波达方向估计 中图分类号：TP319文献标识码：A文章编号：1672-7800（2013）006-0130-04 作者简介：顾小路（1989-），男，山东大学（威海）机电与信息工程学院硕士研究生，研究方向为麦克风阵列信号处理、波束形成；刘若伦（1971-），男，博士，山东大学（威海）机电与信息工程学院副教授，研究方向为音频信号处理、阵列信号处理。 0引言 球面谐波是波动方程在球坐标系下的解中关于角度的函数，它是球面傅里叶变换的基函数，不同阶次的谐波之间相互正交。球面谐波分解在麦克风阵列处理中的应用是基于入射声场的传播与散射，由傅里叶声学原理[1]提出的，这类方法也称为球面谐波域阵列处理。基于球面谐波分解的波束形成就是首先通过球面傅里叶变换将入射信号的声压变换到球面谐波域内，然后对各阶次的谐波进行加权求和得到整个阵列的输出。 Meyer和Elko[2]通过对声场进行球面谐波分解，得到一种简单、灵活、有效的波束形成器结构，输出波束的形状并不会随着方向的改变而变化。Abhayapala和Ward[3]将类似的阵列应用到高阶声场录音和重构中。Abhayapala[4]还将3D球面阵的分解方法扩展到环形阵中，进一步扩展了球面谐波在不同阵形中的应用范围。随着研究的深入，学者们对球面谐波的提取及球面谐波域波束形成进行了更详细的阐述，提出了许多有效的算法。Rafely[5]曾针对远场窄带情况对球面阵的采样和波束形成进行了总结。 本文结合新近出现的算法，从球面谐波的分解原理和提取方法出发，按照空间域波束形成的分类方法[6]，对球面远场窄带波束形成进行分类，并总结了针对近场源和宽带信号提出的新算法。同时也对比分析了球面谐波域中不同阵形波束形成算法的原理和性能，及其在麦克风阵列信号处理中的可能应用，以方便那些从事麦克风阵列应用研究的读者。 1基础理论 球面傅里叶变换是整个算法的基础，它实现了信号从阵元域到球面谐波域的变换，在实际工程应用中，通常利用采样的离散变换来代替连续的变换。本部分将对球面傅里叶变换和空间采样的规则进行介绍。 1.1球面傅里叶变换 1.2空间采样 此外，Vicente[9]提出了两种新的采样方法，将阵元分布在不同的圆环上，每个环可以构成一个子阵，并将圆环按照纬度排列。在保证满足奈奎斯特定理的前提下，增加阵元之间的距离，从而减少阵元的数目。新的方法能够保持较高的方向性和波束响应的对称性。 2球面阵波束形成 2.1远场窄带波束形成 根据Van Veen和Buckly对传统波束形成的分类方法[6]，我们将球面阵远场窄带波束形成也分为三大类：固定波束形成、统计最优波束形成和多约束的最优波束形成。 2.1.1固定波束形成 2.1.2统计最优波束形成 在球面谐波域内，统计最优波束形成与空间域类似，权值依赖于阵列接收数据的统计特性。其实，阵列的输出不仅可以表示成球面谐波域的采样结果，也可以直接在空间域中表示，具有类似的表达形式。如式（6）所示。 因此，可以将空间域中的统计最优算法扩展到球面谐波域中使用，例如，MVDR、LCMV、GSC等经典算法，并且在球面谐波域中可以减少运算量，提高运行效率[10]。 2.1.3多约束最优波束形成 多约束最优波束形成是指在保证期望波束的同时，可以对其它方向的信号或者相关的性能指标进行约束，最后得到一个最优化的权值，达到对波束进行优化的目的。 Li和Duraswami[11]、Meyer和Elko[12]提出了权值优化的球面波束形成算法，引入白噪声增益作为性能约束条件，权值w是整个优化过程的结果，在一定条件下是最优的。Rafely[13]提出了能够同时在空间域和球面谐波域中实现的多零陷算法，针对非对称的普通波束，在保证期望波束的同时可以对多个方向的干扰进行抑制，增强期望信号，并应用到室内脉冲响应的测量中。 Shefeng Yan[17]将二阶锥化与波束形成相结合，在球面谐波域内将多约束的最优算法看成是凸优化的二阶锥化过程，权值和输出的波束是对多项性能指标优化后得到的最优结果。例如，鲁棒的最小旁瓣算法[14]，在保证期望方向的响应和主瓣宽度的同时，使得旁瓣的峰值最小