基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法技术

基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法，包括特征提取、模型参数初始化、参数训练、输出概率计算这四个步骤。(1)特征提取：选取小波函数对振动信号进行3层小波包分解和重构，通过对各频段小波分解系数信号进行分析,实现从各个频段分别提取振动信号所代表的不同故障状态的特征信息。(2)模型参数初始化：将振动信号频带能量作为特征向量进行建模(3)参数训练：根据第二步初始化的参数，用鸟群算法进行重估，(4)输出概率计算：经过上述步骤建好模型后，将监测振动数据特征提取，代入不同故障状态模型，由前向后向算法计算输出概率，概率最大的即为对应的故障类型。

【技术实现步骤摘要】
基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法
本专利技术涉及风机齿轮箱故障诊断领域，尤其涉及到一种基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法。
技术介绍
齿轮传动具有传动平稳、可靠、高效、高寿命、传动比精确、功率范围大等优点，但因结构较复杂，工作环境恶劣等，齿轮较易出现故障，从而导致整个系统出现故障。据有关文献统计，传动机械和旋转机械的故障中分别有80％和10％的故障是齿轮箱故障，随着设备系统的自动化、复杂化、大型化。齿轮的失效给整个工业生产和社会生活造成越来越大的损失，如对齿轮箱进行状态检测及故障诊断，可以从根本上改变齿轮箱事后维修和定期维修的现状，实现视情维修，从而降低事故发生率，减少人员伤亡及不必要的经济损失，创造更多的经济效益和社会效益。目前对齿轮箱故障诊断的方法有很多，如温度测量、油样分析、振动分析、声学测量分析、专家系统等，现有技术缺点：1.温度测量法通过监测箱体处的温度来判断齿轮箱是否工作正常，对于润滑不良而引起的过热现象较敏感，但该方法不适用于齿轮点烛、断齿、裂纹等局部故障；2.油样分析是通过从齿轮箱润滑油中提取出油样，通过收集和分析油样中金属颗粒的大小和形状来判断齿轮工况和故障的方法，但这种方法易受其它非齿轮损坏下的颗粒的影响；3.专家系统需要大量样本进行训练，而且依赖与相关领域的专家经验，对新的故障模式的处理效果不佳。综上，现有诊断方法均有一定的不足，因此，有必要提供一种新的诊断方法。
技术实现思路
本专利技术的目的是提供一种基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法，该方法通过振动分析以及改进的高斯混合隐马尔科夫模型，实现对齿轮箱的故障进行仿真分析，同时为解决隐马尔科夫模型参数训练易陷入局部收敛的问题。本专利技术为解决其技术问题所采用的技术方案是：基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法，包括特征提取、模型参数初始化、参数训练、输出概率计算这四个步骤；(1)特征提取：选取小波函数对振动信号进行3层小波包分解和重构，通过对各频段小波分解系数信号进行分析,实现从各个频段分别提取振动信号所代表的不同故障状态的特征信息。针对提取的每一频段的重构信号求其能量：其中，xjk(j＝0,1...,7,k＝1,2,...,n)表示为重构信号s3j的n个离散点的幅值，为了提高特征向量的聚类性,使其分布离散性更好，便于模型输入和辨识，需要对特征向量进行归一化处理，得到归一化后的特征向量C3j＝[c30,c31,...,c37]，归一化的方法是对根据不同频段的能量特征值序列，按下式进行归一化处理：(2)模型参数初始化：将振动信号频带能量作为特征向量进行建模，则模型中的观测值O＝{o(1),...,o(h),...,o(H)}＝{E(1)，...,E(h),...,E(H)}，其中表示振动信号第h观测值的频带能量，T表示频带数；(3)参数训练：根据第二步初始化的参数，用鸟群算法进行重估，步骤如下：①初始化BSA种群及相关参数，初始化个体X，X＝cell(pi，A，c，u，U，)；②进行迭代，每一次BSA迭代中，进行隐马尔科夫传统训练算法k次迭代，目标函数计算；③当满足结束条件，即输出概率P(O|λ)的增量小于一定阈值，则结束迭代，输出结果；(4)输出概率计算：经过上述步骤建好模型后，将监测振动数据特征提取，代入不同故障状态模型，由前向后向算法计算输出概率，概率最大的即为对应的故障类型；所述步骤(2)中建立的模型为连续隐马尔科夫模型；连续隐马尔可夫模型参数如下：(1)状态个数：设N个状态分别为S1,S2,...,SN,则有t(1≤t≤T)处所处的状态qt∈(S1,S2,...,SN)，N表示所有频带能量的值分成的类数；(2)初始状态分布矢量：π＝{πi}，表示1处状态处于某一状态的概率，即：E1取值为某一类的概率，其中1≤i≤N；(3)状态转移矩阵：A＝{aij},元素为状态之间转移的概率，即：一类E值在下一个频带中变成另一类E值的概率；(4)高斯混合参数和输出概率矩阵:混合高斯元个数M＝{Mj}，混合高斯元权系数c＝{cjl}，均值向量u＝{ujl}，协方差矩阵U＝{Ujl}，观测值矩阵B＝{bj(o)}，其中1≤j≤N，l表示状态Sj的第l个高斯分量，设观测向量为D维，则概率密度函数表示为：初始化时：π、A均匀取值，即：π＝{1/N}N×1；对于N、M、c、u、U，由AP聚类算法确定：N和M的取值通过AP聚类算法确定的过程为：首先，将AP聚类算法所求得的所有数据点的聚类个数确定为N；然后，对得到的N个聚类的每一类继续使用AP算法的第二次聚类结果即可作为向量M的对应数值。C、u、U的取值：其中，ujl＝属于的观测向量的均值向量；Ujl＝属于的观测向量的协方差矩阵。本专利技术的优点在于：针对传统隐马尔科夫模型参数训练易陷入局部收敛的问题，本专利技术提出鸟群连续隐马尔科夫(BSA-CHMM)模型，将鸟群算法(BSA)与隐马尔科夫传统训练算法(Baum-Welch)算法结合，每一次BSA全局搜索后，进行BW算法K次迭代局部搜索，这样能有效解决局部收敛问题。与现有诊断方法相比，所需训练样本更少，诊断速度较快。附图说明图1是本专利技术提出的基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法的结构示意图；图2是该方法中BSA-BW算法流程图；图3是该方法中BSA-BW算法伪码；图4是在一实施例中各个状态特征能量分布图；图5是在一实施例中的正常状态训练曲线；图6是在一实施例中的点蚀状态训练曲线；图7是在一实施例中的磨损状态训练曲线；图8是在一实施例中的断齿状态训练曲线；图9、图10是在一实施例中各个状态检验样本数据图；具体实施方式为了使本专利技术实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合图示与具体实施例，进一步阐述本专利技术。如图1至图10所示，本专利技术提出的基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法包括特征提取、模型参数初始化、参数训练、输出概率计算这四个步骤；(1)特征提取：选取小波函数对振动信号进行3层小波包分解和重构，通过对各频段小波分解系数信号进行分析,实现从各个频段分别提取振动信号所代表的不同故障状态的特征信息。针对提取的每一频段的重构信号求其能量：其中，xjk(j＝0,1...,7,k＝1,2,...,n)表示为重构信号s3j的n个离散点的幅值，为了提高特征向量的聚类性,使其分布离散性更好，便于模型输入和辨识，需要对特征向量进行归一化处理，得到归一化后的特征向量C3j＝[c30,c31,...,c37]，归一化的方法是对根据不同频段的能量特征值序列E3j＝[e30,e31,...,e37]，按下式进行归一化处理：(2)模型参数初始化：将振动信号频带能量作为特征向量进行建模，则模型中的观测值O＝{o(1),...,o(h),...,o(H)}＝{E(1),...,E(h),...,E(H)}，其中表示振动信号第h观测值的频带能量，T表示频带数；(3)参数训练：根据第二步初始化的参数，用鸟群算法进行重估，步骤如下：①初始化BSA种群及相关参数，初始化个体X，X＝cell(pi，A，c，u，U，)；②进行迭代，每一次BSA迭代中，进行隐马尔科夫传统训练算法k次迭代，目标函数计算；③当满足本文档来自技高网...

