基于PMFCC

轨道电路是自动完成列车占用状态检查的铁路信号设备，其工作质量直接影响铁路运输的效率，也是列车安全运行的重要保障。轨道电路故障模式复杂，现场处理故障时间较长，严重影响铁路运输秩序，是铁路信号现场维护的难点。因此，探索各种智能诊断方法以缩短故障处理时间，提升轨道电路故障维护效率，正成为人们关注的热点。

作为我国优秀传统文化的重要内容之一，“为官清廉”极受世人赞誉和推崇。记载先贤“为官清廉”嘉言懿行的古籍可谓汗牛充栋，而我颇为激赏北宋苏轼的“廉为官程”一说。

目前，已有部分学者对ZPW-2000轨道电路故障智能诊断展开研究，如：基于神经网络算法[1-2]构造轨道电路相关变量与轨道电路相关故障的非线性映射关系，继而实现故障的识别分类；基于粗糙集和模糊认知图[3]对监测数据进行属性约简提取特征参数，构建分类器实现故障诊断；通过组合决策树[4]建立组合网络模型获取诊断规则，逐级压缩诊断范围最终完成多故障识别；利用相空间重构技术[5]对一维分路电流包络信号序列进行重构操作，提取补偿电容的故障特征，从而实现补偿电容故障的检测；基于集中监测系统获取的电压值或电流值[1，3-4]、机车或动检车采集的感应信号[2,5-9]，通过对稳态值进行分析处理，从而对轨道电路故障进行诊断。以上方法虽对故障机理明确或有较多历史故障信息的显性故障诊断精度较高，但在稳态条件下可能忽略正常与故障状态之间的暂歇性波动状态[10]，因其故障信息较少或无故障信息。然而如果只对正常和故障两种状态进行判别，很难满足铁路现场对轨道电路设备的诊断要求，因此对轨道电路暂态故障进行诊断尤为必要。

ZPW-2000轨道电路发生故障，其实质是轨道电路参数的突变，从而引发轨道信号幅值由暂态过渡到稳态[11]。若能基于集中监测系统对采集的轨道电路电压信号进行暂态分析，探究其暂态特性与轨道电路故障的关系，实现故障地面实时检测，则可提高诊断的便利性和时效性。

本文首先通过实地测试模拟故障获取数据，并对轨道电路故障时电压暂态波形进行分析，组合利用梅尔频率谱系数和动态时间规整算法建立轨道电路诊断模型；然后通过该模型利用主成分分析约简特征属性，再融入Fisher准则和K-means聚类选择合适故障模板，最终得到的模板具有清晰明显的特点。测试结果表明，基于集中监测系统与模型的使用，能够对轨道电路设备的多状态进行有效识别，提高诊断效率，为轨道电路的状态监测提供合理补充。

本文研究的时态RDF模型是用时间标签来标记RDF数据三元组，且表示有效时间。以下就是时态RDF模型的基本定理。

1 ZPW-2000轨道电路故障类别

1.1 ZPW-2000轨道电路组成

ZPW-2000轨道电路由轨面、发送端和接收端组成，其中发送端包括发送器、传输电缆和防雷单元，接收端包括接收器、衰耗器、传输电缆和防雷单元。发送端至接收端轨面等间距分布着补偿电容，相邻轨道电路区间由调谐区实现电气隔离。ZPW-2000轨道电路结构示意见图1。

图1 ZPW-2000轨道电路结构示意

在集中监测系统中，反映轨道电路性能状态的主要指标是电压，因此本文采用轨道电路电压数据进行故障诊断，以轨道电路的多故障系统作为研究对象。轨入电压在一段时间内的电压状态曲线见图2，其按照时间的先后顺序可以分成三部分：正常、波动、正常。对于ZPW-2000轨道电路，在不同时间、地点、环境下其电压曲线波动情况会略有不同，但是总体上还是一致的。

图2 轨入电压曲线

1.2 故障类型

为了验证轨道电路暂态理论的合理性，通过分析故障产生的原因[12-13]，查阅现场实际案例的相关资料，并结合现场工作人员的历史经验，在西南交通大学轨道电路试验基地进行了常见故障的测试。故障采集系统包括硬件和软件两部分。硬件部分进行通信、自动切换测试点和采集模拟量等；软件部分进行发送指令、数据采集和存储等。测试的核心板卡选用NI-USB 6251，采样频率可达10 MHz，单点定位输出时序差分数据。针对不同故障的暂态过程进行数据采集，输出的数据具有噪声较小、精确度高等特点。使用LabVIEW中的DAQmx控件做上位机，图形化编写简单。测试系统结构见图3。图3中，主机确认采集点后，由DAQmx控件把采集命令传给采集卡，通过采集卡与室内设备通信，切换单元控件自动选择采集点，并通过采集卡将采集的数据传入计算机中存储。轨道电路故障模式见表1。

图3 测试系统结构

表1 轨道电路故障模式

编号故障模式编号故障模式F1发送端故障F5调谐单元短路故障F2防雷单元故障F6主轨故障F3电缆网络模拟故障F7小轨故障F4空心线圈断线故障F8补偿电容断线故障

对表1中8种故障类型每种随机各选一组提取包络，见图4。由图4可以看出,F1～F8故障的表现形式差异较为明显：F1为发送端故障，轨入电压短暂趋于0；F2为防雷单元故障，设备短时间内电压波动幅度较大；F3为电缆网络模拟故障，设备短时间内电压波动趋于0；F4为空心线圈断线故障，轨出1电压瞬间幅度降低；F5为调谐单元短路故障，轨出2电压瞬间波动幅度较大；F6为主轨故障，轨出2电压短时间波动幅度较大；F7为小轨故障，轨出1电压短时间波动幅度较大；F8为单个补偿电容断线故障，主轨电压瞬间降低。

以上两组均在安装假体后，在髌上囊、膝关节两侧、鹅足注入“鸡尾酒”混合液，留置负压引流球后关闭伤口，缝合皮肤前抽尽关节内积血后往关节腔内灌注氨甲环酸（北京双鹤）60 ml。

图4 F1～F8的故障信号包络

图4展示了故障的发生会导致电压值的瞬间波动，工作状态发生跳变，可以较为直观地反映故障信息。但是，实际列车运行条件下设备会受到多方面影响，如天气、人为等因素，使故障波形不能直观精准反映故障位置，难以及时为维修人员提供有效的指导。因此，本文对ZPW-2000轨道电路的多种故障模式进行智能识别，组合使用PMFCC-DTW算法建立轨道电路诊断模型。

在NSCLC发生发展进程中，微环境内的肿瘤相关炎症反应会促进肿瘤细胞转移，而巨噬细胞的趋化聚集在微环境内的肿瘤相关炎症反应中起重要作用[4]。肿瘤相关巨噬细胞(tumor-associated macrophages,TAMs)是从血液或淋巴管迁移至癌巢的巨噬细胞，定位于肿瘤上皮组织、间质、坏死区域和肿瘤浸润的边缘。它们可以通过旁分泌作用产生并释放不同的生长因子、细胞因子、趋化因子等，促进血管生成、组织重建及损伤修复作用，在肿瘤发生发展过程中起促进作用[5]。根据TAMs活化后不同的表型和功能可分为促炎症反应的M1型巨噬细胞和抗炎症反应的M2型巨噬细胞[6]。

2 诊断模型基本原理

2.1 梅尔频率谱系数

梅尔频率倒谱(Mel-Frequency Cepstrum Coefficient, MFCC)在语音处理中表示一个语音的短时功率谱，由对数功率谱在频率的一个非线性梅尔刻度上进行线性余弦转换所得[14]。梅尔频率倒谱fmel与频率f的关系为

(1)

轨道电路在复杂的工况下，暂态故障通常是潜在的故障或早期的微弱故障，具有故障特征不明显、易被噪声影响的特点，使故障难以被发现。MFCC是按临界带宽由密到疏形成一组带通滤波器，对样本信号进行滤波，可以有效抑制非本频带的高频噪声成分。并且通过大量试验发现，当轨道电路发生故障时，不同状态下同一频带的信号频率基本一致，使用MFCC提取得到的不同状态下的特征各不相同，因此可以将MFCC用于轨道电路的故障诊断中。具体诊断过程如下：

Step1 滤波输出为

H(z)=1-αz-1

(2)

Step2 输出对数能量为

(3)

1≤m≤M

Step3 提取梅尔系数为

(4)

式中：n为梅尔系数的数量；Ci(n)为第i帧的第n个MFCC系数；s(m)为对数能量模块的输出；M为Mel滤波器的个数。

由上述过程，依次获取8种故障的总采样点数，分别计算出每帧MFCC参数Sm和差分系数特征集Cn，约简后得到特征矩阵Mmel为

Mmel=|Sm Cn|

(5)

通过MFCC提取的特征能够精准地捕捉信号变化，但是信号特征维度较高，不利于可视化分析，更会导致诊断时间增长。因此，本文利用主成分分析(Principal Component Analysis，PCA)，对高维数据进行压缩，把可能具有相关性的高维变量合成线性无关的低维变量，新的低维变量集会尽可能保留原始信号的特征。组合MFCC和PCA，寻找和保留最有效、最重要的信号特征，从而降低计算复杂度、减少数据存储量，并且提高信号识别的效率。

人才培养模式是民办高校定位和办学特色的关键，是高校根据人才培养目标，选择适合目标人才的培养内容，并在一定培养制度保障下保证实施的人才培养体系，主要包括了培养目标、培养内容、培养制度和培养方式四个方面［5］。独立学院成立之初，多依托母体大学的优质教育资源设置专业、教育教学、培养人才，在发展壮大的过程中呈现出与母体大学的“同质化”现象，难以形成具有自身特色的人才培养模式，进入了发展的瓶颈期。

2.2 动态时间规整

在模式匹配中，通常选用动态时间规整(Dynamic time warping，DTW)求解衡量不同信号间似然度的测度函数，也可以不要求每条故障数据长度严格对齐，允许在时间轴上的拉伸或压缩[15]。本文通过计算模板库和测试集中对应帧之间的失真距离，得出帧匹配距离矩阵，编制DTW表格。表格中的格点是对应的交汇点，按照距离最短原则搜索格点，反复递推得到组最优解，并将累积距离最小的归为此类故障。

将实测数据分为测试故障和故障模板，分别为X=(x1,x2,…,xm)和Y=(y1, y2,…,yn)，长度分别为m和n。其中xi与yj之间的距离dij采用欧式距离计算

dij=d(xi,yj)=|xi-yj|

(6)

计算得出两序列对应点间的距离矩阵D，并在D中确定出一条路径W=(w1,w2,…,wk)让X与Y匹配程度最大。为了求出有效路径，需要符合的约束条件见表2。

本次研究结果中，实验组治疗两周后，NT-ProBNP水平更低，LVEF更高，6分钟步行距离更长，治疗总有效率也比对照组高，说明实验组患者治疗后心脏功能、活动能力都优于对照组患者。

表2 约束条件

属性条件有界性max(m,n)≤K≤m+n-1边界约束起始点(1,1),终点(m, n)连续性单调性wk=(ak,bk), wk+1=(ak+1,bk+1)0≤ak+1-ak≤1, 0≤bk+1-bk≤1

计算得出多条路径，其中最优路径为

(7)

综上所述，本模型的基本原理为，首先针对MFCC算法进行优化，然后结合PCA降维的梅尔频率倒谱(以下简称为“PMFCC算法”)，并组合DTW算法诊断出故障类型，简称为“PMFCC-DTW算法”。

3) 控制段表示数据的字节数，在通用板卡向分油机板卡的数据发送中，LED状态指示的数据以及LCD显示数据是包含在一个数据报文中的，其中液晶显示为20个字符用ASCII码表示，每个字符占用一个字节，另加显示的第几行占用一个字节，10只LED状态指示占用一个字节，因此控制段此时为22个字节。

3 轨道电路特征参数提取

本文将PMFCC特征提取过程分为两部分：梅尔频率系数和主成分分析。通过梅尔频率系数获取MFCC特征值后利用PCA进行降维处理，最终得到轨道电路特征参数。轨道电路特征提取流程见图5。

图5 轨道电路特征提取流程

选用图4中8种不同故障的信号进行MFCC特征提取，得到24维的MFCC信号特征，再进行PCA降维。通过求得各主成分φk的贡献率μk，计算出降维满意度。假如前k个主成分的累计贡献率δk大于阈值(通常为85%)，则选取前k个主成分作为新的样本主成分进行循环。其中，第k个主成分贡献率为

(8)

通过计算得出符合阈值的特征维数并不唯一，为更好地选出最有效的特征维数，本文利用费舍尔判别准则(Fisher准则)对满足累计贡献率δk>85%的特征参数直接进行二分类求解，其最优参数求解公式为

由此可知,基于特征函数法的Merton跳-扩散模型的看涨期权定价公式中含有4个参数,即波动率σ、Poisson过程的强度λ以及随机跳的2个参数μJ和σJ。

(9)

式中：Sa为类内散度矩阵；Sb为类间散度矩阵。

设d维度空间的两类样本为X1和X2，且样本个数为N，当d维度空间的样本映射到一维Y空间时得到y=wtx。其中，类内散度矩阵与类间散度矩阵分别为

(10)

(11)

式中：

为d维度空间样本Xi的均值向量；

为一维Y空间的两类样本均值；w依据拉格朗日乘子求解，

用

作为分界点，第i个样本投影点的yi与y0比较，yi>y0则属于第一类，否则属于第二类。将F1～F8的8种故障两个一组不重复地排列组合形成28组进行分类，统计分类结果的正确率，见图6。

式中：E为油液体积弹性模量；V0为溢流阀下腔容积与连接溢流阀管道容积之和；Q1为流经阻尼孔R1的流量；Qy为主阀口流出的流量。

图6 Fisher准则分类正确率

从图6可以看出，8种故障的特征数据在不同维度下分类的正确性有很大不同。计算不同维数分类正确率的平均值，6维为78.54%、8维为87.39%、10维为81.19%、14维为78.22%、18维为74.09%，可以看出8维的分类正确率最高，随着维度增加或减少，特征数据的分类效果降低，差异性变小。因此，本文以8维数据作为信号特征。

进一步对降维后PMFCC算法得到的8维条件下信号特征与MFCC算法的24维信号特征进行比较分析，利用欧氏距离计算不同故障间的相对距离，通过距离的大小分析不同故障间的差异性。选定F1故障作为参考与F2～F8故障进行比较计算，然后归一化处理得到雷达图，见图7。

图7 F1作为参考的雷达图

从图7中可以看出，8维条件下PMFCC算法的特征数据差异性高于24维条件下MFCC算法，并且F1与F2～F8故障的相对距离较为均匀，不易受到外界干扰。同理，将F2～F8分别作为参考，计算得到的结果与图7类似。8维条件下的PMFCC特征数据间差别明显，为后续的模板训练和故障识别，既减少了计算量又降低了模板库的存储量，同时还将特征值进行最优化，不但提高了信号的识别精度也提高了信号的识别效率。

4 故障诊断模型

4.1 故障模板选取

本文利用DTW计算模板信号与测试信号在每一帧的Mel倒谱距离，通过该距离寻求待检测信号与模板信号的最佳匹配路径，从而实现测试信号与模板信号间的时差失真最小、相似特征最大。诊断流程由增量式训练和匹配识别两个阶段组成，见图8。通过历史故障数据训练出最优的故障模板作为匹配依据，然后将实时故障数据提取特征后与故障模板进行匹配，并将结果在人工复合后作为新的历史数据反馈给故障模式库进行训练形成新的故障模板，从而实现自适应效果。

图8 诊断流程

在实际运用中，增量式的故障模板库能够有效更新报告设备故障情况，提高诊断精准度，因此模板库的选取意义重大。本文使用K-means聚类与最小欧式距离结合的方式找出最优的不同故障的模板库，并且根据实时数据与故障模板进行对比后反馈给故障模式库进行训练形成新的各类故障，从而保证算法实际应用的有效性。K-means聚类主要步骤如下：

Step1 输入故障集{X}，有n个样本和k个故障模式,分类{Gi， i=1,2,…,K}。

什么是语文教育现代化？围绕立德树人，培养创新型人才，满足国家、社会高速发展的需求为目标，建造具有中国传统、中国特色的母语教学课程，构建适应新时代社会所需要的现代语文教学结构，从而培养学生在现代社会所需要的母语交流必备品质、基本能力和正确的价值观，即语文核心素养。

Step2 随机选取k个初始聚类中心{c1,c2,…,ck}，依次计算故障集{X}到聚类核心的距离D(xi,cj)，按照最小距离原则形成k个聚类核心，计算式为

