1.基于退化数据的动量轮可靠性分析与剩余寿命预测方法，具体按照以下步骤实施：步骤1、分析动量轮的性能退化机理，确定关键性能退化参数，并对退化数据进行预处理；

步骤2、分析步骤1中得出的退化参数的特点，建立合适的性能退化数据模型；

步骤3、根据步骤2得到的退化模型，利用极大似然法估计出模型中的未知参数；

步骤4、分析步骤2退化模型的特点以及性质，再根据模型的性质推出关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数；

步骤5、将步骤3估计出的参数值代入步骤4中的各个函数表达式中，得出函数的具体表达式，从而根据函数分析动量轮各个测量时刻的可靠度以及剩余寿命。

2.根据权利要求1所述的基于退化数据的动量轮可靠性分析与剩余寿命预测方法，其特征在于，步骤2中对动量轮的退化建模具体步骤如下：

2.1)将动量轮中的轴系润滑剂作为关键性能参数，通过观测退化数据发现，润滑剂剩余重量具有单调减小的趋势，且具有波动；

2.2)根据维纳过程的性质可知，维纳过程不具有严格单调性，适合对具有单调波动的性能退化过程进行建模且该过程已被广泛的用于退化建模中，而动量轮退化数据正符合维纳过程的性质；

2.3)根据步骤2.2)的分析，建立如下带有测量误差的退化数据模型：

式(1)中：a为漂移参数，满足 X(t)为不含误差的退化量，Y(t)为带有误差的退化量，B(t)用于表征退化过程随机性的标准布朗运动过程，δB为扩散参数，用于表示同类产品的共性，ε为随机测量误差，满足 进一步假设a,B(t)和ε是相互独立的。

3.根据权利要求1所述的基于退化数据的动量轮可靠性分析与剩余寿命预测方法，其特征在于，步骤3中对模型参数估计的具体过程如下：假设Yi(tj)表示第i个动量轮在第j时刻测量得到的性能退化量，则Yi(tj)＝aitj+δBB(tj)+δεεij，1＜i＜n,1＜j＜m  (2)令t＝(t1,t2...tm)',Yi＝(Yi(t1),Yi(t2)...Yi(tm))',Y＝(Y1',...Yn')'，其中，Yi服从维正态分布，其均值和方差为μ0＝μat，∑＝δa2tt'+Ω  (3)其中

Ω＝δB2Q+δε2P

根据极大似然原理写出对数似然函数为：

其中

对式(4)关于μa和δa分别求一阶偏导，并分别令其等于零，得：

将式(5)和(6)中的 代入式(4)中，得到：

通过极大化式(7)，得到δB和δε的极大似然估计值，再将δB和δε的估计值代入式(5)和(6)中得到μa和δa的极大似然估计值，最后估计出模型中的所有未知数。

4.根据权利要求1所述的基于退化数据的动量轮可靠性分析与剩余寿命预测方法，其特征在于，步骤4中关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数各自的具体表达式求解过程如下：

4.1)根据首达时间(First Hitting Time,FHT)的概念来定义寿命，即如果退化过程{X(t),t≥0}首次等于或超过一个预设的失效阈值Df时，认为系统失效，寿命T可以定义为：T＝inf{t:X(t)≥Df|X(0)＜Df}  (8)为不失一般性，令初始退化量X(t)＝0，失效阈值Df＞0；根据维纳过程的性质以及首达时间分布的特点可知，寿命T服从逆高斯(Inverse Gaussian,IG)分布；

a)不考虑漂移a的随机性时，在一元维纳过程{X(t),t≥0}的基础上定义随机过{Z(t),t≥0}为

即任意时刻t≥0，Z(t)取X(t)在时间[0,t]的最大值，取t时刻Z(t)的概率密度函数g(z,t)，由{Z(t),t≥0}的定义可知其是单调随机过程，产品的不失效概率为：

利用Fokker-Planck(Kolmogorov方程)方程即求得g(z,t)的形式为：

将式(11)带入式(10)得：

通过1-F(t)＝R(t)的关系得出寿命T的分布函数：

对式(13)关于时间t求导，得出寿命T的概率密度函数为：

b)考虑漂移系数a为随机变量的时，通过如下的全概率公式得出寿命T的概率密度函数：fT(t)＝∫afT|a(t|a)p(a)da＝Ea[fT|a(t|a)]  (15)其中p(a)为a的概率密度函数，fT|a(t|a)不考虑a的随机性时寿命T概率密度函数，Ea[·]是关于参数a的数学期；

4.2)为了简化推导过程，给出定理：若a～N(μ,δ2)且满足a,b,c∈R+，则

证明：

通过步骤4.1)、4.2)得出考虑漂移系数a为随机变量的时寿命T的概率密度函数为：

4.3)利用 关系，得出可靠性函数为：

4.4)对于当前时刻tk，同样根据首达时间的定义动量轮剩余寿命定义为：Lk＝inf{lk＞0:X(lk+tk)≥Df}  (20)式中，lk为tk时刻动量轮的剩余寿命，根据上述剩余寿命的定义以及漂移维纳过程的性质可知：

1)在不考虑a的随机性时，当获得tk时刻的测量值X(tk)，因为t＞tk，退化过程记作X(t)＝X(tk)+a(t-tk)+δBB(t-tk)，当t是{X(t),t≥tk}的首达时间分布时，lk＝t-tk为tk时刻的剩余寿命且lk＞0，将退化过程变换为：X(lk+tk)-X(tk)＝alk+δBB(lk)  (21)进一步，将tk时刻的剩余寿命记作{Z(lk),lk≥0}通过失效阈值Dfk＝Df-X(tk)的首达时间，记Z(lk)＝X(lk+tk)-X(tk)，且Z(0)＝0，则得如下变换式：Z(lk)＝alk+δBB(lk)  (22)此时参考不考虑a的随机性时的动量轮寿命T的概率密度函数的推理过程可得出tk时刻的剩余寿命概率密度函数为：

2)在考虑a的随机性时，根据步骤a)、4.2)得当前tk时刻的剩余寿命概率密度函数为：