[发明专利]声场分离方法

摘要

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

主权项

一种声场分离方法，其特征在于包括以下步骤：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处的声压表示为p(r)=Σt=1t∫S′σ(rQ)K(r,rQ)dS′---(5);]]>式中，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进行正则化处理，得到：{pE}=TE(αI+THHTH)-1THHpE(H+)---(8);]]>其中，α为正则化参数，I为单位对角矩阵，为TH的共轭转置矩阵，为的逆矩阵；E、利用步骤D获得的声压求解方程，将两侧声源的声场分离，获取两个声源面的声压值。