本发明公开了一种基于趋势分段相似性的机场噪声监测点异常识别方法，属于机场噪声监测点异常分析技术领域。
背景技术：
随着机场噪声环境污染程度日益严重，噪声环境影响投诉逐年增多，机场噪声问题成为影响民航业可持续发展的障碍之一。机场噪声主要是航空器在民航机场内进行起飞、降落和上升、下降等重要操作步骤时产生的噪声，它会对机场尤其是周围环境造成一些负面影响。机场噪声监测既是国家、地方法律法规对机场管理机构的要求，也有助于机场当局了解航空器噪声的影响程度和范围，掌握低噪声程序的执行情况及减噪效果。为了客观反映机场航空噪声影响程度，同时为日后的噪声管理控制提供依据，大部分机场投入大量资金建立了机场噪声监测系统，在关键区域设置噪声监测点。
然而，这些硬件设备不可避免地存在损坏或老化等问题，导致监测点无法准确采集所处区域的噪声数据。现有噪声监测数据文件中发生的多字段连续缺值或出现奇异值正反映了这个问题。另外，由于噪声直接由航空事件产生，起降飞机的设备(消音飞行程序、油门动力系统、襟翼、机身气流摩擦等)异常也会导致噪声监测数值发生异常。因此，当监测点出现异常数值时，如何通过历史监测数据以及其他正常节点采集到的数据来发现异常就成为一个有价值的研究工作。在识别了异常噪声数据后，相关部门可以及时了解监测节点工作状况，追究造成异常的原因(航空器故障或者监测设备故障)，并采取相应措施，避免更大损失。
目前国内的机场噪声监测系统尚处于初级阶段，有关机场噪声异常检测的研究基本空白。考虑到机场周边每个监测点获得的数据是机场噪声在不同时间点的数值，且按时间顺序排列成序列，所以机场噪声异常监测问题可以用时间序列异常检测技术来解决，时间序列相似性度量技术就是其中的一种。
技术实现要素：
本发明针对机场噪声监测中存在的数据异常问题，提出了一种基于趋势分段相似性的机场噪声监测点异常识别方法，以在复杂的机场噪声环境中快速有效地识别异常噪声监测点，以辅助故障定位。
本发明为解决其技术问题采用如下技术方案：
一种基于趋势分段相似性的机场噪声监测点异常识别方法，包括如下步骤：
步骤1：利用机场噪声监测设备获取机场周边多个噪声监测点的原始时间序列；
步骤2：对原始噪声时间序列进行预处理，创建标准的噪声时间序列数据集；
步骤3：用基于趋势分段的时间序列特征表示方法对每个监测点的噪声时序数据进行降维表示；
步骤4：使用基于趋势分段的时间序列相似性度量方法，逐一度量各监测点的噪声时序与其他监测点的噪声时序之间的相似程度，并建立相似性矩阵；
步骤5：根据相似性矩阵，找到每一个监测点的前k个最相似的监测点，创建相似监测点集合；
步骤6：度量各个监测点新的噪声时序与其相似监测点的新噪声时序之间的相似性，若相似性发生明显变化，则判定为异常。
步骤2中所述对原始噪声时间序列进行预处理为将实测噪声数据集作归一化处理。
步骤3中所述基于趋势分段的时间序列特征表示方法具体如下：对噪声数据集中的每一条噪声序列进行了重新表示，运用了迭代终点拟合算法，将各个监测点的噪声时序表示成由多个趋势分段组成的降维序列，每个分段用一个三元组表示。
所述三元组的内容包括：分段的起点、分段均值的符号化表示和分段趋势变量。
步骤4中所述基于趋势分段的时间序列相似性度量方法具体如下：使用符号化距离和欧氏距离分别度量噪声序列中各分段的均值和趋势变量之间的距离，将均值距离与趋势距离结合，定义了一种能度量不同长度序列相似性度量函数。
步骤6中所述相似性发生明显变化为相似性明显偏大。
本发明的有益效果如下：
1、本方法具有智能学习能力，识别所需的机场噪声时间序列数据易于获取。
2、采用基于趋势分段的时序表示方法，结合迭代终点拟合算法，依据机场噪声时间序列的不同趋势对其进行分段，有效保留了噪声时间序列中重要的变化趋势信息，能更加精确地反映相应噪声监测点监测到的噪声变化趋势。
3、采用基于趋势分段的时间序列相似性度量方法，结合了动态时间弯曲距离的优点，能更好地捕获机场噪声时间序列的模式特征，提高异常检测的准确性和可靠性。
4、本方法在复杂的机场噪声环境中，能有效地提高机场噪声监测点异常识别的及时性和有效性。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图对本发明创造作进一步详细说明。
本发明是一种基于时间序列相似性度量的机场噪声监测点异常识别方法流程如图1所示，具体包括以下步骤：
步骤1：利用布置在机场周边的噪声监测设备(丹麦的b&k声学测量仪器)获取机场周边多个噪声监测点的原始时间序列；
监测点i一段时间内监测到的机场噪声时间序列数据di可以表示为：
di＝{di1,di2,...,dim}
其中，m为监测的次数，di1为监测点i在某段时间内第1个时刻监测到的噪声数值，di2为监测点i在第2个时刻监测到的噪声数值，dim为监测点i在第m个时刻监测到的噪声数值。假设一共有n个监测点，那么所有监测点在这段时间内测量到机场噪声时间序列为：
其中：d1为监测点1的噪声时序，d2为监测点2的噪声时序，dn为监测点n的噪声时序，d11为监测点1在第1个时刻监测到的噪声数据，d12为监测点1在第2个时刻监测到的噪声数据，d1m为监测点1在第m个时刻监测到的噪声数据，d21为监测点2在第1个时刻监测到的噪声数据，d22为监测点2在第2个时刻监测到的噪声数据，d2m为监测点2在第m个时刻监测到的噪声数据，dn1为监测点n在第1个时刻监测到的噪声数据，dn2为监测点n在第2个时刻监测到的噪声数据，dnm为监测点n在第m个时刻监测到的噪声数据。
步骤2：对原始噪声时间序列进行预处理，创建标准的噪声时间序列数据集；
对监测到的原始机场噪声时间序列数据d归一化处理，采用方法为线性函数转化，即sij＝(dij-minvaluej)/(maxvaluej-minvaluej),其中，maxvaluej和minvaluej分别是j时刻所有测量值中的最大值和最小值，dij为监测点i在第j个时刻监测到的噪声数值，sij为归一化处理后的dij。从而得到一个标准的机场噪声时间序列数据集。
其中，s1是归一化后的d1，s2是归一化后的d2，sn是归一化后的dn，s是归一化后的d，即归一化的机场噪声时间序列数据集，s11为归一化后的d11，s12为归一化后的d12，s1m为归一化后的d1m，s21为归一化后的d21，s22为归一化后的d22，s2m为归一化后的d2m，sn1为归一化后的dn1，sn2为归一化后的dn2，snm为归一化后的dnm。
步骤3：用基于趋势分段的时间序列特征表示方法对每个监测点的噪声时序数据进行降维表示；
为了对噪声时间序列进行合理降维，保留序列中更有价值的趋势特征，本发明结合已经被成熟应用的迭代终点拟合算法对各监测点的噪声序列si进行降维，先将数据集中的每条噪声时序划分成多个具有不同趋势的分段。迭代终点拟合算法的目的是通过阈值参数ε对机场噪声时间序列进行分割，从而获得一个保留更多有用信息，去除冗余信息的噪声近似序列subi，并用subi中各分段的均值符号化表示
具体地，首先，选定并保留si中的起点si1和终点sim，用欧氏距离计算样本si中剩余的点sij到si1和sim连线
用最终被保留的ti+1个点顺序连成的ti个趋势分段来表示原序列si，那么si被降维表示为subi＝(subi1,subi2,...,
其中：sub1为s1的分段表示，sub2为s2的分段表示，subn为sn的分段表示，p11为s1第1分段的起始位置，p21为s2第1分段的起始位置，pn1为sn第1分段的起始位置，
步骤4：使用基于趋势分段的时间序列相似性度量方法，逐一度量各监测点的噪声时序与其他监测点的噪声时序之间的相似程度，并建立相似性矩阵；
获得监测点噪声时间序列新的表示后，使用基于趋势分段的时间序列相似性度量方法，对基于趋势分段表示下的监测点噪声时间序列进行相似性度量，并建立相应监测点与其他监测点之间的相似度矩阵。
假设要度量监测点i与监测点j的噪声时序数据的相似性，首先计算subi和subj中对应分段的符号化均值距离mdist，计算公式为：
其中，βmax和βmin分别指分段均值对应的高斯分布的划分区间断点的最大值和最小值，
其中，tik为第i个监测点噪声时序数据的第k分段的趋势变量，tjk为第j个监测点噪声时序数据的第k分段的趋势变量。将均值距离与趋势距离结合，定义了一种能度量不同长度序列的相似性度量函数：
其中，rest(subi)表示样本序列subi中除去第一个点subi1之后的剩余序列，rest(subj)表示序列subj中除去第一个点subj1之后的剩余序列。
使用tsdist度量每个监测点的噪声时序与其他监测点噪声时序之间的趋势分段距离，建立相应的相似度矩阵simi：
其中，tsdistij＝tsdist(subi,subj)表示第i个监测点与第j个监测点噪声时序之间的相似距离。
步骤5：根据相似性矩阵，找到每一个监测点的前k个相似性最高的的监测点，创建相似监测点集合；
按相似距离由小到大排列，寻找关联监测点集合中的前k个监测点作为相似监测点，并创建相应的关联监测点集合。
步骤6：度量各个监测点新的噪声时序与其相似监测点的新噪声时序之间的相似性，若相似性发生明显变化，则判定为异常。
逐一度量各个监测点新的噪声时序与其关联监测点的新的噪声时序之间的相似性，若相似性明显偏大，则判断该监测点出现异常，完成监测点异常识别。