CN111998934B - 一种声源声功率测试方法 - Google Patents

一种声源声功率测试方法 Download PDF

Info

Publication number
CN111998934B
CN111998934B CN202010881820.1A CN202010881820A CN111998934B CN 111998934 B CN111998934 B CN 111998934B CN 202010881820 A CN202010881820 A CN 202010881820A CN 111998934 B CN111998934 B CN 111998934B
Authority
CN
China
Prior art keywords
sound
sound source
power
reciprocal
loudspeaker
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.)
Active
Application number
CN202010881820.1A
Other languages
English (en)
Other versions
CN111998934A (zh
Inventor
卢铃
车垚
彭继文
谈翀
徐莲环
胡伟
蔡炜
杨舟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010881820.1A priority Critical patent/CN111998934B/zh
Publication of CN111998934A publication Critical patent/CN111998934A/zh
Application granted granted Critical
Publication of CN111998934B publication Critical patent/CN111998934B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
    • G01H11/06Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H7/00Measuring reverberation time ; room acoustic measurements

Abstract

本发明公开了一种声源声功率测试方法,该方法包括下述步骤:(1)、在一个混响空间中,让声源设备正常运行,测量互易扬声器的开路电压;(2)驱动互易扬声器发声,测量互易扬声器的输入电流,同时将声源设备移除,在声源位置处放置传声器,获得传声器声压;(3)根据公式计算得到声源设备不同频率的声功率;(4)改变互易扬声器的位置,重复步骤1、2、3,并计算不同频率的声功率的多次测量平均值;(5)将不同频率下的声功率的平均值进行求和,获得声源设备总的声功率。该方法能基于电磁互易原理测量得到声源宽频带的声特性,且理论依据充分,系统搭建简单,特别适用于扬声器、变压器、电容器、风机等多种设备的声功率特性测试。

Description

一种声源声功率测试方法
技术领域
本发明属于声学测量领域,具体涉及一种声源声功率测试方法。
背景技术
基于互易原理可以在复杂声场环境中测量得到声源的低频特性(体积速度),而不需要消声室、消声水池等造价昂贵的实验设施,因此具有极好的应用潜力。但是在测试中,传声器的指向性必须与声源的指向性必须完全一致。对于全指向性的点声源或者某些低频声源,传声器可以满足该条件;而对于复杂的声源,其指向性可能十分复杂,找到指向性与之相匹配的传声器几乎不可能。这意味着,基于互易原理测量声源特性仅适用于全指向性的点声源,因此适用范围有限。
对此,基于互易原理的测量声源的高频特性的简易计算方法尚未出现。
发明内容
本发明要解决的技术问题是,针对互易原理测量声源体积速度,只能在低频有效的问题,提出了一种基于互易原理宽频带的声源声功率测试方法。该方法可以获得声源宽频带的声功率,扩展了基于互易原理的声源特性测量应用范围,使得该方法可应用到扬声器、变压器、电容器、风机等各种设备的噪声源特性测试。
为解决上述问题,本发明提供了一种声源声功率测试方法,在混响室中的互易扬声器和传声器满足电声互易条件,所述测试方法具体包含如下步骤:
步骤1:在所述混响室中,使声源设备正常运行发出噪声,测量互易扬声器的开路电压E;
步骤2:驱动互易扬声器发声,测量互易扬声器的输入电流I,同时将声源设备移除,在声源位置处放置传声器,获得传声器声压P0
步骤3:将开路电压E,输入电流I,声压P0变换到频域,根据公式计算声功率随频率的分布:
Figure BDA0002654303640000021
其中,ρ表示介质密度,f表示频率,
Figure BDA0002654303640000022
为频率f下的开路电压E,
Figure BDA0002654303640000023
为频率f下的输入电流,c为声速;
步骤4:改变所述互易扬声器的位置,重复步骤1~3,计算多次试验的得到Pf的平均值
Figure BDA0002654303640000024
步骤5:根据步骤4得到的平均值
Figure BDA0002654303640000025
计算出全声源发声范围频率下的声源总的声功率P。
进一步的,所述步骤2中,驱动互易扬声器发声时,发声频率范围覆盖全部声源发声范围频率。
进一步的,所述步骤4中的多次试验为5次以上的实验。
进一步的,所述步骤5中,对所有频率f的声功率平均值
Figure BDA0002654303640000031
进行离散积分,得到声源总的声功率:
Figure BDA0002654303640000032
进一步的,在计算得到声功率P后,所述声功率P用于小型声源的噪音测试。
进一步的,在计算得到声功率P后,所述声功率P用于水下工作设备、扬声器、变压器、电容器或风机的噪声源特性测试。
此外,本发明还公开了一种声源声功率测试系统,包括:
至少一个处理器以及与所述处理器通信连接的至少一个存储器,其中:所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如上述任一项的声源声功率测量方法。
此外,本发明还公开了一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如上述任一项所述的声源声功率测试方法。
本发明提出了一种声源声功率测试方法,与现有技术相比,本发明具有如下优点和有益效果:
1)本发明为声源声功率测量提供了一种新的方法,利用一个强混响的空间,使用一个满足电声互易条件的互易扬声器、一个传声器,构建一套声源声功率测试系统,在该测试系统使用声源声功率测量方法后,可根据发明人推导出来的简易公式计算得到声源声功率。本发明的技术核心是,利用强混响空间、互易扬声器和传声器构成一个互易声场,并根据声场及电声互易理论分析得到声源声功率与各个物理量之间的关系,从而通过测量所需物理量快速准确的计算出的声源宽频带声功率。
2)此外,该方法能测量得到声源宽频带的声特性(声功率),且该方法理论依据充分扎实,满足电声互易条件的测试系统搭建简单,计算方法简单且快速,能够依赖软件实现,特别适用于扬声器、变压器、电容器、风机等各种设备的噪声声功率特性测试。
附图说明
图1为本发明中声源声功率测试系统的模型示意图。
图2为本发明的利用电声互易原理测量声源辐射声功率的系统示意图。
图3为本发明的声源声功率测试方法的方法流程图。
具体实施方式
下面将结合附图和实施例对本发明进行清楚、完整地描述,同时也叙述了本发明技术方案解决的技术问题及有益效果,需要指出的是,所描述的实施例仅旨在便于对本发明的理解,而对其不起任何限定作用。
本发明基于背景技术中提及的不易发现的技术问题和基于图1的理论推导,专门提供了一种声源声功率测试方法,该测试方法依赖于图2的(a)-(b)中所示的在强混响的空间中满足电声互易条件的互易扬声器和一个传声器构成的声源声功率测试系统,在该系统环境下进行图3所示的声源声功率测试方法后,根据声场及互易理论分析得到声源声功率与各个物理量之间的关系,从而通过测量所需物理量计算出声源宽频带的声功率,从而使得该声功率能够应用到工业噪声的监控场合。
首先,如图1所示,本发明声源声功率测试方法的的理论依据如下:
图1的(a)中,o点为坐标中心,弧线表示半径为a的波振面。假设在o点有一个声源,声源的均方体积速度为
Figure BDA0002654303640000051
那么在波振面上总的均方体积速度也会是
Figure BDA0002654303640000052
如果在波振面上,取面积为dS的面元(dS的维度远远小于波长),那么面元的均方体积速度为U22dS/S。同样,若只在波振面上有总的均方体积速度是
Figure BDA0002654303640000054
的面声源,那么在o点处的声压为:
Figure BDA0002654303640000055
其中,ρ表示介质密度,f表示频率,P0表示自由场声压。若面元均匀分布在波振面上,得到所有分布在波振面上的面元在o点处的总声压为:
Figure BDA0002654303640000056
该公式的意义如图1的(d)所示。由于
Figure BDA0002654303640000057
是均匀分布,尺寸远小于波振面半径a的o点的附近的区域可以认为是近似的扩散声场。
将一个不规则外形的声接收器放置在o点,如图1的(b)所示。在a足够大的情况下,可以认为声接收器被放置在一个扩散声场中。在面元
Figure BDA0002654303640000061
的作用下,声接收器的振膜处的受挡声压为
Figure BDA0002654303640000062
当声接收器变成声源时,如图1的(c)所示,如果其声源均方体积速度为
Figure BDA0002654303640000063
在面元dS处的声压为
Figure BDA0002654303640000064
此时,通过面元的声功率可以表示为:
Figure BDA0002654303640000065
如果整个系统是互易的,那么有:
Figure BDA0002654303640000066
变形得到:
Figure BDA0002654303640000067
将公式(5)代入公式(3),消去
Figure BDA0002654303640000068
得到通过面元dS的声功率为:
Figure BDA0002654303640000069
对公式(6)积分得到通过波振面的总声功率为:
Figure BDA00026543036400000610
公式(7)即为声源声功率表达式,式中
Figure BDA0002654303640000071
物理意义是声接收器放置在一个均方声压为
Figure BDA0002654303640000072
的扩散声场中时,在声接收器的振膜处的受挡声压,如图1的(e)所示。
当声源具有指向性的时候,采用公式(7)进行实验,但是发明人发现在工作中实际上利用公式(7)进行实验几乎不可能,因为对于任意的声源,
Figure BDA0002654303640000073
Figure BDA0002654303640000074
很难直接测得,因为往往不规则的声源均方体积速度
Figure BDA0002654303640000075
的测量依往往需要激光测振仪测得,且
Figure BDA0002654303640000076
测得要求传声器的指向性必须与声源的指向性必须完全一致,这样的条件很难满足。但是发明人发现,基于互易原理,可以巧妙地解决这一问题,故本发明基于电声互易原理,将对未知声源的
Figure BDA0002654303640000077
Figure BDA0002654303640000078
的测量,转换为对辅助扬声器正向实验中开路电压E2和反向实验中输入电流I2的测量(由此可知本发明的技术问题发现和技术手段提出并不常规)。
在混响室中,利用一个辅助扬声器发声(输入电流的均方为I2)形成一个扩散声场。假设该声场的效果正好等效于图1的(e)中总均方体积速度为
Figure BDA0002654303640000079
的波振面,那么在声接收器的振膜处的受挡声压同样为
Figure BDA00026543036400000710
当声接收器发声时,假设其均方体积速度仍为
Figure BDA00026543036400000711
假设此时辅助扬声器的开路电压为E2。如果声接收器与辅助扬声器构成的系统是互易的,那么根据电声互易原理有:
Figure BDA00026543036400000712
即:
Figure BDA0002654303640000081
将公式(9)代入公式(7),可得到:
Figure BDA0002654303640000082
公式(10)中,E2、I2
Figure BDA0002654303640000083
均可以测量得到,c为声速。考虑到多次测量可提高准确性,因此可改变互易扬声器的位置,重复上述操作,计算平均值:
Figure BDA0002654303640000084
进一步对频率进行积分,得到总的声功率:
Figure BDA0002654303640000085
由此,可以实现对声源的辐射声功率测试,且经过发明人多次实验后发现该方案完全可行。
下面结合理论依据和具体实施例对本发明的测试方法作进一步说明。
如图3所示,本发明提出了一种声源声功率测试方法,在如图2所示的混响室中,该方法中使用的互易扬声器和传声器需满足电声互易条件,所述声源声功率测试方法包含如下步骤:
步骤1:在所述混响室中,使声源设备正常运行发出噪声,测量互易扬声器的开路电压E;
步骤2:驱动互易扬声器发声(发声频率范围覆盖全部声源发声范围频率),测量互易扬声器的输入电流I,同时将声源设备移除,在声源位置处放置传声器,获得传声器声压P0
步骤3:将开路电压E,输入电流I,声压P0变换到频域,根据公式计算声功率随频率的分布:
Figure BDA0002654303640000091
其中,ρ表示介质密度,f表示频率,
Figure BDA0002654303640000092
为频率f下的开路电压E,
Figure BDA0002654303640000093
为频率f下的输入电流,c为声速;
步骤4:改变互易扬声器的位置,重复步骤1~3,计算多次(如5次以上)试验的得到Pf的平均值
Figure BDA0002654303640000094
步骤5:根据步骤4得到的平均值
Figure BDA0002654303640000095
计算出全声源发声范围频率下的声源总的声功率P。
进一步的,可以对所有频率f的声功率平均值
Figure BDA0002654303640000096
进行离散积分,得到声源总的声功率:
Figure BDA0002654303640000097
在根据上述测试方法得到声功率P后,可将该声功率P用于工业环境中的扬声器、变压器、电容器、风机等多种设备的噪声源特性测试,也可适用于水下工作设备的水声噪音测试(如水下机器人等),或者更广泛的,该方法也可适用于波长小或者空间小的小型声源的噪音测试,从而对高频率/超低频率的声源噪声进行测试。
需要说明的是,上述的声源声功率测试方法可以作为软件程序或者计算机指令在非暂态计算机可读存储介质中执行或者在带有存储器和处理器的控制系统中执行,且其计算程序简单快速。在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。上述以软件功能单元的形式实现的集成的单元,可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)或处理器(processor)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
最后说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种声源声功率测试方法,其特征在于,在混响室中的互易扬声器和传声器满足电声互易条件,所述测试方法具体包含如下步骤:
步骤1:在所述混响室中,使声源设备正常运行发出噪声,测量互易扬声器的开路电压E;
步骤2:驱动互易扬声器发声,测量互易扬声器的输入电流I,同时将声源设备移除,在声源位置处放置传声器,获得传声器声压P0
步骤3:将开路电压E,输入电流I,声压P0变换到频域,根据公式计算声功率随频率的分布:
Figure FDA0003520997630000011
其中,ρ表示介质密度,f表示频率,Ef为频率f下的互易扬声器的开路电压,If为频率f下的互易扬声器的输入电流,c为声速;
步骤4:改变所述互易扬声器的位置,重复步骤1~3,计算多次试验的得到Pf的平均值
Figure FDA0003520997630000012
步骤5:根据步骤4得到的平均值
Figure FDA0003520997630000013
计算出全声源发声范围频率下的声源总的声功率P。
2.根据权利要求1所述的声源声功率测试方法,其特征在于,所述步骤2中,驱动互易扬声器发声时,发声频率范围覆盖全部声源发声范围频率。
3.根据权利要求1所述的声源声功率测试方法,其特征在于,所述步骤4中的多次试验为5次以上的实验。
4.根据权利要求1所述的声源声功率测试方法,其特征在于,所述步骤5中,对所有频率f的声功率平均值
Figure FDA0003520997630000021
进行离散积分,得到声源总的声功率:
Figure FDA0003520997630000022
5.根据权利要求1所述的声源声功率测试方法,其特征在于,在计算得到声功率P后,所述声功率P用于小型声源的噪音测试。
6.根据权利要求5所述的声源声功率测试方法,其特征在于,在计算得到声功率P后,所述声功率P用于水下工作设备、扬声器、变压器、电容器或风机的噪声源特性测试。
7.一种声源声功率测试系统,其特征在于,包括:
至少一个处理器以及与所述处理器通信连接的至少一个存储器,其中:所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至6任一项所述的声源声功率测试方法。
8.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至6任一项所述的声源声功率测试方法。
CN202010881820.1A 2020-08-28 2020-08-28 一种声源声功率测试方法 Active CN111998934B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010881820.1A CN111998934B (zh) 2020-08-28 2020-08-28 一种声源声功率测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010881820.1A CN111998934B (zh) 2020-08-28 2020-08-28 一种声源声功率测试方法

