高铁钢轨伤损分类方法
【技术领域】
[0001] 本发明涉及高铁钢轨伤损检测领域，尤其涉及一种高铁钢轨伤损分类方法。
【背景技术】
[0002] 随着科技进步，高铁运输技术快速发展，高速铁路的安全运营面临严峻的挑战。忽 略影响高速铁路正常运行的人为因素，列车车辆和钢轨的状况对于车辆安全运行具有重要 影响。事实上，钢轨伤损情况是铁路运输出现安全事故的主要原因。钢轨伤损常见的类型 有核伤、纵向裂纹、水平裂纹等很多种。尤其钢轨在列车运动中受到强烈冲击、挤压等会更 大程度的影响钢轨的健康状况。因此发展快速、准确的钢轨伤损分类技术对于高铁的安全 运行至关重要。
[0003] 超声波技术很早就用于钢轨伤损的检测，包括基于传统超声技术探伤车等设备， 但超声检测技术容易受到钢轨表面情况的影响，其探测速度也无法满足高铁的伤损检测速 度需求。基于超声导波技术的钢轨探伤法采用长波长和低频的超声波来检测轨道，能实现 大范围快速检测，但对微小的伤损效果不好，只有当伤损发展到一定程度才能明显检测。基 于电磁超声技术的钢轨检测法是利用电磁效应车产生超声波来检测钢轨伤损，但容易引入 噪声，且超声转换效率低，其他图像技术等方法也都存在着性价比低和不能对整个轨道的 伤损准确分类等问题。可见，目前已有方法无法满足高铁钢轨伤损分类检测的速度和准确 性的要求。
[0004] 算法用于整个钢轨实时监测系统的不同阶段，而当钢轨出现伤损时利用所提出的 非负张量分解和极限学习机的钢轨伤损分类方法在信息中心对钢轨出现何种伤损进行详 细分析，取得良好速度和识别率且具有一定鲁棒性。
【发明内容】
[0005] 针对目前已有的钢轨伤损分类方法存在处理速度慢和准确性差的问题，提出一种 基于非负张量分解和极限学习机的钢轨伤损分类方法。本发明的基本思路是首先利用小 波分析的方法提取有损信号的时域和频域局部特征，对同一测量点结合不同车厢建立三维 张量信号，将数据扩展到高维空间，得到非负张量，引入奇异值分解对非负张量分解进行改 进，用改进的非负张量分解算法提取隐藏的信息，针对非负张量分解速度慢的问题，引入极 限学习机来对伤损信号进行识别分类。
[0006] 在信号处理时，信息中心收到的振动信号是非稳定、非线性信号，故利用小波分析 的方法提取其特征，对得到的张量引入非负性的限制条件，使在精确度和解释度都能达到 很好效果具有实际意义。利用高维张量模型的TUCKER分解可有效地实现模型降维及特征 提取。但TUCKER分解的结果只有当核张量在超对角线位置上存在非零量时，且满足1 1 = I2= I 3的条件，TUCKER分解转变为CP分解，才具有唯一分解结果，故利用TUCKER分解在 数据压缩上优势，并对其因子矩阵和核张量施加正交性、稀疏性及非负性约束，则得到的CP 分解就会变为非负张量分解（NTD)。本发明采用非负张量分解方法进行振动信号的特征提 取。
[0007] 奇异值分解（SVD)能对数据中变化最大的张量进行识别且提高其收敛速度和鲁 棒性，本发明采用奇异值分解的方法进行非负张量的初始化。
[0008] 本发明提供一种高铁钢轨伤损分类方法，对将时域特征、频域特征和多个车厢特 征综合考虑而利用小波分析方法得到的特征数据进行非负张量分解，选交替最小二乘法作 为其迭代准则，在用SVD算法改进NTD分解，提高其收敛速度，利用改进的NTD特征提取方 法提取特征，接着引入极限学习机算法实现了对钢轨伤损的实时分类，分类识别算法不仅 速度快准确率高且具有鲁棒性。其特征在于：其中对钢轨伤损分类的方法是按照下面的步 骤完成的：
【主权项】
1. 一种高铁钢轨伤损分类方法，对将时域特征、频域特征和多个车厢特征综合考虑而 利用小波分析方法得到的特征数据进行非负张量分解，选交替最小二乘法作为其迭代准 贝1J，在用高阶奇异值分解改进非负张量分解，提高其收敛速度，利用改进的非负张量分解特 征提取算法提取特征，接着引入极限学习机算法实现了对钢轨伤损的实时分类，分类识别 算法不仅速度快准确率高且具有鲁棒性，其特征在于对钢轨伤损分类的方法是按照下面的 步骤完成的： (1) 首先采集不同钢轨伤损类型信号建立数据库，采用小波分析方法处理通过带通 滤波器的频率点，构建时域、频域、车厢的张量信号，针对信号局部特征取函数簇为 b(t)}，函数表达式为
，进行连续小波变换得出Morlet小波函数为
(2) 对所得张量进行TUCKER分解得到的因子矩阵和核张量施加正交性、稀疏性约束， 使核张量只在超对角线的位置上存在非零量且满足I 1= I 2= I 3的条件，对此模型施加非 负性约束，变为非负TUCKER分解，非负TUCKER分解简写为NTD ; (3) 选交替最小二乘算法作为NTD的迭代准则提高其收敛速度，交替最小二乘算法 准则为：
寻找使得误差张量最小的核张量殳以及因子矩阵
(4) 采用高阶奇异值分解来对所得核张量及因子矩阵进行优化，分解过程分为两步； a) 对于η = 1，2, -·，Ν，计算张量A的展开矩阵A(n)，并对其进行奇异值分解，使得A(n) =U(n)S(n)V(n)T，正交矩阵U(n)是A (n)的主左奇异向量，交替地对协方差矩阵进行特征值分解，
b) 利用反演公式来计算核张量，当因子矩阵和核张量都是正交的时候，S= AX1Uw tX2U(2)t…X nU(n)t; (5) 利用NTD进行特征提取，包括训练和测试两个步骤；先将原始数据集A分为两部 分，得到一大一小两个数据集，其中把较大的数据集作为训练数据集B，将较小的数据集当 作测试数据集C ;接着对训练集进行训练，得到特征空间以及训练集特征；然后，将测试集 的数据从原空间投影到训练所得到的特征空间，得到测试集的特征； (6) 将训练集特征和测试集特征作为参数输入极限学习机得到分类结果；前述的基于 NTD的特征提取算法中训练和测试采用下边步骤计算完成： a) 训练是为寻找N个因子矩阵或基矩阵X(n)以及核张量￡(k)，假设数据集有均为N阶 张量的K个样本，记e ，A(N+1)看作联结张量B € N+1维的矩阵化 形式，子张量通过确定张量S的第Ν+1维的k值来得到，_
:)，对N+1维核张 量G '的子张量￡(k)通过确定张量
等价于张量积B形式 PX1XwX2X^…XNX(N);分解得到核张量再将其矩阵化，得出训练数据集的特征； b) 在测试阶段，把基矩阵X(k)固定，再求解核张量变为一个非负约束最小二乘问 题，采用交替最小二乘算法来迭代直到收敛，以得到测试集的特征数据； 前述的极限学习机分类识别部分采用下边步骤计算完成： a) 极限学习机算法简称ELM，算法中隐藏结点参数随机指定，其隐藏层输出用行 向量h(x) = Di1(X),…，hjx)]表示，L为隐藏结点个数，X为输入样本；假定训练 数据集有N个训练样本（ Xi，tj，Xi为输入样本，％为标签，单隐层神经网络数学模 型Ηβ =T中，H为隐藏层的输出矩阵，β为输出权值，T为目标向量；ELM的输出为
^其中a」为第j个隐藏结点的输入权值;b」为第i个隐藏 结点的偏差；为第j个隐藏结点的输出权值；Xi)表示第j个隐藏结点的输出函 数；且h (Xi)为隐藏层关于Xi的输出向量，能够将d维输入空间映射到L维特征空间； b) ELM求解仅需求解一个极小范数的最小二乘问题，采用Moore-Penrose广义逆来求 解；ELM的目标是使训练误差最小及不同类别数据的间距最大，ξ/= [ξ u，…，Ii,Jt是 m个关于训练样本Xi的输出结点的误差，这里C是一个正则规划参数用来平衡训练误差最 小化和边缘距最大化这两个指标；ELM的最优化问题为： 最小化
关于：Ii(Xi) β = tj-ξ/，i = 1，2，…，NELM 的二 次优化问题为：
μ为引入 ' ^ I =1 / =1 j-\ 的拉格朗日乘子，a = [Ci1，…，αΝ]τ，求偏导，可得输出矩阵β为
故ELM面对二分类问题时的决策方程为
同时ELM用一个分类器实现多类别数据分类，arg Hiaxfi(X) ie[l，2, ···，!!!]为其对 应决策方程，m为类别的数目。
【专利摘要】本发明提供一种高铁钢轨伤损分类方法，其主要思想是：首先利用小波分析方法提取有损信号的时域和频域局部特征，对同一测量点结合不同车厢建立三维张量信号，将数据扩展到多维空间得到非负张量，采用交替最小二乘算法作为非负张量分解的迭代准则，接着引入奇异值分解对非负张量的初始化进行改进，利用改进的非负张量分解方法提取隐藏的特征，最后引入极限学习机算法实现对钢轨伤损的实时分类。本发明方法可以准确对钢轨伤损信号进行分类，提高了伤损分类的速度和准确性且具有较好鲁棒性。本发明提出的基于钢轨伤损分类方法优于已有方法，可以获得更好的识别效果，可以在钢轨伤损分类领域广泛应用。
【IPC分类】G06K9-62
【公开号】CN104751169
【申请号】CN201510016567
【发明人】陈玉敏, 马立勇, 孙明健, 王胜利
【申请人】哈尔滨工业大学（威海）
【公开日】2015年7月1日
【申请日】2015年1月10日