Page 78 - 《应该声学》2022年第2期
P. 78

246                                                                                  2022 年 3 月


             扬声器特征矩阵,其分布情况如图4所示。                                   造成扬声器异常声的原因较多,其响应时频图
                                                               特征存在差异性。为进一步减少中间状态样本数量,
                                   10
                                   ∑
                              Z =     Z a ,             (8)    对特征矩阵进行谐波带域及高次谐波域划分,应用
                                   a=1
                                                               式 (9) 提取谐波带域特征点矩阵,并将该矩阵置于
             式(8)中,Z 为扬声器特征矩阵,Z a 为剔除相关性特                      谐波带矩阵N X ,其分布情况如图 5(a)所示,其特征
             征点的合格扬声器特征矩阵。                                     点主要集中于谐波带内,文献[6]即基于扬声器谐波
                                                               带进行分析并得出检测结果,已经初步验证谐波带
                    10
                                                               内有效信息较多。
                     8
                                                                                Z x = N X · Z,            (9)
                    ᮠဋ/kHz  6 4                                式(9)中,Z x 为谐波域特征点矩阵,N X 为谐波带域

                                                               模板矩阵。
                                                                   应用式 (10) 提取其高次谐波特征点矩阵,并对
                     2
                                                               高次谐波域做进一步分割;将其分割为左右共振区
                      0    0.1  0.2   0.3  0.4   0.5           以及中间区域3部分,分布如图5(b)所示,特征点主
                                   ௑ᫎ/s
                                                               要集中于共振区域。
                         图 4  整体特征点分布情况
                                                                                Z g = N G · Z,           (10)
               Fig. 4 The distribution of feature points on the
               time-frequency spectrum                         式 (10) 中,Z g 为高次谐波特征点矩阵,N G 为高次
                                                               谐波域模板矩阵。
                   10                                          2.3  检测模型建立
                                             ࣜЯྲढ़ག
                                             ܱࣜྲढ़ག
                   8                                               将上文所提取扬声器特征矩阵 Z 作为特征提
                                                               取模板,对 3 位专业声音检听员均认证合格的 50 支
                   6
                  ᮠဋ/kHz  4                                    扬声器进行特征提取,将所提取赋值特征点按其在
                                                               矩阵中每列纵坐标从小到大的顺序进行排列,连接
                                                               排列后特征点构建扬声器特征曲线;将50条合格扬
                   2
                                                               声器特征曲线叠加,提取上下边界值得到如图 6 所
                                                               示合格扬声器特征曲线边界,并将其作为检测模型
                    0     0.1   0.2   0.3   0.4   0.5
                                                               的检测域构建全特征点检测模块。若被测扬声器
                                   ௑ᫎ/s
                                                               特征曲线超出检测域点数低于所设定最小阈值,如
                              (a) ៈฉࣜ۫ྲढ़Ѭ࣋
                                                               图 7(a) 所示,则判定该扬声器为合格产品;若超出
                   10
                                                               检测域点数高于所设定最大阈值,如图 7(b) 所示,
                   8
                                                                      15
                                                                                                ೝ฾ʽᬍ
                   6
                  ᮠဋ/kHz  4                                           10
                                                                                                ೝ฾ʾᬍ
                                                                     ᑟ᧚ϙ
                                   ࢻС૝۫
                   2
                                   ˗ᫎӝ۫                               5
                                   ԿС૝۫
                     0    0.1   0.2   0.3   0.4   0.5
                                   ௑ᫎ/s
                                                                        0      50     100     150     200
                            (b) ᰴ൓ៈฉ۫Ѭҟྲढ़Ѭ࣋                                         ྲढ़གऀՂ

                            图 5  特征分布情况                                         图 6  全检测域
                   Fig. 5 The distribution of feature points             Fig. 6 Complete detection area
   73   74   75   76   77   78   79   80   81   82   83