CN110850371A - 一种基于Green函数修正的高分辨率声源定位方法 - Google Patents

Abstract

本发明公开了一种基于Green函数修正的高分辨率声源定位方法，包括：1采集传声器阵列上各个传感器处的声压数据形成测量面声压列向量；2建立测量面声压列向量与聚焦面声源源强列向量之间的Green函数；3通过稀疏重构方法进行求解，获得聚焦面各网格点处的声源源强列向量，并通过其中源强的大小在聚焦面网格上识别与定位出声源的初步位置；4在获得声源的初步位置后，通过一阶泰勒展开，并结合坐标轮换迭代优化方法对Green函数进行修正，利用初步位置处的Green函数迭代估计出实际声源位置对应的修正后的Green函数，从而在聚焦面上的网格点之间实现声源精确位置的亚网格搜索，进而获得声源的高分辨率、高精度定位结果。

Description

一种基于Green函数修正的高分辨率声源定位方法

技术领域

本发明涉及噪声源的识别和定位领域，具体涉及一种基于Green函数修正的高分辨率声源定位方法。

背景技术

噪声源识别技术在机械设备故障诊断、环境噪声治理以及舰船、潜艇减振降噪，高铁、汽车NVH开发等方面有着广泛的应用，并逐步发展出了波束形成、近场声全息、反卷积声源成像等多种先进的噪声源识别技术。近年来，信号处理领域中压缩感知理论的出现为进一步提高噪声源识别精度与空间分辨率，并降低噪声源识别系统的硬件成本提供了理论保证和实现途径，籍此发展出了一些基于压缩感知理论中稀疏重构方法的新兴的噪声源识别技术，例如:基于压缩感知的波束形成算法，基于稀疏等效源和稀疏模态等效源法的近场声全息，基于正交匹配追踪的反卷积声源成像方法等等，这些方法都利用了实际中声源分布通常是空域的先验信息，通过稀疏重构算法求解带有稀疏约束的欠定方程组获得声源分布的稀疏解，实现声源分布的精确重构，从而进一步提高了声源识别的精度与分辨率，降低了测量所需的传感器数目和测量数据量。然而，无论是传统的还是新兴的噪声源识别技术，无论其声源定位精度高低、分辨率如何，其所得到的识别定位结果总是位于聚焦面预先划分好的离散网格点上。这是因为上述方法均假设声源位置恰好落在这些预先划分的离散网格点上，并据此建立相应的Green函数矩阵。然而在实际中声源位置不可能总是正好位于预划分的网格点上，从概率论上来说，声源位置正好在网格点上是小概率事件，其概率接近为0，因此实际中其实很难满足上述假设。实际中对上述声源位置假设的违背，造成算法采用的Green函数与真实的Green函数间存在误差，导致声源识别结果总是落在聚焦面网格点上，使算法无法精确定位位于相邻网格点间的声源，降低了声源识别定位的精度和分辨率。并且由此可见，上述已有的声源识别技术其识别定位精度受聚焦面网格划分密度影响，聚焦面网格越稀疏则算法的识别定位精度越低。当然，加密聚焦面网格可以在一定程度上缓解这个问题，但加密网格同时会导致计算量的成倍增加，影响算法的计算效率和声源识别的实时性。

发明内容

本发明的目的是克服上述现有技术的不足，提出一种基于Green函数修正的高分辨率声源定位方法，以期能实现高精度高分辨率声源识别定位。

为实现上述目的，本发明采用以下技术方案：

本发明一种基于Green函数修正的高分辨率声源定位方法的特点包括以下步骤：

步骤1、在K个声源辐射形成的声场中布置M个传感器，从而形成测量面W；采集各传感器所检测到的声压数据p＝[p(r₁),p(r₂),…p(r_m)]，其中，r_m表示测量面W上第m个测量点的坐标向量，p(r_m)表示第m个测量点的声压数据；

