CN110135121B - 一种基于拉格朗日-粒子群更新算法的故障诊断方法 - Google Patents
一种基于拉格朗日-粒子群更新算法的故障诊断方法 Download PDFInfo
- Publication number
- CN110135121B CN110135121B CN201910512377.8A CN201910512377A CN110135121B CN 110135121 B CN110135121 B CN 110135121B CN 201910512377 A CN201910512377 A CN 201910512377A CN 110135121 B CN110135121 B CN 110135121B
- Authority
- CN
- China
- Prior art keywords
- fault
- particle
- model
- lagrange
- algorithm
- 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
Images
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于拉格朗日‑粒子群更新算法的故障诊断方法。包括:建立故障诊断数学模型、对模型进行拉格朗日松弛、用粒子群更新方法寻找全局最优粒子、将最优粒子代入模型求解诊断结果四个步骤。本方法使用粒子群算法求解出最优的拉格朗日乘子,构造以上下界之差的绝对值为主要参数的适应度函数,基于该适应度函数选取全局最优粒子。为了满足原拉格朗日松弛算法中乘子非负的条件,构造指数函数减缓粒子往负方向的移动速度,并将种群中的粒子限定非负。本发明能够发现系统中的隐含故障,且故障隔离率更高,算法的结果更加准确,提高了算法的精度。
Description
技术领域
本发明涉及一种基于拉格朗日-粒子群更新算法的故障诊断方法,属于故障诊断领域。
技术背景
基于模型和测试结果的故障诊断是目前应用较为广泛的故障诊断方法。该方法基于系统的相关性矩阵,采用贝叶斯理论建立其故障诊断模型并使用人工智能算法进行求解。
拉格朗日松弛算法属于人工智能算法,已成熟应用于求解基于贝叶斯理论建立的故障诊断模型。但是现有的拉格朗日松弛算法理论分析较为复杂,算法涉及最优解与拉格朗日乘子相互迭代的过程,结构层次不分明,算法中判断和循环较多,编程实现难度较大。
目前,尚未发现一种利用粒子群方法更新拉格朗日乘子,对拉格朗日松弛算法进行改进的算法。
发明内容
本发明的目的是采用粒子群算法对拉格朗日松弛算法进行改进,提出一种基于拉格朗日-粒子群更新算法的故障诊断方法。首先,基于故障-测试相关性矩阵,利用贝叶斯理论构造故障诊断模型;其次,对故障诊断模型进行拉格朗日松弛,得到松弛后的模型;最后,采用粒子群算法更新拉格朗日乘子找到其全局最优解,利用该乘子求解拉格朗日模型得到诊断结果。
一种基于拉格朗日-粒子群更新算法的故障诊断方法包括以下4个步骤:
步骤1:利用相关性矩阵构造故障诊断模型。
根据系统结构和信号流向,利用图论分析得到故障-测试相关性矩阵:
其中,fi(i=1,…m)表示第i个故障;tj(j=1,…n)表示第j个测试;相关性矩阵D中元素dij取值为1或0,分别表示第i个故障与第j个测试相关或不相关。
基于该相关性矩阵,利用贝叶斯后验概率公式,将需要求解的最大后验故障概率转化为系统已知的先验故障概率,得到故障诊断模型:
其中,S-为故障集S中去除确认未发生故障的故障集合,A为D矩阵中除去未发生故障的测试列的矩阵再转置得到;pi(i=1,…,m)为故障fi的先验概率;x为故障向量,表示组件的故障状态,当故障fi发生故障时,xi=1;反之xi=0。
步骤2:将故障诊断模型进行拉格朗日松弛,得到松弛后的优化模型。
将约束条件Ax≥e吸收到目标函数ZIP中,得到进行拉格朗日松弛后的模型为:
步骤3:利用粒子群更新方法找到全局最优的拉格朗日乘子λb。
流程图如图2所示,包括以下步骤:
步骤3.1:根据拉格朗日乘子λ的维数和种群数量,构造初始的拉格朗日乘子种群;
步骤3.2:计算种群每个粒子λ的上下界之差的绝对值J;
具体方法为:将每个粒子代入到模型(2)中,求解出下界ZDW和故障x,再根据故障x构造可行解x′,将可行解x′代入模型(1)求解得到上界ZUP,计算上下界之差的绝对值J=|ZUP-ZDW|。
步骤3.4:当本代粒子的适应度函数值Ft大于粒子在历史最佳位置时的适应度函数值Fb时,即Ft>Fb时,更新粒子的历史最佳位置Pb=Pt,否则,不更新历史最佳位置;初始历史最佳位置即种群的初始值。
步骤3.6:惯性权重更新wt=wmax-t(wmax-wmin)/Nmax,其中,t为迭代代数;wmax、wmin为惯性权重的最大值和最小值。
步骤3.7:更新粒子的速度和位置,第i个粒子在第t+1次迭代中第d维的速度为其中,c1、c2为加速度因子;第t代生成的在0-1之间的随机数。若更新后的速度为负,令以减缓粒子向负方向移动的速度(因拉格朗日乘子λ要求非负)。其中,K1为比例系数;vmax为设定的速度最大值。更新后的位置为若令
步骤4:利用得到的最优乘子λb代入到松弛后的模型(2)中,求解得到诊断结果xb。
以上为基于拉格朗日-粒子群更新算法的故障诊断方法具体实施步骤。
本发明的优点在于:
(1)本发明能根据适应度函数值从种群中选出最优粒子,粒子的可选范围更广,选取出的全局最优粒子更好,相比于拉格朗日松弛算法求解出的诊断结果精度更高。
(2)本发明能发现隐含故障,排除所有可能存在故障的组件,提高了故障诊断的安全性和可靠性。
(3)本发明不涉及拉格朗日乘子与故障相互迭代的过程,而是先求解最优拉格朗日乘子,再代入模型求解诊断结果。并且改进后算法只需一次判断,而拉格朗日松弛算法需要三次判断,结构更加复杂。改进后的算法理论分析简单,结构清晰明了,编程易于实现。
附图说明
图1是基于拉格朗日-粒子群更新算法的故障诊断流程图。
图2是拉格朗日-粒子群更新算法的流程图。
图3是拉格朗日松弛算法与拉格朗日-粒子群更新算法的对比图。
图4是实例中故障为1,5时全局最优粒子适应度函数值随迭代次数的变化图像。
具体实施方式
以下结合具体的实施例对基于拉格朗日-粒子群更新算法的故障诊断方法作进一步的说明:
阿波罗飞船发射前系统相关性矩阵如表1所示。
表1阿波罗飞船系统相关性矩阵
选取系统中故障的先验概率相等,均为0.1。若系统故障为1,5时,则检测到故障的测试有2,4,6,7,8,9,10,11,12,13,14。算法应当能够根据测试点未通过(检测到故障)的情况求解出故障1,5。
步骤1:利用相关性矩阵构造故障诊断模型。
构造的故障诊断模型为:
表2化简后得到的的矩阵A
步骤2:将故障诊断模型进行拉格朗日松弛,得到松弛后的优化模型。
进行拉格朗日松弛后的模型为:
步骤3:利用粒子群更新方法找到全局最优的拉格朗日乘子λb。
步骤3.1:构造初始的拉格朗日乘子种群,种群数量选取50,具体编程时构造50行11列的矩阵,每一行相当于一个粒子,初始种群矩阵元素均为0-1之间的随机数;
步骤3.2:计算种群每个粒子λ的上下界之差的绝对值J;
具体方法为:将每个粒子代入到模型(2)中,求解出下界ZDW(即模型求解出的函数值)和故障向量x。构造可行解向量x′,初始值与故障向量x相等,然后,对故障向量x中不满足约束条件Ax≥e的行(即),将对应的行添加到解向量x′中,最后,对解向量x′中的每一个故障进行判断,如果将该故障去除后仍满足约束条件,则将故障去除。将得到的可行解向量x′代入模型(1)求解得到上界ZUP,计算上下界之差的绝对值J=|ZUP-ZDW|。
步骤3.4:根据粒子的适应度函数值决定是否更新粒子的位置,初始历史最佳位置即种群的初始值。
具体方法为:当本代粒子的适应度函数值Ft大于粒子在历史最佳位置时的适应度函数值Fb时,即Ft>Fb时,更新粒子的历史最佳位置Pb=Pt,否则,不更新历史最佳位置。
步骤3.6:惯性权重更新wt=1.2-0.004t,t为迭代代数。
步骤3.7:更新粒子的速度和位置,第i个粒子在第t+1次迭代中第d维(这里速度的与粒子的维数均为11)的速度为其中,第t代生成的在0-1之间的随机数。若更新后的速度为负,令v=-0.2×e-(v+2),减缓向负方向移动的速度(因拉格朗日乘子λ要求非负)。更新后的位置为若令
最后求解得到的全局最优粒子λb如表3所示。
表3故障为1,5时的全局最优粒子λb
故障状态 | <![CDATA[全局最优粒子值λ<sub>b</sub>]]> | 适应度函数值 |
1,5 | 5.6156,0,0.8886,2.9666,0,0,0,4.8169,0.6259,0,0.9080 | 1.0000 |
全局最优粒子在故障状态为1,5时的适应度函数值随迭代次数的变化图像如图4所示。
步骤4:利用得到的最优乘子λb代入到松弛后的优化模型(2)中,求解得到诊断结果向量xb如表4所示。
表4故障为1,5时的故障向量xb
故障状态 | 报故测试的序号 | <![CDATA[x<sub>b</sub>对应的故障状态]]> |
1,5 | 2,4,6,7,8,9,10,11,12,13,14 | 1,5 |
求解出的故障状态与设定的故障状态相同,说明算法能有效地解决故障诊断问题。将该算法应用于其余双故障状态,求解得到的结果如表5所示。
从双故障的诊断结果来看,拉格朗日松弛算法的故障隔离率为88.89%,改进后的算法为91.11%,高于前者的精度。并且当故障为4,6和5,6时,测试的结果相同,改进后的算法能够将故障4,5,6都找出,说明改进后的算法能够发现隐含故障,排除所有故障。
表2拉格朗日-粒子群更新算法双故障诊断结果
Claims (1)
1.一种基于拉格朗日-粒子群更新算法的故障诊断方法,其特征在于:基于系统相关性矩阵和测试结果,利用贝叶斯后验概率公式构造故障诊断模型;将模型中的约束条件吸收到目标函数中,得到拉格朗日松弛后的模型;利用粒子群更新方法找到全局最优拉格朗日乘子,并将该乘子代入到松弛后的优化模型中,求解得到故障诊断结果;
其中,利用相关性矩阵构造故障诊断模型,根据系统结构和信号流向,利用图论分析得到故障-测试相关性矩阵:
其中,fi(i=1,…m)表示第i个故障;tj(j=1,…n)表示第j个测试;相关性矩阵D中元素dij取值为1或0,分别表示第i个故障与第j个测试相关或不相关;
基于该相关性矩阵,利用贝叶斯后验概率公式,将需要求解的最大后验故障概率转化为系统已知的先验故障概率,得到故障诊断模型:
其中,S-为故障集S中去除确认未发生故障的故障集合,A为D矩阵中除去未发生故障的测试列的矩阵再转置得到;pi(i=1,…,m)为故障fi的先验概率;x为故障向量,表示组件的故障状态,当故障fi发生故障时,xi=1;反之xi=0;
其中,将故障诊断模型进行拉格朗日松弛,得到松弛后的优化模型,将约束条件Ax≥e吸收到目标函数ZIP中,得到进行拉格朗日松弛后的模型为:
式中,λi≥0,(i=1,…,|Tf|)为拉格朗日乘子;Tf为测试未通过的测试列;
其中,利用粒子群更新方法找到全局最优的拉格朗日乘子λb,具体步骤为:
步骤1,根据拉格朗日乘子λ的维数和种群数量,构造初始的拉格朗日乘子种群;
步骤2,计算种群每个粒子λ的上下界之差的绝对值J;具体的,将每个粒子代入到模型(2)中,求解出下界ZDW和故障x,再根据故障x构造可行解x′,将可行解x′代入模型(1)求解得到上界ZUP,计算上下界之差的绝对值J=|ZUP-ZDW|;
步骤3,计算每个粒子的适应度函数值其中,K为比例系数,J为粒子的上下界之差的绝对值;
步骤4,当本代粒子的适应度函数值Ft大于粒子历史最佳位置时的适应度函数值Fb时,即Ft>Fb,更新粒子的历史最佳位置Pb=Pt,否则,不更新历史最佳位置;初始历史最佳位置即种群的初始值;
步骤5,根据所有粒子在历史最佳位置时的函数值,找出其中函数值最大的粒子对应的历史最佳位置,即为本代的全局最优粒子
步骤6,惯性权重更新wt=wmax-t(wmax-wmin)/Nmax,其中,t为迭代代数;wmax、wmin为惯性权重的最大值和最小值;
步骤7,更新粒子的速度和位置,第i个粒子在第t+1次迭代中第d维的速度为其中,c1、c2为加速度因子;r1 t、第t代生成的在0-1之间的随机数;若更新后的速度为负,令以减缓粒子向负方向移动的速度;其中,K1为比例系数;vmax为设定的速度最大值;更新后的位置为若令
步骤8,判断t是否达到最大迭代次数Nmax,即t≥Nmax,或者本代的全局最优粒子对应的适应度函数值是否大于等于设定的阈值;若是,执行步骤4;若不是,迭代次数t=t+1,转步骤2;
步骤9,利用得到的最优乘子λb代入到松弛后的模型(2)中,求解得到诊断结果xb。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910512377.8A CN110135121B (zh) | 2019-06-13 | 2019-06-13 | 一种基于拉格朗日-粒子群更新算法的故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910512377.8A CN110135121B (zh) | 2019-06-13 | 2019-06-13 | 一种基于拉格朗日-粒子群更新算法的故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110135121A CN110135121A (zh) | 2019-08-16 |
CN110135121B true CN110135121B (zh) | 2023-04-18 |
Family
ID=67577128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910512377.8A Active CN110135121B (zh) | 2019-06-13 | 2019-06-13 | 一种基于拉格朗日-粒子群更新算法的故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110135121B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007522658A (ja) * | 2004-02-06 | 2007-08-09 | テスト アドバンテージ, インコーポレイテッド | データ解析の方法および装置 |
CN102200759A (zh) * | 2011-05-28 | 2011-09-28 | 东华大学 | 一种非线性核化自适应预测方法 |
CN102426689A (zh) * | 2011-08-23 | 2012-04-25 | 浙江大学 | 基于离散粒子群和拉格朗日乘数算法的机组组合优化方法 |
CN102609593A (zh) * | 2012-03-05 | 2012-07-25 | 浙江大学 | 基于多重先验知识混合模型的聚丙烯熔融指数预报方法 |
CN204027820U (zh) * | 2014-08-26 | 2014-12-17 | 安徽工程大学 | 汽车空调故障检测系统 |
CN104697765A (zh) * | 2014-08-26 | 2015-06-10 | 安徽工程大学 | 汽车空调故障检测系统及方法 |
CN109470477A (zh) * | 2018-09-27 | 2019-03-15 | 昆明理工大学 | 一种基于改进的pso算法优化fsvm的滚动轴承故障诊断方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140257767A1 (en) * | 2013-03-09 | 2014-09-11 | Bigwood Technology, Inc. | PSO-Guided Trust-Tech Methods for Global Unconstrained Optimization |
CN105093095B (zh) * | 2015-07-10 | 2017-12-12 | 四川航天系统工程研究所 | 模拟非线性直流电路的改进粒子群优化诊断方法及系统 |
CN106443297B (zh) * | 2016-09-07 | 2018-05-22 | 江南大学 | 光伏二极管箝位型三电平逆变器的决策树svm故障诊断方法 |
US11355240B2 (en) * | 2017-09-26 | 2022-06-07 | Edge2020 LLC | Determination of health sciences recommendations |
CN107884706B (zh) * | 2017-11-09 | 2020-04-07 | 合肥工业大学 | 基于向量值正则核函数逼近的模拟电路故障诊断方法 |
CN108414923A (zh) * | 2018-02-05 | 2018-08-17 | 武汉大学 | 一种基于深度置信网络特征提取的模拟电路故障诊断方法 |
-
2019
- 2019-06-13 CN CN201910512377.8A patent/CN110135121B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007522658A (ja) * | 2004-02-06 | 2007-08-09 | テスト アドバンテージ, インコーポレイテッド | データ解析の方法および装置 |
CN102200759A (zh) * | 2011-05-28 | 2011-09-28 | 东华大学 | 一种非线性核化自适应预测方法 |
CN102426689A (zh) * | 2011-08-23 | 2012-04-25 | 浙江大学 | 基于离散粒子群和拉格朗日乘数算法的机组组合优化方法 |
CN102609593A (zh) * | 2012-03-05 | 2012-07-25 | 浙江大学 | 基于多重先验知识混合模型的聚丙烯熔融指数预报方法 |
CN204027820U (zh) * | 2014-08-26 | 2014-12-17 | 安徽工程大学 | 汽车空调故障检测系统 |
CN104697765A (zh) * | 2014-08-26 | 2015-06-10 | 安徽工程大学 | 汽车空调故障检测系统及方法 |
CN109470477A (zh) * | 2018-09-27 | 2019-03-15 | 昆明理工大学 | 一种基于改进的pso算法优化fsvm的滚动轴承故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
Xiang Yu.Unit commitment using Lagrangian relaxtion and particle swarm optimization.《International Journal or Electrical Power & Energy Systems》.2014,第61卷第510-522页. * |
张雷.基于PSO-LSSVM的牵引供电系统故障诊断研究.《机车点传动》.2019,(第期),第51-55页. * |
王林军.基于粒子群算法和增广拉格朗日乘子法的混合可靠性分析.三峡大学学报.2019,第41卷(第41期),第108-112页. * |
Also Published As
Publication number | Publication date |
---|---|
CN110135121A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112101530B (zh) | 神经网络训练方法、装置、设备及存储介质 | |
CN110888059B (zh) | 基于改进随机森林联合容积卡尔曼的荷电状态估计算法 | |
CN105825045B (zh) | 多阶段任务系统可修备件需求预测方法 | |
CN108363876B (zh) | 一种考虑关键故障的测试优化选择方法 | |
Zhao et al. | Model-based fault diagnosis methods for systems with stochastic process–a survey | |
CN107909194B (zh) | 系统级测试性设计多目标优化方法 | |
CN110929835B (zh) | 一种新型碳化硅基航空功率变换器故障诊断方法及系统 | |
CN114626307B (zh) | 一种基于变分贝叶斯的分布式一致性目标状态估计方法 | |
CN111339487B (zh) | 基于径向基函数rbf的边坡系统失效概率计算方法 | |
KR101337499B1 (ko) | 신뢰성 기반 최적설계 방법 | |
CN107657311B (zh) | 基于多目标粒子群算法的测试优选方法 | |
CN110457776B (zh) | 一种基于故障决策网络的测试策略快速生成方法 | |
CN114567288B (zh) | 基于变分贝叶斯的分布协同非线性系统状态估计方法 | |
CN110135121B (zh) | 一种基于拉格朗日-粒子群更新算法的故障诊断方法 | |
Aghapour et al. | Performance specified state estimation with minimum risk | |
Liu et al. | Predictive state representations with state space partitioning | |
CN112488295A (zh) | 交叉验证算法优化lstm网络继电器贮存寿命预测方法 | |
CN108090566B (zh) | 基于串并联遗传算法的多目标测试优选方法 | |
CN111412795A (zh) | 测试点设置方案生成方法及装置 | |
CN116595865A (zh) | 一种基于维纳过程的设备剩余寿命预测方法 | |
Iglésis et al. | Simultaneous actuator and sensor faults estimation for aircraft using a Jump-Markov regularized particle filter | |
Qu et al. | Fault tolerant control system based on subspace predictive control and multiple model predictive control | |
Smith et al. | Recurrent neural network ensembles for convergence prediction in surrogate-assisted evolutionary optimization | |
CN110147306A (zh) | 一种考虑关联失效的容错软件可靠性及性能评估方法 | |
Jiang et al. | Robust linear-complexity approach to full SLAM problems: Stochastic variational Bayes inference |
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 |