设备状态分类模型的构建方法、使用方法、计算机可读介质、装置与流程

1.本发明涉及设备状态识别技术领域，具体涉及一种设备状态分类模型的构建方法、一种设备状态分类模型的使用方法、一种存储有设备状态分类程序的计算机可读介质、一种设备状态辨识装置，其可以应用在变压器状态识别上。


背景技术：

2.变压器运行状态改变时，会产生不同于正常工况的异常振动。因此，通过采集、分析变压器运行时产生的振动信息有助于了解变压器的运行状态。由于变压器体积庞大，且带电部位(如套管)不易布置传感器，因此，即使大量布设振动传感器，也只能测量局部的振动信号。由于异常振动会向外辐射噪声，因此，通过采集、分析声纹也有助于也解变压器的运行状态。但是声纹通过空气传播到传声器阵列，存在一定程度的衰减，且易受环境噪音的干扰，与振动信号相比信息完整性较低。
3.专利文献cn114167315a公开了一种变压器智能在线监测系统，包括声纹传感器、振动传感器和处理器，声纹传感器用于采集变压器在工作时的声纹信号，振动传感器用于采集变压器在工作时的振动信号，处理器获取所述声纹传感器采集的声纹信号和振动传感器采集的振动信号，提取对应的频谱特征；并将其获取得到的采集变压器正常工作时的频谱特征作为初始特征；所述处理器构建神经网络模型，并以获取得到的变压器在正常工作时的声纹信号和振动信号及频谱特征作为样本进行训练得到预测模型，并根据所述预测模型对获取得到的变压器在工作时的声纹信号和振动信号进行预测，得到其对应的频谱特征，并与所述初始特征进行比对，并据此判断是否存在异常。其通过神经网络模型对变压器运行时产生的声纹信号、振动信号进行分类，由于声音是由振动引起的，也即声音与振动具有相关性，因此，此种方法准确率低。


技术实现要素：

4.本发明的目的是提供一种设备状态分类模型的构建方法、一种设备状态分类模型的使用方法、一种存储有设备状态分类程序的计算机可读介质、一种设备状态辨识装置，以从新的技术路线的角度实现设备状态分类。
5.本发明的技术方案是：
6.一种设备状态分类模型的构建方法，包括以下步骤：
7.步骤1、获取设备不同运行状态下的声纹数据和振动数据；
8.步骤2、提取声纹数据的短时过零率、谱质心和短时能量，提取振动数据的偏度、峭度、波峰因子和形状因子；
9.步骤3、使用融合分类算法处理短时过零率、谱质心、短时能量、偏度、峭度、波峰因子和形状因子，获得设备状态分类模型。
10.优选的，在所述步骤2之前，还使用稀疏变换算法处理所述声纹数据和/或所述振
动数据。
11.进一步优选的，所述稀疏变换算法包括：
12.取所述步骤1获得的设备不同运行状态下的声纹数据或振动数据，得到样本数据；
13.使用测量矩阵得到压缩测量值y＝af，式中，y为数据的压缩测量值，a为数据的测量矩阵，f为样本数据；其中，测量矩阵采用托普利兹矩阵；
[0014][0015]
式中，b表示稀疏系数矩阵，f表示样本数据，αi表示稀疏系数矩阵b的列向量，表示某一常数，n
*
表示正整数；
[0016]
对于稀疏系数矩阵b的元素α
(n)
，采用最小l1范数优化求解，即有
[0017][0018]
对于稀疏基向量φ，采用moore-penrose广义逆求解方法进行求解，即有
[0019]
φ
(n 1)
＝x(α
(n)
)
t
[α
(n)
(α
(n)
)
t
]-1
ꢀꢀꢀ
(3)式中，“t”表示矩阵转置；
[0020]
声纹数据或振动数据的稀疏系数矩阵为式中，y为数据的压缩测量值，a为数据的测量矩阵，φ为优化后的稀疏基向量。
[0021]
优选的，所述步骤2中获得的短时过零率为
[0022]
所述步骤2中获得的谱质心为
[0023]
所述步骤2中获得的短时能量为
[0024]
所述步骤2中获得的偏度为
[0025]
所述步骤2中获得的峭度为
[0026]
所述步骤2中获得的波峰因子为
[0027]
所述步骤2中获得的形状因子为其中，ξ表示某一常数，n＝1,2,
…
,n，n为信号时域总离散点，xn表示第n个离散点对应的幅值，fn表示信号频率，en表示离散时域信号短时傅里叶变换后所对应频率的谱能量，m表示频域离散总点数，p(
·
)表示概率计算，ωn(
·
)表示汉明窗函数，υ表示帧移长度，表示原始信号幅值的平均值。
[0028]
优选的，在所述步骤3中，所述融合分类算法为梯度提升决策树方法，其映射关系可定义为
[0029][0030]
式中，β表示初始模型，yi表示预测值，loss(yi,f(x))表示损失函数，f(x)表示与模型相关的函数值。
[0031]
优选的，所述设备为变压器。
[0032]
一种设备状态分类模型的使用方法，包括以下步骤：
[0033]
步骤a、采集所述设备的声纹信号和振动信号；
[0034]
步骤b、提取声纹数据的短时过零率、谱质心和短时能量，提取振动数据的偏度、峭度、波峰因子和形状因子；
[0035]
步骤c、将所述短时过零率、所述谱质心、所述短时能量、所述偏度、所述峭度、所述波峰因子和所述形状因子输入所述设备状态分类模型，所述设备状态分类模型输出设备状态。
[0036]
一种存储有设备状态分类程序的计算机可读介质，所述设备状态分类程序被处理器时执行以下步骤：
[0037]
步骤a、获取所述设备的声纹信号和振动信号；
[0038]
步骤b、提取声纹数据的短时过零率、谱质心和短时能量，提取振动数据的偏度、峭度、波峰因子和形状因子；
[0039]
步骤c、将所述短时过零率、所述谱质心、所述短时能量、所述偏度、所述峭度、所述波峰因子和所述形状因子输入所述设备状态分类模型，所述设备状态分类模型输出设备状态。
[0040]
一种设备状态辨识装置，包括处理器和前述的存储有设备状态分类程序的计算机可读介质，所述处理器用于执行前述的设备状态分类程序。
[0041]
优选的，还包括声纹传感器、振动传感器和输出模块，所述声纹传感器用于采集所述设备的声纹信号，所述振动传感器用于采集所述设备的振动信号，所述处理器用于接收所述声纹信号和振动信号，执行所述设备状态分类程序，获得设备状态信息，所述输出模块用于输出所述设备状态信息。
[0042]
本发明的有益效果是：
[0043]
1.本发明的设备状态分类模型的构建方法中，步骤1获取设备不同运行状态下的声纹数据和振动数据，与仅获取设备不同运行状态下的声纹数据或振动数据相比，其携带的信息量更大、更全面，更能准确反映设备工作状态；步骤2中，通过使用短时过零率、谱质心和短时能量表征声纹特征，使用偏度、峭度、波峰因子和形状因子表征振动特征，使其更准确、全面的反映声纹信息和振动信息；步骤3中，使用融合分类算法处理短时过零率、谱质心、短时能量、偏度、峭度、波峰因子和形状因子，获得设备状态分类模型，这样，通过丰富设备状态的信息参量，获得的设备状态分类模型分辨设备工作状态的准确率更高。由于声纹信号具备整个声场的全局特性，而振动信号则对于指定测试点的局部特征更为敏感，两者不同源信号进行稀疏表征与信息融合后的优势在于：相较于之前单一的信号源融合表征方法，所述方法既能够得到测试设备丰富的检测信息参量，同时也显著提升检测状态的辨识度。
[0044]
2.在所述步骤2之前，使用稀疏变换算法处理所述声纹数据和/或所述振动数据。
这样，可以在不丢失原始信号测试信息的情况下，大幅减少了数据通信代价，降低硬件规格要求，提高获得设备状态分类模型的效率。
[0045]
3.通过重新设计稀疏基的优化目标函数，利用l1范数代替原先的l0范数进行优化求解，其优势在于：显著提高了计算效率，避免了l0范数很难优化求解(np难问题)；显著提高了求解精度，更易收敛。
[0046]
4.通过引入短时过零率、谱质心和短时能量、偏度、峭度、波峰因子和形状因子七个特征参量更好地表征声纹、振动稀疏信号，更为全面反映原始信息。
[0047]
5.通过引入梯度提升决策树算法进行信息融合与分类，相较于常规融合分类方法有如下优势：并行计算、显著提高了计算效率；对于采集得到的声纹、振动海量数据集，提升了泛化、表达及数据融合能力；利用所述方法得到的融合分类模型的可解释度与鲁棒性较高。
[0048]
6.本发明的设备状态分类模型的使用方法，可以分辨设备状态。
[0049]
7.本发明的存储有设备状态分类程序的计算机可读介质，设备状态分类程序被处理器时执行时，可以分辨设备状态。
[0050]
8.本发明的设备状态辨识装置，可以分辨设备状态。
附图说明
[0051]
图1为一种设备状态分类模型的构建方法的流程图。
[0052]
图2为一种变压器状态分类模型的构建方法中构造最优稀疏基过程中的损换函数曲线。
[0053]
图3为一种变压器状态分类模型的构建方法中构造最优稀疏基过程中的准确率曲线。
[0054]
图4为一种变压器状态分类模型的构建方法中变压器正常状态反演得到的稀疏系数。
[0055]
图5为一种变压器状态分类模型的构建方法中变压器相间放电状态反演得到的稀疏系数。
[0056]
图6为一种变压器状态分类模型的构建方法中变压器对铁芯放电状态反演得到的稀疏系数。
[0057]
图7为实施例3所述的构建方法构建的变压器状态分类模型的混淆矩阵计算结果。
[0058]
图8为一种设备状态辨识装置的信号框图。
具体实施方式
[0059]
下面结合附图，以实施例的形式说明本发明，以辅助本技术领域的技术人员理解和实现本发明。除另有说明外，不应脱离本技术领域的技术知识背景理解以下的实施例及其中的技术术语。
[0060]
实施例1：一种设备状态分类模型的构建方法，包括以下步骤：
[0061]
步骤1、获取设备不同运行状态下的声纹数据和振动数据；应当明白，信号与数据指同一存在的不同表现形式，数据为信号记录在存储器中的状态，信号为数据传送时的状态。因此，也可以直接从采集装置中获取设备不同运行状态下的声纹信号和振动信号，这也
等同于获取设备不同运行状态下的声纹数据和振动数据。
[0062]
步骤2、提取声纹数据的短时过零率、谱质心和短时能量，提取振动数据的偏度、峭度、波峰因子和形状因子；
[0063]
步骤3、使用融合分类算法处理短时过零率、谱质心、短时能量、偏度、峭度、波峰因子和形状因子，获得设备状态分类模型。
[0064]
本实施例中，融合分类算法采用梯度提升决策树方法。在其它实施例中，融合分类算法还可以使用加权平均法、基于stacking的交叉验证融合分类算法等，但本实施例中所采用的梯度提升决策树方法的优势在于：1、预测阶段的计算速度快，树与树之间可并行化计算；2、在分布稠密的数据集上，泛化、表达及数据融合能力强；3、得到的模型具有较好的解释性和鲁棒性，能够完整诠释特征间的高阶深层关系。
[0065]
实施例2：一种设备状态分类模型的构建方法，包括以下步骤：
[0066]
步骤1、获取设备不同运行状态下的声纹数据和振动数据；
[0067]
步骤2、使用稀疏变换算法处理所述声纹数据和所述振动数据；本步骤可以减少步骤3过程的数据处理量，降低对计算机硬件的要求，提高步骤3的效率。在其它实施例中，也可以仅使用稀疏变换算法处理所述声纹数据，或者仅使用稀疏变换算法处理振动数据。
[0068]
本实施例中，使用稀疏变换算法处理所述声纹数据包括以下分步骤：
[0069]
分步骤2a、取所述步骤1获得的设备不同运行状态下的声纹数据，得到声纹样本数据；
[0070]
分步骤2b、使用测量矩阵得到声纹数据的压缩测量值y
p
＝a
pfp
，式中，y
p
为声纹数据的压缩测量值，a
p
为声纹数据的测量矩阵，f
p
为声纹样本数据；其中，声纹数据的测量矩阵采用托普利兹矩阵；
[0071]
托普利兹矩阵形式为
[0072][0073]
式中，n、m表示元素序号，该矩阵的特点在于构成该测量矩阵的元素完全由其第1行和第1列的2n-1个元素确定；除第一行第一列外，其他每个元素都与左上角的元素相同。由于构造托普利兹测量矩阵需要的自由元素较少、结构简单、循环移位易于在硬件上实现，故该矩阵具有良好的工程应用前景；
[0074]
分步骤2c、采用最佳方向法构造、优化稀疏基向量φ，优化目标函数为
[0075][0076]
式中，b表示稀疏系数矩阵，f表示样本数据，αi表示稀疏系数矩阵b的列向量，表示某一常数，n
*
表示正整数；
[0077]
对于稀疏系数矩阵b的元素α
(n)
，采用最小l1范数优化求解，即有
[0078][0079]
对于稀疏基向量φ，采用moore-penrose广义逆求解方法进行求解，即有
[0080]
φ
(n 1)
＝x(α
(n)
)
t
[α
(n)
(α
(n)
)
t
]-1
ꢀꢀ
(8)
[0081]
式中，“t”表示矩阵转置；
[0082]
分步骤2d、声纹数据的稀疏系数矩阵为式中，y
p
为数据的压缩测量值，a
p
为数据的测量矩阵，φ
p
为优化后的稀疏基向量，也即分步骤2c优化后获得的稀疏基向量φ。
[0083]
本实施例中，使用稀疏变换算法处理所述振动数据包括以下分步骤：
[0084]
分步骤2a、取所述步骤1获得的设备不同运行状态下的振动数据，得到振动样本数据；
[0085]
分步骤2b、使用测量矩阵得到振动数据的压缩测量值yv＝a
vfv
，式中，yv为振动数据的压缩测量值，av为振动数据的测量矩阵，fv为振动样本数据；其中，振动数据的测量矩阵采用托普利兹矩阵；
[0086]
托普利兹矩阵形式为
[0087][0088]
式中，n、m表示元素序号，该矩阵的特点在于构成该测量矩阵的元素完全由其第1行和第1列的2n-1个元素确定；除第一行第一列外，其他每个元素都与左上角的元素相同。由于构造托普利兹测量矩阵需要的自由元素较少、结构简单、循环移位易于在硬件上实现，故该矩阵具有良好的工程应用前景；
[0089]
分步骤2c、采用最佳方向法构造、优化稀疏基向量φ，优化目标函数为
[0090][0091]
式中，b表示稀疏系数矩阵，f表示样本数据，αi表示稀疏系数矩阵b的列向量，表示某一常数，n
*
表示正整数；
[0092]
对于稀疏系数矩阵b的元素α
(n)
，采用最小l1范数优化求解，即有
[0093][0094]
对于稀疏基向量φ，采用moore-penrose广义逆求解方法进行求解，即有
[0095]
φ
(n 1)
＝x(α
(n)
)
t
[α
(n)
(α
(n)
)
t
]-1
ꢀꢀꢀ
(12)式中，“t”表示矩阵转置；
[0096]
分步骤2d、振动数据的稀疏系数矩阵为式中，yv为振动数据的压缩测量值，av为数据的测量矩阵，φv为优化后的稀疏基向量，也即分步骤2c优化后获得的稀疏基向量φ。
[0097]
步骤3、提取处理后的声纹数据的短时过零率、谱质心和短时能量，提取处理后的振动数据的偏度、峭度、波峰因子和形状因子；具体的，提取稀疏系数向量为b
p
的短时过零率、谱质心、短时能量：
[0098]
短时过零率为
[0099]
谱质心为
[0100]
短时能量为
[0101]
提取稀疏系数向量为bv的偏度、峭度、波峰因子和形状因子：
[0102]
偏度为
[0103]
峭度为
[0104]
波峰因子为
[0105]
形状因子为其中，ξ表示某一常数，n＝1,2,
…
,n，n为信号时域总离散点，xn表示第n个离散点对应的幅值，fn表示信号频率，en表示离散时域信号短时傅里叶变换后所对应频率的谱能量，m表示频域离散总点数，p(
·
)表示概率计算，ωn(
·
)表示汉明窗函数，υ表示帧移长度，表示原始信号幅值的平均值。
[0106]
步骤4、使用融合分类算法处理短时过零率、谱质心、短时能量、偏度、峭度、波峰因子和形状因子，获得设备状态分类模型。
[0107]
本实施例中，融合分类算法为梯度提升决策树方法，其映射关系可定义为
[0108][0109]
式中，β表示初始模型，yi表示预测值，loss(yi,f(x))表示损失函数，f(x)表示与模型相关的函数值；
[0110]
这样，获得的设备状态分类模型属于一种随机森林分类器模型。
[0111]
本实施例的方法首先采集变压器设备不同状态的声纹与振动原始数据，其中一部分数据用于获取声纹、振动信号的稀疏变换基，剩余部分数据作为压缩数据y
p
、yv；接着设定残差阈值，随机选取声纹、振动信号数据中的k个元素作为初始稀疏变换基，通过优化迭代至残差满足设定条件，输出最优稀疏变换基φ
p
和φv；然后基于最优稀疏变换基φ
p
和φv对压缩数据y
p
和yv进行稀疏变换，以得到先验信息，并借助l1范数优化求解算法得到两种信号的稀疏系数向量；然后分别提取压缩稀疏声信号的短时过零率、谱质心、短时能量作为特征参量，提取压缩稀疏振动信号的偏度、峭度、波峰因子、形状因子作为特征参量，进而进行特征融合；最后利用随机森林分类器实现变压器不同工况的状态辨识。上述方法在不丢失原始信号测试信息的情况下，大幅减少了数据通信代价，且由于融合了声、振信号的各自特点，因此能够得到更为丰富的测试设备特征信息参量。
[0112]
实施例1和实施例2的方法中，设备可以是变压器，但也可以是其它工作时会发出
振动的设备。
[0113]
实施例3：一种变压器状态分类模型的构建方法，参见图1，包括以下步骤：
[0114]
1)首先采集变压器设备不同状态的声纹与振动原始数据，其中一部分数据用于获取声纹、振动信号的稀疏变换基。
[0115]
测量矩阵采取托普利兹矩阵a，其形式为：
[0116][0117]
式中，n、m表示元素序号，该矩阵的特点在于构成该测量矩阵的元素完全由其第1行和第1列的2n-1个元素确定；除第一行第一列外，其他每个元素都与左上角的元素相同。由于构造托普利兹测量矩阵需要的自由元素较少、结构简单、循环移位易于在硬件上实现，故该矩阵具有良好的工程应用前景。
[0118]
则利用测量矩阵即可得到相应的压缩测量值：
[0119][0120]
式中，f
p
、fv分别表示原始声纹、振动测试信号，a
p
、av分别表示原始声纹、振动测试信号的测量矩阵，y
p
、yv分别表示原始声纹、振动测试信号的压缩测量值。
[0121]
2)进一步地，采用最佳方向法构造、优化稀疏基向量φ，其优化目标函数如下：
[0122][0123]
式中，b表示稀疏系数矩阵，f表示原始测试信号，αi表示稀疏系数矩阵b的列向量，表示某一常数，n
*
表示正整数。
[0124]
显然，该方法的计算过程是稀疏系数求解与稀疏基更新交替进行的过程。对于稀疏系数矩阵b的元素α
(n)
，可采用最小l1范数优化求解，即有：
[0125][0126]
对于稀疏基向量φ计算，可采用moore-penrose广义逆求解方法进行求解：
[0127]
φ
(n 1)
＝x(α
(n)
)
t
[α
(n)
(α
(n)
)
t
]-1
ꢀꢀ
(18)
[0128]
上式中，“t”表示矩阵转置，值得说明的是，声纹与振动信号的稀疏变换基皆由式(16)～式(18)计算可得，限于篇幅，给出通用的计算过程。
[0129]
3)进一步地，利用步骤1得到的压缩值和步骤2得到的稀疏变换基计算声纹和振动信号的稀疏系数向量：
[0130][0131]
式中，φ
p
、φv分别表示声纹和振动信号的稀疏变换基，b
p
、bv分别表示声纹和振动信号的稀疏系数矩阵。
[0132]
4)进一步地，分别提取稀疏声纹信号b
p
的短时过零率、谱质心、短时能量作为特征
参量，提取稀疏振动信号bv的偏度、峭度、波峰因子、形状因子作为特征值作为判别当前变压器状态的依据，并采用特征表达式见下表。
[0133]
表1两种稀疏信号的特征值
[0134][0135]
上表中，ξ表示某一常数，n＝1,2,
…
,n，n为信号时域总离散点，xn表示第n个离散点对应的幅值。fn表示信号频率，en表示离散时域信号短时傅里叶变换后所对应频率的谱能量，m表示频域离散总点数，p(
·
)表示概率计算。ωn(
·
)表示窗函数，本发明方法采用汉明窗，υ表示帧移长度。表示原始信号幅值的平均值。
[0136]
5)进一步地，采用梯度提升决策树方法实现声振特征值的融合与设备的状态分类，其映射关系可定义为：
[0137][0138]
式中，β表示初始模型，yi表示预测值，loss(yi,f(x))表示损失函数，f(x)表示与模型相关的函数值。
[0139]
针对某一实测数据采用本实施例的方法进行计算，算例截取原始声纹、振动数据离散点数为500、根据压缩感知理论设定压缩率为60％，则有观测矩阵的维度为200
×
500，稀疏基维度为500
×
200，稀疏系数的维度为200
×
1。
[0140]
其计算、构造最优稀疏基过程中的性能曲线如图2、3所示。图2中曲线表示对应的损失函数值。其中图3中曲线表示准确率的变化。结果表明，损失函数在经过60次迭代后就能快速达到稳定值。对应的分类准确率也在迭代60次后达到稳定，其准确率几乎达到100％。
[0141]
如图4-6所示为变压器三种状态反演得到的稀疏系数示意，可见经过稀疏技术处理后的三种状态数据在减少采样点数量的情况下波形仍存在一定的差异性。
[0142]
为评判本发明方法对于变压器不同状态下的辨识性能，引入混淆矩阵进行直观地评判。如图7所示为混淆矩阵计算结果。
[0143]
引入混淆矩阵进行分类效果的定量评价，如图7与表2所示即为利用所提方法得到的混淆矩阵计算结果。图中混淆矩阵的对角线分别表示利用所提方法对于每种状态的识别准确率。结合图7和表2可知，在全部300个测试集中只有4个被错误分类，测试样本集的分类正确率为98.7％。由表2可知，所提方法的状态识别精度、召回率和f1值分别为1、0.98和0.98。这一结果表明所提方法具备良好的分类性能，证明了所提方法的有效性。
[0144]
表2评价参数表
[0145][0146][0147]
实施例4：一种设备状态分类模型的使用方法，包括以下步骤：
[0148]
步骤a、采集设备的声纹信号和振动信号；一般地，使用声纹传感器就可以采集设备的声纹信号，使用振动传感器就可以采集设备的振动信号，采集的声纹信号和振动信号应是同步的；
[0149]
步骤b、提取声纹数据的短时过零率、谱质心和短时能量，提取振动数据的偏度、峭度、波峰因子和形状因子；具体可参见实施例1的方法的步骤2；
[0150]
步骤c、将短时过零率、所述谱质心、所述短时能量、所述偏度、所述峭度、所述波峰因子和所述形状因子输入使用实施例1的方法构建的设备状态分类模型，设备状态分类模型输出设备状态。
[0151]
实施例5：一种设备状态分类模型的使用方法，包括以下步骤：
[0152]
步骤a、采集设备的声纹信号和振动信号；一般地，使用声纹传感器就可以采集设备的声纹信号，使用振动传感器就可以采集设备的振动信号，采集的声纹信号和振动信号应是同步的；
[0153]
步骤b、使用稀疏变换算法处理所述声纹数据和所述振动数据；具体可参见实施例2的方法的步骤2；
[0154]
步骤c、提取处理后的声纹数据的短时过零率、谱质心和短时能量，提取处理后的振动数据的偏度、峭度、波峰因子和形状因子；具体可参见实施例2的方法的步骤3；
[0155]
步骤d、将短时过零率、所述谱质心、所述短时能量、所述偏度、所述峭度、所述波峰因子和所述形状因子输入使用实施例2的方法构建的设备状态分类模型，设备状态分类模型输出设备状态。
[0156]
实施例6：一种存储有设备状态分类程序的计算机可读介质，所述设备状态分类程序被处理器时执行以下步骤：
[0157]
步骤a、获取所述设备的声纹信号和振动信号；获取的声纹信号和振动信号应是同
步的；
[0158]
步骤b、提取声纹数据的短时过零率、谱质心和短时能量，提取振动数据的偏度、峭度、波峰因子和形状因子；具体可参见实施例1的方法的步骤2；
[0159]
步骤c、将所述短时过零率、所述谱质心、所述短时能量、所述偏度、所述峭度、所述波峰因子和所述形状因子输入实施例1的方法构建的设备状态分类模型，所述设备状态分类模型输出设备状态。
[0160]
实施例7：一种存储有设备状态分类程序的计算机可读介质，所述设备状态分类程序被处理器时执行以下步骤：
[0161]
步骤a、获取设备的声纹信号和振动信号；获取的声纹信号和振动信号应是同步的；
[0162]
步骤b、使用稀疏变换算法处理所述声纹数据和所述振动数据；具体可参见实施例2的方法的步骤2；
[0163]
步骤c、提取处理后的声纹数据的短时过零率、谱质心和短时能量，提取处理后的振动数据的偏度、峭度、波峰因子和形状因子；具体可参见实施例2的方法的步骤3；
[0164]
步骤d、将短时过零率、所述谱质心、所述短时能量、所述偏度、所述峭度、所述波峰因子和所述形状因子输入使用实施例2的方法构建的设备状态分类模型，设备状态分类模型输出设备状态。
[0165]
实施例8：一种设备状态辨识装置，包括处理器和存储有设备状态分类程序的计算机可读介质，所述处理器用于执行实施例6或实施例7所述的设备状态分类程序。
[0166]
实施例9：一种设备状态辨识装置，包括声纹传感器、振动传感器、处理器、存储有设备状态分类程序的计算机可读介质和输出模块，声纹传感器用于采集所述设备的声纹信号，振动传感器用于采集所述设备的振动信号，处理器用于接收所述声纹信号和振动信号，执行所述设备状态分类程序，获得设备状态信息，输出模块用于输出设备状态信息。输出模块可以是显示器、喇叭、灯、打印机、传真机等。
[0167]
使用时，声纹传感器输出声纹信号，振动传感器输出振动信号，处理器接收声纹信号和振动信号后，在执行设备状态分类程序时，将声纹信号和振动信号作为设备状态分类程序的输入量，设备状态分类程序输出设备状态信息，输出模块输出设备状态信息。
[0168]
上面结合附图和实施例对本发明作了详细的说明。应当明白，实践中无法穷尽地说明所有可能的实施方式，在此通过举例说明的方式尽可能的阐述本发明得发明构思。在不脱离本发明的发明构思、且未付出创造性劳动的前提下，本技术领域的技术人员对上述实施例中的技术特征进行取舍组合、具体参数进行试验变更，或者利用本技术领域的现有技术对本发明已公开的技术手段进行常规替换形成的具体的实施例，均应属于为本发明隐含公开的内容。