步骤2、定义聚焦面T为离散化声源计算平面所形成的网格面；所述聚焦面T包含N个网格点，且每个网格点为聚焦点；

将每个聚焦点作为潜在声源点，从而利用式(1)构造声源源强列向量q：

q＝[q(r₁),q(r₂),...q(r_n)] (1)

式(1)中，r_n表示聚焦面T上第n个聚焦点的坐标向量，n＝1,2,3,…N；

利用式(2)建立所述聚焦面T上各网格点与传感器阵列中各传感器之间的Green函数矩阵G：

式(2)中，g_m(r_n)为第n个聚焦点到第m个传感器之间的Green函数，且

j为虚数单位，k表示声波波数，|r_n-r_m|为第n个聚焦点到第m个传感器之间的距离，m＝1,2,3，...，M；

步骤3、利用式(3)建立声源源强列向量q与传感器阵列声压测量值p之间的关系：

p＝Gq (3)

步骤4、通过稀疏重构算法对(3)式进行求解，获得所述聚焦面T上各网格点处的声源源强列向量，通过所述声源源强列向量中的源强大小确定所述聚焦面上声源的个数K及各声源的初步位置；

根据各声源的初步位置，从所述Green函数矩阵G中抽取相应的列，从而形成如式(4)所示的待修正的Green函数矩阵G_K：

G_K＝[g(r_c1),(r_c2),…,g(r_cK)] (4)

式(4)中，g(r_cK)为第K个声源的初步位置到传感器阵列中各传感器的Green函数所组成的列向量，即g(r_cK)＝[g₁(r_ck),g₂(r_ck),…,g_m(r_cK)]^T；其中g_m(r_cK)为第K个声源的初步位置到第m个传感器的Green函数；

步骤5、对所述待修正的Green函数矩阵G_K进行一阶泰勒展开，并结合坐标轮换迭代优化方法，从而估计出实际声源位置对应的修正后的Green函数矩阵G_t，从而在所述聚焦面T上网格点之间实现亚网格搜索，得到准确的声源位置。

本发明所述的高分辨率声源定位方法的特点也在于：所述步骤5是按如下过程进行：

步骤5.1：利用式(5)所示的一阶泰勒展开方法，得到修正后的Green函数矩阵G_t：

式(5)中，ΔX，ΔY分别为由各声源在x、y方向的位置偏差所组成的对角矩阵，并有：

式中，Δx_K为第K个声源实际位置相对于初步位置在x方向的位置偏差，Δy_K为第K个声源实际位置相对于初步位置在y方向的位置偏差；

步骤5.2：利用式(6)建立声压数据p与K个声源对应的修正后的Green函数矩阵G_t之间的数学关系：

式(6)中，q_K为K个声源的源强列向量；

步骤5.3：定义当前迭代次数为i，最大迭代次数为Imax；定义当前迭代维度为j，并初始化i＝1；

步骤5.4：令第i轮迭代中第j个维度的K个声源的源强列向量为

并初始化为q_K；

步骤5.5：初始化j＝0；

步骤5.6：先进行x方向偏差修正，利用式(7)计算第i轮迭代中由于声源在x方向存在位置偏差造成的传感器阵列中各传感器处的声压偏差：

式(7)中，

为第i轮迭代中由于声源在x方向存在位置偏差造成的第m个传感器处的声压偏差；ΔYⁱ为第i轮迭代中K个声源的y方向位置偏差组成的对角矩阵；当i＝1时，ΔYⁱ为0；

步骤5.7：利用式(8)计算第i轮迭代中声源在x方向的位置偏差：

式(8)中，

为第i轮迭代中第K个声源在x方向的位置偏差；

步骤5.8：将j+1赋值给j；

步骤5.9：利用式(9)计算第i轮迭代中第j个维度的修正后的Green函数矩阵

式(9)中，ΔXⁱ为第i轮迭代中K个声源的x方向位置偏差组成的对角矩阵；

