基于最小二乘映射的机械故障特征参量优化方法
【技术领域】
[0001 ] 本发明涉及一种基于最小二乘映射的机械故障特征参量优化方法，涉及机械振动 信号处理领域和在线监控及故障诊断领域。
【背景技术】
[0002] 机械故障特征参量优化在信号处理和故障诊断领域具有重要的意义。大都数机械 设备在实际运作状态下，其故障信号是非平稳、非线性的，同时早期故障信号中干扰信号过 多，故障信号能量微弱，故障特征的提取十分困难，并且提取的故障特征参量维数较高，冗 余信息会相互削弱其反映故障特性的能力。在故障诊断领域中故障特征参量直接关系到故 障诊断的准确率和故障早期预报的可靠性，是设备状态监测和故障诊断的关键。选择对设 备状态敏感度高的特征参量，能增强监测诊断的针对性，提高诊断的准确性。
[0003] 现有的技术方法中，基于时域特征参量提取算法计算相对简单，容易实现，但是提 取到的特征参量容易受到噪声的影响，适合于信噪比比较高的环境，其应用范围受到限制。 目前广泛采用的信号特征提取方法有：包络解调、共振解调、形态滤波、小波、经验模态法、 稀疏分解、主分量分析法等。相对于简单时域特征参量提取算法提取的特征参量具有更好 的抗噪性，应用较为广泛，但是计算复杂。
【发明内容】
[0004] 针对上述现有技术存在的不足，本发明提供一种基于最小二乘映射的机械故障特 征参量优化方法。该方法计算简单，通过对提取的特征参量进行优化，提高故障特征参量的 敏感性，实现非平稳早期微弱故障的识别。
[0005] 本发明的技术方案是：
[0006] 基于最小二乘映射的机械故障特征参量优化方法，包括以下步骤：
[0007] 1、根据所要诊断的机械设备计算提取特征参量P，并在特征空间构造矢量Pi_j。
[0008] Pij= {p i.n，pij2，…，PijkI i = 1 ~M，j = 1 ~Ν}τ
[0009] 其中，i表示故障类型的数目为M，j表示特征参量的数目为N，k为特征参量的维 度（一般取2)。矢量端点表不故障类型。
[0010] 2、构造矢量Plj由特征空间映射到判决空间的误差向量ε ^和总均方误差ε。
[0013] 其中，Llj表示矢量P u在判决空间的映像，V 1表示矢量P U投射在判决空间的点。
[0014] 3、计算将矢量Plj由特征空间映射到判决空间的变换Α。
[0015] a.矢量Pu由特征空间映射到判决空间的矢量L AP u。
[0016] b.根据总的均方误差最小原则计算变换矩阵A。
[0017] 由误差向量ε u和总均方误差ε以及矢量L u公式可以推导出：
[0019] 因为选取的变换A使得ε最小，所以A可以通过下面公式求得，
[0021] W表示梯度算子。
[0022] 于是得：
[0024] 又因为：
[0028] 将其带入上述恒等式得：
[0030] 所以要寻找的变换A为：
[0032] 当故障类型的数目M彡2时，则变换A为：
[0034] 4、构造在判决空间内的新特征参量NP。
[0035] NP = AXP
[0036] 本发明的有益效果是：本发明提供一种基于最小二乘映射的机械故障特征参量优 化方法，能够增强特征参量对故障的敏感性，更好地反映故障特性，提高了故障诊断的准确 率和可靠性。
【具体实施方式】
[0037] 以滚动轴承的故障特征参量对内轮伤和滚动体故障这两种故障状态的识别为例 对本发明所提出的技术路线作进一步的解释和说明。此处应用的实例并不只适用于本例， 适用于所有旋转机械。具体的实施过程如下：
[0038] 1、根据滚动轴承内轮伤和滚动体故障这两种故障状态选择以下两个特征参量来 识别：
[0040] 其中^是振动信号的离散数据，无是^的均值，
σ是X1的标准方差，
[0041] 在特征空间构造矢量Pl j
[0042] Pij= {ρ 叩，pij2，…，pijk| i = 1 ~2, j = 1 ~2}τ
[0043] 2、构造矢量Plj由特征空间映射到判决空间的误差向量ε ^和总均方误差ε。
[0046] 其中，Llj表示矢量P u在判决空间的映像，V 1表示矢量P U投射在判决空间的点。
[0047] 3、计算将矢量Plj由特征空间映射到判决空间的变换Α。
[0048] a.矢量Pu由特征空间映射到判决空间的矢量L AP u。
[0049] b.根据总的均方误差最小原则计算变换矩阵A。
[0050] 由误差向量ε u和总均方误差ε以及矢量L u公式可以推导出：
[0052] 因为选取的变换A使得ε最小，所以A可以通过下面公式求得，
[0054] A表示梯度算子。
[0055] 于是得：
[0057] 又因为：
[0061] 将其带入上述恒等式得：
[0063] 所以要寻找的变换A为：
[0065] 4、构造在判决空间内的新特征参量NP。
[0066] NP = AXP
[0067] 辨别任意特征参数P对两种状态的优点，可以用DI (敏感性评价因子）辨别指标 来评价。
[0069] 其中μ η μ 2各自是位于特征参数P的状态1和状态2下的平均值、σ JP σ 2是 其标准偏差。DI越大、使用特征P辨别状态1和状态2时的敏感度就越高。
[0070] 使用DI (敏感性评价因子）来评价优化过的特征参量和原来的特征参量对滚动轴 承内轮伤和滚动体故障反映的敏感性。
[0071] 表1特征参量PjP P 2区别内轮伤和滚动体故障的DI值
[0073] 表2优化的特征参量NPjP NP 2区别内轮伤和滚动体故障的DI值
[0075] 通过比较表1和表2,经过本发明所提出的方法优化过的特征参量比原来的特征 参量对故障的反映敏感性更高。所以通过该方法优化的特征参量能大大提高对故障诊断的 准确性。
【主权项】
1. 一种基于最小二乘映射的机械故障特征参量优化方法，其特征在于：包括以下步 骤： 一、 根据所要诊断的机械设备计算提取特征参量P，并在特征空间构造矢量Pi_j。 P1：i={pin，plj2，…，pljk|i= 1 ~M，j= 1 ~Ν}τ 其中，i表示故障类型的数目为M，j表示特征参量的数目为N，k为特征参量的维度（一 般取2)。矢量P、,的端点表示故障类型。 二、 构造矢量由特征空间映射到判决空间的误差向量ε^和总均方误差ε。 ευ=其中，表示矢量Pu在判决空间的映像，V廣示矢量Pu投射在判决空间的点。 三、 计算将矢量由特征空间映射到判决空间的变换A。 a. 矢量Pu由特征空间映射到判决空间的矢量LAPu。 b. 根据总的均方误差最小原则计算变换矩阵A。 由误差向量％和总均方误差ε以及矢量"公式可以推导出：因为选取的变换A使得ε最小，所以A可以通过下面公式求得， = 〇 ▽A表示梯度算子。 于是得： 又因为：将其带入上述恒等式得：所以要寻找的变换A为：当故障类型的数目Μ彡2时，则变换A为：四、构造在判决空间内的新特征参量NP。NP = AXP〇
【专利摘要】本发明涉及一种基于最小二乘映射的机械故障特征参量优化方法，将用于机械设备在线监测和故障诊断领域。具体包括以下步骤：根据所要诊断的机械设备计算提取特征参量，并在特征空间构造特征矢量，再构造特征矢量由特征空间映射到判决空间的误差向量和总均方误差，根据总均方误差最小原则选择变换矩阵，最后用变换矩阵构造出新的特征参量。
【IPC分类】G06F17/50
【公开号】CN105260582
【申请号】CN201510835989
【发明人】赵宇, 李可, 陈鹏, 王华庆
【申请人】江南大学
【公开日】2016年1月20日
【申请日】2015年11月26日