（二）在地理课堂教学中运用生活实际，提升地理思维能力。针对学生出现的地理思维缺陷，特别对高中文科生而言，由于面临高考的现场，加上许多高考题目都是全新的，没有良好的地理知识思维结构在做题的时候是很难正确解决问题的。所以在地理课堂教学中更要注重加强诱导启发，调动学生的地理思维，注重启发诱导，调动学生一切可以利用的头脑中的地理理论知识，从学生的认知能力角度开启，一步一步深入挖掘，让学生深层次的进行思考形成正确的地理思维。此在此基础上老师就应该深入挖掘，提出问题，然后让学生会提出一系列的假设，最终回归到目前比较易懂好理解的理论知识上，充分的让理论知识生动起来。

Dij=‖xi-cj‖ i=1,2,…,n j=1,2,…,k

(12)

Step3 计算新的聚类中心

计算式为

(13)

式中：Nj为第j个聚类中心的样本个数。

Step4 利用式(14)～式(15)聚类准则数值E判断聚类结果是否合理，若合理迭代终止，否则重复Step2和Step3。

民营企业没有各种发展的支持，拥有的只是野蛮生长的力量和一颗不断向上的心——于晟图机械而言，发展路上事无巨细，将7000多天的始终如一，化为了所谓固执的坚持。而这正是像晟图机械这样的民营企业，能够逐步走向成功、推动行业不断前进的内生动力。

(14)

|E′-E|<ζ

(15)

4.2 结果分析

从采集的故障文件中随机抽取数据63条，其中发送端数据7条，防雷单元数据5条，电缆网络模拟数据9条，空心线圈断线数据6条，调谐单元短路数据10条，主轨出数据8条，小轨出数据7条，补偿电容断线数据11条。将随机抽取的数据使用K-means聚类得出各故障参考模板，再利用PMFCC-DTW故障诊断把t1～t9测试数据与F1～F8模式库中的参考样本进行匹配，求出部分动态匹配距离，见表3。其中匹配计算得到的数值越小，则距离越近，相似度就越高，并将该测试数据归为此类故障。例如表3中圈出的数据，t1所在列中最小的值是与F6相对应，故t1是F6故障；t4所在列中最小的值是与F1相对应，故t4是F1故障。同理，按此计算也可以得到动态匹配结果。

表3 部分动态匹配距离

故障t1t2t3t4t5t6t7t8t9F139 252.9517 976.9317 976.9310 661.8620 163.3417 897.9121 041.8317 557.0916 576.64F238 765.1217 976.9317 976.9317 976.9315 921.8419 937.7618 390.1618 385.3221 887.41F327 796.0117 976.9317 976.9317 976.9316 477.8416 728.044 355.8114 800.4120 057.92F428 882.7817 976.9317 976.9317 976.9317 757.9617 446.8414 187.2416 497.8919 770.62F526 294.9417 976.9317 976.9317 976.9318 135.0716 505.5611 586.5814 372.5019 545.84F617 567.7117 976.9317 976.9317 976.9344 832.4133 575.9736 022.7433 289.4729 648.71F717 976.936 397.327 556.6814 837.8817 976.9317 976.9317 976.9317 976.9318 460.01F828 189.5311 241.6317 976.9315 542.3629 571.8320 194.3726 435.176 016.5318 073.35

为了证明本文采用的PMFCC-DTW算法进行故障诊断的有效性，根据实测ZPW-2000轨道电路的数据建立数据库，选取2 400条数据，做改进效果对比实验分析。其中发送端数据、防雷单元数据、电缆网络模拟数据、空心线圈断线数据、调谐单元短路数据、补偿电容断线数据、主轨出数据和小轨出数据各300条。PMFCC-DTW模型与MFCC-DTW模型的诊断效果对比见表4。对比结果显示，PMFCC-DTW模型单个故障的识别率高达96.67%，而且对每一种故障的识别率都要高于MFCC-DTW模型，正确率有明显提升。

表4 故障识别率 %

故障编号PMFCC-DTW故障识别率MFCC-DTW故障识别率F195.3384.67F289.3485.00F390.3377.33F488.6787.33F589.0086.00F696.6781.33F789.3379.33F890.0085.00

