采用边界元法重建循环平稳声源的方法技术

本发明专利技术涉及的是一种信号重建技术领域的采用边界元法重建循环平稳声源的方法。本发明专利技术通过设置待测声源、传声器阵列和参考源传声器，同步记录采集到的参考信号和全息测点声压，通过从参考信号中获取参考相位，利用参考信号与全息测点声压之间的相位关系，得到传声器阵列所采集到的全息测点声压的相对相位关系最终由边界元法重建计算获得待测声源信号的自谱相关密度向量。本发明专利技术通过对用传声器阵列采集得到的全息测量面上声源信号的重建，可以分析任意外形的声源，得到声源的三维谱相关密度分布，与传统的采用近场声全息技术辨识声源的方法相比较，本发明专利技术能够适用于循环平稳声源，同时突破了对声源外形的要求。

【技术实现步骤摘要】

本专利技术涉及的是一种信号重建
的方法，特别是一种采用边界元法 重建循环平稳声源的方法。技术背景为了有效地控制噪声，在降噪措施实施以前，必须首先进行噪声源诊断， 确定各主要噪声源地位置及其特性。随着现代信号处理技术地发展，谱分析技术， 相干和偏相干分析技术、声强分析技术以及声近场全息技术等得到了迅速地发 展。循环平稳信号是一类特殊的非平稳信号，由于自身独特的周期平稳性，使 得单次采集到的记录具有周期遍历特性，增加了分析这类非平稳信号方法。循环平稳信号在工程应用中有着十分重要的现实意义，例如旋转机械设备由于对称或 近似对称的物理结构和周期性的工作运动模式，其声源有明显的周期时变特征， 声源信号具有循环平稳特性。经对现有技术文献的检索发现，Veronesi W. ， Maynard J.在《The Journal of the Acoustical Society of America》(85(1989): 588-598)上撰文"Digital holographic reconstruction of sources with arbitrarily shaped surfaces"(《美国声 学学会杂志》适用于任意外形声源的数值全息重建)，该文提出利用常值边界元 法实现对任意外形声源的近场声全息重建，该方法可以有效地分析噪声源的辐射 特性，预测外部声辐射。后来，围绕着该项技术，展开了许多研究，包括应用环 境，算法改进，精度提高，测量系统等。但是，所有工作都是基于平稳声场进行 的，因此，有必要提出新的技术，用于非平稳声场的分析。然而，对于一般的非 平稳声场，声信号的统计特性参量是时变的，因而也就无法再用时间平均来代替 集合平均，使得数据采集十分困难，很难分析声场的特性。现有技术文献中的中国专利公开号1487500 ，名称为"采用近场声全息 技术辨识非平稳声源的方法"，该技术对于平面的循环平稳声源进行过分析，但是不能应用于任意外形。在进一步检索中，尚未发现与本专利技术主题相同或者相 类似的文献报道。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种采用边界元法重建循环平 稳声源的方法，本专利技术通过循环平稳声源的近场全息重建，由全息测量面上测量 到的声压信号辨识出噪声源的特性和三维声源分布，适用于任意外形的循环平稳 声源。本专利技术通过以下技术方案实现，本专利技术在多个声源产生的复杂循环平稳声源 的场合，采用循环平稳理论替代传统的傅立叶变换技术，选择谱相关密度函数作 为声源重建的物理量，提出了用于分析循环平稳声源的近场全息技术，首先布置 参考传声器，提取参考信号；再布置传声器阵列对全息测量面进行扫描测量，采 集全息测量面数据；然后，利用基于边界元法的循环平稳近场全息技术，重建得 到声源的三维分布。本专利技术具体包括以下步骤第一步、设置待测声源表面为&，源面&可以是任意形状，全息测量面^ 与待测声源表面&共形，但为测量方便也可以是&附近的一个任意面。第二步、布置传声器阵列，在全息测量面&上采用线形均布方式设置传声 器阵列，以便采集待测声源的高频信息；另外在声源附近布置一个相对于声源位 置固定不变的参考源传声器用以采集参考信号；传声器阵列所含传声器总数为 0。第三步、同步记录参考源传声器和传声器列阵采集到的参考信号和全息测 点声压，并用空间定位仪确定传声器列阵中每个传声器的空间位置， 一并存储在 存储设备中。第四步、从参考信号中获取参考相位，利用参考信号与全息测点声压之间 的相位关系，可以得到传声器阵列所采集到的全息测点声压的相对相位关系。具 体如下首先分析全息测点声压和参考信号，参考信号的谱相关密度函数定义为-其中/为声源特性的频率，上标a为循环频率，下标r代表参考信号，at代 表参考信号的谱相关密度函数，即自谱相关密度函数。然后计算出参考信号和全息测点声压的互谱相关密度向量fe"，/)L,及fe"，/)L',其中^表示传声器阵列上各个传声器的空间位置，0为传声器阵列所含传 声器总数，下标r代表参考信号，下标^代表全息测量信号，^代表参考信号谱 分量(/ + /2)与全息测点声压谱分量(/-"/2)的互谱相关密度函数，f 代表全息测点声压谱分量(/ + "/2)与参考信号谱分量(/-a/2)的互谱相 关密度函数；结合参考信号的谱相关密度函数《(/)，获得全息测点声压的自谱相关密度 向量fe"，/))exi ，其中下标/^代表全息测点声压的谱相关密度函数； 最后得到式(1):fc(/)L,老(/)K(/)W W) (1)其中" "表示点乘，该公式反映了参考信号和全息测点声压的时域相位关 系，用于进行相位同步。第五步、由边界元法重建计算获得待测声源信号的自谱相关密度向量。 设定待测声源表面&上iV个节点发出的待测声源信号的自谱相关密度向量为(s:(。，/))潟，其中^表示待测声源表面&上n个节点各自的空间位置，iV为自然常数；设定参考信号和待测声源信号的互谱相关密度向量为fc"，/)L及fe(。，/)L,由于传声器列阵上各个全息测点声压的空间相位关系的损失，s^(^,/)无法 直接由s;(&，/)重建。但是，s;(^/)和s;(^/)仍然保留着空间相位关系，可以直接用于重建。设定r(/)是参考信号的频谱分量，p(/)为全息测点声压的频谱分量，v(/)为 待测声源表面&上JV个节点法向振速的频谱分量，则88^(^/)为吖/ + "/2)和？'(/-a/2)的函数； S;(^/)为P(Z + "/2)和a/2)的函数；其中S^(^/)及S;(:，/)为参考信号和全息测点声压的互谱相关密度向量；<formula>formula see original document page 9</formula>的函数；其中上标*为共轭转置，K(^/)及^^(^,/)为参考信号和表面节点法向速度的互谱相关密度向量，s^^,/)和^(^,/)应该由s)(^,/)在频率(/-"/2)上重建得到，而S;(r、.，/)和^k，/)应该由S;(^/)在频率 (/ + /2)上重建得到。因此，可以得到如下以互谱相关密度作为变量的Helmholtz公式(2)-(3):<formula>formula see original document page 9</formula>其中G(rA，^，/-a/2)和G(:，r,，Z + "/2)是Green函数，p表示介 质密度，C是和测点所处位置相关的系数。利用边界元法，将式(2)和(3)离散为矩阵形式<formula>formula see original document page 9</formula>其中Dfo和Mfo是系数Cd时，全息测量面为A到待测声源表面& 的传递矩阵。利用待测声源表面&的互谱相关密度Helmholtz公式作为限制方程，也可9以得到一对离散的矩阵方程<formula>formula see original document page 10</formula>从方程(4本文档来自技高网...

【技术保护点】
一种采用边界元法重建平稳声源的方法，其特征在于，包括以下步骤：第一步、设置待测声源表面为Ｓ↓［ｓ］，全息测量面为Ｓ↓［ｈ］，全息测量面Ｓ↓［ｈ］是待测声源表面Ｓ↓［ｓ］附近的一个任意面；第二步、布置传声器阵列，在声源附近布置 一个相对于声源位置固定不变的参考源传声器用以采集参考信号；第三步、同步记录参考源传声器和传声器列阵采集到的参考信号和全息测点声压，并用空间定位仪确定传声器列阵中每个传声器的空间位置，一并存储在存储设备中；第四步、从参考信号中 获取参考相位，利用参考信号与全息测点声压之间的相位关系，得到传声器阵列所采集到的全息测点声压的相对相位关系；第五步、由边界元法重建计算获得待测声源信号的自谱相关密度向量。

【技术特征摘要】

【专利技术属性】
技术研发人员：蒋伟康，张海滨，万泉，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：31[中国|上海]