python写出DTW算法

以下是Python代码实现DTW算法：

import numpy as np

def dtw(s, t):
    n = len(s)
    m = len(t)
    dtw_matrix = np.zeros((n+1, m+1))
    for i in range(n+1):
        for j in range(m+1):
            dtw_matrix[i, j] = np.inf
    dtw_matrix[0, 0] = 0
    
    for i in range(1, n+1):
        for j in range(1, m+1):
            cost = abs(s[i-1] - t[j-1])
            dtw_matrix[i, j] = cost + min(dtw_matrix[i-1, j], dtw_matrix[i, j-1], dtw_matrix[i-1, j-1])
    
    return dtw_matrix[n, m]

其中，`s`和`t`分别代表两个时间序列，`dtw_matrix`是DTW矩阵，用于存储两个序列之间的距离。在代码实现中，首先创建一个全0的矩阵，然后将第一行和第一列的值设为无穷大，因为第一行和第一列的值代表其中一个序列为空序列的情况，所以它们与任何一个非空序列的距离都是无穷大。接着，通过双重循环计算DTW矩阵中每个位置的值，最后返回矩阵的右下角元素，即为两个序列的DTW距离。