[0034] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
[0035] 本发明通过实时监测涡喷发动机工作状态，采集振动信号、发动机转速、发动机排气温度、燃油流量、发动机推力5种信号，对采集的到的多种模态下的信号采用深度波尔曼兹机进行特征提取之后，采取数据融合技术，融合多种模态下的状态特征，作为故障诊断的依据。
[0036] 下面结合附图对本发明的应用原理作详细的描述。
[0037] 如图1所示，本发明实施例提供的航空发动机轴承故障诊断方法包括以下步骤：
[0038] S101：通过采集涡喷发动机的振动信号、发动机转速、发动机排气温度、燃油流量、发动机推力5种信号，首先提取各种信号的经典的时域统计学特征、频域特征，以及基于小波包分解获得时频域特征；
[0039] S102：然后分别对各种信号的特征参数采用深度波尔曼兹机进行学习；接下来采用数据融合技术将深度波尔曼兹机学习到的各源特征进行融合；
[0040] S103：最后将融合后的特征参数作为分类器支持向量机的输入，对发动机故障进行分类识别。
[0041] 下面结合附图对本发明的应用原理作进一步的描述。
[0042] 本发明基于多源特征数据融合和深度波尔曼兹机的航空发动机轴承故障诊断方法，流程图如图2所示。具体步骤如下：
[0043] (1)采集振动信号、发动机转速、发动机排气温度、燃油流量、发动机推力5种信号，分别在涡喷航空发动机上安装加速度传感器、转速传感器、温度传感器、流量传感器、压力传感器，在线采集其运行过程中的五种模态下的各种信号；
[0044] (2)提取各种模态下信号特征，对采集到的振动信号等五种信号，提取其时域统计学特征，基于快速傅里叶变换提取频域特征，基于小波包分解提取时频域特征。
[0045] 主要时域特征包括有量纲参数有：平均值，均方根值(有效值)，方差，方根幅值，峰值，峭度，偏度。无量纲参数：波形指标，峰值指标，脉冲指标，裕度指标，峭度指标。
[0046] 频域特征提取包括：应用快速傅里叶变换时域信号，获得频域参数。频域特征参数被考虑，如均值、方差、标准方差、偏度、峭度、中心频率、均方根值。其中均方根值，指的是在指定频带的频率幅值均方根值。
[0047] 基于小波包分解提取时频域特征：每个小波包分解中，分解级数为6，也就是说获得26个小波系数。每个小波系数的能量被计算，每个母小波共计26个小波系数的能量组成一个特征向量，表征故障状态特征。
[0048] (3)分别就每个模态下提取的特征，作为深度神经网络输入，采用深度波尔曼兹机进行学习。图3所示为一个有3个隐藏层的深度波尔曼兹机网络结构。
[0049] 深度波尔曼兹机的能量函数定义如下：
[0050]
[0051] 式中，v表示可视层；h(1)，…，h(L)表示L个隐藏层；模型参数θ＝{W，a，b}，其中W表示网络间的连接权值矩阵，a表示可视层的偏值，b表示隐藏层的偏值；nv表示可视层神经元数，nhl表示第l隐藏层的神经元数。
[0052] 整个网络分配到可视层的概率如下：
[0053]
[0054] 可视层和各隐藏层间的条件分布分别定义如下：
[0055]
[0056]
[0057]
[0058]
[0059] 式 (3) ～ (6) 中， g( · ) 为逻辑函数 。可选取sigmoid函数或者双正切函数。
[0060] 每种模态下，特征学习过程如下：
[0061] a)将特征数据输入第一个受限波尔曼兹机(Restricted Boltzmann machines，RBM)RBM1，假设其可视层为v，采用一步对比散度算法(One-step Contrastive Divergence),基于重构的可视层矢量中场值(mean-field reconstructions of RBM1’s visiblevectors)，学习计算得到隐藏层第一层h(1)、层间自顶向下的连接权值W(1)，层间自底向上的连接权值矩阵为2W(1)。
[0062] b)冻结RBM1的自底向上的连接权值矩阵为2W(1)，按照条件概率 P(h(1)|v；2W(1))抽样，提取特征h(1)作为第一个受限波尔曼兹机RBM2的输入。 RBM2自顶向下和自底向上两个方向的连接权值都为2W(2)，基于重构的可视层矢量中场值，采用一步对比散度算法，训练RBM2。
[0063] c)冻结RBM2的连接权值矩阵为2W(2),按照条件概率P(h(2)|v；2W(1),2W(2))抽样提取特征h(2)，作为RBM3的输入，训练方式同步骤b)。
[0064] d)如图步骤c)进行递归处理，直到L-1层，这里L指深度波尔曼兹机的隐藏层数。
[0065] e)使用一步对比散度算法训练顶层的受限波尔曼兹机RBM，约束条件为：自底向上(L) (L)的连接权重为W ，自顶向下的连接权重为2W 。
[0066] f)使用连接权值{W(1)，W(2)，…,W(L)},构建一个深度波尔曼兹机。
[0067] g)利用构建好的深度波尔曼兹机进行特征提取。
[0068] (4)采用多特征参数数据融合算法，对步骤(3)提取的各种模态特征进行特征融合。
[0069] (5)基于步骤(4)融合的特征参数，采用支持向量机进行故障分类识别。
[0070] 下面结合具体实施例对本发明的应用效果作详细的描述。
[0071] 本发明的应用试验平台，该试验平台为20公斤推力的通用航空涡喷发动机，采用发明的方法，对发动机轴承故障进行诊断。试验方案如下：
[0072] (1)采集信号
[0073] (a)振动加速度信号：a(t)，用振动加速度传感器(1000mV/g)测量，该传感器采用透明胶带缠绕在发动机上面。发动机表面应该磨平(因为加速度传感器的底部是平的)。采样频率：200kHz。
[0074] (b)温度信号：T(t)，由厂家提供，输出方式0-5V，采样频率：1Hz。
[0075] (c)转速信号：r(t)，由厂家提供，输出方式0-5V，采样频率：200kHz。
[0076] (d)推力信号：F(t)。
[0077] (e)油耗信号：O(t)。
[0078] (2)发动机轴承故障诊断
[0079] (a)正常部件下，设定转速为6000、10000、12000rpm，采集信号，长度为 10s，采集20次。
[0080] (b)轴承内圈故障(直径分别为0.05mm、0.1mm、0.2mm、0.3mm、0.5mm)，设定转速为6000、10000、12000rpm，采集信号，长度为10s，采集20次。
[0081] (c)轴承外圈故障(直径分别为0.05mm、0.1mm、0.2mm、0.3mm、0.5mm)，设定转速为6000、10000、12000rpm，采集信号，长度为10s，采集20次。
[0082] (d)轴承球故障(直径分别为0.05mm、0.1mm、0.2mm、0.3mm、0.5mm)，设定转速为6000、10000、12000rpm，采集信号，长度为10s，采集20次。
[0083] 通过分析本发明提出的一种基于多源特征数据融合和深度波尔曼兹机的航空发动机轴承故障诊断方法在20公斤推力的涡喷发动机上应用的试验结果，可以得出如下结论：
[0084] (1)本发明采用多源特征融合数据，较之单源特征状态数据，效果明显，提升了涡喷发动机的轴承故障的识别精度和稳定性。
[0085] (2)区别于传统的时域、频域、时频域的信号特征提取方法，本发明利用深度波尔曼兹机方法直接对各种模态原始数据进行重构，在保留数据原始信息的同时，有效的挖掘了数据的本质特征。
[0086] (3)本发明将无监督学习到的多源特征，进行数据融合之后，输入到支持向量机模型，并对深度波尔曼兹机进行微调，极大的提高了涡喷发动机轴承故障精度。
[0087] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
