基于FPGA的改进局部线性嵌入风机轴承故障诊断方法

505 编辑：管理员来源：傲拓科技股份有限公司

2024-03-12 17:24:14

权利要求书： 1.一种基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，其特征在于，所述方法包括以下步骤：S1、采集风机运行状态下的齿轮箱轴承信号数据；

S2、对信号数据进行预处理，将高维信号数据向量化，构建风机齿轮箱轴承的原始数据空间，利用余弦相似度计算原始空间中每个样本点的前p个近邻点；

S3、计算每个样本点周围的前p个近邻点的距离之和并排序，根据排序结果对每个样本点的邻域分配近邻点；

S4、挖掘局部邻域内的局部线性结构和近邻结构，对局部线性嵌入的权重进行重构；

S5、根据重构后的权重，构建最小误差函数，计算低维映射矩阵；

S6、将得到的低维映射集合输入至支持向量机中，根据支持向量机的分类结果对风机齿轮箱进行故障诊断，得到并输出故障诊断结果；

其中，步骤S2进一步包括：

S21、将高维数据向量化，得到原始数据空间，其中，N为样本点数，D为高维数据维数；

S22、利用余弦相似度计算每个样本点周围的前p个近邻点，余弦相似度的计算公式为；

其中， , 表示的是样本的第t维, 的值越大，与之间的相似度越高，表示样本和之间的余弦相似度，q表示样本个数。

2.根据权利要求1所述的基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，其特征在于，步骤S3进一步包括：S31、计算每个样本点周围的前p个近邻点的距离之和，计算公式为；

其中表示第i个样本点与其前p个近邻点的距离之和；

S32、对计算的近邻点之和进行排序，给定近邻点k的范围[kmin,kmax]；

S33、通过公式为每个样本点选择邻域内近邻点的个数，表示第1个样本点与其前p个近邻点的距离之和，表示最后一个样本点与其前p个近邻点的距离之和。

3.根据权利要求2所述的基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，其特征在于，步骤S4的具体过程包括：S41、计算样本点邻域内局部线性结构权重，计算公式为 ;其中表示和之间的重构权值，如果是的邻居，则可由计算，反之 , 表示的权重系数向量（表示和之间的重构权值）；

S42、计算样本点邻域内的近邻结构权重，计算公式为;

S43、对局部线性结构权重和近邻结构权重进行融合，公式为，其中表示调节参数，；

其中，为重构权重，为原始权重，为包含近邻结构权重。

4.根据权利要求3所述的基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，其特征在于，步骤S5的具体过程包括：S51、计算低维映射Y，保持重构权重矩阵不变，最小化低维重构误差；

将该式进一步写成：

;

式中，且约束函数矩阵化为：，构建Lagrange函数：；

对Y求导并令该式为0，计算得到低维映射Y；

其中，L表示构建的拉格朗日函数，为重构权重，为原始权重，为包含近邻结构权重，表示的是样本xi和xj之间的原始权重，表示的是样本xi和xj之间的重构权重，表示的是样本xi和xj之间的近邻结构权重，I是单位矩阵，M，m是求解拉格朗日函数中的协方差矩阵，tr表示的是矩阵的迹，是构建拉格朗日函数必要的系数。

5.根据权利要求4所述的基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，其特征在于，步骤S6进一步包括：将得到的低维映射集合选择50%样本作为训练集，50%样本作为测试集，根据支持向量机的分类结果对风机齿轮箱进行故障诊断，得到故障诊断结果。

说明书：基于FPGA的改进局部线性嵌入风机轴承故障诊断方法技术领域[0001] 本发明涉及风机齿轮箱故障诊断技术领域，尤其涉及一种基于FPGA的改进局部线性嵌入风机轴承故障诊断方法。背景技术[0002] 在风力发电的风机主控系统中，风机齿轮箱的运行可靠性至关重要，因此常通过设置FPGA模块和传感器来分析风机齿轮箱运行数据从而得到故障诊断结果，再反馈给CPU，振动分析法和机器学习方法是目前风机轴承故障诊断的主要方法，两种方法都是通过传感器测取风机轴承的振动信号，然后对轴承振动信号分析得出风机轴承的故障类型，而实现风机轴承故障诊断。[0003] 振动分析法和机器学习法都对故障特征的选择和提取有着较高的要求。但振动分析法在实际使用的过程中过于依赖人工选取特征，有一定的主观性。机器学习法在处理实时产生的大量高维风机信号时容易造成信号滞后，不能及时发现风机运行的故障，实时性不佳。同时由于风机轴承的实际工况复杂多变，但两种方法在模式识别过程中，通常通过经验确定算法参数，且各个算法的参数固定不变，难以适应实时变换的复杂工作环境，难以保证对轴承故障特征的识别准确可靠。发明内容[0004] 有鉴于此，本发明的目的在于提出一种基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，以解决目前风机齿轮箱轴承故障诊断提取故障特征不准确、实时性不强的问题。[0005] 基于上述目的，本发明提供了一种基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，包括以下步骤：S1、采集风机运行状态下的齿轮箱轴承信号数据；

