一种中心时间序列的动态求解方法技术

本发明专利技术公开了一种中心时间序列的动态求解方法，包括以下步骤：标注时间序列；计算动态匹配；输出中心时间序列。本发明专利技术给出了一种新的中心时间序列动态匹配距离，保持了动态形状特征的相似性，比欧几里德平均距离法求中心时间序列的方法更好地呈现形态特征。本发明专利技术所确定的一种新的中心时间序列动态匹配累计方式，保证了中心时间序列到相关的时间序列动态匹配距离最小，而且时间复杂度为O(n(p+q)/2)3)，比动态时间弯曲距离法求解中心时间序列的方法的时间复杂度好两个数量级。本发明专利技术给出的所有预设中点中获取最小误差的方法，对于时间序列聚类而言，使时间序列聚类更准确。

【技术实现步骤摘要】

本专利技术涉及一种数据处理技术，特别是一种中心时间序列的动态求解方法。
技术介绍
时间序列是指按时间顺序排列的一组数据，是一类重要的复杂数据对象。它广泛存在于各种不同领域中，如股票价格、各种汇率、销售数量、产品的生产能力、天气数据、网购行为和观点评述等。大量时间序列数据真实地记录了各个时刻的所有重要信息，依时间呈现动态性，若能得到有效地分析，发现其中各时间序列之间的相互关系，必将大大提高人们的认识和理解，进而有效地预测和控制，在大数据时代显得尤为重要。人们发现这些时间序列的中心时间序列直接反映了它们的通用模式，有效把握和分析时间序列，至关重要的就是要有效地求解这些中心时间序列，但不幸的是，在人们致力于研发时间序列分析和挖掘系统工具中，中心时间序列求解作为其核心程序，一直是个瓶颈。观察时间序列数据，可以看作是具有一定动态特性的系统的外在表现形式，这种动态特性主要表现在于它因时间变化而形成的全局形态特征和局部形态特征，其中全局形态特征描述了时间序列的起伏变化等，局部形态特征则表现为时间序列局部时间点上的异常观测值，如不连续点、极值点、突变点和转折点等。在求解两个时间序列的中心时间序列过程中，中心时间序列必须具有与参考时间序列相似的动态特性，保持相似性。然而，由于需要保持其特定的形状特征,使得目前常用的一些相似性度量和聚类方法失去了原有的优越性。求解中心时间序列，当前主要采用欧几里德距离法和动态时间弯曲距离法。欧几里德平均距离法求中心时间序列，是直接通过时间序列间对应的两个节点求欧几里德距离上的平均中值，计算复杂度为O(n)，但这种欧几里德二值平均会直接抵消原时间序列在时间弯曲上的形态特征，导致所求的中心时间序列不支持原时间序列依形态特征。动态时间弯曲距离法求解中心时间序列的方法中，先通过随机产生初始中心时间序列，然后通过迭代产生由每个单独的时间序列和暂时的中心时间序列之间的动态时间弯曲距离，逐步精炼，并生产新的临时中心时间序列，直到迭代到临时中心时间序列是稳定的，与欧几里德平均距离法求到的中心时间序列相比，在时间弯曲上的表现较好，但计算复杂度为O(n3)，而且因为不是通过最小匹配而是通过迭代逼近法求中心时间序列，所以不易保持误差最小化。
技术实现思路
为解决现有技术存在的上述问题，本专利技术要提出一种既要在时间弯曲上的表现较好、又要降低时间复杂度的中心时间序列的动态求解方法。为了实现上述目的，本专利技术的技术方案如下：一种中心时间序列的动态求解方法，包括以下步骤：A、标注时间序列：对于多个时间序列，先完成每两个时间序列之间的求解，最终求出多个时间序列的中心时间序列。不失一般性，假定现有两个时间序列x[1:m]=(x1,x2,…,xm)(m>1),y[1:n]=(y1,y2,…,yn)(n>1)，需要求出它们的中心时间序列c[1:t]=(c1,c2,…,cq)(q>1)。B、计算动态匹配：时间序列x与c之间存在动态时间匹配关系W=[[w1,1,w1,2]T,[w2,1,w2,2]T,…,[wk,1,wk,2]T,…,[wL,1,wL,2]],k∈{1,2,…,L本文档来自技高网...

【技术保护点】
一种中心时间序列的动态求解方法，其特征在于：包括以下步骤：A、标注时间序列：对于多个时间序列，先完成每两个时间序列之间的求解，最终求出多个时间序列的中心时间序列；不失一般性，假定现有两个时间序列x[1:m]=(x1,x2,…,xm)(m>1),y[1:n]=(y1,y2,…,yn)(n>1)，需要求出它们的中心时间序列c[1:t]=(c1,c2,…,cq)(q>1)；B、计算动态匹配：时间序列x与c之间存在动态时间匹配关系W=[[w1,1,w1,2]T,[w2,1,w2,2]T,…,[wk,1,wk,2]T,…,[wL,1,wL,2]],k∈{1,2,…,L},wk,1∈{1,2,…,m},wk,2∈{1,2,…,q}；动态匹配距离为：ρ(W,x,y)=Σk=1L(xwk,1-ywk,2)2---(1)]]>因为时间序列间存在形状变化特征，动态时间匹配关系并不是唯一的，而是存在动态时间匹配簇(H>1)，其中，最小动态匹配距离为：在路径上的中心时间序列候选集t中，最小动态时间匹配上的中心时间序列为：C、输出中心时间序列c有以下7个子步骤：C1、输入时间序列x,y，读出x序列的长度赋给m，读出y的长度赋给n；C2、初始化伴随矩阵M[(m+1)*(n+1)]，是所有行的值为0,而误差error值设为无穷大；C3、读取x序列，检索伴随矩阵M，依公式(2)计算规整映射路径到x序列；C4、读取y序列，检索伴随矩阵M，规整映射路径到y序列；C5、从节点i到1倒序累计预设中心点到x序列的距离和，以及从节点j到1倒序累计预设中心点到y序列的距离和，依公式(4)计算预设第k个中心点距离Disti，j，k(x，y)：Disti，j，k(x，y)=DT(c[1:k],x[1:i])+DT(c[1:k],y[1:j])    (4)这里c＝Centi，j，k(x，y);Centi，j，k[k]=avg(x[p:i],y[q:j]);C6、依公式(5)更新伴随矩阵值mvalue：mvalue=Distp，q，k‑1+err(x[p:i],y[q:j])      (5)这里err(x,y)依公式(6)在所有预设中点坐标λ中获取最小误差：err(x,y)=minλ(Σi=1m(xi-λ)2+Σj=1n(yj-λ)2)---(6)]]>C7、依最小误差输出中心时间序列c。...

【技术特征摘要】
1.一种中心时间序列的动态求解方法，其特征在于：包括以下步骤：
A、标注时间序列：对于多个时间序列，先完成每两个时间序列之间的求解，
最终求出多个时间序列的中心时间序列；不失一般性，假定现有两个时间序列
x[1:m]=(x1,x2,…,xm)(m>1),y[1:n]=(y1,y2,…,yn)...

【专利技术属性】
技术研发人员：刘洪波，孙焘，周亮，孙野青，
申请(专利权)人：大连海事大学，
类型：发明
国别省市：辽宁;21