DTW MATLAB代码

以下是一个简单的 MATLAB DTW 代码示例：

function [D, w] = dtw(x, y)

% x: 第一个序列
% y: 第二个序列
% D: DTW距离矩阵
% w: 最佳路径

% 计算两个序列的长度
Nx = length(x);
Ny = length(y);

% 初始化DTW距离矩阵
D = inf(Nx+1,Ny+1);
D(1,1) = 0;

% 计算每个格子的距离
for i = 1:Nx
    for j = 1:Ny
        dist = (x(i)-y(j))^2;
        D(i+1,j+1) = dist + min([D(i,j+1), D(i+1,j), D(i,j)]);
    end
end

% 获取最佳路径
w = [];
i = Nx+1;
j = Ny+1;
while i > 1 || j > 1
    w = [i-1, j-1; w]; %#ok<AGROW>
    if i == 1
        j = j - 1;
    elseif j == 1
        i = i - 1;
    else
        [~, idx] = min([D(i-1,j), D(i,j-1), D(i-1,j-1)]);
        if idx == 1
            i = i - 1;
        elseif idx == 2
            j = j - 1;
        else
            i = i - 1;
            j = j - 1;
        end
    end
end
w = [1, 1; w];
end

使用该函数，你可以计算两个序列之间的 DTW 距离。例如：

x = [1, 3, 4, 2, 1];
y = [1, 1, 3, 2, 0, 1];
[D, w] = dtw(x, y);
disp(D(end,end)) % 输出 DTW 距离
plot(x, 'b'); hold on;
plot(y, 'r'); hold on;
plot(w(:,1), w(:,2), 'g'); hold off;

这将输出 DTW 距离并绘制最佳路径。