1.一种基于优化变分模态分解的滚动轴承故障诊断方法，其特征在于，包括以下步骤：步骤S1，选取原始振动信号4096个采样点作为变分模态分解的输入信号，采用改进的蝙蝠算法以最小平均包络熵作为优化目标对所述变分模态分解的模态数和二次惩罚因子进行优化；

步骤S2，以优化后的模态数和二次惩罚因子对所述原始振动信号分解，并求取分解后的分量的能量熵以及能谱熵；

步骤S3，以峭度、相关系数、边际谱熵作为筛选准则对所述分量进行筛选，求取保留分量的主频分布特征；

步骤S4，将所述能量熵、所述能谱熵、所述主频分布特征作为特征向量输入到支持向量机中，以实现故障诊断。

2.根据权利要求1所述的基于优化变分模态分解的滚动轴承故障诊断方法，特征在于，根据以下公式(1)定义所述平均包络熵：其中，a(j)表示第i个分量的包络线，K表示模态数，N表示采样点数，Pij表示第i个分量的归一化包络线，Eij表示第i个分量的包络熵，Ep表示平均包络熵，i和j为正整数。

3.根据权利要求2所述的基于优化变分模态分解的滚动轴承故障诊断方法，特征在于，采用改进的蝙蝠算法以最小平均包络熵作为优化目标对所述变分模态分解的模态数和二次惩罚因子进行优化，包括：步骤S101，初始化蝙蝠个体的位置、速度、响度A、脉冲率r、响度衰减系数Af、脉冲频度增强系数Rf、种群数量N、最大迭代次数Tmax，其中，所述蝙蝠个体的初始位置采用混沌映射的方法生成，随机生成一个2维的[0,1]上的向量Zi＝(Zi1，Zi2)，根据以下公式(2)迭代生成N-1个向量，Zi+1＝4Zi(1-Zi)   (2)；

步骤S102，根据以下公式(3)-(6)更新所述蝙蝠个体的速度和位置，

fi＝fmin+(fmax-fmin)β (6)；

其中， 以及 分别表示蝙蝠个体i在t和t-1代的飞行速度， 和 分别表示蝙蝠个体i在t和t-1代的位置，x*表示蝙蝠种群的当前最优解，pfitness(i)表示个体的适应度值，best min表示当前全局的最小适应度值，best max表示当前全局的最大适应度值，ω为权重因子，β为一个[0,1]上的随机向量，fmax和fmin分别是最大频率以及最小频率；

步骤S103，判断rand1是否大于rit-1，若满足，则通过以下公式(7)进行局部寻优：其中，rand1为一个[0,1]上的随机数， 表示蝙蝠个体i在t代的位置，x*表示蝙蝠种群的当前最优解，rit为蝙蝠个体i在t代的脉冲率，ε为一个[0,1]上的随机向量，At是t代所有蝙蝠响度的平均值；

步骤S104，如果 且根据式(7)更新后的蝙蝠个体位置求得的适应度值小于蝙蝠个体i在t-1代时的适应度值，那么根据式(8)和式(9)更新所述响度A以及脉冲率r；

其中，rand2表示一个[0,1]上的随机数， 表示蝙蝠个体i在t代的响度，Af表示响度衰减系数，rit表示蝙蝠个体i在t代的脉冲率，ri0表示初始脉冲频率， 表示脉冲频度增强因子，Rf表示脉冲频度增强系数；

返回步骤103，直至达到最大迭代次数，将所述最小平均包络熵所对应的蝙蝠个体位置作为变分模态分解的最优模态数和二次惩罚因子进行优化。

4.根据权利要求3所述基于优化变分模态分解的滚动轴承故障诊断方法，其特征在于，步骤S2，以优化后的模态数和二次惩罚因子对所述原始振动信号分解，并求取分解后的分量的能量熵以及能谱熵，包括：步骤S201，选取10种故障状态的信号的4096个采样点，使用优化后的模态数和二次惩罚因子对10个所述原始振动信号进行变分模态分解，每种故障状态生成K个分量，其中，10种故障状态分别为正常、内圈轻度故障、内圈中度故障、内圈重度故障、滚动体轻度故障、滚动体中度故障、滚动体重度故障、外圈轻度故障、外圈中度故障、外圈重度故障；

步骤S202，根据以下公式(10)和(11)求取所述能量熵和所述能谱熵，式中，Ei表示第i个分量的能量，Pi表示第i个分量的归一化能量，H表示能量熵，ci(t)表示第i个分量，m表示采样点数；

式中，Pi的定义与式(10)中相同，VE表示能谱熵，k表示分量的总个数。

5.根据权利要求4所述基于优化变分模态分解的滚动轴承故障诊断方法，其特征在于，步骤S3，以峭度、相关系数、边际谱熵作为筛选准则对所述分量进行筛选，求取保留分量的主频分布特征，包括：步骤S301，分别根据以下公式(12)、(14)、(14)求取所述分量的峭度、相关系数、边际谱熵的倒数；

式中，Kr表示分量的峭度，xi表示分量信号，N表示分量的采样点数；

式中，ρ表示分量的关系数，X以及Y分别表示分量与原始振动信号，μX与μY分别表示分量与原始信号的均值，σ(X)以及σ(Y)分别表示分量与原始信号的标准差，式中，Hp表示边际谱熵的倒数，Pi表示第i个分量的边际谱的归一化，h(i)表示分量的边际谱；

步骤S302，将所述分量以所述峭度、相关系数、边际谱熵作为筛选准则的值映射到三维直角坐标系中，其中，点的坐标设为：(0，峭度值，0)，(10*边际谱熵的倒数，0，0)，(0，0，10*相关系数)，将三个点相连后得到三角形，选取三角形面积最大的分量；

S303，根据以下公式(15)求取筛选后的分量的主频分布特征；

式中，K表示频谱线数，fk表示第k条频谱线的频率值，s(k)表示频谱。

6.根据权利要求5所述的基于优化变分模态分解的滚动轴承故障诊断方法，其特征在于，步骤S4，将所述能量熵、所述能谱熵、所述主频分布特征作为特征向量输入到支持向量机中，以实现故障诊断，包括：每种故障状态选取80个样本，每个样本有1024个点，获取10种故障状态的所有样本的特征向量，其中，每种故障状态选取60个样本作为训练样本，20个样本则作为测试样本，将所述训练样本输入到所述支持向量机中进行模型的训练，训练完成后，将所述测试样本输入到所述模型中对10种故障状态进行分类识别。