【技术保护点】
基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法，包括特征提取、模型参数初始化、参数训练、输出概率计算这四个步骤，其特征在于：(1)特征提取：选取小波函数对振动信号进行3层小波包分解和重构，通过对各频段小波分解系数信号进行分析,实现从各个频段分别提取振动信号所代表的不同故障状态的特征信息。针对提取的每一频段的重构信号求其能量：

【技术特征摘要】
1.基于鸟群算法和隐马尔科夫模型的齿轮箱故障诊断方法，包括特征提取、模型参数初始化、参数训练、输出概率计算这四个步骤，其特征在于：(1)特征提取：选取小波函数对振动信号进行3层小波包分解和重构，通过对各频段小波分解系数信号进行分析,实现从各个频段分别提取振动信号所代表的不同故障状态的特征信息。针对提取的每一频段的重构信号求其能量：其中，xjk(j＝0,1...,7,k＝1,2,...,n)表示为重构信号s3j的n个离散点的幅值，为了提高特征向量的聚类性,使其分布离散性更好，便于模型输入和辨识，需要对特征向量进行归一化处理，得到归一化后的特征向量C3j＝[c30,c31,...,c37]，归一化的方法是对根据不同频段的能量特征值序列E3j＝[e30,e31,…,e37]，按下式进行归一化处理：(2)模型参数初始化：将振动信号频带能量作为特征向量进行建模，则模型中的观测值O＝{o(1),...,o(h),...,o(H)}＝{E(1),...,E(h),...,E(H)}，其中表示振动信号第h观测值的频带能量，T表示频带数；(3)参数训练：根据第二步初始化的参数，用鸟群算法进行重估，步骤如下：①初始化BSA种群及相关参数，初始化个体X，X＝cell(pi，A，c，u，U，)；②进行迭代，每一次BSA迭代中，进行隐马尔科夫传统训练算法k次迭代，目标函数计算；③当满足结束条件，即输出概率P(O|λ)的增量小于一定阈值，则结束迭代，输出结果；(4)输出概率计算：经过上述步骤建好模型后，将监测振动数据特征提取，代入不同故障状态模型，由前向后向算法计算输出概率，概率最大...

【专利技术属性】
技术研发人员：丁超然，刘三明，王致杰，王帅，潘昭旭，
申请(专利权)人：上海电机学院，
类型：发明
国别省市：上海,31