声场分离方法技术

本发明专利技术旨在提供一种声场分离方法，包括以下步骤：A、测量测量面上的声压值B、根据振动体产生的声压值生成灰度图，根据实际声源点的位置计算出各点的灰度值，然后将各点灰度值作为各自的权值计算质心球虚拟源强球心坐标及非质心球虚拟源强球心坐标；C、对测量面上测定的声压值进行补零扩展，得到测量面扩展面上的声压值；D、建立测量面扩展面、测量面上声压值与多虚拟球等效源强声压之间的传递矩阵，获得测量面上的声压求解方程；E、结合、声压求解方程，获取两个声源面的声压值。该声场分离方法克服现有技术中波叠加法必须要求虚拟源强配置域与测量面共形的缺陷，同时还具有较高的计算效率和计算精度。

【技术实现步骤摘要】
声场分离方法
本专利技术涉及噪声领域，具体涉及一种声场分离方法。
技术介绍
在实际测量时，通常会遇到测量面两侧都有声源，或是测量面的一侧存在反射或散射的情况。这些干扰声源所产生的声场影响了对目标声源所辐射声场的准确测量。在实际工程中，为了更加准确的研究目标声源的声辐射特性或反射面的反射特性，需要将来自测量面两侧的辐射声分开。现有分离方法包括：(1)基于Fourier变换的单面声场分离技术(《物理学报》(2009年58卷12期))。该方法能够去除奇异性，可以采用单面测量进行声场分离，该方法的缺点是受到傅里叶变换算法的影响，测量面孔径至少是重构面大小的数倍大时才能得到准确分离结果，因此测量量较大，无法对大型声源进行分离，即便能够得到声场分离结果，仍需要继续对声源面进行声场重构；(2)基于双面振速测量的声场分离技术(《声学学报》(2010年35卷6期))。该方法首先测量两个平行等间距测量面上的法向质点振速，再采用傅里叶法分离入射和辐射声场，该方法能够获得更高的法向质点振速精度，但却仍然受到测量孔径大小的限制；(3)基于声压和速度测量的统计最优声场分离方法。Jacobsen等人在《J.Acoust.Soc.Am》(2007年121卷3期)的文章中提出基于声压和速度测量的统计最优声场分离方法,该方法采用p-u声强探头，在单测量面上对局部声压和质点振速信息进行测量，再采用建立的联合求解公式实现来自测量面两侧的辐射声场分离和声源面的重建，该方法解决了测量空间大小的限制，但其缺点是计算量大、效率低、计算时间长。(4)基于球面波叠加的声场分离方法：1956年，J.Pachner采用球面波叠加法实现了任意波场中行波和驻波声场的分离，G.Weinreich等在1980年对J.Pachner提出的方法作了进一步改进，建立了基于双球面测量的声场分离理论。现有声场分离技术仍然存在2个缺点：(1)波叠加法必须要求虚拟源强配置域与测量面共形；(2)计算量大、效率低、计算时间长，适应性不强。
技术实现思路
本专利技术旨在提供一种声场分离方法，该声场分离方法克服现有技术中波叠加法必须要求虚拟源强配置域与测量面共形的缺陷，同时还具有较高的计算效率和计算精度。本专利技术的技术方案如下：一种声场分离方法，包括以下步骤：A、在测量声场中两个声源面之间设置有测量面，测量面上呈网格设置声压振速传感器，相邻网格点之间的距离小于半个波长，声压振速传感器测量测量面上的声压值；B、对振动体形状进行有限元仿真，建立多虚拟球等效源强，将多虚拟球等效源强分为质心球虚拟源强和非质心球虚拟源强两类，将声压图像的灰度值作为加权重心法的权重值，利用测量出的声压值生成灰度图，根据振动体上实际声源的位置坐标计算出各点的灰度值，然后将灰度值作为权值求系统重心，将系统重心作为质心球虚拟源强球心坐标，同时任选两个以上实际声源，利用各个实际声源的灰度值作为权值计算出与这些实际声源对应的非质心球虚拟源强球心坐标；C、对测量面上测定的声压值进行补零扩展，得到测量面扩展面上的声压值；D、建立测量面扩展面声压值与多虚拟球等效源强声压之间的传递矩阵，建立测量面上的声压值与多虚拟球等效源强声压之间的传递矩阵，对测量面扩展面声压值与多虚拟球等效源强声压之间的传递矩阵进行正则化处理，获得测量面上的声压求解方程，求解测量面上的声压值；E、利用步骤D获得的声压求解方程，将两侧声源的声场分离，获取两个声源面的声压值。所述步骤B具体步骤如下：B1、对振动体形状进行有限元仿真，根据振动体产生的声压值生成灰度图，根据实际声源点的位置计算出各个声源的灰度值f(xi,yi,zi)，(xi,yi,zi)为各实际声源点的坐标，设置的多虚拟球等效源强分为质心球虚拟源强和非质心球虚拟源强两类；B2、利用式(1)算出整个系统的重心作为质心球虚拟源强球心坐标，其中v为不规则振动体内部声源的个数；B3、从振动体声源中任选w个声源，参照式(2)计算出一定数量的非质心球虚拟源强球心坐标其中2≤w<v；所述步骤C具体步骤如下：依照式(3)对测量面上的声压值进行补零扩展得到测量面扩展面上的声压值：pE(H+)＝D·pE(H)(3)；其中，pE(H)为测量面H上的声压测量值；pE(H+)为测量面的扩展面上的声压值；D＝diag[D11,D22,…,DNN]，D为采样算子，Dii为矩阵D上对角线上的值；N为测量面扩展面上网格的点数；测量面H(x＝xH)与(y,z)坐标面平行，测量面的法向为x方向。所述步骤D具体步骤如下：D1、建立测量面扩展面声压值与虚拟球等效源强声压值之间的传递矩阵如式(4)：pE(H+)＝THQ(4)；其中pE(H+)为测量面扩展面上的声压值，Q为虚拟球等效源强双向傅里叶分解后的系数矩阵；TH为虚拟球等效源强与测量面扩展面上声压之间的传递矩阵；其中TH由以下步骤得到：根据等效源强理论，假设S′是振动体内某一虚拟源强分布表面，则可将外场中某点r处的声压表示为式中，rQ是虚源面上的某一点，t为多虚拟球的个数，σ(rQ)为待求的源强密度函数，K(r,rQ)为积分核函数，K(r,rQ)＝g(r,rQ)＝(1/4πR)eikR，k＝ω/c为波数,R为两点间的距离；将未知源强密度函数进行双向Fourier级数展开，并利用二维快速FFT计算积分，采用梯形公式离散格林函数，将与距离有关的矩阵规整成一个矩阵，即为TH；D2、建立测量面声压值与虚拟球等效源强声压值之间的传递矩阵如下：{pE}＝TEQ(6)；其中，{pE}为测量面H上的声压值，TE为虚拟球等效源强与测量面上声压之间的传递矩阵，TE的计算参照TH的计算进行；Q为虚拟球等效源强双向傅里叶分解后的系数矩阵；D3、联立式(4)、(6)，求解测量面上的声压：{pE}＝TETH-1pE(H+)(7)；对TH进行正则化处理，得到：其中，α为正则化参数，I为单位对角矩阵，为TH的共轭转置矩阵，为的逆矩阵。所述步骤E具体步骤如下：测量面H上的声压和法向质点振速为：pE(H)＝p1E(H)+p2E(H)；vE(H)＝v1E(H)-v2E(H)(9)；式中p1E(H)为声源1在测量面H上所辐射的声压、p2E(H)为声源2在测量面H上所辐射的声压、v1E(H)为声源1在测量面H上所辐射的法向质点振速、v2E(H)为声源2在测量面H上所辐射的法向质点振速；对式(9)进行补零扩展得到：pE(H+)＝p1E(H+)+p2E(H+)；vE(H+)＝v1E(H+)-v2E(H+)(10)；与式(8)相结合，得到其中p1为声源面S1上的声压、p2为声源面S2上的声压,WE为虚拟球等效源强与测量面上质点振速之间的传递矩阵，WH为虚拟球等效源强与测量面扩展面上质点振速之间的传递矩阵；将式(11)和式(12)相加，将式(13)和式(14)相减，可得：将式17和式18相加，可得：将式17和式18相减，可得：本专利技术技术方案采用波叠加法进行声场分离，并采用灰度重心法对虚拟源强进行配置，灰度重心法通过将各实际声源点的灰度值作为权值来计算相应的虚拟源坐标，不仅克服了常规波叠加法必须要求虚拟源强配置域与测量面共形的缺陷，并且通过灰度权值进行配置的虚拟源强更贴近实际声源点的分布情况，提升声场的重建与分离的准确度；同时本文档来自技高网...

【技术保护点】
一种声场分离方法，其特征在于包括以下步骤：A、在测量声场中两个声源面之间设置有测量面，测量面上呈网格设置声压振速传感器，相邻网格点之间的距离小于半个波长，声压振速传感器测量测量面上的声压值；B、对振动体形状进行有限元仿真，建立多球形虚拟源，将多球形虚拟源分为质心球虚拟源强和非质心球虚拟源强两类，将声压图像的灰度值作为加权重心法的权重值，利用测量出的声压值生成灰度图，根据振动体上实际声源的位置坐标计算出各点的灰度值，然后将灰度值作为权值求系统重心，将系统重心作为质心球虚拟源强球心坐标，同时任选两个以上实际声源，利用各个实际声源的灰度值作为权值计算出与这些实际声源对应的非质心球虚拟源强球心坐标；C、对测量面上测定的声压值进行补零扩展，得到测量面扩展面上的声压值；D、建立测量面扩展面声压值与多虚拟球等效源强声压之间的传递矩阵，建立测量面上的声压值与多虚拟球等效源强声压之间的传递矩阵，对测量面扩展面声压值与多虚拟球等效源强声压之间的传递矩阵进行正则化处理，获得测量面上的声压求解方程，求解测量面上的声压值；E、利用步骤D获得的声压求解方程，将两侧声源的声场分离，获取两个声源面的声压值。

