CN112836346B - 基于cn和pca的电机故障诊断方法、电子设备及介质 - Google Patents
基于cn和pca的电机故障诊断方法、电子设备及介质 Download PDFInfo
- Publication number
- CN112836346B CN112836346B CN202110020431.4A CN202110020431A CN112836346B CN 112836346 B CN112836346 B CN 112836346B CN 202110020431 A CN202110020431 A CN 202110020431A CN 112836346 B CN112836346 B CN 112836346B
- Authority
- CN
- China
- Prior art keywords
- data
- original data
- channel
- bit binary
- binary data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000003745 diagnosis Methods 0.000 title claims abstract description 56
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000013507 mapping Methods 0.000 claims abstract description 31
- 238000012545 processing Methods 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 20
- 230000009467 reduction Effects 0.000 claims abstract description 7
- 238000004590 computer program Methods 0.000 claims description 13
- 238000010606 normalization Methods 0.000 claims description 11
- 238000012549 training Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 6
- 238000000513 principal component analysis Methods 0.000 description 34
- 230000008569 process Effects 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000013527 convolutional neural network Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000002759 z-score normalization Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 229910002056 binary alloy Inorganic materials 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000009347 mechanical transmission Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/34—Testing dynamo-electric machines
- G01R31/343—Testing dynamo-electric machines in operation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Computational Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Control Of Electric Motors In General (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明实施例涉及故障检测技术领域,公开了基于CN和PCA的电机故障诊断方法、电子设备及介质。该方法包括:接收传感器采集的电机的原始数据;对原始数据进行归一化处理,并转换为15位的二进制数据;利用CN算法对二进制数据进行空间映射,得到10维的空间映射矩阵;利用PCA算法对空间映射矩阵进行降维处理,获取3维的特征矩阵;将特征矩阵输入预先训练的网络模型,输出诊断结果。实施本发明实施例,可以提升故障类型的空间可分离性,提升故障诊断算法的检测效率。
Description
技术领域
本发明涉及故障检测技术领域,具体涉及基于CN和PCA的电机故障诊断 方法、电子设备及介质。
背景技术
先进制造业是创新驱动发展和经济社会高质量发展的主引擎和主推手。高 端装备与智能制造领域中,电动机可直接将电能转化为直线运动的机械能,由 于其推力大、力密度高、长行程、低惯量、动态响应快、机械结构简单等优点。
电机直接驱动运动设备,省去了机械传动机构,完全消除机械传动元件的 速度和加速度的物理极限限制,已广泛应用于往复式伺服系统、工业机器人和 高精度定位直接驱动系统。
故障诊断技术包括的内容,是指在设备运行状态或称工作状态下,通过各 种监测手段,判别其工作是否正常;如果不正常,经过分析与判断,指出发生 了什么故障,便于管理人员维修;或者在故障未发生之前,提出可能发生故障 的预报,便于管理人员尽早采取措施,避免发生故障,或避免发生重大故障, 因而造成停机停产,给工程带来重大经济损失。这就是故障诊断技术的任务,也是发展设备故障诊断技术的目的。
在大数据这一时代背景下,伴随而来的是极大的计算复杂度和建模复杂度, 数据驱动的智能故障诊断方法更为适用,表现在其对海量、多源、高维数据进 行统计分析和信息提取的直接性和有效性。该技术以采集到的不同来源、不同 类型的监测数据作为基底,利用各种数据挖掘技术获取其中隐含的有用信息, 表征系统运行的正常模式以及故障模式,进而达到检测与诊断的目的。
智能故障诊断方法的性能极大地依赖于所提取的特征质量,包括特征的实 时变化、阶段性变化、趋势性变化及故障模式等,即数据的表示学习是智能故 障诊断技术的核心。传统的特征表示学习方法存在以下问题:
(1)需要领域的先验信息、专业知识和深厚的数学基础才能设计出合适的特 征提取方法;
(2)提取的特征多是浅层特征,针对复杂分类问题其泛化能力受到一定制约;
(3)受限于机械系统的物理特性,组件或故障条件变化均可能显著改变特征 提取方法或其评估标准;
(4)特征提取依赖于原有特征和评估标准,对于新特征的挖掘具有一定局限 性。
针对生产过程和设备运行中的海量状态数据及监控变量,信息采集通常为 多维向量形式,例如:旋转设备运行过程中基座的振动信息、驱动端的振动信 息、设备运行过程中的电流、电压信息等。针对多通道的传感器信息,如何将 传感器采集到的故障信息融合为多通道的矩阵形式是实现基于卷积神经网络的 故障诊断的关键。基于矩阵形式输入的卷积神经网络故障诊断方法中,传统的故障数据拼接实现的矩阵形式样本的处理方法具有一定的随机性,必然要求故 障诊断的深度网络具有较强的特征提取鲁棒性和较强的分类泛化能力。较强的 特征提取鲁棒性和较强的泛化能力必然要求相对复杂的网络结构满足要求。复 杂的卷积网络必然要求具有足够的样本进行模型预训练、实时故障诊断过程中较为强大的硬件支持。
发明内容
针对所述缺陷,本发明实施例公开了一种基于CN和PCA的电机故障诊断 方法、电子设备及介质,应用CN(Color Names)多颜色空间模型方法将原始 数据映射到高维空间,将高维空间映射结果通过PCA的方法进行降维,最终获 得故障数据高维空间主元,提升故障类型的空间可分离性,提升故障诊断算法 的检测效率。
本发明实施例第一方面公开一种基于CN和PCA的电机故障诊断方法,所 述方法包括:
接收传感器采集的电机的原始数据;
对所述原始数据进行归一化处理,并转换为15位的二进制数据;
利用CN算法对所述二进制数据进行空间映射,得到10维的空间映射矩阵;
利用PCA算法对所述空间映射矩阵进行降维处理,获取3维的特征矩阵;
将所述特征矩阵输入预先训练的网络模型,输出诊断结果。
作为一种可选的实施方式,在本发明实施例第一方面中,接收传感器采集 的电机的原始数据,包括:
接收一个或多个类型的传感器采集的电机的单通道或多通道原始数据;
所述传感器为振动传感器、电压互感器、电流互感器、加速度传感器中的 任一种或多种。
作为一种可选的实施方式,在本发明实施例第一方面中,所述原始数据为 单通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述单通道原始数据规范化至0-32,即每个单通道原始数据均采用5位 二进制数据进行表示;
选择相邻的3个5位二进制数据依次连接,组成15位的二进制数据:
Si={Mi-1,Mi,Mi+1}
其中,Si为组成的第i个15位的二进制数据,Mi为单通道的第i个原始数 据规范化后的5位二进制数据;0≤i≤L,L为单通道原始数据的长度;当i=0时, 使用线性差值方式获取Mi-1,当i=L时,使用线性差值方式获取Mi+1。
作为一种可选的实施方式,在本发明实施例第一方面中,所述原始数据为 双通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述双通道原始数据分别规范化至0-127以及0-255,即其中一通道的原 始数据采用7位二进制数据进行表示,另一通道的原始数据采用8位二进制数 据进行表示;
选择其中一通道规范后的7位二进制数据与另一通道对应的规范后的8位 二进制数据连接,组成15位的二进制数据:
Si={Mi,Ni}
其中,Si为组成的第i个15位的二进制数据,Mi为其中一通道的第i个原 始数据经规范化后的7位二进制数据;Ni为另一通道的第i个原始数据经规范 化后的8位二进制数据;0≤i≤L,L为每个通道原始数据的长度。
作为一种可选的实施方式,在本发明实施例第一方面中,所述原始数据为 三通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述三通道原始数据均规范化至0-32,即三通道原始数据中每个通道的原始数据均采用5位二进制数据进行表示;
选择三个通道相对应的原始数据进行连接,组成15位的二进制数据:
Si={Mi,Ni,Oi}
其中,Si为组成的第i个15位的二进制数据,Mi、Ni、Oi分别为三个通 道的第i个原始数据经规范化后的5位二进制数据;0≤i≤L,L为每个通道原始 数据的长度。
作为一种可选的实施方式,在本发明实施例第一方面中,利用CN算法对所 述二进制数据进行空间映射,得到10维的空间映射矩阵,包括:
使用转换矩阵32768×10对所述15位的二进制数据组成的数据集进行高维 空间映射,得到映射后的10为空间映射矩阵L×10,其中,L为原始数据的长度, 也即是15位的二进制数据组成的数据集的个数。
作为一种可选的实施方式,在本发明实施例第一方面中,预先训练的网络 模型,包括:
获取多组不同故障类型的电机样本数据以及未故障的电机样本数据,构建 样本集;
创建网络初始模型,并使用样本集中的各个电机样本数据对所述网络初始 模型进行训练,得到训练后的网络模型。
作为一种可选的实施方式,在本发明实施例第一方面中,所述原始数据的 长度L满足:
其中,β为传感器的采样频率,n为电机转速。
本发明实施例第二方面公开一种电子设备,包括:存储有可执行程序代码 的存储器;与所述存储器耦合的处理器;所述处理器调用所述存储器中存储的 所述可执行程序代码,用于执行本发明实施例第一方面公开的一种基于CN和 PCA的电机故障诊断方法。
本发明实施例第三方面公开一种计算机可读存储介质,其存储计算机程序, 其中,所述计算机程序使得计算机执行本发明实施例第一方面公开的一种基于 CN和PCA的电机故障诊断方法。
本发明实施例第四方面公开一种计算机程序产品,当所述计算机程序产品 在计算机上运行时,使得所述计算机执行本发明实施例第一方面公开的一种基 于CN和PCA的电机故障诊断方法。
本发明实施例第五方面公开一种应用发布平台,所述应用发布平台用于发 布计算机程序产品,其中,当所述计算机程序产品在计算机上运行时,使得所 述计算机执行本发明实施例第一方面公开的一种基于CN和PCA的电机故障诊断方法。
与现有技术相比,本发明实施例具有以下有益效果:
基于CN和PCA的故障诊断方法具有较高的诊断精度,同时模型收敛较快, 能够提升不同类别数据之间的可分离性,进而提升算法的故障诊断效率。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要 使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一 些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还 可以根据这些附图获得其他的附图。
图1是本发明实施例公开的基于CN和PCA的电机故障诊断方法的流程示 意图;
图2是本发明实施例公开的基于CN和PCA的电机故障诊断装置的结构示 意图;
图3是本发明实施例公开的电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清 楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全 部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性 劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书中的术语“第一”、“第二”、“第 三”、“第四”等是用于区别不同的对象,而不是用于描述特定顺序。本发明实施 例的术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,示 例性地,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、 方法、产品或设备固有的其它步骤或单元。
一种基于CN和PCA的电机故障诊断方法、电子设备及介质,用于根据检 测数据判断电机是否故障以及故障类型,基于CN和PCA的故障诊断方法具有 较高的诊断精度,同时模型收敛较快,能够提升不同类别数据之间的可分离性, 进而提升算法的故障诊断效率,以下对其进行详细描述。
实施例一
请参阅图1,图1是本发明实施例公开的一种基于CN和PCA的电机故障 诊断方法的流程示意图。其中,本发明实施例所描述的方法的执行主体通过软/ 硬件搭建而成,可以使用计算机、服务器等具有处理和存储功能的设备,在数据量较小的情况下,也可以使用手机、平板电脑等。如图1所示,该基于CN和 PCA的电机故障诊断方法包括以下步骤:
110,接收传感器采集的电机的原始数据。
传感器可以采集电机的内部参数数据,例如电压、电流以及加速度等,也 可以采集电机在运行过程中基座、驱动端等的振动信息。传感器可以为一个或 多个,为多个时,优选采用不同类型的传感器,例如,使用振动传感器采集基 座的振动信息,使用电压互感器采集电机的电压信息等,当然,在一些场景中,也可以采用相同类型传感器采集电机不同位置的参数信息。
每一个传感器采集的电机的原始数据称为一个通道的原始数据。可以理解 的是:当使用一个传感器采集电机数据时,则可以称为单通道原始数据;当使 用二个传感器采集电机数据时,则可以称为双通道原始数据等等。
原始数据的长度(或称数量)应该满足:
其中,L为每个通道的原始数据的长度,β为传感器的采样频率,n为电机 转速。在本发明较佳的实施例中,如果使用多个传感器对电机的原始数据进行 采集,则该多个传感器的采样频率相同,且采集时间相等,则保证各个通道的 原始数据的长度之间的相等的。
120,对所述原始数据进行归一化处理,并转换为15位的二进制数据。
在CN变换中,因为CN变换的变换矩阵为32768×10,即215×10,因此, 将原始数据构建成15位的二进制数据,则一方面利用后续的变换,另一方面因 为数据并不是孤立存在的,其与前后相邻的一个或多个数据或者不同通道相同时间的数据之间必然存在一定的关联,如果可以将通道内或通道间的原始数据 进行关联,则可以使得诊断更准确。
而本步骤120则是用于上述的转换操作,对原始数据进行归一化处理,并 将通道内或通道间的原始数据进行关联。
具体地:
对于单通道原始数据而言,即通过单一传感器采集的电机的原始数据,则 可以使用以下方式实现通道内的原始数据的关联。
将所述单通道原始数据规范化至0-32,即每个单通道原始数据均可以采用5 位二进制数据进行表示;选择相邻的3个5位二进制数据依次连接,从而就可 以拼接组成15位的二进制数据:
Si={Mi-1,Mi,Mi+1}
其中,Si为组成的第i个15位的二进制数据,Mi为单通道的第i个原始数 据规范化后的5位二进制数据;0≤i≤L,L为单通道原始数据的长度;当i=0时, 使用线性差值方式获取M-1,当i=L时,使用线性差值方式获取ML+1。
例如,如果第i-1、i、i+1个单通道原始数据归一化处理得到的数据转换成 二进制分别为10001、01110、11010、则得到的第i个15位的二进制数据为 100010111011010。
可以理解的是:单通道原始数据经过归一化以及拼接后得到的15位的二进 制数据的数量与单通道原始数据的数量是相等的,均具有L个数据。归一化处 理可以使用Z-score标准化方法。
类似的方法,针对双通道原始数据,则可以使用以下方法实现通道间的原 始数据进行关联:
将所述双通道原始数据分别规范化至0-127以及0-255,即其中一通道的原 始数据采用7位二进制数据进行表示,另一通道的原始数据采用8位二进制数 据进行表示;
选择其中一通道规范后的7位二进制数据与另一通道对应的规范后的8位 二进制数据连接,组成15位的二进制数据:
Si={Mi,Ni}
其中,Si为组成的第i个15位的二进制数据,Mi为其中一通道的第i个原 始数据经规范化后的7位二进制数据;Ni为另一通道的第i个原始数据经规范 化后的8位二进制数据;0≤i≤L,L为每个通道原始数据的长度。
例如,如果第1通道原始数据的第i个原始数据和第2通道的第i个原始数 据经过归一化处理得到的数据转换成二进制分别为1000110、11011100、则得到的第i个15位的二进制数据为100011011011100。
同样可以理解的是:双通道原始数据经过归一化以及拼接后得到的15位的 二进制数据的数量与双通道每个通道的原始数据的数量是相等的,均具有L个 数据。归一化处理可以使用Z-score标准化方法。
同样地,对于三通道原始数据而言,则可以使用以下方式实现通道间的原 始数据的关联。
将所述三通道每个通道的原始数据均规范化至0-32,即三通道原始数据中 每个通道的原始数据均可以采用5位二进制数据进行表示;
选择三个通道相对应的原始数据进行连接,组成15位的二进制数据:
Si={Mi,Ni,Oi}
其中,Si为组成的第i个15位的二进制数据,Mi、Ni、Oi分别为三个通 道的第i个原始数据经规范化后的5位二进制数据;0≤i≤L,L为每个通道原始 数据的长度。
例如,如果第1通道原始数据的第i个原始数据、第2通道的第i个原始数 据以及第3通道的第i个原始数据经过归一化处理得到的数据转换成二进制分别为10001、01110、11010、则得到的第i个15位的二进制数据为100010111011010。
可以理解的是:三通道原始数据经过归一化以及拼接后得到的15位的二进 制数据的数量与单通道原始数据的数量是相等的,均具有L个数据。归一化处 理可以使用Z-score标准化方法。
对于四通道或以上的原始数据而言,可以采用类似于双通道或三通道原始 数据的处理方式,例如,对于四通道而言,则可以将各个通道的原始数据分别 规范至0-9、0-9、0-9、0-16,即分别可以采用三位二进制、三位二进制、三位 二进制、四位二进制表示,然后将相同位置的归一化数据进行拼接,得到15位 的二进制数据。五通道则均可以规范至0-9。
如果超过十五通道,则不再本发明的保护范围内。
130,利用CN算法对所述二进制数据进行空间映射,得到10维的空间映射 矩阵。
近年来,CN(Color Names颜色命名)在目标检测、图像识别和动作识别 领域中得到了广泛的应用。它们是语言颜色标签,根据Berlin和Kay的研究所 得出的结论,颜色命名包括11个基本组成:黑色、蓝色、棕色、灰色、绿色、 橙色、粉红色、紫色、红色、白色和黄色。颜色名称是人为分配的,用于在现 实世界中呈现颜色。基于这种颜色属性方法,本文Google图像搜索检索的图像学习的映射,原本的RGB图像经过映射后,可以转换为概率和为1的11维颜 色空间:
X0(i,j)=W2C32768×10·X1(i,j)
在计算机视觉领域,常规对象跟踪器中的灰度值始终需要在[-0.5,0.5][6]内 进行归一化。颜色名称通过一种技术进行规范化以获得更好的性能。该技术是 通过将颜色名称投影到此10维子空间的正交标准进行的归一化处理。因此,11 维空间可以减少到10维。同时,此投影可以使颜色特征具有中心化的特点。
将CN算法应用于本发明实施例时,计算获取的15位二进制数据对应的十 进制值,确定十进制值对应与CN转换矩阵32768×10中的索引值,进而构建出 与十进制值对应的高维空间映射矩阵,最终获得结果映射矩阵MR[0:L,0:10]。
140,利用PCA算法对所述空间映射矩阵进行降维处理,获取3维的特征 矩阵。
在大多数研究的研究领域,需要大量的数据来找寻它们之间的规则。毫无 疑问,大量数据将提供大量信息,并使其更易于分析。但是,许多变量之间可 能存在相关性,因此有必要丢弃无用的数据并减小数据量以减少变量进而减少 计算量。如果随机减少待分析的数据或其维度,将不可避免地导致有用信息的 丢失PCA则被用来解决上述问题。
PCA(Principle Component Analysis主元分析)的原理是将原始样本数据投 影到新空间中,也就是将一组矩阵映射到另一个坐标系。在新的空间或坐标中, 不需要所有的原始样本,而是仅需要与原始样本的最大线性独立组的特征值相 对应的空间坐标。特征值及其对应的特征向量的计算是PCA算法的关键部分, 它是指此处与原始数据相对应的协方差矩阵的特征值。
本文在CN运算后变为L×10的10维数据。如果要使用PCA减小尺寸,则 首先需要计算10×10的协方差矩阵。其次,通过计算得到协方差矩阵的特征值 及其对应的特征向量。如果前3个特征值已占所有特征值的99%以上,则仅提取与前9个特征值相对应的特征向量,所选择的特征向量组成10×3的变换矩阵。 最后,将L×10的经过CN投影的数据乘以10×3的变换矩阵,即可得到新特征 空间中原始样本数据的对应坐标,10维数据在不丢失有用信息的基础上被成功 地减少到了3维。
150,将所述特征矩阵输入预先训练的网络模型,输出诊断结果。
因为一般的网络模型是针对图像数据的处理,因此可以将特征矩阵转换为 三通道图像,示例性地,将3维特征矩阵每一维特征矩阵中的特征均归一化至 0-255,然后3维特征矩阵分别代表R、G、B三个通道的图像的像素值,从而构 建三通道的R、G、B像素点值矩阵,然后选取满足网络模型的预设数个像素点 例如100个像素点构建10*10的图像,从而得到三通道的R、G、B图像。
网络模型优选可以采用MobileNetV3 small的网络结构。
选取多个电机的样本数据组成的样本集对模型进行训练。优选地,样本数 据包括多种数据来源,例如包括多个单通道样本数据、多个双通道样本数据、 多个三通道样本数据、以及多个四通道样本数据等,而且,样本数据应当包括 多个不同故障类型例如失磁、气隙变化、定子绕组故障等样本数据,也应当包 括多个未故障时的电机样本数据。
样本数据集中每个样本数据的长度也应该满足:
对样本数据同样需要上述步骤120-140的处理过程,然后将得到的3维特征 矩阵转换成三通道图像,输入网络初始模型进行训练,训练的标签为样本的故 障类型,例如,将未故障样本数据的标签定义为0,将失磁故障的样本数据的标 签定义为1,将气隙变化故障的标签定义为2……,使用标签对网络初始模型的 输出结果进行反向传播,从而确定网络初始模型的各个参数的合适函数或值, 得到最终的网络模型。
将原始数据经过处理后的三通道图像输入训练好的网络模型,则可以得到 该原始数据对应的电机的状态,如果为故障状态,则还可以确定故障类型。
为了验证CN+PCA能够提升多类原始数据的空间可分离性,本发明实施例 应用CN+PCA的方法对西储大学的滚动轴承故障数据进行处理,构建训练和测 试样本集。同时依据西储大学原始数据直接构建训练样本和测试样本数据集。两种不同的数据集在轻量级的CNN网络(整体结构为五层)上进行仿真验证, 仿真实验的结果表明,基于CN+PCA故障数据处理的故障诊断方法具有较高的 诊断精度,同时模型收敛较快。识别精度和误差迭代关系如表1所示。西储数 据实验结果已经证明了应用CN实现的高维空间映射,应用PCA进行降维构建 样本数据的方法能够提升不同类别样本数据之间的可分离性。进而提升算法的故障诊断效率。
表1:基于CN和PCA的故障诊断训练、测试的损失和准确率
迭代次数 | 训练损失 | 训练精度 | 测试损失 | 测试精度 |
1 | 0.0217 | 99.5438 | 0.0004 | 100.000 |
2 | 0.0011 | 99.9983 | 0.0001 | 100.000 |
3 | 0.0007 | 99.9983 | 0.0006 | 100.000 |
4 | 0.0005 | 100.000 | 0.0000 | 100.000 |
5 | 0.0004 | 100.000 | 0.0000 | 100.000 |
6 | 0.0005 | 99.9983 | 0.0000 | 100.000 |
7 | 0.0003 | 99.9983 | 0.0000 | 100.000 |
8 | 0.0003 | 100.000 | 0.0000 | 100.000 |
9 | 0.0002 | 99.9983 | 0.0000 | 100.000 |
10 | 0.0002 | 100.000 | 0.0000 | 100.000 |
11 | 0.0001 | 100.000 | 0.0000 | 100.000 |
12 | 0.0001 | 100.000 | 0.0000 | 100.000 |
13 | 0.0001 | 100.000 | 0.0000 | 100.000 |
14 | 0.0001 | 100.000 | 0.0000 | 100.000 |
15 | 0.0001 | 100.000 | 0.0000 | 100.000 |
16 | 0.0001 | 100.000 | 0.0000 | 100.000 |
17 | 0.0001 | 100.000 | 0.0000 | 100.000 |
18 | 0.0001 | 100.000 | 0.0000 | 100.000 |
19 | 0.0001 | 100.000 | 0.0000 | 100.000 |
20 | 0.0001 | 100.000 | 0.0000 | 100.000 |
实施例二
请参阅图2,图2是本发明实施例公开的一种基于CN和PCA的电机故障 诊断装置的结构示意图。如图2所示,该装置可以包括:
接收单元210,用于接收传感器采集的电机的原始数据;
归一化单元220,用于对所述原始数据进行归一化处理,并转换为15位的 二进制数据;
映射单元230,用于利用CN算法对所述二进制数据进行空间映射,得到10 维的空间映射矩阵;
降维单元240,用于利用PCA算法对所述空间映射矩阵进行降维处理,获 取3维的特征矩阵;
诊断单元250,用于将所述特征矩阵输入预先训练的网络模型,输出诊断结 果。
作为一种可选的实施方式,接收传感器采集的电机的原始数据,包括:
接收一个或多个类型的传感器采集的电机的单通道或多通道原始数据;
所述传感器为振动传感器、电压互感器、电流互感器、加速度传感器中的 任一种或多种。
作为一种可选的实施方式,所述原始数据为单通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述单通道原始数据规范化至0-32,即每个单通道原始数据均采用5位 二进制数据进行表示;
选择相邻的3个5位二进制数据依次连接,组成15位的二进制数据:
Si={Mi-1,Mi,Mi+1}
其中,Si为组成的第i个15位的二进制数据,Mi为单通道的第i个原始数 据规范化后的5位二进制数据;0≤i≤L,L为单通道原始数据的长度;当i=0时, 使用线性差值方式获取Mi-1,当i=L时,使用线性差值方式获取Mi+1。
作为一种可选的实施方式,所述原始数据为双通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述双通道原始数据分别规范化至0-127以及0-255,即其中一通道的原 始数据采用7位二进制数据进行表示,另一通道的原始数据采用8位二进制数 据进行表示;
选择其中一通道规范后的7位二进制数据与另一通道对应的规范后的8位 二进制数据连接,组成15位的二进制数据:
Si={Mi,Ni}
其中,Si为组成的第i个15位的二进制数据,Mi为其中一通道的第i个原 始数据经规范化后的7位二进制数据;Ni为另一通道的第i个原始数据经规范 化后的8位二进制数据;0≤i≤L,L为每个通道原始数据的长度。
作为一种可选的实施方式,所述原始数据为三通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述三通道原始数据均规范化至0-32,即三通道原始数据中每个通道的 原始数据均采用5位二进制数据进行表示;
选择三个通道相对应的原始数据进行连接,组成15位的二进制数据:
Si={Mi,Ni,Oi}
其中,Si为组成的第i个15位的二进制数据,Mi、Ni、Oi分别为三个通 道的第i个原始数据经规范化后的5位二进制数据;0≤i≤L,L为每个通道原始 数据的长度。
作为一种可选的实施方式,利用CN算法对所述二进制数据进行空间映射, 得到10维的空间映射矩阵,包括:
使用转换矩阵32768×10对所述15位的二进制数据组成的数据集进行高维 空间映射,得到映射后的10为空间映射矩阵L×10,其中,L为原始数据的长度, 也即是15位的二进制数据组成的数据集的个数。
作为一种可选的实施方式,预先训练的网络模型,包括:
获取多组不同故障类型的电机样本数据以及未故障的电机样本数据,构建 样本集;
创建网络初始模型,并使用样本集中的各个电机样本数据对所述网络初始 模型进行训练,得到训练后的网络模型。
作为一种可选的实施方式,所述原始数据的长度L满足:
其中,β为传感器的采样频率,n为电机转速。
实施例三
请参阅图3,图3是本发明实施例公开的一种电子设备的结构示意图。如图 3所示,该电子设备可以包括:
存储有可执行程序代码的存储器310;
与存储器310耦合的处理器320;
其中,处理器320调用存储器310中存储的可执行程序代码,执行实施例 一中的基于CN和PCA的电机故障诊断方法中的部分或全部步骤。
本发明实施例公开一种计算机可读存储介质,其存储计算机程序,其中, 该计算机程序使得计算机执行实施例一中的基于CN和PCA的电机故障诊断方 法中的部分或全部步骤。
本发明实施例还公开一种计算机程序产品,其中,当计算机程序产品在计 算机上运行时,使得计算机执行实施例一中的基于CN和PCA的电机故障诊断 方法中的部分或全部步骤。
本发明实施例还公开一种应用发布平台,其中,应用发布平台用于发布计 算机程序产品,其中,当计算机程序产品在计算机上运行时,使得计算机执行实施例一中的基于CN和PCA的电机故障诊断方法中的部分或全部步骤。
在本发明的各种实施例中,应理解,所述各过程的序号的大小并不意味着 执行顺序的必然先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应 对本发明实施例的实施过程构成任何限定。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为 单元显示的部件可以是或者也可以不是物单元,即可位于一个地方,或者也可 以分布到多个网络单元上。可根据实际的需要选择其中的部分或全部单元来实 现本实施例方案的目的。
另外,在本发明各实施例中的各功能单元可以集成在一个处理单元中,也 可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。 所述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式 实现。
所述集成的单元若以软件功能单元的形式实现并作为独立的产品销售或使 用时,可以存储在一个计算机可获取的存储器中。基于这样的理解,本发明的 技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或者 部分,可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器 中,包括若干请求用以使得一台计算机设备(可以为个人计算机、服务器或者 网络设备等,具体可以是计算机设备中的处理器)执行本发明的各个实施例所述方法的部分或全部步骤。
在本发明所提供的实施例中,应理解,“与A对应的B”表示B与A相关联, 根据A可以确定B。但还应理解,根据A确定B并不意味着仅仅根据A确定B, 还可以根据A和/或其他信息确定B。
本领域普通技术人体可以理解所述实施例的各种方法中的部分或全部步骤 是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存 储介质中,存储介质包括只读存储器(Read–Only Memory,ROM)、随机存储 器(Random Access Memory,RAM)、可编程只读存储器(Programmable Read–only Memory,PROM)、可擦除可编程只读存储器(Erasable Programmable Read–Only Memory,EPROM)、一次可编程只读存储器(One–time Programmable Read–Only Memory,OTPROM)、电子抹除式可复写只读存储器(Electrically–Erasable Programmable Read–Only Memory,EEPROM)、只读光盘(CompactDisc Read–Only Memory,CD–ROM)或其他光盘存储器、磁盘存储器、磁带存储器、或者能够用于携带或存储数据的计算机可读的任何其他介质。
以上对本发明实施例公开的一种基于CN和PCA的电机故障诊断方法、电 子设备及介质进行了详细介绍,本文中应用了具体个例对本发明的原理及实施 方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心 思想;同时,对于本领域的一般技术人体,依据本发明的思想,在具体实施方 式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发 明的限制。
Claims (6)
1.一种基于CN和PCA的电机故障诊断方法,其特征在于,包括:
接收传感器采集的电机的原始数据;
对所述原始数据进行归一化处理,并转换为15位的二进制数据;
利用CN算法对所述二进制数据进行空间映射,得到10维的空间映射矩阵;
利用PCA算法对所述空间映射矩阵进行降维处理,获取3维的特征矩阵;
将所述特征矩阵输入预先训练的网络模型,输出诊断结果;
接收传感器采集的电机的原始数据,包括:
接收一个或多个传感器采集的电机的单通道或多通道原始数据;
所述传感器为振动传感器、电压互感器、电流互感器、加速度传感器中的任一种或多种;
当所述原始数据为单通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述单通道原始数据规范化至0-32,即每个单通道原始数据均采用5位二进制数据进行表示;
选择相邻的3个5位二进制数据依次连接,组成15位的二进制数据:
Si={Mi-1,Mi,Mi+1}
其中,Si为组成的第i个15位的二进制数据,Mi为单通道的第i个原始数据规范化后的5位二进制数据;0≤i≤L,L为单通道原始数据的长度;当i=0时,使用线性差值方式获取Mi-1,当i=L时,使用线性差值方式获取Mi+1;
当所述原始数据为双通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述双通道原始数据分别规范化至0-127以及0-255,即其中一通道的原始数据采用7位二进制数据进行表示,另一通道的原始数据采用8位二进制数据进行表示;
选择其中一通道规范后的7位二进制数据与另一通道对应的规范后的8位二进制数据连接,组成15位的二进制数据:
Si={Mi,Ni}
其中,Si为组成的第i个15位的二进制数据,Mi为其中一通道的第i个原始数据经规范化后的7位二进制数据;Ni为另一通道的第i个原始数据经规范化后的8位二进制数据;0≤i≤L,L为每个通道原始数据的长度;
当所述原始数据为三通道原始数据;
所述对所述原始数据进行归一化处理,并转换为15位的二进制数据,包括:
将所述三通道原始数据均规范化至0-32,即三通道原始数据中每个通道的原始数据均采用5位二进制数据进行表示;
选择三个通道相对应的原始数据进行连接,组成15位的二进制数据:
Si={Mi,Ni,Oi}
其中,Si为组成的第i个15位的二进制数据,Mi、Ni、Oi分别为三个通道的第i个原始数据经规范化后的5位二进制数据;0≤i≤L,L为每个通道原始数据的长度。
2.根据权利要求1所述的基于CN和PCA的电机故障诊断方法,其特征在于,利用CN算法对所述二进制数据进行空间映射,得到10维的空间映射矩阵,包括:
使用转换矩阵32768×10对所述15位的二进制数据组成的数据集进行高维空间映射,得到映射后的10为空间映射矩阵L×10,其中,L为原始数据的长度,也即是15位的二进制数据组成的数据集的个数。
3.根据权利要求1所述的基于CN和PCA的电机故障诊断方法,其特征在于,预先训练的网络模型,包括:
获取多组不同故障类型的电机样本数据以及未故障的电机样本数据,构建样本集;
创建网络初始模型,并使用样本集中的各个电机样本数据对所述网络初始模型进行训练,得到训练后的网络模型。
4.根据权利要求1所述的基于CN和PCA的电机故障诊断方法,其特征在于,所述原始数据的长度L满足:
其中,β为传感器的采样频率,n为电机转速。
5.一种电子设备,其特征在于,包括:存储有可执行程序代码的存储器;与所述存储器耦合的处理器;所述处理器调用所述存储器中存储的所述可执行程序代码,用于执行权利要求1至4任一项所述的一种基于CN和PCA的电机故障诊断方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储计算机程序,其中,所述计算机程序使得计算机执行权利要求1至4任一项所述的一种基于CN和PCA的电机故障诊断方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110020431.4A CN112836346B (zh) | 2021-01-07 | 2021-01-07 | 基于cn和pca的电机故障诊断方法、电子设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110020431.4A CN112836346B (zh) | 2021-01-07 | 2021-01-07 | 基于cn和pca的电机故障诊断方法、电子设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112836346A CN112836346A (zh) | 2021-05-25 |
CN112836346B true CN112836346B (zh) | 2024-04-23 |
Family
ID=75928342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110020431.4A Active CN112836346B (zh) | 2021-01-07 | 2021-01-07 | 基于cn和pca的电机故障诊断方法、电子设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112836346B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114236386B (zh) * | 2021-12-16 | 2023-10-31 | 河北工业大学 | 一种开关磁阻电机功率变换器故障诊断方法 |
CN117849617B (zh) * | 2024-01-10 | 2024-06-14 | 关兴 | 一种星形连接电机阻感故障检测方法、系统、设备及介质 |
CN117991093B (zh) * | 2024-04-03 | 2024-06-18 | 成都航天凯特机电科技有限公司 | 一种永磁同步电机故障诊断方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108994088A (zh) * | 2018-06-12 | 2018-12-14 | 武汉科技大学 | 基于pca降维与dbn网络的hagc缸故障诊断方法及装置 |
CN109272044A (zh) * | 2018-09-19 | 2019-01-25 | 郑州云海信息技术有限公司 | 一种图像相似度确定方法、装置、设备及存储介质 |
CN110895701A (zh) * | 2019-06-12 | 2020-03-20 | 河南理工大学 | 基于cn和fhog的森林火灾在线识别方法及装置 |
CN111191725A (zh) * | 2019-12-31 | 2020-05-22 | 北京市燃气集团有限责任公司 | 一种基于pca数据处理的rbf故障诊断方法、系统、终端及计算机存储介质 |
-
2021
- 2021-01-07 CN CN202110020431.4A patent/CN112836346B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108994088A (zh) * | 2018-06-12 | 2018-12-14 | 武汉科技大学 | 基于pca降维与dbn网络的hagc缸故障诊断方法及装置 |
CN109272044A (zh) * | 2018-09-19 | 2019-01-25 | 郑州云海信息技术有限公司 | 一种图像相似度确定方法、装置、设备及存储介质 |
CN110895701A (zh) * | 2019-06-12 | 2020-03-20 | 河南理工大学 | 基于cn和fhog的森林火灾在线识别方法及装置 |
CN111191725A (zh) * | 2019-12-31 | 2020-05-22 | 北京市燃气集团有限责任公司 | 一种基于pca数据处理的rbf故障诊断方法、系统、终端及计算机存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112836346A (zh) | 2021-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111721535B (zh) | 一种基于卷积多头自注意力机制的轴承故障检测方法 | |
Yao et al. | A survey of transfer learning for machinery diagnostics and prognostics | |
CN112836346B (zh) | 基于cn和pca的电机故障诊断方法、电子设备及介质 | |
US12228476B2 (en) | Fault signal locating and identifying method of industrial equipment based on microphone array | |
CN112200244B (zh) | 一种基于递阶对抗训练的航天发动机异常智能检测方法 | |
CN114359283B (zh) | 基于Transformer的缺陷检测方法和电子设备 | |
CN113780412A (zh) | 故障诊断模型训练方法、系统及其诊断方法、系统 | |
CN112487949B (zh) | 一种基于多模态数据融合的学习者行为识别方法 | |
CN110705525A (zh) | 一种诊断滚动轴承故障的方法及装置 | |
CN112257263B (zh) | 基于自注意力机制的设备剩余寿命预测系统 | |
CN116503399B (zh) | 基于yolo-afps的绝缘子污闪检测方法 | |
CN118090211A (zh) | 一种基于时频特征融合的电梯曳引机轴承故障诊断方法 | |
CN115859142A (zh) | 一种基于卷积变压器生成对抗网络的小样本滚动轴承故障诊断方法 | |
CN114067286A (zh) | 基于序列化可形变注意力机制的高位相机车辆重识别方法 | |
Bao et al. | A transformer model-based approach to bearing fault diagnosis | |
CN115392301A (zh) | 换流变压器状态识别方法、装置、设备、介质和程序产品 | |
CN113435321A (zh) | 一种主轴轴承状态评估方法、系统、设备及可读存储介质 | |
CN118779807A (zh) | 基于多模态声图融合的电力变压器缺陷诊断方法 | |
CN118918424A (zh) | 一种基于双通道多尺度注意力特征融合的轴承故障诊断方法 | |
CN116503398B (zh) | 绝缘子污闪检测方法、装置、电子设备及存储介质 | |
Xie et al. | Multi-scale deep neural network for fault diagnosis method of rotating machinery | |
CN112733706B (zh) | 基于双线性lbp的电机故障诊断方法、电子设备及介质 | |
CN116089863A (zh) | 一种基于尺度自适应度量的滚动轴承故障诊断方法 | |
CN117439800B (zh) | 一种网络安全态势预测方法、系统及设备 | |
CN118470468B (zh) | 一种射电天文信号的射频干扰检测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |