[发明专利]采用振速测量和局部近场声全息法的噪声源识别方法

权利要求书

1.一种采用振速测量和局部近场声全息法的噪声源识别方法，其特征是包括如下步骤：

a、测量测量面H上法向质点振速

在声源的辐射声场中，距声源d处为测量面H，在测量面上呈网格式分布有测量点，相邻网格点之间的距离小于半个波长；测量所述测量面上各网格点处的法向质点振速幅值和相位信息获得测量面上的法向质点振速；

b、对测量面H上的法向质点振速补零扩展得到测量面H+上法向质点振速

vE(H+)=vH(H)(x,y)∈H0(x,y)∉H]]>

v_E(H+)＝D·v_H(H+)，其中

D＝diag[D₁₁,…,D_MM]

Dii=1(x,y)∈H0(x,y)∉H]]>

v_E(H+)=[v_E1，…，v_EM]^T

v_H(H+)=[v_H1，…，v_HM]^T

v_H(H)为测量面H上的法向质点振速、

v_E(H+)为测量面H+上由v_H(H)补零后得到的法向质点振速、

v_H(H+)为测量面H+上的法向质点振速、

D为采样算子、D_ii为矩阵D中对角线上的值、

M为测量面H+上测量网格的点数、T为矩阵转置、

（x,y,z）为三维笛卡尔坐标，测量面H（z＝z_H）与（x,y）坐标面平行，测量面的法向为z方向；

c、计算测量面H+上的法向质点振速与声源面S上声压和法向质点振速之间的传递矩阵

WD=DBkcF-1GDF]]>

WN=DBkcF-1GNF]]>其中

GD=ejkzd,]]>GN=kzejkzd/ρck,]]>

Bkc=F-1LkcF,]]>d＝z_H-z_S

kz=k2-kx2-ky2kx2+ky2≤k2-jkx2+ky2-k2kx2+ky2>k2,]]>

Lkc=diag[L11,...,LMM],]]>Lii=1kx2+ky2<kc0kx2+ky2>kc;]]>

W_D为测量面H+上法向质点振速与声源面S上法向质点振速之间的传递矩阵、W_N为测量面H+上法向质点振速与声源面S上声压之间的传递矩阵、

G_D为测量面H+上法向质点振速与声源面S上法向质点振速的格林函数、

G_N为测量面H+上法向质点振速与声源面S上声压的格林函数、

F为二维空间傅里叶变换因子、F^-1为二维空间傅里叶变换因子的逆变换、

为波数域内带限算子、d为测量面与声源面间的距离、

ρ为介质密度、c为介质中声速、k为波数、

为低通滤波器、L_ii为矩阵中对角线上的值、

k_c为截止波数；

声源面S（z＝z_S）与（x,y）坐标面平行，声源面的法向为z方向；

d、求解得到声源面S上的声压和法向质点振速

vS=(αI+WDHWD)-1WDHvE(H+),]]>

pS=(αI+WNHWN)-1WNHvE(H+),]]>其中

v_S为声源面S上的法向质点振速、p_S为声源面S上的声压、

α为正则化参数、I为单位对角矩阵、

为W_N的共轭转置矩阵、

(αI+WNHWN)-1]]>为矩阵(αI+WNHWN)]]>的逆矩阵。