步骤5.10：利用式(10)计算第i轮迭代中第j个维度的K个声源的源强列向量

步骤5.11：再进行y方向偏差修正，利用式(11)计算第i轮迭代中由于声源在y方向存在位置偏差造成的传感器阵列中各传感器处的声压偏差：

式(11)中，

为第i轮迭代中由于声源在y方向存在位置偏差造成的第m个传感器处的声压偏差；

步骤5.12：利用式(12)计算第i轮迭代中声源在y方向的位置偏差：

式(12)中，

为第i轮迭代中第K个声源在y方向的位置偏差；

步骤5.13：将j+1赋值给j；

步骤5.14：利用式(13)计算第i轮迭代中第j个维度的修正后的Green函数矩阵

步骤5.15：利用式(14)计算第i轮迭代中第j个维度的K个声源的源强列向量

步骤5.16：判断i>Imax是否成立，若成立，则停止迭代，并得到准确的声源位置，否则，将赋值给

后，将i+1赋值给i，同时返回步骤5.5顺序执行。

与现有技术相比，本发明的有益效果在于：

1、本发明通过对Green函数的精确修正，实现了聚焦面上网格点之间的亚网格声源搜索与定位，使得声源识别的结果更加准确，显著提高了实际位置不在聚焦面网格点上的声源的识别定位精度。

2、本发明方法对实际位置在网格点上的声源识别也很精确，因此本方法可用于聚焦面上任意位置声源的识别定位，且定位精度相对于传统方法得到了明显提升。

3、本发明方法在聚焦面网格密集度仅为传统方法网格密集度的几分之一时即可获得优于传统方法的声源定位精度和空间分辨率。因此，在实际中可配合较为稀疏的聚焦面网格点使用，显著降低了计算规模，可实现高效率的声源识别定位。

附图说明

图1为本发明声源识别示意图；

图2为声源频率为1600Hz时，传统稀疏重构声源定位方法和本发明方法对两个实际位置不在预划分网格点上的声源的定位效果图；

图3为声源频率为2000Hz时，传统稀疏重构声源定位方法和本发明方法对两个实际位置不在预划分网格点上的声源的定位效果图；

图4为声源频率为2400Hz时，传统稀疏重构声源定位方法和本发明方法对两个实际位置不在预划分网格点上的声源的定位效果图；

图5a为传统稀疏重构声源定位方法对两个实际位置在预划分网格点上的声源的定位效果图；

图5b为本发明方法对两个实际位置在预划分网格点上的声源的定位效果图；

图6为本发明方法在聚焦面网格被划分为11×11个聚焦点时对两个实际位置不在预划分网格点上的声源的定位效果图；

图7为传统稀疏重构声源定位方法在聚焦面网格被划分为21×21个聚焦点时对两个实际位置不在预划分网格点上的声源的定位效果图；

图8为传统稀疏重构声源定位方法在聚焦面网格被划分为41×41个聚焦点时对两个实际位置不在预划分网格点上的声源的定位效果图；

具体实施方式

本实施例中，一种基于Green函数修正的高分辨率声源定位方法是利用稀疏重构算法所得到的初步定位结果，通过一阶泰勒展开的方式和利用坐标轮换法实现了Green函数的精确修正，从而能识别与定位出实际位置不在聚焦面网格点上的声源的精确位置，实现了声源的亚网格精确定位。该方法不但可以进一步提高声源的定位精度与分辨率，实际应用时还可配合较为稀疏的聚焦面网格点使用，从而降低了计算规模，提高了计算效率，因此该方法具有广泛的运用范围和实际意义。具体的说，如图1所示，是按如下步骤进行：

步骤1、在K个声源辐射形成的声场中布置M个传感器，从而形成测量面W；采集各传感器所检测到的声压数据p＝[p(r₁),p(r₂),…p(r_m)]，其中，r_m表示测量面W上第m个测量点的坐标向量，p(r_m)表示第m个测量点的声压数据；传感器的数量M大于声场中潜在的声源个数K。

步骤2、定义聚焦面T为离散化声源计算平面所形成的网格面；聚焦面T包含N个网格点，且每个网格点为聚焦点；

将每个聚焦点作为潜在声源点，从而利用式(1)构造声源源强列向量q：

q＝[q(r₁),q(r₂),...q(r_n)] (1)

式(1)中，r_n表示聚焦面T上第n个聚焦点的坐标向量，n＝1,2,3,…N；

利用式(2)建立聚焦面T上各网格点与传感器阵列中各传感器之间的Green函数矩阵G：

式(2)中，g_m(r_n)为第n个聚焦点到第m个传感器之间的Green函数，且

j为虚数单位，k表示声波波数，|r_n-r_m|为第n个聚焦点到第m个传感器之间的距离，m＝1,2,3，...，M；

步骤3、利用式(3)建立声源源强列向量q与传感器阵列声压测量值p之间的关系：

p＝Gq (3)

步骤4、通过稀疏重构算法对(3)式进行求解，获得聚焦面T上各网格点处的声源源强列向量，通过声源源强列向量中的源强大小确定聚焦面上声源的个数K及各声源的初步位置；

根据各声源的初步位置，从Green函数矩阵G中抽取相应的列，从而形成如式(4)所示的待修正的Green函数矩阵G_K：

式(4)中，g(r_cK)为第K个声源的初步位置到传感器阵列中各传感器的Green函数所组成的列向量，即g(r_cK)＝[g₁(r_ck),g₂(r_ck),…,g_m(r_cK)]^T；其中g_m(r_cK)为第K个声源的初步位置到第m个传感器的Green函数；

步骤4中，采用稀疏重构算法求得聚焦面上声源数量与初步位置。具体可以采用的稀疏重构算法有很多选择，可以是正交匹配追踪算法、基追踪算法、L1范数最小化算法、稀疏贝叶斯方法等等。其中正交匹配追踪算法具有收敛速度快，计算量小的优点；同时为了便于说明和仿真例计算，本发明仿真验证将采用正交匹配追踪算法作为选用的稀疏重构算法，其他稀疏重构算法在实际中也同样适用。

步骤5、对待修正的Green函数矩阵G_K进行一阶泰勒展开，并结合坐标轮换迭代优化方法，从而估计出实际声源位置对应的修正后的Green函数矩阵G_t，从而在聚焦面T上网格点之间实现亚网格搜索，得到准确的声源位置。

步骤5.1：利用式(5)所示的一阶泰勒展开方法，得到修正后的Green函数矩阵G_t：

式(5)中，ΔX，ΔY分别为由各声源在x、y方向的位置偏差(实际声源位置与声源初步位置的坐标差值)所组成的对角矩阵，并有：

式中，Δx_K为第K个声源实际位置相对于初步位置在x方向的位置偏差，Δy_K为第K个声源实际位置相对于初步位置在y方向的位置偏差；

步骤5.2：利用式(6)建立声压数据p与K个声源对应的修正后的Green函数矩阵G_t之间的数学关系：

式(6)中，q_K为K个声源的源强列向量；

步骤5.3：定义当前迭代次数为i，最大迭代次数为Imax；定义当前迭代维度为j(j＝1代表进行x方向位置修正，j＝2代表进行y方向位置修正)，并初始化i＝1；

步骤5.4：令第i轮迭代中第j个维度的K个声源的源强列向量初始化为q_K；

步骤5.5：初始化j＝0；

步骤5.6：先进行x方向偏差修正，利用式(7)计算第i轮迭代中由于声源在x方向存在位置偏差造成的传感器阵列中各传感器处的声压偏差：

式(7)中，

为第i轮迭代中由于声源在x方向存在位置偏差造成的第m个传感器处的声压偏差；ΔYⁱ为第i轮迭代中K个声源的y方向位置偏差组成的对角矩阵；当i＝1时，ΔYⁱ为0；

步骤5.7：利用式(8)计算第i轮迭代中声源在x方向的位置偏差：

式(8)中，

为第i轮迭代中第K个声源在x方向的位置偏差；

步骤5.8：将j+1赋值给j；

步骤5.9：利用式(9)计算第i轮迭代中第j个维度的修正后的Green函数矩阵

式(9)中，ΔXⁱ为第i轮迭代中K个声源的x方向位置偏差组成的对角矩阵矩阵；

步骤5.10：利用式(10)计算第i轮迭代中第j个维度的K个声源的源强列向量

步骤5.11：再进行y方向偏差修正，利用式(11)计算第i轮迭代中由于声源在y方向存在位置偏差造成的传感器阵列中各传感器处的声压偏差：

式(7)中，

为第i轮迭代中由于声源在y方向存在位置偏差造成的第m个传感器处的声压偏差；

步骤5.12：利用式(12)计算第i轮迭代中声源在y方向的位置偏差：

式(12)中，

为第i轮迭代中第K个声源在x方向的位置偏差；

步骤5.13：将j+1赋值给j；

步骤5.14：利用式(13)计算第i轮迭代中第j个维度的修正后的Green函数矩阵

步骤5.15：利用式(14)计算第i轮迭代中第j个维度的K个声源的源强列向量

步骤5.16：判断i>Imax是否成立，若成立，则停止迭代，并得到准确的声源位置，否则，令

并将i+1赋值给i，同时返回步骤5.5顺序执行。

综上所述，本方法先通过稀疏重构方法定位与识别出声源的初步位置，再通过一阶泰勒展开，结合坐标轮换的迭代优化方法对Green函数进行修正，利用初步位置处的Green函数迭代估计出实际声源位置对应的更加精确的Green函数，从而在在聚焦面上网格点之间实现了声源的亚网格搜索，提高了声源定位识别的效果，使得声源定位更加精确。

为进一步的描述本发明，此处给出如下仿真例：

仿真例1

仿真目的：验证本发明算法相比传统稀疏重构声源定位方法能更好地识别与定位实际位置不在聚焦面预划分的网格点上的声源。

仿真过程：数值仿真中，测量面W、聚焦面T的分布如图1所示，假设空间中有二个等强度的点声源——声源1和声源2，分别位于聚焦面上(-0.23，-0.23)m，(0.02，0.02)m的位置。测量面W位于Z＝0m平面，W的平面尺寸为1m×1m，其上沿x轴、y轴均匀划分11×11个测量点，测量点间距为0.1m。聚焦面T位于Z＝1m平面，T的平面尺寸也为1m×1m，并划分21×21个聚焦网格点，聚焦点间距为0.05m。根据聚焦点的划分方法可知，两个声源均不在聚焦面的网格点上。为使仿真与实际实施中存在测量噪声的情况更加一致，测量面声压添加了高斯白噪声，信噪比为30dB。

图2、图3、图4则分别给出了声源频率为1600Hz、2000Hz、2400Hz时传统稀疏重构声源定位方法和本发明算法对2个实际位置不在预划分网格点上的声源的识别效果图。声源识别结果的图像中‘o’表示声源实际所在位置，‘▽’表示通过传统稀疏重构声源定位方法定位的声源位置，‘+’表示本发明所提方法识别的声源位置。

从图2、图3、图4中可以看出，当2个声源的实际位置不在预划分的网格点上时，无论声源频率为1600Hz，2000Hz还是2400Hz，传统稀疏重构声源定位方法识别到的定位结果总在实际声源位置临近的网格点上，两个声源的识别位置与其实际位置存在偏移，传统稀疏重构声源定位方法识别与定位误差大。为定量的描述声源定位精度，定义声源定位相对误差如下：

