1. 引言
随着制造业的飞速发展,电子设备的内部构造变得越来越复杂,这导致设备发生故障的几率也相应增加。一旦电子设备出现故障,企业将面临重大的经济损失。一般情况下,电子设备的故障并不会直接显现,而是需要通过分析相关数据来进行判断,这无疑增加了电子设备的维修时间难度。若能提前对电子设备的故障进行诊断,不仅能保障电子设备的安全稳定运行,而且还可延长电子设备的使用寿命,使得企业的运行成本得到有效控制。因此,对电子设备的故障进行有效诊断一直都是重要的研究课题。
目前,国内外已有许多学者利用不同数学模型诊断电子设备故障,取得了较为丰富的研究成果。例如,文献[1]给出了一种基于小波分析和神经网络的传感器故障诊断方法;文献[2]提出了一种BP神经网络的故障诊断模型;文献[3]使用卷积神经网络提取振动信号特征,利用EEMD和近似熵构建信号特征向量,并通过训练BP神经网络创建了水电机组故障识别器;文献[4]提出了一种基于时间频率分析(时频分析)的故障诊断方法;文献[5]为了精确地识别变转速电机轴承的故障,提出了一种优化的时频分析策略;文献[6]使用自动编码器对模拟电路的初始故障信号进行降噪处理,然后将处理后的信号输入到深度信念网络中以提取特征;文献[7]为了有效地检测电子电路中的问题,研发了一种故障检测系统;文献[8]提出了一种基于小波分析的电磁阀故障诊断方法;文献[9]针对航空电子设备故障诊断问题,采用了聚类方法构建了多个模型,并开发了一个实用的故障诊断系统;文献[10]提出了一种MCKD方法,该方法结合了粒子群优化算法,能够更准确、更高效地诊断设备的运行状态;文献[11]提出了一种创新的基于序列平均变化率的灰色诊断模型;文献[12]采用动态贝叶斯网络作为主要工具识别出故障征兆并获取相关训练数据;文献[13]采用加权灰色关联分析,提出了一种改进的复杂电子设备故障诊断方法;文献[14]运用短时傅立叶变换(STFT)和连续小波变换(CWT)进行机电系统故障诊断等。
综合分析上述文献发现,利用不同数学模型诊断电子设备故障各有优势,但有些数学模型的诊断精度还可进一步提高。例如,利用单一的小波分析模型或神经网络模型对电子设备故障进行诊断虽然能取得较好的效果,但其诊断精度可能并不令人满意。为此,本文将考虑结合小波分析模型和神经网络模型的优点,构建一种适用于电子设备故障诊断的组合模型,以期提高电子设备的诊断精度,为电子设备的故障诊断提供一种新选择。
2. 预备知识
2.1. 小波分析的基本原理及其MATLAB实现
小波分析[1]是一种强大的信号处理技术,它的核心原理在于将复杂的信号分解成不同尺度的频率成分,以便于进行精细的时频分析和特征提取。这种方法的关键在于多尺度分析,它仅对信号的低频部分进行进一步的细分,而对高频部分则不作进一步处理。下面简要介绍小波多分辨率分析。
函数
的连续小波变换定义为:
(1)
式中,
为尺度因子;
为平移因子。
在式(1)中令
,
,得离散小波变换为:
(2)
式中,
为尺度指标;k为平移指标;
为离散小波系数。
通常情况下,离散化参数可取为
,
,此时得二进离散小波变换为:
(3)
式中,
为二进离散小波变换系数,小波分析的放大倍数
会随着尺度指标j改变,这也是小波分析的多分辨率特性。
利用MATLAB进行小波分析的步骤为:
Step 1:准备工作。确保已经将“Wavelet Toolbox”工具箱添加到MATLAB的路径中,可使用“waveletdenoise”等函数进行小波分析,且准备好需要分析的数据。
Step 2:选择小波基函数。MATLAB的小波函数库中有许多不同类型的小波基函数可供选择,如Haar小波、Daubechies小波、Symlet小波等。
Step 3:执行小波变换。使用“wavedec”函数执行小波变换,将时间序列数据转换为小波系数。该函数的参数包括时间序列数据、小波基函数和分解层数,调用格式为:
[C, L] = wavedec(signal, level, “waveletname”)
其中,“signal”是需要分解的信号,“level”是分解的层数,“waveletname”是选择的小波基函数的名称。“C”是一个包含近似系数和细节系数的向量,“L”是一个向量,包含近似系数和细节系数对应的长度。
Step 4:提取近似系数和细节系数。使用“appcoef”函数提取近似系数,“detcoef”函数提取细节系数,调用格式为:
ca = appcoef(C, L, “waveletname”, level); %提取近似系数
cd = detcoef(C, L, level); %提取细节系数
其中,“C”和“L”是“wavedec”函数的输出,“waveletname”是小波基函数的名称,“level”是提取系数的层数。
2.2. RBF的基本原理及其MATLAB实现
神经网络[15]是一种基于人工智能的高效模型,它通过自动学习数据,实现对数据的深度分析和诊断。神经网络可以被视为一个由众多神经元模型相互连接构成的复杂网络,具有强大的学习、模式识别和诊断能力。神经网络由输入层、隐藏层和输出层三个结构组成。神经网络的组成结构如图1所示。
这里主要介绍RBF神经网络。RBF神经网络[16]是一种模拟人脑结构和学习机制的先进算法。其核心理念在于采用径向基函数作为隐藏层的激活函数,每个输入特征均拥有独立的径向基中心。在训练过程中,RBF网络着重调整网络中心和半径,力求最小化网络输出与实际输出之间的差异,从而提升诊断精度和整体性能。
Figure 1. Composition structure of neural networks
图1. 神经网络的组成结构
RBF神经网络学习算法需要求解的参数有三个:基函数的中心向量W、方差σ以及隐含层到输出层的权值。根据选取的径向基函数中心的方法不同,RBF神经网络有不同的学习方式。
RBF神经网络学习方法由两个阶段组成:一是自组织学习阶段,此阶段为无导师学习过程,求解隐藏层基函数的中心向量与方差;二是有导师学习阶段,此阶段求解隐藏层到输出层之间的权值。
RBF神经网络的可表示为:
(4)
式中,y为网络的输出;
为第i个RBF神经元的权重;
为径向基函数,它表示输入向量x与第i个神经元的中心
之间的欧氏距离
的径向对称函数;
表示对所有RBF神经元进行求和。
径向基函数φ通常选择高斯函数,其表达式为:
(5)
式中,
为输入向量x与神经元中心
之间的欧氏距离;
为高斯函数的宽度参数,它决定了函数的形状和宽度。
RBF神经网络是一种三层前馈网络,其输入层将输入向量x传递给隐藏层,隐藏层使用径向基函数对输入进行变换并计算输出值,输出层将隐藏层的输出进行线性组合得到最终输出y。通过调整权重
和中心
,RBF神经网络可以学习并逼近任意非线性函数。
在MATLAB中训练和使用RBF神经网络的步骤为:
Step 1:准备数据。准备输入数据矩阵“X”和目标输出数据矩阵“T”。
Step 2:创建RBF网络。使用“newrb”或“newrbe”函数创建一个RBF神经网络,调用规则为:
net = newrb(X, T, GOAL, SPREAD, MN, DF);
其中,“X”为输入数据矩阵,“T”为目标输出数据矩阵,“GOAL”为网络性能目标,“SPREAD”为径向基函数的扩展速度,“MN”为神经元数量的最小值,“DF”为分布密度函数。
Step 3:训练网络。使用“train”函数训练RBF网络,调用格式为:
[net, tr] = train(net, X, T);
其中,“net”为待训练的RBF神经网络,“X”为输入数据矩阵,“T”为目标输出数据矩阵,“tr”为训练记录信息。
Step 4:测试网络。使用“sim”函数测试训练好的网络的性能,调用格式为:
Y = sim(net, Xtest);
其中,“net”为训练好的RBF神经网络,“Xtest”为测试输入数据矩阵,“Y”为网络的预测输出。
Step 5:使用网络进行分析。使用训练好的RBF神经网络对新数据进行分析,调用格式为:
prediction = sim(net, newdata);
其中,“net”为训练好的RBF神经网络,“newdata”为新输入数据,“prediction”为网络的输出。
3. 基于小波分析与神经网络的电子故障诊断方法
为精确高效地对电子设备故障进行诊断,本文结合小波分析和神经网络模型的优点,构建出适用于电子设备故障诊断的组合模型。利用该组合模型进行电子设备故障诊断的基本流程如图2示。
Figure 2. Basic process of using composite models for electronic device fault diagnosis
图2. 利用组合模型进行电子设备故障诊断的基本流程
3.1. 数据预处理
为了使得故障样本数据被限定在一定的范围内,从而消除奇异样本数据导致的不良影响,需将收集到的电子设备故障样本数据进行归一化预处理。由于电子设备故障样本数据是离散的,一般不符合正态分布,故使用最小–最大归一化方法进行数据与处理。最小–最大归一化的计算公式为:
(6)
其中,y为标准化后的数据;x为原数据;min为数据最小值;max为数据最大值。
3.2. 对故障样本数据进行小波分析提取小波系数
小波分析能够通过对故障信号的分解,提取出蕴含故障特征的小波系数。这些系数独一无二地反映了故障的类型、程度和潜在的发展趋势。因此,小波分析不仅捕捉到了隐藏在复杂信号中的关键信息,还为故障的早期识别、预警和精准修复提供了有力支持。它可以大大提升电子设备故障样本数据分析的准确性和效率。
使用MATLAB对电子设备故障样本数据进行小波分析并提取小波系数的主要步骤为:
Step 1:准备故障样本数据。导入归一化后的电子设备故障样本数据。
Step 2:选择小波函数。选择一个合适的小波函数作为基函数。由于Daubechies小波具有紧支撑性、对称性和正交性等特点,并且能够有效地进行多分辨率分析,故Daubechies小波可以更加准确地提取故障特征、识别故障类型并预测故障发展趋势,从而为电子设备的故障诊断和预防性维护提供有力支持。因此,本文选取Daubechies小波作为基函数。
Step 3:进行小波分解。将归一化后的电子设备故障样本数据、Daubechies小波和尺度参数作为输入。本文将尺度参数设置为2,也就是对故障样本数据进行小波分解的层数为2。
Step 4:提取小波系数。从小波系数矩阵中提取需要的小波系数。本文提取的小波系数为近似系数、第一层细节系数和第二层细节系数。
3.3. 用神经网络诊断故障类型
RBF神经网络的核心思想是通过训练样本集来建立一个模型,该模型能够将输入特征映射到输出类别或故障类型上。在电子设备故障诊断中,这意味着可以将监测到的样本数据作为输入向量x,通过RBF神经网络进行处理和分析,最终输出诊断结果,即故障类型或异常状态。
利用MATLAB通过RBF神经网络进行故障诊断的主要步骤为:
Step 1:数据导入。将故障样本数据的小波系数导入。
Step 2:创建RBF神经网络。使用“newrb”函数创建RBF神经网络。
Step 3:训练RBF神经网络。使用故障样本数据的小波系数训练RBF神经网络。
Step 4:测试RBF神经网络。使用测试数据对训练好的RBF神经网络进行测试。
Step 5:使用RBF神经网络进行故障诊断。将待诊断的设备数据输入到训练好的RBF神经网络中,神经网络将输出对应的故障类型。
4. 仿真实验
文献[2]以某电子设备为诊断对象,给出了6个测试点10个板卡故障的样本数据,如表1所示,其中Vi/V (i = 1, 2, 3, 4)表示四个时间段的电压。
Table 1. Fault sample data
表1. 故障样本数据
编号 |
V1/V |
V2/V |
V3/V |
V4/V |
故障类型 |
1 |
10.253 |
14.135 |
10.624 |
11.311 |
MA-1 |
2 |
0.536 |
14.782 |
0.325 |
0.259 |
MA-2 |
3 |
5.126 |
8.395 |
13.698 |
14.687 |
MB-1 |
4 |
0.568 |
4.865 |
9.658 |
14.568 |
MB-2 |
5 |
14.325 |
11.459 |
0.265 |
5.368 |
MC-1 |
6 |
13.684 |
12.084 |
0.657 |
4.865 |
MC-2 |
7 |
13.956 |
13.652 |
9.657 |
11.527 |
MD |
8 |
9.576 |
12.651 |
12.357 |
0.826 |
ME |
9 |
5.367 |
4.795 |
12.657 |
0.264 |
MF |
10 |
14.965 |
10.532 |
4.976 |
9.853 |
MG |
11 |
5.013 |
14.965 |
10.021 |
4.986 |
正常 |
由于表1的数据量太少,不利于RBF神经网络训练。为此,本文将表1中的电压数据分别扩大1.5倍和缩小0.5倍,得到故障的模拟数据如表2所示。
对表1与表2中的数据进行归一化处理,结果如表3所示。
将表3中的数据用MATLAB进行小波分析,提取近似系数、第一层细节系数和第二层细节系数,其中第一层细节系数有两个,第二节细节系数有一个,结果如表4所示。
将表4中的前25组数据用于训练RBF神经网络,剩余8组数据用于检验故障诊断的准确度。剩余8组数据的故障诊断检验结果如表5所示。
由表5可知,RBF神经网络的诊断准确率达到了75%,表现良好。在具体样本中,编号26、28、29、30、32和33的诊断结果完全正确,这进一步印证了RBF神经网络在处理特定数据集时的有效性。然而,编号27和编号31的诊断结果错误,这可能源于多种因素,如网络训练不足、数据中存在噪声或样本分布不平衡等。
Table 2. Simulated data of faults
表2. 故障的模拟数据
编号 |
V1/V |
V2/V |
V3/V |
V4/V |
编号 |
V1/V |
V2/V |
V3/V |
V4/V |
12 |
15.3795 |
21.2025 |
15.936 |
16.9665 |
23 |
5.1265 |
7.0675 |
5.312 |
5.6555 |
13 |
0.804 |
22.173 |
0.4875 |
0.3885 |
24 |
0.268 |
7.391 |
0.1625 |
0.1295 |
14 |
7.689 |
12.5925 |
20.547 |
22.0305 |
25 |
2.563 |
4.1975 |
6.849 |
7.3435 |
15 |
0.852 |
7.2975 |
14.487 |
21.852 |
26 |
0.284 |
2.4325 |
4.829 |
7.284 |
16 |
21.4875 |
17.1885 |
0.3975 |
8.052 |
27 |
7.1625 |
5.7295 |
0.1325 |
2.684 |
17 |
20.526 |
18.126 |
0.9855 |
7.2975 |
28 |
6.842 |
6.042 |
0.3285 |
2.4325 |
18 |
20.934 |
20.478 |
14.4855 |
17.2905 |
29 |
6.978 |
6.826 |
4.8285 |
5.7635 |
19 |
14.364 |
18.9765 |
18.5355 |
1.239 |
30 |
4.788 |
6.3255 |
6.1785 |
0.413 |
20 |
8.0505 |
7.1925 |
18.9855 |
0.396 |
31 |
2.6835 |
2.3975 |
6.3285 |
0.132 |
21 |
22.4475 |
15.798 |
7.464 |
14.7795 |
32 |
7.4825 |
5.266 |
2.488 |
4.9265 |
22 |
7.5195 |
22.4475 |
15.0315 |
7.479 |
33 |
2.5065 |
7.4825 |
5.0105 |
2.493 |
Table 3. Normalized data
表3. 归一化后的数据
编号 |
V1/V |
V2/V |
V3/V |
V4/V |
故障类型 |
编号 |
V1/V |
V2/V |
V3/V |
V4/V |
故障类型 |
1 |
0.4536 |
0.6275 |
0.4702 |
0.5010 |
MA-1 |
18 |
0.9322 |
0.9118 |
0.6432 |
0.7689 |
MD |
2 |
0.0182 |
0.6565 |
0.0088 |
0.0058 |
MA-2 |
19 |
0.6378 |
0.8445 |
0.8247 |
0.0497 |
ME |
3 |
0.2239 |
0.3704 |
0.6080 |
0.6523 |
MB-1 |
20 |
0.3549 |
0.3165 |
0.8449 |
0.0119 |
MF |
4 |
0.0196 |
0.2122 |
0.4269 |
0.6469 |
MB-2 |
21 |
1.0000 |
0.7021 |
0.3286 |
0.6564 |
MG |
5 |
0.6361 |
0.5076 |
0.0061 |
0.2347 |
MC-1 |
22 |
0.3311 |
1.0000 |
0.6677 |
0.3293 |
正常 |
6 |
0.6073 |
0.5356 |
0.0236 |
0.2122 |
MC-2 |
23 |
0.2239 |
0.3109 |
0.2322 |
0.2476 |
MA-1 |
7 |
0.6195 |
0.6059 |
0.4269 |
0.5107 |
MD |
24 |
0.0062 |
0.3254 |
0.0015 |
0.0000 |
MA-2 |
8 |
0.4233 |
0.5610 |
0.5479 |
0.0312 |
ME |
25 |
0.1090 |
0.1823 |
0.3011 |
0.3232 |
MB-1 |
9 |
0.2347 |
0.2090 |
0.5613 |
0.0060 |
MF |
26 |
0.0069 |
0.1032 |
0.2106 |
0.3206 |
MB-2 |
10 |
0.6647 |
0.4661 |
0.2172 |
0.4357 |
MG |
27 |
0.3151 |
0.2509 |
0.0001 |
0.1145 |
MC-1 |
11 |
0.2188 |
0.6647 |
0.4432 |
0.2176 |
正常 |
28 |
0.3008 |
0.2649 |
0.0089 |
0.1032 |
MC-2 |
12 |
0.6833 |
0.9442 |
0.7082 |
0.7544 |
MA-1 |
29 |
0.3069 |
0.3000 |
0.2105 |
0.2524 |
MD |
13 |
0.0302 |
0.9877 |
0.0160 |
0.0116 |
MA-2 |
30 |
0.2087 |
0.2776 |
0.2710 |
0.0127 |
ME |
14 |
0.3387 |
0.5584 |
0.9148 |
0.9813 |
MB-1 |
31 |
0.1144 |
0.1016 |
0.2778 |
0.0001 |
MF |
15 |
0.0324 |
0.3212 |
0.6433 |
0.9733 |
MB-2 |
32 |
0.3295 |
0.2302 |
0.1057 |
0.2149 |
MG |
16 |
0.9570 |
0.7644 |
0.0120 |
0.3550 |
MC-1 |
33 |
0.1065 |
0.3295 |
0.2187 |
0.1059 |
正常 |
17 |
0.9139 |
0.8064 |
0.0384 |
0.3212 |
MC-2 |
|
|
|
|
|
|
Table 4. Wavelet coefficients
表4. 小波系数
编号 |
近似 系数 |
细节 系数1a |
细节 系数1b |
细节 系数2 |
故障 类型 |
编号 |
近似 系数 |
细节 系数1a |
细节 系数1b |
细节 系数2 |
故障 类型 |
1 |
1.0262 |
−0.1230 |
−0.0218 |
0.0549 |
MA-1 |
18 |
1.6281 |
0.0144 |
−0.0889 |
0.2160 |
MD |
2 |
0.3447 |
−0.4513 |
0.0021 |
0.3301 |
MA-2 |
19 |
1.1784 |
−0.1462 |
0.5480 |
0.3040 |
ME |
3 |
0.9273 |
−0.1036 |
−0.0313 |
−0.3330 |
MB-1 |
20 |
0.7641 |
0.0272 |
0.5890 |
−0.0927 |
MF |
4 |
0.6528 |
−0.1362 |
−0.1556 |
−0.4210 |
MB-2 |
21 |
1.3436 |
0.2106 |
−0.2318 |
0.3586 |
MG |
5 |
0.6923 |
0.0909 |
−0.1616 |
0.4515 |
MC-1 |
22 |
1.1641 |
−0.4730 |
0.2393 |
0.1671 |
正常 |
6 |
0.6894 |
0.0507 |
−0.1334 |
0.4536 |
MC-2 |
23 |
0.5073 |
−0.0615 |
−0.0109 |
0.0275 |
MA-1 |
7 |
1.0815 |
0.0096 |
−0.0593 |
0.1439 |
MD |
24 |
0.1666 |
−0.2257 |
0.0011 |
0.1651 |
MA-2 |
8 |
0.7817 |
−0.0974 |
0.3654 |
0.2026 |
ME |
25 |
0.4578 |
−0.0518 |
−0.0156 |
−0.1665 |
MB-1 |
9 |
0.5055 |
0.0182 |
0.3927 |
−0.0618 |
MF |
26 |
0.3207 |
−0.0681 |
−0.0778 |
−0.2106 |
MB-2 |
10 |
0.8919 |
0.1404 |
−0.1545 |
0.2390 |
MG |
27 |
0.3403 |
0.0454 |
−0.0809 |
0.2257 |
MC-1 |
11 |
0.7722 |
−0.3153 |
0.1595 |
0.1114 |
正常 |
28 |
0.3389 |
0.0254 |
−0.0667 |
0.2268 |
MC-2 |
12 |
1.5451 |
−0.1845 |
−0.0327 |
0.0825 |
MA-1 |
29 |
0.5349 |
0.0049 |
−0.0296 |
0.0720 |
MD |
13 |
0.5228 |
−0.6771 |
0.0031 |
0.4952 |
MA-2 |
30 |
0.3850 |
−0.0487 |
0.1826 |
0.1013 |
ME |
14 |
1.3966 |
−0.1554 |
−0.0470 |
−0.4995 |
MB-1 |
31 |
0.2470 |
0.0091 |
0.1964 |
−0.0310 |
MF |
15 |
0.9851 |
−0.2042 |
−0.2333 |
−0.6315 |
MB-2 |
32 |
0.4402 |
0.0702 |
−0.0772 |
0.1196 |
MG |
16 |
1.0442 |
0.1362 |
−0.2425 |
0.6772 |
MC-1 |
33 |
0.3803 |
−0.1577 |
0.0798 |
0.0557 |
正常 |
17 |
1.0400 |
0.0760 |
−0.2000 |
0.6804 |
MC-2 |
|
|
|
|
|
|
Table 5. Fault diagnosis test results
表5. 故障诊断检验结果
编号 |
RBF |
诊断故障 |
实际故障 |
26 |
4.467 |
MB-2 |
MB-2 |
27 |
6.1969 |
MC-2 |
MC-1 |
28 |
6.7419 |
MC-2 |
MC-2 |
29 |
7.3042 |
MD |
MD |
30 |
8.4776 |
ME |
ME |
31 |
10.9565 |
MG |
MF |
32 |
10.0514 |
MG |
MG |
33 |
11.5465 |
正常 |
正常 |
5. 结语
电子设备在持续运行过程中常常面临各种故障的挑战。因此,如何快速、准确地对电子设备的故障进行诊断具有重要的应用价值。本文提出了一种基于小波分析和RBF神经网络的电子设备故障诊断方法,该方法结合了小波分析与神经网络在数据分析中的优点,能有效地实现电子设备故障诊断。利用MATLAB进行仿真实验,结果表明,本文提出的方法在对电子设备进行故障诊断时具有较高的准确率,具有较高的应用价值。
基金项目
曲阜特霓电子科技有限公司横向课题“电子设备数据的分析与可视化”资助。
NOTES
*通讯作者。