本文介绍一种方法,帮助我们了解一个时间序列是否可以预测，或者说了解可预测能力有多强。

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。