根据声源定位相对误差公式，声源频率为1600Hz时，通过传统稀疏重构声源定位方法定位，声源1的x轴和y轴坐标的相对误差分别为8.7％和8.7％，声源2的x轴和y轴坐标的相对误差分别为100％和100％。通过本发明算法进行定位时，声源1的x轴和y轴坐标的相对误差分别为0.65％和0.48％，声源2的x轴和y轴坐标的相对误差分别为3％和13％；声源频率为2000Hz时，通过传统稀疏重构声源定位方法定位，声源1的x轴和y轴坐标的相对误差分别为8.7％和8.7％，声源2的x轴和y轴坐标的相对误差分别为100％和100％。通过本发明算法进行定位时，声源1的x轴和y轴坐标的相对误差分别为0.91％和0.39％，声源2的x轴和y轴坐标的相对误差分别为2％和5％；声源频率为2400Hz时，通过传统稀疏重构声源定位方法定位，声源1的x轴和y轴坐标的相对误差分别为8.7％和8.7％，声源2的x轴和y轴坐标的相对误差分别为100％和100％。通过本发明算法进行定位时，声源1的x轴和y轴坐标的相对误差分别为0.63％和0.91％，声源2的x轴和y轴坐标的相对误差分别为6％和5.5％；本发明算法的声源定位相对误差减小到了传统稀疏重构声源定位方法的十分之一左右甚至更小。由此可以看出，在不同频率下，本发明算法对声源的实际位置不在预划分的网格点上时定位更精确。

仿真例2

仿真目的：验证对于实际位置落在预划分的网格点上的声源，本发明算法和传统稀疏重构声源定位方法有着相同的识别定位精度。

仿真过程：数值仿真中，测量面W、聚焦面T的分布如图1所示，假设空间中有二个等强度的点声源——声源1和声源2，分别位于聚焦面上(-0.25，-0.25)m，(0.05，0.05)m的位置，声源频率为2000Hz。测量面W位于Z＝0m平面，W的平面尺寸为1m×1m，其上沿x轴、y轴均匀划分11×11个测量点，测量点间距为0.1m。聚焦面T位于Z＝1m平面，T的平面尺寸也为1m×1m，并划分21×21个聚焦网格点，聚焦点间距为0.05m。根据聚焦点的划分方法可知，两个声源均在聚焦面的网格点上。

图5a、图5b则分别展示了传统稀疏重构声源定位方法和本发明算法对2个实际位置在预划分网格点上的声源的识别效果图。声源识别结果的图像中‘o’表示声源实际所在位置，‘▽’表示通过传统稀疏重构声源定位方法定位的声源位置，‘+’表示本发明所提方法所识别的声源的位置。通过对比图5a和图5b，可以看出两图都能准确识别出声源位置，从而得到结论：当声源位置在离散网格点上时，两种方法都具有很好的性能，均能精确的识别与定位声源。

仿真例3

仿真目的：验证对于实际位置不在预划分的网格点上的声源，本发明算法可以在较稀疏的聚焦面网格上获得比传统稀疏重构声源定位方法在更密集的网格上更好的定位精度和更高的计算效率。

仿真过程：数值仿真中，测量面W、聚焦面T的分布如图1所示，假设空间中有二个等强度的点声源——声源1和声源2，分别位于聚焦面上(-0.23，-0.23)m，(0.02，0.02)m的位置，声源频率为2000Hz。测量面W位于Z＝0m平面，W的平面尺寸为1m×1m，其上沿x轴、y轴均匀划分11×11个测量点，测量点间距为0.1m。聚焦面T位于Z＝1m平面，T的平面尺寸也为1m×1m，聚焦面分别划分为11×11，21×21，41×41个聚焦网格点，聚焦点间距分别为0.1m，0.05m，0.025m。根据聚焦点的划分方法可知，两个声源均不在聚焦面的网格点上。

