【Matlab实现】动态时间规划调整算法(DTW算法)——计算两个序

【Matlab实现】动态时间规划调整算法(DTW算法)——计算两个序列之间的相似度

【Matlab实现】动态时间规划调整算法(DTW算法)——计算两个序列之间的相似度

概述

算法原理与步骤

算法的实现

概述

DTW (Dynamic time warping)算法是可以度量两个独立时间序列的相似度的一种方法。曾被广泛应用在单词音频的匹配上。该方法主要用来解决在两段序列的时长不同的情况下，进行相似度的判断。

上图中，左侧时长相等，可以逐一进行欧式距离的计算，右侧则是时长不等，经过DTW之后得到的结果，可以看出来两个序列并不是一一对应的。

再比如上面左图，要得到蓝色序列与红色序列的相似度，因为可以看出来两个序列有经过平移的迹象，直接用一一匹配的方法显然是不合理的。要得到左图的对应效果，就需要用DTW方法。

算法原理与步骤

① 计算两个特征点之间的欧氏距离。构成一个 n*m 矩阵，距离矩阵。

②计算累计距离 得到DP矩阵。

计算后的的值，放到DP矩阵中，为了更加直观的理解，把这两个序列绘图如下：

其实在计算过程中，计算的顺序其实是有方向的。网上有很多的博客说的也非常清楚，博主在这里不再赘述。为了更好的理解计算过程，列举一个非常非常非常非常简单的例子来帮助理解，如下图： A B为带有两个特征值的序列，右边是其对应的DP矩阵的求解步骤。

③ 当计算完整个DP矩阵 后，右上角的值(不一定是右上角，就是最终的得到的那个矩阵角上的值)即为两个序列的累计距离。

④从右上角往左下角回溯，找到累计距离最短的路径，根据路径可以得到各个点之间的对应关系。

算法的实现

博主是利用matlab实现的这个算法，只是因为利用matlab可以很方便的查看矩阵和画图，检查算法的正确性，但是没有调用matlab中成形的函数，所以利用这个思路，用C/C++也是可以实现的，便于移植。

首先要写好两个函数。

一个是Get Min()；用来得到三个值中的最小值，在计算 DP矩阵 时用得到。

function min = GetMin(a,b,c)

if(a <= b && a <= c)

min = a;

elseif(b <= a && b <= c)

min = b;

elseif(c <= b && c <= a)

min = c;

end

end

另一个是GetMinIndex()；这是用来在得到 DTW 结果之后，方便显示特征点匹配的结果，返回两个序列对应特征点的索引。

function [index_i,index_j] = GetMinIndex(a,b,c,i,j)

%a 是相邻左上角，b 是相邻正上方，c说相邻正左方

%i 是当前的x坐标 j 是当前 y坐标

if(a <= b && a <= c)

index_i = i-1;

index_j = j-1;

elseif(b <= a && b <= c)

index_i = i-1;

index_j = j;

elseif(c <= b && c <= a)

index_i = i;

index_j = j-1;

end

end

接下来就是主函数了

%生成两个有明显平移性质的时间序列

x = zeros(1,50);

for i = 1:50

x(i) = sin(i*2*pi/50)+2;

end

y = zeros(1,50);

for i = 1:50

y(i) = sin(i*2*pi/50 + pi/6)+2;

end

x_len = length(x);

y_len = length(y);

plot(1:x_len,x);hold on

plot(1:y_len,y);hold on

%计算两序列每个特征点的距离矩阵

distance = zeros(x_len,y_len);

for i = 1:x_len

for j=1:y_len

distance(i,j) = (x(i)-y(j)).^2;

end

end

%计算两个序列

DP = zeros(x_len,y_len);

DP(1,1) = distance(1,1);

for i=2:x_len

DP(i,1) = distance(i,1)+DP(i-1,1);

end

for j=2:y_len

DP(1,j) = distance(1,j)+DP(1,j-1);

end

for i=2:x_len

for j=2:y_len

DP(i,j) = distance(i,j) + GetMin(DP(i-1,j),DP(i,j-1),DP(i-1,j-1));

end

end

%回溯，找到各个特征点之间的匹配关系

i = x_len;

j = y_len;

while(~((i == 1)&&(j==1)))

plot([i,j],[x(i),y(j)],'b');hold on %画出匹配之后的特征点之间的匹配关系

if(i==1)

index_i = 1

index_j = j-1

elseif(j==1)

index_i = i-1

index_j = 1

else

[index_i,index_j] = GetMinIndex(DP(i-1,j-1),DP(i-1,j),DP(i,j-1),i,j)

end

i = index_i;

j = index_j;

end

最终效果如下图，可以看出来是考虑了平移之后的匹配。

【Matlab实现】动态时间规划调整算法(DTW算法)——计算两个序相关教程

数据库实现-索引文件(Index File)BPTree 版本

数据库实现-索引文件(Index File)BPTree 版本 Project 1 - Index file 这个项目呢，是尝试复刻一下数据库中的索引系统，只实现了内存版本的BPTree，索引文件并未建立成功(本博客持续更新ing) 文章目录 Project 1 - Index file 1.Abstract 2.Tasks Require

Cube实现天空盒子，看了你自己也能做

Cube实现天空盒子，看了你自己也能做 大家可能都了解了引擎自带的天空盒子，菜鸟也用过，感觉很棒。但是菜鸟在使用的过程中感觉天空盒子太大了，而场景比较小，所以想缩放一下天空盒子，奈何没找见，所以私下自己使用Cube实现了一下，分享给大家 1.坐标空间

02-链表

02-链表 链表 1.定义 链表是一种线性的动态数据结构，与动态数组，栈，队列的“动态”不同(这三种数据结构的底层都是静态数组，依托resize解决容量问题)，链表是真正意义上的动态数据结构 链表与数组相比，失去了随机访问的能力，但却拥有动态的优点(不需

涉及简历

涉及简历 反射 反射概念：反射就是动态调用一个类/一个对象的成员变量/属性、方法、构造器的功能 定义测试的User：看到原类的构造器、成员属性、方法，想着怎么使用反射生成 package reflect;public class User { private int id=1; private String name=张三

ATM机项目实现：实现登录，注册操作(登录后可以存款，取款，查

ATM机项目实现：实现登录，注册操作(登录后可以存款，取款，查询余额，转账) 提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档 文章目录 前言 一、分析需求 二、建立TestDemo.java 三、建立ATM.java 四、建立Bank.java 五、建立BankCard.

GPU中实现反距离加权插值(IDW)

GPU中实现反距离加权插值(IDW) 利用GPU实现IDW(反距离加权插值) IDW的实现比较简单，已知插值点位比较少的情况下，可以直接遍历所有插值点，来获取临近的几个点，进行插值运算。插值点较多时，需要可以使用 kd-tree 来加速临近点的查找。本次仅对 小数据量 情

前端实现一个小说分页的功能

前端实现一个小说分页的功能 来源 | https://juejin.im/post/6886418644381728776 先让我找找我的思路在哪里？ 在小说读书APP中，都会有分页的功能，那么前端如何实现这个功能呢？ 因为没有什么思路，那就只能在前辈的项目中寻找思路了。 这不，直接打开起点

centos7 docker容器内实现动态maven打包

centos7 docker容器内实现动态maven打包 背景：业务需要动态更新maven依赖，目前想到的办法是动态打包执行jar实现。已有centos7镜像且已安装java sdk。 linux下安装maven环境 官网下载安装包：http://maven.apache.org/download.cgi cd /mkdir optcd opt# 下