S2、对信号数据进行预处理，将高维信号数据向量化，构建风机齿轮箱轴承的原始数据空间，利用余弦相似度计算原始空间中每个样本点的前p个近邻点；

S3、计算每个样本点周围的前p个近邻点的距离之和并排序，根据排序结果对每个样本点的邻域分配近邻点；

S4、挖掘局部邻域内的局部线性结构和近邻结构，对局部线性嵌入的权重进行重

构；

S5、根据重构后的权重，构建最小误差函数，计算低维映射矩阵；

S6、将得到的低维映射集合输入至支持向量机中，根据支持向量机的分类结果对

风机齿轮箱进行故障诊断，得到并输出故障诊断结果；

其中，步骤S2进一步包括：

S21、将高维数据向量化，得到原始数据空间，其中，N为样

本点数，D为高维数据维数；

S22、利用余弦相似度计算每个样本点周围的前p个近邻点，余弦相似度的计算公

式为

；

其中，，表示的是样本的第t维，的值越大，与之间的相

似度越高，表示样本和之间的余弦相似度，q表示样本个数。

[0006] 优选地，步骤S3进一步包括：S31、计算每个样本点周围的前p个近邻点的距离之和，计算公式为

；

其中表示第i个样本点与其前p个近邻点的距离之和；

S32、对计算的近邻点之和进行排序，给定近邻点k的范围[kmin,kmax]；

S33、通过公式为每个样本点选择

邻域内近邻点的个数，表示第1个样本点与其前p个近邻点的距离之和，表示最后一个样本点与其前p个近邻点的距离之和。

[0007] 优选地，步骤S4的具体过程包括：S41、计算样本点邻域内局部线性结构权重，计算公式为

；

其中表示和之间的重构权值，如果是的邻居，则可由

计算，反之 , 表示的权重系数向量（表示和之间

的重构权值）；

S42、计算样本点邻域内的近邻结构权重，计算公式为

；

S43、对局部线性结构权重和近邻结构权重进行融合，公式为

，其中表示调节参数，；

其中，为重构权重，为原始权重，为包含近邻结构权重。

[0008] 优选地，步骤S5的具体过程包括：S51、计算低维映射Y，保持重构权重矩阵不变，最小化低维重构误差

;

将该式进一步写成：

；

式中，且约束函数矩阵化为：，构建Lagrange函

数：；

对Y求导并令该式为0，计算得到低维映射Y；

其中，L表示构建的拉格朗日函数，为重构权重，为原始权重，为包含近

邻结构权重，WLij表示的是样本xi和xj之间的原始权重，WEij表示的是样本xi和xj之间的重构权重，WNij表示的是样本xi和xj之间的近邻结构权重，I是单位矩阵，M，m是求解拉格朗日函数中的协方差矩阵，tr表示的是矩阵的迹，是构建拉格朗日函数必要的系数。

[0009] 优选地，步骤S6进一步包括：将得到的低维映射集合选择50%样本作为训练集，50%样本作为测试集，根据支持

向量机的分类结果对风机齿轮箱进行故障诊断，得到故障诊断结果。

[0010] 本发明的有益效果：1、本发明对原始高维数据的实时处理速度有一定地提高，有效地增强了风机齿轮箱轴承故障诊断的实时性。

[0011] 2、本发明通过对局部线性嵌入算法的改进，能够有效挖掘原始高维数据的局部结构，进一步提取故障数据的本质特征，提高了风机齿轮箱轴承故障诊断的准确度。附图说明[0012] 为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。[0013] 图1为本发明实施例的基于FPGA的改进局部线性嵌入风机轴承故障诊断方法流程示意图；图2为本发明实施例的支持向量机得到的不同种类风机齿轮箱轴承故障的分类结

果示意图；

图3为使用现有的LLE算法得到的不同种类风机齿轮箱轴承故障的分类结果示意

图。

