[发明专利]一种基于最小平均距离的动态时间规整方法

说明书

技术领域

本发明涉及一种动态时间规整方法。

背景技术

距离是度量两段时间序列相异(似)性的常用方法。由于序列值在时间轴上分布会有平移、伸缩，从而导致直接的点对点计算会有偏离常识性判断的风险。例如，对于两个序列(a,a,b,b,a)和(a,b,b,a,a)模式，在实际应用或工程实践中需要判断为一致，但使用距离计算其数值为2|a-b|。使用传统的动态时间规整算法计算距离所得值为0，这是比较符合实际工程的需要的。但是传统动态时间规整算法也有弊端，例如，考虑使用动态时间规整距离计算以下三段直流信号的匹配情况：

A＝(0.1,0.1)

B＝(0.25,0.25,0.25,0.25)

C＝(0.28,0.28,0.28)

信号A与B的动态时间规整距离为0.6，动态时间规整路径的长度为4；而信号C与A的动态时间规整距离0.54，动态时间规整路径的长度为3。从我们直觉判断信号B与A比信号C与A更接近，但从动态时间规整距离决定的相异性来看，反而是C比B更接近A。造成这个反直觉的效果是我们仅仅考虑了距离，而忽视了路径长度，假如以动态时间规整距离相对路径长度的平均值而言，B比C更接近A。

发明内容

为了克服已有技术动态时间规整方式的误识率较高的不足,本发明提供了一种误识率较低的基于最小平均距离的动态时间规整方法。

本发明解决其技术问题所采用的技术方案是：

一种基于最小平均距离的动态时间规整方法，包括以下步骤：

1)对于两段时间序列

S＝(s₁,s₂,...,s_n)

T＝(t₁,t₂,...,t_m)

创建平均距离矩阵其中是点(1,1)到点(i,j)的最小平均距离，是起点(1,1)到终点(n,m)的平均距离，也是时间序列S和T的基于最小平均距离的最优匹配距离。从起点(1,1)到终点(n,m)的最优匹配路径P＝(p₁,p₂,......,p_l)，其中p_i＝(q_i,r_i)，而1＜＝q_i＜＝n,1＜＝r_i＜＝m，与传统动态时间规整算法一样，满足如下的条件：

单调性：

q_i-1≤q_i且r_i-1≤r_i

连续性：

q_i-1＝q_i或q_i-1+1＝q_i且r_i-1＝r_i或r_i-1+1＝r_i

边界条件：

p₁＝(1,1)p_l＝(n,m)

再创建路径长度矩阵l＝(l(i,j))_n×m，其中l(i,j)点(1,1)到点(i,j)的最小平均距离的路径长度，自然而l(n,m)＝l为路径长度矩阵满足的单调性、连续性、和边界条件且为S和T的匹配距离(路径长度平均距离)，它是所有的从起点(1,1)到终点(n,m)的可能的路径里平均距离最小的；这样的路径长度l满足max(n,m)≤l≤n+m-1，任意的路径数目的一个下界值为：

num(P)＞2·3^min(m，n)-3；

2)求的动态规划的过程如下：

2.1)初始化

2.2)对平均距离矩阵进行首行首列的赋值，首行取与其同行前列的近邻进行加权组合计算算术平均，首列取与其同列前行的近邻进行加权组合计算算术平均；

2.3)对平均距离矩阵进行非首行、非首列的赋值，先将当前点取与其前行前列的三邻点分别进行组合并计算算术平均，在得到的三个平均值中选取最小的那个，作为当前点的平均距离值；并将当前点的路径长度值在被选中的那个邻点的长度值基础上增加1；

2.4)输出和l(n,m)。