python - Python - SciPy：从一组周期性点进行预测

Question

我正在学习 SciPy，如果找到正确的函数，我应该很容易解决一个问题：我有一组涵盖 10 个周期的点（2D）。我想预测下一个时期的价值。

我想我必须创建一个与我的点相对应的周期函数，然后在这个模型上取点，但我不知道该怎么做！

你可以帮帮我吗 ？

提前致谢

Answer 1

外推绝非易事。它几乎总是很差，除非你对数据有一些强有力的假设。

在您的情况下，您可以尝试一下，但我认为没有立即可用的东西。

我会尝试：

• 确定自相关的第一个最大值
• 通过将此值的倍数移位来扩展您的信号

如果需要，之后进行插值。

例子：

import numpy as np
import matplotlib.pyplot as plt

def autocorr(x):
    result = np.correlate(x, x, mode='full')
    return result[result.size/2:]

data = np.sin(np.linspace(0,30,300)) + np.random.random((300)) * 0.1

plt.subplot(3,1,1)
#plt.plot(data,"b-")
plt.plot(data,"bx")

acorr = autocorr(data)
acorr_diff = np.diff(acorr)

maxima = [i+1 for i in range(acorr_diff.shape[0]-1) 
          if acorr_diff[i]>=0 and acorr_diff[i+1]<0]


plt.subplot(3,1,2)
plt.plot(acorr)
for m in maxima:
    plt.axvline(m, color="b", alpha=0.5)

first_max = maxima[0]
new_data = np.hstack([data[:4*first_max],data])

plt.subplot(3,1,3)
plt.plot(data)
#plt.plot(data,"b-", alpha=0.1)
plt.plot(data,"bx")
plt.plot(new_data,"r-")
#plt.plot(new_data,"rx")

plt.show()

这只是一个非常基本的实现。当然，它有局限性，但原则应该是明确的。