Publications (2)

Publication Number Publication Date
CN111998934A CN111998934A (zh) 2020-11-27
CN111998934B true CN111998934B (zh) 2022-07-01

Family

ID=73471349

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010881820.1A Active CN111998934B (zh) 2020-08-28 2020-08-28 一种声源声功率测试方法

Country Status (1)

Country Link
CN (1) CN111998934B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115683323B (zh) * 2022-12-29 2023-05-02 杭州兆华电子股份有限公司 一种发声设备的声功率测量系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106501795A (zh) * 2016-11-09 2017-03-15 哈尔滨工程大学 一种利用混响水池进行水声换能器互易校准的方法
CN107436190A (zh) * 2016-05-27 2017-12-05 中国科学院声学研究所 一种电声互易装置的非互易补偿方法
WO2018164438A1 (en) * 2017-03-10 2018-09-13 Samsung Electronics Co., Ltd. Method and apparatus for in-room low-frequency sound power optimization
CN108627240A (zh) * 2018-03-23 2018-10-09 南京大学 一种恒体积速度声源低频声功率的估算方法
CN111586546A (zh) * 2020-03-26 2020-08-25 中国船舶重工集团公司第七一五研究所 一种低频换能器谐振点发射响应的测量方法和系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107436190A (zh) * 2016-05-27 2017-12-05 中国科学院声学研究所 一种电声互易装置的非互易补偿方法
CN106501795A (zh) * 2016-11-09 2017-03-15 哈尔滨工程大学 一种利用混响水池进行水声换能器互易校准的方法
WO2018164438A1 (en) * 2017-03-10 2018-09-13 Samsung Electronics Co., Ltd. Method and apparatus for in-room low-frequency sound power optimization
CN108627240A (zh) * 2018-03-23 2018-10-09 南京大学 一种恒体积速度声源低频声功率的估算方法
CN111586546A (zh) * 2020-03-26 2020-08-25 中国船舶重工集团公司第七一五研究所 一种低频换能器谐振点发射响应的测量方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种用于声源低频性能测试的互易装置;卢铃 等;《应用声学》;20161130;第35卷(第6期);第520-526页 *
混响场中水声换能器互易校准研究;孙俊东 等;《声学技术》;20160831;第35卷(第4期);第44-48页 *

