CN116206154B - 一种非常见工况下的污水处理过程的故障分类方法 - Google Patents
一种非常见工况下的污水处理过程的故障分类方法 Download PDFInfo
- Publication number
- CN116206154B CN116206154B CN202310130653.0A CN202310130653A CN116206154B CN 116206154 B CN116206154 B CN 116206154B CN 202310130653 A CN202310130653 A CN 202310130653A CN 116206154 B CN116206154 B CN 116206154B
- Authority
- CN
- China
- Prior art keywords
- data
- fault
- source domain
- target domain
- stationary
- 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 130
- 230000008569 process Effects 0.000 title claims abstract description 63
- 239000010865 sewage Substances 0.000 title claims abstract description 45
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 238000012706 support-vector machine Methods 0.000 claims abstract description 26
- 238000012549 training Methods 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 14
- 238000010219 correlation analysis Methods 0.000 claims abstract description 10
- 238000013508 migration Methods 0.000 claims abstract description 7
- 230000005012 migration Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 63
- 239000013598 vector Substances 0.000 claims description 54
- 230000006870 function Effects 0.000 claims description 34
- 238000001514 detection method Methods 0.000 claims description 15
- 238000013507 mapping Methods 0.000 claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 9
- 238000013112 stability test Methods 0.000 claims description 9
- 238000013145 classification model Methods 0.000 claims description 8
- 238000012544 monitoring process Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 230000009977 dual effect Effects 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 230000007774 longterm Effects 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 230000002195 synergetic effect Effects 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 239000000758 substrate Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 101000984533 Homo sapiens Ribosome biogenesis protein BMS1 homolog Proteins 0.000 description 3
- 102100027057 Ribosome biogenesis protein BMS1 homolog Human genes 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000006911 nucleation Effects 0.000 description 2
- 238000010899 nucleation Methods 0.000 description 2
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000013526 transfer learning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
- G06V10/765—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects using rules for classification or partitioning the feature space
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0224—Process history based detection method, e.g. whereby history implies the availability of large amounts of data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/761—Proximity, similarity or dissimilarity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Multimedia (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种非常见工况下的污水处理过程的故障分类方法,首先分离出非平稳变量的过程数据,基于邻域保持成分分析的方法提取非平稳特征;然后基于典型相关分析的迁移学习方法减小源域和目标域正常数据的分布差异,获得投影后的源域和目标域数据;画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图,对两者做相似度分析,实现对目标域的第一次故障分类;如果目标域单类未带标签的故障数据的贡献图与源域多类带标签的故障数据的贡献图的相似度相同,将对应的投影后的源域多类带标签的故障数据训练支持向量机模型;根据对应的投影后的目标域少量带标签的故障数据更新支持向量机模型,实现对目标域的第二次故障分类。
Description
技术领域
本发明涉及工业过程故障分类技术领域,尤其涉及一种非常见工况下的污水处理过程的故障分类方法。
背景技术
污水处理过程是一个具有多工况特点的复杂工业过程。增强故障检测的能力是保证污水处理过程安全稳定运行的关键。由于天气和季节的影响,污水处理过程一般具有非线性、多工况、多变量等特征。一般来说,污水处理过程的工况可以分为常见工况和非常见工况。常见工况的特点是带标签数据容易获取且有着大量的历史数据。非常见工况的特点是带标签数据不容易获取,而且只有少量的历史数据。非常见工况下的污水处理过程的故障分类模型不容易被建立,这是因为传统的数据驱动的故障分类方法如支持向量机、神经网络等方法均依赖于大量的历史数据。当历史数据不足时,就无法训练模型。此外,传统的数据驱动的故障分类模型会受到工况变化的影响,这是因为不同工况下提取到的数据特征是有差异的,而大多数数据驱动的故障分类方法没有考虑到这一点。除此之外,大多数故障分类方法如支持向量机、神经网络等都是先离线训练好分类模型,再进行故障分类。在数据量和故障类型比较多的情况下,这些方法训练模型所耗费的时间是比较长的。
目前,数据驱动的故障分类方法在污水处理过程的非常见工况的故障分类问题上的应用有限,这些方法大多数应用于历史数据较多的常见工况中。传统方法建立的故障分类模型会受到工况变化的影响,这是因为常见工况和非常见工况提取到的数据特征是有差异的,而大多数方法在特征提取时没有考虑到工况变化的影响。另外,现有的大部分故障分类算法没有考虑从故障检测模型的基础上去开发算法。
发明内容
为了解决污水处理过程的非常见工况下因历史数据不足无法建立故障分类模型的问题,本发明提供一种非常见工况下的污水处理过程的故障分类方法,针对多工况的污水处理过程,采用迁移学***稳特征从而降低工况变化对***模型的影响,采用基于邻域保持成分分析的特征提取方法,最后采用贡献图和支持向量机结合的方法实现对非常见工况下的故障分类。
为解决上述技术问题,本发明的实施例提供如下方案:
一种非常见工况下的污水处理过程的故障分类方法,包括以下步骤:
S1,采集污水处理过程数据,包括源域和目标域的历史数据,数据预处理后分离出非平稳变量的过程数据,基于邻域保持成分分析的方法提取非平稳特征;
S2,基于典型相关分析的迁移学习方法减小源域和目标域正常数据的分布差异,获得投影后的源域和目标域的数据;
S3,根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析,对目标域进行第一次故障分类;
S4,如果目标域单类未带标签的故障数据的贡献图与源域多类带标签的故障数据的贡献图的相似度相同,将对应的投影后的源域多类带标签的故障数据训练支持向量机模型;根据对应的投影后的目标域少量带标签的故障数据更新支持向量机模型,最后将目标域单类未带标签的故障数据输入支持向量机模型,对目标域进行第二次故障分类。
优选地,所述步骤S1具体包括以下步骤:
S11,采集污水处理过程数据,包含两类工况数据:常见工况即源域的正常数据集Xsc和带标签的故障数据集Xsf,非常见工况即目标域的少量的正常数据集Xtc和少量带标签的故障数据集Xtf;
S12,分离出非平稳变量的过程数据;首先建立平稳性检验模型,对每组单变量的数据进行平稳性检验;根据平稳性检验结果,分离出非平稳变量的过程数据,获得常见工况即源域的非平稳变量的正常数据集和带标签的非平稳变量的故障数据集/>获得非常见工况即目标域的非平稳变量的少量的正常数据集/>和少量带标签的非平稳变量的故障数据集/>
S13,基于邻域保持成分分析的方法提取非平稳特征;对S12获得的数据进行归一化处理,得到源域的非平稳变量的正常数据集为:
其中表示/>的列向量,ns1表示源域的非平稳变量的正常数据集的样本数,ms表示源域的非平稳变量的正常数据集的变量数;
源域的非平稳变量的带标签的故障数据集为:
其中表示/>的列向量,ns2表示源域的非平稳变量的带标签的故障数据集的样本数;
目标域的非平稳变量的正常数据集为:
其中表示/>的列向量,nt1表示目标域的非平稳变量的正常数据集的样本数,mt表示目标域的非平稳变量的正常数据集的变量数;
目标域的非平稳变量的带标签的故障数据集为:
其中表示/>的列向量,nt2表示目标域的非平稳变量的带标签的故障数据集的样本数;
以提取源域的正常数据集的非平稳特征为例,根据给定的核函数k,计算正定核矩阵K:
其中,Kij表示正定核矩阵中的元素,表示/>中的任意两列样本,/>表示从原始特征空间到高维特征空间的一个核映射,/>表示/>的转置;
计算中心核矩阵
其中,Ik表示全为的ms阶方阵;
计算矩阵所有特征值,对这些特征值进行排序,选择前p个特征值λ1≥λ2≥…≥λp,然后计算它们所对应的特征向量v1,v2,…,vp。特征提取矩阵W的第i个列向量wi为:
其中λi和vi分别表示第i个特征值和对应的第i个特征向量,W=[w1,w2,…,wp];
因此,的输出特征集Ysc的列向量ysci为:
其中Ysc=[ysc1,ysc2,…,yscp],表示/>中的所有列向量,p表示保留的特征值的个数;
利用K最近邻算法构建邻域图;Ysc=[ysc1,ysc2,…,yscp],则邻域图共有p个结点,其中ysci表示第i个结点;如果yscj是ysci的k个最近邻居之一,那么将这两点连接,反之不连接;设权重矩阵为U,其中元素Uij代表结点i和结点j之间的边的权重,如果两点之间没有边,则对应的矩阵元素为0;矩阵U的元素值Uij主要通过最小化如下目标函数得到:
其中U应满足归一化约束:
邻域保持算法的思想是降维后的特征空间与原始高维空间具有相似的局部结构;通过求解广义特征向量问题来计算映射矩阵A:
其中矩阵M=(IU-U)T(IU-U),矩阵IU=diag(1,…,1),是特征值,a是特征向量,所有特征值按照从小到大的次序排列并保留前d个特征值/>对应的前d个特征值对应的特征向量为a1,a2,…,ad,因此邻域保持算法降维后的数据为:
其中A=(a1,a2,…,ad),降维后的数据矩阵为其中d表示矩阵/>的变量的个数;
因此,源域的正常数据的非平稳特征为源域的带标签的故障数据集的非平稳特征为/>目标域的少量的正常数据的非平稳特征为/>目标域的少量带标签的故障数据的非平稳特征为/>
优选地,所述步骤S2具体包括以下步骤:
根据步骤S1获得源域的正常数据的非平稳特征为目标域的少量的正常数据的非平稳特征为/>将典型相关分析扩展到求解一对投影矩阵Ps,Pt的情况,因此有以下目标函数:
Zsc和Ztc有如下约束:
ZscZsc T=Id,ZtcZtc T=Id
其中,Zsc和Ztc分别是和/>的投影,zsc和ztc分别是Zsc和Ztc的元素,trace表示矩阵的迹,Id表示全为1的d阶方阵;
因此,上述目标函数表示为如下优化问题:
s.t.ZscZsc T=Id
ZtcZtc T=Id.
源域的正常数据的非平稳特征变为源域的带标签的故障数据集的非平稳特征变为/>目标域的少量的正常数据的非平稳特征变为/>目标域的少量带标签的故障数据的非平稳特征变为/>
优选地,所述步骤S3具体包括以下步骤:
S31,根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图;
源域的正常数据的非平稳特征变为Zsc,给定非平稳时间序列Zsc,计算参考均值和参考标准差:
其中i表示第i个变量,μi是均值,σi是标准差,ns1是源域的正常数据的样本数,Iz是具有适当维数的全为1的向量,t是采样时刻,zsc,i是第i个变量,zsc,t,i是第t个采样时刻的第i个变量的值,因此,原始数据Zsc被归一化为:
其中表示归一化的数据,μ=(μ1,μ2,…,μd),Σ=diag{σ1,σ2,…,σd};
利用Johansen多变量协整分析方法对选取的数据进行协整关系检验;如果一个时间序列经过一次差分变成平稳的,则称原序列是1阶单整的,记为I(1);假设所有时间序列/>都是1阶单整的。矢量误差修正模型描述为:
式中
其中 是第t个时刻的时间序列,/>是第t-1个时刻的时间序列;/> 是第t-i个时刻的时间序列,/>是第t-i-1个时刻的时间序列;εt为白噪声向量并服从高斯分布N(0,Ξ),Ξ是高斯分布的方差,p是模型的阶次;Πi和Πj为系数矩阵,Im为随机矩阵,Im可以分解为Υ和Bf两个列满秩矩阵,Γ=ΥBf T,Bf是协积分矩阵,Bf的列向量是协积分向量,Υ是一个辅助矩阵;
根据Johansen多变量协整分析方法,Bf可以通过下面的似然函数L求解:
利用特征值分解可以得到Bf的最大似然估计:
其中E为预测误差, 为特征值,系数矩阵Θi和Φi通过最小二乘方法得到;
在线监测时,q+1时刻的长期的动态平稳关系表示为:
其中EB为反映了长期动态平稳关系的误差,Bf,q为第q个时刻的协积分矩阵,表示第q+1个时刻的时间序列;
因此,监测统计量L2可以表示为:
L2=EB TEB
其中,μL和sL分别为L2的均值和方差,表示卡方分布,g和h是参数,置信水平为α;
监测统计量L2还可以表示为:
第i个变量对L2的贡献可以表示为:
其中eB,i表示EB的第i个元素,d表示源域正常数据的非平稳特征Zsc的变量个数;
利用上述的故障检测模型画出源域带标签的故障数据的贡献图,作为先验知识储存;利用上述的故障检测模型再获取目标域未带标签的故障数据,进一步画出目标域未带标签的故障数据的贡献图;
S32,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析;
源域和目标域的贡献图可以看作是很多副不同的直方图图像,计算目标域未带标签的故障数据的贡献图与源域所有带标签的故障数据的贡献图之间的图像相似度,找出与目标域贡献图的图像相似度最大的那一类源域的贡献图,此时对应的源域的故障标签即为目标域的故障类型;
通过对图像进行处理,利用欧式距离表示图像之间的相似度,欧氏距离越小,相似度越大;n维空间中的欧式距离公式为:
其中DEuc表示欧式距离,x和y表示n维空间中的两个点,x=(x1,…,xn),y=(y1,…,yn);
通过对图像进行处理,也可以用余弦相似度表示图像之间的相似度,余弦相似度越大,相似度越大;余弦相似度把图片表示成一个向量,通过计算向量之间的余弦值来表征两张图片的相似度,余弦值接近1,夹角趋于0,表明两个向量越相似,余弦值接近于0,夹角趋于90度,表明两个向量越不相似;n维空间中的余弦相似度公式为:
其中DCos表示余弦相似度,x和y表示n维空间中的两个点,x=(x1,…,xn),y=(y1,…,yn);
因此,基于欧式距离和余弦相似度的图像相似度的计算公式为:
其中DCE表示图像相似度,DCE越大,图像之间的相似度越大。
优选地,所述步骤S4具体包括以下步骤:
支持向量机进行非线性分类时通过高维空间变换,将非线性分类变成高维空间的线性分类问题;假设大小为l的训练样本集{(hi,vi),i=1,2,...,l}由两个类别组成,通过非线性映射Φ,将原始空间的样本映射到高维的特征空间中,对应的对偶问题变为:
其中hi表示原始空间的样本,Φ(hi)表示hi映射到高维特征空间后的样本,vi是原始空间的样本对应的标签,F(β)是对偶函数,β=[β1,…,βl]T是拉格朗日乘子,<Φ(hi),Φ(hj)>表示内积,C为惩罚因子;
设β的最优解为则/>其中ω*是一个实数向量,从而得到最优的分类模型为:
其中是f(h)分类函数,sgn表示阶跃函数,b*是一个常数,<Φ(h),Φ(hi)>表示内积,表示核函数,这里采用高斯核函数。
本发明实施例提供的技术方案带来的有益效果至少包括:
1)针对多工况的污水处理过程,基于邻域保持成分分析的特征提取方法提取过程数据的非平稳特征,可以减小常见工况和非常见工况下提取到的特征的差异,降低工况变化对***模型的影响。
2)针对多工况的污水处理过程,基于典型相关分析的迁移学习方法,解决了非常见工况因历史数据不足而无法建立故障诊断模型的问题。
3)针对多工况的污水处理过程,开发了一种非常见工况下的污水处理过程的故障分类方法。贡献图分类是第一次故障分类,支持向量机分类是第二次故障分类,对第一次分类过程中未成功分类的故障进行分类。两者结合大大提高了在线监测时故障分类的效率,降低了模型复杂度且缩短了离线训练的时间。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种非常见工况下的污水处理过程的故障分类方法的流程图;
图2是根据本发明一实施例的污水处理过程的结构示意图;
图3是根据本发明一实施例的雨天模式下带标签的故障A的贡献图;
图4是根据本发明一实施例的雨天模式下带标签的故障B的贡献图;
图5是根据本发明一实施例的雨天模式下带标签的故障C的贡献图;
图6是根据本发明一实施例的雨天模式下带标签的故障D的贡献图;
图7是根据本发明一实施例的目标域未带标签的故障1和故障2的支持向量机故障分类结果。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例的附图,对本发明实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的实施例提供了一种非常见工况下的污水处理过程的故障分类方法,如图1所示,所述方法包括以下步骤:
S1,采集污水处理过程数据,包括源域和目标域的历史数据,数据预处理后分离出非平稳变量的过程数据,基于邻域保持成分分析的方法提取非平稳特征。
所述步骤S1具体包括:
S11,采集污水处理过程数据,包含两类工况数据:常见工况即源域的正常数据集Xsc和带标签的故障数据集Xsf,非常见工况即目标域的少量的正常数据集Xtc和少量带标签的故障数据集Xtf;
S12,分离出非平稳变量的过程数据;首先建立平稳性检验模型,对每组单变量的数据进行平稳性检验;根据平稳性检验结果,分离出非平稳变量的过程数据,获得常见工况即源域的非平稳变量的正常数据集和带标签的非平稳变量的故障数据集/>获得非常见工况即目标域的非平稳变量的少量的正常数据集/>和少量带标签的非平稳变量的故障数据集/>
S13,基于邻域保持成分分析的方法提取非平稳特征;对S12获得的数据进行归一化处理,得到源域的非平稳变量的正常数据集为:
其中表示/>的列向量,ns1表示源域的非平稳变量的正常数据集的样本数,ms表示源域的非平稳变量的正常数据集的变量数;
源域的非平稳变量的带标签的故障数据集为:
其中表示/>的列向量,ns2表示源域的非平稳变量的带标签的故障数据集的样本数;
目标域的非平稳变量的正常数据集为:
其中表示/>的列向量,nt1表示目标域的非平稳变量的正常数据集的样本数,mt表示目标域的非平稳变量的正常数据集的变量数;
目标域的非平稳变量的带标签的故障数据集为:
其中表示/>的列向量,nt2表示目标域的非平稳变量的带标签的故障数据集的样本数;
以提取源域的正常数据集的非平稳特征为例,根据给定的核函数k,计算正定核矩阵K:
其中,Kij表示正定核矩阵中的元素,表示/>中的任意两列样本,/>表示从原始特征空间到高维特征空间的一个核映射,/>表示/>的转置;
计算中心核矩阵
其中,Ik表示全为的ms阶方阵;
计算矩阵所有特征值,对这些特征值进行排序,选择前p个特征值λ1≥λ2≥…≥λp,然后计算它们所对应的特征向量v1,v2,…,vp。特征提取矩阵W的第i个列向量wi为:
其中λi和vi分别表示第i个特征值和对应的第i个特征向量,W=[w1,w2,…,wp];
因此,的输出特征集Ysc的列向量ysci为:
其中Ysc=[ysc1,ysc2,…,yscp],表示/>中的所有列向量,p表示保留的特征值的个数;
利用K最近邻算法构建邻域图;Ysc=[ysc1,ysc2,…,yscp],则邻域图共有p个结点,其中ysci表示第i个结点;如果yscj是ysci的k个最近邻居之一,那么将这两点连接,反之不连接;设权重矩阵为U,其中元素Uij代表结点i和结点j之间的边的权重,如果两点之间没有边,则对应的矩阵元素为0;矩阵U的元素值Uij主要通过最小化如下目标函数得到:
其中U应满足归一化约束:
邻域保持算法的思想是降维后的特征空间与原始高维空间具有相似的局部结构;通过求解广义特征向量问题来计算映射矩阵A:
其中矩阵M=(IU-U)T(IU-U),矩阵IU=diag(1,…,1),是特征值,a是特征向量,所有特征值按照从小到大的次序排列并保留前d个特征值/>对应的前d个特征值对应的特征向量为a1,a2,…,ad,因此邻域保持算法降维后的数据为:
其中A=(a1,a2,…,ad),降维后的数据矩阵为其中d表示矩阵/>的变量的个数;
因此,源域的正常数据的非平稳特征为源域的带标签的故障数据集的非平稳特征为/>目标域的少量的正常数据的非平稳特征为/>目标域的少量带标签的故障数据的非平稳特征为/>
S2,基于典型相关分析的迁移学习方法减小源域和目标域正常数据的分布差异,获得投影后的源域和目标域的数据。
所述步骤S2具体包括:
根据步骤S1获得源域的正常数据的非平稳特征为目标域的少量的正常数据的非平稳特征为/>将典型相关分析扩展到求解一对投影矩阵Ps,Pt的情况,因此有以下目标函数:
Zsc和Ztc有如下约束:
ZscZsc T=Id,ZtcZtc T=Id
其中,Zsc和Ztc分别是和/>的投影,zsc和ztc分别是Zsc和Ztc的元素,trace表示矩阵的迹,Id表示全为1的d阶方阵;
因此,上述目标函数表示为如下优化问题:
s.t.ZscZsc T=Id
ZtcZtc T=Id.
源域的正常数据的非平稳特征变为源域的带标签的故障数据集的非平稳特征变为/>目标域的少量的正常数据的非平稳特征变为/>目标域的少量带标签的故障数据的非平稳特征变为/>
S3,根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析,对目标域进行第一次故障分类。
所述步骤S3具体包括:
S31,根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图;
源域的正常数据的非平稳特征变为Zsc,给定非平稳时间序列Zsc,计算参考均值和参考标准差:
其中i表示第i个变量,μi是均值,σi是标准差,ns1是源域的正常数据的样本数,Iz是具有适当维数的全为1的向量,t是采样时刻,zsc,i是第i个变量,zsc,t,i是第t个采样时刻的第i个变量的值,因此,原始数据Zsc被归一化为:
其中表示归一化的数据,μ=(μ1,μ2,…,μd),∑=diag{σ1,σ2,…,σd};
利用Johansen多变量协整分析方法对选取的数据进行协整关系检验;如果一个时间序列经过一次差分变成平稳的,则称原序列是1阶单整的,记为I(1);假设所有时间序列/>都是1阶单整的。矢量误差修正模型描述为:
式中
其中 是第t个时刻的时间序列,/>是第t-1个时刻的时间序列;/> 是第t-i个时刻的时间序列,/>是第t-i-1个时刻的时间序列;εt为白噪声向量并服从高斯分布N(0,Ξ),Ξ是高斯分布的方差,p是模型的阶次;Πi和Πj为系数矩阵,Im为随机矩阵,Im可以分解为Υ和Bf两个列满秩矩阵,Γ=ΥBf T,Bf是协积分矩阵,Bf的列向量是协积分向量,Υ是一个辅助矩阵;
根据Johansen多变量协整分析方法,Bf可以通过下面的似然函数L求解:
利用特征值分解可以得到Bf的最大似然估计:
其中E为预测误差, 为特征值,系数矩阵Θi和Φi通过最小二乘方法得到;
在线监测时,q+1时刻的长期的动态平稳关系表示为:
其中EB为反映了长期动态平稳关系的误差,Bf,q为第q个时刻的协积分矩阵,表示第q+1个时刻的时间序列;
因此,监测统计量L2可以表示为:
其中,μL和sL分别为L2的均值和方差,表示卡方分布,g和h是参数,置信水平为α;
监测统计量L2还可以表示为:
第i个变量对L2的贡献可以表示为:
其中eB,i表示EB的第i个元素,d表示源域正常数据的非平稳特征Zsc的变量个数;
利用上述的故障检测模型画出源域带标签的故障数据的贡献图,作为先验知识储存;利用上述的故障检测模型再获取目标域未带标签的故障数据,进一步画出目标域未带标签的故障数据的贡献图;
S32,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析;
源域和目标域的贡献图可以看作是很多副不同的直方图图像,计算目标域未带标签的故障数据的贡献图与源域所有带标签的故障数据的贡献图之间的图像相似度,找出与目标域贡献图的图像相似度最大的那一类源域的贡献图,此时对应的源域的故障标签即为目标域的故障类型;
通过对图像进行处理,利用欧式距离表示图像之间的相似度,欧氏距离越小,相似度越大;n维空间中的欧式距离公式为:
其中DEuc表示欧式距离,x和y表示n维空间中的两个点,x=(x1,…,xn),y=(y1,…,yn);
通过对图像进行处理,也可以用余弦相似度表示图像之间的相似度,余弦相似度越大,相似度越大;余弦相似度把图片表示成一个向量,通过计算向量之间的余弦值来表征两张图片的相似度,余弦值接近1,夹角趋于0,表明两个向量越相似,余弦值接近于0,夹角趋于90度,表明两个向量越不相似;n维空间中的余弦相似度公式为:
其中DCos表示余弦相似度,x和y表示n维空间中的两个点,x=(x1,…,xn),y=(y1,…,yn);
因此,基于欧式距离和余弦相似度的图像相似度的计算公式为:
其中DCE表示图像相似度,DCE越大,图像之间的相似度越大。
S4,如果目标域单类未带标签的故障数据的贡献图与源域多类带标签的故障数据的贡献图的相似度相同,将对应的投影后的源域多类带标签的故障数据训练支持向量机模型;根据对应的投影后的目标域少量带标签的故障数据更新支持向量机模型,最后将目标域单类未带标签的故障数据输入支持向量机模型,对目标域进行第二次故障分类。
所述步骤S4具体包括:
支持向量机进行非线性分类时通过高维空间变换,将非线性分类变成高维空间的线性分类问题;假设大小为l的训练样本集{(hi,vi),i=1,2,…,l}由两个类别组成,通过非线性映射Φ,将原始空间的样本映射到高维的特征空间中,对应的对偶问题变为:
其中hi表示原始空间的样本,Φ(hi)表示hi映射到高维特征空间后的样本,vi是原始空间的样本对应的标签,F(β)是对偶函数,β=[β1,…,βl]T是拉格朗日乘子,<Φ(hi),Φ(hj)>表示内积,C为惩罚因子;
设β的最优解为则/>其中ω*是一个实数向量,从而得到最优的分类模型为:/>
其中是f(h)分类函数,sgn表示阶跃函数,b*是一个常数,<Φ(h),Φ(hi)>表示内积,表示核函数,这里采用高斯核函数。
本发明实施例中,针对多工况的污水处理过程,在非常见工况缺乏大量历史数据的条件下,利用污水处理过程的常见工况的历史数据,通过迁移学习实现对非常见工况的故障分类。
下面以收集来自污水处理过程的BMS1模型的运行数据为例,对本发明提供的一种非常见工况下的污水处理过程的故障分类方法的有效性进行详细的说明,具体步骤包括:
步骤1、建立污水处理过程的BMS1模型,采集污水处理过程中的实际运行数据,一共有晴天、雨天、暴雨三种工况,三种工况可以分为常见工况和非常见工况,假设常见工况为源域,非常见工况为目标域,对数据集进行预处理。在预处理后的过程数据中分离出非平稳变量构成非平稳数据集,基于邻域保持成分分析的方法提取非平稳特征。
步骤2、基于典型相关分析的迁移学习方法减小源域和目标域正常数据的分布差异,获得投影后的源域和目标域的正常和故障数据集。
步骤3、根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析。
步骤4、如果目标域单类未带标签的故障数据的贡献图与源域多类带标签的故障数据的贡献图的相似度相同,将对应的投影后的源域多类带标签的故障数据训练支持向量机模型;根据对应的投影后的目标域少量带标签的故障数据集更新支持向量机模型,最后将目标域单类未带标签的故障数据输入支持向量机模型。
采用污水处理过程的BMS1模型的运行数据进行仿真。如图2所示,污水处理模型由1个生物反应器和1个二次沉淀池组成,生物反应器包括5个混合完全的小单元。整个污水处理过程共收集27个过程变量,在数据预处理时分离出了共7个非平稳变量,整个污水处理过程共采集了14天的运行数据,数据采样间隔为4次/小时。
所有数据根据工况可以划分为晴天、雨天、暴雨三种数据集,每一种工况的运行数据包括正常数据和故障数据。三种工况可以分为常见工况和非常见工况,假设常见工况雨天为源域,非常见工况暴雨为目标域,源域含有大量带标签的故障数据和正常数据,目标域含有少量带标签的故障数据和正常数据,假设源域和目标域的故障的类型相同,共有四类故障。
源域带标签的四类故障(故障A、故障B、故障C、故障D)的贡献图与目标域未带标签的四类故障(故障1、故障2、故障3、故障4)的贡献图的图像相似度分析的结果如表1所示,根据表格结果得到,目标域中只有故障3和故障4可以通过贡献图分类确定故障类型,目标域故障3的类型是源域故障C,目标域故障4的类型是源域故障D,目标域的故障1与源域的故障A和故障B的图像相似度基本一致,目标域的故障2也与源域的故障A和故障B的图像相似度基本一致。
表1源域和目标域四类贡献图图像相似度分析结果
图像相似度 | 源域故障A | 源域故障B | 源域故障C | 源域故障D |
目标域故障1 | 49.2799 | 50.6801 | 2.0672 | 2.7584 |
目标域故障2 | 55.9722 | 50.3952 | 5.7263 | 2.3676 |
目标域故障3 | 5.0351 | 5.2944 | 185.7077 | 4.5461 |
目标域故障4 | 1.4448 | 1.5135 | 2.3305 | 53.4451 |
因此,目标域的故障1和故障2需要在支持向量机模型中进行分类。首先利用源域带标签的故障A和故障B的数据训练支持向量机模型,再利用目标域少量带标签的故障A和故障B的数据训练模型,微调模型参数,直到模型分类的准确率最高。支持向量机模型的训练样本集为190,测试集样本为82。本实施例对基于贡献图和支持向量机的非常见工况下的故障分类方法的实验结果如图3、图4、图5、图6和图7所示。其中,图3是根据本实施例画出的源域带标签的故障A的贡献图,图4是根据本实施例画出的源域带标签的故障B的贡献图,图5是根据本实施例画出的源域带标签的故障C的贡献图,图6是根据本实施例画出的源域带标签的故障D的贡献图,图7是根据本实施例,假设常见工况雨天为源域,非常见工况暴雨为目标域,把源域训练好的支持向量机模型迁移到目标域,画出的支持向量机模型对目标域未带标签的故障1和故障2的分类结果。
与现有技术相比,本发明提供的一种非常见工况下的污水处理过程的故障分类方法,具有以下有益效果:
1)针对多工况的污水处理过程,基于邻域保持成分分析的特征提取方法提取过程数据的非平稳特征,可以减小常见工况和非常见工况下提取到的特征的差异,降低工况变化对***模型的影响。
2)针对多工况的污水处理过程,基于典型相关分析的迁移学习方法,解决了非常见工况因历史数据不足而无法建立故障诊断模型的问题。
3)针对多工况的污水处理过程,开发了一种非常见工况下的污水处理过程的故障分类方法。贡献图分类是第一次故障分类,支持向量机分类是第二次故障分类,对第一次分类过程中未成功分类的故障进行分类。两者结合大大提高了在线监测时故障分类的效率,降低了模型复杂度且缩短了离线训练的时间。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。
在说明书中提到“一个实施例”、“实施例”、“示例性实施例”、“一些实施例”等指示所述的实施例可以包括特定特征、结构或特性,但未必每个实施例都包括该特定特征、结构或特性。另外,在结合实施例描述特定特征、结构或特性时,结合其它实施例(无论是否明确描述)实现这种特征、结构或特性应在相关领域技术人员的知识范围内。
本发明涵盖任何在本发明的精髓和范围上做的替代、修改、等效方法以及方案。为了使公众对本发明有彻底的了解,在以下本发明优选实施例中详细说明了具体的细节,而对本领域技术人员来说没有这些细节的描述也可以完全理解本发明。另外,为了避免对本发明的实质造成不必要的混淆,并没有详细说明众所周知的方法、过程、流程、元件和电路等。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于计算机可读取存储介质中,如:ROM/RAM、磁碟、光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种非常见工况下的污水处理过程的故障分类方法,其特征在于,包括以下步骤:
S1,采集污水处理过程数据,包括源域和目标域的历史数据,数据预处理后分离出非平稳变量的过程数据,基于邻域保持成分分析的方法提取非平稳特征;
所述步骤S1具体包括以下步骤:
S11,采集污水处理过程数据,包含两类工况数据:常见工况即源域的正常数据集Xsc和带标签的故障数据集Xsf,非常见工况即目标域的少量的正常数据集Xtc和少量带标签的故障数据集Xtf;
S12,分离出非平稳变量的过程数据;首先建立平稳性检验模型,对每组单变量的数据进行平稳性检验;根据平稳性检验结果,分离出非平稳变量的过程数据,获得常见工况即源域的非平稳变量的正常数据集和带标签的非平稳变量的故障数据集/>获得非常见工况即目标域的非平稳变量的少量的正常数据集/>和少量带标签的非平稳变量的故障数据集/>
S13,基于邻域保持成分分析的方法提取非平稳特征;对S12获得的数据进行归一化处理,得到源域的非平稳变量的正常数据集为:
其中表示/>的列向量,ns1表示源域的非平稳变量的正常数据集的样本数,ms表示源域的非平稳变量的正常数据集的变量数;
源域的非平稳变量的带标签的故障数据集为:
其中表示/>的列向量,ns2表示源域的非平稳变量的带标签的故障数据集的样本数;
目标域的非平稳变量的正常数据集为:
其中表示/>的列向量,nt1表示目标域的非平稳变量的正常数据集的样本数,mt表示目标域的非平稳变量的正常数据集的变量数;
目标域的非平稳变量的带标签的故障数据集为:
其中表示/>的列向量,nt2表示目标域的非平稳变量的带标签的故障数据集的样本数;
以提取源域的正常数据集的非平稳特征为例,根据给定的核函数k,计算正定核矩阵K:
其中,Kij表示正定核矩阵中的元素,表示/>中的任意两列样本,/>表示从原始特征空间到高维特征空间的一个核映射,/>表示/>的转置;
计算中心核矩阵
其中,Ik表示全为的ms阶方阵;
计算矩阵所有特征值,对这些特征值进行排序,选择前p个特征值λ1≥λ2≥···≥λp,然后计算它们所对应的特征向量v1,v2,···,vp;特征提取矩阵W的第i个列向量wi为:
其中λi和vi分别表示第i个特征值和对应的第i个特征向量,W=[w1,w2,···,wp];
因此,的输出特征集Ysc的列向量ysci为:
其中Ysc=[ysc1,ysc2,…,yscp],表示/>中的所有列向量,p表示保留的特征值的个数;
利用K最近邻算法构建邻域图;Ysc=[ysc1,ysc2,…,yscp],则邻域图共有p个结点,其中ysci表示第i个结点;如果yscj是ysci的k个最近邻居之一,那么将这两点连接,反之不连接;设权重矩阵为U,其中元素Uij代表结点i和结点j之间的边的权重,如果两点之间没有边,则对应的矩阵元素为0;矩阵U的元素值Uij主要通过最小化如下目标函数得到:
其中U应满足归一化约束:
邻域保持算法的思想是降维后的特征空间与原始高维空间具有相似的局部结构;通过求解广义特征向量问题来计算映射矩阵A:
其中矩阵M=(IU-U)T(IU-U),矩阵IU=diag(1,…,1),是特征值,a是特征向量,所有特征值按照从小到大的次序排列并保留前d个特征值/>对应的前d个特征值对应的特征向量为a1,a2,…,ad,因此邻域保持算法降维后的数据为:
其中A=(a1,a2,…,ad),降维后的数据矩阵为其中d表示矩阵/>的变量的个数;
因此,源域的正常数据的非平稳特征为源域的带标签的故障数据集的非平稳特征为/>目标域的少量的正常数据的非平稳特征为/>目标域的少量带标签的故障数据的非平稳特征为/>
S2,基于典型相关分析的迁移学习方法减小源域和目标域正常数据的分布差异,获得投影后的源域和目标域的数据;
S3,根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析,对目标域进行第一次故障分类;
S4,如果目标域单类未带标签的故障数据的贡献图与源域多类带标签的故障数据的贡献图的相似度相同,将对应的投影后的源域多类带标签的故障数据训练支持向量机模型;根据对应的投影后的目标域少量带标签的故障数据更新支持向量机模型,最后将目标域单类未带标签的故障数据输入支持向量机模型,对目标域进行第二次故障分类。
2.根据权利要求1所述的故障分类方法,其特征在于,所述步骤S2具体包括以下步骤:
根据步骤S1获得源域的正常数据的非平稳特征为目标域的少量的正常数据的非平稳特征为/>将典型相关分析扩展到求解一对投影矩阵Ps,Pt的情况,因此有以下目标函数:
Zsc和Ztc有如下约束:
ZscZsc T=Id,ZtcZtc T=Id
其中,Zsc和Ztc分别是和/>的投影,zsc和ztc分别是Zsc和Ztc的元素,trace表示矩阵的迹,Id表示全为1的d阶方阵;
因此,上述目标函数表示为如下优化问题:
s.t.ZscZsc T=Id
ZtcZtc T=Id.
源域的正常数据的非平稳特征变为源域的带标签的故障数据集的非平稳特征变为/>目标域的少量的正常数据的非平稳特征变为/>目标域的少量带标签的故障数据的非平稳特征变为/>
3.根据权利要求2所述的故障分类方法,其特征在于,所述步骤S3具体包括以下步骤:
S31,根据投影后的源域数据集训练故障检测模型,画出源域带标签的故障数据的贡献图和目标域未带标签的故障数据的贡献图;
源域的正常数据的非平稳特征变为Zsc,给定非平稳时间序列Zsc,计算参考均值和参考标准差:
其中i表示第i个变量,μi是均值,σi是标准差,ns1是源域的正常数据的样本数,Iz是具有适当维数的全为1的向量,t是采样时刻,zsc,i是第i个变量,zsc,t,i是第t个采样时刻的第i个变量的值,因此,原始数据Zsc被归一化为:
其中表示归一化的数据,μ=(μ1,μ2,…,μd),∑=diag{σ1,σ2,…,σd};
利用Johansen多变量协整分析方法对选取的数据进行协整关系检验;如果一个时间序列经过一次差分变成平稳的,则称原序列是1阶单整的,记为I(1);假设所有时间序列/>都是1阶单整的;矢量误差修正模型描述为:
式中
其中 是第t个时刻的时间序列,/>是第t-1个时刻的时间序列; 是第t-i个时刻的时间序列,/>是第t-i-1个时刻的时间序列;εt为白噪声向量并服从高斯分布N(0,Ξ),Ξ是高斯分布的方差,p是模型的阶次;Πi和Πj为系数矩阵,Im为随机矩阵,Im可以分解为Υ和Bf两个列满秩矩阵,Γ=ΥBf T,Bf是协积分矩阵,Bf的列向量是协积分向量,Υ是一个辅助矩阵;
根据Johansen多变量协整分析方法,Bf可以通过下面的似然函数L求解:
利用特征值分解可以得到Bf的最大似然估计:
其中E为预测误差, 为特征值,系数矩阵Θi和Φi通过最小二乘方法得到;
在线监测时,q+1时刻的长期的动态平稳关系表示为:
其中EB为反映了长期动态平稳关系的误差,Bf,q为第q个时刻的协积分矩阵,表示第q+1个时刻的时间序列;
因此,监测统计量L2可以表示为:
其中,μL和sL分别为L2的均值和方差,表示卡方分布,g和h是参数,置信水平为α;
监测统计量L2还可以表示为:
第i个变量对L2的贡献可以表示为:
其中eB,i表示EB的第i个元素,d表示源域正常数据的非平稳特征Zsc的变量个数;
利用上述的故障检测模型画出源域带标签的故障数据的贡献图,作为先验知识储存;利用上述的故障检测模型再获取目标域未带标签的故障数据,进一步画出目标域未带标签的故障数据的贡献图;
S32,利用欧式距离和余弦距离对源域和目标域的贡献图做相似度分析;
源域和目标域的贡献图可以看作是很多副不同的直方图图像,计算目标域未带标签的故障数据的贡献图与源域所有带标签的故障数据的贡献图之间的图像相似度,找出与目标域贡献图的图像相似度最大的那一类源域的贡献图,此时对应的源域的故障标签即为目标域的故障类型;
通过对图像进行处理,利用欧式距离表示图像之间的相似度,欧氏距离越小,相似度越大;n维空间中的欧式距离公式为:
其中DEuc表示欧式距离,x和y表示n维空间中的两个点,x=(x1,…,xn),y=(y1,…,yn);
通过对图像进行处理,也可以用余弦相似度表示图像之间的相似度,余弦相似度越大,相似度越大;余弦相似度把图片表示成一个向量,通过计算向量之间的余弦值来表征两张图片的相似度,余弦值接近1,夹角趋于0,表明两个向量越相似,余弦值接近于0,夹角趋于90度,表明两个向量越不相似;n维空间中的余弦相似度公式为:
其中DCos表示余弦相似度,x和y表示n维空间中的两个点,x=(x1,…,xn),y=(y1,…,yn);
因此,基于欧式距离和余弦相似度的图像相似度的计算公式为:
其中DCE表示图像相似度,DCE越大,图像之间的相似度越大。
4.根据权利要求3所述的故障分类方法,其特征在于,所述步骤S4具体包括以下步骤:
支持向量机进行非线性分类时通过高维空间变换,将非线性分类变成高维空间的线性分类问题;假设大小为l的训练样本集{(hi,vi),i=1,2,...,l}由两个类别组成,通过非线性映射Φ,将原始空间的样本映射到高维的特征空间中,对应的对偶问题变为:
其中hi表示原始空间的样本,Φ(hi)表示hi映射到高维特征空间后的样本,vi是原始空间的样本对应的标签,F(β)是对偶函数,β=[β1,…,βl]T是拉格朗日乘子,<Φ(hi),Φ(hj)>表示内积,C为惩罚因子;
设β的最优解为则/>其中ω*是一个实数向量,从而得到最优的分类模型为:
其中是f(h)分类函数,sgn表示阶跃函数,b*是一个常数,<Φ(h),Φ(hi)>表示内积,表示核函数,这里采用高斯核函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310130653.0A CN116206154B (zh) | 2023-02-17 | 2023-02-17 | 一种非常见工况下的污水处理过程的故障分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310130653.0A CN116206154B (zh) | 2023-02-17 | 2023-02-17 | 一种非常见工况下的污水处理过程的故障分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116206154A CN116206154A (zh) | 2023-06-02 |
CN116206154B true CN116206154B (zh) | 2024-02-13 |
Family
ID=86510733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310130653.0A Active CN116206154B (zh) | 2023-02-17 | 2023-02-17 | 一种非常见工况下的污水处理过程的故障分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116206154B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414226A (zh) * | 2017-12-25 | 2018-08-17 | 哈尔滨理工大学 | 基于特征迁移学习的变工况下滚动轴承故障诊断方法 |
CN111914888A (zh) * | 2020-06-13 | 2020-11-10 | 宁波大学 | 一种多工况识别与故障检测一体化的化工过程监测方法 |
CN112784872A (zh) * | 2020-12-25 | 2021-05-11 | 北京航空航天大学 | 一种基于开放集联合迁移学习的跨工况故障诊断方法 |
CN112964469A (zh) * | 2021-02-28 | 2021-06-15 | 哈尔滨理工大学 | 一种迁移学习的变负载下滚动轴承在线故障诊断方法 |
CN114004251A (zh) * | 2021-10-12 | 2022-02-01 | 北京科技大学 | 基于二维灰度图与迁移学习的冷连轧断带故障诊断方法 |
CN114842383A (zh) * | 2022-04-29 | 2022-08-02 | 合肥工业大学 | 一种基于dvrl-vst的电熔镁炉工况视频识别方法 |
CN115439710A (zh) * | 2022-08-23 | 2022-12-06 | 北京数慧时空信息技术有限公司 | 基于组合迁移学习的遥感样本标注方法 |
CN115563536A (zh) * | 2022-09-28 | 2023-01-03 | 北京工业大学 | 基于子域自适应的滚动轴承故障诊断方法 |
CN115563494A (zh) * | 2022-12-01 | 2023-01-03 | 北京科技大学 | 基于迁移学习的多模式流程工业过程的数据驱动故障检测方法 |
CN115687925A (zh) * | 2022-10-31 | 2023-02-03 | 中国船舶重工集团公司第七一九研究所 | 针对不平衡样本的故障类型识别方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112414714A (zh) * | 2020-11-05 | 2021-02-26 | 苏州大学 | 基于自适应流形概率分布的轴承故障诊断方法 |
-
2023
- 2023-02-17 CN CN202310130653.0A patent/CN116206154B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414226A (zh) * | 2017-12-25 | 2018-08-17 | 哈尔滨理工大学 | 基于特征迁移学习的变工况下滚动轴承故障诊断方法 |
CN111914888A (zh) * | 2020-06-13 | 2020-11-10 | 宁波大学 | 一种多工况识别与故障检测一体化的化工过程监测方法 |
CN112784872A (zh) * | 2020-12-25 | 2021-05-11 | 北京航空航天大学 | 一种基于开放集联合迁移学习的跨工况故障诊断方法 |
CN112964469A (zh) * | 2021-02-28 | 2021-06-15 | 哈尔滨理工大学 | 一种迁移学习的变负载下滚动轴承在线故障诊断方法 |
CN114004251A (zh) * | 2021-10-12 | 2022-02-01 | 北京科技大学 | 基于二维灰度图与迁移学习的冷连轧断带故障诊断方法 |
CN114842383A (zh) * | 2022-04-29 | 2022-08-02 | 合肥工业大学 | 一种基于dvrl-vst的电熔镁炉工况视频识别方法 |
CN115439710A (zh) * | 2022-08-23 | 2022-12-06 | 北京数慧时空信息技术有限公司 | 基于组合迁移学习的遥感样本标注方法 |
CN115563536A (zh) * | 2022-09-28 | 2023-01-03 | 北京工业大学 | 基于子域自适应的滚动轴承故障诊断方法 |
CN115687925A (zh) * | 2022-10-31 | 2023-02-03 | 中国船舶重工集团公司第七一九研究所 | 针对不平衡样本的故障类型识别方法及装置 |
CN115563494A (zh) * | 2022-12-01 | 2023-01-03 | 北京科技大学 | 基于迁移学习的多模式流程工业过程的数据驱动故障检测方法 |
Non-Patent Citations (4)
Title |
---|
A New Deep Transfer Learning Method for Bearing Fault Diagnosis Under Different Working Conditions;Jun Zhu 等;《IEEE SENSORS JOURNAL》;8394-8402 * |
Domain Adaptation-Based Transfer Learning for Gear Fault Diagnosis Under Varying Working Conditions;Chao Chen 等;《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》;第70卷卷;1-10 * |
数据驱动的污水处理过程监测与故障诊断研究;程洪超;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;B027-183 * |
数据驱动的燃煤发电装备运行工况监控 —— 现状与展望;赵春晖 等;《自动化学报》;第第48卷卷(第第11期期);2611-2633 * |
Also Published As
Publication number | Publication date |
---|---|
CN116206154A (zh) | 2023-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111142501B (zh) | 基于半监督自回归动态隐变量模型的故障检测方法 | |
CN109766950B (zh) | 一种基于形态聚类和LightGBM的工业用户短期负荷预测方法 | |
CN111695456B (zh) | 一种基于主动判别性跨域对齐的低分辨人脸识别方法 | |
CN110728694B (zh) | 一种基于持续学习的长时视觉目标跟踪方法 | |
CN109492750B (zh) | 基于卷积神经网络和因素空间的零样本图像分类方法 | |
CN114818996B (zh) | 基于联邦域泛化的机械故障诊断方法及*** | |
CN107403618B (zh) | 基于堆叠基稀疏表示的音频事件分类方法及计算机设备 | |
CN114722892A (zh) | 基于机器学习的持续学习方法及装置 | |
CN115393631A (zh) | 基于贝叶斯层图卷积神经网络的高光谱图像分类方法 | |
CN115563494A (zh) | 基于迁移学习的多模式流程工业过程的数据驱动故障检测方法 | |
CN114358157A (zh) | 一种基于时序数据LSTM特征的K-Shape聚类方法 | |
CN113988415A (zh) | 一种中长期电力负荷预测方法 | |
CN117611838A (zh) | 一种基于自适应超图卷积网络的多标签图像分类方法 | |
CN115705393A (zh) | 一种基于持续学习的雷达辐射源分级识别方法 | |
CN116206154B (zh) | 一种非常见工况下的污水处理过程的故障分类方法 | |
CN116527346A (zh) | 基于深度学习图神经网络理论的威胁节点感知方法 | |
CN114036947B (zh) | 一种半监督学习的小样本文本分类方法和*** | |
CN111401519B (zh) | 一种基于物体内和物体间相似性距离的深层神经网络无监督学习方法 | |
Mirzaeifard et al. | Dynamic graph topology learning with non-convex penalties | |
CN113296947A (zh) | 基于改进XGBoost模型的资源需求预测方法 | |
Li et al. | A BYY scale-incremental EM algorithm for Gaussian mixture learning | |
Xu-Dong et al. | Control Chart Recognition Method Based on Transfer Learning | |
CN112364892B (zh) | 一种基于动态模型的图像识别方法及装置 | |
CN117878933B (zh) | 一种基于突变气象环境识别的风电功率预测方法 | |
CN111832427B (zh) | 基于欧氏对齐和Procrustes分析的EEG分类的迁移学习方法和*** |
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 |