DTW matlab代码

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

function [dist,D] = dtw(x,y,w)
% x: NxD feature matrix, N是序列长度，D是特征维度
% y: MxD feature matrix, M是序列长度，D是特征维度
% w: DTW窗口大小
% dist: DTW距离
% D: D(i,j)是x(i,:)和y(j,:)之间的欧几里得距离

N = size(x,1);
M = size(y,1);
if nargin < 3
    w = max([N,M]);
end

% 计算点之间的距离矩阵
D = zeros(N,M);
for i = 1:N
    for j = 1:M
        D(i,j) = norm(x(i,:)-y(j,:));
    end
end

% 初始化距离矩阵
C = inf(N+1,M+1);
C(1,1) = 0;

% 动态规划计算最短路径
for i = 2:N+1
    for j = max(2,i-w):min(M+1,i+w)
        C(i,j) = D(i-1,j-1) + min([C(i-1,j),C(i,j-1),C(i-1,j-1)]);
    end
end

% 返回DTW距离
dist = C(N+1,M+1);
end

这个函数可以计算两个时间序列之间的DTW距离。在这个函数中，我们首先计算了输入时间序列之间的所有点之间的欧几里得距离，然后使用动态规划计算了最短路径，并返回最终DTW距离。我们还可以通过设置窗口大小w来限制DTW算法的搜索范围。