利用matlab实现dtw

function [distpath] = dtw（abr）
%% This function is to compute the two time series distance using DTW
% input:  a-- a time serise
%         b-- a time serise
%         r-- warping window
% output: dist-- the distance of a and b
%         path-- the warping path matrix which align points

%% initializing
a                    = zscore（a）;
b                    = zscore（b）;
n                    = length（a）;
m                    = length（b）;
r                    = ceil（min（[nm]）*r）;
Dist（1:n1:m）        = inf;

%% construct the distance matrix
for i = 1 : n
    for j = 1 : m %max（[1i-r]）:min（[mi+r]）%这里是为了减少计算，但是我做的实验几乎都是要使用r=1的。（并且发现一个问题，当两条时间序列长度差别特别大的时候，这里的r限制不起作用，计算都得inf，说明会导致过度弯曲）
%             Dist（ij）=sqrt（（a（i）-b（j））^2）; 
        Dist（ij） = （a（i）-b（j））^2;
    end
end

%% construct the accumulative matrix
for i = 1 : n
    for j = max（[1i-r]） : min（[mi+r]）
        if i == 1 && j == 1
            continue;
        end
        if i == 1   %第一行的累积距离