[发明专利]一种基于质点振速的部分场分解方法

说明书

技术领域

本发明涉及非完全相干声场的分解技术，尤其是涉及一种基于质点振速的部分场分解方法。

背景技术

近场声全息技术是一种非常有效的噪声源识别和声场可视化技术，但该技术仅能适用于完全相干声场。然而在实际工程中，待研究的声场往往是非完全相干声场，通常需要借助多参考源互谱测量方法将非完全相干声场分解成若干个独立的部分场，使得每个部分场都是完全相干声场，从而采用近场声全息技术对每个部分场分别进行研究。

目前，已有很多文献对部分场分解方法进行了研究，发展了不同的算法，比如基于奇异值分解和主谱能量分析技术的虚相干法，基于偏相干函数的偏相干法，并将其应用扩展至非平稳声场、循环平稳声场以及柱面坐标下的气动声学声场。但是，这些研究都选择声压作为参考。由于声压衰减慢，不能很好地感应和区别不同的潜在声源，导致部分场分解不彻底，一定程度上限制了部分场的分解精度。虽然，有文献选用了声压梯度作为参考，对部分场的分解精度有所提高，但声压梯度是通过两个传声器测量声压，并结合有限差分近似法而获得的，导致该方法的分解精度依赖于两传声器之间的距离，且增加了测量工作和计算难度。

所以，现有的部分场分解方法有两个缺陷和不足：(1)声压衰减慢，限制了基于声压参考的部分场分解方法的分解精度；(2)采用双传声器测量，获得精确的声压梯度成本较高。

发明内容

本发明的目的在于克服上述技术不足，提出一种基于质点振速的部分场分解方法，解决现有技术中基于声压信号的部分场分解方法的分解精度低、分解成本高的技术问题。

为达到上述技术目的，本发明的技术方案提供一种基于质点振速的部分场分解方法，包括如下步骤：

(1)在非完全相干声场的的多个不相干声源附近分别一一对应布置多个质点振速传感器，获取多个质点振速传感器形成的质点振速信号，设定非完全相干声场中有多个测量场点；

(2)对每个质点振速传感器检测的质点振速信号及每个测量场点的场点声压信号进行取样，并将所有的质点振速信号与其中一个测量场点的场点声压信号进行有序排列；

(3)计算有序排列的多个信号的初始条件谱，相继消除有序排列的多个信号中的部分信号的线性作用，并计算迭代条件谱；

(4)计算有序排列的最后两个信号的条件自谱及其之间的条件互谱，并获取偏相干函数；

(5)根据有序排列的最后两个信号的条件互谱和偏相干函数计算有序排列的最后一个质点振速信号相对应的部分场在相对应的测量场点的声压；

(6)将所有的质点振速信号与其中一个测量场点的场点声压信号重新进行有序排列，重复步骤(3)～(5)，计算任一质点振速信号相对应的部分场在任一测量场点的声压。

优选的，所述不相干声源和质点振速传感器均为N个，非完全相干声场中的测量场点为M个，M、N均为自然数，则：

v＝[v₁,…,v_n,…v_N]^T(1)

p＝[p₁,…,p_m,…p_M]^T(2)

其中，v为质点振速信号，v_n为第n个质点振速传感器测量的质点振速信号，T为转置，p为测量场点的场点声压信号，p_m为第m个测量场点的场点声压信号。

优选的，所述步骤(2)中取样次数k，并采用第m个测量场点的场点声压信号进行有序排列，则有序排列顺序为：且其中，是第n个质点振速信号的第k次采样，是第m个场点声压信号的第k次采样。

优选的，所述步骤(3)中初始条件谱和迭代条件谱的计算公式分别为：