专利名称：采用偏相干技术辨识非平稳声源特性的方法
技术领域：
本发明涉及的是一种辨识声源特性的方法，特别是一种利用循环平稳声场的偏相干部分声场分离技术辨识非平稳声源特性的方法，属于物理类中的噪声领域。
背景技术：
对噪声源特性的了解是控制噪声的前提，因此，为了有效地控制噪声，在降噪措施实施以前，必须首先进行噪声源诊断，确定各主要噪声源地位置，及其特性。随着现代信号处理技术地发展，谱分析技术，相干和偏相干分析技术、声强分析技术以及声近场全息技术等得到了迅速地发展。经文献检索发现，M.A.Tomlinson在《Applied Acoustics》(57(1999)243-261)上撰文“Partial source discrimination innear field acoustic holography”(《应用声学》，近场声全息中的部分源辨识)，该文提出了近场声全息偏相干分析的部分声源分离技术，可以有效地诊断噪声源。但是，这些技术只能分析平稳声场，因此，有必要提出新的技术，用于非平稳声场的分析。然而，对于一般的非平稳声场，声信号的统计特性参量是时变的，因而也就无法再用时间平均来代替集合平均，使得数据采集十分困难，很难分析声场的特性。
循环平稳信号是一类特殊的非平稳信号，由于自身独特的周期平稳性，使得单次采集到的记录具有周期遍历特性，这种特性提供了非平稳信号分析可能性。循环平稳信号在工程应用中有着十分重要的现实意义，例如旋转机械设备由于对称或近似对称的物理结构和周期性的工作运动模式，其声场有明显的周期时变特征，声信号具有循环平稳特性。在进一步检索中，尚未发现与本发明主题相同或者相类似的文献报道。
发明内容
本发明的目的在于克服现有技术中的不足，提供一种采用偏相干技术辨识非平稳声源特性的方法，使其实现得到各个声源分别形成的部分声场，通过循环平稳声场的近场全息重建，还可以由平面上的测量数据得到三维声场分布，实现各声源传播途径可视化。
本发明是通过以下技术方案实现的，本发明在多个声源产生的复杂循环平稳声场的场合，采用循环平稳理论替代传统的傅立叶变换，提出适用于循环平稳声场的偏相干分析技术，并将之用于循环平稳声场的部分声场分离中，首先确定参考源的数量和位置，并布置参考源传声器阵列，提取参考源信号，再设计传声器阵列对全息面进行扫描测量，采集全息面数据，然后，利用循环平稳偏相干分析技术，分离各个声源的部分声场。
以下对本发明方法作进一步的限定，方法步骤如下1、确定参考源的数量和位置，并布置参考源传声器阵列，提取参考源信号。为此，需要确定参考源的数目K和位置，可以采用平稳声场的近场综合试验分析(或称偏奇异值分析)得到。
2、布置全息面测量传声器阵列，如果通道数足够多，可以在整个全息平面上布置传声器；如果通道数不够，可以将传声器布置成线阵列，在整个全息面上进行扫描测量；在全息面数据采集的同时，利用参考源阵列采集参考源声信号。测量时，同步采集全部传声器通道的声信号，记录存储在计算机、磁带记录仪或其他设备中，设全息面通道数为Q。
3、分析所采集的声信号，选取能够反映声场特性的频率和循环频率，可以分析各个参考源信号的谱相关密度函数来选择。
4、利用采集到的全息面时域声信号数据和参考源的时域声信号数据，计算所选择的频率f与循环频率a上的参考源信号自谱相关密度矩阵 ，以及参考源信号和全息面上传感器信号的互谱相关密度矩阵 与 。可以通过这三个谱相关密度矩阵，计算获得全息面上声信号的谱相关密度矩阵 ，计算方法如式(1)[Sppα(f)]Q×Q=[Sprα(f)]Q×K[Srrα(f)]K×K-1[Srpα(f)]K×Q---(1)]]>5、利用采集到的全息面数据和参考源数据，进行循环偏相干分析，分离得到各个声源的部分声场，方法如下设下标r表示第r个参考源信号，下标p表示测量得到的全息平面上测点p的信号，则在除去r-1个参考源的影响后，第r个参考源在全息面测点处产生的输出信号的谱相关密度为Spr(r-1)!α(f)=γrp(r-1)!α2(f)Spp(r-1)!α(f)---(2)]]>其中γrp(r-1)!α2(f)=Srp.(r-1)!α(f)Spr(r-1)!α(f)Srr.(r-1)!α(f)Spp(r-1)!α(f)]]>为循环偏相干函数。
本发明具有实质性特点和显著进步，本发明利用循环平稳声场的独特的周期性，在偏相干分析技术的基础上，提出了循环平稳偏相干分析技术，并在循环偏相干技术和声近场全息技术的基础上，发明了循环偏相干分析部分声场分离技术，可以对多个循环平稳声源形成的复杂循环平稳声场进行分析，得到各个声源分别形成的部分声场；通过循环平稳声场的近场全息重建，还可以从全息平面上的测得的声压信号中推算出单个声源形成的三维声场分布，实现各声源传播途径可视化。
具体实施例方式
结合本发明方法内容提供以下实施例1、假设已知声源数目和位置，采用两个扬声器发声形成两个循环平稳声源，布置两个传声器，提取参考源信号。
扬声器的激励源为Vsource1＝Acos(2πf1t)*noise(t)Vsource2＝B(1+Ccos(2πfbt))*cos(2πfat)其中，A＝10，B＝C＝1，f1＝600，fa＝600，fb＝200，noise为带通白噪声。
2、将14个传声器布置成线阵，在整个全息面上进行扫描测量，形成14×14的全息面阵列；在全息面数据采集的同时，利用两个参考源传声器采集参考源声信号。测量时，同步采集全部传声器通道的声信号，记录存储在计算机或其他设备中。
3、在实验室回放分析所采集的声信号，选取能够反映声场特性的频率f＝200Hz和循环频率alpha＝1200Hz。
4、利用采集到的全息面数据和参考源数据，分析所选择的频率f与循环频率a上的参考源信号自谱相关密度矩阵 ，参考源信号和全息面上传感器信号的互谱相关密度矩阵 及 5、分离得到各个声源的部分声场。
6、利用循环平稳声场的近场全息重建公式进行声场重建，以得到各个声源部分声场的三维信息。
在该实施例中发现当得到的两个参考源信号相干性比较小时，分离效果非常好，可以明显地分析得到两个声源各自的声场，当得到的两个参考源信号相干性很强时，分离效果仍然十分理想，说明采用循环偏相干分析分离部分声场时在很大程度上可以削弱参考源阵列之间的内在相干性。
权利要求
1.一种采用偏相干技术辨识非平稳声源特性的方法，其特征在于，在多个声源产生的复杂循环平稳声场的场合，采用循环平稳理论替代传统的傅立叶变换，提出适用于循环平稳声场的偏相干分析技术，并将之用于循环平稳声场的部分声场分离中，首先确定参考源的数量和位置，并布置参考源传声器阵列，提取参考源信号，再设计传声器阵列对全息面进行扫描测量，采集全息面数据，然后，利用循环平稳偏相干分析技术，分离各个声源的部分声场。
2.根据权利要求1所述的采用偏相干技术辨识非平稳声源特性的方法，其特征是，以下对本发明作出进一步限定，具体的方法步骤如下(1)采用平稳声场的近场综合试验分析或称偏奇异值分析确定参考源的数量和位置，并布置参考源传声器阵列，提取参考源信号；(2)布置全息面测量传声器阵列，如果通道数足够多，在整个全息平面上布置传声器，否则将传声器布置成线阵列，在整个全息面上进行扫描测量，在全息面数据采集的同时，利用参考源阵列采集参考源声信号，测量时，同步采集全部传声器通道的声信号，记录存储在计算机、磁带记录仪或其他设备中，设全息面通道数为Q；(3)分析所采集的声信号，通过分析各个参考源信号的谱相关密度函数来选取反映声场特性的频率和循环频率；(4)利用采集到的全息面时域声信号数据和参考源的时域声信号数据，计算所选择的频率f与循环频率a上的参考源信号自谱相关密度矩阵 以及参考源信号和全息面上传感器信号的互谱相关密度矩阵 与 ，通过这三个谱相关密度矩阵，计算获得全息面上声信号的谱相关密度矩阵 ，计算方法如下[Sppα(f)]Q×Q=[Sprα(f)]Q×K[Srrα(f)]K×K-1[Srpα(f)]K×Q;]]>(5)利用采集到的全息面数据和参考源数据，进行循环偏相干分析，分离得到各个声源的部分声场。
3.根据权利要求2所述的采用偏相干技术辨识非平稳声源特性的方法，其特征是，步骤(5)中，利用采集到的全息面数据和参考源数据，进行循环偏相干分析，分离得到各个声源的部分声场，方法如下设下标r表示第r个参考源信号，下标p表示测量得到的全息平面上测点p的信号，则在除去r-1个参考源的影响后，第r个参考源在全息面测点处产生的输出信号的谱相关密度为Spr.(r-1)!α(f)=γrp.(r-1)!α2(f)Spp.(r-1)!α(f)]]>其中γrp.(r-1)!α2(f)=Srp.(r-1)!α(f)Spr(r-1)!α(f)Srr.(r-1)!α(f)Spp.(r-1)!α(f)]]>为循环偏相干函数。
4.根据权利要求1或2所述的采用偏相干技术辨识非平稳声源特性的方法，其特征是，通过循环平稳声场的近场全息重建，从全息平面上的测得的声压信号中获得单个声源形成的三维声场分布，实现各声源传播途径可视化。
全文摘要
一种采用偏相干技术辨识非平稳声源特性的方法，属于物理中的噪声类领域。本发明在多个声源产生复杂循环平稳声场的场合，采用循环平稳理论替代傅立叶变换，提出适用于循环平稳声场的偏相干分析技术，并将之用于循环平稳声场的部分声场分离中，方法如下首先确定参考源的数量和位置，并布置参考源传声器阵列，提取参考源信号，再设计传声器阵列对全息面进行扫描测量，采集全息面数据，然后，利用循环平稳偏相干分析技术，分离各个声源部分声场。本发明可对多个循环平稳声源形成的复杂循环平稳声场进行分析，得到各个声源分别形成的部分声场，还可从全息平面上测得的声压信号中推算出单个声源形成的三维声场分布，实现各声源传播途径可视化。
文档编号G01H17/00GK1472515SQ0312940
公开日2004年2月4日 申请日期2003年6月19日 优先权日2003年6月19日
发明者蒋伟康, 万泉 申请人:上海交通大学