[发明专利]基于逆边界元法机械噪声远场声压预测方法

摘要

本发明公开了基于逆边界元法机械噪声远场声压预测方法，属于声压预测领域。在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量得到传声器阵列与参考传声器互谱后的场点复声压；基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递矩阵；对传递矩阵进行奇异值分解，得到法向振动速度；根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与法向振动速度的关系py＝ATMyvn，ATMy为对应远场声压的传递矩阵。本发明能够适用于复杂结构表面的，具有高精度的优点。

主权项

基于逆边界元法机械噪声远场声压预测方法，其特征在于：包括以下几个步骤，步骤一：在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量得到传声器阵列与参考传声器互谱后的场点复声压；步骤二：基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递矩阵；步骤三：对传递矩阵进行奇异值分解，得到法向振动速度；步骤四：根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与法向振动速度的关系py＝ATMyvn，ATMy为对应远场声压的传递矩阵；所述的场点复声压与法向振动速度的传递关系为：pf＝ATMvn其中，pf场点声压,vn为法向振动速度,ATM为传递矩阵；所述的法向振动速度为：vn=Σi=1nfiuiTPfσivi]]>其中ui、vi为单位正交矩阵U、V中的列向量，λ为正则化参数，σ1＞σ2＞…＞σnf＞0为传递矩阵的奇异值，对传递矩阵奇异值分解：[ATM]T＝UΣVT其中U、V为单位正交矩阵，即UUT＝I，VVT＝I，正则化参数λ利用广义交叉检验法选取：g(λ)=||[ATM]+xλ-b||22(trace(In-ATM[ATM]+))2]]>其中，[ATM]+为ATM的广义伪逆矩阵，xλ为正则化后得到的正则解，trace()表示矩阵的迹，In单位矩阵，当g(λ)取得最小值对应的λ作为正则化参数。