一种扬声器异常音检测方法[发明专利]

(12)发明专利申请

(10)申请公布号(10)申请公布号 CN 1043635 A(43)申请公布日 2015.02.18

(21)申请号 201410511623.5(22)申请日 2014.09.29

(71)申请人嘉善恩益迪电声技术服务有限公司

地址314100 浙江省嘉兴市嘉善县晋阳东路

568号科创中心1号楼1层申请人浙江中科电声研发中心(72)发明人徐楚林 温周斌 李宏斌(74)专利代理机构杭州君度专利代理事务所

(特殊普通合伙) 33240

代理人沈志良(51)Int.Cl.

H04R 29/00(2006.01)

权利要求书2页 说明书7页 附图1页权利要求书2页 说明书7页 附图1页

()发明名称

一种扬声器异常音检测方法(57)摘要

本发明提出了一种扬声器异常音检测方法。该方法只要激励信号等条件不改变，每次检测只需重复进行第三阶段即可，因此本算法的计算量非常低，可以实现实时处理。为了进一步提高估计精度可以定期执行第四阶段的更新操作，利用检测空闲时间进行更新操作或使用异步运算方式可以在保持实时处理的同时更新参数，从而提高本算法性能。本发明所述方法的主要优势是：一次性准确得到所有阶次谐波；抗噪声性能好，虚警率低，计算量小；可实时处理，实时发现问题；对激励信号形式和频率范围没有算法上的。由于有这些特点，采用本方法可以开发出性能良好的扬声器异常音自动检测设备，以快速准确检测发现扬声器的异常音故障，还可同时检测扬声器频率响应、阻抗曲线等常规参数。

CN 1043635 A CN 1043635 A

权 利 要 求 书

1/2页

1.一种扬声器异常音检测方法，其特征在于所述的该方法包含以下阶段：

第一阶段：根据被测扬声器特点确定激励信号，进一步得到每个采样点处激励信号的瞬时频率、瞬时相位

、当前最高谐波次数

和整个过程的最高谐波次数

；根据需要

选择系统阶次M、平滑距离S、组成异常音曲线的高阶谐波阶次范围；确定量测噪声方差和状态噪声协方差矩阵的估计值；

第二阶段：首先准备状态向量协方差矩阵的初始值，对应激励信号的每一个采样点顺序按算法迭代构造为：

初始化状态协方差矩阵P：对整数的0值方阵；

准备状态修正向量;

对从1到N的每个n顺序迭代进行：测量向量

保存

和整数

，赋值为边长

和更新

，并在此基础上算出并保存状态修正向量

，具体操作步骤

对、：（其中

是组合数）

更新并保存状态修正向量：对保存

，

2

CN 1043635 A

权 利 要 求 书

2/2页

更新状态协方差矩阵：

；

对对

，和

，

；

第三阶段：对每个待测扬声器，将待测扬声器和麦克风装在消声箱内指定位置，在待测扬声器两输入端加载激励信号，用声卡或数据采集卡接收麦克风采集信号，确保接收信号与激励信号时间对齐；从初始值开始，对应每个接收信号采样点顺序迭代更新状态估计向量

，并根据估计结果

得到异常音曲线;

将异常音曲线与根据经验确定的阈值比较得出是否存在异常音故障的判断；第四阶段：基于历史数据用统计方法估计、、和新进行第二阶段的计算，得到更符合测试条件的

更准确的值，利用这些参数重

序列，从而在之后的检测中进一步提高

检测性能。

2.权利要求1所述的扬声器异常音检测方法，其特征在被测扬声器包括各种动圈式电动扬声器。

3.权利要求1所述的扬声器异常音检测方法，其特征在于在主要参数不改变的情况下对任一扬声器的检测只需进行第三阶段操作。

3

CN 1043635 A

说 明 书

一种扬声器异常音检测方法

1/7页

技术领域

[0001]

本发明属于信号处理领域，涉及扬声器异常音故障的一种检测方法。

背景技术

扬声器是一种重要的电声器件。由于制造工艺问题，部分扬声器产品会出现音圈摩擦磁路部件、音圈碰撞磁路底部、振膜或纸盆由于粘接不牢固出现漏气、内部混入颗粒杂质，等等诸多问题，这些问题统称为异常音故障。相对于扬声器的线性和正常的非线性输出而言，扬声器异常音是十分微弱的，异常音的自动检测方法是一个世界性难题，目前对于异常音问题电声企业主要采用人工检听的方法。

[0003] 经过多年研究也出现了一些异常音自动检测方法，例如基于对数扫频信号和Volterra模型的算法、基于短时傅里叶变换的算法、基于小波的方法、基于经验模态分解（EMF）的算法等。但是这些方法也有一些不足，例如：多数方法需要得到完整的响应信号才能开始计算，不能实时处理，不利于尽快发现问题；有些方法的计算量比较大；某些方法只适用于特定种类的激励信号；某些方法对激励和采样的频率范围有较高要求等。

[0002]

发明内容

[0004] 本发明的目的是针对上述问题，设计出一种扬声器异常音故障的一种检测方法。它采用信号处理方法根据激励信号和麦克风采集的声信号自动判定该扬声器是否存在异常音故障。本发明基于自适应Kalman滤波方法，实现准确、抗噪声、高效地检测异常音故障。

[0005] 本发明中提及的一种扬声器异常音检测方法，具体步骤是：

第一阶段：参数设定

1)根据被测扬声器特点选定激励信号频率范围，这个范围应包括该扬声器容易产生异常音故障的主要频率点，选定任一种形式的连续扫频信号作为激励信号，写为

的形式，其中

是激励信号的幅值序列；是激励信号的相位序

一般通过加窗使初始列；激励信号的总采样点数为N。为了避免幅值跳变产生爆音的现象，幅值

并逐步增大到异常音测试需要的水平；

2)根据激励信号形式得到在第个采样点处激励信号的瞬时频率输出信号的采样频率得到每一采样时刻最高谐波次数

记为

；

，又根据声卡

，和整个过程的最高谐波次数

4

CN 1043635 A

说 明 书

2/7页

3)根据需要选择系统阶次M和平滑距离一般选择M=1或2，。越大对谐波幅

值估计越准确，同时计算量和所需存储空间增大；

4)根据被测扬声器类型和特点选择构成异常音曲线的高阶谐波阶次范围；一般可取所有高于7阶的谐波作为异常音曲线；

5)根据环境噪声和声信号本底噪声估计量测噪声方差行

列的状态噪声协方差矩阵第二阶段：准备状态修正向量6)初始化状态协方差矩阵P：对方阵。

[0006]

；根据扬声器特性估计

；

和，赋值为边长的0值

7)准备状态修正向量。对从1到N的每个n顺序迭代进行：

a)测量向量保

存

b)c)对

、（其中

：

是组合数）

d)更新并保存：对保存e)更新：

；

5

，

CN 1043635 A

说 明 书

，和

，

；

3/7页

对对

第三阶段：实时检测

8) 将待测扬声器和麦克风装在消声箱内指定位置，在待测扬声器两端加载激励信号，用声卡或数据采集卡接收麦克风采集信号，并采用一般的信号处理方法找到与激励信号时间对齐的部分，记为响应；

9)初始化状态向量） 其中别是，当激励的

是幅值为

： 对频率为

，

（向量总长度为

，即确保对每个整数

，是扬声器对

的

的激励信号作用于该型扬声器测得的声信号幅值。特

行1列的0向量；

时，初始化为

10)对从1到N的每个n顺序迭代进行：a)

b)更新状态向量：

（其中

对c)当

，时，向量

是根据接收信号采样点

估计的

采样时刻，幅值为

是组合数）；

的系统状态。 其中，第阶谐波为

；根据预先确定的高阶谐波阶次范围，将范围内的谐波幅值求

平方和得到异常音曲线的第个样点；

11)如果待检测扬声器的异常音曲线超过预先确定的该型扬声器异常音曲线阈值则说明待检测扬声器存在异常音故障，反之说明待检测扬声器合格； 异常音曲线阈值可按上述方法预先处理多个同型号无异常音的扬声器，根据这些扬声器的异常音曲线包络的最大值加一定余量构成；

6

CN 1043635 A

说 明 书

4/7页

第四阶段：更新统计量

12)基于上述处理结果可以利用统计方法估计、、和

等参数更准确的初始值。

例如，统计无信号时的噪声方差作为的初始值、系统方程的误差方差作为的初始值、作为的

的初始值等。定期利用这些参数重新进行第二阶段的计算，可得到更符合测试条件序列，从而在之后的检测中进一步提高检测性能。

13)第二阶段所述更新

、计算状态修正向量

和第三阶段所述新状态估计向

[0007]

量的方法是在两种公开技术：“自适应Vold-Kalman阶次跟踪方法”和“固定延迟平滑

Kalman滤波方法”的基础上改进结合形成的。[0008] 从上述步骤中可以知道，只要激励信号等条件不改变，每次检测只需重复进行第三阶段即可，因此本算法的计算量非常低，可以实现实时处理。为了进一步提高估计精度可以定期执行第四阶段的更新操作，利用检测空闲时间进行更新操作或使用异步运算方式可以在保持实时处理的同时更新参数，从而提高本算法性能。[0009] 本发明所述方法的主要优势是：1、一次性准确得到所有阶次谐波；2、抗噪声性能好，虚警率低，3、计算量小；4、可实时处理，实时发现问题；5、对激励信号形式和频率范围没有算法上的。由于有这些特点，采用本方法可以开发出性能良好的扬声器异常音自动检测设备，以快速准确检测发现扬声器的异常音故障，还可同时检测扬声器频率响应、阻抗曲线等常规参数。

附图说明

[0010] 图1为本发明检测流程框图。

具体实施方式

[0011] 下面以某款微型扬声器异常音故障检测过程为例说明本发明的具体实施方式：

1)扬声器异常音故障往往体现在共振频率附近。该微型扬声器共振频率约700Hz，经验表明，选择300Hz到2000Hz频率的激励信号可以激励出有故障的扬声器的异常音现象。本例使用对数扫频信号，激励时间0.5秒，采用参数为0.03的Tukey窗，声卡采样频率为48000Hz，因此激励信号长度

其中

，激励信号可写为：

7

CN 1043635 A

说 明 书

5/7页

2)根据对数扫频信号性质，上述激励信号的频率序列可写作：

整运算），最大谐波次数

[0012] [0013] [0014]

每个采样点处谐波次数。

），平滑距离

。

（为下取

3)本例使用2阶系统模型（

4)取高于7阶的谐波组成异常音曲线。5)估计

，

，其中，所有

和

是边长为160的单位矩阵。都赋值为边长160的0值方阵。

[0015] 6)初始化P： 对7)准备

[0016] 。 对从1到N=24000的每个n顺序迭代进行：

a)保存

b)c)对

：

d)更新并保存：对保存e)更新

：

8

，

CN 1043635 A

说 明 书

；；

6/7页

对，

以上是进行具体检测前的准备工作。针对某一个扬声器的检测在以上工作保存的内容的基础上进行以下操作步骤：

8)将待测扬声器和麦克风装在消声箱内指定位置，在待测扬声器两端加载激励信号，用声卡或数据采集卡接收麦克风采集信号，并采用一般的信号处理方法找到与激励信号时间对齐的部分，记为（即确保对每个整数

[0017]

，是扬声器对的响应）。

，全部初

9)初始化： 由于加窗使激励信号幅值从0开始，因此，对

行1列的0向量。

10)对从1到N的每个n顺序迭代进行：

始化为

[0018]

a)

b)更新：

；

对c)当

，时，向量

是根据接收信号采样点

估计的

采样时，幅值为

刻的系统状态。其中，第阶谐波为

。因为是采用所有高于7阶的谐波组成异常音曲线，因此异常

音曲线对应序列第

[0019]

个样点处的值为。

11)异常音曲线阈值可按上述方法预先处理多个同型号无异常音的扬声器，根据这些扬声器的异常音曲线包络的最大值加一定余量构成。 如果待检测扬声器的异常音曲线超过预先确定的该型扬声器异常音曲线阈值则说明待检测扬声器存在异常音故障，反之说明待检测扬声器合格。从以上步骤可以看到本方法在发送激励信号的同时得到声压采样序列，随着声压采样序列的取得可以同步地（仅有S-1个样点的延迟）计算得到异常音曲线

9

CN 1043635 A

说 明 书

7/7页

并与阈值进行比较。因此本方法可以实现实时判定异常音故障。[0020] 在经过一段时间检测后，在已有数据的基础上可以进行以下步骤的更新操作。根据所使用电脑的性能调整进行更新操作的频率，实现在不影响实时检测的同时提高检测性能。

[0021]

12)基于上述处理结果可以用统计方法改善对、、和等的估计值：

a)每个样点处的值（步骤10.a）构成余量序列，通过统计方法估计这个序列的协方差作为的值。

[0022]

b)在的情况下，对每一个满足，向量是状态噪

声，通过统计方法估计这个向量的协方差矩阵作为的值。

[0023]

c)对本例而言，由于初始状态输入谐波幅值为0，输出的各次谐波必然为0，因此是准确的，不需要调整估计。 对于初始时刻输入幅值不为0的情况，可以统计各次的值，这些向量的均值是对

更好的估计，协方差矩阵是对

更好的估计。

和检测中

[0024]

在得到这些初始参数更好的估计值后，可以在空闲时段或在后台线程中进行第二

序列。在更新的

序列基础上进行的第三阶段检测工作将更

阶段操作，从而得到更新的

适应实际检测环境和当前扬声器产品的特点，因此会有更好的检测性能。

[0025] 以上实施实例及实施步骤仅用以说明本发明而并非本发明所描述的技术方案。尽管本说明书参照上述实施实例对本发明进行了详细的说明，但本领域的普通技术人员应该理解，在不脱离所附权利要求的所限定的精神和范围的情况下，可做出许多修改、变化或等效替换；而一切不脱离本发明的精神和范围的技术方案及其改进，其均应涵盖于本发明的权利保护范围中。

10

CN 1043635 A

说 明 书 附 图

1/1页

图1

11