CN102944867A - 一种无需相位信息的点声源识别方法 - Google Patents
一种无需相位信息的点声源识别方法 Download PDFInfo
- Publication number
- CN102944867A CN102944867A CN2012103806684A CN201210380668A CN102944867A CN 102944867 A CN102944867 A CN 102944867A CN 2012103806684 A CN2012103806684 A CN 2012103806684A CN 201210380668 A CN201210380668 A CN 201210380668A CN 102944867 A CN102944867 A CN 102944867A
- Authority
- CN
- China
- Prior art keywords
- sound source
- vector
- point
- face
- sound
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明涉及基于传声器阵列的声源识别领域,是一种基于柯西施瓦茨不等式数学定理、无需测量声压相位信号的声源识别方法。
背景技术
声源识别是噪声控制工程的重要内容之一,在产品的低噪声设计和声学故障诊断等方面有着重要的意义。通过声源识别获得声源的位置、强度等信息,从声源环节降低噪声是最有效的噪声控制方式;通过分析声信号来获得机械设备运转的信息,能够及时发现故障避免损失。
传统的噪声测量方法如覆盖法、分别运转法、声强法,这些传统的声源测量方法都存在着一定的局限性,获得的只是测量点处的声压或声强信息。随着计算机技术和信号处理技术的发展,现已采用传声器阵列对声源的信号进行测量,然后采用一定的算法进行声源识别,以获得声源的信息。
使用一定数量的传声器构成一个测量平面,放置于测量面上的传声器捕捉声源辐射到该点的声压信号,测量获得声压信号的幅值和相位信息。目前使用的有声全息方法和波束形成方法。在这些方法中相位信息是必不可少的,能否准确地获得传声器的相位信息对声源识别的结果具有很大的影响。使用到相位信息以后,必然存在各传声器之间相位匹配的问题,实现各传声器直接的相位匹配需要额外的投入资金和精力,并且在实际测量中相位比幅值更容易受到噪声信号的干扰而发生偏差,影响测量精度。目前还没有一种基于传声器阵列测量的声源识别方法不需要测量相位信息。
发明内容
本发明的目的在于克服上述现有技术存在的问题而提出一种无需相位信息的点声源识别方法,只需要测量声压的幅值信息,不存在各传声器之间相位匹配的问题,能实现点声源的准确识别。
本发明为实现上述目的采用的技术方案是按如下步骤:(1) 在平行于声源面的测量面上,用 规则矩形网格传声器阵列测量声源面上处点声源辐射到传声器阵列面上的声压,获取声压幅值向量为,;(2)构造声源面上任意位置的某一个点声源在测量面形成的声压幅值向量为; (3)对向量归一化处理,得到归一化向量;(4)将向量与归一化向量做内积运算,有:,,;当和线性相关时,和线性相关;(5)搜索最大值,当获得最大值时,归一化向量中对应的点声源的位置就是点声源s在声源面的位置; (6)计算得到声源的强度。
本发明利用柯西施瓦茨不等式原理,基于内积进行相关识别,识别过程中不需要测量传声器声压信号相位,只需要测量声压幅值便可进行点源的准确识别;仅将点声源信号的幅值信息与关于声源位置的向量做内积运算,根据柯西施瓦茨定理,应用单纯搜索法优化搜索内积模的最大值,此时基向量中声源位置就是实际声源所处的位置;识别出声源的位置后,通过内积运算计算出声源强度。方法简单,实际工程测量中更容易实现,能够实现点声源的准确识别。声源识别精度高。测量传声器幅值是传声器最基本的功能,不需要测量相位就不存在相位匹配的问题,因此对仪器设备的要求降低,经济实用。
附图说明
图1是声源识别测量示意图;
图2是本发明的流程图;
具体实施方式
设空间某处有一点声源,该点声源向空间中辐射的声压为:
则点声源在空间声场某一位置辐射的声压可以表示为:
(3)
如附图1所示的面上有一点声源,点声源位置用坐标表示为,点声源所在的平面称之为声源面,在平行于声源面的测量面上用一个传声器阵列测量点声源辐射到传声器阵列面上的声压,测量面距离点声源的距离为d。测量面某一位置接收到点声源辐射的声压信号,经过傅里叶变换可以表示为,其中, ,是测量面上某一点到点声源的距离。测量面上某一点的声压幅值函数是:
为某一常数,式(7)中等式左右两端可进一步表示为:
(8)
上述过程是通过函数理论分析获得,工程实现过程需要经过离散化,如图2所示,实施步骤如下:
1、用传声器阵列获取声源的声压幅值向量
整个传声器阵列测得的声压向量为:
(10)
2、构造声源面上某一点源在测量面形成的声压幅值向量
在传声器阵列上产生的声压幅值向量为:
3、的归一化处理
4、声压幅值向量和归一化的向量做内积运算
5、优化计算的最大值获取声源位置
6、计算声源的强度
根据式(16),计算出点声源的强度:
通过上述的理论计算可以识别声源的位置和强度。
以下提供本发明的一个实施例:
实施例
设声源面上存在一个点声源,点源的频率为1200Hz,声源强度为,点源的位置位于(0.12m,-0.18m)处.。在距离声源面处布置传声器阵列,对该点声源辐射的声压信号进行数据采集,传声器之间的间距为0.1m,传声器阵列的数目为。对传声器测量的时域声压信号进行傅里叶变换,经频谱分析后获得声压幅值(声压的相位信息不需要),根据公式(10),其中,获取测量阵列的声压幅值向量。由公式(12),其中,构造声源面上某一点源辐射到测量面上的声压幅值向量。对进行归一化处理,获得归一化后的向量。将传声器阵列获得声压幅值向量和归一化后的向量进行内积运算,并计算内积的模。使用单纯搜索法求取内积模取最大值,内积模最大值时向量中对应声源位置,根据柯西施瓦茨不等式,此时就是实际声源所在的位置。根据公式(17)识别出声源的强度值。识别结果见表1:
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210380668.4A CN102944867B (zh) | 2012-10-10 | 2012-10-10 | 一种无需相位信息的点声源识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210380668.4A CN102944867B (zh) | 2012-10-10 | 2012-10-10 | 一种无需相位信息的点声源识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102944867A true CN102944867A (zh) | 2013-02-27 |
CN102944867B CN102944867B (zh) | 2014-06-25 |
Family
ID=47727828
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210380668.4A Expired - Fee Related CN102944867B (zh) | 2012-10-10 | 2012-10-10 | 一种无需相位信息的点声源识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102944867B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765038A (zh) * | 2015-03-27 | 2015-07-08 | 江苏大学 | 一种基于内积相关性原理追踪运动点声源轨迹的方法 |
CN107064877A (zh) * | 2017-05-11 | 2017-08-18 | 江苏大学 | 一种运用向量正交化识别多点声源的方法 |
CN109100727A (zh) * | 2018-06-06 | 2018-12-28 | 江苏大学 | 一种使用运动虚拟声源识别声源运动轨迹的方法 |
CN109696659A (zh) * | 2019-01-21 | 2019-04-30 | 柳州市展虹科技有限公司 | 一种自由声场批量传声筒幅值灵敏度与相位测量装置智能控制系统 |
CN109696660A (zh) * | 2019-01-21 | 2019-04-30 | 柳州市展虹科技有限公司 | 一种用于自由声场传声筒幅值灵敏度与相位检测的精确测量方法 |
CN109975812A (zh) * | 2017-12-27 | 2019-07-05 | 中国科学院声学研究所 | 一种爆炸声的声源深度和声源级的计算方法 |
CN110459220A (zh) * | 2019-08-26 | 2019-11-15 | 杭州涂鸦信息技术有限公司 | 一种语音唤醒方法和系统以及可读存储介质、计算机 |
CN111537058A (zh) * | 2020-04-16 | 2020-08-14 | 哈尔滨工程大学 | 一种基于Helmholtz方程最小二乘法的声场分离方法 |
CN113281727A (zh) * | 2021-06-02 | 2021-08-20 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006044051A2 (en) * | 2004-10-13 | 2006-04-27 | Wayne State University | Farfield analysis of noise sources |
CN102508204A (zh) * | 2011-11-24 | 2012-06-20 | 上海交通大学 | 基于波束形成和传递路径分析的室内噪声源定位方法 |
CN102680071A (zh) * | 2012-05-21 | 2012-09-19 | 哈尔滨工程大学 | 采用振速测量和局部近场声全息法的噪声源识别方法 |
-
2012
- 2012-10-10 CN CN201210380668.4A patent/CN102944867B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006044051A2 (en) * | 2004-10-13 | 2006-04-27 | Wayne State University | Farfield analysis of noise sources |
CN102508204A (zh) * | 2011-11-24 | 2012-06-20 | 上海交通大学 | 基于波束形成和传递路径分析的室内噪声源定位方法 |
CN102680071A (zh) * | 2012-05-21 | 2012-09-19 | 哈尔滨工程大学 | 采用振速测量和局部近场声全息法的噪声源识别方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765038A (zh) * | 2015-03-27 | 2015-07-08 | 江苏大学 | 一种基于内积相关性原理追踪运动点声源轨迹的方法 |
CN107064877A (zh) * | 2017-05-11 | 2017-08-18 | 江苏大学 | 一种运用向量正交化识别多点声源的方法 |
CN107064877B (zh) * | 2017-05-11 | 2019-12-03 | 江苏大学 | 一种运用向量正交化识别多点声源的方法 |
CN109975812A (zh) * | 2017-12-27 | 2019-07-05 | 中国科学院声学研究所 | 一种爆炸声的声源深度和声源级的计算方法 |
CN109100727A (zh) * | 2018-06-06 | 2018-12-28 | 江苏大学 | 一种使用运动虚拟声源识别声源运动轨迹的方法 |
CN109696659A (zh) * | 2019-01-21 | 2019-04-30 | 柳州市展虹科技有限公司 | 一种自由声场批量传声筒幅值灵敏度与相位测量装置智能控制系统 |
CN109696660A (zh) * | 2019-01-21 | 2019-04-30 | 柳州市展虹科技有限公司 | 一种用于自由声场传声筒幅值灵敏度与相位检测的精确测量方法 |
CN109696660B (zh) * | 2019-01-21 | 2024-03-22 | 广西科技大学 | 一种检测自由声场传声筒幅值灵敏度与相位精确测量方法 |
CN110459220A (zh) * | 2019-08-26 | 2019-11-15 | 杭州涂鸦信息技术有限公司 | 一种语音唤醒方法和系统以及可读存储介质、计算机 |
CN111537058A (zh) * | 2020-04-16 | 2020-08-14 | 哈尔滨工程大学 | 一种基于Helmholtz方程最小二乘法的声场分离方法 |
CN113281727A (zh) * | 2021-06-02 | 2021-08-20 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
CN113281727B (zh) * | 2021-06-02 | 2021-12-07 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102944867B (zh) | 2014-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102944867B (zh) | 一种无需相位信息的点声源识别方法 | |
Kuang et al. | A complete characterization and solution to the microphone position self-calibration problem | |
CN101403634B (zh) | 压力速度法声场分离方法 | |
McBrearty et al. | Earthquake arrival association with backprojection and graph theory | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN109341845B (zh) | 一种海洋环境稳态声场空间实时仿真的方法及装置 | |
CN113064147B (zh) | 一种低信噪比下新型匹配场被动定位方法 | |
CN106596088A (zh) | 基于近场声源聚焦定位的碰摩声发射故障位置识别方法 | |
Ammari et al. | Shape identification and classification in echolocation | |
JP2008249702A (ja) | 音響測定装置及び音響測定方法 | |
Kita et al. | Fundamental study on sound source localization inside a structure using a deep neural network and computer-aided engineering | |
Sun et al. | Array geometry calibration for underwater compact arrays | |
CN113159447B (zh) | 一种激光雷达电磁环境效应预测方法及系统 | |
Nouri et al. | Providing a simple and efficient algorithm for passive sound localization | |
CN112503399B (zh) | 一种地下管线位置的确定方法、装置及其存储介质 | |
CN107797091B (zh) | 一种基于子空间的新型纯方位目标定位方法 | |
CN107064877B (zh) | 一种运用向量正交化识别多点声源的方法 | |
Flückiger et al. | Optimization of receiver arrangements for passive emitter localization methods | |
CN106124044A (zh) | 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 | |
Pang et al. | Impact localization method for composite plate based on low sampling rate embedded fiber bragg grating sensors | |
Zhayida et al. | Minimal solutions for dual microphone rig self-calibration | |
US20150293244A1 (en) | Estimation of a far field signature in a second direction from a far field signature in a first direction | |
Larsson et al. | Fast and robust stratified self-calibration using time-difference-of-arrival measurements | |
Yang et al. | Passive localization based on multipath time-delay difference with two hydrophones in deep ocean | |
CN108594163B (zh) | 用单个传声器移动结合内积运算识别声源的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140625 Termination date: 20171010 |