CN108718438B - 一种电声产品检测系统的响应信号时间起点的判断方法 - Google Patents

一种电声产品检测系统的响应信号时间起点的判断方法 Download PDF

Info

Publication number
CN108718438B
CN108718438B CN201810470005.9A CN201810470005A CN108718438B CN 108718438 B CN108718438 B CN 108718438B CN 201810470005 A CN201810470005 A CN 201810470005A CN 108718438 B CN108718438 B CN 108718438B
Authority
CN
China
Prior art keywords
response signal
starting point
signal
sliding
acoustic
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
CN201810470005.9A
Other languages
English (en)
Other versions
CN108718438A (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.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic Technology
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 Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN201810470005.9A priority Critical patent/CN108718438B/zh
Publication of CN108718438A publication Critical patent/CN108718438A/zh
Application granted granted Critical
Publication of CN108718438B publication Critical patent/CN108718438B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R29/00Monitoring arrangements; Testing arrangements

Landscapes

  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Otolaryngology (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及电声产品参数检测技术领域,具体涉及一种电声产品检测系统的响应信号时间起点的判断方法,对响应信号加窗处理,并查找加窗处理结果达到最大值时对应的矩形窗移动次数,从而获得电流响应信号和声响应信号的时间起点。本发明通过软件滑动窗处理实现电声产品响应信号时间起点判断,可准确得到时间起点,精度可根据实际情况调整滑动步长进行选择,适用于电声产品生产线使用;且本发明对电声产品的电流响应信号和声响应信号时间起点判断采用同一种方法,若忽略声传播时延,可仅计算电流响应时间起点,简化了软件处理流程,本发明计算简单、实用、有效、易实现,且不会显著增加系统运算量和存储量。

Description

一种电声产品检测系统的响应信号时间起点的判断方法
技术领域
本发明涉及电声产品参数检测技术领域,具体涉及一种电声产品检测系统的响应信号时间起点的判断方法。
背景技术
扬声器、传声器、麦克风、咪头等电声产品的参数检测是产品研发、生产、品质控制等环节不可或缺的技术保障。通过参数检测,可使研发人员了解设计缺点,使品质控制发现生产过程问题,调整生产环节。
现有电声产品的参数检测手段均是发送测试信号激励受话器,通过受话器的电流响应和声响应测试产品参数(齐娜,孟子厚.声频声学测量技术原理[M].北京:国防工业出版社,2008.7)。在此过程中,如何判断响应信号时间起点、准确快速高效的捕获响应信号是参数检测的基础,直接影响到检测性能和实用效果。值得注意的是,硬件检测电路和软件处理均存在短暂响应时延,使得发送和接收同时动作难以实现;同时受话器等电声器件本质上属于非线性换能系统,对激励信号同样存在响应时延,且随着激励频率不同时延也有所不同。这些都影响了响应信号捕获的准确性。
已有的电声产品检测系统的响应信号时间起点判断方法,通常是在发送正式测试信号之前添加单频等约定信号(郭庆,杨东奇,徐翠锋.扬声器电声参数自动测试系统的设计[J].科学技术与工程,2015,15(34):56-62),或是使用硬件电路触发响应信号采集(祁宇明,许增朴,王红星,胡立祥.扬声器激励与响应信号同步采集的一种设计方法[J].机械设计与制造,2010(6):185-187)。添加单频等约定信号实现时间起点判断,增加了测试的时间和计算开销,降低了测试效率;同时扬声器等电声器件本质是非线性换能系统,对频率信号存在时滞响应,且信号频率不同引起的时滞也不相同,根据单频信号并不能够准确判断响应信号起始时间。使用硬件电路触发信号采集增加了测试硬件的复杂度,引入了触发电路和采集电路的响应延时,时间起点判断精度不能保证。可见,现有时间起点判断方法不能够保证准确捕获响应信号,必须采用恰当的处理方法判断响应信号的时间起点。
发明内容
为了解决上述问题,本发明提供了一种电声产品检测系统的响应信号时间起点的判断方法,具体技术方案如下:
一种电声产品检测系统的响应信号时间起点的判断方法包括以下步骤:
(1)发送数字激励信号s(n)激励受话器,设数字激励信号s(n)的采样率为fs,持续时间为T,总点数为N;
(2)采用取样电阻R获得电流响应信号,采用送话器接收声响应信号;
对通过取样电阻获得的电流响应信号和送话器接收到的声响应信号以采样率为fs做同步采集,得到长度为M的电流响应信号c(n)和声响应信号r(n);
(3)将电流响应信号c(n)加长度为N点的滑动矩形窗
Figure BDA0001663054760000021
将声响应信号r(n)加长度为N点的滑动矩形窗
Figure BDA0001663054760000022
其中wci(n)和wri(n)偶对称于原点,i为滑动次数(i=0,1,2…),k为电流响应信号c(n)对应矩形窗滑动步进长度,m为声响应信号r(n)对应矩形窗滑动步进长度;
(4)将
Figure BDA0001663054760000023
Figure BDA0001663054760000024
分别与本地保存的数字激励信号s(n)逐点相乘并求和得yci、yri,求和可表示为:
yci=c(i*k)s(0)+c(i*k+1)s(1)+…+c(i*k+N-1)s(N-1);
yri=r(i*m)s(0)+r(i*m+1)s(1)+…+r(i*m+N-1)s(N-1);
(5)判断滑动矩形窗下降沿位置,若(i+1)*k+N-1≤M则转到步骤(3)继续对电流响应信号c(n)进行滑动处理,直至(i+1)*k+N-1>M;
若(i+1)*m+N-1≤M则转到步骤(3)继续对声响应信号r(n)进行滑动处理直至(i+1)*m+N-1>M;
(6)查找yci和yri达到最大值时对应的矩形窗移动次数a和b,则电流响应信号c(n)的时间起点tc和声响应信号r(n)的时间起点tr为:
tc=a*k;
tr=b*m。
进一步,所述步骤(2)中的同步采集是在数字激励信号s(n)开始发送前开始同步采集,在数字激励信号s(n)停止发送后停止同步采集。
本发明的有益效果为:
(1)时间起点判断所需的数字激励信号与电声产品性能参数测试信号可共用,不需增加额外硬件和测试时间开销;
(2)通过软件滑动窗处理实现电声产品响应信号时间起点判断,可准确得到时间起点,且精度可根据实际情况调整滑动步长进行选择,适用于电声产品生产线使用;
(3)对电声产品的电流响应信号和声响应信号时间起点判断采用同一种方法,若忽略声传播时延,可仅计算电流响应时间起点,简化了软件处理流程;
(4)本发明所述时间起点判断方法计算简单、实用、有效,不会显著增加系统运算量和存储量,工程易实现。
附图说明
图1为本发明的步骤流程图;
图2为本发明的实施例中发送数字激励信号以及采集响应信号的原理框图;
图3为发送数字激励信号s(n)时域图;
图4为采集到的电流响应信号c(n)时域图;
图5为采集到的声响应信号r(n)时域图;
图6为电流响应信号与发送信号的滑动处理求和结果yci示意图;
图7为声响应信号与发送信号的滑动处理求和结果yri示意图;
图8为时间起点判断后获得的电流响应信号;
图9为时间起点判断后获得的声响应信号;
图10为根据捕获的响应信号测量得到的扬声器频响曲线、总谐波失真曲线、相位曲线、阻抗曲线。
具体实施方式
为了更好的理解本发明,下面结合附图和具体实施例对本发明作进一步说明:
如图1所示,一种电声产品检测系统的响应信号时间起点的判断方法包括以下步骤:
(1)发送数字激励信号s(n)激励受话器,设数字激励信号s(n)的采样率为fs,持续时间为T,总点数为N;数字激励信号与电声产品性能参数测试信号共用,具体可采用正弦信号、白噪声信号、脉冲信号、线性正弦扫频信号、对数正弦扫频信号、指数正弦扫频信号、最长序列伪随机噪声中的一种。上述7种信号是电声测试国家标准(如:国家标准GB/T12060.5-2011等)和现有电声产品测试系统(齐娜,孟子厚.声频声学测量技术原理[M].北京:国防工业出版社,2008.7)常用的性能参数测试信号。本申请所用的数字激励信号与电声产品性能参数测试信号可共用,不需增加额外硬件和测试时间开销。
在本实施例中,数字激励信号s(n)采用对数正弦扫频信号,可表达为:
Figure BDA0001663054760000041
其中U为信号幅度,f1为扫频起始频率,f2为扫频终止频率,扫频速度
Figure BDA0001663054760000042
fs为信号采样频率,N为信号的总点数。
将数字激励信号s(n)经过功率放大器放大后送至受话器,如图2所示,取样电阻R与受话器串联。取fs为192000Hz,f1为20Hz,f2为20000Hz,U为0.089V,N为19200,即信号时长为0.1s,发送数字激励信号s(n)的时域图如图3所示。
(2)采用取样电阻R获得电流响应信号,采用送话器接收声响应信号;
对通过取样电阻R获得的电流响应信号和送话器接收到的声响应信号以采样率为fs通过采集设备做同步采集,得到长度为M的电流响应信号c(n)和声响应信号r(n)。
取样电阻R的电阻值影响采集到的电流响应信号的幅度,若取样电阻R太小,则影响电流响应信号c(n)的采集,而因为受话器不是一个纯阻抗元件,取样电阻R与受话器串联,若取样电阻取值过大则会影响受话器的输出功率。因此在本实施例中,选择取样电阻R为0.25欧姆,常用受话器的阻抗范围4-32欧姆,在本实施例中的测试中,受话器的电阻值为4欧姆,由于受话器的电阻值是取样电阻R的电阻值的15倍以上,所以可将取样电阻R等效为短路,不影响受话器的功率输出和测试结果。
同步采集是在数字激励信号s(n)开始发送前开始同步采集,在数字激励信号s(n)停止发送后停止同步采集,如此可得到完整的电流响应信号c(n)和声响应信号r(n)。采集到的电流响应信号c(n)和声响应信号r(n)分别如图4和图5所示。从图4和图5可以看出,由于取样电阻R较小,所以电流响应信号c(n)的幅度微弱;声响应信号r(n)失真较大,且随着频率不同失真明显不同。同时电流响应信号c(n)和声响应信号r(n)在有用信号前均存在一段空白,如若不去寻找响应信号时间起点而直接采用此信号测量电声器件参数,势必会对测量性能造成影响。
(3)将电流响应信号c(n)加长度为N=19200点的滑动矩形窗
Figure BDA0001663054760000043
将声响应信号r(n)加长度为N点的滑动矩形窗
Figure BDA0001663054760000051
其中wci(n)和wri(n)偶对称于原点,i为滑动次数(i=0,1,2…),k为电流响应信号c(n)对应矩形窗滑动步进长度,m为声响应信号r(n)对应矩形窗滑动步进长度;矩形窗滑动步进长度越小,时间起点判断精度越高,随之计算量也相应增加;可根据实际应用情况选择矩形窗滑动步进长度,若追求时间起点判断精度且计算能力足够,可取步进长度为1个采样点。另,wci(n)和wri(n)的滑动步进长度k和m可以相同。
(4)将
Figure BDA0001663054760000052
Figure BDA0001663054760000053
分别与本地保存的数字激励信号s(n)逐点相乘并求和得yci、yri,求和可表示为:
yci=c(i*k)s(0)+c(i*k+1)s(1)+…+c(i*k+N-1)s(N-1); (2)
yri=r(i*m)s(0)+r(i*m+1)s(1)+…+r(i*m+N-1)s(N-1); (3)
其中,
Figure BDA0001663054760000054
表示加矩形窗后的电流响应信号;
Figure BDA0001663054760000055
表示加矩形窗后的声响应信号;
取滑动步进长度k和m均为1个点,滑动处理后,yci、yri求和结果如图6和图7所示。其中,yci为矩形窗第i次滑动后加窗电流响应信号与本地保存的数字激励信号逐点相乘后求和的结果;yri为矩形窗第i次滑动后加窗声响应信号与本地保存的数字激励信号逐点相乘后求和的结果。
(5)判断滑动矩形窗下降沿位置,若(i+1)*k+N-1≤M则转到步骤(3)继续对电流响应信号c(n)进行滑动处理,直至(i+1)*k+N-1>M;
若(i+1)*m+N-1≤M则转到步骤(3)继续对声响应信号r(n)进行滑动处理直至(i+1)*m+N-1>M。
(6)查找yci和yri达到最大值时对应的矩形窗移动次数a和b,则电流响应信号c(n)的时间起点tc和声响应信号r(n)的时间起点tr为:
tc=a*k; (4)
tr=b*m; (5)
查找yci最大值对应滑动次数a=3851,yri最大值对应滑动次数b=3905,根据公式(4)和公式(5)计算得到电流响应信号c(n)的时间起点tc=3851,声响应信号r(n)的时间起点tr=3905。可以发现,计算得到的电流响应信号c(n)和声响应信号r(n)的起点不同,这是因为声响应信号r(n)存在声波传播延时,同时电声器件属于非线性换能系统,存在系统响应延时。实际应用若可忽略送话器到受话器之间的声音传播延时,则可只计算电流响应信号c(n)时间起点,并将电流响应信号c(n)的时间起点作为声响应信号r(n)的时间起点。
从起点处截取长度为N=19200点的信号,得到捕获后的电流响应信号c(n)和声响应信号r(n)分别如图8和图9所示。从图中可以看出按照本方法可以很准确的判断电流响应信号和声响应信号的时间起点。
图10为根据捕获的电流响应信号和声响应信号测量得到的扬声器频率响应曲线、总谐波失真曲线、相位曲线、阻抗曲线。从图中看出,采用本发明提供的判断电流响应信号和声响应信号的时间起点的方法,对扬声器的性能参数测试可以得到很准确的测试结果。
本发明不局限于以上所述的具体实施方式,以上所述仅为本发明的较佳实施案例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种电声产品检测系统的响应信号时间起点的判断方法,其特征在于:包括以下步骤:
(1)发送数字激励信号s(n)激励受话器,设数字激励信号s(n)的采样率为fs,持续时间为T,总点数为N;
(2)采用取样电阻R获得电流响应信号,采用送话器接收声响应信号;
对通过取样电阻获得的电流响应信号和送话器接收到的声响应信号以采样率为fs做同步采集,得到长度为M的电流响应信号c(n)和声响应信号r(n);所述同步采集是在数字激励信号s(n)开始发送前开始同步采集,在数字激励信号s(n)停止发送后停止同步采集;
(3)将电流响应信号c(n)加长度为N点的滑动矩形窗
Figure FDA0002385137000000011
将声响应信号r(n)加长度为N点的滑动矩形窗
Figure FDA0002385137000000012
其中wci(n)和wri(n)偶对称于原点,i为滑动次数, i=0,1,2…,k为电流响应信号c(n)对应矩形窗滑动步进长度,m为声响应信号r(n)对应矩形窗滑动步进长度;
(4)将
Figure FDA0002385137000000013
Figure FDA0002385137000000014
分别与本地保存的数字激励信号s(n)逐点相乘并求和得yci、yri,求和可表示为:
yci=c(i*k)s(0)+c(i*k+1)s(1)+…+c(i*k+N-1)s(N-1);
yri=r(i*m)s(0)+r(i*m+1)s(1)+…+r(i*m+N-1)s(N-1);
(5)判断滑动矩形窗下降沿位置,若(i+1)*k+N-1≤M则转到步骤(3)继续对电流响应信号c(n)进行滑动处理,直至(i+1)*k+N-1>M;
若(i+1)*m+N-1≤M则转到步骤(3)继续对声响应信号r(n)进行滑动处理直至(i+1)*m+N-1>M;
(6)查找yci和yri达到最大值时对应的矩形窗移动次数a和b,则电流响应信号c(n)的时间起点tc和声响应信号r(n)的时间起点tr为:
tc=a*k;
tr=b*m。
CN201810470005.9A 2018-05-16 2018-05-16 一种电声产品检测系统的响应信号时间起点的判断方法 Active CN108718438B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810470005.9A CN108718438B (zh) 2018-05-16 2018-05-16 一种电声产品检测系统的响应信号时间起点的判断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810470005.9A CN108718438B (zh) 2018-05-16 2018-05-16 一种电声产品检测系统的响应信号时间起点的判断方法

Publications (2)

Publication Number Publication Date
CN108718438A CN108718438A (zh) 2018-10-30
CN108718438B true CN108718438B (zh) 2020-05-01

Family

ID=63900062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810470005.9A Active CN108718438B (zh) 2018-05-16 2018-05-16 一种电声产品检测系统的响应信号时间起点的判断方法

Country Status (1)

Country Link
CN (1) CN108718438B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110519680B (zh) * 2019-10-28 2020-04-07 展讯通信(上海)有限公司 音频器件测试方法及装置
CN111049596B (zh) * 2019-12-18 2022-04-26 公安部第三研究所 一种对讲终端音频激活声压的测试方法及测试系统
CN112911483B (zh) * 2021-01-22 2022-08-09 苏州华兴源创科技股份有限公司 反馈信号有效数据起点的截取方法、设备及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105827550A (zh) * 2015-01-04 2016-08-03 电信科学技术研究院 一种利用滑动窗确定目标参数的方法及装置
CN107844289A (zh) * 2017-10-27 2018-03-27 周燕红 一种响度控制方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060078139A1 (en) * 2003-03-27 2006-04-13 Hilmar Meier Method for adapting a hearing device to a momentary acoustic surround situation and a hearing device system
US7522736B2 (en) * 2004-05-07 2009-04-21 Fuji Xerox Co., Ltd. Systems and methods for microphone localization
TW201406171A (zh) * 2012-07-27 2014-02-01 Primax Electronics Ltd 麥克風檢測方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105827550A (zh) * 2015-01-04 2016-08-03 电信科学技术研究院 一种利用滑动窗确定目标参数的方法及装置
CN107844289A (zh) * 2017-10-27 2018-03-27 周燕红 一种响度控制方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"扬声器电声参数自动测试系统的设计";郭庆 等;《科学技术与工程》;20151231;第15卷(第34期);全文 *

Also Published As

Publication number Publication date
CN108718438A (zh) 2018-10-30

Similar Documents

Publication Publication Date Title
CN108718438B (zh) 一种电声产品检测系统的响应信号时间起点的判断方法
CN106331951B (zh) 一种基于连续对数扫频的谐波失真测量方法及系统
CN101426169B (zh) 一种快速检测发声体声响应参数的时域跟踪滤波器及系统
CN105371945B (zh) 去除厅堂脉冲响应测量声信号中噪声的方法及装置
CN104050964A (zh) 音频信号还原度检测方法及系统
CN106768266B (zh) 一种剔除背景噪声能量的厂房混响时间的测量方法
WO2004054319A1 (ja) スピーカとマイクロホン間の音波伝搬時間測定方法およびその装置
CN109991590B (zh) 一种在有限空间压力罐内测试换能器低频发射特性的系统与方法
EP2562745A3 (en) Method and apparatus for impulse response measurement and simulation
US11976961B2 (en) Transducer transfer impedance calibration device based on spatial frequency domai smoothing technology
CN107785025B (zh) 基于房间脉冲响应重复测量的噪声去除方法及装置
CN106710602B (zh) 一种声学混响时间估计方法和装置
Policardi MLS and Sine-Sweep measurements
CN112414483A (zh) 一种超声波水表计量信号采集自适应方法
JP4960838B2 (ja) 距離測定装置、距離測定方法、距離測定プログラム、および記録媒体
CN106644049B (zh) 一种电容声的测量装置及测量方法
WO2019056760A1 (zh) 一种电子设备及其声波测距方法
CN104394028A (zh) 一种测量通信物理层信道频率幅度响应的方法和系统
Nguyen et al. Methods to Estimate the Channel Delay Profile and Doppler Spectrum of Shallow Underwater Acoustic Channels
Kwon et al. Time-selective windowing technique in free-field microphone reciprocity calibration
JP5611393B2 (ja) 遅延時間測定装置、遅延時間測定方法及びプログラム
Ma et al. Measurement of echo reduction for passive-material samples using sparse Bayesian learning and least squares estimation
CN107389803B (zh) 一种液体与固体延迟材料之间声反射系数的测量方法
NL2029967B1 (en) System for measuring objective acoustic quality parameters of hall
CN107219393B (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