首页/文章/ 详情

中科院一区Top | 用于多目标域故障诊断的模糊跨域矩阵机方法研究

3月前浏览1431

论文基本信息

论文题目:A fuzzy cross domain matrix machine for fault diagnosis under multi-objective domain

论文期刊:Engineering Applications of Artificial Intelligence

论文日期:2026

论文链接:https://doi.org/10.1016/j.engappai.2025.113368

作者:Haiyang Pan [a,b,*], Chunan Chen [a], Tong Chen [a], Jian Cheng [a], Jinde Zheng [a], Shuchao Deng [a,c]

机构:

[a] School of Mechanical Engineering, Anhui University of Technology, Ma’anshan, 243032, China

[b] Anhui Key Laboratory of Mine Intelligent Equipment and Technology, Anhui University of Science and Technology, Huainan, 232001, China

[c] Anhui Gongda Information Technology Co., Ltd, Ma’anshan, 243002, China

作者简介:

潘海洋,安徽工业大学副教授,湖南大学机械工程博士,西安交通大学访问学者,安徽省优秀青年研究生导师,连续4年入选“全球前2%顶尖科学家榜单”(2022-2025),主持与参与国家重点研发计划、国家自然科学基金、安徽省自然科学基金、安徽省优秀青年教师培育重点项目等20余项;以第一作者或通讯作者在TII、MSSP、《机械工程学报》等期刊发表论文80余篇,领域内高被引5篇;获安徽省科学技术奖(自然科学)二等奖、安徽省科学技术奖(科技进步)三等奖、中国振动工程学会科学技术奖二等奖各1项。

陈春安,安徽工业大学在读硕士研究生,在Engineering Applications of Artificial Intelligence、Mechanis m and Machine Theory、IEEE Transactions on Instrumentation and Measurement、《机械工程学报》等期刊发表高水平学术论文10余篇,获2025年度研究生国家奖学金,安徽省优秀毕业生等荣誉。

程健,安徽工业大学资格教授,硕士生导师,湖南大学机械工程博士毕业,参与国家重点研发计划项目、国家自然科学基金、安徽省自然科学基金、校企合作课题等8项;担任MSSP、ESWA、RESS、TIM和ISA等20多个期刊审稿人;以第一作者或通讯作者发表SCI/EI论文27篇。

郑近德,博士生导师,安徽工业大学科研院副院长,安徽省学术和技术带头人后备人选,安徽省领军人才特聘教授,安徽省“青年皖江学者”,省级优秀青年研究生导师;主持国家自然科学基金(3项),安徽省高校杰出青年科研项目、安徽省自然科学基金等20余项;在国内外期刊发表学术论文120余篇,授权发明专利5项,领域内高被引论文4篇,获得安徽省科学技术奖二等奖2项(第一完成人和第六完成人),中国振动工程学会科技奖二等奖(第一完成人)。

目录

摘要

1 引言

2 相关工作

    2.1 支持矩阵机

    2.2 直接多任务支持向量机

模糊跨域矩阵机原理

    3.1 FCDMM模型

    3.2 参数选择指南

实验分析

    4.1 基于FCDMM的机械故障诊断流程

    4.2 实验1

    4.3 实验2

    4.4 计算复杂性分析

    4.5 应用分析

5 结论

摘要

由于传统矩阵分类器对异常样本敏感且无法利用跨域信息,机械故障诊断面临重大挑战。为了解决这些局限性,本文提出了一种新颖的模糊跨域矩阵机(Fuzzy Cross Domain Matrix Machine, FCDMM),它将模糊隶属度建模与多任务学习相结合,以增强鲁棒性并实现跨相关任务的知识迁移。FCDMM构建了基于模糊的决策超平面,并设计了自适应边界因子以减轻异常值的影响。在多个故障数据集上进行的大量实验表明,FCDMM在轴承故障和齿轮故障上分别取得了98.86%和99.57%的显著准确率,显著优于最先进的方法。

关键词:故障诊断、模糊跨域矩阵机器、多任务学习、异常样本

1 引言

随着机械设备在现代工业中扮演着越来越不可或缺的角色,机械设备已渗透到各个行业,成为不可或缺的核心力量。不幸的是,在长期连续运行过程中,机械设备不可避免地受到多种内外因素交织的影响,进而影响工业生产的连续性和稳定性,并对人身安全构成潜在威胁。因此,开展机械故障诊断研究以实现早期检测和诊断已成为一个紧迫的问题。

近年来,随着人工智能技术的快速发展,机器学习作为一种强大的数据分析工具,为机械故障诊断提供了新的研究思路和方法。通过从海量设备运行数据中提取特征并建立智能诊断模型,机器学习可以有效监控设备运行状态。例如,Ali等人提出了一种基于动态机器学习的新型可解释融合方法,以解决工业和化工过程系统中的过程监控问题。此外,他们还提出了一种基于鲁棒机器学习和低频生成单变量控制图的创新分布式典型相关分析,利用小波变换监控工业和制造系统。另外,一些学者的研究集中于电机故障诊断,例如Ahmadpour等人通过三维建模电机并使用希尔伯特-黄变换来监测故障特征。然而,当前主流的机械诊断方法是基于单任务诊断,其侧重于当前特征信息的利用当面对多种运行工况时,单一的机器学习模型可能表现不佳。原因是不同类型的故障可能具有相似的特征,或者同一故障在不同运行工况下可能表现出不同的模式。在复杂系统故障诊断中,需要同时完成故障识别、定位和严重程度评估等多个相关任务。然而,工业现场标注数据通常稀缺且成本高昂,因此实现多任务故障诊断尤为必要。为了解决这个问题,多任务学习( Multi-Task Learning,MTL)作为一种新兴的机器学习范式逐渐受到关注。与传统的单任务学习不同,MTL充分利用任务之间的相关性,通过共享表示捕获共同特征,从而降低过拟合风险并提高学习效率。在MTL中,任务相关性是基础,它可以是数据分布、特征表示或决策边界上的相似性。共享表示是实现多任务学习的关键,它可以是低级的原始特征或高级的抽象特征。同时,MTL也具有一定的正则化效果,有助于防止模型在单个任务上过拟合。因此,MTL能够通过同时学习多个相关任务并利用任务间的共享信息来提高泛化能力和鲁棒性。

目前,基于MTL的机器学习方法已经被提出并得到广泛应用,例如多任务高斯过程(Multi-Task Gaussian Process,MTGP)、v-软间隔多任务学习逻辑回归( v-Soft Margin Multi-Task Learning logistic regression,v-S MMTL-LR)等。此外,Du等人提出了多任务深度强化学习,通过MTL优化应用深度强化学习方法的效率。Zhang等人提出了信息最大化多任务聚类(Information Maximization Multi-task Clustering,IMMC),以实现任务内聚类和跨任务信息传递。Xie等人通过将孪生支持向量机纳入MTL框架,提出了直接多任务孪生支持向量机(Direct Multitask Twin SVM,DMTSVM)。鉴于支持向量机(Support Vector Machines, SVM)在多任务分类中的有效性,一些学者开始关注SVM中的MTL算法。例如,Mei等人提出了多任务最小二乘孪生支持向量机(Multi-Task Least Squares Twin SVM,MTLS-TWSVM),通过引入最小二乘损失提高了模型的计算速度。Liu等人在DMTSVM中引入了不敏感损失,提出了多任务非平行支持向量机(Multi-Task Non-Parallel SVM, MTNPSVM),有效避免了矩阵求逆运算。An等人鉴于DMTSVM仅考虑经验风险最小化,提出了多任务孪生有界支持向量机(Multi-Task Twin Bound SVM,MT-TBSVM),通过加入正则化项来减少结构风险。Pang等人考虑将TSVM的两个正负超平面设计为超球面,并提出了多任务孪生超球面支持向量机(Multi-Task Twin Hypersphere SVM,MTTHSVM),以进一步提高模型的分类性能。然而,上述方法都使用向量形式的数据作为模型输入。当提取的特征信号表示为矩阵数据时,通常采取的措施是将矩阵重塑为向量。这种重塑可能导致矩阵特征数据中固有结构信息的丢失,并导致维数灾难。

为了避免向量化的局限性,学者们提出了一些针对矩阵数据的模式识别方法,例如Wolf等人通过约束低秩矩阵提出了低秩支持向量机(Low Rank SVM,LRSVM)。Luo提出了支持矩阵机(Support Matrix Machine,S MM),其通过设计谱弹性网络来利用矩阵数据行与列之间的强相关性。鉴于S MM的优越性,相关学者提出了S MM的衍生算法,例如稀疏支持矩阵机(Sparse S MM,SS MM)、最优间隔分布矩阵机(Optimal-margin Distribution Matrix Machine,ODMM)和非平行有界支持矩阵机(Non-Parallel Bounded Support Matrix Machine,NPBS MM)。尽管上述方法能够从矩阵样本中提取结构信息,但它们都集中于当前任务的特征信息,无法实现多信息利用。基于此,研究开始研究多任务矩阵机,并相继提出了多任务协同增强矩阵机(Multi-Task Collaborative Enhancement Matrix Machine,MTCEMM)和多任务孪生弹球核矩阵分类器(Multi-Task Twin Pinball Kernel Matrix Classifier,MTTPKMC)。MTCEMM设计了维度增强项以保留原始信号的结构完整性,而MTTPKMC引入了弹球约束项以减轻对噪声的敏感性。然而,上述方法在建模过程中未能考虑异常样本对超平面的重要性,以及特征空间数据分布对模型的影响。基于此分析,我们提出了一种新的分类方法,称为模糊跨域矩阵机(Fuzzy Cross Domain Matrix Machine,FCDMM)。与多任务支持向量机不同,FCDMM不仅可以充分利用矩阵行与列之间的强相关信息,还可以更灵活地解决异常样本的误分类问题。实验表明,FCDMM通过使用两个数据集展现出卓越的分类性能。

本文的主要贡献概述如下:

(1) 设计了模糊跨域项,用于量化类别样本点的模糊隶属关系,并基于模糊逻辑构建决策超平面,以实现数据分布的非线性建模和优化。

(2) 引入了一对边界因子来控制间隔误差,从而降低了模型对异常样本的敏感性。

(3) FCDMM已在多个故障数据集上实现,结果表明它能有效减轻异常样本对超平面的影响。

本文其余部分安排如下。第二部分介绍了S MM和DMTSVM。第三部分具体介绍了FCDMM方法,第四部分通过大量实验验证了所提方法的优越性。第五部分给出结论。

2 相关工作

2.1 支持矩阵机

支持矩阵机S MM通过设计谱弹性网络来利用矩阵样本的内部信息,具有更鲁棒的分类性能。S MM通过结合铰链损失和核范数构建目标函数,如下示所示。

 

其中,    为权重矩阵,b为阈值,    为核范数,    为核范数系数,    为折衷参数。

S MM虽然保护了原始数据的结构信息,但是其构造的是一对平行超平面限制了其处理复杂数据的分类性能,并且其可利用的特征信息过于单一。

2.2 直接多任务支持向量机

DMTSVM通过建立多任务学习框架,使得原本的向量机模型能够实现多信息利用。假设存在T个任务,所有任务的正类样本表示为a,所有的任务的负类样本表示为b。第t个任务的正类样本与负类样本分别表示为    和    。DMTSVM假设所有任务的公共信息为    和    ,    和    表示任务t与公共信息的偏差。DMTSVM的目标函数如下:

 
 

其中,A=(a,e),B=(b,e),    ,    ,    和    为调整参数,为松弛参数。

DMTSVM虽然能够实现多信息的利用,但是其构建的铰链损失函数与S MM类似,都容易受离群点的影响。此外,DMTSVM在处理矩阵式数据集时,会破坏原始数据的结构信息。

3 模糊跨域矩阵机原理

3.1 FCDMM模型

在第二节分析的基础上,提出了一种能够实现跨域信息利用的矩阵形式的分类方法FCDMM。FCDMM通过构造非平行超平面来提高模型的分类性能和计算速度,同时设计FCT来实现跨域学习,减少异常样本对模型的影响。FCT量化类别样本点之间的模糊隶属关系,并基于模糊逻辑构造决策超平面,实现数据分布的非线性建模和优化。此外,FCDMM还引入了一对边界因子来控制边际误差,以降低模型对异常样本的敏感性。

假设存在T个任务,每个任务的训练集为  ,所有任务的正类样本表示为a,所有的任务的负类样本表示为b。第t个任务的正类样本表示为  ,第t个任务的负类样本表示为  。为了充分利用原始信号的结构信息,FCDMM定义了一种矩阵核技巧如下。

定义:假设存在矩阵样本    ,则  ,其中k为偏移参数,  。因此,定义        ,其中c为所有任务的总样本。此外,为了方便,对公共信息设定,        。对任务t设定,        。FCDMM的目标函数构建如下:

   

其中,  为折衷参数,  为自适应参数,用于控制误差边界,为尺度参数。上式的前两项为FCT,第三项为正则项,避免矩阵计算过于奇异的问题。    为属于任务的模糊隶属度值的对角矩阵,  其中为类别的中心,  为样本点距离中心的最大距离,  为非常小的正值。

由于上式互不影响且都是不等式约束,其构造拉格朗日函数如下式所示。

 

其中,    为非负拉格朗日乘子。对上式分别求        的次梯度如下式所示。

 

由上式可以得到

 

将上式代入,上式的对偶问题可以被确定为下式所示:

 

其中,      。 最后,通过求解上述两个二次规划问题,任务t的未知样本的标签由下式确定。

 

3.2 参数选择指南

从上式可以看出,FCDMM的训练结果与参数    和    直接相关。因此,参数    和    直接影响整个模型的分类精度。    是尺度参数,用于缩放边界,对模型的分类精度影响较小。    和    主要用于避免奇异矩阵计算,对模型的分类精度影响相对较小。在核技巧中,d通常取2,k是一个小的自然数,主要用于核矩阵样本,不影响模型的分类结果。综上所述    和    是模型学习性能最重要的因素。

基于以上分析,首先设置模型的核技巧参数。此处设 d=2 和 k=4。其次,设置尺度参数    (其他值也可接受)。由于上式分别建立了正负超平面且互不影响,因此设    和    此时,只有    和    两个参数影响模型。因此,在训练模型时,应主要考虑    和    的值。    和    的取值范围分别为    和    

模型主要受参数    和    影响,其他参数可由范围内的某个值替代。论文中设置的参数可进行微调,以进一步提高模型的分类结果。

4 实验分析

4.1 基于FCDMM的机械故障诊断流程

为了验证FCDMM的诊断性能,使用多目标故障数据集进行实验。同时,选择DMTSVM、MTNPSVM、MTCEMM和MTTPKMC方法与FCDMM进行比较,并使用准确率、Kappa系数、召回率、F1分数和精确率等评估指标来评估FCDMM的分类性能。FCDMM的诊断流程如下,图1显示了其诊断流程图。所有实验均在Windows 11系统上的MATLAB R2022a上运行,计算机硬件配置为第12代英特尔酷睿i5-12450H处理器,16 GB内存,NVIDIA GeForce RTX 3060 GPU。

步骤一:利用不同的加速度传感器采集多故障部位(轴承和齿轮)不同故障状态下的振动信号;

步骤二:采用辛几何相似变换(SGST)预处理振动信号构建多任务矩阵数据集,并设定训练集和测试集;

步骤三:将多个任务的训练集输入至FCDMM算法中,训练模型;

步骤四:将测试集输入FCDMM模型中,对输出的结果进行分析。

图1 FCDMM方法的流程图

4.2 实验1

为了验证FCDMM在多目标故障诊断中的有效性,使用山东科技大学的轴承和齿轮故障数据集进行了实验验证。实验数据来源于旋转机械传动部件故障植入试验台,轴承和齿轮的不同故障类型如图2所示,试验部件包括6205轴承和行星变速箱。轴承和齿轮不同故障状态和试验条件的具体参数如表1和表2所示。在实验中,轴承和齿轮每种工况的1024个采样点作为一个样本,每种故障类型构造200个样本,其中100个样本作为训练集,100个样本作为测试集。

图2 旋转机械传动部件故障植入试验台

表1 任务1轴承故障数据

表2 任务2齿轮故障数据

实验前,对维度为1024的原始信号通过SGST方法进行处理,得到3232的辛几何系数矩阵作为模型的输入样本。FCDMM模型中的矩阵核技巧的参数设定为k=1,d=2,尺度参数设定为    。为了方便,折衷参数设定为        ,并且分别从集 合和        中选取。图3给出了FCDMM在轴承数据集和齿轮数据集的参数寻优结果,可以看出当        的值分别为0.1和0.001时,FCDMM在两个数据集中均能取得最佳的分类结果。因此,FCDMM的参数选取为        

图3 FCDMM在不同参数下的分类结果

在优化的参数下,对FCDMM模型进行训练,并与DMTSVM、MTNPSVM、MTCEMM和MTTPKMC等方法的实验结果进行了比较。结果如表3和表4所示,图4给出了所提出方法的分类混淆矩阵。从表3和表4来看,FCDMM方法对滚动轴承故障数据和齿轮故障数据的识别准确率分别达到98.86%和99.57%,识别精度明显优于其他方法。图4示出FCDMM在任务1和任务2中分别有3个和2个错误分类的类别,并且错误分类的样本的数量相对较少。这是因为FCT能够量化类别样本点在特征空间中的模糊隶属关系,使得不确定样本更加清晰,从而提高分类精度。MTCEMM和MTTPKMC方法在任务1和任务2上的识别率分别达到97%和96%,但总体分类性能略逊于FCDMM,因为这两种方法只考虑了经验风险最小化的原则。DMTSVM和MTNPSVM在任务1中的分类正确率达到96%,但在任务2中的分类性能并不理想。这是因为DMTSVM和MTNPSVM虽然也可以利用任务之间的共同信息,但都破坏了样本的结构信息。FCDMM方法各评价指标的标准差最小,说明该方法具有更稳定的分类性能。

表3 任务1分类结果

表4 任务2分类结果

图4 FCDMM方法的分类混淆矩阵

为了避免测试集和训练集之间的紧密分布对测试结果的影响,将多任务样本置乱,重新选择新的样本进行训练和测试。验证进行了五次,每种方法的平均性能指标如图5所示。在两个故障数据集中,结果表明FCDMM在所有评估指标中获得了最高的分数,证明了其在多任务故障诊断中的有效性。

图5 五次验证的平均值

为了验证FCDMM对异常样本的分类性能,将随机振动信号添加到两个故障数据集的原始信号中来模拟异常样本,图6显示了正常样本和异常样本的辛几何系数矩阵样本的比较。每个故障类别增加10个异常样本,训练集和测试集各包含5个异常样本,实验结果如图7所示。显然,在异常样本分类的情况下,FCDMM的分类结果最好。该方法的优点在于它能够建立多任务模糊超平面,并利用FCT中的模糊隶属函数对异常样本进行加权,使其更接近模糊超平面。其他比较方法在建模过程中没有考虑异常样本对模型的影响,导致分类性能较差。

图6 正常样本与异常样本(异常样本有随机振动信号)的对比

图7 异常样本下的评价指标结果

4.3 实验2

为了进一步验证FCDMM的诊断性能,使用安徽工业大学的故障数据集进行实验验证。该数据集在轧机振动特性实验平台上采集,如图8所示。为了构建具有相关性的多任务故障数据,分别采集了不同故障状态下变速箱轴承和轧辊轴承的振动信号。具体参数见表5和表6。

图8 轧钢机振动特性实验台

表5 任务1:齿轮箱轴承的详细参数

表6 任务6:滚动轴承的详细参数

同样,在利用SGST方法对采集的振动信号进行预处理后,为每个故障状态选择200个样本进行分析,其中一半作为训练集,另一半作为测试集。将训练样本输入到所提出的方法和比较方法中,然后用测试样本对它们进行测试。结果如图9所示,FCDMM在两个数据集中的分类结果都是最好的,每个评价度量值都达到了98%以上。该方法的优点在于其基于模糊逻辑的FCT决策超平面实现了数据分布的非线性建模和优化,从而进一步提高了模型分类的稳健性。MTCEMM和MTTPKMC在两个任务中都表现出了良好的分类性能,但它们的建模没有考虑不确定样本对模型的影响,导致分类性能不如所提出的方法。虽然MTNPSVM和DMTSVM都建立了多任务学习框架,但它们都没有考虑矩阵样本的结构信息,而FCDMM设计的矩阵核技术可以有效地利用原始信号中嵌入的内在信息。

图9 每项任务的不同评价指标值

为了消除实验的随机性,对两个数据集进行置乱,并为五次验证运行随机生成新的训练样本和测试样本。表7和表8列出了各评价指标的均值和标准差。结果表明,在两个故障数据集中,FCDMM的评价准则的结果最高,标准差也最小,充分证明了FCDMM具有良好的跨域诊断性能和良好的分类稳定性。

表7 任务1验证结果

表8 任务2验证结果

同样,在上述两个故障数据集的每个故障类别中添加10个异常样本,以验证FCDMM的分类性能,两个任务分类的准确性如图10所示。从图10可以看出,FCDMM在两个故障数据集中仍然表现出最好的分类性能,进一步验证了该方法的优越性和有效性。这是因为FCDMM引入参数来控制边际误差,允许一些样本点位于区间内或错误分类,从而提高了模型对包含异常样本的数据集的稳健性。

图10 包含异常样本的评估指标结果

为了进一步分析FCDMM方法的优越性,采用Friedman检验和Nemenyi检验进行统计分析。采用Friedman检验作为初步分析,用来比较不同方法的分类性能。Nemenyi检验是一种后验方法,用于确定不同算法之间的差异,如下式所示。

 
 

其中,N为实验的数量,    为第i个算法的平均秩,n为算法的数量,    为Nemenyi检验的临界值。

对上述五种方法进行排名排序,获得的平均排名如表9所示。由上式可得    ,当自由度为(4,95)时,临界值为2.005。显然,    的实际值远大于其临界值,因此这五种方法之间存在显著差异。接下来,使用Nemenyi检验进行后续验证,由上式可计算得CD=1.23。图11显示了Nemenyi检验编码图,结果表明所提方法与其他方法之间存在显著差异。

表9 五种方法的平均秩

图11 Nemenyi测试编码结果

4.4 计算复杂性分析

FCDMM的主要计算成本是执行二次规划问题以求解上式,计算复杂度为    ,其中 N 是样本总数。此外,多任务下矩阵乘法的计算复杂度为    。因此,FCDMM的计算复杂度为    

实验2中的数据集用于时间成本验证,表10显示了每种方法的培训和测试时间。FCDMM、MTCEMM和MTTPKMC方法的分类时间较长,因为这些方法处理的是二维矩阵样本,而MTNPSVM和DMTSVM的目标是一维向量样本。然而,FCDMM的分类精度明显好于其他方法。因此,综合分析FCDMM方法的分类性能最好。FCDMM方法采用一对一的策略,因此其所需时间会随着样本类别的增加而增加。

表10 各方式计算时间成本

4.5 应用分析

在实际的工厂环境中,机械设备的故障诊断通常依赖于监控和数据采集(SCADA)系统。将FCDMM集成到这类系统中可以作为高级分析模块部署在云中或本地服务器上,从底层传感器网络(如振动加速度传感器)接收实时或定期上传的数据快照。具体地说,工厂现有SCADA系统的振动传感器连续采集设备运行数据,然后使用SGST等方法对原始振动信号进行预处理。将生成的特征矩阵输入训练好的FCDMM模型。该模型基于从多个任务中学习到的知识进行联合推理,并输出故障分类结果。诊断结果可以返回到SCADA系统的人机界面,触发报警或向维护人员提供维护建议,从而实现预测性维护。 在FCDMM和更简单的单任务方法(如标准支持向量机或单任务矩阵机)之间进行选择,本质上涉及到平衡诊断精度、计算资源和管理复杂性。当需要同时诊断变速箱中的轴承和齿轮故障,或者需要诊断同一生产线上同一类型但不同位置的多个设备时,FCDMM可以通过跨域信息共享显著提高每项任务的诊断准确率和稳健性。如果只需要诊断单个组件或不相关的多个任务,则使用单个任务分类器更合适。目前的FCDMM是一种批量学习模式,不适合在线学习。工业设备的运行条件和退化模式可能会随着时间的推移而变化(概念漂移),这需要使用新数据定期对模型进行再训练。因此,将FCDMM扩展到在线学习是未来重要的研究方向。

5 结论

提出了一种新的用于稳健多任务故障诊断的模糊跨域矩阵机器。实验结果表明,FCDMM在轴承故障数据集和齿轮故障数据集上的分类准确率分别达到了98.86%和99.57%,达到了最好的分类效果。这比表现最好的基线方法MTTPKMC提高了大约1.8-2.5个百分点,MTTPKMC在各自的任务中分别达到了97.02%和97.11%。此外,FCDMM表现出了卓越的稳定性,与所有其他方法相比,其在五次验证运行中的最低标准偏差(低于0.5%)。 然而,这种性能的提高伴随着计算需求的增加。训练过程涉及两个二次规划问题的求解和参数优化带来的时间消耗。因此,未来的工作将集中在开发自适应参数选择算法,如贝叶斯优化框架,以实现调整过程的自动化,并减少计算开销。此外,将FCDMM框架扩展到针对有限样本故障诊断的单一学习也将是一个重点研究方向。



编辑:曹希铭

校核:李正平、陈凯歌、赵栓栓、赵学功、白亮、任超、Tina、陈宇航、海洋、Kira、王金、赵诚、肖鑫鑫、张优

该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除



来源:故障诊断与python学习
ACTMechanical振动非线性化学旋转机械ANSAMATLAB海洋UM参数优化电机传动控制人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-05-08
最近编辑:3月前
故障诊断与python学习
硕士 签名征集中
获赞 87粉丝 152文章 327课程 0
点赞
收藏
作者推荐

中科院一区Top | 基于决策树启发模拟电路分类器的机械信号实时分类与故障诊断方法研究

在数字信号处理器上运行的传统机器学习算法需要进行模拟信号与数字信号之间的转换,这会导致不可避免的延迟和大量的能耗。本文作者研究了一种基于决策树启发的新型模拟电路分类器,用于实时机器信号分类和故障诊断,该方法具有低延迟响应时间和低计算成本等优势,供大家交流讨论。论文基本信息论文题目:Real-Time Classification and Fault Diagnosis of Machine Signals via a Novel Decision-Tree-Inspired Analog Circuit Classifier 论文期刊:Mechanical Systems and Signal Processing论文日期:2026论文链接:https://doi.org/10.1016/j.ymssp.2026.114074作者:Siliang Lu [a], Junyan Zhu [a], Ling Zheng [a], Juncai Song [b], Zhiyong Hu [a], Xiaoxian Wang [c,*], Dong Wang [d]机构:a College of Electrical Engineering and Automation, Anhui University, Hefei, Anhui 230601, PR Chinab College of Internet, Anhui University, Hefei, Anhui 230601, PR Chinac College of Electronics and Information Engineering, Anhui University, Hefei, Anhui 230601, PR Chinad Department of Industrial Engineering and Management, School of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai 200240, PR China作者简介:朱俊炎,安徽大学电气工程与自动化学院机械专业在读硕士生,在Mechanical Systems and Signal Processing、IEEE Transactions on Instrumentation and Measurement上发表论文,获批国家发明专利,获全国大学智能车竞赛、“西门子杯”中国智能制造挑战赛全国二等奖。王骁贤,安徽大学光电信息获取与防护技术全国重点实验室副教授、硕士生导师,于2024年获得中国科学技术大学工学博士学位。研究方向包括信息处理、在线检测、智能运维等。主持国家自然科学基金青年项目1项,以第一作者/通信作者发表SCI期刊论文20余篇。陆思良,安徽大学电气工程与自动化学院教授、博士生导师、副院长,研究方向为机电系统信号处理与智能运维,主持国家自然科学基金4项,国家电网、奇瑞汽车等企业委托开发项目10余项。入选安徽省优青、安徽省青拔、全球前2%顶尖科学家终身榜单,曾获上海市科学技术二等奖、安徽省自然科学二等奖、中国电工技术学会科学技术奖一等奖、基金委工材学部机械工程学科“优秀结题项目”。担任《IEEE Transactions on Instrumentation and Measurement》、《IEEE Sensors Journal》期刊Associate Editor。王冬,国家特聘青年专家,上海交大机动学院长聘副教授,博导。致力于智能运维与大数据分析研究,成果应用于航天五院、中国核动力、中国船舶、联想、华为、上海核工程、上海电气、西部数据等合作研究中。获中国振动工程学会青年科技奖、中国科技产业化促进会科学技术奖一等奖、中国发明协会“发明创业奖”金奖、上海教学成果一等奖等。担任第九、第十届中国振动工程学会副秘书长、IEEE Sensors Journal数据科学领域主编、JRSE预测与健康管理领域主编、MSSP、IEEE TIM、JDMD和Measurement副主编。指导学生获首届博士生国自然青年学生基础项目、上交“学术之星”提名奖2人、机动学院重点推荐A级科创全国特等奖和一等奖5次、会议Best Paper最佳论文和口头汇报奖5次等。目录摘要1 引言2 基于模拟电路的决策树分类器设计 2.1. 使用模拟电路实现的决策树算法 2.2. 模拟电路分类器的设计 2.3. 模拟电路分类器的参数优化3 不同数据集上的实验验证 3.1 模拟电路分类器的拓扑结构 3.2 无刷直流电机(BLDCM)轴承数据集 (1)实验设置及信号预处理 (2)性能验证 (3)抗噪性能验证 (4)分类精度的定量分析 (5)信号持续时间的影响 3.3 凯斯西储大学(CWRU)轴承数据集 3.4 开关磁阻电机(SRM)数据集 3.5 数据集的关键信息摘要4 电机故障信号实时分类 4.1 实验装置 4.2 实时信号处理的性能验证 4.3 变速条件下的性能验证 4.4 验证系统的响应时间和能耗变速条件下的性能验证 4.5 基于模拟电路的分类器的比较5 讨论 5.1 信号类型数量与滤波器数量之间的关系 5.2 集成电路实时故障诊断和显示 5.3 滤波器特性对系统输出的影响 5.4 噪声干扰在实际场景中的影响6 结论摘要在数字信号处理器(Digital Signal Processor, DSP)上运行的传统机器学习算法需要进行模拟信号与数字信号之间的转换,这会导致不可避免的延迟和大量的能耗。本文研究了一种基于决策树的新型模拟电路分类器,用于实时机器信号分类和故障诊断。与传统的基于DSP的故障诊断方法相比,所设计的基于模拟电路的方法具有低延迟响应时间和低计算成本等优势。决策树算法的主要思想被整合到由滤波矩阵和幅值量化模块组成的模拟电路设计中。经过条件判断,传感器输出的模拟信号流经不同的分支通路。这种设计将传统决策树的分层判断机制创新性地转化为模拟电路实现,其中滤波器阵列对应于决策树的子节点,幅值量化模块实现信号分支选择。最终的输出结果是相应叶节点的分类结果,并且这些输出是根据要分类的信号类型进行编码。这种全时域模拟信号处理方法避免了模拟到数字的转换延迟,保持了决策树的可解释性,并且具有低延迟、易于实现和可扩展性的优点,在实时智能机器信号处理、状态监测和故障诊断等领域展现出良好的应用潜力。关键词:模拟电路分类器、决策树、机器信号处理、故障诊断、滤波矩阵1 引言近年来,人工智能(AI)技术发展迅速,其在特征提取和模式识别方面的强大能力在多个领域展现出卓越的表现。借助大规模数据集、先进的算法架构和高性能计算硬件,深度学习不断突破传统算法的局限,成为科学研究和工程创新的重要推动力。然而,这些模型的快速扩展和复杂性增加导致了显著的能耗问题,使计算效率成为现代人工智能系统中的主要瓶颈。目前,深度神经网络的训练和推理所需的计算需求和能耗远远超过了传统电子硬件的能力。这种资源密集型特性不仅限制了模型的可扩展性,也对可持续目标构成挑战。为解决这一问题,研究人员开始探索物理神经网络(PNNs)作为替代方案。这些日益严峻的挑战重新点燃了对替代硬件平台的兴趣。基于这些非常规计算平台的人工智能硬件统称为PNNs。这些硬件系统利用物理过程,如光学、机械或电子工艺,来执行计算任务,实现更高效、更低功耗的机器学习。PNNs在推理阶段展现出显著优势,充分利用物理系统的固有特性。例如,光学神经网络利用光学器件执行高效的矩阵运算,而机械神经网络则利用物理形变进行计算。这种基于物理过程的计算范式在降低能耗和提升处理速度方面具有巨大潜力。尽管PNN在信号推断方面具有巨大潜力,但在硬件中直接实现高效且稳定的训练过程仍是一大挑战。当前主流的物理硬件训练方法主要依赖反向传播算法或其变体。该算法需要高度精确的数学建模和数值求解过程,但由于物理系统中常见的非线性、噪声和动态特性,这些过程较为复杂,不仅大幅增加了整体计算时间和能耗,还降低了物理硬件在推理阶段的高效率。这些挑战严重限制了传统培训方法在工业环境中的实际应用。因此,提出一种能够克服传统训练方法局限、进一步探索实体硬件潜力的新方法至关重要。在此背景下,模拟电路的特性为实现机器学习提供了新的视角。与传统机器学习算法不同,受决策树启发的模拟电路系统可以直接利用电压和电流等物理量来实现分类和决策的能力。该方法不仅将机器学习的核心概念嵌入电路设计中,还实现了硬件层面的特征提取和分类任务。通过这种方法,可以实现振动信号的实时处理,这对于监测和诊断工业关键部件(如轴承故障)尤为重要。通过利用模拟电路的快速响应能力,振动信号可以高效地提取特征并准确分类,从而实现对不同类型轴承故障的实时监测和诊断。从另一个角度看,电机被广泛应用于工业自动化、精密制造和航空航天行业,而轴承是电机中最重要的部件之一。然而,在长时间运行过程中,电机轴承容易因负载冲击、磨损不均匀和润滑不足等因素而出现故障。这些故障不仅会导致异常振动和电机效率下降,还可能导致设备故障,甚至严重的安全事故和经济损失。因此,实现电机轴承故障的高效、精确实时监测和诊断已成为确保设备稳定运行和延长使用寿命的关键技术。目前,传统的轴承故障诊断方法主要依赖数字信号处理技术和复杂的算法模型。尽管这些方法能够实现较高的故障诊断准确性,但它们常常存在实时性能较差、高能耗和系统复杂度较高的问题。这些限制不仅增加了实施成本,也限制了其在对时间延迟敏感的工业环境中的适用性。此外,随着监控系统的规模和复杂度不断增长,传统方法相关的计算负担和延迟进一步加剧,使得实现高效且实时的故障诊断变得越来越困难。相比之下,基于模拟电路的故障诊断技术在响应速度快、功耗低和实施成本低方面具有显著优势,并有效避免了传统方法的延迟和高能耗问题。尤其在嵌入式和低功耗场景下,模拟电路故障诊断技术展现出极大的应用潜力。研究空白、挑战及主要贡献总结如下:(1)传统故障诊断依赖于模拟信号到数字信号转换,该转换存在固有延迟和高能耗,阻碍满足嵌入式低功耗的需求。受决策树启发的模拟电路分类器完全避免了数字转换,实现了实时信号处理,同时大幅降低了延迟和功耗。(2)鉴于对模拟电路分类的研究不足,该模拟电路设计了滤波矩阵和均方根值计算模块。重现了传统决策树的节点划分逻辑,分类方法完全基于模拟电路进行信号分类。 所提出的决策树启发模拟电路分类器可以直接处理传感器采集的原始连续振动信号,这在很大程度上满足了工业场景中低功耗嵌入式系统的需求。研究的其余部分如下。第二节介绍基于模拟电路的决策树分类器原理及模拟电路各模块的功能。第三部分讨论了模拟电路的拓扑图、数据集和实验验证结果。第四节介绍了用于验证实时信号分类实验的实际实验平台。第五节讨论了该方法的可能改进,第六部分总结了全文。2 基于模拟电路的决策树分类器设计2.1. 使用模拟电路实现的决策树算法决策树被应用于各种场景,尤其是在分类任务中。该算法通常使用某些标准,如信息增益、增益比或基尼指数,来为每个分割选择最佳属性。例如,利用信息增益,通过选择信息增益最大的属性进行数据划分,可以逐渐降低数据中的无序程度,从而实现分类预测,如下式: 其中, 表示数据集D的子集,被属性A的值v划分,Entropy (D)是数据集D的熵(熵越高,数据的不确定性越大;熵越低,数据的纯度越高),如上式所示, 是数据集D中属于 类数据的比例。与线性判别分析(LDA)分类器和基于规则的电路分类器的模拟电路实现相比,受决策树算法启发的模拟电路分类器具有明显优势。LDA电路在实际中需要高度精确的权重值才能实现正确的分类,这使得确保分类边界的稳定性变得困难。而基于规则的电路则随着规则数量的增加,电路复杂度和调试难度显著增加,且缺乏对信号变化的适应能力。相比之下,受决策树启发的模拟电路分类器执行分层阈值比较,这些比较可以自然地在硬件中实现为多级比较器结构,从而实现可解释性和稳定性之间的平衡。树状结构直观地展示了整个决策过程,使用户能够轻松理解每个决策的基础。具体来说,本研究设计的电路采用了受决策树启发的架构,滤波矩阵作为电路每个子模块中的决策节点。从每个子模输出信号计算出的均方根值对应于特定的分支准则,类似于传统决策树中的基于阈值的决策。实现受决策树启发的模拟电路系统的关键在于选择合适的模拟组件来替代决策树的子节点,并确定分支路径。如图1所示,设计了一个受决策树启发的模拟电路分类器的示意图。在这些组件中,由于滤波器的特性,同一信号通过具有相同参数的滤波器时,不同工作模式下的输出会有所变化。例如,信号经过低通和高通滤波模式时,输出会有所不同。因此,滤波矩阵电路取代了传统决策树中的子节点。然而,与传统的决策树不同,传统决策树是数字化实现并处理离散数据的,模拟电路处理的是连续信号。在传统决策树中,最佳分段阈值是基于信息增益选择的,而在模拟电路中,通过滤波电路提取物理量之间的关系,如均方根值,以实现信号分类。最终,电路根据信号的有效值关系对信号进行分类。如图1所示,不同类型的故障对应不同的RMS输出,从而实现分类。图1 受决策树启发的模拟电路分类器示意图2.2 模拟电路分类器的设计该电路使用Butterworth滤波器(MAX261, Maxim Inc.)作为滤波电路模块。MAX261 是由Maxim推出的CMOS双二阶通用开关电容有源滤波器,包含两个二阶开关电容有源滤波器和一个自由运算放大器。每个滤波单元由两个级联积分器和一个加法放大器组成,如图2所示,其中HP、LP和BP分别代表高通滤波器、低通滤波器和带通滤波器的输出。SCN指的是开关电容网络。通过反馈控制,每个滤波器的中心频率f0和品质因数Q被调整。其特性与连续时间滤波器极为接近,最大工作频率范围为0.4 Hz至57 MHz。为了实现理想的二阶状态变量响应,必须确保输入时钟频率与滤波器中心频率比(fCKL/f0)较大。此外,滤波器功能可由微处理器精确控制,带通、低通、高通、带阻和全通滤波器可在无外部元件的情况下构建。滤波器内部包含两个二阶滤波单元,每个滤波器的中心频率、Q值和工作模式可编程设置。该电路采用HP和LP滤波模式。为了控制所有滤波器,该设计在模式3使用MAX261实现LP和HP滤波模式。为了确定滤波器的最优参数,采用基于遗传算法(GA)的优化策略。该算法用于自动搜索能够最小化滤波器输出RMS差值的截止频率。这两种不同滤波模式的模型可以用二阶传递函数来描述。例如,LP滤波模式和HP滤波模式的传递函数G(s)如下所示: 其中, 是直流处的LP输出增益, 是HP增益, 近似为 。图2 二阶滤波器的内部框图该模拟电路的另一个重要组成部分是RMS计算模块(AD637,Analog Devices, Inc.)。RMS计算模块是一个专门为高精度RMS测量设计的硬件模块。当供电电压为±12V至±15V时,其测量范围可达0 Vrms至7 Vrms,输出电压纹波小于2mV。AD637模块采用了RMS方程的隐式解,克服了简单RMS计算的固有局限。图3(a)展示了AD637模块的功能实现图,图3(b)展示了AD637模块及其外设电路的示意图。AD637模块分为四个主要部分:绝对值电路、平方电路、滤波电路和缓冲放大器电路。AD637模块的输入电压可以是交流或直流信号,通过有源整流器A1和A2将其转换为单级电流,驱动平方电路。平方电路的输出电流驱动A4,A4与外部平均电容协同形成LP滤波器,而A3则利用该滤波器的放大器输出驱动最终的均方根计算。此外,为了获得更准确的有效值值,AD637模块必须包含两个运算放大器模块:(AD825, Analog Devices, Inc.),作为信号输入端的运算放大器,(OP07CR,Texas Instruments, Inc.),作为输出端的运算放大器,分别负责放大输入和输出信号。图3(a) AD637模块的功能实现图,以及(b)AD637模块及其外围电路的示意图由于AD637模块响应时间短,滤波后轴承振动信号的均方根值可以被准确测量。AD637模块实际执行的计算见下式: 其中, 是输入连续模拟信号, 是 在均方根计算后输出的连续模拟信号。除了上述两个模块外,该模拟电路还包括一个微控制器单元(MCU)和一个电源模块。MCU控制不同滤波器的工作模式,无论是LP滤波还是HP滤波器。电源模块为滤波器和AD637模块供电。原始模拟信号的输出经过不同滤波器后,作为AD637模块的输入。经过AD637模块处理后,滤波后的信号会从混沌信号转变为易于观测的均方根值信号作为输出。2.3 模拟电路分类器的参数优化在新型模拟电路参数优化中,首先描述参数优化问题。由于模拟电路使用滤波模块的LP和HP滤波模式,确定每种模式的参数成为关键问题。例如,一组包含两种轴承故障的模拟信号,当通过不同参数的滤波器时,会产生不同的输出。同样,这些信号在经过不同工作模式的滤波器后,其输出也会存在差异。然而,如果LP和HP滤波器的参数相同,则总会存在一个参数使得模拟信号经过LP和HP滤波器后的均方根值相等。将实际参数设置为该值,有助于区分各类轴承故障对应的模拟信号分别经过LP与HP滤波器后的均方根值,从而实现信号分类。因此,优化的目标函数表达为下式: 其中, 和 分别是截止频率为 的通过HP滤波器和LP滤波器输出信号S的有效值;S是一个包含n个数据点的信号段。因此,优化的目标函数是 的最小值。滤波器所需的参数是对应于 最接近零位置的频率值。在参数优化领域,存在多种智能算法,如GA、粒子群体优化、蜂群算法和神经网络算法。对于本文研究的单变量全局优化问题,选择GA作为求解工具,以获得全局最优解。GA是一种受自然选择和生物遗传过程启发的方向性、半随机搜索算法。在全能遗传学中,染色体代表具有独特遗传特征的个体,基因对应一组需要优化的特定参数。为了评估人群中的染色体,GA利用适应度函数为每条染色体分配一个分数,每一代中,优先选择得分较高的染色体进入下一代。经过多代进化,保留的高分染色体将代表优化问题的全局最优解或近似最优解。基于上述目标函数,算法中的适应度函数如下: 其中, 是信号段,即原始信号被划分为一百个长度相等的信号段,每个信号段有2000个数据点,因此n的值为2000; 是对应不同信号段的截止频率,通过持续调整和优化;最后,得到最接近0的 值,然后对应的 是最佳截止频率。图4 参数优化和信号流程图这里用三类示例来说明整个参数优化过程。在本例中,使用了1000、1500和2000赫兹的正弦波以及随机噪声。这三种正弦波和噪声被组合起来作为输入。通过GA,进行第一层优化以获得第一个参数值,如图4所示。第一层优化的结果约为1750Hz。该参数随后用作巴特沃斯HP和LP滤波器的截止频率。当1000、1500和2000Hz的正弦波通过HP和LP滤波器时,2000Hz的正弦波产生的输出是HP滤波后的均方根值大于LP滤波后的均方根值。相反,1000Hz和1500Hz的正弦波表现出相反的结果,因此对2000Hz的正弦波实现了分类。同样,通过输入1000Hz、1500Hz正弦波和随机噪声,第二层优化利用GA得到第二个参数,如图4所示。第二层优化的结果约为1400Hz。该参数用作巴特沃斯HP和LP滤波器的截止频率,而滤波后1000和1500Hz正弦波的均方根值关系则相反;因此,三种不同的信号都成功被区分。为避免仅使用正弦信号所带来的潜在偏差,并进一步验证所提出方法在非平稳和非线性振动信号中的适用性,本研究又增加了对三组具有不同故障特征频率的模拟滚动轴承故障信号进行了测试。所采用的信号模型旨在代表具有局部故障的滚动轴承的典型动态行为,其中故障位置与滚动元件之间的周期性接触会产生瞬态脉冲激励,进而激发系统的结构共振响应。对应的数学表述见下式: 其中, 和 表示第k次脉冲激发的幅度和时间, 是系统的固有角频率, 是阻尼比,n(t)表示背景噪声。具体来说,如图5所示,信号1、信号2和信号3的故障特征频率分别设置为100Hz、400Hz和800Hz。系统截止频率会根据这三种信号重新优化,然后利用训练好的系统对非固定和非线性振动信号进行分类。图5右侧面板的结果表明,设计的基于模拟电路的分类方法仍能有效区分不同类型的轴承故障信号,从而验证系统处理不同类型信号的性能。处理实际振动信号的方法性能将在第三节介绍。图5 系统在处理模拟非固定和非线性振动信号时的性能验证在所提方法中,GA被用于根据不同机器信号的特性调整滤波器频率。由于故障信号的频谱分布在不同机器间存在差异,需要进行GA重新优化以确保分类准确性。GA优化的核心目标是确定滤波器模块的截止频率 ,而 的最优值直接依赖于故障信号的频谱分布。因此,当电路应用于新机器时,必须再次进行GA优化,以匹配该条件下故障信号的频谱特性,确保滤波模块能够有效区分具有不同特征的信号。3 不同数据集上的实验验证3.1 模拟电路分类器的拓扑结构该电路旨在通过新型模拟电路分类器,利用轴承振动的模拟信号对轴承故障类型进行分类。在此过程中,电路设计尤为重要。真实电路设计与仿真之间存在一定差异。二阶巴特沃斯滤波器在仿真实验中与实际电路中产生的效果不同。巴特沃斯二阶滤波器在实电路中的滤波效应未达到仿真滤波器的性能。为了更好地使仿真与实际电路实验的滤波效果保持一致,对实际的巴特沃斯二阶滤波器进行改进,通过级联的方式提高其阶数,级联后滤波器阶数由二阶提升至四阶。由于实际巴特沃斯滤波器对滤波后的输出信号存在放大作用,实际模拟电路中级联滤波器所产生的实验结果与仿真结果基本相似,仅幅值略有放大,其余效果基本保持一致。本文以三分类轴承故障类型为例,介绍新型模拟电路的整体工作流程。对轴承故障类型进行分类时,轴承信号需要经过两层滤波器矩阵,对轴承故障类型进行分类时,电路需要通过两层滤波矩阵,每个滤波矩阵由一个HP滤波器和一个LP滤波器组成。每个过滤器还需要级联处理。因此,每个滤波矩阵需要四个MAX261模块。每层的滤波器成对工作,配置相同。因此,当使用MCU控制滤波器工作模式时,除时钟引脚外的所有引脚都可以使用相同的引脚,而时钟引脚则必须分配到不同的引脚,以设定每个滤波矩阵的不同工作参数。当轴承信号通过滤波矩阵的第一层时,会产生两个输出:一个来自LP滤波器,另一个来自HP滤波器。这两个输出各自分成两条相同的路径:其中一路作为AD637模块的输入,另一路作为加法器的输入。由于经过LP滤波器后的信号会显著削弱高频特征,而经过HP滤波器后的信号则会明显削弱低频特征,因此,滤波矩阵第一层输出的信号必须通过加法器进行叠加,以恢复信号的全频率特性。恢复后的信号将作为滤波矩阵第二层的输入。经过滤波矩阵第二层滤波后的输出被送入AD637模块作为输入。此时,AD637模块的输出为一种近似直线的模拟信号,它表示信号更加直观、易于观察的均方根值。若LP滤波器的输出大于HP滤波器的输出,则编码为 [1, 0];反之,若HP滤波器的输出大于LP滤波器的输出,则编码为 [0, 1],如表1所示。其中, 表示经过LP滤波后由AD637模块输出的有效值,而 表示经过HP滤波后由AD637模块输出的有效值。表1 编码方法新型模拟电路的原理图如图6所示。模拟信号进入整体电路,依次经过第一层滤波矩阵、加法器电路、第二层滤波矩阵,最终,第一层滤波矩阵的两个输出以及第二层滤波矩阵的两个输出共同作为AD637模块的输入。AD637模块的四个输出在图中表现为四条近似直线的模拟信号,可以通过示波器进行观察,并根据其幅值之间的关系判断故障类型。图6 模拟电路分类器原理图:(a)原理图;(b)详细电路接线图6(b)展示了该新型模拟电路的详细接线图。按照前文给出的器件参数,并注意滤波器级联时的相关设计要求,即可依据图6(b)所示的接线图复现实验系统。3.2 无刷直流电机(BLDCM)轴承数据集(1) 实验设置及信号预处理在本小节中,使用BLDCM的数据集来验证所设计系统的性能。实验装置如图7所示。图7(a)展示了不同故障类型BLDCM的测试平台,包括健康状态。图7(b)展示了三种故障类型。BLDCM及其故障类型的详细参数分别列于表2和表3中。图7 实验设置:(a)电机轴承振动信号采集系统;(b)外圈故障、内圈故障和健康轴承的说明表2 BLDCM 参数表3 BLDCM故障类型在搭建实际电路之前,首先在计算机上进行了仿真实验。在仿真实验中,所使用的数据集为通过NI数据采集卡采集得到的BLDCM轴承数据集。该数据集包含300万个数据点,采样频率为20 kHz。在本次仿真实验中,选取200,000个数据点用于参数优化,因为仿真过程需要模拟真实的模拟信号。因此,这200,000个数据点必须是连续选取的。在优化之前,需要对数据进行归一化处理,如下式所示,以确保所有信号的数据范围均处于 [0, 1] 区间。随后,通过计算信号的平均值并去除直流分量,可以使信号更加平滑,从而有利于后续的滤波与分析。 将200,000个数据点划分为100个片段,每个片段包含2,000个数据点。信号经滤波器后,对每个片段的输出计算其均方根值,因此从这200,000个数据点中可以得到100个RMS值。对于具有相同故障类型的轴承信号,当其通过不同滤波器时,输出的RMS值会有所不同。通过分析这些RMS值之间的关系,可以实现不同信号的分类。(2) 性能验证在本实验中,选取BLDCM健康信号、2 mm内圈故障信号以及2 mm外圈故障信号进行分析。仿真实验首先将这三种信号混合为一个复合信号进行优化。通过优化得到的最优截止频率分别为:第一层滤波器2690 Hz,第二层滤波器2400 Hz。实际上,如果滤波器参数设置在最优参数附近的一定范围内,系统仍然可以完成分类任务,如补充材料[34]所示。在确定滤波器的最佳参数后,将三种不同信号分别输入滤波器进行处理。每种信号在经过LP滤波器和HP滤波器后,会得到不同的RMS值关系。图8(a)展示了健康信号通过第一层LP和HP滤波器后的RMS值。其中HP滤波后的RMS值大于LP滤波后的RMS值,因此在编码中表示为 [0, 1]。图8(c)展示了2 mm内圈故障信号,其编码表示为 [1, 0]。 图8(e)展示了2 mm外圈故障信号,其编码同样表示为 [1, 0]。整体编码结果如表4所示。其中,只有健康信号的编码方式与另外两种信号不同。因此,在第一层滤波后即可成功识别健康信号。表4 不同信号经过全层滤波后的编码值对于2 mm内圈故障信号和2 mm外圈故障信号,在经过LP和HP滤波后,将两个输出信号进行相加。随后,将合成信号输入到经过优化的滤波器中以获得最佳参数。在为不同滤波器设置最佳参数后,将经过第一层滤波并相加后的内圈2 mm故障信号和外圈2 mm故障信号再次通过滤波器进行处理。该方法 会产生不同的RMS值关系,与第一层滤波的情况类似。图8(b)展示了健康信号通过第二层LP和HP滤波器后的RMS值,其编码为 [0, 1]。图8(d)展示了2 mm内圈故障信号,其编码同样为 [0, 1]。图8(f)展示了2 mm外圈故障信号经过第二层滤波后的结果,其编码为 [1, 0]。整体编码结果同样列于表4中。图8 (a)第一层和(b)第二层健康轴承信号的滤波矩阵输出;(c)第一层和(d)第二层的内圈故障信号;(e)第一层和(f)第二层外圈故障信号。(3) 抗噪声性能验证同时,为了验证系统在工业环境中的稳定性,在 BLDCM 健康信号、外圈故障信号和内圈故障信号的分类系统中加入了一定程度的合成混合噪声。考虑到实际应用场景中可能存在的噪声类型,混合噪声包括高斯白噪声、脉冲噪声和50 Hz工频噪声。其中,高斯白噪声占总能量的50%,而脉冲噪声和50 Hz工频噪声各占25%。通过改变混合噪声的信噪比(SNR),测试系统在干扰条件下的鲁棒性,实验结果如表5所示。实验结果表明,在加入噪声后系统性能有所下降,但在35dB噪声水平下系统仍然能够保持较高的分类准确率。此外,还进行了长期运行测试。系统连续运行24小时,并每4小时采集一次输出数据。结果表明,在整个测试过程中系统的编码输出基本保持稳定。表5 加入噪声后的实验结果(4) 分类精度的定量分析为了进一步定量分析系统在信号分类中的性能,引入了多种评价指标,其定义如下: 其中,TP、TN、FP和FN分别表示真正例、真反例、假正例和假反例。当系统正确识别健康状态信号,即两层滤波矩阵均满足 大于 时,该结果被计为真正例;当系统正确识别内圈故障信号时,则计为真反例。相应地,将健康信号误判为内圈故障计为假正例,将内圈故障误判为健康状态计为假反例。实验共测试400段健康信号和400段内圈故障信号,实验结果如表6和表7所示。结果表明,该系统能够较好地区分正样本与负样本,并具有较低的总体分类错误率。这说明该系统在当前分类任务中表现良好,具有较高的可靠性。表6 轴承状态预测结果统计表7 指标的结果(5) 信号持续时间的影响同时,考虑到分类准确率可能受到输入信号持续时间的影响,还进行了实验以验证实现准确分类所需的最短信号时长。实验中将信号持续时间设置为0–25ms,步长为2.5ms,实验结果如图9所示。结果表明,随着信号持续时间的增加,分类准确率逐渐提高。当信号持续时间达到10ms时,分类准确率达到87.50%;当信号持续时间增加到20ms时,分类准确率稳定达到99.50%。图9 不同持续时间输入信号的精度3.3 凯斯西储大学(CWRU)轴承数据集为了进一步验证所提出系统在不同数据集上的有效性,本研究使用公开的CWRU轴承数据集,在不同负载和故障严重程度条件下评估其多类别分类能力。本研究使用驱动端轴承振动数据,其故障类型、故障等级和负载情况列于表8中。其中IRF和ORF分别表示内圈故障和外圈故障。在本研究中,由于需要对4种信号进行分类,因此重新设置并优化系统参数,输出编码如表9所示。实验结果总结如图10所示。可以看出,所设计的系统能够成功分类4种不同工况下的轴承振动信号,结果验证了所提出的方法能够在不同数据集之间有效应用。表8 CWRU公共轴承数据集故障类型图10系统输出不同故障类型、故障大小、负载的CWRU信号表9 不同信号经过全层滤波后的编码值3.4 开关磁阻电机(SRM)数据集除了机械故障之外,电机还可能出现电气故障。为了进一步验证该方法在混合故障诊断中的性能,本研究还处理了来自SRM实验平台的另一组数据。SRM的故障类型如表10所示,其中混合故障为轴承机械故障、高电阻连接(HRC)电气故障以及霍尔传感器故障的组合。SRM的详细信息可参考文献[36]。在表10中,I1表示轴承内圈存在宽度为1 mm的缺陷,O1表示轴承外圈存在宽度为1 mm的缺陷,H表示霍尔故障,A1表示A相出现电阻值为1 Ω的HRC故障,B0.5表示B相出现0.5 Ω的HRC故障。表10 SRM的混合轴承和HRC故障分类结果如图11所示。经过第一层滤波矩阵后,系统将轴承内圈故障与HRC混合故障信号分离;随后第二层滤波矩阵进一步区分轴承外圈与HRC混合故障信号以及霍尔与HRC 混合故障信号。实验结果表明,即使在另一种类型电机的复杂混合故障信号条件下,所提出的系统仍然能够实现准确的故障分类。图11 SRM中混合轴承、HRC和霍尔故障的分类结果3.5 数据集的关键信息摘要上述实验结果表明,所设计系统能够对具有不同特征的多种机械信号进行分类。为了便于快速参考,数据集的关键信息总结于表11。表11 实验数据集关键信息4 电机故障信号实时分类4.1 实验装置实验装置如图12所示。从右到左,三个BLDCM分别代表不同故障状态:健康状态、内圈2 mm故障以及外圈2 mm故障。每个电机下方对应一个BLDCM驱动器用于驱动电机。电机左侧为直流稳压电源,为电机驱动器提供电源。电机右侧为实际设计的模拟电路模型,其中包括两层滤波矩阵、四AD637模块、一个加法器模块、一个主控制模块以及三个电源模块。电路板下方为IEPE信号调理放大器(LZ-1106E,Langzhen 公司),其为加速度传感器(CT1010L,Chengko 公司)提供恒流源驱动。图12 实时信号处理与故障诊断实验装置原始加速度传感器信号经过信号调理放大器后输出为未进行额外处理的归一化信号。加速度传感器通过磁铁固定在电机上方以获取轴承振动信号。信号调理放大器的输出连接至新型模拟电路的输入端。电机上方为示波器,实际模拟电路的四个输出连接至示波器的四个输入通道,通过观察示波器上的四个输出通道即可直接获得对应编码,从而判断相应的故障类型。4.2 实时信号处理的性能验证图 12 所示首先将加速度传感器放置在具有2 mm外圈故障的电机上,随后开启恒流源和电机电源,此时加速度传感器输出的电压信号表示电机在2 mm外圈故障条件下的轴承振动,该信号直接输入至设计电路板的输入端。电路板四个通道的输出从右到左依次对应为第一层滤波矩阵的LP输出、第一层滤波矩阵的HP输出、第二层滤波矩阵的LP输出以及第二层滤波矩阵的HP输出。这四个通道的输出连接至示波器对应的输入通道,结果如图13(a)所示。示波器的通道1和通道2分别表示第一层滤波矩阵的LP和HP输出,而通道3和通道4分别表示第二层滤波矩阵的LP和HP输出。由于示波器四个通道的零点均设置在同一位置,因此它们在屏幕上的相对位置代表其幅值大小。图13(a)显示通道1位于通道2上方,因此第一层滤波矩阵的输出编码为 [1, 0];同样,通道3位于通道4上方,因此第二层滤波矩阵的输出编码也为 [1, 0],因此整体输出编码为 [1, 0, 1, 0]。类似地,当将加速度传感器放置在具有2 mm内圈故障的电机上时,可得到图13(b)所示结果,其中第一层滤波矩阵输出编码为 [1, 0],第二层滤波矩阵输出编码为 [0, 1],整体输出编码为 [1, 0, 0, 1]。当加速度传感器放置在健康轴承电机上时,输出结果如图13(c)所示,其中第一层滤波矩阵输出编码为 [0, 1],第二层滤波矩阵输出编码为 [0, 1],因此整体输出编码为 [0, 1, 0, 1]。真实实验结果与仿真实验结果一致,验证了所提出方法的有效性。图13 实时信号处理和故障诊断的实验验证:(a)轴承外圈故障,(b)轴承内圈故障,(c)健康轴承4.3 变速条件下的性能验证本小节在电机转速变化条件下验证所提出系统的性能。如图14所示,图中左列展示系统能够成功分类信号的情况。图14(a)给出了电机转速变化曲线,其在1300 rpm到1600 rpm之间以5s为周期进行正弦变化。图14(b)和图14(c)分别展示第一层和第二层滤波器的RMS输出。图14 不同转速范围下的转速变化及相应的系统输出:(a)小转速范围条件下电机转速变化曲线;(b)和(c)成功分类的第1层和第2层滤波器输出;(d)大转速范围工况下电机转速变化曲线;(e)和(f)第1层和第2层未成功分类的滤波器输出在相同时间区间内,模拟电路系统输出的变化趋势与电机转速变化高度一致。在该转速范围内,经过短暂的瞬态调整后,系统能够输出正确编码,与之前在恒定转速条件下得到的结果一致。图14右列子图展示系统无法 正确分类信号的情况。图14(d)显示电机转速在1000 rpm到1600 rpm之间波动。在较高转速时系统能够保持正确输出编码,如图14(e)左侧粉色框所示;然而在最低转速处,如图14(e)红色框所示,输出编码出现偏差。这是因为转速变化会改变振动信号的频率分布,从而影响系统输出。因此,所提出系统更适用于恒速或小幅转速变化的运行工况。4.4 验证系统的响应时间和能耗本小节验证系统的推理时间和能耗。为了评估系统推理时间,在BLDCM的两种工况下进行实验,健康状态和霍尔传感器故障状态。霍尔传感器故障通过断开其信号线实现[37]。如图15左侧蓝色背景区域所示,系统输出对应电机健康状态,其编码为 [0, 1, 0, 1]。随后引入突发霍尔传感器故障,系统对该变化作出响应,如黄色背景区域所示。图15右侧粉色背景区域表示系统稳定检测到霍尔故障后的输出,其编码为 [0, 1, 1, 0]。通过测量黄色 区域宽度估算突发故障的检测时间约为78 ms,结果表明系统能够以较小延迟高效检测故障。图15 模拟电路分类系统对突发故障检测的响应时间此外,将所设计模拟系统的推理时间与在Raspberry Pi 嵌入式系统上实现的两种经典方法SVM和CNN 进行比较。所有实验均使用相同的BLDCM 信号。如表12和图16所示,尽管所提出系统的分类准确率略低于SVM和CNN,但其推理时间显著低于对比方法。这些结果表明,所提出系统在保持较高分类准确率的同时,在实时性能和计算延迟方面具有明显优势。图16 不同分类器在树莓派嵌入式系统上的推理结果表12 不同分类方法的准确率和推理时间的比较此外,通过对所有电路组件能耗进行求和评估系统功耗。在给定电源条件下,MCU、MAX261芯片以及RMS计算电路的功耗分别约为 89 mW、200 mW 和 120 mW,总系统功耗约为2.1 W。相比之下,使用Raspberry Pi 4B进行基于深度学习的故障信号训练与分类时功耗约为7 W,而Jetson Orin NX平台约为15 W。详细对比结果见表13。未来工作中,可使用电阻、电容和电感等基本电子元件构建固定截止频率滤波器,以进一步降低系统整体功耗。表13 功耗对比4.5 基于模拟电路的分类器的比较基于模拟电路的分类器近年来受到研究人员广泛关注,因为其无需ADC即可实现信号分类并具有较高处理速度。在本小节中,将所设计的决策树分类器与其他类型的模拟电路分类器进行比较,包括LDA、SVM和神经网络,结果如表14所示。LDA、SVM 以及神经形态分类器的模拟实现具有低功耗和高并行处理能力等显著优势,但这些方法通常需要定制的大规模集成电路(VLSI)技术,目前仍处于实验室研究阶段,在实现大规模生产之前VLSI成本极高。相比之下,本文所设计系统采用基本离散电子元件构建,具有良好的灵活性。当然,如果需要进行大规模生产,本文提出的方法和系统同样可以采用VLSI技术实现。此外,本研究从信号采集、信号处理到结果展示的完整信号流程角度验证了所设计系统的工程实用性。在参考文献39–43的对比方法中,仅使用离线数据集进行分类验证。因此,所设计系统更适用于部署成本低、系统可靠性高且工程实用性强的工业应用场景。表14 基于模拟电路的不同分类器的比较5 讨论本研究提出了一种基于模拟电路的实时信号处理与机器故障诊断系统。尽管所提出方法的性能已从多个方面得到验证,但类似主题的研究仍然较少。因此,为了使该方法在实际应用中更加灵活,本文对所提出方法的改进方向进行了讨论。5.1 信号类型数量与滤波器数量之间的关系本文提出了一种基于模拟电路的决策树分类器,该分类器利用多级滤波矩阵实现特征提取与分类。实验结果表明,对于本文所涉及的三分类轴承故障任务,需要两个内部决策节点,每个节点对应一层滤波矩阵,即HP和LP滤波模块。每个模块由四个滤波模块组成,因此完成三分类任务需要八个滤波模块。进一步分析表明,分类数量与滤波器数量之间呈线性增长关系,这一关系来源于决策树的二叉树结构,其具体表达如下式: 其中,K为待分类故障个数,L为所需滤波器个数。 此外,在当前系统实现中,每一层滤波矩阵均配置为双输出结构,即LP和HP滤波器。对于三分类任务,需要两层滤波矩阵,总计涉及八个MAX261模块,这在一定程度上引入了冗余。为了解决这一问题,未来可以研究一种结合分层复用策略和差分放大机制的协同方案。在优化各层滤波器的截止频率之后,可利用MCU精确控制时间间隔并动态调节MAX261器件的截止频率,使算法优化能够确定每一层的最佳截止频率。在这种架构下,仅需要一对LP和HP滤波器即可完成分类任务。5.2 集成电路实时故障诊断和显示当前电路设计依赖离散模块,包括1个MCU、8个模拟滤波器、4个RMS计算模块以及1个求和电路,总系统功耗约为2.1 W。未来可以利用专用集成电路(ASIC)技术将整个系统集成到单一芯片中。在ASIC实现中,这些模拟功能可以通过低功耗CMOS跨导单元和乘法单元实现,其功耗有望降低至约50 mW。此外,大规模集成不仅可以降低单芯片成本,同时还能提高系统可靠性,使系统更加适用于工业实时故障诊断应用。图17 未来ASIC技术后的模拟电路模型示意图最后,采用ASIC技术后的模拟电路模型结构如图17所示。传感器采集的信号可以直接连接到电路板预设的信号输入接口,这些信号随后由集成多种功能模块的ASIC进行处理,处理后的输出信号直接连接到电路板上方的LED灯组。通过引入基于差分放大的编码方案,输出可以表示为数字电压电平的形式:当输出电压为正时编码为 1,反之编码为 0。因此可以利用LED来识别编码输出,当编码输出为1时LED点亮,当编码输出为0时LED熄灭。通过观察LED状态即可直接识别故障类型,无需额外的显示设备如示波器。在这种情况下,可以实现设备状态信号的实时处理并降低功耗,从而减少能耗、延长设备运行时间并降低运行成本。5.3 滤波器特性对系统输出的影响本研究中使用四阶集成模拟滤波器构建决策树节点。滤波器特性如滤波类型、阶数和截止频率都会影响系统输出。为了进一步评估滤波器特性的影响,采用由离散电子元件构建的二阶滤波器替代原有四阶滤波器进行对比。所构建二阶LP滤波器的原理图和实物照片分别如图18(a)和图18(b)所示。该滤波器由电阻、电容和运算放大器构成,其详细信息见图18(a)。图18 由分立电子元件构成的二阶低电压滤波器:(a)原理图,(b)照片具有相同截止频率2670 Hz的二阶和四阶滤波器的频率响应如图19所示。可以看出,两种滤波器在衰减比、过渡带斜率以及整体频率响应形状方面均存在差异,这进一步导致相位失真和信号延迟的不同。图19 二阶和四阶LP滤波器的频率响应一段振动信号被输入到这两个LP滤波器中,其输出如图20所示。为了评估信号通过滤波器后的波形变化,在输入信号波形中标记了一个参考点。如图20(b)和图20(c)所示,随着滤波器阶数增加,系统群延迟表现出可预测的累积特性,这与巴特沃斯无限冲激响应(IIR)滤波器的理论相位响应特性一致。二阶和四阶滤波器的时间延迟分别为0.14 ms和0.32 ms。图20 (a)输入信号、(b)二阶和(c)四阶LP滤波器输出信号的波形比较将原有四阶LP滤波器替换为二阶LP滤波器后的最终编码输出如图21所示。与四阶滤波器输出相比,二阶滤波器输出(红线)中可以观察到明显的波形波动。然而第一层LP滤波器的幅值仍然高于第一层HP滤波器,因此仍可得到正确编码 [1, 0]。由此可见,不同滤波器类型和结构不会影响本研究中的正确输出,因此为了简化系统,仅对截止频率进行优化以构建基于决策树的模拟电路分类器。此外,与集成四阶滤波器(MAX261)相比,由离散电子元件构建的二阶滤波器将功耗从200 mW降低到21 mW。从另一方面来看,尽管四阶滤波器结构具有更高功耗和更大相位延迟,但其输出信号质量和稳定性更好,这在当前应用场景中具有重要的工程意义。总体而言,所设计系统具有较高鲁棒性,在考虑复杂度、噪声条件和成本等因素时,其滤波器选择具有较好的灵活性,适用于不同潜在应用场景。图21 二阶和四阶滤波器的输出结果5.4 噪声干扰在实际场景中的影响在工业机电系统中,信号采集和处理性能容易受到多种干扰源的影响。为了增强噪声测试的工程代表性,第3.2节采用由高斯白噪声、脉冲噪声和50 Hz工频干扰组成的混合噪声模型。具体而言,高斯白噪声用于模拟传感器热噪声和模拟前端电路中的电子噪声,脉冲噪声用于模拟由机械冲击、开关动作或强电磁干扰引起的瞬态扰动,而50 Hz工频噪声则反映工业环境中常见的干扰源,例如电源耦合、接地不良以及线路感应干扰[45–47]。这些噪声成分的组合旨在模拟连续干扰和脉冲干扰同时存在的复杂工业噪声环境。需要指出的是,尽管该混合噪声模型在一定程度上捕捉了工业环境中主要的干扰机制,但模拟噪声仍无法完全再现真实工业场景中复杂的电磁耦合条件。为了进一步验证所提出系统在实际工业干扰环境下的鲁棒性,本文在真实电磁干扰环境中进行了补充实验,如图22所示。在这些实验中,利用电磁线圈产生的交变磁场和永磁体产生的静态磁场作为干扰源,通过在电路板上方均匀移动干扰源来模拟工业环境中常见的磁场扰动和电磁耦合效应。图22 静磁场和交变磁场干扰下的系统鲁棒性实验两种磁场干扰条件下的系统输出与无干扰条件下的结果进行对比,如图23所示。信号输出波形可分为三个阶段:第一段、第二段和第三段分别对应电磁线圈干扰、无干扰以及永磁体干扰条件下的系统输出。实验结果表明,在外部磁场干扰作用下,系统输出仅出现轻微波动,但从输出编码角度来看,分类结果未受到影响。这进一步验证了所提出基于模拟电路的分类系统在复杂噪声干扰条件下的鲁棒性。图23 磁场干扰下的系统输出6 结论本研究提出了一种受决策树启发的模拟电路分类器,该分类器能够直接处理由加速度传感器采集的连续振动信号,而无需使用传统的ADC和DSP。该分类器由离散电子元件构建,并根据待处理信号对系统参数进行优化。通过仿真信号以及来自三种电机实验平台的真实振动信号验证了所设计系统的有效性和优势。与传统DSP方法相比,该方法避免了数字信号转换与预处理所固有的延迟和能耗,因此特别适用于低功耗和实时嵌入式应用场景。需要强调的是,所提出的方法并非旨在替代当前最先进的高精度模式分类方法,如新兴的深度学习模型。该方法更适用于信号类别数量有限且需要实时处理的应用场景。未来的研究工作将重点扩展模拟特征提取能力,并提高系统在更广泛信号类型和运行条件下的适应能力。编辑:曹希铭校核:李正平、陈凯歌、赵栓栓、赵学功、白亮、任超、Tina、陈宇航、海洋、Kira、王金、赵诚、肖鑫鑫、张优该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