我在MATLAB中进行后续时间序列匹配时遇到了困难(我刚开始使用MATLAB)。
我有两个时间序列:A(长度为a)和B(长度为b)。假设a远大于b。任务是根据欧几里得度量找到A中与B最接近的窗口。
为了实现这一点,我构建了一个额外的矩阵C,用于存储A中长度为b的所有子序列,然后使用pdist2(C, B)。显然,这种方法运行速度慢且需要大量内存。
因此,我有几个问题:
-
如何在不使用循环的情况下获得C(实际上是重塑A)?
-
解决这个问题的常见方法有哪些?(最好是MATLAB中的方法,但其他环境的方法也可以)
感谢您的帮助!
回答:
对于第一个问题,您可以尝试
tmp = repmat(A,1,b);C = reshape([tmp zeros(1,b)],a,b);C = C(1:(a-b+1),:);
此外,与这个非常好的解决方案相比,pdist2
的速度非常慢:在MATLAB中高效计算成对平方欧几里得距离