Also Published As

Publication number Publication date
CN111998934A (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
EP3232686B1 (en) Neural network-based loudspeaker modeling with a deconvolution filter
KR20070072518A (ko) 소음원의 원거리-장 분석
CN108683982B (zh) 一种音箱性能测试方法、系统及音箱共振性测试方法
KR20140053283A (ko) 노이즈를 제어하는 전자 디바이스들
CN111998934B (zh) 一种声源声功率测试方法
Ju et al. Near-field characteristics of the parametric loudspeaker using ultrasonic transducers
Melon et al. Evaluation of a method for the measurement of subwoofers in usual rooms
Abel et al. Synchronization of sound sources
CN108229030B (zh) 一种主动降噪系统控制器参数的设计方法
Shi et al. The use of wideband acoustical holography for noise source visualization
WO2017220816A1 (en) Method for simulating total harmonic distortion of a loudspeaker
Andersen et al. 3D shape optimization of loudspeaker cabinets for uniform directivity
CN111159945A (zh) 一种基于主辐射模态的水下圆柱壳低频声辐射预报方法
Koyama et al. Sound field decomposition in reverberant environment using sparse and low-rank signal models
Sanalatii et al. Estimation of loudspeaker frequency response and directivity using the radiation-mode method
Blanco et al. Experimental verification of projection algorithms and optimization routines for acoustic field uniformity enhancement in MIMO direct field acoustic control
Bradley et al. Numerical prediction of sound scattering from surfaces with fractal geometry: A preliminary investigation
Sagers et al. An extended lumped-element model and parameter estimation technique to predict loudspeaker responses with possible surround-dip effects
Kemp et al. Pressure fields in the vicinity of brass musical instrument bells measured using a two dimensional grid array and comparison with multimodal models
Chojnacki et al. Coupled speakers directivity measurements for small acoustic omnidirectional source development
Skudrzyk Sound radiation of a system with a finite or an infinite number of resonances
RU2782354C2 (ru) Устройство для градуировки электроакустических преобразователей
Murata et al. Sparse sound field decomposition with parametric dictionary learning for super-resolution recording and reproduction
CN117496997B (zh) 基于惩罚机制的声源检测方法、装置及存储介质
Cheng et al. Evaluation of free-field transfer functions in anomalous reverberant fields

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