dd/PA/code/stationary season.py
2025-05-15 16:32:52 +04:00

12 lines
446 B
Python

# Ñîçäàåì ñåçîííûé ðÿä
seasonal_ts = 5 * np.sin(np.linspace(0, 10*np.pi, n)) + np.linspace(0, 5, n) + np.random.normal(size=n)
# Ñåçîííîå äèôôåðåíöèðîâàíèå ñ ïåðèîäîì 30 äíåé
seasonal_diff = seasonal_ts[30:] - seasonal_ts[:-30]
plt.figure(figsize=(12, 6))
plt.plot(time, seasonal_ts, label='Èñõîäíûé ðÿä')
plt.plot(time[30:], seasonal_diff, label='Ñåçîííî äèôôåðåíöèðîâàííûé ðÿä')
plt.legend()
plt.title('Ñåçîííîå äèôôåðåíöèðîâàíèå')
plt.show()