具体实施方式[0014] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，对本发明进一步详细说明。[0015] 需要说明的是，除非另外定义，本发明使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也可能相应地改变。

[0016] 如图1所示，本说明书实施例提供一种基于FPGA的改进局部线性嵌入风机轴承故障诊断方法，包括以下步骤：S1、采集风机运行状态下的齿轮箱轴承信号数据，如传感器采集齿轮箱轴承信号

数据传输给FPGA模块；

S2、对信号数据进行预处理，将高维信号数据向量化，构建风机齿轮箱轴承的原始数据空间，利用余弦相似度计算原始空间中每个样本点的前p个近邻点，具体来说，包括以下子步骤：

S21、将高维数据向量化，得到原始数据空间，其中，N为样

本点数，D为高维数据维数。

[0017] S22、利用余弦相似度计算每个样本点周围的前p个近邻点，余弦相似度的计算公式为；

其中，，表示的是样本的第t维，的值越大，与之间的相似

度越高。通过使用余弦相似度来度量每个样本之间的距离，相较于传统欧式距离来说，能够更好地反应样本之间的关系。

[0018] S3、计算每个样本点周围的前p个近邻点的距离之和并排序，根据排序结果对每个样本点的邻域分配近邻点，举例来说，进一步包括以下子步骤：S31、计算每个样本点周围的前p个近邻点的距离之和，计算公式为

;

其中表示第i个样本点与其前p个近邻点的距离之和；

S32、对计算的近邻点之和进行排序，给定近邻点k的范围[kmin,kmax]；

S33、通过公式为每个样本点选择邻域

内近邻点的个数。

[0019] 传统LLE算法人工指定近邻点个数，无法保证算法的特征提取能力，本文采用根据样本分布的稀疏稠密程度自适应地选取样本近邻点，能够让算法更好地自动地选择合适的近邻点。[0020] S4、挖掘局部邻域内的局部线性结构和近邻结构，对局部线性嵌入的权重进行重构，举例来说，步骤S4的具体过程包括：S41、计算样本点邻域内局部线性结构权重，计算公式为

；

其中表示和之间的重构权值，如果是的邻居，则可由式（1）计

算，反之，。表示的权重系数向量；

S42、计算样本点邻域内的近A邻结构权重，计算公式为

；

S43、对局部线性结构权重和近邻结构权重进行融合，公式为

；

其中，为重构权重，为原始权重，为包含近邻结构权重。

[0021] 通过在局部线性结构的基础之上加入邻域内近邻序列结构，能够更好地保持原始流形结构。[0022] S5、根据重构后的权重，构建最小误差函数，计算低维映射矩阵，举例来说，步骤S5的具体过程包括：S51、计算低维映射Y，保持重构权重矩阵不变，最小化低维重构误差

；

将该式进一步写成：

；

式中，且约束函数矩阵化为：，构建Lagrange函

数：

；

对Y求导并令该式为0，计算得到低维映射Y。

[0023] S6、将得到的低维映射集合输入至支持向量机中，根据支持向量机的分类结果对风机齿轮箱进行故障诊断，得到故障诊断结果并进行输出，如将得到的低维映射集合选择50%样本作为训练集，50%样本作为测试集，通过支持向量机进行故障识别。

[0024] 本发明提出的基于样本点分布密集程度的自适应邻域方法能有效选取样本点的邻域，融合局部线性结构权重和邻域内近邻结构权重能够有效挖掘高维数据的本质结构，提高风机齿轮箱轴承故障诊断的准确度，提出改进局部线性嵌入的风机齿轮箱轴承故障诊断方法，解决了提取故障特征不准确、实时性不强的问题。[0025] 本发明实施例中，通过输入我们从东北石油大学数据库中获取的风机齿轮箱数据集，支持向量机得到的不同种类风机齿轮箱轴承故障的分类结果如图2所示，根据分类结果可以得出风机齿轮箱轴承发生的故障类别。[0026] 图3是使用同样的数据集，通过现有的LLE算法得到的风机齿轮箱轴承故障的分类结果，从结果可以看出，本发明提供的方法，分类准确度更高。[0027] 所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本发明的范围（包括权利要求）被限于这些例子；在本发明的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，步骤可以以任意顺序实现，并存在如上所述的本发明的不同方面的许多其它变化，为了简明它们没有在细节中提供。[0028] 本发明旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此，凡在本发明的精神和原则之内，所做的任何省略、修改、等同替换、改进等，均应包含在本发明的保护范围之内。