CN102208028B - 一种适用于动态复杂系统的故障预测和诊断方法 - Google Patents
一种适用于动态复杂系统的故障预测和诊断方法 Download PDFInfo
- Publication number
- CN102208028B CN102208028B CN 201110144382 CN201110144382A CN102208028B CN 102208028 B CN102208028 B CN 102208028B CN 201110144382 CN201110144382 CN 201110144382 CN 201110144382 A CN201110144382 A CN 201110144382A CN 102208028 B CN102208028 B CN 102208028B
- Authority
- CN
- China
- Prior art keywords
- data
- fault
- prediction
- dynamic complex
- time sequence
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000003745 diagnosis Methods 0.000 claims abstract description 43
- 238000013528 artificial neural network Methods 0.000 claims abstract description 39
- 230000004927 fusion Effects 0.000 claims abstract description 28
- 238000001514 detection method Methods 0.000 claims abstract description 18
- 238000007781 pre-processing Methods 0.000 claims abstract description 13
- 238000011058 failure modes and effects analysis Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims description 29
- 238000012544 monitoring process Methods 0.000 claims description 28
- 238000005070 sampling Methods 0.000 claims description 21
- 238000012549 training Methods 0.000 claims description 17
- 238000005315 distribution function Methods 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 10
- 238000012952 Resampling Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000013519 translation Methods 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 108091092919 Minisatellite Proteins 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000012986 modification Methods 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 241001123248 Arma Species 0.000 claims 1
- 238000003860 storage Methods 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000001976 improved effect Effects 0.000 description 4
- 239000008186 active pharmaceutical agent Substances 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000002955 isolation Methods 0.000 description 3
- 108091092878 Microsatellite Proteins 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000003446 memory effect Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000008093 supporting effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出一种适用于动态复杂系统的故障预测和诊断方法,可应用与航天器等动态复杂系统故障预测及诊断领域。本发明方法首先对动态复杂系统进行FMEA分析,得到主要故障模式和相应的性能检测参数,并性能检测参数划分为缓变数据和快变数据,再对性能检测参数进行预处理,针对缓变数据建立ARMA模型进行时序预测,针对快变数据建立多分辨率小波神经网络进行时序预测,然后建立预测区间模型对时序预测结果进行故障预警,最后通过建立基于D-S证据理论的多信号融合模型进行故障诊断。本发明方法能够对动态复杂系统进行高精度的故障预测及诊断,具有较强的通用性。
Description
技术领域
本发明属于航天器等动态复杂系统故障预测及诊断领域,为该领域内具有结构复杂且具有大量性能监测参数数据的系统提供了一套故障预测和诊断方法。
背景技术
故障预测是指根据系统过去和现在的状态,针对不确定事件,利用已有的知识,采用预测推理的方法预测系统在未来时间段出现故障的时间和发生故障的种类,找出故障原因和发生故障部件,为有计划地安排维修和维护提供理论依据。
故障预测技术最初是为机械的应用发展起来的,并在直升机和复杂的机械系统与HUMS子系统里所采用。对于故障预测技术,从其实际研究中应用的理论、方法和技术路线来看,可分为四类:基于模型的故障预测技术;基于知识的故障预测技术;基于数据的故障预测技术;基于传统可靠性的故障预测技术。
1.基于模型的故障预测技术对于大多数工业应用来说,基于物理模型的故障预测可能不是最切合实际的解决方案,因为对于不同的部件来说,某部件的当前故障类型通常是唯一的,并且很难在不终止运行的情况下鉴别出该故障。
2.基于知识的故障预测技术基于知识的方法不需要对象精确的数学模型,同时能够有效的表达对象相关的领域专家的经验知识,最大优势是能够充分利用对象系统有关的领域专家经验知识。但是,由于基于知识的故障预测技术本身更适合于定性推理而不太适合于定量计算,因而限制了其实际应用。
3.基于数据的故障预测技术基于数据的故障预测技术不需要对象系统的先验知识(数学模型和专家经验),以采集的数据为基础,通过各种数据分析处理方法挖掘其中的隐含信息进行预测操作,从而避免了基于模型和基于知识的故障预测技术的缺点,成为了一种较为实用的故障预测方法。但是,实际应用中一些关键设备的典型数据(历史工作数据、故障注入数据以及仿真实验数据)的获取代价通常十分高昂;而且即使对于所获得的数据来说,往往其具有很强的不确定性和不完备性。这些问题都增加了故障预测技术的实现难度。
4.基于传统可靠性的故障预测技术传统的可靠性评估方法都是基于同类部件/设备/系统的事件记录的分布。目前的可靠性方法基本上均是采用历史失效数据去估计对象的整体特性,如:平均故障间隔时间(MTBF)、平均故障前时间(MTTF)、可靠运行概率等。然而,这些方法仅提供了对同类对象的整体评估。此类可靠性评估通常对于进行大批量生产的制造商比较有价值,而对于最终的用户却没有多少价值可言。
自20世纪90年代以来,预测区间技术在国际上发展较为迅速。作为故障预测中的关键技术,预测区间技术是通过建立遥测数据预测误差的统计模型,进一步得到未来状态参数的预测区间,以提高故障预测的精度和工程实用性。但是,由于预测区间技术的研究还不甚成熟,现有方法大多用来解决预测残差分布已知的情形,针对在实际工程中常常遇到的残差分布未知的情况研究不多。
证据理论全称为Dempster-Shafer(D-S)证据理论,它首先由Dempster提出,并由Shafer进一步发展起来的一种不确定推理理论。证据理论是经典概率论的一种扩充形式,它基于对某一假设之上的各个片断的证据组合来估计假设的真实性。证据理论的基本策略就是将证据集合划分成多个不相关的部分,利用它们对识别框架独立进行判断,然后用证据合成规则,将多个证据的判断结果组合起来。这样,经过D-S证据结合,多个具体证据组合成一个抽象证据,该抽象证据综合了具体证据的信息,聚焦了具体证据的共同支持点,但形式上抽象证据与具体证据有相同的表达,这使得证据理论具有融合的开放性和处理的兼容性等优点。
动态复杂系统是复杂性科学的基本研究对象,广泛存在于工业、航空航天、生物、化学、物理、经济、管理等众多领域。整体涌现性和动态性是复杂系统的基本特征,即动态复杂系统具有不同于各组分特征线性加和的特征,如航天器、工业流水线、数控机床等。
动态复杂系统常常在国民经济中扮演者重要的角色,它们的健康稳定运行至关重要。但是由于动态复杂系统具有组成结构复杂、工作状态多变和易受运行环境影响的特点,前述四类主要的故障预测和诊断方法难以适用于此类系统。因此,亟需一种针对这类重要系统的故障预测和诊断方法。
发明内容
本发明为了解决由于系统组成结构、运行环境复杂等因素所造成的难以建立精度较高故障预测模型,以及现有方法无法直接、准确解决动态复杂系统的故障预测和诊断问题,提供了一种适用于动态复杂系统的故障预测和诊断方法。它能够实现三大功能:实现对动态复杂系统的故障预警;能够对可能发生故障的单元进行诊断隔离;能够给出故障预测的置信区间。
本发明一种适用于动态复杂系统的故障预测和诊断方法,具体包括以下步骤:
步骤一、对动态复杂系统的结构、性能和历史故障数据进行FMEA(Failure Mode andEffect Analysis,失效模式和效果分析)分析,得到动态复杂系统的主要故障模式和各故障模式相应的性能检测参数,并将性能检测参数划分为缓变数据和快变数据;
步骤二、对动态复杂系统的性能检测参数进行预处理,所述的预处理包括:剔除奇异值和滤波降噪;
步骤三、建立单参数时序预测模型:针对缓变数据,通过ARMA模型建立缓变数据的时序预测模型;针对快变数据,通过多分辨率小波神经网络建立快变数据的时序预测模型;将预处理的性能检测参数输入相应的时序预测模型后输出时序预测结果;
步骤四、建立预测区间模型,通过预测区间对时序预测结果进行故障预警,当时序预测结果超过预测区间的范围时进行故障预警;
步骤五、通过建立基于D-S证据理论的多信号融合模型进行故障诊断。
所述步骤四中的预测区间模型的建立包括步骤:构建k个基于向量化观测值的抽样块;对k个抽样块进行重抽样;确定各统计量;确定概率α的分位点值;得到预测区间。其中,对k个抽样块进行重抽样的次数可以进行自动调整,具体是:首先定义数据变动率Di:Si表示第i个数据块,max(Si)是抽样块Si中的最大值,min(Si)是抽样块Si中的最小值;然后按照数据变动率的定义分别求出各抽样块的数据变动率,若相邻抽样块的数据变动率之比Di/Di+1大于1,则对抽样块Si+1的抽样次数进行调整,得到抽样块Si+1新的重抽样次数M′:M′=M/(Di/Di+1),M为对数据块进行重抽样的次数的初始值。
所述步骤五中进行故障诊断具体是:构建动态复杂系统的多信号模型;将动态复杂系统的所有可能故障单元组合组成了识别框架;构造基本概率分配函数,并结合每个性能监测参数报警与对应故障模式的可能性进行修正;按照证据融合规则,确定证据融合结果及各种故障模式的置信区间;对动态复杂系统是否存在的故障模式进行判断。
本发明的优点与积极效果在于:
(1)本发明方法克服了众多预测和诊断方法限制条件苛刻的弊端,针对动态复杂系统可以进行高精度的故障预测及诊断,且针对具有动态复杂特点的一类对象的故障预测及诊断,具有较强的通用性;
(2)本发明方法在故障预测及诊断过程中,无需建立故障影响因素、预测对象性能表征参数和故障之间的解析函数关系等,相比现有方法操作简单;
(3)本发明方法中采用的多分辨率小波神经网络结合了小波多分辨理论和神经网络理论,有机融合了两者的优点,即小波变换的时频域局部特性和变焦特性,神经网络的自学习、自适应、鲁棒性、容错性和推广能力,通过建立多分辨率小波神经网络实现了对快变数据的时序预测;
(4)本发明方法中预测区间采用改进的块状自助法来建立,改进了传统区间预测的方式并弥补了块状自助法的不足和缺陷,解决了预测误差分布未知时的区间预测问题,提高了故障预测的实用性和精度;
(5)本发明提出的利用D-S证据理论来对故障隔离模糊组进行分解的方法,很好了解决了如何将专家经验融入故障诊断的问题,进而实现故障诊断技术从理论研究到工程实际应用的一次跨越。
附图说明
图1为本发明的故障预测和诊断方法的步骤流程图;
图2为本发明方法中步骤三中建立ARMA模型的流程示意图;
图3为本发明实施例中蓄电池放电初压经小波降噪的效果对比图;
图4为本发明实施例中ARMA模型预测出的温度预测曲线和预警示意图;
图5为本发明实施例中多分辨率小波神经网络预测出的蓄电池放电初压的曲线图;
图6为本发明实施例中多分辨率小波神经网络预测出的母线电压的曲线图;
图7为本发明实施例母线电压预测区间的曲线示意图图;
图8为本发明实施例HY-1B小卫星电源系统的多信号模型。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明方法适用于动态复杂系统的故障预测和诊断,其主要原理如下:首先通过性能监测参数的时序预测和区间分析来实现故障预警,然后对相应的故障预警单元利用基于D-S证据理论的多信号融合诊断模型进行诊断推理,预测出未来可能发生的故障模式及其概率。本发明方法包含两大核心内容:时序预测方法和诊断推理方法。下面结合图1说明本发明适用于动态复杂系统的故障预测和诊断方法。
步骤一:对动态复杂系统的结构、性能和历史故障数据进行FMEA分析,得到动态复杂系统的主要故障模式以及与各故障模式相关联的性能监测参数。
通过传感器搜集与动态复杂系统的结构、性能相关的数据以及历史故障数据并保存在数据库中,根据数据库中的数据确定动态复杂系统的关重件以及不同任务剖面下各部件的工作状态,然后对动态复杂系统进行FMEA分析得到该动态复杂系统的主要故障模式和各故障模式相关联的性能监测参数如遥测数据。由于性能监测参数的特点决定了数据预处理和预测方法的选取,因此在步骤一中把性能监测参数按照参数采样的频率分为缓变数据和快变数据。根据惯例采样间隔时间高于32秒的为缓变数据,反之为快变数据。
所述的FMEA是分析系统中每一产品可能产生的故障模式及其对系统造成的所有可能影响,并按每一个故障模式的严重程度及其发送概率予以分类的一种归纳分析方法,FMEA可采用国军标GJB1391标准。
步骤二:对动态复杂系统的性能监测参数进行预处理。
性能监测参数是了解动态复杂系统运行状态的基础,它包含了动态复杂系统的运行和故障信息。但是,由于动态复杂系统具有结构复杂、运行状态受运行环境影响较大的特点,如振动、电磁辐射以及运行状态切换指令的影响,因此,得到的性能监测参数通常含有噪声。由于这些噪声严重干扰了有用信号,所以,为了得到更精确的故障信息、提高故障预测的精度,就首先需要在对动态复杂系统进行故障预测和诊断前,对各故障模式相应的性能监测参数进行滤波降噪。有用信号通常表现为低频信号或是一些比较平稳的信号,而噪声信号则通常表现为高频信号。
本发明的预处理主要包含两部分内容:剔除奇异值和滤波降噪。奇异点的剔除,主要通过将严重超出正常范围的个别点剔除实现。滤波降噪主要做如下处理:首先对原始性能检测参数进行小波分解,噪声通常包含在高频系数中;然后对小波分解的高频系数以门限阈值等形式进行量化处理;最后再对性能检测参数重构。
步骤三:建立单参数时序预测模型,针对缓变数据,通过ARMA模型建立缓变数据的时序预测模型;针对快变数据,通过多分辨率小波神经网络建立快变数据的时序预测模型;将预处理的性能检测参数输入相应的时序预测模型后输出时序预测结果。
单参数时序预测是故障预警的基础,因此需要预测方法具有较高的预测精度。预测方法的选取主要根据性能监测参数的特点进行。在步骤一中,将性能监测参数划分为两类:缓变数据和快变数据。
缓变数据如温度属于平稳数据,选用适合平稳数据预测的ARMA模型(Auto-Regressiveand Moving Average Model,自回归滑动平均模型)进行预测。实际的时间序列如航天数据、股票等,在演化过程中并不具有固定不变的均值,除去局部水平和趋势不同以外,时间序列显示出具有某种意义上的同质性,即时间序列的某一部分与任何其他部分极其相似,对该时间序列作适当的差分便可使之平稳化,就得到描述同质非平稳特性的模型,也就是ARMA模型。具体如图2所示,对非平稳时间序列(遥测数据)xt作差分得到变化量如果一阶差分还是非平稳序列,再进一步对差分结果作差分,直到得到平稳序列为止,得到平稳序列后,通过参数估计、定阶、拟合检验最终建立ARMA模型,图2中d=1表示一阶差分,表示d阶差分后的时间序列xt。对缓变数据经如图2所示的ARMA模型处理后得到时序预测结果。
快变数据如电压属于非平稳数据,这类数据的预测要求预测方法具有较快的速度和较强的跟踪能力,因此本发明以神经网络预测模型为基础,利用小波多分辨理论对其进行改进,通过建立多分辨率小波神经网络来实现快变数据的时序预测模型。快变数据的时序预测模型采用的是多分辨率小波神经网络来建立的,建立的步骤如下:
步骤3.1.建立多分辨率小波神经网络,主要包括网络的层数、每层的节点数等等。由于多分辨率小波神经网络自身的特点,多分辨率小波神经网络不仅要设置网络权值、闭值的初始值,还要设置尺度参数m和平移参数n的初始值。m初始值的选取可以提高网络的收敛性能,一般m选取大于1的值。
一般网络的层数和每层节点的个数主要是靠经验和反复实验所得,而且预先并不固定小波函数尺度参数m和平移参数n的大小,而是让它们无约束地自由取值。本发明给出了确定小波函数尺度参数m和平移参数n范围的方法:
尺度参数m∈[mmin,mmax]:
mmin=int-∞((lnωmin-lnω1)/ln2)
mmax=int+∞((lnωmax-lnω0)/ln2)
其中,int-∞(·)和int+∞(·)分别表示向-∞和+∞方向取整,ωmin和ωmax为通过对训练数据作频率带宽估计获得的最小频率带宽和最大频率带宽,ω0和ω1分别训练数据频率区间的下界频率带宽和上界频率带宽。
确定了小波的尺度参数后,就可以根据待处理信号的时域空间区域来确定每个伸缩参数m∈[mmin,mmax]下的平移参数n了。对于每个尺度参数m,只有选择小波函数ψmn位于[tmin,tmax]的平移尺度才会对逼近产生作用。由此,可得到对于每个m下覆盖待处理信号整个时域空间的平移参数n的取值范围为[nmin,nmax]:
nmin=int-∞(2mtmin-t1)
nmax=int+∞(2mtmax-t0)
其中,tmin表示时域空间的起点时刻,tmax表示时域空间的终点时刻,t0和t1分别表示时域空间的下界时刻和上界时刻。由m、n值决定的小波函数Ψmn加权和描述了待处理信号的时频局部特性,同时,也就得到了多分辨率小波神经网络的结构。
步骤3.2.进行输入向量的预处理。多分辨率小波神经网络的输入向量初始值在取[0,1]或者[-1,1]之间时,具有较好的收敛性能,因此,才为了提高多分辨率小波神经网络的收敛速度,对多分辨率小波神经网络的输入数据进行预处理,将输入向量的每个分量化归到[0,1]或者[-1,1]之间。下面以将输入向量的每个分量化归到[0,1]之间为例说明。本发明中对输入向量进行预处理的方法是采用的归一化处理方法,首先找到输入向量的最大最小值,然后对输入向量值归一化,最后得到的结果是:
其中,ximin为输入向量的最小值,ximax为输入向量的最大值,xik为第k个输入向量的值,为第k个输入向量归一化后得到的值。通过归一化处理,将所有输入向量的分量将被归一化到[0,1]之间。
对于多分辨率小波神经网络的输出向量,通过下式进行反归一化变换:
步骤3.3.将预处理后的输入向量输入多分辨率小波神经网络,对多分辨率小波神经网络进行训练。为了提高网络训练的速度,本发明引入加速因子。
加速因子使网络在修正其权值时,不仅考虑误差在梯度上的作用,而且考虑在误差曲面上变化趋势的影响。在没有加速因子的作用下,网络可能陷入浅的局部极小值,利用加速因子的作用有可能滑过这些极小值。
带有加速因子的权值调节公式为:
Δwij(k+1)=(1-a)ηδi+aΔwij(k)
其中Δwij(k)为权值变化量,k为训练次数,α为加速因子,一般取0.95左右,η为学习步长,δi为第i个隐层节点的误差。
步骤3.4.利用训练好的多分辨率小波神经网络对快变数据进行时序预测。
步骤3.5.根据得到的时序预测结果,判断预测精度是否达到要求,若预测精度未达到要求,则对网络结构和参数如网络层数、权值初值等作修改,转步骤3.3重新训练;若预测精度达到了要求,则建立完成了快变数据的时序预测模型。
通过上述两种时序预测模型对动态复杂系统的故障模式相应的性能监测参数进行预测分析,以预测该动态复杂系统的某子系统在未来是否会出现故障。如果预测可能出现故障,则发出故障预警,以便及时对故障进行诊断,避免故障扩大和传播。
步骤四、建立预测区间模型,对经步骤三得到的时序预测结果通过预测区间的上、下限比较进行故障预警,当时序预测结果超过预测区间的范围时进行故障预警。
利用步骤三的预测结果虽然可以进行故障预测,但是在实际生产生活中,动态复杂系统的运行过程往往会因受到外界复杂环境的影响而表现为性能监测参数的波动,而单点预测往往无法识别这种波动,因此无法满足工程实际需求。所以故障预测的精度不可能达到100%,真实值与预测值总会存在误差,所以依靠单点预测难以保证动态复杂系统故障预警的准确性。本发明通过建立的预测区间模型,利用具有一定置信度的预测区间取代单点预测值,然后在此基础上通过预测区间的上、下限进行故障预警,这样就降低了不确定性因素对故障预警精度的影响,进一步增强了预警的工程实用性。
将经步骤三所建立的单参数时序预测模型进行预测后得到的时序预测结果作为一个连续观测值,针对一个连续观测值的向量Y=(X1,…,Xn),n为向量的维数,具体如下构建预测区间模型:
步骤4.1.构建基于向量化观测值的抽样块(Block Bootstrap)。
首先构建k个连续交叉的抽样块S1=(Xm,…,Xm+l-1),S2=(Xm+1,…,Xm+l),...,Sk=(Xn-l+1,…,Xn);其中m是指抽样起点值Xm在向量Y中的序号,l是抽样块的长度,一般取l=n1/3,n为观测值向量的维数,k取正整数。
步骤4.2.对k个抽样块进行重抽样。
对数据块进行重抽样的次数M一般设定为50次以上。但抽样次数确定后,缺少根据数据特点变动的能力,抽样难以跟随数据的特点,会降低抽样的效率,为此本发明对抽样次数的确定方法进行了改进,实现了重抽样次数根据数据波动大小进行自动调整,既提高了所获取数据变动趋势的精度,又提高了抽样效率。
然后按照数据变动率的定义分别求出各抽样块的数据变动率,若相邻抽样块的数据变动率之比Di/Di+1大于1,则表明数据的趋势发生了较大的变化,因此需要对抽样块Si+1的抽样次数进行调整,得到新的重抽样次数M′:M′=M/(Di/Di+1)。
对k个抽样块S1=(Xm,…,Xm+l-1),S2=(Xm+1,…,Xm+l),...,Sk=(Xn-l+1,…,Xn)进行非独立的替换抽样,抽样得到的样本称为Bootstrap复制样本对得到的每个Bootstrap复制样本内的数据进行N次重抽样,N取大于50的整数,得到N个Bootstrap样本对每个样本进行抽样得到标准差的估计值
步骤五:通过建立基于D-S证据理论的多信号融合模型进行故障诊断。
在进行故障预测时,虽然通过对性能监测参数的预测可以预测故障演化的过程,但是仅仅依赖单一性能监测参数的预测信息进行故障诊断并不可信。要获得可信度较高的预测诊断结果,需要融合所有性能监测参数的预测信息进行综合全面的预测诊断。因此,本发明提出一种改进的带证据权值的D-S证据理论的多信号融合模型来解决上述的诊断问题。
利用多信号融合模型进行故障诊断隔离时常存在系统组件与性能监测参数之间的依赖矩阵并不完全确定的问题,即参数与故障并不一定存在一一对应关系,这就很难将故障隔离到独立的组件,而只能隔离到某个模糊组(即包含多个可能故障组件的集合)。因此,本发明利用证据权值的DS证据理论对多信号融合模型进行改进,借助专家经验来解决了上述故障诊断中的融合诊断问题,最终给出故障隔离结果。
基于DS证据理论的多信号融合模型进行故障诊断具体步骤如下:根据信任区间按拟定的决策规则得到最终决策。下面针对每一步骤进行详述。
步骤5.1.构建动态复杂系统的多信号模型。
利用测试性工程和维护系统(TEAMS)软件,结合动态复杂系统的结构原理图,构建其多信号模型。
步骤5.2.形成识别框架。
在证据理论中,一个样本空间称为一个识别框架,用Ω表示。识别框架Ω由一系列对象θi组成:Θ={θ1,θ2,...,θi},θi称为Θ的一个单子,只含一个单子的集合称为单子集合。证据理论的基本问题是:已知识别框架Θ,判明Θ中一个先验的未定元素θi属于Θ中某一个子集A的程度。在动态复杂系统故障预测中,需要判断哪些单元发生了故障,因此动态复杂系统的所有可能故障单元组合组成了识别框架。
步骤5.3构造基本概率分配函数。令识别框架Θ为一论域集合,2Θ为该Θ的所有子集构成的集合,则基本概率分配函数为:m:2Θ→[0,1],且基本概率分配函数满足如下公理:
其中,P(Θ)为幂集,m(A)表示子集A的概率分配。
专家根据经验知识、动态复杂系统的结构与地面测试的各性能监测参数的预测时序结果,给出每个性能监测参数报警与对应故障模式的可能性,相当于给出了D-S证据推理中的mass值。根据每个性能监测参数报警的可能性不同,它对证据支持程度不一样。把性能监测参数报警的可能性μ(ti)当作证据权值,ti表示第i个故障模式,对原始D-S证据基本概率分配函数进行修正,得到的修正后的基本概率分配函数为:
ma(A)表示修正后的该子集A的概率分配。
可以看出,修正的结果是使得传感器对具体的故障诊断的基本概率分配值减小,也就是降低了它的确定性故障诊断,增加了不确定故障诊断。这样就可以减少权值小的传感器提供的确定性信息,减小错误结果对整个系统的故障诊断的影响。
步骤5.4.按照证据融合规则,确定证据融合结果及各种故障模式的置信区间。证据理论的一个基本策略是将证据集合划为2个或多个不相关的部分,并利用它们分别对识别框架Θ独立进行判断,然后用Dempster组合规则将它们融合起来,其融合规则为:
根据融合结果,确定各种故障模式的置信区间,可采用计算证据推理理论中的信任区间方法得到。所述的信任区间表示为(Bel(A),Pl(A)),其中,Bel(A)为信息测度,Pl(A)为似然测度,信任测度与似然测度分别表示对判断信任程度估计的下限估计(悲观估计)与上限估计(乐观估计)。
步骤5.5.故障诊断。经D-S证据理论的融合规则融合后,可得到各个故障模式的置信区间,可按如下规则对动态复杂系统是否存在该故障模式进行判断:
规则一、故障元件具有最大的可信度,且Bel(A)>ε1,Pl(A)>ε2;
规则二、不确定的可信度m(θ)≤ε3。
式中ε1、ε2和ε3是用户定义的阈值,阈值的选取主要依照对象特点和经验。
实施例
本实施例以我国小卫星电源系统为预测对象,小卫星电源系统是典型的动态复杂系统。由于小卫星电源系统结构复杂、缺少统一的物理模型、故障机理复杂,符合本发明所需要解决的动态复杂系统的故障预测问题。通过本实施例的详细阐述,进一步说明本发明的实施过程及工程应用过程。
本发明实施例对小卫星电源系统应用本发明提出的故障预测和诊断方法的步骤如下:
步骤一:对小卫星电源系统进行FMEA分析,得到小卫星电源系统的主要故障模式和与各故障模式相关联的性能监测参数。
对于小卫星电源系统,它包括蓄电池、太阳电池阵和电源控制器三大部分,故障预测和诊断主要针对这三大部分。通过资料以及和卫星专家沟通,目前小卫星电源系统的太阳电池阵、电源控制器可靠性较好,基本不会出现故障,蓄电池是主要的故障子系统。
小卫星蓄电池的主要故障模式有四种:插接件短路,单体电池短路,单体电池性能衰减,单体电池记忆效应。与这四种故障模式相关的性能监测参数都为:蓄电池放电初压、蓄电池温度、母线电压。蓄电池处于正常工作状态时,性能监测参数的正常阈值为:蓄电池温度阈值5-8℃,蓄电池放电初始电压阈值24-25v,母线电压阈值22-25v。
步骤二:故障预测相关状态参数的预处理。预处理主要为剔除奇异值和滤波降噪。滤波降噪采用小波降噪,小波分解层次为4,采用软阈值方法来确定降噪阈值,如图3所示为放电初压在经过小波降噪后的示意图。
步骤三:建立单参数时序预测模型。
(1)缓变数据的时序预测模型,本发明实施例中蓄电池温度属于缓变数据。
针对蓄电池温度的遥测数据,建立如下7阶ARMA模型:
x(t)=1.4662×x(t-1)-0.2452×x(t-2)-0.2297×x(t-3)-0.1531×x(t-4)
-0.2868×x(t-5)+0.6451×x(t-6)-0.1558×x(t-7)+ε(t)
ε(t)~WN(0,2.9208),
x(t)表示t时刻的蓄电池温度的遥测数据值,ε(t)表示误差,WN表示白噪声。
该7阶ARMA模型通过了白噪声检验。
根据如图4所示的预测结果,卫星继续运行170圈左右时,蓄电池温度超出阈值范围,将会进行故障预警。HY-1B小卫星每圈周期为100.7分钟,即大约11天后将可能发生故障,需要提前采取措施预防。这个预警可能是由于单体电池短路造成,也可能是用电载荷功率大引起。
(2)快变数据的时序预测模型。本发明实施例中,蓄电池放电初压与母线电压这两种属于快变数据。
步骤3.1.建立多分辨率小波神经网络模型:输入层节点为8个,输出层节点为1个,隐层节点为6个,网络核函数采用morlet小波。
步骤3.2.进行输入向量的预处理:归一化蓄电池放电初压;归一化母线电压。
步骤3.3.将预处理后的输入向量输入多分辨率小波神经网络,对多分辨率小波神经网络进行训练。利用808个数据进行训练,192个数据作为验证,采用梯度下降法进行训练。
步骤3.4.利用训练好的多分辨率小波神经网络进行仿真,得到预测的蓄电池放电初压和预测的母线电压。
步骤3.5.多分辨率小波神经网络训练参数及预测结果如下:
训练时间:1.0238s
预测精度:均方误差MSE=0.0213
如图5和图6所示,为多分辨率小波神经网络预测得出的蓄电池放电初压和母线电压曲线的对比示意图,通过图中与BP神经网络的预测曲线对比可以看出多分辨率小波神经网络预测精度较高。
步骤四:建立预测区间模型。
每个数据段的长度为l=n1/3=1251/3=5,数据段个数k=25,设置抽样次数初始值M=200。母线电压得到的预测区间曲线如图7所示。
步骤五:基于DS证据理论的多信号融合模型的故障诊断。
如图8所示为建立的HY-1B小卫星电源系统的多信号模型,该多信号模型中的故障模式,包括:插接件短路F1,单体电池短路F2,单体电池性能衰减F5,单体电池记忆效应F4,电池片参数变化F3。t1、t2和t3表示三个性能监测参数,t1表示蓄电池放电初压,t2表示母线电压,t3表示蓄电池温度。由于缺乏有效的信息,单体电池性能衰减F5与参数t1、t2+和t3之间的依赖矩阵不能确定,图8中用虚线表示故障模式与参数间的依赖关系不确定,但可由专家根据经验和已有的信息,给出当参数报警时电源系统故障模式发生的概率,即证据理论中的基本概率赋值(BPA)。
通过和东方红卫星公司的电源专家进行沟通,利用专家经验知识、系统结构与遥测参数预测结果,给出每个预警参数与相应故障模式发生的概率,这相当于D-S证据推理中的mass值。
当参数t1预警时,各故障模式发生的可能性为:
m1(F2)=0.4,m1(F5)=0.3,m1(F3)=0.05
m1(F1)=0.1,m1(F3,F4)=0.05,m1(θ)=0.1,
其中m1(F3,F4)=005表示当参数t1预警时,卫星电源系统故障模式F3,F4同时发生故障的概率为0.05,其他参数含义相似,不再赘述。
当参数t2预警时,各组件发生故障的可能性为:
m2(F4)=05,m2(F5)=03,m2(F3,F5)=01,m2(θ)=01。
当参数t3预警时,各元件发生故障的可能性为:
m3(F5)=05,m3(F2)=04,m3(θ)=01。
有前面的步骤对参数进行预测预警,得到参数预警可能性为μ(t1)=1.0,μ(t2)=1.0,下面介绍一下如何采用D-S证据推理判断将要发生的故障模式和相应的概率。
根据融合规则合成t1和t2,因为预测参数预警可靠度μ(t1)=1.0,μ(t2)=1.0,所以t1和t2的证据可信度不受影响,两参数的融合过程如下所示:
进行故障诊断,设置故障模式F5可信度最大,且满足
Bel(F5)=0.669>ε1=0.1,pl(F4)=0.756>ε2=0.5,m(θ)=0.024<ε3=0.5。
而其他故障模式则不满足故障条件,因此只有模式F5预测会发生,其他模式将不会发生,预测F5故障模式发生的时间为未来第11天,概率为0.669。
Claims (9)
1.一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,具体包括以下步骤:
步骤一、对动态复杂系统的结构、性能和历史故障数据进行FMEA分析,得到动态复杂系统的主要故障模式和各故障模式相应的性能检测参数,并将性能检测参数划分为缓变数据和快变数据;
步骤二、对动态复杂系统的性能检测参数进行预处理,所述的预处理包括:剔除奇异值和滤波降噪;
步骤三、建立单参数时序预测模型:针对缓变数据,采用ARMA模型建立缓变数据的时序预测模型;针对快变数据,采用多分辨率小波神经网络建立快变数据的时序预测模型;将预处理的性能检测参数输入相应的时序预测模型后输出时序预测结果;
步骤四、建立预测区间模型,通过预测区间对时序预测结果进行检测,当时序预测结果超过预测区间的范围时进行故障预警;
步骤五、通过建立基于D-S证据理论的多信号融合模型进行故障诊断。
2.根据权利要求1所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,步骤一中所述将性能检测参数划分为缓变数据和快变数据的方法是:当性能检测参数的采样间隔时间高于32秒时,将该性能监测参数划分为缓变数据,否则将该性能监测参数划分为快变数据。
3.根据权利要求1所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,步骤三中所述的采用多分辨率小波神经网络建立快变数据的时序预测模型,具体步骤如下:
步骤3.1、建立多分辨率小波神经网络,设置网络的层数、每层的节点数和网络参数的初始值,所述的网络参数包括网络权值、网络闭值、尺度参数和平移参数;
所述的尺度参数p∈[mmin,mmax]:
mmin=int-∞((lnωmin-lnω1)/ln2)(1)
mmax=int+∞((lnωmax-lnω0)/ln2)(2)
其中,int-∞(·)和int+∞(·)分别表示向-∞和+∞方向取整,ωmin和ωmax分别表示通过对训练数据作频率带宽估计获得的最小频率带宽和最大频率带宽,ω0和ω1分别表示训练数据频率区间的下界频率带宽和上界频率带宽;
对于每个尺度参数p下的覆盖待处理信号的整个时域空间的平移参数n的取值范围为[nmin,nmax]:
nmin=int-∞(2ptmin-t1)(3)
nmax=int+∞(2ptmax-t0)(4)
其中,tmin表示所述时域空间的起点时刻,tmax表示所述时域空间的终点时刻,t0和t1分别表示所述时域空间的下界时刻和上界时刻;
步骤3.2、对输入向量进行预处理,将输入向量归一化;
步骤3.3、将预处理后的输入向量输入多分辨率小波神经网络,对多分辨率小波神经网络进行训练,训练过程使用加速因子进行权值调节,带有加速因子的权值调节公式为:
Δwij(k+1)=(1-a)ηδi+aΔwij(k)(5)
其中,Δwij(k+1)与Δwij(k)分别为第k+1次训练与第k次训练的第i个隐层节点与第j个输出层节点之间的权值变化量,k表示训练次数,a表示加速因子,η表示学习步长,δi为第i个隐层节点的误差;
步骤3.4、利用训练好的多分辨率小波神经网络对快变数据进行时序预测,得到时序预测结果;
步骤3.5、根据得到的时序预测结果,判断预测精度是否达到要求,若达到要求,则建立完成快变数据的时序预测模型,若没有,则对多分辨率小波神经网络的层数和网络参数作修改,转步骤3.3重新训练。
4.根据权利要求3所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,步骤3.1中所述的尺度参数选取大于1的值。
5.根据权利要求3所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,步骤3.3中所述的加速因子为0.95。
6.根据权利要求1所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,步骤四中的预测区间模型具体通过下面步骤得到:
步骤4.1.构建基于向量化观测值的抽样块;
首先构建k个连续交叉的抽样块:S1=(Xq,…,Xq+l-1),S2=(Xq+1,…,Xq+l),...,Sk=(Xn-l+1,…,Xn);其中,n为连续观测值的向量Y=(X1,…,Xn)的维数,q是指抽样起点值Xq在向量Y中的序号,l是抽样块的长度;
步骤4.2.对k个抽样块进行重抽样;
然后按照数据变动率的定义分别求出各抽样块的数据变动率,若相邻抽样块的数据变动率之比Di/Di+1大于1,则对抽样块Si+1的抽样次数进行调整,得到抽样块Si+1新的重抽样次数M′:M′=M/(Di/Di+1),M为对数据块进行重抽样的次数的初始值;
对k个抽样块S1=(Xq,…,Xq+l-1),S2=(Xq+1,…,Xq+l),...,Sk=(Xn-l+1,…,Xn)进行非独立的替换抽样,抽样得到Bootstrap复制样本对得到的每个Bootstrap复制样本内的数据进行N次重抽样,得到N个Bootstrap样本对每个Bootstrap样本进行抽样得到标准差的估计值
7.根据权利要求6所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,步骤4.2中所述的对数据块进行重抽样的次数的初始值为大于50的整数,所述的对得到的每个Bootstrap复制样本内的数据进行重抽样的次数N取大于50的整数。
8.根据权利要求1所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,所述的步骤五具体包括以下步骤:
步骤5.1.构建动态复杂系统的多信号模型;
步骤5.2.形成识别框架:将动态复杂系统的所有故障单元组合组成识别框架Θ;
步骤5.3.构造基本概率分配函数:
其中,2Θ为识别框架Θ的所有子集构成的集合,A表示识别框架Θ中一个子集,m(A)表示该子集A的概率分配,ma(A)表示修正后的该子集A的概率分配,ti表示第i个故障模式;
步骤5.4.按照证据融合规则,确定证据融合结果及各种故障模式的置信区间;
所述的融合规则为:
根据融合规则进行证据融合得到证据融合结果,根据融合结果,确定各种故障模式的置信区间(Bel(A),Pl(A)),其中,Bel(A)表示信息测度,Pl(A)表示似然测度;
步骤5.5.按照下面两条规则对动态复杂系统是否存在该故障模式进行判断:
规则一、故障元件具有最大的可信度,且Bel(A)>ε1,Pl(A)>ε2;
规则二、不确定的可信度m(θ)≤ε3;
ε1、ε2和ε3是用户根据实际情况定义的阈值,θ表示虚警。
9.根据权利要求8所述的一种适用于动态复杂系统的故障预测和诊断方法,其特征在于,所述的步骤5.5中针对HY-1B小卫星电源系统,设置ε1为0.1,ε2为0.5,ε3为0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110144382 CN102208028B (zh) | 2011-05-31 | 2011-05-31 | 一种适用于动态复杂系统的故障预测和诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110144382 CN102208028B (zh) | 2011-05-31 | 2011-05-31 | 一种适用于动态复杂系统的故障预测和诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102208028A CN102208028A (zh) | 2011-10-05 |
CN102208028B true CN102208028B (zh) | 2013-06-19 |
Family
ID=44696850
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110144382 Expired - Fee Related CN102208028B (zh) | 2011-05-31 | 2011-05-31 | 一种适用于动态复杂系统的故障预测和诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102208028B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI551981B (zh) * | 2014-01-09 | 2016-10-01 | 富士通股份有限公司 | 分析方法、分析程式及分析裝置 |
Families Citing this family (60)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064340B (zh) * | 2011-10-21 | 2014-12-03 | 沈阳高精数控技术有限公司 | 一种面向数控机床的故障预测方法 |
CN102495949A (zh) * | 2011-11-23 | 2012-06-13 | 西北工业大学 | 一种基于飞行数据的故障预报方法 |
CN102426334B (zh) * | 2011-11-28 | 2013-07-24 | 北京航空航天大学 | 一种放大电路的贮存性能表征参数确定方法 |
CN102542159B (zh) * | 2011-12-08 | 2014-10-08 | 北京空间飞行器总体设计部 | 一种在轨航天器状态预测方法 |
CN102662325B (zh) * | 2012-05-04 | 2014-10-15 | 重庆邮电大学 | 一种改进自适应学习树电源管理方法 |
CN102707708B (zh) * | 2012-05-25 | 2014-04-09 | 清华大学 | 多模型飞行控制系统故障诊断方法及装置 |
CN103678918A (zh) * | 2013-12-16 | 2014-03-26 | 中国人民解放军63791部队 | 航天器飞行时序最小差值处理方法 |
CN104731664A (zh) * | 2013-12-23 | 2015-06-24 | 伊姆西公司 | 用于故障处理的方法和装置 |
CN104361026B (zh) * | 2014-10-22 | 2017-09-19 | 北京航空航天大学 | 一种fmea分析过程中的故障知识存储和推送方法 |
CN104462826B (zh) * | 2014-12-11 | 2017-07-18 | 云南师范大学 | 基于矩阵奇异值分解的多传感器证据冲突检测与度量方法 |
CN104598745B (zh) * | 2015-01-27 | 2017-07-07 | 北京空间飞行器总体设计部 | 一种月地高速再入返回飞行品质评估方法 |
CN105424224B (zh) * | 2015-12-14 | 2018-04-06 | 国家电网公司 | 一种隔离开关状态监测方法及装置 |
CN107024915B (zh) * | 2016-02-02 | 2019-10-01 | 同济大学 | 一种电网控制器板卡故障检测系统及检测方法 |
CN105825271B (zh) * | 2016-03-21 | 2018-10-19 | 南京邮电大学 | 基于证据推理的卫星故障诊断与预测方法 |
CN107222322A (zh) * | 2016-03-22 | 2017-09-29 | 中国移动通信集团陕西有限公司 | 一种通信故障诊断方法及装置 |
CN107065897B (zh) * | 2016-06-20 | 2019-12-24 | 浙江工业大学 | 三自由度直升机显式模型预测控制方法 |
CN106446506A (zh) * | 2016-08-26 | 2017-02-22 | 航天东方红卫星有限公司 | 一种多方法组合的小卫星状态自主智能预测方法 |
CN106779135A (zh) * | 2016-11-11 | 2017-05-31 | 上海海事大学 | 一种混合动力船舶负载功率预测方法 |
CN106742068B (zh) * | 2016-12-07 | 2019-01-04 | 中国人民解放军国防科学技术大学 | 一种诊断卫星姿态控制系统未知故障的方法 |
WO2018105104A1 (ja) * | 2016-12-09 | 2018-06-14 | 三菱電機株式会社 | 故障リスク指標推定装置および故障リスク指標推定方法 |
CN106598029B (zh) * | 2016-12-21 | 2019-02-26 | 北京交通大学 | 基于性能退化的列控车载设备可靠性预测方法 |
CN108241114B (zh) * | 2016-12-26 | 2020-05-26 | 北京邮电大学 | 混合电路芯片的在线故障检测方法、装置及故障检测电路 |
CN109150564B (zh) * | 2017-06-19 | 2021-09-21 | 中国移动通信集团广东有限公司 | 一种用于小区故障告警的预测方法及装置 |
CN107505922B (zh) * | 2017-08-23 | 2019-11-29 | 北京赛普泰克技术有限公司 | 工厂预警实现方法、装置及系统 |
US11188518B2 (en) * | 2017-09-30 | 2021-11-30 | Siemens Aktiengesellschaft | Method and apparatus for generating fault diagnosis information base of numerical control machine tool |
CN107862375A (zh) * | 2017-10-30 | 2018-03-30 | 北京计算机技术及应用研究所 | 一种两阶段的设备故障诊断方法 |
CN107862126B (zh) * | 2017-11-02 | 2020-11-27 | 中国科学院数学与系统科学研究院 | 一种部件级信息多样性条件下的系统可靠性评估方法 |
CN108563142B (zh) * | 2017-12-11 | 2020-07-24 | 南京航空航天大学 | 一种基于sha的航天器acs任务寿命预测方法 |
CN108304661B (zh) * | 2018-02-05 | 2021-05-07 | 南京航空航天大学 | 基于tdp模型的诊断预测方法 |
CN108470222B (zh) * | 2018-02-08 | 2022-11-11 | 山东蓝光软件有限公司 | 一种用于典型事件预测预报的多征兆增信方法 |
CN108398939B (zh) * | 2018-03-01 | 2020-01-10 | 西北工业大学 | 一种基于ds证据理论的故障诊断方法 |
CN108196514A (zh) * | 2018-03-29 | 2018-06-22 | 武汉理工大学 | 一种数控机床运行状态远程监控方法 |
CN108764520B (zh) * | 2018-04-11 | 2021-11-16 | 杭州电子科技大学 | 一种基于多层循环神经网络和d-s证据理论的水质参数预测方法 |
CN110163403A (zh) * | 2018-04-24 | 2019-08-23 | 北京机电工程研究所 | 基于灰色理论的战术导弹故障预测方法 |
CN108920341B (zh) * | 2018-04-26 | 2021-08-10 | 航天东方红卫星有限公司 | 一种基于蒙特卡洛仿真的小卫星及其星座可用度评估方法 |
CN108681747B (zh) * | 2018-05-11 | 2021-05-04 | 武汉理工大学 | 基于深度学习的旋转机械故障诊断与状态监测系统及方法 |
CN108920426B (zh) * | 2018-07-04 | 2019-08-09 | 西北工业大学 | 一种基于幂均算子和ds证据理论的故障诊断方法 |
CN109189089B (zh) * | 2018-08-29 | 2021-04-13 | 北京控制工程研究所 | 用于卫星地面飞控任务的智能辅助决策方法 |
CN110895625A (zh) * | 2018-09-11 | 2020-03-20 | 湖南银杏可靠性技术研究所有限公司 | 性能退化产品可靠度置信区间估计数值仿真方法 |
CN111064433A (zh) | 2018-10-17 | 2020-04-24 | 太阳能安吉科技有限公司 | 光伏系统故障和警报 |
CN109657835A (zh) * | 2018-10-31 | 2019-04-19 | 中国电力科学研究院有限公司 | 一种配电网区域故障次数预测方法及系统 |
CN110069810B (zh) * | 2019-03-11 | 2023-04-07 | 北京百度网讯科技有限公司 | 电池故障预测方法、装置、设备和可读存储介质 |
CN110135629B (zh) * | 2019-04-24 | 2021-07-02 | 中车青岛四方机车车辆股份有限公司 | 一种列车平均故障周期的预测方法及装置 |
CN110263944A (zh) * | 2019-05-21 | 2019-09-20 | 中国石油大学(华东) | 一种多变量故障预测方法和装置 |
CN110515781B (zh) * | 2019-07-03 | 2021-06-22 | 北京交通大学 | 一种复杂系统状态监测及故障诊断方法 |
CN110503217B (zh) * | 2019-08-29 | 2022-08-26 | 新誉轨道交通科技有限公司 | 一种空调冷媒慢泄漏故障预测方法、装置、设备及系统 |
US11544422B2 (en) * | 2019-09-16 | 2023-01-03 | Palo Alto Research Center Incorporated | Machine learning based systems and methods for real time, model based diagnosis |
CN112290885B (zh) * | 2020-09-30 | 2024-06-07 | 国网浙江安吉县供电有限公司 | 一种户用光伏系统直流侧故障诊断方法 |
CN112329828B (zh) * | 2020-10-26 | 2024-07-23 | 北京旋极信息技术股份有限公司 | 一种故障关联分析方法和装置 |
CN112565593B (zh) * | 2020-11-23 | 2021-10-15 | 中国科学院西安光学精密机械研究所 | 一种针对空间调焦机构的动态可靠性分析方法 |
CN112949683B (zh) * | 2021-01-27 | 2023-02-07 | 东方红卫星移动通信有限公司 | 一种低轨卫星智能故障诊断及预警方法及系统 |
CN113218667B (zh) * | 2021-04-06 | 2022-06-10 | 南京航空航天大学 | 一种超磁致伸缩线控制动系统故障诊断装置及方法 |
CN113255965B (zh) * | 2021-04-26 | 2024-09-13 | 大连海事大学 | 一种雷达发射机退化故障预后智能处理系统 |
CN113391622B (zh) * | 2021-06-22 | 2022-09-02 | 中国人民解放军国防科技大学 | 利用多元多步预测技术的航天器姿态系统异常检测方法 |
CN113626285A (zh) * | 2021-07-30 | 2021-11-09 | 平安普惠企业管理有限公司 | 基于模型的作业监控方法、装置、计算机设备和存储介质 |
CN114154583B (zh) * | 2021-12-08 | 2024-05-24 | 深圳博沃智慧科技有限公司 | 一种小波分析耦合lstm神经网络的水质预测方法 |
CN115131943B (zh) * | 2022-07-07 | 2023-10-31 | 杭州申昊科技股份有限公司 | 一种基于深度学习的声光联动预警方法 |
CN116389229B (zh) * | 2023-06-07 | 2023-08-15 | 南京科羿康光电设备有限公司 | 基于rs485总线的自愈环网络系统 |
CN116956458A (zh) * | 2023-06-26 | 2023-10-27 | 北京空间飞行器总体设计部 | 一种航天器可重构性的表征、判定和量化方法 |
CN118037280B (zh) * | 2024-04-12 | 2024-06-11 | 山东鑫林纸制品有限公司 | 一种瓦楞纸生产线维护与故障诊断系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807047A (zh) * | 2010-03-19 | 2010-08-18 | 北京航空航天大学 | 基于模糊奇偶方程及ar模型的非线性系统故障预测方法 |
-
2011
- 2011-05-31 CN CN 201110144382 patent/CN102208028B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807047A (zh) * | 2010-03-19 | 2010-08-18 | 北京航空航天大学 | 基于模糊奇偶方程及ar模型的非线性系统故障预测方法 |
Non-Patent Citations (2)
Title |
---|
唐建等.小卫星电源系统遥测数据的区间预测技术.《华中科技大学学报( 自然科学版)》.2009,第37卷(第1期), |
小卫星电源系统遥测数据的区间预测技术;唐建等;《华中科技大学学报( 自然科学版)》;20090831;第37卷(第1期);210-212 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI551981B (zh) * | 2014-01-09 | 2016-10-01 | 富士通股份有限公司 | 分析方法、分析程式及分析裝置 |
Also Published As
Publication number | Publication date |
---|---|
CN102208028A (zh) | 2011-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102208028B (zh) | 一种适用于动态复杂系统的故障预测和诊断方法 | |
Kordestani et al. | Failure prognosis and applications—A survey of recent literature | |
CN111666982B (zh) | 一种基于深度神经网络的机电设备故障诊断方法 | |
Heydarzadeh et al. | A two-stage fault detection and isolation platform for industrial systems using residual evaluation | |
US20130116996A1 (en) | Method for integrating models of a vehicle health management system | |
Chang et al. | Remaining useful life prediction for rolling bearings using multi-layer grid search and LSTM | |
Cui et al. | A novel fault measure and early warning system for air compressor | |
CN112082639A (zh) | 发动机状态的诊断方法及其诊断建模方法 | |
Lu et al. | Joint training of a predictor network and a generative adversarial network for time series forecasting: A case study of bearing prognostics | |
Mathew et al. | Regression kernel for prognostics with support vector machines | |
Wang et al. | A remaining useful life prediction model based on hybrid long-short sequences for engines | |
Chen et al. | A novel Bayesian-optimization-based adversarial TCN for RUL prediction of bearings | |
Kefalas et al. | Automated machine learning for remaining useful life estimation of aircraft engines | |
CN111241629B (zh) | 基于数据驱动的飞机液压泵性能变化趋势智能预测方法 | |
CN112488171A (zh) | 一种变分模式分解的sru故障预测方法 | |
Sadoughi et al. | A deep learning approach for failure prognostics of rolling element bearings | |
Lv et al. | A new performance degradation evaluation method integrating PCA, PSR and KELM | |
Dang et al. | seq2graph: Discovering dynamic non-linear dependencies from multivariate time series | |
Wang et al. | Assessing the Performance Degradation of Lithium‐Ion Batteries Using an Approach Based on Fusion of Multiple Feature Parameters | |
Liu et al. | Causality-based PCA Methods for Condition Modeling of Mechatronic Systems | |
CN116933003A (zh) | 基于DaNet的无人机发动机的剩余使用寿命预测方法 | |
Zhou et al. | Degradation indicator construction using dual-class component feature fusion recalibration for bearing performance evaluation | |
Li | Remaining useful life prediction of bearings using fuzzy multimodal extreme learning regression | |
Chen et al. | Remaining useful life interval estimation-based maintenance decision-making for key equipment | |
Jang | A deep learning-based conditional system health index method to reduce the uncertainty of remaining useful life prediction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130619 Termination date: 20190531 |
|
CF01 | Termination of patent right due to non-payment of annual fee |