[0083] 下面结合具体附图对本发明做进一步分析。
[0084] 如图1所示,基于阈值判断的GA-BP的磁悬浮列车故障检测方法,步骤如下:
[0085] 步骤一、通过在列车上安装多组加速度、电流、间隙传感器获得所需数据,对原始数据进行筛选和预处理。对40000+数据筛选,选取其中的18000+的数据,针对选取的数据进行漏检测数据和过检测数据使用五点平均法处理数据。再针对处理过的数据进行切片(t时间段)。
[0086] 步骤二、将步骤一处理后的每种数据使用信号处理技术和统计学习方法,提取t时段内的特征参数。其中特征参数包括时域指标和频域指标以及时频特征。
[0087] 2.1时域指标提取
[0088] 时域指标包括方差(Var)、均值 均方值(X2rms)、峰值(Xp),峰值指标(Ip)、偏度(Cp)、峭度(Cq)、波形指标(W)、脉冲指标(Cf)。时域特征根据公式(1)-(9)获得输入向量[0089] 均值公式:
[0090]
[0091] 方差公式:
[0092]
[0093] 均方根值公式:
[0094]
[0095] 峰值:
[0096] Xp=max{x1,x2,x3...xn} 公式(4)
[0097] 峰值指标:
[0098]
[0099] 偏度指标:
[0100]
[0101] 峭度指标:
[0102]
[0103] 波形指标:
[0104]
[0105] 脉冲指标:
[0106]
[0107] n为t时段内数据个数,xi表示t时段内数据中的第i个,|xi|表示第i个数据的模值,其中i=1,2,3……n。
[0108] 2.2频域指标提取
[0109] 先将步骤一处理后的数据进行快速傅里叶变换得到不同时刻的频谱,然后根据公式(10)-(14)得到t时段内频域指标。频域指标包括重心频率(FC)、均方频率(MSF)、均方根频率(RMSF)、频率方差(VF)、频率标准差(RVF)。频域特征输入向量xp=[FC,MSF,RMSF,VF,RVF]。
[0110] 重心频率FC计算公式:
[0111]
[0112] 均方频率MSF计算公式:
[0113]
[0114] 均方根频率RMSF:
[0115]
[0116] 频率方差VF的计算公式:
[0117]
[0118] 频率标准差RVF计算公式:
[0119]
[0120] fi为i时刻的功率谱对应的频率值,s(fi)表示i时刻的功率谱值,其中i=1,2,3……n。
[0121] 2.3时频特征提取
[0122] 时频特征提取通过希尔伯特黄变换获取。时频特征输入xsp=[imf1,imf2,imf3,imf4],步骤如下:
[0123] 2.3.1初始化i=1,步骤一处理后的数据x(t),x(t)表示t时刻的数据组;
[0124] 2.3.2找到上述数据组中极大值点和极小值点,利用三次样条插值,求出信号上下包络线;
[0125] 2.3.3计算上下包络线的平均值m(t);
[0126] 2.3.4根据公式(15)得到信号分量h(t);
[0127] h(t)=x(t)-m(t) 公式(15)
[0128] 2.3.5判断信号分量h(t)是否满足以下两个条件,若满足则认为h(t)为imf分量,记为imfi,并进入步骤2.3.6;反之则跟新x(t)=h(t),重复步骤2.3.1-2.3.5;
[0129] (3)h1(t)的极值点个数和过零点数相等或相差一个;
[0130] (4)h1(t)上包络线和下包络线的平均值为0;
[0131] 2.3.6将上述imf分量从原始信号x(t)中分离出来,根据公式(16)获得残差分量h'(t)i+1跟新,x(t)=h'(t)i+1。重复2.3.1-2.3.6步骤,最终获得4个IMF分量:imf1,imf2,imf3,imf4。
[0132] h'(t)i+1=x(t)-imfi,i=1,2,3,4 公式(16)
[0133] 由于基于时域的特征提取可以对故障定位的优点,基于频域的特征提取具有更直观的特征信息优点,基于快速傅里叶变换的能够得到不同时刻的频谱的优点,而emd方法具有结合时间和频率的变化规律的优点,基于以上优点,选取这三种特征参数用于检测输入。
[0134] 步骤三、对以上所求时域、频域、时频域特征参数分别进行归一化处理,将数据的值转换到[0,1]区间内,使用的归一化函数为:
[0135]
[0136] x为步骤二计算得到的时域、频域、和时频特征,即[xs,xp,xsp]T作为输入数据,T为转置,max为输入数据中的最大值,min为输入数据中的最小值;y表示归一化后的数据,作为神经网络输入向量。
[0137] 步骤四、对步骤一处理后的数据初步提取未见异常的数据;根据公式(18)-(20)获得自协方差的高低阈值和变化率阈值,利用阈值判断对训练数据进行初步异常判断,若训练数据自协方差落在自协方差的高低阈值范围内且变化率落在零至变化率阈值范围内则为正常数据,反之则为异常数据;上述正常数据标0,异常数据标1,作为BP神经网络训练数据的标签;
[0138] 上述异常数据通过数据波动变化判断,为公知技术。
[0139] 自协方差的高阈值公式:
[0140] 自协方差的低阈值公式:
[0141] 变化率阈值公式:
[0142] 其中 为t时段初步数据的平均值,Varj为第j个t时间段内方差值,α、β、γ为调节权重参数,这里取值3、0.8、1.2,Varmax为各时段的方差中的最大值,ΔΦmax为各t时段变化率中的最大变化率,VΔΦ为各个t时段内变化率ΔΦ的方差值,μΔΦ表示各个t时段内变化率方差的平均值。
[0143] 步骤五、使用遗传算法,经过种群初始化、选择操作、交叉操作、变异操作和复制操作,通过适应度的计算,把获得的最优解作为BP神经网络的初始值权值和阈值。
[0144] 具体步骤如下:
[0145] 4.1、种群初始化,初始化种群规模为10,遗传代数30代,交叉概率为0.2,变异概率为0.1,随机初始生成一个可行解;
[0146] 4.2、根据上一步确定的BP神经网络结构,输入层15个节点,隐含层10个节点,输出层1个节点,确定种群个体的染色体个数为15*10+10*1+10+15=185;
[0147] 4.3、确定适应度函数,要求预测值和输出值尽可能小,使用的是期望值和输出值误差矩阵的范数作为适应度函数;
[0148] 4.4、选择操作,通过计算个体适应度选择适应度大的遗传给下一代;
[0149] 4.5、交叉操作,染色体选择通过计算个体适应度和适应度之和的比值确定,交叉位置是随机选择,但是交叉操作通过交叉概率决定。
[0150] 4.6、变异操作,变异是以一定概率0.1,用随机方法选择染色体中的一点进行变异,变异出更有秀的个体;
[0151] 4.7、复制操作,每次进化都保留上一代适应度最高的基因,所以需要原封不动的复制给下一代;
[0152] 4.8、适应度计算,通过计算适应度均值,当符合要求或者达到最大迭代次数,就进入下一步,否则继续执行4.2-4.7步;
[0153] 4.9、把获得的最优解作为BP神经网络的权值和阈值初始值。
[0154] 步骤六、把步骤五获取的权值和阈值用于训练BP神经网络。BP神经网络隐含层使用S型正切函数:
[0155]
[0156] 而输出为线性传输函数(purelin):
[0157] y=k·x+c 公式(22)
[0158] k为斜率,c为偏移值,这里k取1,c取0,既y=x。
[0159] 神经网络训练函数采用最速下降算法,设置允许迭代最大次数为3000,神经网络训练目标误差为0.001,学习速率lr=0.01。BP神经网络训练分为三个步骤:正向传播、误差反向传播、权重更新。
[0160] 权重更新:
[0161] 权值更新首先是从输出层开始反向传播误差,然后调整输入层权值,通过不断的更新,最终使得网络输出值和期望输出值误差不断减小,直到达到所需精度要求。每次调整都是沿着负梯度方向调整α个单位,因为负梯度方向是误差函数值下降最快的方向。权值调整公式:
[0162]
[0163] 通过训练优化后的遗传BP神经网络,得到一个检测识别率高的遗传算法优化BP神经网络。
[0164] 图2为本发明方法和未采用阈值判断的传统GA-BP方法对磁悬浮列车传感器测试数据的处理结果,本发明故障识别率明显高于传统GA-BP方法,识别率可达92.5%。
[0165] 上述实施例并非是对于本发明的限制,本发明并非仅限于上述实施例,只要符合本发明要求,均属于本发明的保护范围。