【技术特征摘要】
1.一种声场分离方法，其特征在于包括以下步骤：A、在测量声场中两个声源面之间设置有测量面，测量面上呈网格设置声压振速传感器，相邻网格点之间的距离小于半个波长，声压振速传感器测量测量面上的声压值；B、对振动体形状进行有限元仿真，建立多虚拟球等效源强，将多虚拟球等效源强分为质心球虚拟源强和非质心球虚拟源强两类，将声压图像的灰度值作为加权重心法的权重值，利用测量出的声压值生成灰度图，根据振动体上实际声源的位置坐标计算出各点的灰度值，然后将灰度值作为权值求系统重心，将系统重心作为质心球虚拟源强球心坐标，同时任选两个以上实际声源，利用各个实际声源的灰度值作为权值计算出与这些实际声源对应的非质心球虚拟源强球心坐标；C、对测量面上测定的声压值进行补零扩展，得到测量面扩展面上的声压值；所述步骤C具体步骤如下：依照式(3)对测量面上的声压值进行补零扩展得到测量面扩展面上的声压值：pE(H+)＝D·pE(H)(3)；其中，pE(H)为测量面H上的声压测量值；pE(H+)为测量面的扩展面上的声压值；D＝diag[D11,D22,…,DNN]，D为采样算子，Dii为矩阵D上对角线上的值；N为测量面扩展面上网格的点数；测量面H(x＝xH)与(y,z)坐标面平行，测量面的法向为x方向；D、建立测量面扩展面声压值与多虚拟球等效源强声压之间的传递矩阵，建立测量面上的声压值与多虚拟球等效源强声压之间的传递矩阵，对测量面扩展面声压值与多虚拟球等效源强声压之间的传递矩阵进行正则化处理，获得测量面上的声压求解方程，求解测量面上的声压值；所述步骤D具体步骤如下：D1、建立测量面扩展面声压值与虚拟球等效源强声压值之间的传递矩阵如式(4)：pE(H+)＝THQ(4)；其中pE(H+)为测量面扩展面上的声压值，Q为虚拟球等效源强双向傅里叶分解后的系数矩阵；TH为虚拟球等效源强与测量面扩展面上声压之间的传递矩阵；其中TH由以下步骤得到：根据等效源强理论，假设S′是振动体内某一虚拟源强分布表面，则可将外场中某点r处的声压表示为式中，rQ是虚源面上的某一点，t为多虚拟球的个数，σ(rQ)为待求的源强密度函数，K(r,rQ)为积分核函数，K(r,rQ)＝g(r,rQ)＝(1/4πR)eikR，k＝ω/c为波数,R为两点间的距离；将未知源强密度函数进行双向Fourier级数展开，并利用二维快速FFT计算积分，采用梯形公式离散格林函数，将与距离有关的矩阵规整成一个矩阵，即为TH；D2、建立测量面声压值与虚拟球等效源强声压值之间的传递矩阵如下：{pE}＝TEQ(6)；其中，{pE}为测量面H上的声压值，TE为虚拟球等效源强与测量面上声压之间的传递矩阵，TE的计算参照TH的计算进行；Q为虚拟球等效源强双向傅里叶...

【专利技术属性】
技术研发人员：王玉江，向宇，
申请(专利权)人：广西科技大学，
类型：发明
国别省市：广西;45