本文介绍一种方法,帮助我们了解一个时间序列是否可以预测，或者说了解可预测能力有多强。
Sample Entropy （样本熵）
Sample Entropy是Approximate Entropy(近似熵)的改进，用于评价波形前后部分之间的混乱程度， 熵越大，乱七八糟的波动越多，越不适合预测；熵越小，乱七八糟的波动越小，预测能力越强。
具体思想和实现如下：
- 思想 Sample Entropy最终得到一个 -np.log(A/B) ，该值越小预测难度越小，所以A/B越大，预测难度越小。 A：从0位置开始，取m 1个元素构成一个向量，然后移动一步，再取m 1个元素构成一个向量，如此继续直到最后得到一个向量集合Xa,看有多少向量彼此的距离小于容忍度r(即有多少向量彼此相似，又称自相似个数)。 B：从0位置开始，取m个元素构成一个模板向量，然后移动一步，再取m个元素构成一个模板向量，如此继续直到最后得到一个向量集合Xb，看有多少向量彼此的距离小于容忍度r(即有多少向量彼此相似，又称自相似个数)。 而实际上A总是小于等于B的，所以A/B越接近1，预测难度越小，直觉上理解，应该就是波形前后部分之间的变化不大，那么整个时间序列的波动相对来说会比较纯(这也是熵的含义，熵越小，信息越纯，熵越大，信息越混乱)，或者说会具有一定的规律，而如果A和B相差很大，则时间序列波动不纯，或者说几乎没有规律可言。 比如：U = [0.2, 0.6, 0.7, 1.2, 55, 66]，m=2, 那么可以计算得到： Xa = [[0.2, 0.6, 0.7], [0.6, 0.7, 1.2], [0.7, 1.2, 55.0], [1.2, 55.0, 66.0]] Xb = [[0.2, 0.6], [0.6, 0.7], [0.7, 1.2], [1.2, 55.0], [55.0, 66.0]] 假设Xa中相似向量的个数比Xb多，那么应该出现Xa满足r，但是Xb不满足r的情况，但是拿Xa和Xb的前两个向量来分析，如果Xa满足r，则0.2-0.6 ，0.6-0.7，0.7-1.2中的最大值应该<=r,也就是说0.2-0.6 ，0.6-0.7肯定<=r,如此推断，Xb肯定也满足r， 所以只有可能出现Xb满足r，Xa不满足r的情况。
- python实现
def SampEn(U, m, r):
"""
用于量化时间序列的可预测性
:param U: 时间序列
:param m: 模板向量维数
:param r: 距离容忍度，一般取0.1~0.25倍的时间序列标准差，也可以理解为相似度的度量阈值
:return: 返回一个-np.log(A/B)，该值越小预测难度越小
"""
def _maxdist(x_i, x_j):
"""
Chebyshev distance
:param x_i:
:param x_j:
:return:
"""
return max([abs(ua - va) for ua, va in zip(x_i, x_j)])
def _phi(m):
x = [[U[j] for j in range(i, i m - 1 1)] for i in range(N - m 1)]
C = [len([1 for j in range(len(x)) if i != j and _maxdist(x[i], x[j]) <= r]) for i in range(len(x))]
return sum(C)
N = len(U)
return -np.log(_phi(m 1) / _phi(m))
if __name__ == '__main__':
_U = [0.2, 0.6, 0.7, 1.2, 55, 66]
rand_small = np.random.randint(0, 100, size=120)
rand_big = np.random.randint(0, 100, size=136)
m = 2
print(SampEn(_U, m, r=0.2 * np.std(_U)))
print(SampEn(rand_small, m, r=0.2 * np.std(rand_small)))
print(SampEn(rand_big, m, r=0.2 * np.std(rand_big))) ok,本篇就这么多内容啦~，感谢阅读O(∩_∩)O。
