CN110530510A - 一种利用线性声阵列波束形成的声源辐射声功率测量方法 - Google Patents
一种利用线性声阵列波束形成的声源辐射声功率测量方法 Download PDFInfo
- Publication number
- CN110530510A CN110530510A CN201910907324.6A CN201910907324A CN110530510A CN 110530510 A CN110530510 A CN 110530510A CN 201910907324 A CN201910907324 A CN 201910907324A CN 110530510 A CN110530510 A CN 110530510A
- Authority
- CN
- China
- Prior art keywords
- sound source
- sound
- wave beam
- beam forming
- array
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提出一种利用线性声阵列波束形成的声源辐射声功率测量方法,首先选择间距为声源辐射频点对应半波长的传声器线阵置于正对着待测声源区域的测量位置处,采集待测声源的辐射声压信号;对声压信号进行波束形成,得到波束形成空间谱;将波束形成空间谱与转化因子α相乘得到转化后的空间谱;再将转化后得到的空间谱在主瓣区域进行积分计算,从而得到声源辐射声功率。本发明不受稳态工况的限制,不需要布置大量的测点,也不需要进行测量网格面的划分,操作简便、测试时间短,效率高,重复性好,测量结果准确可靠,测试成本低,并且能够对测量区域存在多声源的条件下得到其中的每个声源的辐射声功率,能够快速适用于不同的工程问题进行现场测量。
Description
技术领域
本发明涉及声源辐射声功率测量技术领域,具体为一种利用线性声阵列波束形成的声源辐射声功率测量方法。
背景技术
声功率作为对噪声源辐射强度与特性的定量描述,是评价降噪效果的重要参量,也是声源识别的主要内容。因此,声功率的准确测量具有十分重要的意义。传统的声功率测量方法如声压法(GB/T 6882-2016等),声强法(GB/T 16404-1996等),虽然能够获得目标声源的辐射声功率,但是存在诸多弊端。
声强法是通过声强探头直接进行噪声源近场数据的测量,需要进行测量面的网格划分。虽然能够获得声源的声强数据从而计算出相应的辐射声功率,但是往往需要大量的时间进行测试,并且费用高昂,仅适合在稳态工况下进行等等,这些不足使得声强法在实际工程应用中会受到很大的制约。
声压法是指通过测量声压值获取声功率的测量方法,可以分为自由场和混响场法两类。混响场法则是把噪声源放在混响室内,通过测得室内平均声压级和混响时间,结合声功率与声压级以及混响室的吸声量的换算关系,来计算噪声源的辐射声功率。由于混响法对测量环境的要求限制,使得该方法难以应用于实际工程中的现场测量。自由场测试原理基本上是通过测量声场中所划分测点处的声压级,结合对应的声压与声功率的关系,将测试的声压级结果换算成声功率级。当声源具有指向性时,需要在声源远场处的假想球面上测出多点的声压级,需要在空间布置大量的测点,在实际工程应用的现场测量中,测点的选取也会受到测量条件的限制,操作复杂,测量繁琐,工作量大。
无论是声压法还是声强法在测量区域存在多声源的情况下,只能测量出该区域总声功率,并不能得到该区域其中每个声源单独的辐射声功率。
发明内容
针对现有技术存在的问题,本发明提出一种利用线性传声器阵列波束形成方法对测量区域的声源进行测试来获得声源的辐射声功率的方法,简单易行,重复性好,方便效率,大大减少了测试工作量,并且能够在多声源条件下获得其中每个声源的辐射声功率,能够更好地用于实际的工程应用。
本发明的技术方案为:
所述一种利用线性声阵列波束形成的声源辐射声功率测量方法,其特征在于:包括以下步骤:
步骤1:根据声源的辐射频谱特性,选择间距为声源辐射频点对应半波长的传声器线阵;然后将传声器线阵置于正对着待测声源区域的测量位置处,采集待测声源的辐射声压信号;
步骤2:将步骤1中传声器线阵所采集到的声压信号进行波束形成,得到波束形成空间谱;然后将波束形成空间谱与转化因子α相乘得到转化后的空间谱;再将转化后得到的空间谱在主瓣区域进行积分计算,从而得到声源辐射声功率;其中转化因子α根据公式
计算,L为传声器线阵与声源的法向距离,ρ0为空气密度,c0为空气声速,k为波数,λ为声源辐射频点对应波长,M为阵元个数,d为阵元间距,R1为主瓣区域的主瓣圆半径,根据公式计算得到。
进一步的优选方案,所述一种利用线性声阵列波束形成的声源辐射声功率测量方法,其特征在于:当声源为宽带声源时,对宽带声源频带范围内各频点分别按照步骤1和步骤2的过程进行计算,得到宽带声源各个频点处的辐射声功率的计算结果,将这些计算结果叠加得到宽带声源的声功率结果。
有益效果
本发明所提供的声功率测量方法不受稳态工况的限制,不需要布置大量的测点,也不需要进行测量网格面的划分,操作简便、测试时间短,效率高,重复性好,测量结果准确可靠,测试成本低,并且能够对测量区域存在多声源的条件下得到其中的每个声源的辐射声功率,能够快速适用于不同的工程问题进行现场测量。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1本发明的流程图;
图2测量系统示意图;
图3实施例中声源1的测量结果对比;
图4实施例声源2的测量结果对比。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
如图1所示,本发明中利用线性声阵列波束形成的声源辐射声功率测量方法,包括以下步骤:
步骤1:根据声源的辐射频谱特性,选择间距为声源辐射频点对应半波长的传声器线阵;然后将传声器线阵置于正对着待测声源区域的测量位置处,采集待测声源的辐射声压信号。
步骤2:将步骤1中传声器线阵所采集到的声压信号进行波束形成,得到波束输出结果即波束形成空间谱。
同时通过理论推导单极子点声源的辐射声功率与波束输出结果之间的换算关系建立起声源辐射声功率的计算模型:
积分的主瓣圆半径可以由式(1)计算得到,令阵列的波束输出结果经转化后在主瓣区域的积分结果与单极子理论计算的点声源辐射声功率值相等,得到波束输出结果与声源辐射声功率的转化因子α如式(2)。
其中L为传声器线阵与声源的法向距离,ρ0为空气密度,c0为空气声速,k为波数,λ为声源辐射频点对应波长,M为阵元个数,d为阵元间距,R1为主瓣区域的主瓣圆半径。
这样在实际测试中根据相应的阵列参数,测试距离,频谱等信息,确定积分主瓣圆半径以及计算相应的转化因子。
然后将波束形成空间谱与转化因子α相乘得到转化后的空间谱;再将转化后得到的空间谱在主瓣区域进行积分计算,从而得到声源辐射声功率。
当声源为宽带声源时,对宽带声源频带范围内各频点分别按照步骤1和步骤2的过程进行计算,得到宽带声源各个频点处的辐射声功率的计算结果,将这些计算结果叠加得到宽带声源的声功率结果。
下面以消声室中利用传声器阵列波束形成对双声源的辐射声功率测量为例,对本发明做详细说明。
(1)如图2所示,发射信号经由功率放大器B&K2716固定增益后加载到扬声器上发声,B&K3560配合PULSE软件进行数据的采集。声源1置于正对阵列中心位置,声源2偏离约0.6m,声源到阵列的测试距离约为1.8m。
(2)声源1发射单频信号,声源2发射宽带信号,通过各阵元接收的声信号数据,分别进行波束形成得到各自相应的空间谱,从而确定声源的位置,然后按照波束输出结果与声源辐射声功率的换算关系,分别计算出两声源在阵列侧的辐射声功率。
(3)利用声压法测量声功率,分别获得两声源单独发声时的辐射声功率,并作为对比参考值。对比结果如图3,图4所示。
结果分析:采用本发明提供的利用线性声阵列波束形成来计算声源辐射声功率方法,在消声室中对两位置处声源的声源辐射声功率计算结果与在声压法测量结果都能极好的吻合,偏差小于2dB,因此采用本发明提供的方法进行声源辐射声功率的测量是可信的。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (2)
1.一种利用线性声阵列波束形成的声源辐射声功率测量方法,其特征在于:包括以下步骤:
步骤1:根据声源的辐射频谱特性,选择间距为声源辐射频点对应半波长的传声器线阵;然后将传声器线阵置于正对着待测声源区域的测量位置处,采集待测声源的辐射声压信号;
步骤2:将步骤1中传声器线阵所采集到的声压信号进行波束形成,得到波束形成空间谱;然后将波束形成空间谱与转化因子α相乘得到转化后的空间谱;再将转化后得到的空间谱在主瓣区域进行积分计算,从而得到声源辐射声功率;其中转化因子α根据公式
计算,L为传声器线阵与声源的法向距离,ρ0为空气密度,c0为空气声速,k为波数,λ为声源辐射频点对应波长,M为阵元个数,d为阵元间距,R1为主瓣区域的主瓣圆半径,根据公式计算得到。
2.根据权利要求1所述一种利用线性声阵列波束形成的声源辐射声功率测量方法,其特征在于:当声源为宽带声源时,对宽带声源频带范围内各频点分别按照步骤1和步骤2的过程进行计算,得到宽带声源各个频点处的辐射声功率的计算结果,将这些计算结果叠加得到宽带声源的声功率结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910907324.6A CN110530510B (zh) | 2019-09-24 | 2019-09-24 | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910907324.6A CN110530510B (zh) | 2019-09-24 | 2019-09-24 | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110530510A true CN110530510A (zh) | 2019-12-03 |
CN110530510B CN110530510B (zh) | 2021-01-05 |
Family
ID=68669871
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910907324.6A Active CN110530510B (zh) | 2019-09-24 | 2019-09-24 | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110530510B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111405456A (zh) * | 2020-03-11 | 2020-07-10 | 费迪曼逊多媒体科技(上海)有限公司 | 一种网格化3d声场采样方法及系统 |
CN112161814A (zh) * | 2020-09-21 | 2021-01-01 | 西北工业大学 | 一种平面叶栅尾缘噪声识别测量方法与装置 |
CN112179656A (zh) * | 2020-09-21 | 2021-01-05 | 西北工业大学 | 一种移动式线性传声器阵列声源指向性测量方法与装置 |
CN113548134A (zh) * | 2021-06-28 | 2021-10-26 | 神龙汽车有限公司 | 一种空调通风噪音计算的系统和方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183150A (zh) * | 2007-11-30 | 2008-05-21 | 浙江大学 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
JP2008219884A (ja) * | 2007-03-02 | 2008-09-18 | Honda Motor Co Ltd | 移動音源からの音の抽出方法および装置 |
CN101398482A (zh) * | 2008-10-24 | 2009-04-01 | 西北工业大学 | 声接收基阵被动宽带检测中的噪声场数值计算方法 |
CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
CN105263085A (zh) * | 2011-01-13 | 2016-01-20 | 高通股份有限公司 | 用移动平台进行的可变波束成形 |
CN106023996A (zh) * | 2016-06-12 | 2016-10-12 | 杭州电子科技大学 | 基于十字形声阵列宽带波束形成的声识别方法 |
CN106096184A (zh) * | 2016-06-24 | 2016-11-09 | 中国电力科学研究院 | 一种基于声强法的变压器多点声源模型的噪声计算方法及系统 |
CN106813774A (zh) * | 2017-01-10 | 2017-06-09 | 中国计量大学 | 一种高强度聚焦超声的声强和声功率测量方法 |
CN106910493A (zh) * | 2017-02-22 | 2017-06-30 | 青岛理工大学 | 一种多环同轴圆环阵稳健超指向性波束形成方法 |
JP2017207399A (ja) * | 2016-05-19 | 2017-11-24 | 株式会社小野測器 | 音源探査装置及び音源探査方法 |
CN108627240A (zh) * | 2018-03-23 | 2018-10-09 | 南京大学 | 一种恒体积速度声源低频声功率的估算方法 |
CN109238436A (zh) * | 2018-09-01 | 2019-01-18 | 哈尔滨工程大学 | 混响水池中瞬态声源特性的测量方法 |
CN110244288A (zh) * | 2019-07-01 | 2019-09-17 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于聚焦归位原理的水声阵列信号处理方法 |
-
2019
- 2019-09-24 CN CN201910907324.6A patent/CN110530510B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008219884A (ja) * | 2007-03-02 | 2008-09-18 | Honda Motor Co Ltd | 移動音源からの音の抽出方法および装置 |
CN101183150A (zh) * | 2007-11-30 | 2008-05-21 | 浙江大学 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
CN101398482A (zh) * | 2008-10-24 | 2009-04-01 | 西北工业大学 | 声接收基阵被动宽带检测中的噪声场数值计算方法 |
CN105263085A (zh) * | 2011-01-13 | 2016-01-20 | 高通股份有限公司 | 用移动平台进行的可变波束成形 |
CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
JP2017207399A (ja) * | 2016-05-19 | 2017-11-24 | 株式会社小野測器 | 音源探査装置及び音源探査方法 |
CN106023996A (zh) * | 2016-06-12 | 2016-10-12 | 杭州电子科技大学 | 基于十字形声阵列宽带波束形成的声识别方法 |
CN106096184A (zh) * | 2016-06-24 | 2016-11-09 | 中国电力科学研究院 | 一种基于声强法的变压器多点声源模型的噪声计算方法及系统 |
CN106813774A (zh) * | 2017-01-10 | 2017-06-09 | 中国计量大学 | 一种高强度聚焦超声的声强和声功率测量方法 |
CN106910493A (zh) * | 2017-02-22 | 2017-06-30 | 青岛理工大学 | 一种多环同轴圆环阵稳健超指向性波束形成方法 |
CN108627240A (zh) * | 2018-03-23 | 2018-10-09 | 南京大学 | 一种恒体积速度声源低频声功率的估算方法 |
CN109238436A (zh) * | 2018-09-01 | 2019-01-18 | 哈尔滨工程大学 | 混响水池中瞬态声源特性的测量方法 |
CN110244288A (zh) * | 2019-07-01 | 2019-09-17 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于聚焦归位原理的水声阵列信号处理方法 |
Non-Patent Citations (2)
Title |
---|
WANG QIAN: "《Research on the Methods of Underwater VehicleNoise Source Identification》", 《2015 IEEE INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING, COMMUNICATIONS AND COMPUTING (ICSPCC)》 * |
杨洋: "《基于波束形成的发动机噪声源识别及声功率计算》", 《内燃机工程》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111405456A (zh) * | 2020-03-11 | 2020-07-10 | 费迪曼逊多媒体科技(上海)有限公司 | 一种网格化3d声场采样方法及系统 |
CN111405456B (zh) * | 2020-03-11 | 2021-08-13 | 费迪曼逊多媒体科技(上海)有限公司 | 一种网格化3d声场采样方法及系统 |
CN112161814A (zh) * | 2020-09-21 | 2021-01-01 | 西北工业大学 | 一种平面叶栅尾缘噪声识别测量方法与装置 |
CN112179656A (zh) * | 2020-09-21 | 2021-01-05 | 西北工业大学 | 一种移动式线性传声器阵列声源指向性测量方法与装置 |
CN113548134A (zh) * | 2021-06-28 | 2021-10-26 | 神龙汽车有限公司 | 一种空调通风噪音计算的系统和方法 |
CN113548134B (zh) * | 2021-06-28 | 2022-04-29 | 神龙汽车有限公司 | 一种空调通风噪音计算的系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110530510B (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110530510A (zh) | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 | |
CN106501795B (zh) | 一种利用混响水池进行水声换能器互易校准的方法 | |
CN106680376B (zh) | 一种基于三维声强阵列的隔声测量系统与测量方法 | |
CN104501939B (zh) | 一种利用单水听器测量非消声水池混响时间的反演方法 | |
WO2019214570A1 (zh) | 一种阵列天线总辐射功率的测量方法、装置和系统 | |
CN105823935B (zh) | 一种半电波暗室场均匀性测试装置及其方法 | |
CN109302667A (zh) | 一种水声发射换能器宽带发送响应的快速测量方法及装置 | |
CN110186546A (zh) | 基于粉红噪声的水听器灵敏度自由场宽带校准方法 | |
CN113868583B (zh) | 一种子阵波束聚焦的声源距离计算方法及系统 | |
CN109238436B (zh) | 混响水池中瞬态声源特性的测量方法 | |
CN109375198B (zh) | 一种低频声呐阵阻抗特性校准方法 | |
CN104501938A (zh) | 一种在置于空气中的矩形混响水槽内测量水下声源低频辐射声功率的方法 | |
CN110243521A (zh) | 一种薄板应力测量方法及薄板应力测量系统 | |
CN109324113B (zh) | 一种混响水池壁面平均吸声系数的测量方法 | |
CN206563731U (zh) | 一种基于三维声强阵列的隔声测量系统 | |
CN109324320A (zh) | 一种利用混响水池进行水听器批量校准的方法 | |
CN110082431A (zh) | 一种用于材料表面声阻抗测量的方法及装置 | |
CN109916497B (zh) | 一种在混响水槽测量水下声源甚低频辐射特性的方法 | |
CN101398482B (zh) | 声接收基阵被动宽带检测中的噪声场数值计算方法 | |
CN115656926A (zh) | 一种用于电力设备异响监测的定位识别方法和装置 | |
CN112327281B (zh) | 一种高声扩散混响水池中窄带声信号的声源级测量方法 | |
CN108344497A (zh) | 水下声基阵互辐射阻抗的光学测试方法 | |
RU2782354C2 (ru) | Устройство для градуировки электроакустических преобразователей | |
CN216132987U (zh) | 一种基于近场声全息的隔声测量系统 | |
CN110332986A (zh) | 一种湍流控制屏声学校准系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |