CN107169243A - 一种基于hsmm和经验模型的燃料电池故障预测方法 - Google Patents

一种基于hsmm和经验模型的燃料电池故障预测方法 Download PDF

Info

Publication number
CN107169243A
CN107169243A CN201710499198.6A CN201710499198A CN107169243A CN 107169243 A CN107169243 A CN 107169243A CN 201710499198 A CN201710499198 A CN 201710499198A CN 107169243 A CN107169243 A CN 107169243A
Authority
CN
China
Prior art keywords
mrow
msub
hsmm
data
fuel cell
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.)
Granted
Application number
CN201710499198.6A
Other languages
English (en)
Other versions
CN107169243B (zh
Inventor
吴小娟
叶倩文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710499198.6A priority Critical patent/CN107169243B/zh
Publication of CN107169243A publication Critical patent/CN107169243A/zh
Application granted granted Critical
Publication of CN107169243B publication Critical patent/CN107169243B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Abstract

本发明公开了一种基于HSMM和经验模型的燃料电池故障预测方法,先采集一组燃料电池的全寿命电压数据,作为训练数据,再通过Welch‑Baum算法利用训练数据训练HSMM模型;再采集一组燃料电池的电压退化数据,作为测试数据;将测试数据输入HSMM中,通过前向算法估计当前健康状态,并根据状态持续时间计算其剩余寿命;根据测试数据建立经验模型,并估计参数及预测未来走势,根据未来电压走势计算其剩余寿命;以训练数据和测试数据的电压梯度值的相似度作为标准,将HSMM得到的剩余寿命和经验模型计算出的剩余寿命结合起来,估算出燃料电池的最终剩余寿命。

Description

一种基于HSMM和经验模型的燃料电池故障预测方法
技术领域
本发明属于燃料电池技术领域,更为具体地讲,涉及一种基于HSMM和经验模型的燃料电池故障预测方法。
背景技术
燃料电池是一个通过电化学反应把燃料中化学能量直接转化为电能的装置,它具有环境友好、能量转换率高等优点,然而,可靠性和持久性是阻碍燃料电池发展的一个主要原因。
故障预测技术可以在燃料电池故障早期及时的诊断它的健康状态、预测其剩余可用寿命,使得我们可以及时的对它进行维修,以此延长燃料电池寿命。近几十年来,故障预测技术在解决燃料电池的寿命相对较短的问题上表现出了很大的潜力。
基于隐式半马尔科夫模型(HSMM)的故障预测方法是一种基于数据的预测方法,不需要对燃料电池的内部结构和物理退化过程进行建模,直接对系统的退化数据进行分析即可预测其剩余寿命。然而当燃料电池的电压退化趋势与HSMM全寿命训练数据的电压退化趋势相似度较低时,基于HSMM的预测结果较差。
为了解决这个问题,本发明提出了将HSMM与经验模型相结合的混合故障预测方法,若燃料电池的当前运行电压退化趋势与HSMM训练电压退化趋势相似度较高,则利用HSMM对燃料电池的剩余寿命进行预测;若两者的相似度较低,则利用经验模型估计燃料电池剩余运行寿命。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于HSMM和经验模型的燃料电池故障预测方法,根据燃料电池的退化数据,对燃料电池的剩余寿命进行预测。
为实现上述发明目的,本发明一种基于HSMM和经验模型的燃料电池故障预测方法,其特征在于,包括以下步骤:
(1)、设置初始状态概率矩阵π0、初始状态转移概率矩阵A0、初始观测值概率矩阵B0、初始状态持续时间概率矩阵P0,以及隐含状态数N和观测值个数M,构建HSMM模型λ0=(N,M,A0,B00,P0);
(2)、采集一组燃料电池的全寿命电压数据作为训练数据,利用Welch-Baum算法更新HSMM模型中A0、B0、π0和P0的值,得到更新后的HSMM模型λ=(N,M,A,B,π,P);
(3)、采集当前燃料电池的一段电压数据作为测试数据,将测试数据输入至训练好的HSMM模型中,利用前向算法计算当系统处于N种不同隐含状态时的P(O|λk),k=1,2,…,N;其中,λ1=(N,M,A,B,π1,P),λ2=(N,M,A,B,π2,P),……,λN=(N,M,A,B,πN,P);
αT(i)=P(o1,o2…oT,St=i|St+1≠i,λk)
其中,t=1,2,…,T,表示T个观测时刻,O=(o1,o2,o3,…oT)表示从初始时刻t=1到T时刻的测试数据,St=i表示在t时刻时系统处于i状态;
(4)、从N个不同的P(O|λk)中选出最大值P(O|λm);
P(O|λm)=max(P(O|λ1),P(O|λ2),…,P(O|λN)),1≤m≤N
将P(O|λm)所对应的HSMM模型标记为当前燃料电池的模型λm=(N,M,A,B,πm,P),当前燃料电池的健康状态处于隐含状态m;
(5)、从训练数据中提取处于隐含状态m且与测试数据相同长度的观测数据On,n表示观测数据组数,将观测数据On输入至λm=(N,M,A,B,πm,P),通过前向后向算法分别计算出n个不同的P(Ojm),j=1,2,…,n;将P(O|λm)分别与P(Ojm)进行一一比较,选择大小最相近的一组观测数据,以该组观测数据离状态转移点的时间作为当前测试数据距离状态转移点的时间ttransfer
(6)、利用HSMM算法,计算燃料电池的剩余寿命Rulhsmm
其中,a表示从状态m转移到状态μ的概率;
(7)、根据测试数据建立经验模型:
其中,a,b,c,d为待估计参数,可以通过梯度下降法估算;
(8)、利用经验模型,计算燃料电池剩余寿命Rulempirial
Rulempirial=arg(Vfailure(t))-tcurrent
其中,arg(Vfailure(t))表示经验模型到预设电压阈值的时间,tcurrent为测试数据的当前时间;
(9)、计算训练数据和测试数据的电压梯度值gradζ(t):
gradζ(t)=Vζ(t)-Vζ(t-1)
其中,ζ=1,2,V1(t)和grad1(t)分别表示t时刻训练数据的电压值和电压梯度值;V2(t)和grad2(t)分别表示t时刻测试数据的电压值和电压梯度值;
(10)、计算指标β:
其中,tcurrent表示当前时间;Tdef为预设值,用来决定用几个单位的电压梯度值来计算相似度;
(11)、估算燃料电池的剩余寿命Rulfinal
Rulfinal=α·Rulhsmm+(1-α)·Rulempirial
本发明的发明目的是这样实现的:
本发明基于HSMM和经验模型的燃料电池故障预测方法,先采集一组燃料电池的全寿命电压数据,作为训练数据,再通过Welch-Baum算法利用训练数据训练HSMM模型;再采集一组燃料电池的电压退化数据,作为测试数据;将测试数据输入HSMM中,通过前向算法估计当前健康状态,并根据状态持续时间计算其剩余寿命;根据测试数据建立经验模型,并估计参数及预测未来走势,根据未来电压走势计算其剩余寿命;以训练数据和测试数据的电压梯度值的相似度作为标准,将HSMM得到的剩余寿命和经验模型计算出的剩余寿命结合起来,估算出燃料电池的最终剩余寿命。
同时,本发明基于HSMM和经验模型的燃料电池故障预测方法还具有以下有益效果:
(1)、故障预测所需数据仅为燃料电池的输出电压数据,对燃料电池内部的数据没有要求,因此在采集数据方面成本较低,较方便,也不会在采集过程中对燃料电池的内部结构造成影响。
(2)、该方法相较于基于经验模型的故障预测方法,在预测速度上有很大的提高。
(3)、本发明相较于目前提出的其它燃料电池故障预测方法,有较高的预测精度。
附图说明
图1是本发明基于HSMM和经验模型的燃料电池故障预测方法流程图;
图2是HSMM训练数据状态划分结果;
图3是HSMM健康状态诊断结果;
图4是基于HSMM剩余寿命预测结果;
图5是经验模型各个时刻预测的电压值;
图6是经验模型预测的剩余寿命
图7是基于HSMM和经验模型的燃料电池的预测剩余寿命。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明基于HSMM和经验模型的燃料电池故障预测方法流程图。
在本实施例中,如图1所示,本发明一种基于HSMM和经验模型的燃料电池故障预测方法,包括以下步骤:
S1、随机设置初始状态概率矩阵π0、初始状态转移概率矩阵A0、初始观测值概率矩阵B0、初始状态持续时间概率矩阵P0,以及隐含状态数N为3和观测值个数M为10,构建HSMM模型λ0=(N,M,A0,B00,P0)。
PAI0:[3×1double]
A0:[3×3double]
B0:[3×10double]
P0:[3×300double]
S2、采集一组燃料电池的全寿命电压数据作为训练数据,利用Welch-Baum算法更新HSMM模型中A0、B0、π0和P0的值,得到更新后的HSMM模型λ=(N,M,A,B,π,P);
Welch-Baum算法是极大似然估计的一个应用,它采用了一种多次迭代优化的算法,用已知的训练数据和HSMM模型构造函数,该函数中包括了HSMM中所有的六个参数N,M,A0,B00,P0,然后令各偏导为零,得到新的HSMM各参数的估计值,用新旧HSMM模型参数之间的函数关系反复进行迭代运算,直到HSMM的各个参数估计都稳定不再发生明显变化为止,得到更新后的HSMM模型λ=(N,M,A,B,π,P)。其中更新后的模型中,
另外,还有一个3*10的矩阵B和一个3*300的矩阵P,由于矩阵过大在此处就没有展示。
根据状态持续时间概率矩阵P,利用公式:
计算得到各个状态的平均持续时间为:
根据训练结果,可作出如图2所示的状态划分图。
S3、沿用上面的训练结果,采集当前燃料电池的一段电压数据作为测试数据,将测试数据输入至3个不同的HSMM模型:λ1=(N,M,A,B,π1,P),λ2=(N,M,A,B,π2,P),λ3=(N,M,A,B,π3,P),HSMM模型通过前向算法分别计算出3个不同的P(O|λk),k=1,2,3,如表1所示。
表1 是健康状态诊断时各样本的P(O|λ)值。
正常状态 退化状态1 退化状态2
样本1 -Inf -Inf -0.2344
样本2 -0.1860 -5.3669 -5.2402
样本3 -Inf -0.1846 -5.5593
样本4 -Inf -Inf -0.2511
样本5 -Inf -0.1503 -4.8085
样本6 -Inf -Inf -0.1732
样本7 -0.0700 -2.6605 -2.5971
样本8 -Inf -0.1735 -5.5482
样本9 -Inf -Inf -0.1876
样本10 -Inf -0.1120 -4.4118
样本11 -Inf -Inf -0.1375
样本12 -0.1001 -3.8007 -3.7101
样本13 -Inf -0.0560 -2.2059
样本14 -Inf -Inf -0.0952
表1
S4、以第二个样本为例,P(O|λ1)为-0.1860,P(O|λ2)也为-5.3669,P(O|λ3)为-5.2402。因P(O|λ1)的值最大,我们认为此时的燃料电池系统处于第一个隐含状态。
在本实施例中,我们取了14组数据来做当前健康状态的估计,估计结果如图3所示,从图3中可以看出,对于实验所选取的14个样本点,HSMM都能精确的估计其健康状态。
S5、以第二个样本为例,该样本处于第一个隐含状态,根据上述发明内容中第五步描述的方法,计算得到的ttransfer为20。
S6、利用HSMM算法,计算燃料电池的剩余寿命Rulhsmm。以第二个样本为例,该样本处于第一个隐含状态,根据上述式子,可计算得到:
Rul3=72h
Rul2=44h+1*Rul3=116h
Rulhsmm=20h+1*Rul2+0*Rul3=136h
根据该方法计算所得的RULhsmm如图4所示。图中,横坐标表示的是进行寿命预测时的时刻;纵坐标代表各时刻寿命预测的结果。粗体实线表示的是HSMM对于阳极氯气浓度为的燃料电池寿命预测曲线;粗体虚线为各个时刻所对应的真实寿命值;上侧细虚线为真实寿命值的置信上限;下侧细点线为真实寿命值的置信下限,两者的置信概率均为15%。
S7、在本实施例中,根据测试数据设置50组不同的初始值,限制Voltage(t)是关于时间t的单调递减函数,并且限制在一定时间内电压不会退化到失效,再用梯度下降法分别估计出参数a,b,c,d的值;然后根据估计出的参数值,分别预测燃料电池电压在未来的走势,取这50组不同预测曲线中的中位数所对应的模型,作为我们所求的经验模型。以燃料电池前10个小时运行的电压数据为例,根据梯度下降法训练得到的经验模型模型为:
S8、利用经验模型,计算燃料电池剩余寿命Rulempirial
Rulempirial=arg(Vfailure(t))-tcurrent
其中,arg(Vfailure(t))表示经验模型到预设电压阈值的时间,tcurrent为测试数据的当前时间;以燃料电池前10个小时运行的电压数据为例,我们预设的电压阈值为0.89V,测试数据的当前时间为10小时,根据上述式子计算可得:
Rulempirial=144h-10h=134h
在10、20、……140时刻所选的经验模型电压预测结果如图5所示。经验模型的剩余寿命预测结果如图6所示。
S9、根据上述发明内容中第九步展示的电压梯度计算公式,分别计算HSMM训练数据和当前测试数据的电压梯度值。
S10、同样以燃料电池前10个小时运行的电压数据为例,设置Tdef=10,根据上述发明内容中第十步展示的β计算公式,可得β=0.26。其中,指标β用来衡量训练数据和测试数据之间的相似度,β值越小,则相似度越高;β值越大,相似度越低。
S11、同样以燃料电池前10个小时运行的电压数据为例,根据上述发明内容中第十一步展示的α计算公式及Rulfinal,此时由于β=0.26,那么α=1,因此Rulfinal=Rulhsmm=136h。
根据该结合方法得到的最终剩余寿命结果如图7所示。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

Claims (2)

1.一种基于HSMM和经验模型的燃料电池故障预测方法,其特征在于,包括以下步骤:
(1)、设置初始状态概率矩阵π0、初始状态转移概率矩阵A0、初始观测值概率矩阵B0、初始状态持续时间概率矩阵P0,以及隐含状态数N和观测值个数M,构建HSMM模型λ0=(N,M,A0,B00,P0);
(2)、采集一组燃料电池的全寿命电压数据作为训练数据,利用Welch-Baum算法更新HSMM模型中A0、B0、π0和P0的值,得到更新后的HSMM模型λ=(N,M,A,B,π,P);
(3)、采集当前燃料电池的一段电压数据作为测试数据,将测试数据输入至训练好的HSMM模型中,利用前向算法计算当系统处于N种不同隐含状态时的P(O|λk),k=1,2,…,N;其中,λ1=(N,M,A,B,π1,P),λ2=(N,M,A,B,π2,P),……,λN=(N,M,A,B,πN,P);
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>O</mi> <mo>|</mo> <msub> <mi>&amp;lambda;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>&amp;alpha;</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow>
αT(i)=P(o1,o2…oT,St=i|St+1≠i,λk)
其中,t=1,2,…,T,表示T个观测时刻,O=(o1,o2,o3,…oT)表示从初始时刻t=1到T时刻的测试数据,St=i表示在t时刻时系统处于状态;
(4)、从N个不同的P(O|λk)中选出最大值P(O|λm);
P(O|λm)=max(P(O|λ1),P(O|λ2),…,P(O|λN)),1≤m≤N
将P(O|λm)所对应的HSMM模型标记为当前燃料电池的模型λm=(N,M,A,B,πm,P),当前燃料电池的健康状态处于隐含状态m;
(5)、从训练数据中提取处于隐含状态m且与测试数据相同长度的观测数据On,n表示观测数据组数,将观测数据On输入至λm=(N,M,A,B,πm,P),通过前向后向算法分别计算出n个不同的P(Ojm),j=1,2,…,n;将P(O|λm)分别与P(Ojm)进行一一比较,选择大小最相近的一组观测数据,以该组观测数据离状态转移点的时间作为当前测试数据距离状态转移点的时间ttransfer
(6)、利用HSMM算法,计算燃料电池的剩余寿命Rulhsmm
<mrow> <msub> <mi>Rul</mi> <mrow> <mi>h</mi> <mi>s</mi> <mi>m</mi> <mi>m</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>s</mi> <mi>f</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>&amp;mu;</mi> <mo>=</mo> <mi>m</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>a</mi> <mrow> <mi>m</mi> <mi>&amp;mu;</mi> </mrow> </msub> <msub> <mi>Rul</mi> <mi>&amp;mu;</mi> </msub> </mrow>
<mrow> <msub> <mi>Rul</mi> <mi>&amp;mu;</mi> </msub> <mo>=</mo> <msub> <mi>Dur</mi> <mi>&amp;mu;</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>&amp;delta;</mi> <mo>=</mo> <mi>&amp;mu;</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>a</mi> <mrow> <mi>&amp;mu;</mi> <mi>&amp;delta;</mi> </mrow> </msub> <msub> <mi>Rul</mi> <mi>&amp;delta;</mi> </msub> <mo>,</mo> <mi>&amp;mu;</mi> <mo>&gt;</mo> <mi>m</mi> </mrow>
其中,a表示从状态m转移到状态μ的概率;
(7)、根据测试数据建立经验模型:
<mrow> <mi>V</mi> <mi>o</mi> <mi>l</mi> <mi>t</mi> <mi>a</mi> <mi>g</mi> <mi>e</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>a</mi> <mo>&amp;CenterDot;</mo> <mi>t</mi> <mo>+</mo> <mi>b</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mn>1</mn> <mi>t</mi> </mfrac> <mo>+</mo> <mi>c</mi> <mo>&amp;CenterDot;</mo> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> </mrow>
其中,a,b,c,d为待估计参数,可以通过梯度下降法估算;
(8)、利用经验模型,计算燃料电池剩余剩余寿命Rulempirial
Rulempirial=arg(Vfailure(t))-tcurrent
其中,arg(Vfailure(t))表示经验模型到预设电压阈值的时间,tcurrent为测试数据的当前时间;
(9)、计算训练数据和测试数据的电压梯度值gradζ(t):
gradζ(t)=Vζ(t)-Vζ(t-1)
其中,ζ=1,2,V1(t)和grad1(t)分别表示t时刻训练数据的电压值和电压梯度值;V2(t)和grad2(t)分别表示t时刻测试数据的电压值和电压梯度值;
(10)、计算指标β:
<mrow> <mi>&amp;beta;</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>u</mi> <mi>r</mi> <mi>r</mi> <mi>e</mi> <mi>n</mi> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>d</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> </mrow> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>u</mi> <mi>r</mi> <mi>r</mi> <mi>e</mi> <mi>n</mi> <mi>t</mi> </mrow> </msub> </munderover> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>grad</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>-</mo> <msub> <mi>grad</mi> <mn>2</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>u</mi> <mi>r</mi> <mi>r</mi> <mi>e</mi> <mi>n</mi> <mi>t</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>d</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> </mrow> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>u</mi> <mi>r</mi> <mi>r</mi> <mi>e</mi> <mi>n</mi> <mi>t</mi> </mrow> </msub> </munderover> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>grad</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow> </mfrac> </mrow>
其中,tcurrent表示当前时间;Tdef为预设值,用来决定用几个单位的电压梯度值来计算相似度;
(11)、估算燃料电池的剩余寿命Rulfinal
Rulfinal=α·Rulhsmm+(1-α)·Rulempirial
<mrow> <mi>&amp;alpha;</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>&amp;beta;</mi> <mo>&amp;GreaterEqual;</mo> <mn>0.7</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>&amp;beta;</mi> <mo>&lt;</mo> <mn>0.7</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
2.根据权利要求1所述的基于HSMM和经验模型的燃料电池故障预测方法,其特征在于,所述的指标β用来衡量训练数据和测试数据之间的相似度,β值越小,则相似度越高;β值越大,相似度越低。
CN201710499198.6A 2017-06-27 2017-06-27 一种基于hsmm和经验模型的燃料电池故障预测方法 Active CN107169243B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710499198.6A CN107169243B (zh) 2017-06-27 2017-06-27 一种基于hsmm和经验模型的燃料电池故障预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710499198.6A CN107169243B (zh) 2017-06-27 2017-06-27 一种基于hsmm和经验模型的燃料电池故障预测方法

Publications (2)

Publication Number Publication Date
CN107169243A true CN107169243A (zh) 2017-09-15
CN107169243B CN107169243B (zh) 2019-08-20

Family

ID=59826187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710499198.6A Active CN107169243B (zh) 2017-06-27 2017-06-27 一种基于hsmm和经验模型的燃料电池故障预测方法

Country Status (1)

Country Link
CN (1) CN107169243B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615003A (zh) * 2018-12-06 2019-04-12 哈尔滨工业大学 一种基于elm-chmm的电源故障预测方法
CN110909455A (zh) * 2019-11-06 2020-03-24 电子科技大学 一种延缓固体氧化物燃料电池性能衰退的方法
CN111983475A (zh) * 2020-08-24 2020-11-24 哈尔滨理工大学 一种基于隐马尔科夫的锂离子动力电池安全度评估方法及装置
CN112687926A (zh) * 2019-10-18 2021-04-20 中国科学院大连化学物理研究所 一种燃料电池衰减趋势预测、剩余使用寿命预测方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070005311A1 (en) * 2005-04-25 2007-01-04 Smartsignal Corporation Automated model configuration and deployment system for equipment health monitoring
CN103698709A (zh) * 2013-12-11 2014-04-02 清华大学 一种燃料电池剩余寿命预测方法
CN105893208A (zh) * 2016-03-31 2016-08-24 城云科技(杭州)有限公司 基于隐半马尔可夫模型的云计算平台系统故障预测方法
CN106371030A (zh) * 2016-08-29 2017-02-01 丹阳亿豪电子科技有限公司 一种非确定性推理的新能源汽车电池故障诊断方法
CN106599920A (zh) * 2016-12-14 2017-04-26 中国航空工业集团公司上海航空测控技术研究所 一种基于耦合隐半马尔可夫模型的航空轴承故障诊断方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070005311A1 (en) * 2005-04-25 2007-01-04 Smartsignal Corporation Automated model configuration and deployment system for equipment health monitoring
CN103698709A (zh) * 2013-12-11 2014-04-02 清华大学 一种燃料电池剩余寿命预测方法
CN105893208A (zh) * 2016-03-31 2016-08-24 城云科技(杭州)有限公司 基于隐半马尔可夫模型的云计算平台系统故障预测方法
CN106371030A (zh) * 2016-08-29 2017-02-01 丹阳亿豪电子科技有限公司 一种非确定性推理的新能源汽车电池故障诊断方法
CN106599920A (zh) * 2016-12-14 2017-04-26 中国航空工业集团公司上海航空测控技术研究所 一种基于耦合隐半马尔可夫模型的航空轴承故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WU, XJ 等: "A hybrid prognostic model applied to SOFC prognostics", 《INTERNATIONAL JOURNAL OF HYDROGEN ENERGY》 *
彭颖: "基于退化隐式半马尔科夫模型的设备健康预测及系统性维护策略研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615003A (zh) * 2018-12-06 2019-04-12 哈尔滨工业大学 一种基于elm-chmm的电源故障预测方法
CN109615003B (zh) * 2018-12-06 2022-11-08 哈尔滨工业大学 一种基于elm-chmm的电源故障预测方法
CN112687926A (zh) * 2019-10-18 2021-04-20 中国科学院大连化学物理研究所 一种燃料电池衰减趋势预测、剩余使用寿命预测方法及系统
CN112687926B (zh) * 2019-10-18 2021-10-26 中国科学院大连化学物理研究所 一种燃料电池衰减趋势预测、剩余使用寿命预测方法及系统
CN110909455A (zh) * 2019-11-06 2020-03-24 电子科技大学 一种延缓固体氧化物燃料电池性能衰退的方法
CN110909455B (zh) * 2019-11-06 2023-04-18 电子科技大学 一种延缓固体氧化物燃料电池性能衰退的方法
CN111983475A (zh) * 2020-08-24 2020-11-24 哈尔滨理工大学 一种基于隐马尔科夫的锂离子动力电池安全度评估方法及装置
CN111983475B (zh) * 2020-08-24 2022-12-30 哈尔滨理工大学 一种基于隐马尔科夫的锂离子动力电池安全度评估方法及装置

Also Published As

Publication number Publication date
CN107169243B (zh) 2019-08-20

Similar Documents

Publication Publication Date Title
CN107169243A (zh) 一种基于hsmm和经验模型的燃料电池故障预测方法
CN112765772B (zh) 一种基于数据驱动的动力电池剩余寿命预测方法
CN108734360B (zh) 一种基于修正的elm预测模型多维遥测数据智能判读方法
Zhou et al. Remaining useful life prediction of individual units subject to hard failure
CN104699606A (zh) 一种基于隐马尔科夫模型的软件系统状态预测方法
KR102439041B1 (ko) 신경망 기반의 배터리 셀 불량 및 화재 사전 진단 방법 및 장치
CN113722985B (zh) 航空发动机健康状态评估和剩余寿命预测方法及系统
CN110412470B (zh) 电动汽车动力电池soc估计方法
CN110222431A (zh) 基于门控循环单元神经网络和卡尔曼滤波模型融合的锂离子电池剩余寿命预测方法
CN111680848A (zh) 基于预测模型融合的电池寿命预测方法及存储介质
CN112327188A (zh) 一种模型-数据混合驱动的锂离子电池剩余寿命预测方法
CN102542155A (zh) 基于加速退化数据的粒子滤波剩余寿命预测方法
CN107887023A (zh) 一种基于相似性和双随机游走的微生物‑疾病关系预测方法
CN109949929A (zh) 一种基于深度学习大规模病历的辅助诊断系统
CN102509023A (zh) 航天驱动组件综合应力加速寿命试验损伤累积模型的建模方法
CN104462757A (zh) 基于监测数据的Weibull分布可靠性序贯验证试验方法
CN109523021A (zh) 一种基于长短时记忆网络的动态网络结构预测方法
CN113095367B (zh) 一种基于深度强化学习的压缩机故障诊断方法
CN107133400A (zh) 一种飞机结构疲劳可靠度贝叶斯组合预测方法
CN107220500B (zh) 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法
Xia et al. Lithium-ion battery ageing behavior pattern characterization and state-of-health estimation using data-driven method
CN103678886B (zh) 一种基于地面测试数据的卫星贝叶斯网络健康确定方法
CN111369120A (zh) 一种基于设备迁移学习的非侵入式负荷监测方法
CN106896219A (zh) 基于变压器油中溶解气体数据的变压器亚健康状态识别及平均剩余寿命估计方法
CN104914850A (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