另一方面利用混淆矩阵(Confusion Matrix)计算实际样本值和模型诊断值的正确率、精准率和召回率，并统计模型诊断耗时。

脑子里所想的与课堂能讲的“表里不一”的困惑；先进的特色理论与复杂的社会现实的差距；自由与纪律的约束，等等，由此产生的心里纠结，直接或间接地诉之于学生。

正确率为

(16)

精准率为

Precision=TP/(FP+TP)

(17)

召回率为

Recall=TP/(FN+TP)

(18)

式中：TP为样本的真实类别是某类故障，并且模型诊断的结果是同类故障的个数；FN为样本的真实类别是某类故障，但是模型将其诊断成为其他故障的个数；FP为样本的真实类别是正常，但是模型将其诊断成为故障的个数；TN为样本的真实类别是正常，并且模型诊断的结果也是正常的个数。

基于改进效果对比实验的2 400条数据以及1 000条正常数据，仿真了HHT-DBWT算法和SVM算法，并与本文改进的PMFCC-DTW算法进行对比分析，评价结果见表5。

表5 评价结果

算法正确率/%精准率/%召回率/%时间/sPMFCC-DTW92.2497.7691.08116.48HHT-DBWT90.0697.4788.21134.56SVM84.3596.9880.33128.31

由表5可见，PMFCC-DTW算法的正确率较高，而且精准率高达97.76%，表示本模型诊断故障的准确度很高。通过召回率分析出模型诊断故障可信性强，且用时短。

在实际应用中，通过安装轨道电路室内外设备信号采集装置，结合信号集中监测系统获取室内外监测信息，利用本文方法对轨道电路信号实时监控实时诊断，若能获取较完整的监测数据和突发故障数据以便故障模式更深入研究，对提高轨道电路故障诊断及日常维护效率具有重要意义。

5 结论

(1)本文针对高速铁路轨道电路常见的故障，通过实地测试采集到的暂态电压数据按时域、频域进行智能化区段划分，利用MFCC提取出表征轨道电路状态的特征参数，并基于主成分分析进行特征降维，使用Fisher准则选出最优特征集，通过K-means聚类进行故障集训练，得到轨道电路暂态数据的初始模板，形成轨道电路的模板库，最后通过DTW模型计算出测试集与每种模型的匹配距离，从而进行故障诊断。

(2)通过实验以及与其他诊断方法进行对比可以发现，基于PMFCC-DTW的轨道电路故障诊断的识别时间大为缩短，正确率达到90%以上，为轨道电路设备健康状态监测提供了一种新思路。

参考文献：

[1] 曹雲梦, 边冰, 张天奇. 基于BP神经网络的轨道电路故障诊断[J]. 华北理工大学学报, 2018, 40(1): 78-82.

CAO Yunmeng, BIAN Bing, ZHANG Tianqi. Track Circuit Fault Diagnosis Based BP Neural Network[J]. Journal of North China University of Science and Technology, 2018, 40(1): 78-82.

[2] 陈姝姝，田慕琴，宋建成．基于PSE-ANFIS的ZPW-2000A轨道电路故障诊断方法 [J]．太原理工大学学报，2018，49(6)：898-902.

CHEN Shushu，TIAN Muqin，SONG Jiancheng．Fault Diagnosis Method of ZPW-2000A Track Circuit Based on PSE-ANFIS [J]．Journal of Taiyuan University of Technology，2018，49(6)：898-902.

[3] 董昱，陈星．基于粗糙集和模糊认知图的ZPW-2000轨道电路故障诊断 [J]．铁道学报，2018，40(6)：83-89.

DONG Yu，CHEN Xing．Intelligent Fault Diagnosis for ZPW-2000 Track Circuit Based on Rough Set Theory and Fuzzy Cognitive Map [J]．Journal of the China Railway Society，2018，40(6)：83-89.

[4] 朱文博，王小敏．基于组合决策树的无绝缘轨道电路故障诊断方法研究 [J]．铁道学报，2018，40(7)：74-79.

ZHU Wenbo，WANG Xiaomin．Research on Fault Diagnosis of Railway Jointless Track Circuit Based on Combinatorial Decision Tree [J]．Journal of the China Railway Society，2018，40(7)：74-79．

[5] 孙上鹏，赵会兵．基于相空间重构的无绝缘轨道电路补偿电容故障检测方法 [J]．铁道学报，2012，34(10)：79-84.

SUN Shangpeng，ZHAO Huibing．The Method of Fault Detection of Compensation Capacitor in Jointless Track Circuit Based on Phase Space Reconstruction [J]．Journal of the China Railway Society，2012，34(10)：79-84.

[6] 王秋实，刘扬，王小敏．基于车载单元和IPSO-SVM的轨道电路分路不良诊断方法 [J]．系统仿真技术，2016，12(4)：291-296，306.

WANG Qiushi，LIU Yang，WANG Xiaomin．Diagnosis of Shunt Malfunction of Track Circuit Based on TCR and IPSO-SVM [J]．System Simulation Technology，2016，12(4)：291-296，306.

[7] 傅佳伟，王小敏，郭进．ZPW-2000轨道电路的多轮对动态分路建模研究 [J]．铁道科学与工程学报，2018，15(9)：2374-2384.

FU Jiawei，WANG Xiaomin，GUO Jin．Research on Dynamic Multi-wheel Sets Shunted Track Circuit Modeling [J]．Journal of Railway Science and Engineering，2018，15(9)：2374-2384.

[8] 李亚兰，董昱，南接龙．基于CEEMD特征提取的无绝缘轨道电路补偿电容故障诊断 [J]．计算机测量与控制，2015，23(2)：378-381.

LI Yalan，DONG Yu，NAN Jielong．Jointless Track Circuit Compensation Capacitor Fault Diagnosis Based on CEEMD Feature Extraction [J]．Computer Measurement & Control，2015，23(2)：378-381.

[9] 苏丽娜，张亚东，王梓丞，等．基于变分模态分解和能量谱的轨道电路分路不良故障监测 [J]．铁道标准设计，2019，63(2)：155-160.

SU Lina，ZHANG Yadong，WANG Zicheng，et al．Prediction of Shunt Malfunction of Track Circuit Based on Variational Mode Decomposition and Energy Spectrum [J]．Railway Standard Design，2019，63(2)：155-160.

[10] 韦达明．轨道电路红光带分析方法新思路 [J]．铁道运营技术，2017，23(3)：40-42.

WEI Daming．A New Method of Red Band Analysis for Track Circuits[J].Railway Operation Technology，2017，23(3)：40-42.

[11] 赵斌，张友鹏．基于傅里叶变换结合Q-D算法的轨道电路暂态分析[J]．铁道学报，2016，38(3)：78-83.

ZHAO Bin，ZHANG Youpeng．Transient Analysis of Track Circuit Based on FFT & Q-D Algorithm[J]．Journal of the China Railway Society，2016，38(3)：78-83.

[12] 付淳川, 朱文博. 基于粗糙集—C4.5的轨道电路故障诊断方法研究[J]. 铁路通信信号工程技术, 2018, 15(3): 22-27.

Fu Chunchuan, Zhu Wenbo. Research on Railway Track Circuit Fault Diagnosis Based on Rough set and C4.5 Decision Tree Algorithm[J]. Railway Signalling & Communication Engineering, 2018, 15(3): 22-27.

[13] 陈欣．基于粒子群支持向量机的轨道电路故障诊断[J]．铁路计算机应用，2016，25(8)：56-63.

CHEN Xin．Track Circuit Fault Diagnosis Based on Particle Swarm Optimization and Support Vector Machine[J]．Railway Computer Application，2016，25(8)：56-63.

[14] 樊晓鹤，赵鹤鸣，陈雪勤，等．倒谱参数稀疏分解下的汉语音谎言检测[J]．声学学报，2018，43(1)：121-128.

FAN Xiaohe，ZHAO Heming，CHEN Xueqin，et al．Deceptive Chinese Speech Detection Based on Sparse Decomposition of Cepstral Feature [J]．Acta Acustica，2018，43(1)：121-128.

[15] 曹源，王玉珏，马连川，等．基于DTW的车辆轴温监测方法 [J]．交通运输工程学报，2015，15(3)：78-84，100.

CAO Yuan，WANG Yujue，MA Lianchuan，et al．Monitoring Method of Vehicle Axle Temperature Based on Dynamic Time Warping [J]．Journal of Traffic and Transportation Engineering，2015，15(3)：78-84，100.