CN107436190B - 一种电声互易装置的非互易补偿方法 - Google Patents
一种电声互易装置的非互易补偿方法 Download PDFInfo
- Publication number
- CN107436190B CN107436190B CN201610366092.4A CN201610366092A CN107436190B CN 107436190 B CN107436190 B CN 107436190B CN 201610366092 A CN201610366092 A CN 201610366092A CN 107436190 B CN107436190 B CN 107436190B
- Authority
- CN
- China
- Prior art keywords
- receptance function
- electroacoustic
- reciprocity device
- reciprocity
- compensation
- 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
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
技术领域
本发明涉及声学测量技术领域,特别涉及一种可用于互易测试的电声互易装置的非互易补偿方法。
背景技术
在任意声场环境中的低频声源测试是一个困难的问题,目前已有的方法往往存在低频极限。互易原理提供了一种解决该难题的途径。在基于互易原理的测量中,互易换能器是最重要的装置。
专利号为ZL 201520753800.0的中国专利《一种用于低频声源测试的双端面的力声互易装置》和申请号为201510622864.1的中国专利申请《一种用于低频声源测试的双端面的力声互易装置和方法》,采用了对称腔体结构来实现互易装置的设计。进一步地,申请号为201610286795.6以及申请号为201620390803.7的中国专利申请《一种用于水声低频声源测试的单端面互易装置》,针对双端面的力声互易装置存在对称腔体结构的加工和精度要求非常高的问题,提出了一种单端面的力声互易装置设计方法。
在理论上,许多电声装置是互易的,但是实际情况中存在非线性等因素的影响,会导致电声互易装置不互易。在现有技术中尚不存在对电声互易装置进行补偿的方法。
发明内容
本发明的目的在于克服现有技术中的某些电声互易装置在实际使用中不互易的缺陷,从而提供一种电声互易装置的非互易补偿方法。
为了达到上述目的,本发明提供了一种电声互易装置的非互易补偿方法,包括:
在全消声实验室内测量待测电声互易装置的声中心以及指向性,进而获得其表现为单极子特性的频率范围;随后在全消声室中,测量该电声互易装置的发射响应函数和接收响应函数;如果发射响应函数和接收响应函数相等,那么电声互易装置满足互易的要求;如果发射响应函数和接收响应函数不相等,那么该电声互易装置不满足互易的要求;如果发射响应函数和接收响应函数不相等并且不相等的频率范围位于单极子特性所在的频率范围之内,那么进行补偿,补偿的范围限于单极子特性所在的频率范围之内,补偿方法为:补偿量等于发射响应函数除以接收响应函数。
上述技术方案中,该方法进一步包括:
(1)在全消声实验室内,测得电声互易装置的声中心;
(2)在全消声实验室中,在电声互易装置前方不同角度测量该互易装置的指向性,确定其开始出现指向性的频率f0;
(3)在全消声实验室内,给电声互易装置输入电流I(f),在其正前方一定距离处放置传声器,测得声压pt(f),根据公式1计算该电声互易装置的发射响应函数,即单位电流激励下的体积速度:
其中,pt表示声压,r表示传声器距离电声互易装置的声中心的距离,I表示输入电流,ρ表示空气密度,f表示频率,绝对值表示响应函数的幅度;
(4)在全消声实验内,在距离电声互易装置一定距离处放置一个声源,在电声互易装置的声中心处放置传声器;让声源发声,测得传声器的声压,以及互易装置中扬声器的开路电压E(f);根据公式2计算该电声互易装置的接收响应函数:
其中,pr(f)表示声压,E(f)表示互易装置中扬声器的开路电压;
(5)比较电声互易装置的发射响应函数SG(f)和接收响应函数SR(f),如果传递函数完全相同,说明该电声互易装置具有较好的互易性;如果发射响应函数和接收响应函数不同,并且发射响应函数和接收响应函数不相等并且不相等的频率范围位于单极子特性所在的频率范围之内,需要进行非互易补偿;非互易补偿量为:
其中,Δ表示非互易补偿量。
本发明的优点在于:针对电声互易装置非互易的问题,本发明提供一种电声互易装置的非互易补偿方法,不仅提高了测试结果的准确信,同时降低了电声互易装置的获取难度。
附图说明
图1是本发明的电声互易装置的非互易补偿方法流程图;
图2(a)是基于电声互易原理的声源测试方法中正向实验的示意图;
图2(b)是基于电声互易原理的声源测试方法中反向实验的示意图;
图3是本发明的一个电声互易装置的发射响应函数和接收电流响应函数;
图4是本发明的一个电声互易装置的非互易补偿结果;
图5是本发明的基于电声互易装置的测试结果及经过非互易补偿的结果。
具体实施方式
现结合附图对本发明作进一步的描述。
如图1所示,本发明在全消声实验室内测量该电声互易装置的声中心以及指向性,进而获得其表现为单极子特性的频率范围;随后在全消声室中,测量该电声互易装置的发射性质(发射响应函数)和接收性质(接收响应函数)。如果发射响应函数和接收响应函数相等,那么电声互易装置满足互易的要求;如果发射响应函数和接收响应函数不相等,那么该电声互易装置不满足互易的要求。如果发射响应函数和接收响应函数不相等并且不相等的频率范围位于单极子特性所在的频率范围之内,那么进行补偿,补偿的范围限于单极子特性所在的频率范围之内,补偿方法为:补偿量等于发射响应函数除以接收响应函数。
本发明的非互易补偿方法具体包括以下步骤:
(1)在全消声实验室内,根据ISO3745:2003或者GB/T 6882-2008中的方法测得电声互易装置的声中心。
(2)在全消声实验室中,在电声互易装置前方不同角度测量该互易装置的指向性,确定其开始出现指向性的频率(f0)。
(3)在全消声实验室内,给电声互易装置输入电流I(f),在其正前方一定距离处放置传声器,测得声压pt(f),根据公式1计算该电声互易装置的发射响应函数,即单位电流激励下的体积速度:
其中,pt表示声压,r表示传声器距离电声互易装置的声中心的距离,I表示输入电流,ρ表示空气密度,f表示频率,绝对值表示响应函数的幅度。本步骤是根据单极子在自由场中的辐射特性推导而来,因此也被称为自由场法。
(4)在全消声实验内,在距离电声互易装置一定距离处放置一个声源,在电声互易装置的声中心处放置传声器。让声源发声,测得传声器的声压,以及互易装置中扬声器的开路电压E(f)。根据公式2即可以计算该电声互易装置的接收响应函数:
其中,pr(f)表示声压,E(f)表示互易装置中扬声器的开路电压。
(5)比较电声互易装置的发射响应函数SG(f)和接收响应函数SR(f),如果传递函数完全相同,说明该电声互易装置具有较好的互易性;如果发射响应函数和接收响应函数不同,并且发射响应函数和接收响应函数不相等并且不相等的频率范围位于单极子特性所在的频率范围之内,需要进行非互易补偿。非互易补偿量为:
Δ表示非互易补偿量。
该非互易补偿量的应用背景如下:
图2表示的是基于电声互易原理的声源测试方法示意图。图中S1表示待测声源,S2表示作为辅助传感器的电声互易装置。声场环境满足线性、无源、稳定的条件。
如图2(a)所示,正向实验中,给S1输入电流I1,测得S2的开路电压E1;如图2(b)所示,反向实验中,给S2输入电流I2,可测得S1的开路电压E2和传声器声压P2(传声器放置在S1的声中心处)。
为了以后的推导方便,对即将引入的物理量进行说明:
S1G:未知声源的发射响应,等于体积速度除以输入电流,S1G=Q1/I1;
S1R:未知声源的接收响应,等于声中心处的声压除以开路电压,S1R=E2/P2;
S2G:辅助扬声器的发射响应,等于体积速度除以输入电流,S2G=Q2/I2;
S2R:辅助扬声器的接收响应,等于声中心处的声压除以开路电压,S2R=E1/P1;
TF1:未知声源声中心到辅助扬声器声中心处的传递函数,等于辅助扬声器声中心处声压除以未知声源的体积速度,TF1=P1/Q1;
TF2:辅助扬声器声中心到未知声源声中心处的传递函数,等于未知声源声中心处声压除以辅助扬声器的体积速度,TF2=P2/Q2。
对于图2所示的电声互易实验系统,根据互易原理的测试方法可知:
根据图2(a)可知:
E1=Q1·TF1·S2R (公式5)
根据图2(b)可知:
P2=I2·S2G·TF2 (公式6)
定义体积速度响应函数
体积速度响应函数的单位为m3/sA,其物理意义为单位电流输入的体积速度响应。
结合公式(公式4)(公式5)(公式6)(公式7)可知:
如果声场和辅助扬声器S2都是互易的,即
那么(公式8)可以表示为:
显然,S1G和具有同样的物理意义,都表示待测声源S1的体积速度响应函数。但是S1G是S1的固有属性,而是通过电声互易实验测量得到的值,只有在理想测量情况下,两者才能完全相等。
如果辅助扬声器非互易,即S2G≠S2R,可知这说明,当辅助扬声器非互易时,利用电声互易原理测量的结果会出现误差。此时,(公式8)可简化为:
变形可知,电声实验的相对误差为:
同时可知:
其中:Δ=S2G/S2R,表示电声互易装置的发射响应函数除以接收响应函数。
(公式12)揭示,互易测试的结果需要经过电声互易装置的非互易补偿,才能得到正确的结果。在只考虑幅度的情况下,将(公式12)转化为dB形式,更为直观:
由(公式13)可以看出,非互易补偿量相当于在互易测试结果的基础上,增加一个修正量。
为了更好地说明本发明的非互易补偿方法,下面结合具体的实例做进一步说明。
电声互易装置一般由电动式扬声器和后腔组成,在一个实施例中,所采用的电声互易装置的后腔为有机玻璃材质,称为有机玻璃后腔电声互易装置,该装置的互易性质较差。下文以有机玻璃后腔电声互易装置为例结合本发明提供的方法进行说明。
在全消声实验室内,通过ISO3745:2003或者GB/T 6882-2008中的方法测量出有机玻璃后腔的电声互易装置的声中心的位置。经过测量,其声中心位置在其振膜中心前方3cm处。在全消声实验室中,在电声互易装置前方不同角度测量该互易装置的指向性,确定其开始出现明显指向性的频率。发现,该电声互易装置在600Hz以下具有良好的单极子特性。
在全消声实验室内,给电声互易装置输入电流,在其正前方2m处放置传声器,测得声压。根据公式1计算该电声互易装置的发射响应函数。注意此时公式1中的r的值应该是1.97m。发射响应函数如图3中实线所示。在全消声实验室内放置一个声源,在电声互易装置的声中心处放置传声器。让声源发声,测得传声器的声压,以及电声互易装置的开路电压,并根据公式2计算接收响应函数。接收响应函数如图3中虚线所示。
比较测量的发射响应函数和接收响应函数(图3中的实线和虚线),发现两条曲线不完全相等,尤其在低频有较大差异。这说明该电声互易装置并不互易。考虑到发射响应函数和接收响应函数差异较大的区域在350Hz以下,利用测量的发射响应函数和接收响应函数,根据公式3计算非互易补偿量,该非互易补偿量如图4所示。
如图5所示,是一个待测声源,利用有机后腔电声互易装置进行互易测试的结果,在该图中分别给出了不经过非互易补偿、经过非互易补偿的互易测试结果,以及参考方法——自由场法的测试结果。从中可以看出,经过非互易补偿之后,明显提高了测量的准确性。这从侧面说明了非互易补偿的效果。
最后应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种电声互易装置的非互易补偿方法,包括:
在全消声实验室内测量待测电声互易装置的声中心以及指向性,进而获得其表现为单极子特性的频率范围;随后在全消声室中,测量该电声互易装置的发射响应函数和接收响应函数;如果发射响应函数和接收响应函数相等,那么电声互易装置满足互易的要求;如果发射响应函数和接收响应函数不相等,那么该电声互易装置不满足互易的要求;如果发射响应函数和接收响应函数不相等并且不相等的频率范围位于单极子特性所在的频率范围之内,那么进行补偿,补偿的范围限于单极子特性所在的频率范围之内,补偿方法为:补偿量等于发射响应函数除以接收响应函数。
2.根据权利要求1所述的电声互易装置的非互易补偿方法,其特征在于,该方法进一步包括:
(1)在全消声实验室内,测得电声互易装置的声中心;
(2)在全消声实验室中,在电声互易装置前方不同角度测量该互易装置的指向性,确定其开始出现指向性的频率f0;
(3)在全消声实验室内,给电声互易装置输入电流I(f),在其正前方一定距离处放置传声器,测得声压pt(f),根据公式1计算该电声互易装置的发射响应函数,即单位电流激励下的体积速度:
其中,pt(f)表示声压,r表示传声器距离电声互易装置的声中心的距离,I(f)表示输入电流,ρ表示空气密度,f表示频率,绝对值表示响应函数的幅度;
(4)在全消声实验内,在距离电声互易装置一定距离处放置一个声源,在电声互易装置的声中心处放置传声器;让声源发声,测得传声器的声压,以及互易装置中扬声器的开路电压E(f);根据公式2计算该电声互易装置的接收响应函数:
其中,pr(f)表示声压,E(f)表示互易装置中扬声器的开路电压;
(5)比较电声互易装置的发射响应函数SG(f)和接收响应函数SR(f),如果传递函数完全相同,说明该电声互易装置具有较好的互易性;如果发射响应函数和接收响应函数不同,并且发射响应函数和接收响应函数不相等并且不相等的频率范围位于单极子特性所在的频率范围之内,需要进行非互易补偿;非互易补偿量为:
其中,Δ表示非互易补偿量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610366092.4A CN107436190B (zh) | 2016-05-27 | 2016-05-27 | 一种电声互易装置的非互易补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610366092.4A CN107436190B (zh) | 2016-05-27 | 2016-05-27 | 一种电声互易装置的非互易补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107436190A CN107436190A (zh) | 2017-12-05 |
CN107436190B true CN107436190B (zh) | 2019-10-22 |
Family
ID=60453312
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610366092.4A Active CN107436190B (zh) | 2016-05-27 | 2016-05-27 | 一种电声互易装置的非互易补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107436190B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109443514B (zh) * | 2018-09-13 | 2020-12-22 | 中国船舶重工集团公司第七一五研究所 | 一种声压灵敏度校准腔及其测试方法 |
CN111998934B (zh) * | 2020-08-28 | 2022-07-01 | 国网湖南省电力有限公司 | 一种声源声功率测试方法 |
CN112214838B (zh) * | 2020-10-10 | 2024-02-09 | 南京大学 | 一种基于互易定理的汽车低音单元最优布放方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4418116B2 (ja) * | 2001-02-01 | 2010-02-17 | 株式会社小松製作所 | 防音性能測定方法 |
CN102523057B (zh) * | 2011-12-27 | 2015-01-07 | 中国船舶重工集团公司第七一五研究所 | 一种低频声波自由场声压校准方法 |
JP2015077393A (ja) * | 2013-09-10 | 2015-04-23 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置、及び超音波測定方法 |
-
2016
- 2016-05-27 CN CN201610366092.4A patent/CN107436190B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107436190A (zh) | 2017-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107436190B (zh) | 一种电声互易装置的非互易补偿方法 | |
US10568543B2 (en) | Method and apparatus for testing vital capacity | |
US9368118B2 (en) | Voice analyzer, voice analysis system, and non-transitory computer readable medium storing program | |
WO2015109894A1 (zh) | 电涡流检测钞票厚度的方法以及装置 | |
CN102650545B (zh) | 一种高声压源结合相位控制的高声强级校准装置及方法 | |
EP2633275A2 (en) | Ultrasonic flow meter with zero impedance measuring electronics | |
KR101825363B1 (ko) | 진동막의 고유 주파수 획득 방법 및 시스템 | |
CN106568838B (zh) | 脉冲法声学材料性质测量装置、测量系统及测量方法 | |
US9949055B2 (en) | Signal control apparatus | |
US2416557A (en) | Electroacoustic transducer | |
CN105323694B (zh) | 密封验证方法、密封检测装置以及密封品质测量系统 | |
CN106813767B (zh) | 一种usp探头的灵敏度测量方法 | |
CN105578375B (zh) | 一种麦克风测试装置和测试方法 | |
CN208026374U (zh) | 在流动介质环境中实时测量静态声速的装置 | |
JP2006329718A (ja) | 膜スチフネス測定装置および膜スチフネス測定方法 | |
US9140619B2 (en) | Piezoelectric vacuum gauge and measuring method thereof | |
JPWO2015186177A1 (ja) | 液量計測装置 | |
CN202485795U (zh) | 一种高声压源结合相位控制的高声强级校准装置 | |
CN107621495A (zh) | 一种柔性薄膜固有频率的测试装置及方法 | |
WO2014187382A1 (zh) | 一种测距方法及设备 | |
Rosenfeld et al. | A modified Zisman apparatus for measuring contact potential differences in air | |
RU2563603C1 (ru) | Способ определения чувствительности по полю гидроакустического приемника | |
CN103234493B (zh) | 一种非线性系数自相关声学体积测量方法及系统 | |
RU157156U1 (ru) | Эффузионный детектор газов | |
US10859421B2 (en) | Electro acoustic volume measurement |
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 |