LCSS算法 matlab

以下是 MATLAB 实现的 LCSS 算法示例代码：

function [lcss, s1_match_idx, s2_match_idx] = LCSS(s1, s2, epsilon, window_size)
    % 计算两个序列 s1 和 s2 的 LCSS
    % epsilon：允许的误差
    % window_size：滑动窗口大小

    % 初始化 lcss 矩阵
    lcss = zeros(length(s1)+1, length(s2)+1);

    % 遍历 s1 和 s2，计算 lcss 矩阵
    for i = 2:length(s1)+1
        for j = 2:length(s2)+1
            % 如果当前位置上的元素相等，则 lcss(i,j) = lcss(i-1, j-1) + 1
            if abs(s1(i-1) - s2(j-1)) <= epsilon
                lcss(i, j) = lcss(i-1, j-1) + 1;
            % 如果当前位置上的元素不相等，则 lcss(i,j) = max(lcss(i-1,j), lcss(i,j-1))
            else
                lcss(i, j) = max(lcss(i-1, j), lcss(i, j-1));
            end
        end
    end

    % 计算匹配的索引
    s1_match_idx = [];
    s2_match_idx = [];
    i = length(s1) + 1;
    j = length(s2) + 1;
    while i > 1 && j > 1
        if abs(s1(i-1) - s2(j-1)) <= epsilon
            s1_match_idx = [i-1, s1_match_idx];
            s2_match_idx = [j-1, s2_match_idx];
            i = i - 1;
            j = j - 1;
        else
            if lcss(i-1, j) > lcss(i, j-1)
                i = i - 1;
            else
                j = j - 1;
            end
        end
    end

    % 如果滑动窗口大小不为 0，则对匹配的索引进行滑动窗口处理
    if window_size > 0
        s1_match_idx = slide_window(s1_match_idx, window_size);
        s2_match_idx = slide_window(s2_match_idx, window_size);
    end
end

function [result] = slide_window(arr, window_size)
    % 对数组 arr 进行滑动窗口处理，窗口大小为 window_size
    result = [];
    i = 1;
    while i <= length(arr)
        j = i + 1;
        while j <= length(arr) && arr(j) - arr(i) < window_size
            j = j + 1;
        end
        result = [result, round(mean(arr(i:j-1)))];
        i = j;
    end
end

该算法的输入为两个序列 `s1` 和 `s2`，以及两个参数 `epsilon` 和 `window_size`。其中 `epsilon` 表示允许的误差，`window_size` 表示滑动窗口大小（如果为 0，则不进行滑动窗口处理）。输出为 LCSS 长度、s1 和 s2 的匹配索引。