1.基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：包括以下步骤：
(1)行星齿轮箱不同于传统定轴齿轮箱，根据其结构的特点及诊断的难度，选取合适的振动测点布置位置和传感器类型，制定出合理的行星齿轮箱振动信号采集方案；
(2)根据阶比重构技术，将处于变工况的风电行星齿轮箱传感器所采集的振动信号进行预处理，将非线性、非平稳的时域信号转化为具有平稳性的角域信号，避免使用硬件方式实现等角度采样的昂贵成本；采用线性插值方法的非平稳振动时域信号的阶比重构技术，将等时间间隔采样的非平稳振动时域信号转化为具有平稳特性的角域振动信号，保证行星齿轮箱振动角域信号的整周期性；
(3)多传感器信息的融合；传感器振动信号中含有重要的故障诊断信息，而无量纲指标可以标征行星齿轮箱的运行状态，不会受工况变化的影响，在步骤(2)的基础上，对传感器所采集的振动信号经阶比重构后的角域信号提取出传统无量纲指标和新无量纲指标，作为故障诊断的特征参数；
(4)根据行星齿轮箱的结构特点，将行星齿轮箱的故障按分布故障和局部故障进行划分，通过行星齿轮箱的结构和故障机理分析，推出不同故障模式的故障特征频率，对行星齿轮箱处于正常、一级行星轮系太阳轮局部故障、一级行星轮系太阳轮分布故障三个状态下的振动信号进行简单介绍；
(5)故障模式识别的方法；由于行星齿轮箱的结构复杂，简单的故障诊断方法难以进行精确诊断；传统的BP神经网络广泛的应用于齿轮的故障诊断中，但在输入信息不精确或不确定时，神经网络的精确性大大降低，因此采用遗传算法优化BP神经网络作为行星齿轮箱故障诊断方法；
(6)在步骤(2)、(3)、(4)、(5)的基础上，以华锐SL1500风电行星齿轮箱为研究对象，进行验证，表明：不同传感器测点包含的故障信息不够完全，很难以单一传感器测点的故障信息进行行星齿轮箱的故障诊断，多传感器信息融合的方法较单一传感器的方法的诊断效果更好，准确率更高。
2.根据权利要求1所述的基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：步骤(1)所述的行星齿轮箱振动信号采集方案：
行星齿轮箱一般由太阳轮、行星轮、齿圈和行星架组成；一般情况下，齿圈固定不动，太阳轮绕自身的中心轴线旋转，行星轮不仅自转，还围绕太阳轮公转；行星轮既和太阳轮 啮合，又和齿圈啮合；行星齿轮箱中多个齿轮的复合运动造成了振动信号的复杂性；在振动的监测中，传感器一般安装在齿圈或与之相连的箱体上采集振动信号，太阳轮-行星轮以及行星轮-齿圈啮合副的啮合点对传感器的位置随行星架旋转变化，使得啮合点至传感器之间的振动传递路径发生变化，这种时变的传递路径对振动测试信号产生调幅调制效应；
以华锐SL1500行星齿轮箱为研究对象，为了测得行星齿轮箱运作状态的全部信息，在行星齿轮箱布置了四处传感器测点，分别为一级行星轮大齿圈处、二级行星轮大齿圈处以及输出低速端轴承座处、输出高速端轴承座处，4个传感器测点所测得的径向振动信号包含齿轮箱的全部振动信息，其中传感器所测振动信号的长度为齿轮箱输入轴转一圈时所测得的振动信号。
3.根据权利要求1所述的基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：步骤(2)所述的振动信号进行预处理：
行星齿轮箱处于变转速、变工况的工作环境下，其采集的振动信号为非平稳信号，这对齿轮箱的故障诊断产生很大的困难，采用阶比重采样的方法对振动信号进行角域重采样，将非平稳的时域振动信号转化为平稳或准平稳特性的角域信号，便于对振动信号的分析；
阶次分析技术的核心在于获得相对参考轴的恒定角增量采样数据，因此需要能准确获得阶次采样的时刻及相应的基准转速，即实现阶次跟踪；常见的阶次跟踪方法有硬件阶次跟踪法、计算阶次跟踪法和基于瞬时频率估计的阶次跟踪法；这里采用计算阶次跟踪法实现振动信号的重采样计算；实际的齿轮箱振动信号一般情况下都含有多种干扰成分，这就使得其故障特征的提取变得比较困难；EMD分解方法可以根据信号的局部时变特性，自适应的将任意一个复杂信号分解为一系列分量，通过相关系数法则对信号进行重构，剔除原始信号中的干扰成分；采用阶比分析与EMD分解方法可以有效的对振动信号进行预处理，对下一步的分析做好了准备。
4.根据权利要求1所述的基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：步骤(3)所述的提取出传统无量纲指标和新无量纲指标，作为故障诊断的特征参数：
无量纲诊断是一种将无量纲参数用于设备故障诊断的技术方法，与有量纲幅域诊断参数(如方根幅值、平均幅值)不同，无量纲幅域诊断参数对故障有一定的敏感，而对信号的幅值和频率的变化不敏感，即和监测对象的工况关系不大；这里从预处理过的振动信号中提取以下时域无量纲指标；
无量纲指标定义如下：
式(1)中,x表示振动幅值，p(x)表示振动幅值的概率密度函数； 为积分符号，dx为微分符号；l、m为常数，不同无量纲指标值会有所不同；目前，最常用的无量纲指标有波形指标、脉冲指标、裕度指标、峰值指标和峭度指标；根据上式有：
若l＝2，m＝1，则有波形指标：
式(2)中：x表示振动幅值，p(x)表示振动幅值的概率密度函数； 为积分符号，dx为微分符号；
若l＝∞，m＝1，则有脉冲指标：
式(3)中：x表示振动幅值，p(x)表示振动幅值的概率密度函数； 为积分符号，dx为微分符号；Xrms表示均方根平均值，lim表示极限符号， 为平均幅值；
若l＝∞，m＝1/2，则有裕度指标：
式(4)中：x表示振动幅值，p(x)表示振动幅值的概率密度函数； 为积分符号，dx为微分符号；lim表示极限符号，Xmax为最大值；Xr为方根幅值；
若l＝∞，m＝2，则有峰值指标：
式(5)中：x表示振动幅值，p(x)表示振动幅值的概率密度函数； 为积分符号，dx为微分符号；Xrms表示均方根平均值，lim表示极限符号，Xmax为最大值；
若l＝4，m＝2，则有定义峭度指标：
式(6)中：x表示振动幅值，p(x)表示振动幅值的概率密度函数； 为积分符号，dx为微分符号；Xrms表示均方根平均值，β为峭度；
计算信号y＝ax的无量纲指标：
其中: d(ax)＝ad(x)；
式(7)中：x表示振动幅值，a为常数，随机变化，y表示信号幅值变化；p(y)表示振动幅值的概率密度函数,dy为微分符号, 为积分符号，l、m为常数，不同无量纲指标值会有所不同；
也就是说，无量纲指标不随随机过程的幅值变化，即不会受工作状况变化的影响，它只对概率密度函数p(x)的形状变化敏感；无量纲指标的优点具体归纳为：1)完全具有反应故障特征的能力；2)几乎与振动信号绝对水平无关；3)对不同类型故障存在不同敏感度；4)基本不受工况、载荷、转速等变化的影响；
上述传统无量纲特征参数在一定程度上仍受能量和工况变化的影响，为了克服该影响，基于整周期采样新的无量纲特征参数被提出，主要包括重复性因子、相似性因子以及跳跃性因子；对于整周期采样n×m则可得数据为：
{x11,x12,…x1m；…xn1,xn2,…xnm}       (8)；
式(8)中：x表示振动幅值，{xn1,xn2,…xnm}表示第n段长度为m的数据，其中n＝1,2,
3,…，m为常数，随具体周期采样而定；
计算上述数据的各段差分为：
{Δx11,…Δx1m-1；…；Δxn1,…Δxnm-1}         (9)；
式(9)中：{Δxn1,Δxn2,…Δxnm-1}表示各段整周期采样数据的差分序列；n＝1,2,3,…，m为常数；
计算一段平均重复波形：
式(10)中： 表示m个整周期采样数据相同序列号所对应数据的平均值，共有m个；
再计算重复性波形平均差分：
式(11)中：Δxm表示每列数据的平均值与平均重复波形之间的差值，m＝1,2,3,…；
为了表示出波形的趋势，规定下降为0，不变为1，上升为2，则k表示与平均重复波形不同的点，则重复性因子为：
式(12)中：k表示与平均重复波形不同的点的个数；m表示整周期采样数据的个数；
相似性因子是由分形维度得到的，为了降低幅值变化的敏感性，首先将数据标准化 在这里，k是一个大于1的整数；计算分形锥度—盒维数，盒维数要求覆盖单元具
有自相似性，并要求曲线具有严格的自相似性，基于这些原理该方法已经在振动中得到应用；设F是Rn中任意一非空有限子集，记N(F,δ)表示最大直径δ为且能覆盖F集合最小数，则F的盒维数定义为的盒维数定义：
式(13)中：N(F,δ)表示最大直径δ为且能覆盖F集合最小数，In表示对数，δ表示最大直径；
最常用N(F,δ)的取法是将F划分网格，取相交于F的边长为δ的网络块数；可以在In(1/δ)～InN(F,δ)关系图中，确定线性好的一段直线，拟合该段直线的斜率就是对应的分形盒维数；相似性因子为
Ff＝dimBF       (14)；
式(14)中：F表示Rn中任意一非空有限子集，dimB表示直线的斜率；
与相似性因子一样，对所得数据进行标准化，对获取的标准化后的分段波形数据求方差，取跳跃性因子为方差值：
式(15)中：x表示振动幅值；xip标准化后的分段波形数据； 表示振动数据的平均值，n为振动数据的长度。
5.根据权利要求1所述的基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：步骤(4)所述的将行星齿轮箱的故障按分布故障和局部故障进行划分，通过行星齿轮箱的结构和故障机理分析：
风电行星齿轮箱一般分为行星轮系、平行级齿轮，对于行星轮系和平行级齿轮，其故障可分为局部故障和分布式故障；在单级行星齿轮箱中，太阳轮–行星轮和行星轮–齿圈两种啮合副的啮合频率相同；通常齿圈固定不动，太阳轮、行星轮和行星架旋转，在这种情况下，啮合频率：
式(16)中：Zr和Zs分别为齿圈和太阳轮的齿数；fm为啮合频率；fc为行星架的旋转频率；
为太阳轮的绝对旋转频率；
太阳轮局部故障特征频率为：
式(17)中：fm为啮合频率；Zs为太阳轮齿数；fs为太阳轮局部故障特征频率；N为行星轮数量；
行星轮局部故障特征频率为：
式(18)中：fm为啮合频率；Zp为行星轮齿数，fp为行星轮局部故障特征频率；齿圈局部故障特征频率为：
式(19)中：fm为啮合频率；fr为齿圈局部故障特征频率；N为行星轮数量，Zr为齿圈齿数；
行星齿轮箱中各种齿轮的分布式故障特征频率等于齿轮相对行星架(太阳轮和齿圈故障)或齿圈(行星轮故障)的旋转频率；已知行星齿轮箱的啮合频率fm和某个齿轮的齿数Zg，则该齿轮相对行星架(太阳轮和齿圈故障)或齿圈(行星轮故障)的旋转频率：
fg＝fm/Zg        (20)；
式(20)中：fg该齿轮相对行星架(太阳轮和齿圈故障)或齿圈(行星轮故障)的旋转频率，Zg为某个齿轮的齿数；
则太阳轮、行星轮和齿圈分布式故障的特征频率分别为：
fs'＝fm/Zs            (21)；
fp'＝fm/Zp            (22)；
fr'＝fm/Zr            (23)；
式(21)-(23)中：fm为啮合频率；fs'、fp'、fr'太阳轮、行星轮和齿圈分布式故障的特征频率；Zs为太阳轮齿数；Zp为行星轮的齿数；Zr为齿圈齿数；fm为啮合频率；
以与主轴相连接的一级行星轮系的行星架为参考转速，对行星齿轮箱中各级的各个齿轮的局部故障和分布故障的特征阶比进行计算；
这里的信号为matlab仿真信号，信号中忽略齿轮箱中轴承以及各齿轮之间的影响，假设各齿轮箱传动级之间的振动影响不存在，对该方法进行仿真验证针对行星轮系处于正常、分布故障、局部故障时的振动信号模型参考文献，这里就不再赘述；
由此，可以仿真行星齿轮箱各级处于正常工况时，其振动信号模型为:
行星齿轮箱一级行星轮系太阳轮发生局部故障时，其振动信号模型为:
行星齿轮箱一级行星轮系太阳轮发生分布故障时，其振动信号模型为:
式(24)～(26)中：x1(t)、x2(t)、x3(t)为行星齿轮箱处于正常一级行星轮系太阳轮发生局部故障、分布故障时的振动信号序列；t为时间序列；θ1、θ2、θ3、φ、 为初始相位；fm1、fm2、fm3为各级啮合频率；fc1、fc2为一、二级行星架的旋转频率； 为一级太阳轮的绝对旋转频率；fs1、fs1'为一级太阳轮发生局部故障和分布故障时的特征频率；A、B、C为无量纲常数，行星齿轮箱各个状态时值会有所不同，这里就不再详述：各个振动信号采用频率为8192HZ。
6.根据权利要求1所述的基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：步骤(5)所述的故障模式识别的方法：
BP神经网络广泛的应用于齿轮的故障诊断中，具有自学习和数处理能力，但在输入信息不精确或不确定时，神经网络的精确性大大降低，针对这一问题，这里采用遗传算法优化的BP神经网络齿轮故障诊断方法；在BP神经网络的拓扑结构确定的情况下，用遗传算法训练BP神经网络的权值和阈值，经过若干代的交叉、变异后得到稳定的权值和阈值，再将它们赋值给BP神经网络，加快网络收敛速度，克服易陷入局部极小的问题，改善故障诊断的精度和速度；
遗传算法，从数学上看是一种概率性搜索算法；从工程学角度上看，是一种自适应的迭代寻优过程；遗传算法从某一随机产生的或是特定的初始群体出发(作为父本)，按照一定的操作规则，如选择、复制、交叉、变异等不断的迭代计算，并根据每一个个体的适应度值，保留优良品种，淘汰次品，引导搜索过程向最优解逼近；遗传算法的主要特点是直接对结构对象进行操作，不存在求导和函数连续性的限定；具有内在的隐并行性和和较好的全局搜索能力；采用概率化的寻优方法，能自动获取搜索过程中的有关知识并用于指导优化，自适应的调整搜索方向，不需要确定的规则；
传统遗传算法的实现步骤如下：
(1)随机产生一定数目的初始个体(染色体)；些随机产生的染色体组成一个种群，种群中染色体数目称为种群的规模或大小；
(2)用评价函数来评价每个染色体的优劣；染色体对环境的适应程度(称为适应度)，用作以后遗传操作的依据；
(3)基于适应值的选择策略；从当前种群中选取一定的染色体作为新一代的染色体，染色体适应度越高，其被选择的机会越大；
(4)对于新生成的种群进行交叉、变异操作；变异操作的目的是使种群中的个体具有多样性，防止陷入局部最优解，这样产生的染色体群称为后代；
(5)对新生的种群重复进行选择、交叉、变异操作过程，经过给定次数的迭代后，把最好的染色体作为优化问题的最优解；
遗传算法用于神经网络的训练，利用遗传算法全局搜索的特性，缩小搜索范围，得到一个初始的权值矩阵和初始的阈值向量，将初步得到的权值矩阵和阈值向量赋给尚未开始训练的BP网络，再利用BP算法进行精确求解；这种GA和BP网络相结合的方法，能显著地提高BP神经网络的性能，可以达到全局寻找和快速高效的目的；
以3层BP网络为例，用遗传算法学习BP网络的过程如下：
(1)初始化种群P，包括交叉规模、交叉概率Pc，突变概率己以及对任一wih、who初始化；在编码中，采用实数进行编码，初始种群取n(视实际应用选择数值大小)；
(2)计算每一个个体评价函数，并将其排序，按下式概率值选择网络个体：
式(27)中，f为个体的适配值，N表示个体评价函数的个数，fi为个体i的适配值，可用误差平方和E来衡量，即：
式(29)中：i—染色体数，i＝1,2,…,n；
o—输出层结点数，D＝1，2，…，g；
k—学习样本数，k＝1，2，…，m；
yo—网络的实际输出；
d—期望输出；
m—网络的实际输出的个数；
q—期望输出；
(3)交叉概率pc对个体Gi和Gi+1，进行交叉操作，产生新个体G′i和G′i+1，没有进行交叉操作的个体直接进行复制；
(4)利用变异概率pm突变个体Gj产生新的个体G'j；
(5)将产生新个体插入到种群p中，并计算新个体的评价函数；
(6)判断算法是否结束；如果找到了满意的个体，则结束，否则转3进入下一轮运算；
算法结束，如达到预先设定的性能指标后，将最终最优个体解码即可得到优化后的网络连接权值系数，阈值的优化过程类似。
7.根据权利要求1所述的基于多传感器信息融合的变工况行星齿轮箱故障诊断方法，其特征在于：步骤(6)所述的进行验证，包括：
按结构把行星齿轮箱分为三级，按照齿轮故障分为局部故障和分布式故障划分，总结归为15类故障；对这四个测点的振动信号进行阶比重采样以及EMD分解重构后的平稳角域信号进行无量纲指标的计算；
故障诊断模型的遗传优化的BP神经网络由输入层、隐含层及输出层组成，第一层是输入层，输入层节点由齿轮的故障特征参数组成，分别为峰值指标、峭度指标、脉冲指标、裕度指标、波形指标、重复性因子、相似性因子、跳跃性因子；第二层是隐含层，隐含层节点用于提取信号中的相关特征量；第三层为输出层，输出层节点对应着齿轮箱的故障类型；
选择三层网络进行故障识别，即包含输入层、隐含层、输出层；网络的每一个输入节点对应样本一个特征，选择四个测点的径向振动的时域指标中的重复性因子、相似性因子、跳跃性因子、峭度指标、裕度指标、波形指标、脉冲指标、峰值指标，共32个特征值作为网络的输入参数；网络的输出层节点数等于类别数，一个输出节点对应一个类别，输出层节点数为
16，分别对应故障类型为正常、一级太阳轮局部故障、一级行星轮局部故障、一级齿圈局部故障、一级太阳轮分布式故障、一级行星轮分布式故障、一级齿圈分布式故障、二级太阳轮局部故障、二级行星轮局部故障、二级齿圈局部故障、二级太阳轮分布式故障、二级行星轮分布式故障、二级齿圈分布式故障、高速级小齿轮局部故障、高速级大齿轮局部故障、高速级分布式故障；对于隐含层数的确定参考经验公式 nH、nI、l分别为隐含
层、输入层和输出层神经元数目，为1-10之间的整数；采用“1-0”编码模式对网络的输出类型进行编码；在识别阶段，当一个未知类别的样本作用到输入端时，考察各输出节点的输出；
matlab神经网络工具箱提供了多种神经网络改进算法，BP网络选择‘trainrp’函数作为遗传网络的训练函数；传递函数是BP网络的重要组成部分，BP网络经常采用‘tansig’、‘logsig’及‘pureline’函数；输出层的传递函数一般选择‘pureline’函数，选择‘tansig’函数作为隐含层传递函数；另外取net.trainParam.epochs＝10000，LP.1r＝0.05，net.trainParam.goal＝0.001；种群规模popu＝50，遗传代数gen＝100；其它参数选择默认；根据经过BP神经网络训练后进行测试的结果，看出，有两组测试数据(14、16)难以识别，准确率为87.5％，测试集的均方误差为1.0251，训练集的均方误差为0.772376；
与BP网络相同的训练函数、传递函数及相同的参数设计建立的遗传优化BP网络，以提取四个传感器测点的振动信号无量纲指标为遗传优化BP网络的输入量，进行故障模式的识别，能够看出，行星齿轮箱全部故障模式都可识别出来，准确率为100％，测试集的均方误差为0.16624，训练集的均方误差为0.27746；
将在无量纲指标作为网络输入量在不同神经网络模型下的测试结果列出，进行分析比较BP网络的测试结果与遗传BP网络的测试结果，发现遗传BP网络的测试准确率及网络测试速度明显优于BP网络，并且具有更高的网络稳定性；
以单一传感器测点振动信号的特征参数为输入量的测试结果与多传感器测点振动信号的特征参数为输入量的的测试结果进行对比，能够得出，不同传感器测点包含的故障信息不够完全，很难以单一传感器测点的故障信息进行行星齿轮箱的故障诊断，多传感器信息融合的方法较单一传感器的方法的诊断效果更好，准确率更高。