matlab写dtw算法

下面是MATLAB代码实现DTW算法的示例：

function [dist, D] = dtw(s, t, w)
% DTW算法的MATLAB实现
% s：序列1
% t：序列2
% w：约束窗口大小（默认为无限制）

if nargin < 3
    w = inf;
end

% 初始化距离矩阵
[n, m] = size(s, t);
D = zeros(n+1, m+1) + inf;
D(1, 1) = 0;

% 计算动态时间规整距离
for i = 2:n+1
    for j = max(2, i-w):min(m+1, i+w)
        cost = norm(s(i-1, :) - t(j-1, :));
        D(i, j) = cost + min([D(i-1, j), D(i, j-1), D(i-1, j-1)]);
    end
end

% 计算DTW距离
dist = D(n+1, m+1);

其中：

- s和t是两个数据序列
- w是约束窗口的大小，用于限制两个序列在时间轴上的扭曲程度。如果w等于无限制，默认为无限制。
- dist是DTW距离
- D是距离矩阵，其中D(i,j)表示在序列1的第i个位置和序列2的第j个位置匹配的最小距离。

使用示例：

% 生成两个长度为10的随机序列
s = randn(10, 2);
t = randn(10, 2);

% 计算DTW距离
dist = dtw(s, t);