一种中心时间序列的动态求解方法
【专利摘要】本发明公开了一种中心时间序列的动态求解方法，包括以下步骤：标注时间序列；计算动态匹配；输出中心时间序列。本发明给出了一种新的中心时间序列动态匹配距离，保持了动态形状特征的相似性，比欧几里德平均距离法求中心时间序列的方法更好地呈现形态特征。本发明所确定的一种新的中心时间序列动态匹配累计方式，保证了中心时间序列到相关的时间序列动态匹配距离最小，而且时间复杂度为O(n(p+q)/2)3)，比动态时间弯曲距离法求解中心时间序列的方法的时间复杂度好两个数量级。本发明给出的所有预设中点中获取最小误差的方法，对于时间序列聚类而言，使时间序列聚类更准确。
【专利说明】—种中心时间序列的动态求解方法
【技术领域】
[0001]本发明涉及一种数据处理技术，特别是一种中心时间序列的动态求解方法。
【背景技术】
[0002]时间序列是指按时间顺序排列的一组数据，是一类重要的复杂数据对象。它广泛存在于各种不同领域中，如股票价格、各种汇率、销售数量、产品的生产能力、天气数据、网购行为和观点评述等。大量时间序列数据真实地记录了各个时刻的所有重要信息，依时间呈现动态性，若能得到有效地分析，发现其中各时间序列之间的相互关系，必将大大提高人们的认识和理解，进而有效地预测和控制，在大数据时代显得尤为重要。人们发现这些时间序列的中心时间序列直接反映了它们的通用模式，有效把握和分析时间序列，至关重要的就是要有效地求解这些中心时间序列，但不幸的是，在人们致力于研发时间序列分析和挖掘系统工具中，中心时间序列求解作为其核心程序，一直是个瓶颈。
[0003]观察时间序列数据，可以看作是具有一定动态特性的系统的外在表现形式，这种动态特性主要表现在于它因时间变化而形成的全局形态特征和局部形态特征，其中全局形态特征描述了时间序列的起伏变化等，局部形态特征则表现为时间序列局部时间点上的异常观测值，如不连续点、极值点、突变点和转折点等。在求解两个时间序列的中心时间序列过程中，中心时间序列必须具有与参考时间序列相似的动态特性，保持相似性。然而，由于需要保持其特定的形状特征，使得目前常用的一些相似性度量和聚类方法失去了原有的优越性。
[0004]求解中心时间序列，当前主要采用欧几里德距离法和动态时间弯曲距离法。欧几里德平均距离法求中心时间序列，是直接通过时间序列间对应的两个节点求欧几里德距离上的平均中值，计算复杂度为O (η)，但这种欧几里德二值平均会直接抵消原时间序列在时间弯曲上的形态特征，导致所求的中心时间序列不支持原时间序列依形态特征。动态时间弯曲距离法求解中心时间序列的方法中，先通过随机产生初始中心时间序列，然后通过迭代产生由每个单独的时间序列和暂时的中心时间序列之间的动态时间弯曲距离，逐步精炼，并生产新的临时中心时间序列，直到迭代到临时中心时间序列是稳定的，与欧几里德平均距离法求到的中心时间序列相比，在时间弯曲上的表现较好，但计算复杂度为0(η3)，而且因为不是通过最小匹配而是通过迭代逼近法求中心时间序列，所以不易保持误差最小化。
【发明内容】
[0005]为解决现有技术存在的上述问题，本发明要提出一种既要在时间弯曲上的表现较好、又要降低时间复杂度的中心时间序列的动态求解方法。
[0006]为了实现上述目的，本发明的技术方案如下:一种中心时间序列的动态求解方法，包括以下步骤:
[0007]Α、标注时间序列:对于多个时间序列，先完成每两个时间序列之间的求解，最终求出多个时间序列的中心时间序列。不失一般性，假定现有两个时间序列x[l:m] = (xi，x2，…，xm) (m>l), y[l:n] = (y1, y2,…，yn) (η>1),需要求出它们的中心时间序列(^[Ι:!:] =、，c2,…，cq) (q>l)。
[0008]B、计算动态匹配:时间序列X与c之间存在动态时间匹配关系W=Ewm, Wl，2]T.[w2；1, w2；2]T,..., [wk；1, wk；2]T,..., [wL；1, wL；2] ], k ∈ {1，2，...，L}，wk’i e {1,2，...，m}，wk，2e {1,2, -,q}0动态匹配距离为:
【权利要求】
1.一种中心时间序列的动态求解方法，其特征在于:包括以下步骤: A、标注时间序列:对于多个时间序列，先完成每两个时间序列之间的求解，最终求出多个时间序列的中心时间序列；不失一般性，假定现有两个时间序列x[l:m] = (xi，x2，…，xm)(m>l), y[l:n] = (y1, y2,…，yn) (η>1),需要求出它们的中心时间序列(^[Ι:!:] =、，C2,…，cq)(q>l); B、计算动态匹配:时间序列X与c之间存在动态时间匹配关系T.[w2；1, w2；2]T,..., [wk；1, wk；2]T,..., [wL；1, wL；2] ], k e {1，2，...，L}，wk’i e {1,2，...，m}，wk，2 e {I, 2，…，q};动态匹配距离为:
【文档编号】G06F17/30GK103942300SQ201410151135
【公开日】2014年7月23日 申请日期:2014年4月15日 优先权日:2014年4月15日
【发明者】刘洪波, 孙焘, 周亮, 孙野青 申请人:大连海事大学