基于回归分析和RBF神经网络的退化数据缺失插补方法技术

本发明专利技术公开了一种基于回归分析和RBF神经网络的退化数据缺失插补方法，包括以下几个步骤：步骤一、已观测的退化数据趋势回归建模；步骤二、计算已观测退化数据的残差序列；步骤三、建立RBF神经网络，并利用已观测数据的残差序列训练网络；步骤四、通过训练好的RBF神经网络估计缺失数据的残差序列；步骤五、合并缺失数据的趋势项与残差序列的估计结果为退化数据插补结果。本发明专利技术将回归分析方法与RBF神经网络方法结合，提出了一种退化数据缺失插补方法，解决了加速退化试验中缺失性能退化数据的插补问题。

【技术实现步骤摘要】
基于回归分析和RBF神经网络的退化数据缺失插补方法
本专利技术涉及一种基于回归分析RBF神经网络的退化数据缺失插补方法，属于加速退化试验

技术介绍
由于检测仪器故障、记录人员过失等原因，常常会造成在加速退化试验中采集到的性能退化数据缺失的情况。而在加速退化试验的数据处理与评估、故障预测或寿命预测中，都需要完整的数据作为输入。数据缺失给性能退化数据的处理带来了很大的困难，很多传统的性能退化数据处理方法无法对有缺失数据的数据进行统计分析，比如一些关于时间序列数据分析的方法就要求所分析的数据为完整的等间距数据集。在故障预测或寿命预测中，性能退化数据的缺失会造成预测结果的偏差。这时要么改进传统的数据处理方法，使之能够适应有缺失的性能退化数据，要么针对有缺失的数据进行处理，使之转化为完整的数据。前一种方案，在实现上有难度，后一种方案的实现更切合实际，且有缺失的数据转化为完整的数据后可以应用到更多的数据处理方法中，适用性更强。缺失数据的插补方法是如今应用最多、发展最快的缺失数据处理方法。采用有效的方法处理缺失数据可以提高对数据资源的有效利用，减少数据资源的浪费；并会有助于性能退化数据的处理，有助于加速退化试验评估，有助于故障预测或寿命预测等工作，甚至可能影响研究结论，减少如寿命预测结果不准确等情况所带来的隐患。因此，对有缺失的性能退化数据采用正确和有效的方法进行处理是能否得出正确结论的关键之一。常用的缺失数据插补法有均值插补法、热卡插补法、冷卡插补法、最邻近插补法、回归插补法、EM（ExpectationMaximization）算法等。这些方法有一个共同的缺点：扭曲样本分布。比如均值插补法会降低变量间的相关关系，回归插补法则会人为的增大变量间的相关关系，而且这些插补法忽略了缺失数据插补的不确定性。尽管由于回归插补引入随机误差项，能够缓解这一问题，但是随机误差项的确定是比较困难的。回归分析插补是一种非常常用的插补方法。回归分析插补的基本思想是利用辅助变量与已观测数据的关系，建立回归模型，利用已知的辅助变量的信息，对缺失的数据进行估计。回归分析插补通常使用服从正态分布的残差作为随机项，以使插补数据保持原始数据的离散性，但有些情况下正态分布随机项并不能够反映原始数据的离散性。径向基函数（RadialBasisFunction,RBF）神经网络是目前应用最广泛的神经网络模型之一，相比于以BP神经网络拓扑结构为基础的小波神经网络，RBF神经网络在函数逼近能力、分类能力和学习速度等方面均有优势。RBF神经网络是以函数逼近理论为基础而构造的一类前向网络，这类网络的学习等价于在多维空间中寻找训练数据的最佳拟合平面，RBF神经网络的每个隐层神经元的函数都构成了拟合平面的一个基函数。RBF神经网络结构简单、训练简洁而且学习收敛速度快，能够逼近任意非线性函数。
技术实现思路
本专利技术的目的是为了解决性能退化数据的缺失插补问题，提出了一种通用性较强的基于回归分析和RBF神经网络的退化数据缺失插补方法。本专利技术综合采用回归分析方法与RBF神经网络方法，利用已观测退化数据的信息来实现缺失退化数据的插补。本专利技术是一种基于回归分析和RBF神经网络的退化数据缺失插补方法，包括以下几个步骤：步骤一、已观测的退化数据趋势回归建模；步骤二、计算已观测退化数据的残差序列；步骤三、建立RBF神经网络，并利用已观测退化数据的残差序列训练网络；步骤四、通过训练好的RBF神经网络估计缺失数据的残差序列；步骤五、合并缺失数据的趋势项与缺失数据的残差序列的估计结果为退化数据插补结果。本专利技术的优点在于：（1）对已观测退化数据的趋势进行回归建模可以使插补数据的趋势与原始数据保持一致。（2）使用RBF神经网络估计缺失数据的残差序列可以使插补数据保持原始数据的离散性，使插补数据更加接近真实。（3）RBF神经网络可以应对不同类型的残差，适用范围广泛。附图说明图1是本专利技术所述基于回归分析和RBF神经网络的退化数据缺失插补方法的流程图；图2是本专利技术实施例有缺失的退化数据曲线；图3是本专利技术实施例缺失数据的退化趋势估计结果；图4是本专利技术实施例缺失数据的残差序列估计结果；图5是本专利技术实施例缺失数据的最终插补结果。具体实施方式下面将结合附图和实施例对本专利技术作进一步的详细说明。本专利技术提供一种基于回归分析和RBF神经网络的退化数据缺失插补方法，所述方法针对有缺失的性能退化数据进行插补处理，假设完整的退化数据为Y，其对应的时间为T，记Y=(Yobs,Ymis)，T=(Tobs,Tmis)，其中已观测退化数据Yobs、已观测退化数据对应的时间Tobs和缺失数据对应的时间Tmis为已知数据，具体为：式中，n为已观测退化数据的数据量，m为缺失数据的数据量。通过已知数据，综合利用回归分析方法与RBF神经网络方法来估计缺失数据Ymis：Ymis＝(ymis_1,ymis_2,…,ymis_m)(2)并最后得到插补完整的退化数据Y=(Yobs,Ymis)。本专利技术是一种基于回归分析和RBF神经网络的退化数据缺失插补方法，方法流程图1所示，包括以下几个步骤：步骤一、已观测退化数据趋势回归建模；根据已观测退化数据的趋势特点选择回归函数，本专利技术中选择指数函数作为回归函数。再根据已观测退化数据Yobs=(yobs_1，yobs_2，…，yobs_n)与对应的时间Tobs=(tobs_1，tobs_2，…，tobs_n)，利用最小二乘法来估计回归模型参数，得到退化趋势的函数表达式：式中，a、b与c为回归模型参数，tobs_i表示已观测退化数据对应的时间。通过回归分析得到的退化趋势模型f(t)，将缺失数据对应的时间Tmis=(tmis_1，tmis_2，…，tmis_m)作为输入，计算缺失数据的趋势序列Qmis=(qmis_1,qmis_2,…,qmis_m)：步骤二、计算已观测退化数据的残差序列；将已观测退化数据对应的时间Tobs作为输入，通过趋势模型f(t)计算已观测退化数据的趋势序列Qobs=(qobs_1,qobs_2,…,qobs_n)，并用相应的已观测退化数据的真实值Yobs减去已观测数据的趋势序列Qobs，得到已观测退化数据的残差序列Eobs=(eobs_1,eobs_2,…,eobs_n)：eobs_i＝yobs_i-qobs_i,i＝1,2,…,n(5)步骤三、建立RBF神经网本文档来自技高网...

