CN102158793A - 一种使用激光传感器测量扬声器参数的方法及系统 - Google Patents
一种使用激光传感器测量扬声器参数的方法及系统 Download PDFInfo
- Publication number
- CN102158793A CN102158793A CN2011100838481A CN201110083848A CN102158793A CN 102158793 A CN102158793 A CN 102158793A CN 2011100838481 A CN2011100838481 A CN 2011100838481A CN 201110083848 A CN201110083848 A CN 201110083848A CN 102158793 A CN102158793 A CN 102158793A
- Authority
- CN
- China
- Prior art keywords
- signal
- module
- frequency
- digital
- loud speaker
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000004044 response Effects 0.000 claims abstract description 67
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 15
- 238000005086 pumping Methods 0.000 claims description 55
- 238000006243 chemical reaction Methods 0.000 claims description 36
- 238000005259 measurement Methods 0.000 claims description 28
- 238000012545 processing Methods 0.000 claims description 17
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 16
- 230000005284 excitation Effects 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000006872 improvement Effects 0.000 abstract description 2
- 238000012942 design verification Methods 0.000 abstract 1
- 238000001514 detection method Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 20
- 230000008859 change Effects 0.000 description 5
- 238000013461 design Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000978750 Havardia Species 0.000 description 1
- HBBOZFUQJDYASD-LPHOMBEVSA-N alpha-L-Fucp-(1->3)-[beta-D-Galp-(1->4)]-beta-D-GlcpNAc Chemical compound O[C@H]1[C@H](O)[C@H](O)[C@H](C)O[C@H]1O[C@H]1[C@H](O[C@H]2[C@@H]([C@@H](O)[C@@H](O)[C@@H](CO)O2)O)[C@@H](CO)O[C@@H](O)[C@@H]1NC(C)=O HBBOZFUQJDYASD-LPHOMBEVSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提出一种使用激光位移传感器测量扬声器参数的方法，并给出了实现该方法的测试系统。该方法采用步进扫频信号激励被测器件，同时测量了被测器件振膜的位移响应信号和流过被测器件的电流响应信号。将上述信号代入扬声器的线性模型和大信号沃特拉模型，通过系统辨识得到扬声器的线性参数和非线性参数。由于扬声器的线性参数和非线性参数可以表征扬声器在不同幅度输入下的工作特性，该方法不但可以应用于扬声器的设计验证、产品质量检测，还可以应用于扬声器产品的音质改善。
Description
技术领域
本发明属于扬声器的参数测量方法及测量系统，具体地说是通过对振膜振动幅度和负载电流进行系统辨识，获得电声产品多个参数的方法及系统。可广泛应用于电声企业的扬声器研发及产品质量检测。
背景技术
扬声器参数测量的意义：
扬声器参数对扬声器的设计、制造、和质量控制具有重要的作用。自1970年以来，澳大利亚的A. N. Thiele和R. H. Small提出了扬声器系统的T/S参数，使通过滤波器综合法实现扬声器系统的设计成为可能。其中的小信号参数是在小信号激励、系统输出中失真可以忽略的情况下测得，故又称线性参数。该参数可用于表征扬声器的低频性能，应于扬声器的设计和生产的质量控制。
1980年以来，研究人员开始对扬声器的大信号性能进行各种量化分析。许多模型被用于扬声器的非线性特性分析，如沃特拉模型(Volterra)、状态空间模型、有限元及边界元模型、Hammerstein与Wiener模型等等。通过对扬声器建模和大信号参数测量，可以建立被测扬声器的计算机模型，在扬声器的设计阶段可以了解该设计的预期电声性能，在小批量样品试制阶段可以测量样品的参数以验证和修改设计，在量产阶段可以用于故障分析和工艺改进。
扬声器的参数测量，特别是非线性参数的测量是电声测量界的难题，现有的线性参数和非线性参数的测量技术和方法如下：
1)通过改变振膜质量、腔体体积的方法：测量扬声器的阻抗曲线，从而得到f s 、Q ts 等参数，之后利用附加质量或附加容积法再次测量阻抗曲线，从而得到f s ’、Q ts ’等，并计算得到所有的线性参数。这种方法缺点为：1、一般需要二次测量，操作麻烦；2、附加质量或附加容积法在一定程度上破坏了扬声器原有工作状态，故结果不准确；3、二次测量导致振膜顺性、音圈温度发生变化，对测量结果有影响。
2) 通过测量声压、电流响应的方法：美国专利6269318 (Method For determining Transducer linear operation parameters)提出了使用电压－声压传递函数和电流－声压传递函数测量扬声器线性参数的方法。该方法使用的传声器易受环境噪声、环境温湿度影响，同时该方法需要提供一个已知容积的测试箱体作为准确测量的条件。
3) 通过改变平衡位置测量：已公开发明专利CN200510041525.0(基于系统辨识的扬声器单元非线性参数测量方法)使用馈给可调直流电流的方法使音圈偏移一定的位置，在不同音圈位移下测量阻抗并根据扬声器单元非线性电阻抗模型辨识得到非线性参数。该方法需要使用可调直流源向扬声器输出偏置电流，而该电流会使音圈温度上升。此外由于振膜的平衡位置改变，扬声器工作在特殊的状态。亦有学者提出使用腔体和气泵改变振膜平衡位置的方法，但所需设备较为特殊，同时也存在改变扬声器正常工作状态的问题。
发明内容
本发明针对上述不足，提出一种使用激光传感器测量扬声器线性及非线性参数的方法及系统。
本发明提出的使用激光传感器测量扬声器线性及非线性参数的方法包括以下步骤：
1）设于计算机中的激励信号生成模块生成激励信号，该激励信号为数字信号，该数字信号经数模模数转换模块转换为模拟信号后输出；
2）功率放大器将数模模数转换模块输出的模拟信号放大后输出至被测器件；
3）激光传感器测量被测器件的位移响应信号，电流传感器测量流过被测器件的电流响应信号，得到的电流信号和位移信号输入数模模数转换模块，通过数模模数转换模块转换为数字信号后送入响应信号采集及处理模块，再把经处理后的信号输入到线性参数辨识模块和非线性参数辨识模块，经线性参数辨识模块的线性参数辨识模型辨识处理和非线性参数辨识模块的大信号沃特拉模型辨识处理，得到扬声器的线性参数和非线性参数。
扬声器的线性参数包括R e、l 0、L ces 、R es 和C mes ；扬声器的非线性参数数包括Bl(x)、
K(x)和Le(x)，其中：
Bl(x) = b 0 + b 1 x + b 2 x 2
K(x) = k 0 + k 1 x + k 2 x 2
L e (x) = l 0 + l 1 x + l 2 x 2 。
被测器件包括动圈扬声器单元、或扬声器闭箱系统、或扬声器倒相箱系统。
由于扬声器在大信号激励时可视为非线性器件，使用2个以上的频率激励扬声器，将产生多个谐波和互调响应项，这些响应有可能在频率上重合，导致测量失效，故本发明的激励信号选用单频或双频信号。
单频激励信号形式有多种，本发明采用如式（1）的步进对数扫频信号：
其中，U为激励信号幅度。w i为激励频率，w i+1一般按oct =1/3、1/6、1/12、1/24及其他倍频程标准或按对数关系选取，即w i+1/w i =2 oct 。t i+1 -t i 为激励信号某一激励频率w i 的持续时间。为保证两个频率切换时幅度不产生冲击地平缓过渡，需要考虑相位的连续性。故需要根据上一个信号的终止相位 u(t)为阶跃信号。
若采用双频信号，其形式为：
其中，第1项为幅度U 1固定频率为w 0 的激励信号，用于提供足够的振膜位移，称为位移激励信号。第2项为幅度U 1的步进对数扫频信号。
当流过被测器件的电流及其两端的电压已知时，在已知被测器件的等效振动质量M ms 的前提下，通过阻抗特性Z(s)公式，辨识得到线性参数；通过Volterra核函数关系式，采用非线性最小二乘法辨识得到非线性参数；
阻抗特性Z(s)为：
Volterra核函数关系式为：
数模模数转换模块将激励信号输出至功率放大器的同时，响应信号采集及处理模块采集扬声器的电流及位移瞬时值并将采集到的信号按激励信号的频率拆分，并做快速傅立叶变换，得到各个激励频率对应的基频响应以及谐波或互调响应。
采用单频信号时需测量基频的2倍频率和3倍频率处的谱线，得到2次谐波、3次谐波的幅度和相位。采用双频信号时需测量扫频信号为中心频率的2阶和3阶互调响应。
本发明所述的一种使用电流传感器测量扬声器参数的系统包括计算机和测试支架。计算机与数模模数转换模块连接、数模模数转换模块与功率放大器、激光传感器、电流传感器连接，激光传感器安装于测试支架上。
计算机内设有生成激励信号的激励信号生成模块、响应信号采集及处理模块、线性参数辨识模块、非线性参数辨识模块、用户接口模块和校准模块，激励信号生成模块、响应信号采集及处理模块和数模模数转换模块连接，用户接口模块和激励信号生成模块、线性参数辨识模块、非线性参数辨识模块连接，线性参数辨识模块、非线性参数辨识模块和响应信号采集及处理模块连接，校准模块和激励信号生成模块连接。
本发明优点是：首先，激光测量是非接触测量，无需改变振膜的平衡位置、质量等，扬声器工作在正常状态；其次，激光测量不受噪声影响，激光传感器受温湿度影响小。目前，激光三角法型位移传感器精度可以满足测量需求，不受被测物的材料、质地、 形状、反射率的限制，售价和传声器相当，因此本发明可以满足扬声器的生产线测试需求。
附图说明
图1为使用激光传感器测量扬声器参数的方法模块图。
图2为使用激光传感器测量扬声器参数的系统装置模块图。
图3为扬声器的大信号机电等效电路图。
图4为扬声器的大信号Volterra模型方框图。
图5为实施例被测扬声器的位移基频响应、2次谐波及3次谐波响应幅度图。
图6为实施例被测扬声器的位移基频响应、2次谐波及3次谐波响应相位图。
图7为实施例被测扬声器的阻抗特性测量结果和拟合结果对比图。
图8为实施例辨识结果与实验结果对比图。
具体实施方式
下面结合附图及实施例，对本发明作进一步的说明。
如图所示，本发明提出的使用激光传感器测量扬声器参数的方法包括以下步骤：
1）设于计算机7中的激励信号生成模块1生成激励信号，该激励信号为数字信号，该数字信号经数模模数转换模块8转换为模拟信号后输出；
2）功率放大器9将数模模数转换模块8输出的模拟信号放大后输出至被测器件13；
3）激光传感器10测量被测器件13的位移响应信号，电流传感器11测量流过被测器件的电流响应信号，得到的电流信号和位移信号输入数模模数转换模块8，通过数模模数转换模块8转换为数字信号后送入响应信号采集及处理模块2，再把经处理后的信号输入到线性参数辨识模块3和非线性参数辨识模块4，经线性参数辨识模块3的线性参数辨识模型辨识处理和非线性参数辨识模块4的大信号沃特拉模型辨识处理，得到扬声器的线性参数和非线性参数。
扬声器的线性参数包括R e、l 0、L ces 、R es 和C mes ；扬声器的非线性参数数包括Bl(x)、
K(x)和Le(x)，其中：
Bl(x) = b 0 + b 1 x + b 2 x 2
K(x) = k 0 + k 1 x + k 2 x 2
L e (x) = l 0 + l 1 x + l 2 x 2 。
被测器件13包括动圈扬声器单元、或扬声器闭箱系统、或扬声器倒相箱系统。
由于扬声器在大信号激励时可视为非线性器件，使用2个以上的频率激励扬声器，将产生多个谐波和互调响应项，这些响应有可能在频率上重合，导致测量失效，故本发明的激励信号选用单频或双频信号。
单频激励信号形式有多种，本发明采用如式（1）的步进对数扫频信号：
其中，U为激励信号幅度。w i为激励频率，w i+1一般按oct =1/3、1/6、1/12、1/24及其他倍频程标准或按对数关系选取，即w i+1/w i =2 oct 。t i+1 -t i 为激励信号某一激励频率w i 的持续时间。为保证两个频率切换时幅度不产生冲击地平缓过渡，需要考虑相位的连续性。故需要根据上一个信号的终止相位
u(t)为阶跃信号。
若采用双频信号，其形式为：
其中，第1项为幅度U 1固定频率为w 0 的激励信号，用于提供足够的振膜位移，称为位移激励信号。第2项为幅度U 1的步进对数扫频信号。
由于位移的基频响应在信号频率大于机械共振频率后，以每倍频程12dB的速度衰减，故测量双频信号的互调可获得较谐波更高的信噪比。然而双频信号也有缺点。一方面，由于基频响应谐波的存在，扫频信号的起始频率至少需要大于位移激励信号频率的3倍；另一方面，位移激励信号的频率不可能太低，否则互调频率与中心频率太过接近，受到傅立叶变换中频率分辨率的限制，无法准确取得互调响应。由于扬声器的小信号模型和大信号Volterra模型适用于扬声器的低频段，故对于共振频率前的测试频率点更为关注。考虑到上述问题，应根据被测扬声器的共振频率点选取激励信号。
当流过被测器件13的电流及其两端的电压已知时，在已知被测器件的等效振动质量M ms 的前提下，通过阻抗特性Z(s)公式，辨识得到线性参数；通过Volterra核函数关系式，采用非线性最小二乘法辨识得到非线性参数；
阻抗特性Z(s)为：
Volterra核函数关系式为：
数模模数转换模块8将激励信号输出至功率放大器9的同时，响应信号采集及处理模块2采集扬声器的电流及位移瞬时值并将采集到的信号按激励信号的频率拆分，并做快速傅立叶变换，得到各个激励频率对应的基频响应以及谐波或互调响应。
采用单频信号时需测量基频的2倍频率和3倍频率处的谱线，得到2次谐波、3次谐波的幅度和相位。采用双频信号时需测量扫频信号为中心频率的2阶和3阶互调响应。
线性参数辨识模块3建立线性参数辨识模型
线性参数辨识模块3建立在扬声器集中参数机电类比等效电路的基础上。由等效电路可建立如下微分方程：
当扬声器工作在小振幅状态时，对式（3）与式（4）进行Laplace变换，可写为
由于振膜位移x很小，Bl(x)、Le(x)和K(x)可视为常数，分别记为b 0、l 0和k 0，式（5）和式（6）中的dLe(x)/dx项为0。数模模数转换模块的增益和功率放大器的增益已知，故可由激励信号幅度计算出扬声器两端电压值U(s)。采集激光传感器得到电流响应I(s)，可计算得到扬声器的阻抗特性Z(s)和电压-位移线性传递函数H u1 (s)为
测量得到的阻抗特性即功率放大器9输出电压比流经扬声器的电流的模值。使用非线性最小二乘拟合，可以求得R e、l 0、k 0/b 0 2、R ms/b 0 2和M ms/b 0 2共5个值。由位移线性传递函数的辨识结果，可以求出b 0 。或者，M ms是一个不易受环境变化影响的值，若事先已知，可以求出k 0、R ms和b 0。
其中后3个值也可以分别表示为3个电等效参数：
L ces = b 0 2/k 0
R es =b 0 2/R ms
C mes =M ms/b 0 2
非线性参数辨识模块4建立非线性参数辨识模型
当扬声器工作在大振幅状态时，由于参数Bl(x)、Le(x)和K(x)随着振膜位移x的变化而变化，其输入与输出不再是线性关系。通常使用Volterra级数求解式（3）与式（4）组成的方程组，得到输入电压与输出位移的关系。根据Volterra级数理论，非线性系统的响应可以表示为各阶核函数与激励信号的卷积的形式，如式（9），其框图如图4所示。
式中，u(t)为输入，h 1表示系统的1阶核函数即线性冲激响应，h 2与h 3分别表示系统的2阶和3阶核函数。
当输入信号x(t) = exp(s 1 t)时，系统线性响应表达式为
y(t) = H 1(s1)exp(s 1 t) (10)
除线性响应外，输出中还具有输入信号的2次以上谐波响应。
输入信号x(t) = exp(s 1 t)+exp(s 2 t)时，只考虑到2阶响应输出，有
y (t) = H 1 (s1) exp(s 1 t) + H 1(s2) exp(s 2 t) +
H 2 (s 1,s 1) exp(2s 1 t) + H 2(s 2,s 2) exp(2s 2 t) + 2H 2(s 1,s 2) exp(s 1 t+s 2 t) (11)
式中前两项为线性响应。第3、4项对应着系统的2次谐波失真。H 2(s 1,s 2)是对称的，有H 2(s 1,s 2) = H 2(s 2,s 1)，故第5项系数为2，对应着系统的互调失真。
输入信号x(t) = exp(s 1 t)+exp(s 2 t)+exp(s 3 t)时，只考虑到3阶响应输出，有
y(t) = y 1(t) + y 2(t) + y 3(t)+…… (12)
其中，y 1(t)是1阶核函数的响应即线性响应。y 2(t)与y 3(t)分别是2、3阶核函数的响应，由谐波响应和互调响应构成。其具体形式如下：
y 1(t) = H 1(s 1) exp(s 1 t) + H 1(s 2) exp(s 2 t) + H 1(s 3) exp(s 3 t) (13)
y 2(t) = H 2(s 1,s 1) exp(2s 1 t)+ H 2(s 2,s 2) exp(2s 2 t) + H 2(s 3,s 3) exp(2s 3 t) +
2H 2(s 1,s 2)exp(s 1 t+s 2 t) + 2H 2(s 1,s 3)exp(s 1 t+s 3 t) + 2H 2(s 2,s 3) exp(s 2 t+s 3 t) (14)
y 3(t) = H 3(s 1,s 1,s 1) exp(3s 1 t) + H 3(s 2,s 2,s 2) exp(3s 2 t) + H 3(s 3,s 3,s 3) exp(3s 3 t) +
3H 3(s 1,s 1,s 2)exp(2s 1 t+s 2 t) +3H 3(s 1,s 1,s 3)exp(2s 1 t+s 3 t) + 3H 3(s 2,s 2,s 1)exp(s 1 +2s 2 t) +
3H 3(s 2,s 2,s 3)exp(2s 2 t+s 3 t) +3H 3(s 3,s 3,s 1)exp(s 1 +2s 3 t) + 3H 3(s 3,s 3,s 2)exp(s 2 t+2s 3 t) +
6H 3(s 1,s 2,s 3) exp(s 1 t+s 2 t+s 3 t) (15)
式（14）的前3项为2次谐波响应，后3项为2次互调响应。与之类似，3次谐波响应参见y 3(t)的前3项，其余项为3次互调响应项。可见，随着分析阶数的增加，互调响应项的数量显著增加。为了便于阐述，说明书中仅对扬声器进行1～3阶Volterra建模。
将Bl(x)、K(x)与L e (x)近似表示为2阶幂级数的形式，如下：
Bl(x) = b 0 + b 1 x + b 2 x 2 (16)
K(x) = k0 + k1x + k2x2 (17)
Le(x) = l0 + l1x + l2x2 (18)
令输入U = exp(s 1 t) + exp(s 2 t) + exp(s 3 t)，且位移为如式（12）的形式，代入式（5）和式（6），由谐波平衡法可求出线性核函数如式（8），及2、3阶核函数如下
。
由上述推导可以看出，如果已知扬声器的线性参数，且已知谐波或互调响应，可以使用最小二乘非线性拟合法进行曲线拟合，从而求出非线性参数的幂级数系数b 1、k 1、l 1和b 2、k 2、l 2。
本发明所述的一种使用电流传感器测量扬声器参数的系统包括计算机7和测试支架12。计算机7与数模模数转换模块8连接、数模模数转换模块8与功率放大器9、激光传感器10、电流传感器11连接，激光传感器10安装于测试支架12上。
计算机7内设有生成激励信号的激励信号生成模块1、响应信号采集及处理模块2、线性参数辨识模块3、非线性参数辨识模块4、用户接口模块6和校准模块5，激励信号生成模块1、响应信号采集及处理模块2和数模模数转换模块8连接，用户接口模块6和激励信号生成模块1、线性参数辨识模块3、非线性参数辨识模块4连接，线性参数辨识模块3、非线性参数辨识模块4和响应信号采集及处理模块2连接，校准模块5和激励信号生成模块1连接。
本发明的数模模数转换模块8是已被校准的部件，故可以用于对系统中功率放大器9、激光传感器10、电流传感器11的校准。
在测量扬声器参数前，校准模块控制数模模数转换模块8发出激励信号，并测量功率放大器9两端带负载情况下的电压相应，从而计算出功率放大器9的频响和失真。通过测量测试支架12上的已知规格的梯形台阶的数个距离对应的电压值，可由直线拟合计算得到位移—电压比，实现激光传感器的校准。将已知阻值的校准电阻接入功率放大器的输出端，由输出电压比校准电阻值得到输出回路电流，再由电流传感器11两端电压比输出回路电流得到电流传感器的等效电阻值。
本发明的用户接口模块6功能是为用户接口模块提供测试界面、为用户输入参数判断及保存、为测量结果显示等。
计算机7中运行着包含本发明所述测试方法的软件，控制数模模数转换模块8，实现校准和测量、参数辨识和结果显示等功能。
数模模数转换模块8实现了数字信号和模拟信号的相互转换。计算机7中的激励信号生成模块1生成的激励信号为数字信号，经该模块转换后输出为模拟信号。
功率放大器9将数模模数转换模块8的模拟信号放大，输出至被测扬声器。功率放大器的频响、失真等参数可由校准模块实现校准。功率放大器的输出回路中串接着电流传感器，电流传感器11两端的电压采用四线测量法接入数模模数转换模块8。
激光传感器10使用的是激光三角法测量位移。用一束激光以某角度聚焦在扬声器振膜表面，物体表面上的激光光斑成像于光电传感器上。振膜位移不同，表面光斑的位置高度不同，所接收散射或反射光线的角度也不同。光电传感器可测出光斑像的位置，计算出主光线的角度，从而计算出振膜的瞬时位移。
测试支架12用于激光传感器10和被测扬声器13的固定。通过调整扬声器的水平位置和激光头的垂直位置，可以把激光光束打在防尘帽的中心位置，使光束与振动表面垂直以获得最佳的信噪比。
现按本发明各模块的顺序说明其测试方法：
实施例选用的激励信号为单频信号。测试前，用户需要向用户接口模块6提供被测产品的基本参数，包括激励信号的电压幅度U、扫频的起止频率f 1、f 2及扫频频率步进。比如U=1V，f 1=20Hz，f 2=20kHz，扫频频率按1/3倍频程关系由起始频率递增。
上述激励信号由数模模数转换模块8输出至功率放大器9。功率放大器9推动被测扬声器13发声。首次测量时，校准模块5控制数模模数转换模块8采集功率放大器9两端带负载情况下的电压相应，用于得到功放的频响。之后的测量中，数模模数转换模块8采集位移传感器10及电流传感器11的输出，并存储到计算机6的内存中。通过计算，得到各个激励频率对应的基频响应、2次谐波和3次谐波响应的幅度和相位，如图5和图6。
流过被测扬声器13的电流及其两端的电压、振膜位移已知，故可以求出电压-位移传递函数H u1 (s)和阻抗特性Z(s)，之后辨识得到线性参数。类似地，可以辨识得到非线性参数。
为了验证辨识结果，上述参数被代回扬声器的微分方程中。通过数值方法，可计算出另一电压下扬声器的线性响应、及谐波失真，并且将其与测量结果进行比较。为了衡量预测结果与测量结果的一致性，可使如下用误差函数
(21)
其中w表示权重函数，形式为
使用激光辨识非线性参数，并用于预测位移线性响应、2阶响应和3阶响应的误差分别为6.2%，11.6%和11.9%。
实施例中，数模模数转换模块8支持4通道输入、4通道输出，采样位数24位，采样率最高可达192kHz。功率放大器9输出功率可达100W，失真小于0.001%。激光传感器10的量程40mm，精度4mm，测量位移的频率范围为直流～1000Hz。电流传感器11为一只测量用的4线功率电阻，阻值0.100Ω。被测扬声器13为直径6.5寸的汽车扬声器。
以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案；因此，尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明，但是，本领域的普通技术人员应当理解，仍然可以对本发明进行修改或等同替换；而一切不脱离发明的精神和范围的技术方案及其改进，其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种使用激光传感器测量扬声器参数的方法，其特征在于其特征在于该方法包括如下步骤：
1）设于计算机中的激励信号生成模块生成激励信号，该激励信号为数字信号，该数字信号经数模模数转换模块转换为模拟信号后输出；
2）功率放大器将数模模数转换模块输出的模拟信号放大后输出至被测器件；
3）激光传感器测量被测器件的位移响应信号，电流传感器测量流过被测器件的电流响应信号，得到的电流信号和位移信号输入数模模数转换模块，通过数模模数转换模块转换为数字信号后送入响应信号采集及处理模块，再把经处理后的信号输入到线性参数辨识模块和非线性参数辨识模块，经线性参数辨识模块的线性参数辨识模型辨识处理和非线性参数辨识模块的大信号沃特拉模型辨识处理，得到扬声器的线性参数和非线性参数。
2.根据权利要求1所述的使用激光传感器测量扬声器参数的方法，其特征在于扬声器的线性参数包括R e、l 0、L ces 、R es 和C mes ；扬声器的非线性参数线性参数包括Bl(x)、K(x)和Le(x)，其中：
Bl(x) = b 0 + b 1 x + b 2 x 2
K(x) = k 0 + k 1 x + k 2 x 2
L e (x) = l 0 + l 1 x + l 2 x 2 。
3.根据权利要求1所述的使用电流传感器测量扬声器参数的方法，其特征在于被测器件包括动圈扬声器单元、或扬声器闭箱系统、或扬声器倒相箱系统。
4.根据权利要求1所述的使用电流传感器测量扬声器参数的方法，其特征在于所述的激励信号选用单频或双频信号；
若激励信号采用单频信号，即采用步进对数扫频信号，其形式为：
其中，U为激励信号幅度，w i 为激励频率，w i+1一般按oct =1/3、1/6、1/12、1/24及其他倍频程标准或按对数关系选取，即w i+1/w i =2 oct ，t i+1 -t i 为激励信号某一激励频率w i 的持续时间，每个信号的初相为f i ，u(t)为阶跃信号；
若激励信号采用双频信其形式为
其中，第1项为幅度U 1固定频率为w 0 的激励信号，用于提供足够的振膜位移，称为位移激励信号，第2项为幅度U 1的步进对数扫频信号；
若激励信号采用单频信号，需测量基频的2倍频率和3倍频率处的谱线，得到2次谐波、3次谐波的幅度和相位；
若激励信号采用双频信号，需测量扫频信号为中心频率的2阶和3阶互调响应。
5.根据权利要求1所述的使用电流传感器测量扬声器参数的方法，其特征在于当流过被测器件的电流及其两端的电压已知时，在已知被测器件的等效振动质量M ms 的前提下，通过阻抗特性Z(s)公式，辨识得到线性参数；通过Volterra核函数关系式，采用非线性最小二乘法辨识得到非线性参数；
阻抗特性Z(s)为：
Volterra核函数关系式为：
。
6.根据权利要求1所述的一种使用电流传感器测量扬声器参数的系统，包括计算机和测试支架，其特征在于计算机与数模模数转换模块连接、数模模数转换模块与功率放大器、激光传感器、电流传感器连接，激光传感器安装于测试支架上。
7.根据权利要求1所述的一种使用电流传感器测量扬声器参数的系统，其特征在计算机内设有生成激励信号的激励信号生成模块、响应信号采集及处理模块、线性参数辨识模块、非线性参数辨识模块、用户接口模块和校准模块，激励信号生成模块、响应信号采集及处理模块和数模模数转换模块连接，用户接口模块和激励信号生成模块、线性参数辨识模块、非线性参数辨识模块连接，线性参数辨识模块、非线性参数辨识模块和响应信号采集及处理模块连接，校准模块和激励信号生成模块连接。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110083848 CN102158793B (zh) | 2011-04-02 | 2011-04-02 | 一种使用激光传感器测量扬声器参数的方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110083848 CN102158793B (zh) | 2011-04-02 | 2011-04-02 | 一种使用激光传感器测量扬声器参数的方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102158793A true CN102158793A (zh) | 2011-08-17 |
| CN102158793B CN102158793B (zh) | 2013-12-18 |
Family
ID=44439920
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201110083848 Expired - Fee Related CN102158793B (zh) | 2011-04-02 | 2011-04-02 | 一种使用激光传感器测量扬声器参数的方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102158793B (zh) |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102968543A (zh) * | 2012-12-13 | 2013-03-13 | 嘉善恩益迪电声技术服务有限公司 | 一种扬声器音圈及磁路的温度特性数值模拟方法 |
| CN103841505A (zh) * | 2014-02-21 | 2014-06-04 | 歌尔声学股份有限公司 | 声学产品的ccd声阻测试方法及系统 |
| CN104640051A (zh) * | 2013-11-06 | 2015-05-20 | 亚德诺半导体股份有限公司 | 估计扬声器的振膜偏移的方法 |
| TWI508576B (zh) * | 2013-05-15 | 2015-11-11 | Lite On Opto Technology Changzhou Co Ltd | 揚聲器異音檢測方法及裝置 |
| CN105101032A (zh) * | 2015-08-27 | 2015-11-25 | 深圳市音沃仕科技有限公司 | 一种电声转换器复合天线模组的声学测试的方法 |
| CN106068007A (zh) * | 2016-06-07 | 2016-11-02 | 瑞声科技(新加坡)有限公司 | 扬声器非线性系统辨识方法 |
| CN106303882A (zh) * | 2016-08-13 | 2017-01-04 | 厦门傅里叶电子有限公司 | 定义扬声器工作极限的方法 |
| CN106454679A (zh) * | 2016-11-17 | 2017-02-22 | 矽力杰半导体技术(杭州)有限公司 | 扬声器振膜状态估计方法及应用其的扬声器驱动电路 |
| CN106851514A (zh) * | 2017-02-25 | 2017-06-13 | 中山市天键电声有限公司 | 扬声器非线性参数测试装置系统及测试方法 |
| CN104170404B (zh) * | 2012-03-05 | 2018-01-26 | 奥音科技（北京）有限公司 | 音频系统、具有电‑声换能器的设备、测量装置及方法 |
| CN108614431A (zh) * | 2018-06-08 | 2018-10-02 | 河海大学常州校区 | 一种基于夹角的Hammerstein-Wiener系统多模型分解及控制方法 |
| CN109655735A (zh) * | 2018-12-07 | 2019-04-19 | 武汉市聚芯微电子有限责任公司 | 功放芯片评估板以及功放芯片评估板系统 |
| CN109738797A (zh) * | 2019-01-16 | 2019-05-10 | 珠海市精实测控技术有限公司 | 一种基于大电压模型的线性马达bl方法 |
| CN110022522A (zh) * | 2019-01-29 | 2019-07-16 | 浙江中科电声研发中心 | 采用激振器激励的扬声器振动部件共振频率测量系统及测量方法 |
| CN110769350A (zh) * | 2019-09-29 | 2020-02-07 | 惠州市杰德创新科技有限公司 | 喇叭响度的调节方法及喇叭 |
| CN111448804A (zh) * | 2018-01-17 | 2020-07-24 | 三星电子株式会社 | 用于扬声器驱动器的运动的非线性控制的方法和系统 |
| CN115031863A (zh) * | 2022-08-11 | 2022-09-09 | 成都凯天电子股份有限公司 | 一种阻滞式总温传感器 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1761366A (zh) * | 2005-08-18 | 2006-04-19 | 南京大学 | 基于系统辨识的扬声器单元非线性参数测量方法 |
| CN101442698A (zh) * | 2008-09-26 | 2009-05-27 | 嘉兴中科声学科技有限公司 | 一种扬声器线性参数测量方法 |
-
2011
- 2011-04-02 CN CN 201110083848 patent/CN102158793B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1761366A (zh) * | 2005-08-18 | 2006-04-19 | 南京大学 | 基于系统辨识的扬声器单元非线性参数测量方法 |
| CN101442698A (zh) * | 2008-09-26 | 2009-05-27 | 嘉兴中科声学科技有限公司 | 一种扬声器线性参数测量方法 |
Non-Patent Citations (3)
| Title |
|---|
| 《华中理工大学学报》 19971231 刘春晖等 "基于Volterra模型的一类通信系统的参数辨识" , * |
| WOLFGANG.KLIPPEL: ""Fast and Accurate Measurement of the Linear Transducer Parameters"", 《第110次[国际]音频工程学会大会》 * |
| 刘春晖等: ""基于Volterra模型的一类通信系统的参数辨识"", 《华中理工大学学报》 * |
Cited By (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104170404B (zh) * | 2012-03-05 | 2018-01-26 | 奥音科技（北京）有限公司 | 音频系统、具有电‑声换能器的设备、测量装置及方法 |
| CN102968543B (zh) * | 2012-12-13 | 2015-09-02 | 嘉善恩益迪电声技术服务有限公司 | 一种扬声器音圈及磁路的温度特性数值模拟方法 |
| CN102968543A (zh) * | 2012-12-13 | 2013-03-13 | 嘉善恩益迪电声技术服务有限公司 | 一种扬声器音圈及磁路的温度特性数值模拟方法 |
| TWI508576B (zh) * | 2013-05-15 | 2015-11-11 | Lite On Opto Technology Changzhou Co Ltd | 揚聲器異音檢測方法及裝置 |
| CN104640051A (zh) * | 2013-11-06 | 2015-05-20 | 亚德诺半导体股份有限公司 | 估计扬声器的振膜偏移的方法 |
| DE102014115719B4 (de) | 2013-11-06 | 2021-08-12 | Analog Devices Global | Verfahren zur Schätzung des Membranhubs eines Lautsprechers |
| CN104640051B (zh) * | 2013-11-06 | 2018-09-11 | 亚德诺半导体股份有限公司 | 估计扬声器的振膜偏移的方法 |
| US9980068B2 (en) | 2013-11-06 | 2018-05-22 | Analog Devices Global | Method of estimating diaphragm excursion of a loudspeaker |
| CN103841505A (zh) * | 2014-02-21 | 2014-06-04 | 歌尔声学股份有限公司 | 声学产品的ccd声阻测试方法及系统 |
| CN105101032A (zh) * | 2015-08-27 | 2015-11-25 | 深圳市音沃仕科技有限公司 | 一种电声转换器复合天线模组的声学测试的方法 |
| CN106068007A (zh) * | 2016-06-07 | 2016-11-02 | 瑞声科技(新加坡)有限公司 | 扬声器非线性系统辨识方法 |
| CN106068007B (zh) * | 2016-06-07 | 2019-09-13 | 瑞声科技(新加坡)有限公司 | 扬声器非线性系统辨识方法 |
| CN106303882A (zh) * | 2016-08-13 | 2017-01-04 | 厦门傅里叶电子有限公司 | 定义扬声器工作极限的方法 |
| CN106454679A (zh) * | 2016-11-17 | 2017-02-22 | 矽力杰半导体技术(杭州)有限公司 | 扬声器振膜状态估计方法及应用其的扬声器驱动电路 |
| CN106454679B (zh) * | 2016-11-17 | 2019-05-21 | 矽力杰半导体技术(杭州)有限公司 | 扬声器振膜状态估计方法及应用其的扬声器驱动电路 |
| CN106851514A (zh) * | 2017-02-25 | 2017-06-13 | 中山市天键电声有限公司 | 扬声器非线性参数测试装置系统及测试方法 |
| CN111448804A (zh) * | 2018-01-17 | 2020-07-24 | 三星电子株式会社 | 用于扬声器驱动器的运动的非线性控制的方法和系统 |
| CN111448804B (zh) * | 2018-01-17 | 2022-04-29 | 三星电子株式会社 | 用于扬声器驱动器的运动的非线性控制的方法和系统 |
| CN108614431A (zh) * | 2018-06-08 | 2018-10-02 | 河海大学常州校区 | 一种基于夹角的Hammerstein-Wiener系统多模型分解及控制方法 |
| CN109655735A (zh) * | 2018-12-07 | 2019-04-19 | 武汉市聚芯微电子有限责任公司 | 功放芯片评估板以及功放芯片评估板系统 |
| CN109655735B (zh) * | 2018-12-07 | 2021-07-13 | 武汉市聚芯微电子有限责任公司 | 功放芯片评估板以及功放芯片评估板系统 |
| CN109738797A (zh) * | 2019-01-16 | 2019-05-10 | 珠海市精实测控技术有限公司 | 一种基于大电压模型的线性马达bl方法 |
| CN110022522A (zh) * | 2019-01-29 | 2019-07-16 | 浙江中科电声研发中心 | 采用激振器激励的扬声器振动部件共振频率测量系统及测量方法 |
| CN110022522B (zh) * | 2019-01-29 | 2023-11-07 | 浙江中科电声研发中心 | 采用激振器激励的扬声器振动部件共振频率测量系统及测量方法 |
| CN110769350A (zh) * | 2019-09-29 | 2020-02-07 | 惠州市杰德创新科技有限公司 | 喇叭响度的调节方法及喇叭 |
| CN115031863A (zh) * | 2022-08-11 | 2022-09-09 | 成都凯天电子股份有限公司 | 一种阻滞式总温传感器 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102158793B (zh) | 2013-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102158793B (zh) | 一种使用激光传感器测量扬声器参数的方法及系统 | |
| CN102118678B (zh) | 一种使用电流传感器测量扬声器参数的方法及系统 | |
| CN101411213B (zh) | 音频系统中的校准方法和装置 | |
| CN102742300B (zh) | 扬声器输出的控制 | |
| CN101442698B (zh) | 一种扬声器线性参数测量方法 | |
| CN101416533B (zh) | 在音频系统中的方法和设备 | |
| CN101600144B (zh) | 采用连续对数扫频信号获得电声产品多个参数的方法及系统 | |
| US9578416B2 (en) | Control of a loudspeaker output | |
| US6269318B1 (en) | Method for determining transducer linear operational parameters | |
| CN102843633A (zh) | 扬声器输出的控制 | |
| CN104837092A (zh) | 用于电声通信设备的回声消除方法和总成 | |
| CN101448180A (zh) | 一种移动电话扬声器测试系统 | |
| EP2456229A1 (en) | Loudspeaker system and control method | |
| CN106105264A (zh) | 用于机电电容式传感器的集成自测试 | |
| CN103096234A (zh) | 一种麦克风测试设备和方法 | |
| CN109655735B (zh) | 功放芯片评估板以及功放芯片评估板系统 | |
| CN202168211U (zh) | 一种使用激光传感器测量扬声器参数的装置 | |
| EP3476135A1 (en) | Method for simulating total harmonic distortion of a loudspeaker | |
| Garg et al. | Measurement accuracy of secondary standards of sound pressure in comparison to primary standards | |
| US20200218841A1 (en) | Prediction method of porous material and the system thereof | |
| CN111796792A (zh) | 一种手势动作判断方法、装置、电子设备和存储介质 | |
| CN100559896C (zh) | 基于电阻抗的扬声器低频频率响应的非消声室测量方法 | |
| US9136814B2 (en) | Method and system of equalization pre-preocessing for sound receivng system | |
| Yeh et al. | Nonlinear modeling of a guitar loudspeaker cabinet | |
| Kim et al. | Practical design of a speaker box with a passive vibrator (February 2018) |
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 |
Granted publication date: 20131218 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |