发明名称 基于回归分析和RBF神经网络的退化数据缺失插补方法
摘要 本发明公开了一种基于回归分析和RBF神经网络的退化数据缺失插补方法,包括以下几个步骤:步骤一、已观测的退化数据趋势回归建模;步骤二、计算已观测退化数据的残差序列;步骤三、建立RBF神经网络,并利用已观测数据的残差序列训练网络;步骤四、通过训练好的RBF神经网络估计缺失数据的残差序列;步骤五、合并缺失数据的趋势项与残差序列的估计结果为退化数据插补结果。本发明将回归分析方法与RBF神经网络方法结合,提出了一种退化数据缺失插补方法,解决了加速退化试验中缺失性能退化数据的插补问题。
申请公布号 CN103793613A 申请公布日期 2014.05.14
申请号 CN201410061310.4 申请日期 2014.02.24
申请人 北京航空航天大学 发明人 孙富强;范晔;李晓阳;姜同敏
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京永创新实专利事务所 11121 代理人 姜荣丽
主权项 1.一种基于回归分析和RBF神经网络的退化数据缺失插补方法,其特征在于,包括以下几个步骤:步骤一、已观测的退化数据趋势回归建模;根据已观测退化数据的趋势特点选择回归函数,再根据已观测退化数据Y<sub>obs</sub>=(y<sub>obs_1</sub>,y<sub>obs_2</sub>,…,y<sub>obs_n</sub>)与对应的时间T<sub>obs</sub>=(t<sub>obs_1</sub>,t<sub>obs_2</sub>,…,t<sub>obs_n</sub>),利用最小二乘法来估计回归模型参数,得到退化趋势的函数表达式;通过回归分析得到的退化趋势模型f(t),将缺失数据数据对应的时间T<sub>mis</sub>=(t<sub>mis_1</sub>,t<sub>mis_2</sub>,…,t<sub>mis_m</sub>)作为输入,计算缺失数据的趋势序列Q<sub>mis</sub>=(q<sub>mis_1</sub>,q<sub>mis_2</sub>,…,q<sub>mis_m</sub>);步骤二、计算已观测退化数据的残差序列;将已观测数据对应的时间T<sub>obs</sub>作为输入,通过趋势模型f(t)计算已观测数据的趋势序列Q<sub>obs</sub>=(q<sub>obs_1</sub>,q<sub>obs_2</sub>,…,q<sub>obs_n</sub>),并用已观测退化数据的真实值Y<sub>obs</sub>减去已观测数据的趋势序列Q<sub>obs</sub>,得到已观测退化数据的残差序列E<sub>obs</sub>=(e<sub>obs_1</sub>,e<sub>obs_2</sub>,…,e<sub>obs_n</sub>):e<sub>obs_i</sub>=y<sub>obs_i</sub>-q<sub>obs_i</sub>,i=1,2,…,n;步骤三、建立RBF神经网络,并利用已观测数据的残差序列训练网络;建立单输入单输出的RBF神经网络:<maths num="0001"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>w</mi><mi>i</mi></msub><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中,y表示一维输出向量,x表示一维输入向量,w<sub>i</sub>是隐层与输出层间权值,k是感知单元的个数,R<sub>i</sub>(x)表示基函数。将已观测退化数据的残差序列E<sub>obs</sub>作为输出向量,相应的时间T<sub>obs</sub>作为输入向量,采用RBF神经网络的学习算法对网络进行训练,得到基函数的中心c<sub>i</sub>和方差σ<sub>i</sub>,以及权值w<sub>i</sub>;步骤四、通过训练好的RBF神经网络估计缺失数据的残差序列;将缺失数据对应的时间T<sub>mis</sub>作为输入,通过训练好的RBF神经网络估计缺失数据的残差序列E<sub>mis</sub>=(e<sub>mis_1</sub>,e<sub>mis_2</sub>,…,e<sub>mis_m</sub>):<maths num="0002"><![CDATA[<math><mrow><msub><mi>e</mi><mrow><mi>mis</mi><mo>_</mo><mi>i</mi></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><msub><mi>R</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>mis</mi><mo>_</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><mrow><mo>|</mo><mo>|</mo><msub><mi>t</mi><mrow><mi>mis</mi><mo>_</mo><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>c</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mrow><mn>2</mn><msup><msub><mi>&sigma;</mi><mi>j</mi></msub><mn>2</mn></msup></mrow></mfrac><mo>]</mo><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>m</mi></mrow></math>]]></maths>在估计缺失数据残差序列的过程中,不断更新RBF神经网络的训练数据,将已经估计得到的残差序列值e<sub>mis_i</sub>与对应的时间t<sub>mis_i</sub>补充到训练数据{E<sub>obs</sub>,T<sub>obs</sub>}中,通过新训练数据训练得到的RBF神经网络去估计下一个缺失数据的残差值e<sub>mis_i+1</sub>;这样不断更新训练数据再估计,直到估计完所有缺失数据的残差序列值;步骤五、合并缺失数据的趋势项与残差序列的估计结果为退化数据插补结果;将由步骤一中得到的缺失数据的趋势序列Q<sub>mis</sub>与由步骤四中得到的缺失数据残差序列E<sub>mis</sub>合并,得到最终的缺失数据插补结果Y<sub>mis</sub>=(y<sub>mis_1</sub>,y<sub>mis_2</sub>,…,y<sub>mis_m</sub>):y<sub>mis_i</sub>=q<sub>mis_i</sub>+e<sub>mis_i</sub>,i=1,2,…,m这样最后得到插补完整的性能退化数据Y=(Y<sub>obs</sub>,Y<sub>mis</sub>),完成了退化数据的缺失插补工作。
地址 100191 北京市海淀区学院路37号