[发明专利]一种具有可调整能力的多工况混合变量异常检测方法在审
申请号: | 202110889880.2 | 申请日: | 2021-08-04 |
公开(公告)号: | CN113505496A | 公开(公告)日: | 2021-10-15 |
发明(设计)人: | 周东华;王敏;陈茂银;纪洪泉;王建东;钟麦英 | 申请(专利权)人: | 山东科技大学 |
主分类号: | G06F30/20 | 分类号: | G06F30/20;G06F17/18;G06K9/62 |
代理公司: | 青岛智地领创专利代理有限公司 37252 | 代理人: | 肖峰 |
地址: | 266590 山东*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 具有 可调整 能力 工况 混合 变量 异常 检测 方法 | ||
1.一种具有可调整能力的多工况混合变量异常检测方法,其特征在于,包括如下步骤:
步骤1:离线建模,具体包括如下步骤:
步骤1.1:采样历史工况构建训练数据集共包含N个样本,为第i个采样并包含d个样本,li为xi对应的标签,i表示采样时间;记Xk=[Xkn,Xkf]为工况k的数据,其中和分别为工况k的正常和故障样本;总共包含K1个工况,记xj为第j个变量,二值变量和连续变量的集合分别记为和
步骤1.2:当xj∈Xc时,假设其服从高斯分布
Pc(xj|l=γ;θj,γ)=N(xj|l=γ;uj,γ,Σj,γ) (1)
其中,uj,γ,Σj,γ为对应的均值和标准差;
步骤1.3:当xj∈Xb时,假设其服从伯努利分布
其中,δj,γ为第j个变量在γ下的响应函数;
步骤1.4:计算条件概率
其中,为标签变量,γ∈{1n,1f,...,K1n,Kfn},参数向量θ={θj,γ,Wj,k},θj,γ为第j个变量在类别γ下的参数,Wj,k为第j个变量在工况k下的权值;
步骤1.5:如果xj={x1j,...,xij,...,xnj}∈Xc,构造二值变量x′j={x′1j,...,x′ij,...,x′nj}
其中所有关于权值的计算都用构造的二值变量代替原始的连续变量进行计算;
步骤1.6:计算条件概率
其中,φ(xj)为示性函数(如果xj=1,φ(xj)=1;否则φ(xj)=0),步骤1.7:计算条件联合概率
P(xj,xj′|k)=P(xj′|k){1-φ(xj)+(2φ(xj)-1)[φ(xj′)η+(1-φ(xj′))η′]} (6)
其中,
步骤1.8:计算条件互信息
步骤1.9:如果xj,xj′∈Xc,则按照下式进行修正
步骤1.10:计算相关系数
其中,M(xj,l|k)为条件互信息;
步骤1.11:计算特征权值
其中,ξj,k为相关系数;
步骤1.12:当xj∈Xb时,计算响应函数
其中,如果li=γ,Ξi,γ=1,否则Ξi,γ=0;
步骤1.13:当xj∈Xc时,计算均值
其中
步骤1.14:当xj∈Xc时,计算标准差
步骤1.15:计算先验概率
步骤1.16:计算条件概率
步骤1.17:新来采样的标签为
其中,P(l=γ)为先验概率;
步骤2:在线检测和更新,具体包括如下步骤:
步骤2.1:当有新采样到来时,计算(16)式判断是否发生故障;
步骤2.3:当有新工况到来时,更新先验概率,补充新到来工况的对应参数;
步骤2.3:当有工况消失时,更新先验概率,删除模型中对应消失工况的参数;
若:预测标签与正常数据标签相同则认定为正常,预测标签与故障数据标签一致则认为发生故障。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于山东科技大学,未经山东科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110889880.2/1.html,转载请声明来源钻瓜专利网。