图6为本发明算法在更稀疏的网格上(聚焦面被划分为11×11个聚焦点)对两个实际位置不在预划分网格点上的声源的声源定位效果图；图7为传统稀疏重构声源定位方法在密集一倍的网格上(聚焦面被划分为21×21个聚焦点)对两个实际位置不在预划分网格点上的声源的声源定位效果图。图8为传统稀疏重构声源定位方法在更密集的网格上(聚焦面被划分为41×41个聚焦点)对两个实际位置不在预划分网格点上的声源的声源定位效果图。声源识别结果的图像中‘o’表示声源实际所在位置，‘▽’表示通过传统稀疏重构声源定位方法定位的声源位置，‘+’表示本发明所提方法所识别的声源的位置。可以看出，本发明算法在更稀疏的网格上定位到的声源位置比传统稀疏重构声源定位方法在更密集的网格上还准确；本发明算法在聚焦面被划分为11×11个聚焦点的网格上定位时，声源1的x轴和y轴坐标的相对误差分别为1.69％和1.26％，声源2的x轴和y轴坐标的相对误差分别为10％，10％，计算时间为0.508s；传统稀疏重构声源定位方法在聚焦面被划分为21×21个聚焦点的网格上定位时，声源1的x轴和y轴坐标的相对误差分别为8.7％和8.7％，声源2的x轴和y轴坐标的相对误差分别为100％，100％，计算时间为0.561s；传统稀疏重构声源定位方法在聚焦面被划分为41×41个聚焦点的网格上定位时，声源1的x轴和y轴坐标的相对误差分别为2.17％和2.17％，声源2的x轴和y轴坐标的相对误差分别为25％，25％，计算时间为1.101s；表1为本发明算法与传统稀疏重构声源定位方法在稀疏程度不同的聚焦面网格上进行声源定位时的相对误差和计算时间，所使用的同一服务器的处理器CPU为Intel Corei5-7200U 2.50GHz、内存RAM为8.0GB。

表1

从表1中可以看出本发明算法采用11×11点的聚焦面网格获得的定位误差甚至比传统稀疏重构方法采用41×41点的聚焦面网格还要小很多，这充分说明了本方法在定位精度与分辨率方面的巨大优势。同时在计算效率上，由于本发明算法可采用更稀疏的聚焦面网格，因此其计算规模远小于传统稀疏重构方法，所需的计算时间也就更少。因此，本方法相比于传统稀疏重构声源定位方法，在识别精度和计算效率上都有明显的提升。

Claims (2)

1.一种基于Green函数修正的高分辨率声源定位方法，其特征包括以下步骤：

步骤1、在K个声源辐射形成的声场中布置M个传感器，从而形成测量面W；采集各传感器所检测到的声压数据p＝[p(r₁),p(r₂),…p(r_m)]，其中，r_m表示测量面W上第m个测量点的坐标向量，p(r_m)表示第m个测量点的声压数据；

步骤2、定义聚焦面T为离散化声源计算平面所形成的网格面；所述聚焦面T包含N个网格点，且每个网格点为聚焦点；

将每个聚焦点作为潜在声源点，从而利用式(1)构造声源源强列向量q：

q＝[q(r₁),q(r₂),...q(r_n)] (1)

式(1)中，r_n表示聚焦面T上第n个聚焦点的坐标向量，n＝1,2,3,…N；

利用式(2)建立所述聚焦面T上各网格点与传感器阵列中各传感器之间的Green函数矩阵G：

式(2)中，g_m(r_n)为第n个聚焦点到第m个传感器之间的Green函数，且

j为虚数单位，k表示声波波数，|r_n-r_m|为第n个聚焦点到第m个传感器之间的距离，m＝1,2,3，...，M；

步骤3、利用式(3)建立声源源强列向量q与传感器阵列声压测量值p之间的关系：

p＝Gq (3)

步骤4、通过稀疏重构算法对(3)式进行求解，获得所述聚焦面T上各网格点处的声源源强列向量，通过所述声源源强列向量中的源强大小确定所述聚焦面上声源的个数K及各声源的初步位置；

根据各声源的初步位置，从所述Green函数矩阵G中抽取相应的列，从而形成如式(4)所示的待修正的Green函数矩阵G_K：

G_K＝[g(r_c1),(r_c2),…,g(r_cK)] (4)

式(4)中，g(r_cK)为第K个声源的初步位置到传感器阵列中各传感器的Green函数所组成的列向量，即g(r_cK)＝[g₁(r_ck),g₂(r_ck),…,g_m(r_cK)]^T；其中g_m(r_cK)为第K个声源的初步位置到第m个传感器的Green函数；

步骤5、对所述待修正的Green函数矩阵G_K进行一阶泰勒展开，并结合坐标轮换迭代优化方法，从而估计出实际声源位置对应的修正后的Green函数矩阵G_t，从而在所述聚焦面T上网格点之间实现亚网格搜索，得到准确的声源位置。

2.根据权利要求1中所述的高分辨率声源定位方法，其特征在于：所述步骤5是按如下过程进行：

步骤5.1：利用式(5)所示的一阶泰勒展开方法，得到修正后的Green函数矩阵G_t：

式(5)中，ΔX，ΔY分别为由各声源在x、y方向的位置偏差所组成的对角矩阵，并有：

式中，Δx_K为第K个声源实际位置相对于初步位置在x方向的位置偏差，Δy_K为第K个声源实际位置相对于初步位置在y方向的位置偏差；

步骤5.2：利用式(6)建立声压数据p与K个声源对应的修正后的Green函数矩阵G_t之间的数学关系：

式(6)中，q_K为K个声源的源强列向量；

步骤5.3：定义当前迭代次数为i，最大迭代次数为Imax；定义当前迭代维度为j，并初始化i＝1；

步骤5.4：令第i轮迭代中第j个维度的K个声源的源强列向量为

并初始化为q_K；

步骤5.5：初始化j＝0；

步骤5.6：先进行x方向偏差修正，利用式(7)计算第i轮迭代中由于声源在x方向存在位置偏差造成的传感器阵列中各传感器处的声压偏差：

式(7)中，

为第i轮迭代中由于声源在x方向存在位置偏差造成的第m个传感器处的声压偏差；ΔYⁱ为第i轮迭代中K个声源的y方向位置偏差组成的对角矩阵；当i＝1时，ΔYⁱ为0；

步骤5.7：利用式(8)计算第i轮迭代中声源在x方向的位置偏差：

式(8)中，

为第i轮迭代中第K个声源在x方向的位置偏差；

步骤5.8：将j+1赋值给j；

步骤5.9：利用式(9)计算第i轮迭代中第j个维度的修正后的Green函数矩阵

式(9)中，ΔXⁱ为第i轮迭代中K个声源的x方向位置偏差组成的对角矩阵；

步骤5.10：利用式(10)计算第i轮迭代中第j个维度的K个声源的源强列向量

步骤5.11：再进行y方向偏差修正，利用式(11)计算第i轮迭代中由于声源在y方向存在位置偏差造成的传感器阵列中各传感器处的声压偏差：

式(11)中，

为第i轮迭代中由于声源在y方向存在位置偏差造成的第m个传感器处的声压偏差；

步骤5.12：利用式(12)计算第i轮迭代中声源在y方向的位置偏差：

式(12)中，为第i轮迭代中第K个声源在y方向的位置偏差；

步骤5.13：将j+1赋值给j；

步骤5.14：利用式(13)计算第i轮迭代中第j个维度的修正后的Green函数矩阵

步骤5.15：利用式(14)计算第i轮迭代中第j个维度的K个声源的源强列向量

步骤5.16：判断i>Imax是否成立，若成立，则停止迭代，并得到准确的声源位置，否则，将

赋值给

后，将i+1赋值给i，同时返回步骤5.5顺序执行。

Images