CN104596636B - 声场分离方法 - Google Patents

声场分离方法 Download PDF

Info

Publication number
CN104596636B
CN104596636B CN201410836637.4A CN201410836637A CN104596636B CN 104596636 B CN104596636 B CN 104596636B CN 201410836637 A CN201410836637 A CN 201410836637A CN 104596636 B CN104596636 B CN 104596636B
Authority
CN
China
Prior art keywords
sound pressure
sound
virtual
source
measuring surface
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201410836637.4A
Other languages
English (en)
Other versions
CN104596636A (zh
Inventor
王玉江
向宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liuzhou Daoyuan Technology Co Ltd
Original Assignee
Guangxi University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangxi University of Science and Technology filed Critical Guangxi University of Science and Technology
Priority to CN201410836637.4A priority Critical patent/CN104596636B/zh
Publication of CN104596636A publication Critical patent/CN104596636A/zh
Application granted granted Critical
Publication of CN104596636B publication Critical patent/CN104596636B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Golf Clubs (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

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

Description

声场分离方法
技术领域
本发明涉及噪声领域，具体涉及一种声场分离方法。
背景技术
在实际测量时，通常会遇到测量面两侧都有声源，或是测量面的一侧存在反射或散射的情况。这些干扰声源所产生的声场影响了对目标声源所辐射声场的准确测量。在实际工程中，为了更加准确的研究目标声源的声辐射特性或反射面的反射特性，需要将来自测量面两侧的辐射声分开。
现有分离方法包括：(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相减，可得：
本发明技术方案采用波叠加法进行声场分离，并采用灰度重心法对虚拟源强进行配置，灰度重心法通过将各实际声源点的灰度值作为权值来计算相应的虚拟源坐标，不仅克服了常规波叠加法必须要求虚拟源强配置域与测量面共形的缺陷，并且通过灰度权值进行配置的虚拟源强更贴近实际声源点的分布情况，提升声场的重建与分离的准确度；同时，通过将多虚拟球波叠加法与多虚拟域的二维快速Fourier变换算法相结合，极大地改善了计算精度、计算速度和适应性。
附图说明
图1为本发明声场分离方法的步骤流程图
图2为本发明声场分离原理图
图3为多虚拟球等效源强配置面与测量面的布置示意图
图4为实施例1振动体内某截面声压幅值分布图的灰度图
图5为实施例1重建声场声压幅值图
图6为实施例1重建声场声压理论幅值图
图7为实施例1声压幅值相对误差图
图8为实施例1声源1声压幅值图
图9为实施例1声源1声压幅值相对误差图
图10为实施例1声源2声压幅值图
图11为实施例1声源2声压幅值相对误差图
图3中各部分名称及标号如下:
1为测量平面，2为传声器，3为虚拟球等效源。
具体实施方式
下面结合具体实施例说明本发明。
实施例1
如图1所示，本实施例声场分离方法包括以下步骤：
A、在测量声场中两个声源面之间设置有测量面，测量面上呈网格设置声压振速传感器，相邻网格点之间的距离小于半个波长，声压振速传感器测量测量面上的声压值；
B、对振动体形状进行有限元仿真，建立多虚拟球等效源强，将多虚拟球等效源强分为质心球虚拟源强和非质心球虚拟源强两类，将声压图像的灰度值作为加权重心法的权重值，利用测量出的声压值生成灰度图，根据振动体上实际声源的位置坐标计算出各点的灰度值，然后将灰度值作为权值求系统重心，将系统重心作为质心球虚拟源强球心坐标，同时任选两个以上实际声源，利用各个实际声源的灰度值作为权值计算出与这些实际声源对应的非质心球虚拟源强球心坐标；
具体步骤如下：
B1、对振动体形状进行有限元仿真，根据振动体产生的声压值生成灰度图，根据实际声源点的位置计算出各个声源的灰度值f(xi,yi,zi)，(xi,yi,zi)为各实际声源点的坐标，设置的多虚拟球等效源强分为质心球虚拟源强和非质心球虚拟源强两类；
B2、利用式(1)算出整个系统的重心作为质心球虚拟源强球心坐标，
其中v为不规则振动体内部声源的个数；
B3、从振动体声源中任选w个声源，参照式(2)计算出一定数量的非质心球虚拟源强球心坐标其中2≤w<v；
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为虚拟球等效源强双向傅里叶分解后的系数矩阵；
D3、联立式(4)、(6)，求解测量面上的声压：
{pE}＝TETH -1pE(H+) (7)；
对TH进行正则化处理，得到：
其中，α为正则化参数，I为单位对角矩阵，为TH的共轭转置矩阵，为的逆矩阵；
E、利用步骤D获得的声压求解方程，将两侧声源的声场分离，获取两个声源面的声压值；
所述步骤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相减，可得：
如图2所示，本实施例设置位于测量面的两侧为两个振动体S1、S2，其中S1为36个点声源组成的振动体，S2为9个点声源组成的振动体，都分布在0.4mх0.4mх0.2m的长方体内，以测量面的中心为坐标原点，测量面指向声源S2方向为x轴正方向，多虚拟球等效源强配置面与测量面的布置如图3所示；
本实施例采用灰度重心法布置多虚拟球的位置：根据测量出的声压值生成灰度图，振动体某截面声压幅值分布图的灰度图中灰度与该灰度点数量的关系如图4所示，根据实际声源点的位置计算出各点的灰度值，然后将各点灰度值作为各自的权值计算系统重心，进而作为质心球虚拟源强球心坐标；选择25个实际声源点，以其各自灰度值作为各自权值，计算出X个非质心球虚拟源强球心坐标，25个实际声源点灰度值如表1所示；
表1选取的实际声源灰度值统计表
其中：对于质心位置的虚拟球的半径为r0＝0.025m；非质心位置的多虚拟球的半径r0＝0.05m；
本实施例重建的声场声压幅值如图5所示，进行声场分离后的声源1声压幅值、声源2声压幅值分别如图8、10所示；
在自由空间任意一场点P的理论声压表示为：
式中，rn为场点P到第n个点源的距离，Qn为第n个点源的强度；ρ为大气密度，c为声速，k＝2πf/c为波数；
声压幅值、幅值相对误差分别为：
式中pT(r)为各点理论声压幅值，pF(r)为本实施例重建的各点声压幅值；
根据式(20)算出测量面处的理论声压幅值如图6所示；
根据式(21)算出测量面处的理论声压幅值相对误差如图7所示；
根据式(20)、(21)结合声源1声压幅值算出声源1声压幅值相对误差如图9所示；
根据式(20)、(21)结合声源2声压幅值算出声源2声压幅值相对误差如图11所示。
由上述图示可以看出，本实施例的声场分离方法的分离出的声源声压幅值相对误差处于可接受的范围内，验证了本实施例声场分离方法的准确性。

Claims (2)

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处的声压表示为
<mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>t</mi> </munderover> <msub> <mo>&amp;Integral;</mo> <msup> <mi>S</mi> <mo>&amp;prime;</mo> </msup> </msub> <mi>&amp;sigma;</mi> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mi>Q</mi> </msub> <mo>)</mo> </mrow> <mi>K</mi> <mrow> <mo>(</mo> <mi>r</mi> <mo>,</mo> <msub> <mi>r</mi> <mi>Q</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>dS</mi> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式中，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进行正则化处理，得到：
<mrow> <mo>{</mo> <msub> <mi>p</mi> <mi>E</mi> </msub> <mo>}</mo> <mo>=</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>p</mi> <mi>E</mi> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，α为正则化参数，I为单位对角矩阵，为TH的共轭转置矩阵，为的逆矩阵；
E、利用步骤D获得的声压求解方程，将两侧声源的声场分离，获取两个声源面的声压值。
2.如权利要求1所述的声场分离方法，其特征在于所述步骤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)相结合，得到
<mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>p</mi> <mrow> <mn>1</mn> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>p</mi> <mrow> <mn>2</mn> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>W</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>W</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>v</mi> <mrow> <mn>1</mn> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>W</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>W</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>v</mi> <mrow> <mn>2</mn> <mi>E</mi> </mrow> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中p1为声源面S1上的声压、p2为声源面S2上的声压,WE为虚拟球等效源强与测量面上质点振速之间的传递矩阵，WH为虚拟球等效源强与测量面扩展面上质点振速之间的传递矩阵；
将式(11)和式(12)相加，将式(13)和式(14)相减，可得：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mrow> <mn>1</mn> <mi>E</mi> </mrow> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>+</mo> <msub> <mi>p</mi> <mrow> <mn>2</mn> <mi>E</mi> </mrow> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>p</mi> <mi>E</mi> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>W</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>W</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mrow> <mn>1</mn> <mi>E</mi> </mrow> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mn>2</mn> <mi>E</mi> </mrow> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mi>W</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>+</mo> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>W</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>v</mi> <mi>E</mi> </msub> <mrow> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
将式(15)和式(16)相加，可得：
<mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>p</mi> <mi>E</mi> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>+</mo> <msub> <mi>W</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>W</mi> <mi>H</mi> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>v</mi> <mi>E</mi> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
将式(15)和式(16)相减，可得：
<mrow> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>T</mi> <mi>H</mi> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>T</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>p</mi> <mi>E</mi> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>-</mo> <msub> <mi>W</mi> <mi>E</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>&amp;alpha;</mi> <mi>I</mi> <mo>+</mo> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>W</mi> <mi>H</mi> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>W</mi> <mi>H</mi> <mi>H</mi> </msubsup> <msub> <mi>v</mi> <mi>E</mi> </msub> <mo>(</mo> <msup> <mi>H</mi> <mo>+</mo> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
CN201410836637.4A 2014-12-29 2014-12-29 声场分离方法 Expired - Fee Related CN104596636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410836637.4A CN104596636B (zh) 2014-12-29 2014-12-29 声场分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410836637.4A CN104596636B (zh) 2014-12-29 2014-12-29 声场分离方法

Publications (2)

Publication Number Publication Date
CN104596636A CN104596636A (zh) 2015-05-06
CN104596636B true CN104596636B (zh) 2018-03-27

Family

ID=53122567

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410836637.4A Expired - Fee Related CN104596636B (zh) 2014-12-29 2014-12-29 声场分离方法

Country Status (1)

Country Link
CN (1) CN104596636B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106017837B (zh) * 2016-06-30 2018-12-21 北京空间飞行器总体设计部 一种等效声模拟源的模拟方法
EP3571663B1 (en) * 2018-03-29 2020-08-19 Leica Microsystems CMS GmbH Apparatus and method, particularly for microscopes and endoscopes, using baseline estimation and half-quadratic minimization for the deblurring of images
CN111537058B (zh) * 2020-04-16 2022-04-29 哈尔滨工程大学 一种基于Helmholtz方程最小二乘法的声场分离方法
CN112146751B (zh) * 2020-09-27 2022-04-01 合肥工业大学 一种基于时域等效源法的实时声场分离方法
CN112415497B (zh) * 2020-12-11 2023-10-13 中国科学院声学研究所 有限元联合波叠加的浅海波导中散射声场分离方法及系统
CN114414039B (zh) * 2022-03-28 2022-06-14 苏州浪潮智能科技有限公司 一种声源近场中声压测试方法、装置、设备及介质
CN117423350A (zh) * 2023-10-18 2024-01-19 安徽建筑大学 一种基于球形测量面的声场分离系统和方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3986785B2 (ja) * 2001-09-20 2007-10-03 日本放送協会 音源分離収音マイクロホン装置および方法
CN102901559A (zh) * 2012-09-27 2013-01-30 哈尔滨工程大学 一种基于单面测量和局部声全息法的分离声场方法
CN103267571A (zh) * 2013-05-10 2013-08-28 合肥工业大学 采用单面声压和振速测量的实时声场分离方法
CN103592022A (zh) * 2013-11-25 2014-02-19 合肥工业大学 采用声压和质点加速度测量的实时声场分离方法
CN103712684A (zh) * 2013-12-25 2014-04-09 广西科技大学 声场重建方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3986785B2 (ja) * 2001-09-20 2007-10-03 日本放送協会 音源分離収音マイクロホン装置および方法
CN102901559A (zh) * 2012-09-27 2013-01-30 哈尔滨工程大学 一种基于单面测量和局部声全息法的分离声场方法
CN103267571A (zh) * 2013-05-10 2013-08-28 合肥工业大学 采用单面声压和振速测量的实时声场分离方法
CN103592022A (zh) * 2013-11-25 2014-02-19 合肥工业大学 采用声压和质点加速度测量的实时声场分离方法
CN103712684A (zh) * 2013-12-25 2014-04-09 广西科技大学 声场重建方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于多球域波叠加法的PNAH的虚源配置研究;王玉江;《沈阳理工大学学报》;20141215;第33卷(第6期);第54-60页 *

Also Published As

Publication number Publication date
CN104596636A (zh) 2015-05-06

Similar Documents

Publication Publication Date Title
CN104596636B (zh) 声场分离方法
CN104237849B (zh) 一种双五元十字阵被动声定位融合方法
CN107576388A (zh) 一种浅海信道下三维结构声源辐射声场预报方法
CN107478325B (zh) 一种非稳态平面声源的自由场还原方法
CN111157096B (zh) 一种基于局部测量的封闭空间面板声学贡献度识别方法
JP5156934B2 (ja) 音響測定装置
CN108646221A (zh) 一种基于声源阵列的空间麦克风定位方法
CN107576275A (zh) 一种利用摄影测量技术对充气结构进行应变场测量的方法
CN110399680A (zh) 一种浅海弹性结构辐射声场计算方法
CN109238444A (zh) 一种采用稀疏测量的声场分离方法
Shen et al. Periodic boundary based FFT-FISTA for sound source identification
IL273294B1 (en) Metrology method and system
CN106017837B (zh) 一种等效声模拟源的模拟方法
CN105004416B (zh) 基于逆边界元法机械噪声远场声压预测方法
CN104105049A (zh) 一种减少传声器使用数量的房间冲激响应函数测量方法
CN103577696B (zh) 旋转声场作用下非规则缺陷散射声场的计算方法
Yan et al. A 3D temperature field reconstruction algorithm based on acoustic travel-time tomography
KR101469761B1 (ko) 방사 패턴과 총 방사파워 분리계산을 통한 함정의 수중 방사 소음 패턴 산출 방법 및 장치
CN117233266A (zh) 一种基于循环神经网络的全波形反演导波层析成像方法
Martínez et al. SHM via topological derivative
Zhou et al. Impact location estimation in anisotropic structures
Li et al. A low-frequency sound power measurement technique based on the separation of near-field acoustic energy flux in a non-anechoic tank
Liu et al. Acoustic tomography of temperature and wind flow fields in a wind power plant
CN115474965B (zh) 一种超声ct阵列探头阵元指向性自适应评估方法及装置
Havranek et al. Comparison of patch acoustic holography methods for confined space

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200520

Address after: Room 501-1, Jingxuan building, Liuzhou Vocational and technical college, No. 30, Shewan Road, Liuzhou City, 545000 Guangxi Zhuang Autonomous Region

Patentee after: Liuzhou Daoyuan Technology Co., Ltd

Address before: 545006 No. 268 East Ring Road, the Guangxi Zhuang Autonomous Region, Liuzhou

Patentee before: GUANGXI UNIVERSITY OF SCIENCE AND TECHNOLOGY

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180327

Termination date: 20191229