sin(x)自相关函数

2024-10-17 12:33:46 发布

您现在位置:Python中文网/ 问答频道 /正文

我有一个时间序列数据,我想显示自相关函数。 (我们知道正弦函数的自相关是余弦函数)

我应用以下几种方法来实现这一点

import numpy as np
import matplotlib.pyplot as plt
from statsmodels.graphics.tsaplots import plot_acf, plot_pacf

x = np.arange(0,30,0.1)  #interval=0.1, 300 samples
y = np.sin(x)
y_cor = np.correlate(y,y,'full')  
lags = np.arange(-x[-1],x[-1]+0.1,0.1)

#sin(x)
plt.figure()
plt.plot(x,y)

#autocorrelation(numpy)
plt.figure()
plt.plot(lags,y_cor)
plt.xlabel('Lag')
plt.ylabel('autocorrelation')

#matplotlib
plt.figure()
plt.acorr(y,maxlags=y.size-1)

#statsmodels
plt.figure()
plot_acf(y,lags=y.size-1)

plt.show()

然而,结果是一个衰减的余弦函数,而不是一个纯的cos(x)。 我看到一些答案说,这是因为在计算自相关时,包在x区域外填充了零,但是如何修正这一点以得到纯cos(x)? enter image description hereenter image description hereenter image description hereenter image description here


Tags: 函数importnumpyplotmatplotlibasnpplt