DTM 相似度计算（matlab）代码记录

 DTM算法，一种衡量两个长度不同的时间序列的相似度的方法。

clc,clear
a=[8 9 1]'; % 输入两个列向量
b=[2 5 4 6]';
[rowa,la]=size(a); 
[rowb,lb]=size(b);
% 计算对应点的欧式距离d,a为列，b为行
for i=1:rowa
    for j=1:rowb
        d(i,j)=abs(a(i)-b(j));
    end
end
% 计算从左上到右下的累加距离D
D(1,1)=d(1,1);
for i=2:rowa     % 第一列累加
    D(i,1)=d(i,1)+D(i-1,1);
end
for i=2:rowb     % 第一行累加
    D(1,i)=d(1,i)+D(1,i-1);
end
for i=2:rowa     % 其他格子的累加
    for j=2:rowb
        D(i,j)=d(i,j)+min(D(i-1,j),min(D(i-1,j-1),D(i,j-1)));%双层的min结构能加快速度
    end
end
D(rowa,rowb)     % 显示两个列向量，消除大小规模后的最小距离