【技术保护点】
一种基于回归分析和RBF神经网络的退化数据缺失插补方法，其特征在于，包括以下几个步骤：步骤一、已观测的退化数据趋势回归建模；根据已观测退化数据的趋势特点选择回归函数，再根据已观测退化数据Yobs=(yobs_1，yobs_2，…，yobs_n)与对应的时间Tobs=(tobs_1，tobs_2，…，tobs_n)，利用最小二乘法来估计回归模型参数，得到退化趋势的函数表达式；通过回归分析得到的退化趋势模型f(t)，将缺失数据数据对应的时间Tmis=(tmis_1，tmis_2，…，tmis_m)作为输入，计算缺失数据的趋势序列Qmis=(qmis_1,qmis_2,…,qmis_m)；步骤二、计算已观测退化数据的残差序列；将已观测数据对应的时间Tobs作为输入，通过趋势模型f(t)计算已观测数据的趋势序列Qobs=(qobs_1,qobs_2,…,qobs_n)，并用已观测退化数据的真实值Yobs减去已观测数据的趋势序列Qobs，得到已观测退化数据的残差序列Eobs=(eobs_1,eobs_2,…,eobs_n)：eobs_i＝yobs_i‑qobs_i,i＝1,2,…,n；步骤三、建立RBF神经网络，并利用已观测数据的残差序列训练网络；建立单输入单输出的RBF神经网络：y=Σi=1kwiRi(x)]]>式中，y表示一维输出向量，x表示一维输入向量，wi是隐层与输出层间权值，k是感知单元的个数，Ri(x)表示基函数。将已观测退化数据的残差序列Eobs作为输出向量，相应的时间Tobs作为输入向量，采用RBF神经网络的学习算法对网络进行训练，得到基函数的中心ci和方差σi，以及权值wi；步骤四、通过训练好的RBF神经网络估计缺失数据的残差序列；将缺失数据对应的时间Tmis作为输入，通过训练好的RBF神经网络估计缺失数据的残差序列Emis=(emis_1,emis_2,…,emis_m)：emis_i=Σj=1kwjRj(tmis_i)=Σj=1kwjexp[-||tmis_i-ci||2σj2],i=1,2,···,m]]>在估计缺失数据残差序列的过程中，不断更新RBF神经网络的训练数据，将已经估计得到的残差序列值emis_i与对应的时间tmis_i补充到训练数据{Eobs,Tobs}中，通过新训练数据训练得到的RBF神经网络去估计下一个缺失数据的残差值emis_i+1；这样不断更新训练数据再估计，直到估计完所有缺失数据的残差序列值；步骤五、合并缺失数据的趋势项与残差序列的估计结果为退化数据插补结果；将由步骤一中得到的缺失数据的趋势序列Qmis与由步骤四中得到的缺失数据残差序列Emis合并，得到最终的缺失数据插补结果Ymis=(ymis_1，ymis_2，…，ymis_m)：ymis_i＝qmis_i+emis_i,i＝1,2,…,m这样最后得到插补完整的性能退化数据Y=(Yobs,Ymis)，完成了退化数据的缺失插补工作。...

