如何在Matlab中识别时间序列数据中的重复模式

Question

我正在使用Matlab计算ENSO指数，一个条件是我必须找到异常的海面温度。条件是，厄尔尼诺事件的特征是海面温度在5个月内高于标准化“0值”0.5度。我已经尽可能使我的月度时间序列数据具有逻辑性（即，“1”是高于0.5的月度数据值，“0”是低于0.5的月度数据值），但我想知道Matlab中是否有命令允许我识别该值何时重复5次或更多。

作为示例代码：

Monthly_data=[0 0 1 1 1 1 1 0 0 0 1 1 0 0 0 0 1 0 1 1 1 1 1 1 1 0]

理想情况下，我需要一个命令，该命令可以发现每个命令后面至少出现五个“1”。这是否存在？

如果需要更多信息，请告诉我，我是matlab的新手，所以我还不确定在这里提问的结构和语法。

非常感谢。

Answer 1

不确定这是你需要的，但可能会给你一些指导。

> x = diff(Monthly_data);
> find(x==-1)-find(x==1)
ans =

   5   2   1   7

这些是序列的长度。您可能需要用0填充数组的前端和末端，以消除缺少一个边界的序列。

要查找大于5的序列的开始索引：

> s=find(x==1);
> s(find(x==-1)-s>5)
ans =  18

或

> s(find(x==-1)-s>=5)
ans =

    2   18

请注意，由于差异滞后，这些比数组索引多一个，或者将其视为从零开始索引的位置。