【技术特征摘要】
2013.12.13 CN 201310684436.21.一种基于回归分析和RBF神经网络的退化数据缺失插补方法，其特征在于，包括以下几个步骤：步骤一、已观测的退化数据趋势回归建模；根据已观测退化数据的趋势特点选择回归函数，再根据已观测退化数据Yobs＝(yobs_1，yobs_2，…，yobs_n)与对应的时间Tobs＝(tobs_1，tobs_2，…，tobs_n)，利用最小二乘法来估计回归模型参数，得到退化趋势的函数表达式；式中，a、b与c为回归模型参数；通过回归分析得到的退化趋势模型f(t)，将缺失数据数据对应的时间Tmis＝(tmis_1，tmis_2，…，tmis_m)作为输入，计算缺失数据的趋势序列Qmis＝(qmis_1,qmis_2,…,qmis_m)；步骤二、计算已观测退化数据的残差序列；将已观测数据对应的时间Tobs作为输入，通过趋势模型f(t)计算已观测数据的趋势序列Qobs＝(qobs_1,qobs_2,…,qobs_n)，并用已观测退化数据的真实值Yobs减去已观测数据的趋势序列Qobs，得到已观测退化数据的残差序列Eobs＝(eobs_1,eobs_2,…,eobs_n)：eobs_i＝yobs_i-qobs_i,i＝1,2,…,n；步骤三、建立RBF神经网络，并利用已观测数据的残差序列训练网络；建立单输入单输出的RBF神经网络：式中，y表示一维输出向量，x表示一维输入向量，wi是隐层与输出层间权值，k是感知单元的个数，Ri(x)表示基函数；将已观测退化数据的残差序列E...

【专利技术属性】
技术研发人员：孙富强，范晔，李晓阳，姜同敏，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11